warpbreaks 분석하는 건데 울 종류에 따라 breaks 개수 차이나는지를 아노바로 확인해야 합니다. 어떻게 해야 하나요??ㅠㅠ
아래 결과를 확인해보시기 바랍니다.
require(stats); require(graphics) summary(warpbreaks) opar <- par(mfrow = c(1, 2), oma = c(0, 0, 1.1, 0)) plot(breaks ~ tension, data = warpbreaks, col = "lightgray", varwidth = TRUE, subset = wool == "A", main = "Wool A") plot(breaks ~ tension, data = warpbreaks, col = "lightgray", varwidth = TRUE, subset = wool == "B", main = "Wool B") mtext("warpbreaks data", side = 3, outer = TRUE) par(opar) summary(fm1 <- lm(breaks ~ wool*tension, data = warpbreaks)) anova(fm1)
아래 결과를 확인해보시기 바랍니다.