在r

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 here的代码中看到该消息。 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)

那应该这样做。