如何在gtsummary的tbl_regression中显示beta估计的标准误？

当然可以实现！在gtsummary::tbl_regression()里把标准误（SE）加入输出表格有两种实用的方式，我结合你的示例代码给你演示：

方法1：将SE与估计值放在同一列 #

如果希望把SE直接跟在估计值后面（比如显示为“系数 (SE)”），可以通过自定义statistic参数来实现：

library(gtsummary)

# 拟合你的示例模型
test <- with(trial, lm(marker~trt+age))

# 自定义统计量格式，包含估计值和SE
tbl_regression(test,
               statistic = ~ "{estimate} ({std.error})",
               conf.int = FALSE) # 若不需要置信区间可关闭，需要的话保留并调整statistic格式

如果想同时保留置信区间，可以把statistic参数改成这样：

statistic = ~ "{estimate} ({std.error}) [{conf.low}, {conf.high}]"

这样列里会显示「估计值 (SE) [置信区间下限, 置信区间上限]」的格式。

方法2：单独添加一列显示SE #

要是你希望SE作为独立的一列展示，可以用add_stat()函数新增统计量列：

tbl_regression(test) %>%
  add_stat(
    # 提取并格式化标准误（这里保留3位小数，可按需调整）
    fns = ~ tibble(stat_0 = paste0("(", round(.std.error, 3), ")")),
    # 指定放在估计值列的后面
    location = "after estimate"
  ) %>%
  # 给新列设置表头
  modify_header(stat_0 = "SE")

小提示 #

不管是线性回归还是其他模型（比如逻辑回归），这个思路都适用。你可以通过查看?tbl_regression文档里的「Statistic Functions」部分，找到对应的统计量占位符（比如.std.error就是标准误的专属占位符）。

内容的提问来源于stack exchange，提问作者West_End_Line