R Survminer Pairwise Survival difference 变长误差

Question

您好，我正在使用 survival 包和 survminer 包进行一些分析。我在使用 survminer 包中的 pairwise_survdiff 时遇到问题。

这里是模仿我自己的一些模拟数据。

set.seed(123)
library(survival)
library(survminer)
size <-  30
treatment <- c(rep(1, size),
               rep(2, size),
               rep(3, size))
deathtime <- c(rexp(size, rate = 1),
               rexp(size, rate = 2),
               rexp(size, rate = 3))
censor <- c(rbinom(size, 1, .1),
            rbinom(size, 1, .1),
            rbinom(size, 1, .1))


df <- data.frame(treatment, deathtime, censor)

进行生存分析

surv.obj<- Surv(df$deathtime, abs(df$censor -1))

survdiff(surv.obj ~ treatment, data = df)

现在尝试成对

pairwise_survdiff(surv.obj ~ treatment, data = df) # from survminer

length(df[,"treatment"]) == length(surv.obj)

我不明白是什么问题。

Answer 1

尝试不使用 df$

pairwise_survdiff(Surv(deathtime, abs(censor -1)) ~ treatment, data = df)

-输出

Pairwise comparisons using Log-Rank test 

data:  df and treatment 

  1      2     
2 0.1389 -     
3 0.0022 0.0364

P value adjustment method: BH

R Survminer Pairwise Survival difference 变长误差

R Survminer Pairwise Survival difference Variable length error

r

survival