如何在R中自定义线性回归显著性水平及代码并生成美观表格

核心思路 #

借助modelsummary包直接自定义显著性阈值与对应代码，无需事后手动判断，生成符合需求的美观回归结果表格。

步骤与代码 #

1. 安装并加载工具包

install.packages("modelsummary")
library(modelsummary)

2. 计算Bonferroni校正后的显著性阈值

p_sig <- 0.05 / 12  # 约等于0.0041

3. 自定义显著性代码规则
   按p值严格程度排序，将校正后的阈值分配给指定符号（如*.*）：

custom_stars <- c(
  "***" = 0.001,       # 最严格阈值
  "*.*" = p_sig,       # Bonferroni校正后的显著性阈值
  "**" = 0.01,         # 次严格阈值
  "*" = 0.05           # 常规阈值
)

4. 拟合模型并生成带自定义显著性的表格

# 拟合线性回归
fit <- lm(y ~ x1 + x2, data = mydata)

# 生成美观表格，自动应用自定义显著性规则
modelsummary(
  fit,
  stars = custom_stars,
  title = "线性回归结果（Bonferroni校正显著性）",
  coef_map = c("x1" = "自变量1", "x2" = "自变量2", "(Intercept)" = "截距")
)

说明 #

• modelsummary会自动匹配你定义的p值阈值，为对应系数标注指定的显著性代码
• 可通过coef_map自定义变量名称、title添加表格标题，进一步优化表格可读性

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