rm(list = ls()) library(tidyverse) source("http://riseki.php.xdomain.jp/index.php?plugin=attach&refer=ANOVA君&openfile=anovakun_486.txt") # Between ----------------------------------------------------------------- # データフレームを組み上げる Example1 <- data.frame( id = rep(1:4, 3), condition = rep(1:3, each = 4), value = c(6, 6, 5, 5, 7, 4, 7, 6, 4, 3, 4, 6) ) # Factor型にする Example1$condition <- factor(Example1$condition, labels = c("control", "exp1", "exp2") ) Example1 result.lm1 <- lm(value ~ condition, data = Example1) summary(result.lm1) ggplot(Example1) + aes(x = condition, y = value, fill = condition) + stat_summary(geom = "bar", fun = "mean") anova(result.lm1) anovakun(Example1[, 2:3], "As", 3, eps = T) # 二要因 --------------------------------------------------------------------- Example2 <- data.frame( id = 1:12, num = rep(1:3, 4), temp = rep(1:2, each = 6), maker = rep(rep(1:2, each = 3), 2), value = c(13, 11, 12, 7, 6, 8, 9, 9, 9, 13, 11, 9) ) Example2$temp <- factor(Example2$temp, labels = c("Hot", "Cold") ) Example2$maker <- factor(Example2$maker, label = c("A", "B") ) Example2 result.lm2 <- lm(value ~ temp * maker, data = Example2) summary(result.lm2) anova(result.lm2) anovakun(Example2[, 3:5], "ABs", 2, 2, eps = T) # Within ------------------------------------------------------------------ Example3 <- data.frame( ID = 1:4, Time1 = c(10, 9, 4, 7), Time2 = c(5, 4, 2, 3), Time3 = c(9, 5, 3, 5) ) Example3 anovakun(Example3[, 2:4], "sA", 3, eps = T)