data(dataMultiv)
dataMultiv$timegap<-dataMultiv$t1-dataMultiv$t0
## Calendar-time ##
MultivSpli <-multivePenal(Surv(timegap,deltaEvent1)~cluster(id)+var1+var2+terminal(deltadc)+event2(deltaEvent2),formula.terminalEvent=~var1,formula2=~var1+var2+var3,data=dataMultiv,Frailty=TRUE, recurrentAG=F,cross.validation=T,n.knots=c(4,4,4),kappa=c(1,1,1),maxit=20,hazard="Splines",initialize=T)
## summary a fit
summary(MultivSpli)
Run the code above in your browser using DataLab