2011-11-30 2 views
7

R을 시도하는 중. 오래된 통계 텍스트의 질문은 다른 건설 현장에서 휴식 시간에 차이가 있는지를 알고 싶습니다. 문제는 각 사이트마다 다른 수의 직원을 고용하기로 한 본문 내용입니다. 그래서, 나는 고집스럽고 샘플 크기가 다른 ANOVA를 사용하여 도움을 요청합니다.동일하지 않은 샘플 크기로 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) 

답변

16

이 하나

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 
+10

멋지게 완료되었습니다. 비록 샘플의 크기가 같거나 다른지 여부는 방법이 완전히 동일하기 때문에 OP가 올라간 불평등 샘플 크기는 중요하지 않다고 말할 수도 있습니다. – Aaron