如何使用 tbl_summary 显示方差分析的检验统计量（F 值）

Question

这是我用来输出方差分析摘要 table 的代码：

hiv %>% 
  select(education, sexfirsttime) %>% 
  mutate(education=
  factor(education, levels= c("no education", "primary","secondary","college"))) %>% 
  tbl_summary(missing="no", 
              by=education,
    statistic = all_continuous() ~"{mean} ({sd})", 
    label = sexfirsttime ~ "Age of first time sex") %>% 
  add_p(test= all_continuous() ~ "aov") %>% 
  modify_header(statistic ~ "**Test Statistic**")

执行代码后，我收到以下错误消息：错误：update= 参数输入错误。 Select 来自 '变量'、'test_name'、'var_type'、'var_label'、'row_type'、'标签'、'stat_1', 'stat_2', 'stat_3', 'stat_4', 'test_result', 'p.value'

当我尝试用 test_result 替换 modify_header 中的 statistic 时，我得到的输出结果很奇怪，如图所示。

我刚开始使用 gtsummary。任何帮助将不胜感激。谢谢。

Answer 1

使用最新版本的 gtsummary 并重试。在最新版本中，“aov”测试的处理与其他测试更加一致，包括返回“统计”列。

library(gtsummary)
packageVersion("gtsummary")
#> [1] '1.5.0'

tbl <- 
  trial %>%
  select(grade, age, marker) %>%
  tbl_summary(
    by = grade, 
    missing = "no"
  ) %>%
  add_p(all_continuous() ~ "aov") %>%
  # add a header (which also unhides a hidden column)
  modify_header(statistic ~ "**Test Statistic**") %>%
  # add a function to format the column
  modify_fmt_fun(statistic ~ style_sigfig)

^{由 reprex package (v2.0.1)}

于 2021-10-17 创建

如何使用 tbl_summary 显示方差分析的检验统计量（F 值）

How do I display test statistic (F value) for anova using tbl_summary

statistics

r

anova

gtsummary