R语言Friedman检验报错：unreplicated complete block design问题如何解决？

错误核心原因 #

Friedman检验要求无重复完全区组设计：每个区组（即你的ID）中，每个处理（TREATMENT）必须仅有1个观测值。你的数据里，每个ID-TREATMENT组合下对应多个TIME点的重复观测，直接违反了这一要求，因此触发报错。此外，你手动按ID, TREATMENT分组排名的方式完全错误——Friedman检验需要在每个区组内对不同处理的观测值排名，而非按区组+处理分组。

针对性解决方案 #

根据你的实验目的，分两种场景处理：

场景1：仅比较不同TREATMENT的整体差异（忽略TIME单独效应） #

将每个ID-TREATMENT组合下的多个TIME观测值合并为一个汇总统计量（非参数分析优先用中位数），确保每个区组内每个处理仅有1个值：

# 按ID和TREATMENT分组，计算中位数（比均值更适配非参数场景）
EXPERIMENT_summary <- EXPERIMENT %>%
  group_by(ID, TREATMENT) %>%
  summarise(SUMMARY_VALUE = median(VALUE), .groups = "drop")

# 运行符合要求的Friedman检验
friedman_result <- friedman.test(SUMMARY_VALUE ~ TREATMENT | ID, data = EXPERIMENT_summary)
print(friedman_result)

场景2：同时分析TREATMENT和TIME的效应（含交互作用） #

Friedman检验仅支持单处理因素+区组的设计，无法直接分析双因素。此时可使用rstatix包的非参数重复测量方差分析，同时纳入TREATMENT和TIME作为重复测量因素，ID作为区组：

# 执行双因素非参数重复测量方差分析
anova_result <- EXPERIMENT %>%
  anova_test(
    dv = VALUE,          # 响应变量
    wid = ID,            # 区组/重复测量单元标识
    within = c(TREATMENT, TIME),  # 重复测量的两个因素
    type = "III",        # 效应类型
    method = "friedman"  # 指定采用Friedman秩转换逻辑
  )

print(anova_result)

关键检查点 #

确认你的ID定义是否符合区组要求：区组应指代同一实验单元（如同一培养皿、同一植株），该单元需接受所有处理的观测。如果你的ID是每个时间点的重复编号（而非实验单元唯一标识），请重新定义ID为实验单元的唯一ID，否则区组划分完全错误。

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