如何在R中进行不等样本大小的单向方差分析？

Question

我想学习R语言。一本旧的统计教材上的问题想要知道在不同施工现场的休息时间是否有差异。问题在于，该教材认为每个施工现场雇用了不同数量的工人。因此，我被卡住了，需要使用具有不均等样本大小的ANOVA方法来寻求帮助。

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)

- tora0515

Answer 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

方差分析（ANOVA）

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 

- MYaseen208

---

做得不错。虽然你可能想提到，让被提问者困扰的样本不均衡并不重要，因为方法在样本大小相等或不相等时完全相同。 - Aaron left Stack Overflow