gtsummary：单变量和多变量

Question

大家，最近我用gtsummary中的coxph用R包gtsummary创建单变量和多变量模型，但是出现如下错误，

add_global_p：已计算变量 add_global_p(include = c("Age", "Gender", "Path_T", "Path_N", "Path_M", "Path_stage", "TP53")) 的全局 p 值 和 car::Anova(x$model_obj, type = "III") x add_global_p() 使用 car::Anova() 计算全局 p 值， 并且函数在计算 p 值时返回错误。 car::Anova() 是否支持您的模型类型？ df.terms.default(mod) 中的错误：模型有别名项； df 模棱两可。

我的代码如下，

library(gtsummary)
data<-read.table("patient_data.txt",sep="\t",header=T)
data1<-na.omit(data)
    library(survival)
      tbl_uvsurv <- select(data1,everything()) %>%
         tbl_uvregression(
            method=coxph,
             y=Surv(Time,Status),
             exponentiate=TRUE
         )  %>%
      add_global_p() %>%
      bold_p(t=0.05)%>%
     bold_labels()

多变量表：

library(survival)
  tbl_mvsurv <-coxph(
   Surv(Time,Status) ~.,
        data=data1
   ) %>%
  tbl_regression(
      exponentiate=TRUE
   ) %>%
   add_global_p() %>%
  bold_p(t=0.05)%>%
 bold_labels()

你能帮帮我吗？提前谢谢你。

Answer 1

当您提供可重现的示例时，我们可以最好地回答您的问题，这意味着提供我们可以使用公开可用的数据集（或模拟数据集）在我们的机器上运行的代码。下面是两个示例：一个是多变量模型，一个是单变量模型表。

library(survival)
#> Warning: package 'survival' was built under R version 4.1.1
library(gtsummary)
packageVersion("gtsummary")
#> [1] '1.4.2'
# multivariable model
tbl1 <-
  coxph(Surv(ttdeath, death) ~ age + trt, trial) %>%
  tbl_regression(exponentiate = TRUE) %>%
  add_global_p()
#> add_global_p: Global p-values for variable(s) `add_global_p(include = c("age",
#> "trt"))` were calculated with
#>   `car::Anova(x$model_obj, type = "III")`

# univariable model
tbl2 <-
  trial %>%
  select(ttdeath, death, age, trt) %>%
  tbl_uvregression(
    y = Surv(ttdeath, death),
    method = coxph,
    exponentiate = TRUE
  ) %>%
  add_global_p()
#> add_global_p: Global p-values for variable(s) `add_global_p(include = c("age",
#> "trt"))` were calculated with
#>   `car::Anova(mod = x$model_obj, type = "III")`

^{由reprex package (v2.0.1) 于 2021-09-19 创建}