我正在尝试使用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模型不支持“mright”生存数据
答案 0 :(得分:15)
是您的因子类事件?它应该是数字或布尔值。因此,如果您的活动类别为0
或1
,请尝试:
df$Event <- as.numeric(df$Event)