Differences

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

--- r:anova [2024/04/16 23:28] – [ANOVA in R: Output] hkimscil
+++ r:anova [2024/04/16 23:30] (current) – [ANOVA in R: Output] hkimscil
@@ Line 377: / Line 377: @@
 <code>
+> # 위에서
+> # ssd라는 function을 만들면
+>
+> ssd <- function(x) {
++     var(x) * (length(x)-1) }
+> ss.a1 <- ssd(A)
+> ss.b1 <- ssd(B)
+> ss.c1 <- ssd(C)
+>
+> ss.a == ss.a1
+[1] TRUE
+> ss.b == ss.b1
+[1] TRUE
+> ss.c == ss.c1
+[1] TRUE
+>
 </code>
 <code>
+> # 그러나 정확히 어떤 그룹에서 차이가 나는지는 판단해주지 않음
+> pairwise.t.test(comb3$values, comb3$group, p.adj = "none")
+	Pairwise comparisons using t tests with pooled SD
+data:  comb3$values and comb3$group
+  a      b
+b 0.4279 -
+c 0.0075 0.0516
+P value adjustment method: none
+> # OR
+> pairwise.t.test(comb3$values, comb3$group, p.adj = "bonf")
+	Pairwise comparisons using t tests with pooled SD
+data:  comb3$values and comb3$group
+  a     b
+b 1.000 -
+c 0.023 0.155
+P value adjustment method: bonferroni
+> pairwise.t.test(comb3$values, comb3$group, p.adj = "holm")
+	Pairwise comparisons using t tests with pooled SD
+data:  comb3$values and comb3$group
+  a     b
+b 0.428 -
+c 0.023 0.103
+P value adjustment method: holm
+>
+> # OR TukeyHSD(anova.output)
+> TukeyHSD(a.res)
+  Tukey multiple comparisons of means
+% family-wise confidence level
+Fit: aov(formula = values ~ group, data = comb3)
+$group
+    diff        lwr       upr     p adj
+b-a   -2  -8.059034  4.059034 0.7049466
+c-a   -7 -13.059034 -0.940966 0.0201250
+c-b   -5 -11.059034  1.059034 0.1238770
+>
 </code>
 ====== Post hoc test ======