R中Friedman检验报not an unreplicated complete block design错误如何解决？

1 错误核心原因 #

首先明确R中friedman.test()的适用前提：要求数据为无重复完全区组设计，即groups（待比较的处理因素）和blocks（待控制的区组因素）的每一个交叉组合仅能存在1个观测值。你的数据中，同一treatment+day的组合有3个重复值，不符合该检验的输入要求，因此触发报错。
另外你的参数设置和研究目标不匹配：你要探究day对score的影响，应该将待比较的因素day设为groups，而不是作为blocks参数传入。

2 修正方案 #

方案一：先合并重复值适配Friedman检验 #

你可以先计算每个treatment+day组合下score的中位数/均值，把数据转换成每个交叉组合仅1个值的格式，再运行检验：

library(dplyr)
# 按treatment和day分组，计算每组score中位数，得到无重复的数据集
fried_data <- sample_data %>%
  group_by(treatment, day) %>%
  summarise(median_score = median(score), .groups = "drop")
# 运行Friedman检验：groups为待比较的day，blocks为控制的treatment
friedman3 <- friedman.test(y = fried_data$median_score, 
                           groups = fried_data$day, 
                           blocks = fried_data$treatment)
print(friedman3)

注意该方案会损失重复测量的变异信息，仅适用于探索性分析。

方案二：改用更适配重复测量设计的非参数检验 #

如果要保留所有重复数据的信息，建议使用对齐秩变换的ANOVA（ART）方法，适配含重复的非参数重复测量设计：

library(ARTool)
# 构建模型，控制treatment的影响，检验day的主效应
art_model <- art(score ~ day + (1|treatment), data = sample_data)
# 输出检验结果
anova(art_model)

额外说明 #

• 不建议使用attach()函数，容易出现变量名冲突导致的错误，直接用$符或者dplyr的管道调用变量更稳妥。
• friedman.test()的输出不需要用summary()包裹，直接print()就可以显示完整检验结果。

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