문제 1
- 데이터 iris 에서 Species 에 따른 Petal.Length값에 대해 일원분산분석을 시행하고 그 결과를 해석하여라.
풀이 1
aov1 = aov(iris$Petal.Length ~ iris$Species)
summary(aov1)
p값이 유의수준보다 작게 나온다. 즉 귀무가설을 기각하고 대립가설을 채택하게 된다. -> 종에 따라 Petal.Length값의 평균이 다르다.
TukeyHSD(aov1)
사후분석을 진행했을 때 3가지 그룹 모두 p값이 유의수준보다 작게 나왔다. 즉 귀무가설을 기각하고 대립가설을 채택하면 된다.
문제 2
- 아래표를 활용하여 이원분산분석을 시행하고 그 결과를 해석하여라.
풀이 2
val = c(700, 820, 710, 830, 540, 680, 530, 710, 450, 590, 470, 590, 460, 600, 470, 610)
group1 = c(rep(c("s", "mm", "m", "b"), c(4, 4, 4, 4)))
group2 = rep(rep(c("A", "B"), c(2, 2)), 4)
dataz = data.frame(val, group1, group2)
aov2 = aov(dataz$val ~ dataz$group1*dataz$group2)
summary(aov2)
p값이 유의수준보다 크게 나온다. ->귀무가설 채택