data("dataDIVAT2")
# The training of the gompertz model with the first 400 patients
formula<-Surv(times,failures) ~ age + hla + retransplant + ecd
data<-dataDIVAT2[1:400,]
model <- LIB_PHgompertz(formula, data=data)
# The prognostic capacities from the same training sample
summary(model)
# The prognostic capacities from a validation of the next 150 patients
# (up to 4 years for several indicators)
#newdata<-dataDIVAT2[401:550,]
#summary(model, pro.time=4, newdata=newdata)
Run the code above in your browser using DataLab