我想知道在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中进行单因素方差分析的不同方式(仅涉及两组)总是会得到一个与等方差t检验相同的p值。那么,有没有一种方法可以进行方差不等的单因素方差分析呢?