生存无法识别右删失数据
Posted
技术标签:
【中文标题】生存无法识别右删失数据【英文标题】:Survival not recognizing right censored data 【发布时间】:2017-07-03 18:21:50 【问题描述】:我正在尝试使用 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)
survdiff(surv.dfs ~ group) 中的错误:仅右删失数据
> coxph(surv.dfs ~ group)
coxph(surv.dfs ~ group) 中的错误: Cox 模型不支持“正确”的生存数据
【问题讨论】:
欢迎来到 ***。你能提供一个reproducible R example吗?dput(surv.dfs)
应该对此有所帮助。
【参考方案1】:
你的事件是一个因素类吗?它应该是数字或布尔值。因此,如果您的 Event 类是 0
或 1
,请尝试:
df$Event <- as.numeric(df$Event)
【讨论】:
谢谢你,成功了。我的事件被归类为一个因素。以上是关于生存无法识别右删失数据的主要内容,如果未能解决你的问题,请参考以下文章
生存分析——KM生存曲线hazard比例PH假定检验非比例风险模型(分层/时变/参数模型)