## ----include = FALSE----------------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

## ----setup--------------------------------------------------------------------
library(autorelevate)

## ----data---------------------------------------------------------------------
data(bladder_cancer)
summary(bladder_cancer)

## ----ttt----------------------------------------------------------------------
ttt_plot(bladder_cancer)

## ----compare_families---------------------------------------------------------
family_table <- compare_families(bladder_cancer)
print(family_table)

## ----fit, warning=FALSE, message=FALSE, fig.width=8, fig.height=4-------------
best_dist <- family_table$Family[1]
fit_best <- fit_autorelevate(bladder_cancer, dist = best_dist, method = "mle")
summary(fit_best)
plot(fit_best)

## ----fit_all_methods----------------------------------------------------------
fit_all_methods(bladder_cancer, dist = best_dist)

## ----density------------------------------------------------------------------
x <- seq(0.01, 5, by = 0.05)
plot(x, dautorelevate(x, dist = "weibull", p1 = 0.5, p2 = 1.5), type = "l",
     ylab = "Density", main = "Autorelevated Weibull density")

## ----hazard-------------------------------------------------------------------
plot(x, haautorelevate(x, dist = "weibull", p1 = 0.5, p2 = 0.8), type = "l",
     ylab = "Hazard rate", main = "Upside-down bathtub hazard (beta = 0.8)")