Differences

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

--- c:ms:2023:w07_anova_note [2023/04/18 22:31] – hkimscil
+++ c:ms:2023:w07_anova_note [2023/06/04 21:38] (current) – [Post hoc test] hkimscil
@@ Line 125: / Line 125: @@
 </code>
-====== Output ======
+===== ANOVA in R: Output =====
 <code>
 > #
@@ Line 433: / 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 467: / Line 474: @@
 </code>
-====== post hoc test: output ======
+<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 =====
 <code>
 > m.a
@@ Line 544: / Line 558: @@
 </code>
+{{:c:ms:2023:pasted:20230418-223608.png}}
 ====== R square or Eta square ======
 SS toal
@@ Line 581: / Line 596: @@
 </code>
-====== R square: output ======
+===== R square: output =====
 <code>
 > ss.tot