/*** a penalized survival model ***/
data(test)
fit.su <- shr(Hist(time=list(l,r),id)~cov,data=test,hazard="Splines")
# no confidence bands
plot(fit.su,conf.bands=FALSE)
# confidence bands
plot(fit.su)
Run the code above in your browser using DataLab