r - 生存无法识别正确的审查数据

Question

我正在尝试使用 survival 包在 R 中运行生存分析。我只有右审查数据(事件时间，然后是事件“1”或审查“0”)。我使用 Surv 函数创建了一个生存对象。但是，survdiff 和 coxph 给出错误，表明它们不将其识别为右删失数据。有什么建议吗？

> surv.dfs <- Surv(DaysFromTx,Event)
> surv.dfs
[1]   99:1  334:1 1024+  1341+   210+  1069+   890+  1242+   255+   228+   349+   300+   717+     2+   657+   995+   491+  1544+   265:1  440+   362:1  845+ 
[23]  669+  1176+   718+   768+  1171+  2276+  1152+   207+  1138+  1002+   942+   644+  1110+   179:1 1535+   841+   923+   904+   367:1  959+   746+  1256+ 
[45]   83:1  439:1   69+   449+   591+   983+   787+   704+   825+   747+    28+    41+   907+   181:1  371+   388+   166:1  702+   647+   944+   903+   797+ 
[67] 1095+   770:1 1118+    63:1 1762+  1662+   127:1  634+   312+   483+ 
> survdiff(surv.dfs ~ group)

Error in survdiff(surv.dfs ~ group) : Right censored data only

> coxph(surv.dfs ~ group)

Error in coxph(surv.dfs ~ group) : Cox model doesn't support "mright" survival data

Answer 1

您的事件属于因子类吗？它应该是数字或 bool 值。因此，如果您的事件类是 0 或 1，请尝试:

df$Event <- as.numeric(df$Event)