# NOT RUN {
# }
# NOT RUN {
data(dataAdditive)
##-- Additive with one covariate --##
modAdd1cov <- additivePenal(Surv(t1,t2,event)~cluster(group)+var1+
slope(var1),data=dataAdditive,n.knots=8,kappa=10000,hazard="Splines")
##-- Additive with two covariates --##
set.seed(1234)
dataAdditive$var2 <- rbinom(nrow(dataAdditive),1,0.5)
modAdd2cov <- additivePenal(Surv(t1,t2,event)~cluster(group)+var1+
var2+slope(var1),data=dataAdditive,n.knots=8,kappa=10000,
hazard="Splines")
##-- Additive with 2 covariates and stratification --##
dataAdditive$var2 <- rbinom(nrow(dataAdditive),1,0.5)
modAddstrat <- additivePenal(Surv(t1,t2,event)~cluster(group)+
strata(var2)+var1+slope(var1),data=dataAdditive,n.knots=8,
kappa=c(10000,10000),hazard="Splines")
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab