This contains a simulated sample of of 800 subjects and 1652 observations. This dataset can be used to illustrate how to fit a joint multivariate frailty model. Two gaussian correlated random effects were generated with mean 0, variances 0.5 and a correlation coefficient equals to 0.5. The coefficients \(\alpha_1\) and \(\alpha_2\) were fixed to 1. The three baseline hazard functions followed a Weibull distribution and right censoring was fixed at 5.
data(dataMultiv)
This data frame contains the following columns:
identification of patient
number of observation for a patient
start of interval
end of interval (death or censoring time)
recurrent of type 1 status (0:no, 1:yes)
recurrent of type 2 status (0:no, 1:yes)
censoring status (0:alive, 1:death)
dichotomous covariate (0,1)
dichotomous covariate (0,1)
dichotomous covariate (0,1)
time to event