如何将lme4重复测量ANOVA结果保存为APA格式表格？

方法1：用apaTables包一键生成（最简便） #

这个包专门适配lmerTest输出的Kenward-Roger自由度ANOVA结果，直接生成符合APA规范的表格，还能导出Word：

1. 安装加载包：

install.packages("apaTables")
library(apaTables)

2. 生成并导出APA格式表格：

# 生成表格并导出为Word文档，table.number设置表格编号
apa_lmer_anova(a1.a, filename = "RM_ANOVA_APA.doc", table.number = 1)

生成的表格会自动处理P值格式（如< .001），并保留Kenward-Roger自由度的计算结果，完全符合APA要求。

方法2：sjstats+flextable自定义格式（灵活度高） #

如果需要自定义表格样式，可以用sjstats提取标准化统计量，再用flextable调整成APA格式：

1. 安装加载依赖包：

install.packages(c("sjstats", "flextable", "dplyr"))
library(sjstats)
library(flextable)
library(dplyr)

2. 提取并整理统计数据：

# 提取包含Kenward-Roger自由度的ANOVA统计量
anova_data <- tidy_anova(a1.a, ddf = "Kenward-Roger")

# 调整P值为APA规范格式
anova_data <- anova_data %>%
  mutate(p.value = case_when(
    p.value < .001 ~ "< .001",
    p.value >= .001 ~ paste0("= ", round(p.value, 3))
  ))

# 重命名列名适配APA术语
colnames(anova_data) <- c("效应", "平方和", "自由度", "均方", "F值", "P值")

3. 渲染并导出表格：

# 生成flextable并美化
apa_table <- flextable(anova_data) %>%
  set_header_labels(
    效应 = "Effect",
    平方和 = "SS",
    自由度 = "df",
    均方 = "MS",
    F值 = "F",
    P值 = "p"
  ) %>%
  bold(part = "header") %>%
  align(align = "center", part = "all") %>%
  autofit()

# 导出为Word文档
save_as_docx(apa_table, path = "Custom_RM_ANOVA_APA.doc")

修复你之前用report包的问题 #

report包输出的P值默认格式不符合APA要求，你可以在转换为数据框后手动调整：

library(dplyr)
stats.table <- as.data.frame(report(a1.a))

# 调整P值列格式
stats.table <- stats.table %>%
  mutate(p = case_when(
    p < .001 ~ "< .001",
    p >= .001 ~ paste0("= ", round(p, 3))
  ))

# 重新渲染表格
flextable(stats.table) %>% autofit()

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