1

拜托,我在尝试适应时变 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 命令?

4

1 回答 1

0

没有看到您的数据就很难说head(hip),但错误消息表明您混淆了开始和结束时间 - 您可以在 github 的代码中看到该消息survSplit如果您没有开始时间,Indeed 会创建一个开始时间并"tstart"默认调用它。除非您使用旧版本的包,否则不要使用eventAndstart参数,因为它们已被弃用并且仅用于遗留使用。给episodeco 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)

那应该这样做。

于 2016-10-02T14:15:06.560 回答