R 中参数化生存模型的校准斜率

Question

我有两个用于weibul分布参数生存分析的数据集：一个用于训练，另一个用于测试。 我想评估测试数据集中的校准。 我尝试了 rms 包，但我不知道如何应用测试数据。

例如，

train <- lung[1:100,]
test <- lung[101:nrow(lung),]
s <- survreg(Surv(time, status) ~ sex, data = train, dist = "weibul")

#for discrimination
rcorr.cens(predict(s, test), Surv(test$time, test$status)) 

#for calibration
??????????