Question

我在等待器官移植的患者名单上有一组患者的数据。事件是死亡时间。该集合包含以下变量：PatientID，时间，状态，年龄，性别，bmi，疾病。我通过使用加速失败时间模型建立了一个候补名单上患者生存时间的模型，该模型是年龄，性别，bmi和疾病的协变量，其中基线分布由对数正态分布给出。

fitLog <- survreg(Surv(time, status)~age+gender+bmi+disease, dist = "lognormal", data = patients)

现在，我需要预测60岁，bmi为23岁的男性纤维化（疾病= 2）患者的生存曲线，并且还需要在该曲线上添加逐点置信区间。

new_data <- with(patients, data.frame(gender = c(1),
                                  age = c(60),
                                  bmi = c(23),
                                  disease = c(2)))

plot(predict(fitLog, newdata=new_data, type="quantile", p=seq(.01,.99,by=.01)), seq(.99,.01,by=-.01), col="red", type="l", xlab = "Time (days)", ylab = "Survival probability (%)", main="Survival time")

因此，当我运行这段代码时，我得到了预测模型，但是现在我必须向该图添加逐点置信区间。我该怎么做？我已经尝试过conf.type，但是什么也没做。

编辑：我使用以下代码来计算和绘制置信区间：

pred <- predict(fitLog, newdata=new_data, 
type="quantile",p=seq(.01,.99,by=.01), se=TRUE)

plot(pred$fit - 1.96*pred$se.fit, seq(.99,.01,by=-.01),col="red", type="l", 
lty=2, xlab = "Time (days)", ylab = "Survival probability (%)", 
main="Survival time")

lines(pred$fit + 1.96*pred$se.fit,seq(.99,.01,by=-.01),col="red", type="l", 
lty=2)

lines(predict(fitLog, newdata=new_data, 
type="quantile",p=seq(.01,.99,by=.01)),seq(.99,.01,by=-.01),col="red", 
type="l")

所以现在我可以绘制置信区间了，但是现在的问题是如何使虚线持续到随访结束？

image of the plot

将逐点置信区间添加到生存曲线图的预测中

0 个答案: