拟合时变cox模型

Question

请，我一直有一些挑战，试图适应一个时变的cox模型。确切地说，在我阅读了使用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指定一个名称，以便在co/ph模型中作为交互使用，那么应该就可以了。

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)

应该可以了