假设参数生存模型:
surv <- Surv(time, event)
surv.reg.t <- survreg(surv ~ 1, dist = "t", data = myData)
可以通过
预测一个给定病例的百分位数pct=1:9000/10000
predict.t = predict(surv.reg.t, newdata=data.frame(F = "X"), type='quantile', p=pct, se=TRUE)
这是How to predict survival probabilities in R和How to make prediction in survival analysis using R中描述的方式。在这些答案中,给出了幸存者函数的值,并预测了事件时间。
但是,如果在一些事件时间x
下预测幸存者功能?有没有办法使用survreg.predict
?如果不是,当使用分布surv.reg.t
时,如何使用dist="t"
的系数执行此操作?