中位生存率

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

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

Answer 1

pkg：survival中没有surv函数。它是Surv。

?Surv

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