如何在 R 中以不相等的样本量进行单向方差分析？答案

【问题标题】：How to do one-way ANOVA in R with unequal sample sizes?如何在 R 中以不相等的样本量进行单向方差分析？
【发布时间】：2012-01-09 08:58:03
【问题描述】：

尝试学习 R。来自旧统计文本的一个问题想知道不同建筑工地的休息时间是否存在差异。麻烦的是，文本决定每个站点雇用不同数量的工人。因此，我被困住并寻求帮助使用具有不等样本量的方差分析。

site1 <- c(34,25,27,31,26,34,21)
site2 <- c(33,35,31,31,42,33)
site3 <- c(17,30,30,26,32,28,26,29)
site4 <- c(28,33,31,27,32,33,40)

【问题讨论】：

标签： r size sample anova

【解决方案1】：

试试这个

site1 <- c(34,25,27,31,26,34,21)
site2 <- c(33,35,31,31,42,33)
site3 <- c(17,30,30,26,32,28,26,29)
site4 <- c(28,33,31,27,32,33,40)

Data <- data.frame(
       Y=c(site1, site2, site3, site4),
       Site =factor(rep(c("site1", "site2", "site3", "site4"), times=c(length(site1), length(site2), length(site3), length(site4))))
       )

数据

  Y  Site
 34 site1
 25 site1
 27 site1
 31 site1
 26 site1
 34 site1
 21 site1
 33 site2
 35 site2
 31 site2
 31 site2
 42 site2
 33 site2
 17 site3
 30 site3
 30 site3
 26 site3
 32 site3
 28 site3
 26 site3
 29 site3
 28 site4
 33 site4
 31 site4
 27 site4
 32 site4
 33 site4
 40 site4

方差分析

fm1 <- aov(Y~Site, data=Data)
anova(fm1)

输出

Analysis of Variance Table

Response: Y
          Df Sum Sq Mean Sq F value  Pr(>F)  
Site       3 212.35  70.782  3.4971 0.03098 *
Residuals 24 485.76  20.240                  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

【讨论】：

做得很好。尽管您可能想提一下，导致 OP 的样本量不等并不重要，因为无论样本量是否相等，方法都是完全相同的。