gtsummary 中的似然比检验 p 值

Question

如何在 gtsummary 输出表中加入似然比检验 p 值？

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") 添加轻轨 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 (v2.0.0) 于 2021-05-12 创建}