3

Rで分散が等しくないANOVAを実行する方法があるかどうか疑問に思っていましたか?

次の例を想像してください。

x <- c(10,11,15,8,16,12,20)
y <- c(10,14,18,25,28,30,35)

d <- c(x,y)
f <- as.factor(c(rep("a",7), rep("b",7)))

# Unequal variance:

t.test(x,y)$p.value
t.test(d~f)$p.value

# Equal variance:

t.test(x,y, var.equal=TRUE)$p.value
t.test(d~f, var.equal=TRUE)$p.value

anova(lm(d~f))[[5]]
summary(aov(lm(d~f)))[[1]][5]
summary(lm(d~f))[[4]][8]

この例からわかるように、R で ANOVA を実行するさまざまな方法は、2 つのグループのみの場合、等分散の t.test によって取得されたものと常に同じ p 値になります。繰り返しますが、分散が等しくない ANOVA を実行する方法はありますか?

4

1 に答える 1

5

この場合には、oneway.test()

R> oneway.test(d~f)

    One-way analysis of means (not assuming equal variances)

data:  d and f 
F = 6.631, num df = 1.000, denom df = 8.339, p-value = 0.03179
于 2012-08-05T13:53:43.363 に答える