r - 竞争风险/生存分析 R

Question

关于处理 R 中的竞争风险的另一个问题，

在两个不同的治疗组中，无论是否考虑竞争风险“局部复发”，我都必须分析“转移或死亡”事件之前的生存函数。我有这个：状态status_meta的变量向量：1如果在时间“t”转移，0如果在其他情况下status_death：1如果在时间“t”死亡，0如果在其他情况下status_loc：1如果在时间“t”局部复发，如果是0

直到事件的时间变量向量：time_meta、time_death、time_loc

我需要合并状态（或时间）的向量，所以我只能有一个状态向量：1 表示转移或死亡，2 表示局部复发，0 表示其他情况（三个基线中的每一个中都存在非信息性审查向量）我认为可以使用 apply() 和 which.min 但我不知道如何使用它（使用时间或状态？）我试过这个：

data$etat_comp<-apply(cbind(data$time_meta,data$time_death,data$time_loc),1,
which.min)-1

但我有这个消息错误

Error in apply(cbind(data$time_meta, data$time_death, data$time_loc), 1,  : 
non numeric argument for binary operator

我也尝试为 cbind na.omit 但也不起作用

我还需要修复时间向量，因此我将有一个向量表示直到 mestatse 或死亡的时间，另一个向量表示直到局部复发的时间。

这些步骤是必要的，因为我必须考虑主要事件并考虑到 Fine and Gray 或 Cox ph 回归中的竞争风险（局部复发）来选择哪种状态（元或死亡）。

谢谢您的帮助

标签： rsurvival-analysis