Differences

This shows you the differences between two versions of the page.

--- c:ms:2023:w07_anova_note [2023/04/18 13:36] – [post hoc test: output] hkimscil
+++ c:ms:2023:w07_anova_note [2023/06/04 12:38] (current) – [Post hoc test] hkimscil
@@ Line 434: / Line 434: @@
 d.bc
-# se
+# mse (ms within) from the a.res.sum output
-# from the a.res.sum output
+# a.res.sum == summary(aov(values ~ group, data=comb3))
 a.res.sum
 # mse = 50
 mse <- 50
+# 혹은 fansy way from comb3 data.frame
+# 15 는 각 그룹의 df
+sse.ch <- sum(tapply(comb3$values, comb3$group, var)*15)
+sse.ch
+mse.ch <- sse.ch/45
+mse.ch
 se <- sqrt(mse/length(A))
@@ Line 470: / Line 476: @@
 <code>
 plot(TukeyHSD(a.res, conf.level=.95), las = 2)
+</code>
+<code>
+pairwise.t.test(comb3$values, comb3$group, p.adj = "bonf")
 </code>
 ===== post hoc test: output =====