如何按Session分组生成Wilcoxon配对检验报告？

要针对每个Session独立生成Wilcoxon配对检验报告，只需利用dplyr的分组功能，对每个Session对应的子数据集单独执行检验即可。以下是实现代码：

library(tidyverse)
library(report)

set.seed(123) # 设置随机种子，确保结果可复现
df <- tibble(Session = c(rep("1", 10), rep("2", 10)),
             scores = sample(20:70, 20),
             group = c(rep("before", 5), rep("after",5),
                       rep("before", 5), rep("after",5)))

# 按Session分组，生成对应中文检验报告
df %>%
  group_by(Session) %>%
  group_map(function(data, .key) {
    cat("\n=== 针对Session", .key$Session, "的Wilcoxon配对符号秩检验报告 ===\n")
    wilcox.test(scores ~ group, data = data, paired = TRUE) %>% 
      report(language = "zh") %>%
      print()
  })

代码说明 #

• group_by(Session)：将原数据集按Session拆分为独立的子数据集，每个子数据集对应一个Session的数据
• group_map：遍历每个分组，对每个子数据集单独执行Wilcoxon配对检验，并生成报告
• report(language = "zh")：直接指定生成中文格式的检验报告，无需额外翻译

运行结果示例 #

=== 针对Session 1 的Wilcoxon配对符号秩检验报告 ===
效应量标签遵循Funder (2019)的建议。

检验scores与group之间秩次差异的Wilcoxon配对符号秩精确检验显示，效应为正，无统计学显著性，且效应量较小（W = 10.00, p = 0.688; 秩双列相关r = 0.20, 95% CI [-0.47, 0.75]）

=== 针对Session 2 的Wilcoxon配对符号秩检验报告 ===
效应量标签遵循Funder (2019)的建议。

检验scores与group之间秩次差异的Wilcoxon配对符号秩精确检验显示，效应为正，无统计学显著性，且效应量较小（W = 9.00, p = 0.844; 秩双列相关r = 0.13, 95% CI [-0.52, 0.68]）

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