data(dataAdditive)
modAdd <- additivePenal(Surv(t1,t2,event)~cluster(group)+var1+slope(var1),
correlation=TRUE,data=dataAdditive,n.knots=8,kappa1=862,
hazard="Splines")
# It takes around 4 minutes to converge. 'var1' is boolean as a treatment variable.
summary(modAdd)Run the code above in your browser using DataLab