Differences

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

--- r:anova [2023/06/05 10:27] – [ANOVA in R] hkimscil
+++ r:anova [2024/04/16 23:30] (current) – [ANOVA in R: Output] hkimscil
@@ Line 283: / Line 283: @@
 [1] 56.7234
 >
+> # within part 구하기
+> # 간단한 방법
+> tapply(comb3$values, comb3$group, var)*15
+  a   b   c
+750 900
+> sse <- sum(tapply(comb3$values, comb3$group, var)*15)
+> sse
+[1] 2250
+>
+> # 이해한 개념대로 얼른 구하기
 > # A, B, C 평균
 > m.a <- mean(A)
@@ Line 300: / Line 310: @@
 >
 > ss.within <- ss.a + ss.b + ss.c
-> ss.within <- ss.within
+> ss.within == sse # check both are the same
+[1] TRUE
 >
+>
+> # between part 구하기
 > 16*((m.a-mean.tot)^2)
 [1] 144
@@ Line 359: / Line 372: @@
 ---
 Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
+>
+>
+</code>
+<code>
 > # 위에서
 > # ssd라는 function을 만들면
@@ Line 365: / Line 383: @@
 +     var(x) * (length(x)-1) }
 > ss.a1 <- ssd(A)
-> ss.b2 <- ssd(B)
+> ss.b1 <- ssd(B)
 > ss.c1 <- ssd(C)
 >
@@ Line 374: / Line 392: @@
 > ss.c == ss.c1
 [1] TRUE
+>
+</code>
+<code>
 > # 그러나 정확히 어떤 그룹에서 차이가 나는지는 판단해주지 않음
 > pairwise.t.test(comb3$values, comb3$group, p.adj = "none")
@@ Line 423: / Line 445: @@
 c-b   -5 -11.059034  1.059034 0.1238770
 >
->
->
 </code>
 ====== Post hoc test ======
 [[:post hoc test]]