Last chance! 50% off unlimited learning
Sale ends in
cluster()
function in the context of survival joint models
for clustered data. This function identifies subject index. It is used on the right hand
side of a 'frailtyPenal' formula. Using num.id()
in a formula implies that
a joint frailty model for clustered data is fitted (Rondeau et al. 2011).num.id(x)
frailtyPenal
data(readmission)
#-- here is a generated cluster (31 clusters of 13 subjects)
readmission <- transform(readmission,group=id%%31+1)
###--- Joint Frailty model for clustered data ---###
#-- splines
joi.clus <- frailtyPenal(Surv(t.start,t.stop,event)~
cluster(group)+num.id(id)+dukes+charlson+sex+chemo+terminal(death),
formula.terminalEvent=~dukes+charlson+sex+chemo,
data=readmission,recurrentAG=TRUE,
n.knots=10,kappa=c(2.11e+08,9.53e+11))
Run the code above in your browser using DataLab