从coxph模型进行预测

Question

我很难使用coxph进行预测。

我希望在样本数据上学习Cox PH模型，然后使用从样本数据中导出的参数，如下所示：

# learn IS params
model.PH <- coxph(Surv(days.IS, outcome.IS) ~ predictor.IS)

# apply IS params to OOS data to make predictions
predictions.raw <- predict(model.PH, newdata = predictor.OS)

# binarise predictions
predictions.OS <- rep(0,length(predictions.raw))    
predictions.OS[which(predictions.raw>0)]<- 1 

# fit survival model
fittedModel <- survdiff(Surv(days.OS, outcome.OS) ~ predictions.OS)

预测器.IS具有维数Y_1 x D，预测器.OS具有维数Y_2 x D

但是，这不起作用，因为predictions.OS中的元素数量是Y_1而不是Y_2。

我做错了什么？

Answer 1

如果使用data = dfrm参数将值传递给``coxph，然后传递newdata参数以预测为具有相同列名的不同数据帧，那么这可能会更干净。我感觉你传递的矩阵可能没有与模型中找到的相同的列名。维度匹配不是predict.coxph函数的唯一要求。