gtsummary的tbl_regression函数能否输出稳健标准误？

gtsummary稳健标准误的计算逻辑与有效性说明 #

结论先行

你得到的稳健标准误输出是完全真实有效的，计算方法就是你预期的HC1类型稳健标准误。

• gtsummary::tbl_regression()的robust参数底层直接调用sandwich包的相关函数实现稳健协方差矩阵计算：当你传入字符串类型的稳健类型（比如"HC1"）时，参数会直接透传给sandwich::vcovHC()的type参数，计算逻辑和你用jtools::summ()调用sandwich的逻辑完全一致，所以二者输出结果没有差异。
• 你可以自行通过以下步骤验证结果一致性：
  1. 手动提取glm模型的HC1稳健协方差矩阵：robust_vcov <- sandwich::vcovHC(model, type = "HC1")
  2. 计算系数的稳健标准误：sqrt(diag(robust_vcov))
  3. 和tbl_regression()输出表中的标准误列做对比，二者数值完全相等
• 关于文档说明：旧版本的gtsummary确实未在公开文档中完整覆盖该参数的说明，v1.6.0及之后的版本已经补充了robust参数的详细说明，该参数除了支持传入sandwich兼容的稳健类型字符串，还支持直接传入用户自定义的协方差矩阵，适配特殊场景的估计需求。

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