如何在R中仅打印带显著性标记的lm模型摘要系数表

方法1：直接捕获并提取summary的输出文本 #

这种方法完全复用summary()的原生输出，不用自行计算显著性标记，步骤如下：

# 构建模型（你的原有代码）
frame      <- data.frame(y=rnorm(100), a=rnorm(100), b=rnorm(100))
frame["c"] <- frame["y"] + 2*rnorm(100)
model <- lm(y ~ a + b + c, data=frame)

# 捕获summary的完整输出文本
summary_text <- capture.output(summary(model))

# 定位系数表的起止位置
coef_start <- grep("Coefficients:", summary_text)
signif_start <- grep("Signif. codes:", summary_text)
coef_end <- signif_start - 1

# 提取并打印系数表及显著性说明
cat(paste(summary_text[coef_start:coef_end], collapse = "\n"), "\n")
cat(summary_text[signif_start:length(summary_text)], sep = "\n")

运行后会输出和summary(model)里完全一致的系数表，包含所有原生显著性标记。

方法2：基于系数表生成标记（复用summary的规则） #

如果需要结构化的结果而非纯文本，可以用coef(summary(model))获取系数数据，再复用summary()内置的显著性标记规则生成标记：

# 获取系数表
coef_data <- coef(summary(model))

# 生成显著性标记（和summary用的规则完全一致）
signif_levels <- c(0, 0.001, 0.01, 0.05, 0.1, 1)
signif_labels <- c("***", "**", "*", ".", " ")
signif_marks <- cut(coef_data[, 4], breaks = signif_levels, labels = signif_labels, include.lowest = TRUE)

# 合并系数表与标记
coef_with_marks <- cbind(coef_data, " " = signif_marks)

# 打印结果
print(coef_with_marks, quote = FALSE)
cat("\n---\nSignif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1\n")

这种方法得到的是带标记的矩阵，既保留了结构化数据特性，又和summary输出的标记逻辑完全一致。

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