用於測試組內均值顯着性的 R 函數

Question

我的問題是關於這個虛構的數據。 我想測試一下一起考慮的三個均值（V1、V2 和 V3）之間是否存在顯着差異。 在 R 中並測試 v1 的平均值是否與 V2 有顯着差異。

id <- c(1,2,3,4,5,6,7,8,9,10)
V1<- c(50,  42, 58, 56, 25, 85, 12, 23, 89, 52)
V2<- c(65,  63, 52, 45, 89, 58, 74, 51, 26, 25)
V3<- c(68,  95, 62, 14, 12, 25, 48, 56, 32, 57)
sex <- c("F","F","F","F","F","M","F","F","M","M")
data<- data.frame(id,V1,V2,V3,sex) 

我嘗試使用方差分析但沒有成功

Answer 1

如果要使用 anova()，則需要使用 lm() 包裝公式。

id <- c(1,2,3,4,5,6,7,8,9,10)
V1<- c(50,  42, 58, 56, 25, 85, 12, 23, 89, 52)
V2<- c(65,  63, 52, 45, 89, 58, 74, 51, 26, 25)
V3<- c(68,  95, 62, 14, 12, 25, 48, 56, 32, 57)
sex <- c("F","F","F","F","F","M","F","F","M","M")
data<- data.frame(id,V1,V2,V3,sex) 
anova(lm(id ~ V1 + V2 + V3, data = data))
Analysis of Variance Table

Response: id
           Df Sum Sq Mean Sq  F value   Pr(>F)  
 V1         1  0.438  0.4382   0.0751  0.79330  
 V2         1 29.750 29.7497   5.0959  0.06478 .
 V3         1 17.285 17.2846   2.9607  0.13610  
 Residuals  6 35.028  5.8379                  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1