在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)

這應該做到這一點。