R Survminer成对生存差可变长度错误

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)

-output

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