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,Frailty=TRUE,joint=TRUE,recurrentAG=TRUE,
n.knots=10,kappa1=2.11e+08,kappa2=9.53e+11)
Run the code above in your browser using DataLab