r - gsummary 中的似然比检验 p 值

Question

如何将似然比检验 p 值合并到 gtsummary 输出表中？

library(gtsummary)

    trial %>% 
      select(response, grade) %>% 
      tbl_uvregression(
        method = glm,
        y = response,
        method.args = list(family = binomial),
exponentiate = TRUE)

Answer 1

您可以使用 add_global_p(test = "LR") 添加 LRT p 值。在后台，该函数使用 car::Anova(mod = x, type = "III", test = "LR") 计算 p 值

library(gtsummary)
#> #BlackLivesMatter

tbl <-
  trial %>% 
  select(response, grade) %>% 
  tbl_uvregression(
    method = glm,
    y = response,
    method.args = list(family = binomial)
  ) %>%
  add_global_p(test = "LR")

#> add_global_p: Global p-values for variable(s) `add_global_p(include = "grade")`
#> were calculated with
#>   `car::Anova(mod = x$model_obj, type = "III", test = "LR")`

^{由 reprex package 创建于 2021-05-12 (v2.0.0)}