survdiff报错“仅含右删失数据”：surv_pvalue调用异常求助

问题详情 #

• 问题与Stack Overflow上《Survival not recognizing right censored data》帖子类似，但该帖的解决方案不适用
• 已确认删失变量为数值型，使用coxph函数时无异常
• 尝试用surv_pvalue函数获取Kaplan-Meier模型的p值用于可视化时，出现报错：

  Error in survival::survdiff(eval(fit$call$formula), data = data) : Right censored data only

• 数据集包含标准化入组日期、周期间隔，用于后续cox模型的时变协变量

所用代码 #

library(readr)
library(survival)
library(survminer)
NOBO_surv_com <- read.csv("WSIdf.csv")

# 创建删失变量（右删失）
NOBO_surv_com$censored[NOBO_surv_com$Death == 1] <- 1
NOBO_surv_com$censored[NOBO_surv_com$Death == 0 | NOBO_surv_com$Death == 2] <- 0
NOBO_surv_com$censored <- as.numeric(NOBO_surv_com$censored)

# 指定带分类协变量（年份和性别）的K-M分析模型
km_fit_Yr <- survfit(Surv(tstart, tstop, censored) ~ Yr,
                  data=NOBO_surv_com,
                  type="kaplan-meier")
print(km_fit_Yr)
surv_pvalue(km_fit_Yr, data = NOBO_surv_com)

内容的提问来源于stack exchange，提问作者Blake Baum