Generation of survival data from a time-homogeneous Markov model.
genTHMM(n, model.cens, cens.par, beta, covar, rate)
Sample size.
Model for censorship. Possible values are "uniform" and "exponential".
Parameter for the censorship distribution. Must be greater than 0.
Vector of three regression parameters for the three transitions: (beta_12,beta_13,beta_23).
Parameter for generating the time-fixed covariate. An uniform distribution is used.
Vector of dimension three. We assume an exponential baseline hazard function with constant hazard rate for each transition.
An object with two classes, data.frame
and THMM
.
For generating survival data from the THMM model, the counting process data structure must contain the following variables:
id
, time
, state
, covariate
. Each patient is identified by id.
The variable time
represents time for each interval of follow-up while variable state
denotes the state of the individual.
Variable covariate
is the (time-fixed) covariate to be studied in the regression model.
Individuals without change in the time dependent covariate are represented by two lines of data,
whereas patients with a change in the time-dependent covariate must be represented by three lines.
Jackson, C. (2011). Multi-State Models for Panel Data: The msm Package for R. Journal of Statistical Software, 38(8), 1<U+2013>28. 10.18637/jss.v038.i08
Meira-Machado, L., Cadarso-Su<U+00E1>rez, C., De U<U+00F1>a- <U+00C1>lvarez, J., Andersen, P.K. (2009). Multi-state models for the analysis of time to event data. Statistical Methods in Medical Research, 18(2), 195-222. 10.1177/0962280208092301
Meira-Machado L., Faria S. (2014). A simulation study comparing modeling approaches in an illness-death multi-state model. Communications in Statistics - Simulation and Computation, 43(5), 929-946. 10.1080/03610918.2012.718841
Meira-Machado, L., Sestelo M. (2019). Estimation in the progressive illness-death model: a nonexhaustive review. Biometrical Journal, 61(2), 245<U+2013>263. 10.1002/bimj.201700200
Therneau, T.M., Grambsch, P.M. (2000). Modelling survival data: Extending the Cox Model, New York: Springer.
# NOT RUN {
thmmdata <- genTHMM( n=100, model.cens="uniform", cens.par=80, beta= c(0.09,0.08,-0.09),
covar=80, rate= c(0.05,0.04,0.05) )
head(thmmdata, n=20L)
# }
Run the code above in your browser using DataLab