2011-11-30 17 views
7

Đang cố gắng tìm hiểu R. Một câu hỏi từ văn bản thống kê cũ muốn biết nếu có sự khác biệt về thời gian nghỉ tại các trang web xây dựng khác nhau. Vấn đề là, văn bản quyết định rằng mỗi trang web sử dụng một số lượng công nhân khác nhau. Vì vậy, tôi đang bị mắc kẹt và tìm kiếm sự trợ giúp bằng cách sử dụng ANOVA với kích thước mẫu không bằng nhau.Làm thế nào để thực hiện ANOVA một chiều trong R với kích thước mẫu không bằng nhau?

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) 

Trả lời

16

Hãy thử điều này một

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)))) 
     ) 

liệu

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) 

Đầu ra

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

Thực hiện tốt. Mặc dù bạn có thể muốn đề cập đến rằng các kích thước mẫu không bằng nhau mà vấp phải OP không quan trọng, vì phương pháp này hoàn toàn giống nhau cho dù kích thước mẫu có bằng nhau hay không. – Aaron