在R？

Question

我用权重的时间相关变量和时间独立的变量拟合了Coxph模型，如下所示：

id start stop death weight age smoke
1   0     1    0     60     60   0
2   0     1    0     55     57   1 
2   1     2    1     60     58   1
...

mod <- coxph(Surv(start, stop, death) ~ weight + age + factor(smoke) + cluster(id), 
             data=data)

我想绘制权重对生存曲线的影响，该图可能有2条生存曲线，一条代表权重60，另一条代表权重80，这表明增加权重对曲线的影响。

如何绘制此图？

Answer 1

您可以指定一个在不同时间点具有权重值的data.frame，并根据此数据集，可以使用 survival 包的survfit()函数来计算生存概率合适的Cox模型 but ，还可以通过适当地指定id参数。

但是，请注意，权重是潜在的和内生的时变协变量，因此从时间上依赖Cox模型在理论上不适合。您可以改为使用联合模型来计算纵向重量结果和事件发生的时间。这些模型例如在R包JM和JMbayes中实现。从拟合的关节模型中，您还可以获得（动态）生存概率。