R中Friedman检验报错‘not an unreplicated complete block design’的解决咨询

报错原因 #

Friedman检验要求无重复完全区组设计：每个区组（此处为Species）内，每个处理（Category_ID）只能对应1个观测值。你的数据中，同一Species+Category_ID组合存在多条重复记录（比如Common hawthorn+UN对应多个Total_visits值），不符合检验的设计前提，因此触发报错。

两种解决方案 #

根据你的分析需求选择对应方法：

方案1：支持重复观测的扩展Friedman检验

如果实验设计本身就是每个物种在每个地点分类下有重复采样，推荐使用rstatix包的friedman_test()函数，它专门适配带重复的区组设计：

• 安装并加载包：

  install.packages("rstatix")
  library(rstatix)
  

• 执行检验：

  friedman_test(Total_visits ~ Category_ID | Species, data = diversity_columns)
  

方案2：汇总数据后使用标准Friedman检验

如果分析目标是比较不同地点分类下物种的访花次数集中趋势，可先对重复数据汇总，再用标准检验：

1. 用dplyr汇总每个组合的中位数（非参数分析更推荐中位数）：

   library(dplyr)
   summary_data <- diversity_columns %>%
     group_by(Species, Category_ID) %>%
     summarise(Median_visits = median(Total_visits), .groups = "drop")
   

2. 转换为宽格式（符合标准Friedman检验的输入要求）：

   wide_data <- summary_data %>%
     tidyr::pivot_wider(names_from = Category_ID, values_from = Median_visits)
   

3. 执行标准Friedman检验：

   friedman.test(as.matrix(wide_data[, -1]))
   

额外说明 #

你之前尝试的as.matrix转换、将变量转为因子的操作，都未解决核心问题——重复观测的存在，因此无法消除报错。

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