A Case Study Using the Beta-Danish Distribution"
In BetaDanish: The Beta-Danish Distribution for Lifetime Data Analysis

knitr::opts_chunk$set(collapse = TRUE, comment = '#>')

Introduction

This vignette demonstrates a typical survival analysis workflow using the BetaDanish package.

library(BetaDanish)
library(survival)
data('remission', package = 'BetaDanish')
head(remission)

Fitting the Beta-Danish model

fit <- fit_betadanish(Surv(time, status) ~ 1, data = remission, n_starts = 1)
summary(fit)

Three-parameter submodel

fit_sub <- fit_betadanish(Surv(time, status) ~ 1, data = remission, submodel = TRUE, n_starts = 1)
compare_models(fit, fit_sub)

Diagnostic plots

plot(fit, type = 'survival')
plot(fit, type = 'hazard')

Interpretation

The fitted model can be used to estimate survival probabilities, hazard behavior, and overall model fit. Users should compare the Beta-Danish model with alternative lifetime distributions and inspect diagnostic plots before drawing final conclusions.