data(readmission)
#-- gap-time
modJoint.gap <- frailtyPenal(Surv(time,event)~cluster(id)+sex+dukes+
charlson+terminal(death),formula.terminalEvent=~sex+dukes+charlson,
data=readmission,n.knots=14,kappa1=9.55e+9,kappa2=1.41e+12,
Frailty=TRUE,joint=TRUE)
#-- calendar time
modJoint.calendar <- frailtyPenal(Surv(t.start,t.stop,event)~cluster(id)+
sex+dukes+charlson+terminal(death),formula.terminalEvent=~sex+dukes+charlson,
data=readmission,n.knots=10,kappa1=9.55e+9,kappa2=1.41e+12,Frailty=TRUE,
joint=TRUE,recurrentAG=TRUE)
#-- It takes around 1 minute to converge
summary(modJoint.gap)
summary(modJoint.calendar)Run the code above in your browser using DataLab