如何在R包lmerTest的lmer()函数结果中显示更多小数位？

我使用lmerTest包的lmer()构建混合效应模型后，调用summary()输出结果时，AIC、BIC，随机效应中的Variance、Std.Dev，以及固定效应的Estimate、Std.Error小数位数不足，无法满足报告撰写需求。尝试过以下方法均无效：

• 修改R全局选项：options(digits=3)
• 在summary()或coef(summary())中指定digits参数：

  summary(lmer, digits=3)
  coef(summary(lmer), digits=3)
  

• 对summary()或coef(summary())的结果使用round()函数：

  round(summary(lmer), digits=3)
  round(coef(summary(lmer)), digits=3)
  

• 在lmerTest包手册中搜索"digit"和"decimal"关键词未找到有效方法

可行解决方案 #

方法1：手动提取并格式化各部分结果 #

直接从模型对象中提取目标内容，手动设置小数位数：

固定效应部分

# 替换model为你的实际模型对象
fixed_coef <- coef(summary(model))
# 为Estimate和Std. Error保留4位小数（可自行调整位数）
fixed_coef[, c("Estimate", "Std. Error")] <- round(fixed_coef[, c("Estimate", "Std. Error")], 4)
print(fixed_coef)

随机效应部分

rand_vc <- VarCorr(model)
# 转换为数据框后格式化方差和标准差
rand_df <- as.data.frame(rand_vc)
rand_df[, c("vcov", "sdcor")] <- round(rand_df[, c("vcov", "sdcor")], 4)
print(rand_df)

AIC与BIC指标

aic_val <- round(AIC(model), 4)
bic_val <- round(BIC(model), 4)
cat(sprintf("AIC: %.4f\nBIC: %.4f\n", aic_val, bic_val))

方法2：使用broom.mixed包统一格式化 #

broom.mixed包专门用于整理混合效应模型的输出，支持灵活设置小数位数：

# 未安装则先安装
install.packages("broom.mixed")
library(broom.mixed)

# 提取固定效应并格式化
tidy(model, effects = "fixed") |> 
  mutate(across(c(estimate, std.error), ~round(.x, 4)))

# 提取随机效应参数并格式化
tidy(model, effects = "ran_pars") |> 
  mutate(across(c(estimate), ~round(.x, 4)))

# 提取模型拟合指标（含AIC、BIC）并格式化
glance(model) |> 
  mutate(across(c(AIC, BIC), ~round(.x, 4)))

方法3：重定义打印函数修改默认小数位数 #

如果希望直接调用summary()时就显示指定小数位数，可以重写打印方法：

# 复制原打印函数并指定默认digits
print.summary.merModLmerTest <- function(x, digits = 4, ...) {
  lmerTest:::print.summary.merModLmerTest(x, digits = digits, ...)
}

# 现在调用summary会自动用4位小数输出
summary(model)

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