生存无法识别右删失数据

Posted

技术标签:

【中文标题】生存无法识别右删失数据【英文标题】:Survival not recognizing right censored data 【发布时间】:2017-07-03 18:21:50 【问题描述】:

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

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

survdiff(surv.dfs ~ group) 中的错误:仅右删失数据

> coxph(surv.dfs ~ group)

coxph(surv.dfs ~ group) 中的错误: Cox 模型不支持“正确”的生存数据

【问题讨论】:

欢迎来到 ***。你能提供一个reproducible R example吗? dput(surv.dfs) 应该对此有所帮助。 【参考方案1】:

你的事件是一个因素类吗?它应该是数字或布尔值。因此,如果您的 Event 类是 01,请尝试:

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

【讨论】:

谢谢你,成功了。我的事件被归类为一个因素。

以上是关于生存无法识别右删失数据的主要内容,如果未能解决你的问题,请参考以下文章

生存分析——KM生存曲线hazard比例PH假定检验非比例风险模型(分层/时变/参数模型)

scipy.stats.weibull_min.fit() - 如何处理右删失数据?

R生存分析AFT

JAGS 泊松计数删失数据

怎么在SPSS中做kaplan-meier生存分析

为啥 Xcode 无法识别我的核心数据实体的属性