使用Cox模型來估計生存

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行」，「有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 

你不應該使用$提供論據監測網。在數據框的環境中構建模型非常重要。