使用考克斯模型来估计生存率

Question

我首先在R中建立一个cox模型：

test1<- test[1:20,]
model.1 <- coxph(Surv(test1$days,test1$status==1) ~ test1$MTT+test1$ADC,data=test1)

当我试图预测下一位患者的生存时：

covs1 <- data.frame(test[21,]$MTT,test[21,]$ADC)
summary(survfit(model.1, newdata= covs1, type ="aalen"))

它给了我太多的生存结果，警告就是 “'newdata'有1行，但找到的变量有20行” fyi，有20个事件，结果包含20个生存结果。

Answer 1

作为预测基础的数据框中列的名称必须与模型公式的RHS中的列名相同。除非你做这样的事情，否则我认为你不会有资格：

test1<- test[1:20,]
model.1 <- coxph( Surv(days, status==1) ~ MTT + ADC, data=test1)

covs1 <- test[21, c("MTT", "ADC")]
 # then do your prediction

您不应该使用$向Surv提供参数。在数据框的环境中构建模型非常重要。