Differences

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

--- r:twoway_anova [2018/10/18 23:35] – [5] hkimscil
+++ r:twoway_anova [2018/10/18 23:39] (current) – [e.g. 5] hkimscil
@@ Line 218: / Line 218: @@
 [{{:r:interaction effect eg04.png|Interaction effect example 4}}]
-===== 5 =====
+===== e.g., 5 =====
 download {{:r:dataset_anova_twoWay_comparisons.csv|dataset_anova_twoWay_interactions.csv}}
@@ Line 268: / Line 268: @@
 </code>
-<code>> pairwise.t.test(d$StressReduction, d$Treatment, p.adj = "none")
+<code>> TukeyHSD(a.mod, which="Treatment")
+  Tukey multiple comparisons of means
+% family-wise confidence level
-	Pairwise comparisons using t tests with pooled SD
+Fit: aov(formula = StressReduction ~ Treatment * Age, data = stressdata)
-data:  d$StressReduction and d$Treatment
+$Treatment
+                 diff        lwr       upr     p adj
+mental-medical      2  0.7968988 3.2031012 0.0013531
+physical-medical    1 -0.2031012 2.2031012 0.1135025
+physical-mental    -1 -2.2031012 0.2031012 0.1135025
-         medical mental
+>
-mental   0.13    -
+</code>
-physical 0.45    0.45
+<code>
+> TukeyHSD(a.mod, which="Age")
+  Tukey multiple comparisons of means
+% family-wise confidence level
-P value adjustment method: none</code>
+Fit: aov(formula = StressReduction ~ Treatment * Age, data = stressdata)
-<code>> pairwise.t.test(d$StressReduction, d$Age, p.adj = "none")
-	Pairwise comparisons using t tests with pooled SD
+$Age
+          diff       lwr       upr    p adj
+old-mid     -3 -4.203101 -1.796899 1.54e-05
+young-mid    3  1.796899  4.203101 1.54e-05
+young-old    6  4.796899  7.203101 0.00e+00
+>
+</code>
-data:  d$StressReduction and d$Age
-      mid     old
-old   2.5e-05 -
-young 2.5e-05 2.3e-10
-P value adjustment method: none
->
-</code>
 <WRAP info>위는 아래의 linear model 을 이용하여도 가능. 사실 모든 ANOVA 테스트는 linear model이기도 함(lm)
@@ Line 309: / Line 314: @@
 >
-</code>
-<code>> pairwise.t.test(d$StressReduction, d$Treatment, p.adj = "none")
-	Pairwise comparisons using t tests with pooled SD
-data:  d$StressReduction and d$Treatment
-         medical mental
-mental   0.13    -
-physical 0.45    0.45
-P value adjustment method: none</code>
-<code>> pairwise.t.test(d$StressReduction, d$Age, p.adj = "none")
-	Pairwise comparisons using t tests with pooled SD
-data:  d$StressReduction and d$Age
-      mid     old
-old   2.5e-05 -
-young 2.5e-05 2.3e-10
-P value adjustment method: none
->
 </code>