library(frailtypack)
data(dataAdditive)
# Additive with one covariate #
modAdd1cov<-additivePenal(Surv(t1,t2,event)~cluster(group)+var1+
slope(var1),data=dataAdditive,n.knots=8,kappa1=10000)
# 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,kappa1=10000)
# 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,kappa1=10000,kappa2=10000)Run the code above in your browser using DataLab