User Tools

Site Tools


r:anova

Differences

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

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision
Previous revision
r:anova [2023/06/05 19:19] – [Post hoc test] hkimscilr:anova [2024/04/17 08:30] (current) – [ANOVA in R: Output] hkimscil
Line 42: Line 42:
 var(comb3$values) var(comb3$values)
  
 +# within part 구하기
 +# 간단한 방법
 +tapply(comb3$values, comb3$group, var)*15
 +sse <- sum(tapply(comb3$values, comb3$group, var)*15)
 +sse
 +
 +# 이해한 개념대로 얼른 구하기
 # A, B, C 평균 # A, B, C 평균
 m.a <- mean(A) m.a <- mean(A)
Line 59: Line 66:
  
 ss.within <- ss.a + ss.b + ss.c ss.within <- ss.a + ss.b + ss.c
-ss.within <- ss.within+ss.within == sse # check both are the same 
  
 +# between part 구하기
 16*((m.a-mean.tot)^2) 16*((m.a-mean.tot)^2)
 16*((m.b-mean.tot)^2) 16*((m.b-mean.tot)^2)
Line 274: Line 283:
 [1] 56.7234 [1] 56.7234
  
 +> # within part 구하기
 +> # 간단한 방법
 +> tapply(comb3$values, comb3$group, var)*15
 +  a     
 +600 750 900 
 +> sse <- sum(tapply(comb3$values, comb3$group, var)*15)
 +> sse
 +[1] 2250
 +
 +> # 이해한 개념대로 얼른 구하기
 > # A, B, C 평균 > # A, B, C 평균
 > m.a <- mean(A) > m.a <- mean(A)
Line 291: Line 310:
  
 > ss.within <- ss.a + ss.b + ss.c > 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) > 16*((m.a-mean.tot)^2)
 [1] 144 [1] 144
Line 350: Line 372:
 --- ---
 Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
 +
 +
 +</code>
 +
 +<code>
 > # 위에서  > # 위에서 
 > # ssd라는 function을 만들면  > # ssd라는 function을 만들면 
Line 356: Line 383:
 +     var(x) * (length(x)-1) } +     var(x) * (length(x)-1) }
 > ss.a1 <- ssd(A) > ss.a1 <- ssd(A)
-> ss.b2 <- ssd(B)+> ss.b1 <- ssd(B)
 > ss.c1 <- ssd(C) > ss.c1 <- ssd(C)
  
Line 365: Line 392:
 > ss.c == ss.c1 > ss.c == ss.c1
 [1] TRUE [1] TRUE
 +
 +</code>
 +
 +<code>
 > # 그러나 정확히 어떤 그룹에서 차이가 나는지는 판단해주지 않음  > # 그러나 정확히 어떤 그룹에서 차이가 나는지는 판단해주지 않음 
 > pairwise.t.test(comb3$values, comb3$group, p.adj = "none") > pairwise.t.test(comb3$values, comb3$group, p.adj = "none")
Line 414: Line 445:
 c-b   -5 -11.059034  1.059034 0.1238770 c-b   -5 -11.059034  1.059034 0.1238770
  
->  +>
->  +
-+
 </code> </code>
- 
 ====== Post hoc test ====== ====== Post hoc test ======
 [[:post hoc test]] [[:post hoc test]]
r/anova.1685960381.txt.gz · Last modified: 2023/06/05 19:19 by hkimscil

Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki