Differences

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

--- c:ms:2017:schedule:week07 [2017/04/14 09:46] – hkimscil
+++ c:ms:2017:schedule:week07 [2017/04/14 12:04] (current) – hkimscil
@@ Line 1: / Line 1: @@
-Read [[:ANOVA]] for the logics of ANOVA
+====== ANOVA ======
+Read [[:ANOVA]] for the logics of ANOVA
+also read: [[:r:oneway_anova|Oneway ANOVA R]]
+{{:anova_egs.xlsx}}
 Eysenck는 (1974) 단어자료를 기억하는 정도는 단어를 외울 때 사용한 처리방법에 따라서 다르게 나타날 것이라고 예측 (p. 426) 55-65세 50명을 다섯 그룹에 무작위배치(random assignment)
   - 의도학습
@@ Line 166: / Line 169: @@
 Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
 >
+> TukeyHSD(a.out)
+  Tukey multiple comparisons of means
+% family-wise confidence level
+Fit: aov(formula = RECALL ~ GROUP, data = a.data)
+$GROUP
+    diff         lwr       upr     p adj
+-1 -0.1 -4.05223799  3.852238 0.9999937
+-1  4.0  0.04776201  7.952238 0.0460196
+-1  6.4  2.44776201 10.352238 0.0003180
+-1  5.0  1.04776201  8.952238 0.0068354
+-2  4.1  0.14776201  8.052238 0.0385792
+-2  6.5  2.54776201 10.452238 0.0002524
+-2  5.1  1.14776201  9.052238 0.0055623
+-3  2.4 -1.55223799  6.352238 0.4291513
+-3  1.0 -2.95223799  4.952238 0.9510451
+-4 -1.4 -5.35223799  2.552238 0.8510119
 </code>
+====== eg. p.473 ======
+<code>> ns <- c(27, 34, 19, 20, 56, 35, 23, 37, 4, 30, 4, 42, 34, 19, 49)
+> ds <- c(48, 29, 34, 6, 18, 63, 9, 54, 28, 71, 60, 54, 51, 25, 49)
+> as <- c(34, 65, 55, 33, 42, 54, 21, 44, 61, 38, 75, 61, 51, 32, 47)
+>
+> smoke <- data.frame(ns,ds,as)
+> smoke
+   ns ds as
+  27 48 34
+  34 29 65
+  19 34 55
+  20  6 33
+  56 18 42
+  35 63 54
+  23  9 21
+  37 54 44
+   4 28 61
+30 71 38
+  4 60 75
+42 54 61
+34 51 51
+19 25 32
+49 49 47
+> smoke.st <- stack(smoke)
+> smoke.st
+   values ind
+      27  ns
+      34  ns
+      19  ns
+      20  ns
+      56  ns
+      35  ns
+      23  ns
+      37  ns
+       4  ns
+     30  ns
+      4  ns
+     42  ns
+     34  ns
+     19  ns
+     49  ns
+     48  ds
+     29  ds
+     34  ds
+      6  ds
+     18  ds
+     63  ds
+      9  ds
+     54  ds
+     28  ds
+     71  ds
+     60  ds
+     54  ds
+     51  ds
+     25  ds
+     49  ds
+     34  as
+     65  as
+     55  as
+     33  as
+     42  as
+     54  as
+     21  as
+     44  as
+     61  as
+     38  as
+     75  as
+     61  as
+     51  as
+     32  as
+     47  as
+> smoke.mod <- aov(values~ind,data=smoke.st)
+> smoke.mod
+Call:
+   aov(formula = values ~ ind, data = smoke.st)
+Terms:
+                      ind Residuals
+Sum of Squares   2643.378 11700.400
+Deg. of Freedom         2        42
+Residual standard error: 16.69074
+Estimated effects may be unbalanced
+> summary(smoke.mod)
+            Df Sum Sq Mean Sq F value Pr(>F)
+ind          2   2643  1321.7   4.744 0.0139 *
+Residuals   42  11700   278.6
+---
+Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
+> TukeyHSD(smoke.mod)
+  Tukey multiple comparisons of means
+% family-wise confidence level
+Fit: aov(formula = values ~ ind, data = smoke.st)
+$ind
+           diff       lwr       upr     p adj
+ds-as  -7.60000 -22.40679  7.206788 0.4327383
+ns-as -18.66667 -33.47345 -3.859878 0.0104480
+ns-ds -11.06667 -25.87345  3.740122 0.1768164
+> </code>