在r中拟合时变cox模型

Question

拜托，我一直在尝试适应时变的考克斯模型。 确切地说，在我阅读了使用了survSplit命令的示例后，我发现很难得到相同的结果。 我最终得到一个错误或另一个错误，这是一个示例和错误消息。

hip2=survSplit(Surv(Time,status)~Operation.Performed,
           data=hip,cut=c(200,500,1000,1600),
           end='Time',event='status',start="start")

Error in survSplit(Surv(Time, status) ~ Operation.Performed, data = hip,  : 
start time must be < stop time

注意：我的数据集最初没有标记为start的列。 但是我看到它会被命令创建，我不知道这是不对的？

拜托，有什么东西没能正确地得到正确的结果??? 或者更好的是有更好的方法来实现survSplit命令？

Answer 1

很难没有看到你的数据说head(hip) ，但该错误消息指出您具有混淆的开始和结束时间-你可以在github上看到消息中的代码在这里 。 如果没有，则确实会创建一个开始时间，默认情况下将其survSplit为"tstart" 。 除非您使用旧版本的软件包，否则请不要使用event并start参数，因为它们已被弃用，仅供传统使用。 给episode命名以在合作模型中用作交互，你应该没问题 。

hip2 <- survSplit(Surv(Time, status) ~., data=hip, cut=c(200,500,1000,1600),episode="tgroup")

model.coxph <- coxph(Surv(tstart, Time, status) ~ Operation.Performed:strata(tgroup), data=hip2)

这应该做到这一点。