中位生存率

Question

我正在尝试使用 R 中的生存包对中位生存率进行对数秩检验 - 这可能吗？我正在使用代码：

survdiff(surv(Age, Survival) ~ Treatment, data = GBM)

其中年龄 = 诊断时的年龄；生存 = 是他们存活了多少年；治疗 = 是否接受治疗（编码为 1 或 2）

如果这是正确的，为什么我会收到错误：

Error in survdiff.fit(y, groups, strata.keep, rho) : 
  There is only 1 group
In addition: Warning message:
In Surv(GBM$Age, GBM$Survival) : Invalid status value, converted to NA

Treatment 是否需要编码为 0 或 1 而不是 1 或 2？

Answer 1

pkg:survival 中没有surv 函数。我是Surv。

?Surv

协变量编码为 0/1 还是 1/2 无关紧要。该估计将描述取决于一个单位差异的风险比率。如果你想设置Surv 调用来处理年龄和持续时间，那么事件时间应该是年龄+生存。当然，您还需要一个事件指示器。