anova_note:output01
Differences
This shows you the differences between two versions of the page.
| Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
| anova_note:output01 [2025/12/07 10:39] – hkimscil | anova_note:output01 [2025/12/09 22:54] (current) – hkimscil | ||
|---|---|---|---|
| Line 48: | Line 48: | ||
| [1] 0.009215657 | [1] 0.009215657 | ||
| > | > | ||
| + | </ | ||
| + | |||
| + | < | ||
| > # | > # | ||
| > comb <- list(o = o, p = p) | > comb <- list(o = o, p = p) | ||
| Line 70: | Line 73: | ||
| 5 102.94545 | 5 102.94545 | ||
| 6 103.89794 | 6 103.89794 | ||
| - | > boxplot(op$values~op$group) | + | > boxplot(op$values~ss.o <- ss(o) |
| + | ss.p <- ss(p) | ||
| + | df.o <- length(o)-1 | ||
| + | df.p <- length(p)-1 | ||
| + | |||
| + | |||
| + | op$group) | ||
| > | > | ||
| > | > | ||
| Line 81: | Line 90: | ||
| + c(" | + c(" | ||
| > | > | ||
| - | > # o = IQ | + | > |
| - | > # p = IQ | + | </code> |
| - | > # 둘을 합하여 구한 집합 op 의 IQ 점수 | + | {{.: |
| - | > # 이 전체집합의 분산은? | + | {{.: |
| + | |||
| + | <code> | ||
| > m.tot <- mean(op$values) | > m.tot <- mean(op$values) | ||
| > m.o <- mean(o) | > m.o <- mean(o) | ||
| Line 101: | Line 112: | ||
| > abline(v=m.tot, | > abline(v=m.tot, | ||
| > | > | ||
| - | > hist(o, breaks=br, | + | </code> |
| - | + col=rgb(1, | + | {{.: |
| - | > abline(v=m.o, col=" | + | |
| - | > abline(v=m.tot, | + | <code> |
| - | > | + | </code> |
| - | > hist(p, breaks=br, | + | |
| - | + col=rgb(.5, | + | <code> |
| - | > abline(v=m.p, | + | |
| - | > abline(v=m.tot, | + | |
| - | > | + | |
| > hist(o, breaks=br, | > hist(o, breaks=br, | ||
| + col=rgb(1, | + col=rgb(1, | ||
| Line 126: | Line 134: | ||
| [1] 5383.588 | [1] 5383.588 | ||
| > | > | ||
| - | > ss.o <- ss(o) | + | </code> |
| - | > ss.p <- ss(p) | + | {{.: |
| - | > df.o <- length(o)-1 | + | |
| - | > df.p <- length(p)-1 | + | <code> |
| - | > | + | |
| > m.tot | > m.tot | ||
| [1] 99.71228 | [1] 99.71228 | ||
| Line 145: | Line 152: | ||
| + col=rgb(1, | + col=rgb(1, | ||
| > abline(v=m.o, | > abline(v=m.o, | ||
| - | > hist(p, add=T, breaks=br, | + | > abline(v=m.tot, col=' |
| + | > | ||
| + | > hist(p, breaks=br, | ||
| + col=rgb(.5, | + col=rgb(.5, | ||
| > abline(v=m.p, | > abline(v=m.p, | ||
| > abline(v=m.tot, | > abline(v=m.tot, | ||
| - | > | + | > |
| - | > ss.bet <- n.o*(m.tot-m.o)^2 + n.p*(m.tot-m.p)^2 | + | > ss.bet <- n.o*(m.tot-m.o)^2 + # m.tot 에서 o그룹공통 까지의 거리를 제곱해서 모두 더한 값 |
| - | > ss.bet | + | # 아래 그림에서 빨간색 선에서 검은색 선까지의 거리를 제곱해서 모두 더한 값 |
| + | + n.p*(m.tot-m.p)^2 | ||
| + | # 아래 그림에서 빨간색 선에서 파란색 선까지의 거리를 제곱해서 모두 더한 값 | ||
| + | # 이것은 그룹 (IV, 독립변인) 때문에 생긴 그룹 간 차이이다 | ||
| + | > ss.bet | ||
| [1] 598.7747 | [1] 598.7747 | ||
| > | > | ||
| + | </ | ||
| + | {{.: | ||
| + | {{.: | ||
| + | |||
| + | < | ||
| > hist(o, breaks=br, | > hist(o, breaks=br, | ||
| + col=rgb(1, | + col=rgb(1, | ||
| > abline(v=m.o, | > abline(v=m.o, | ||
| - | > ss.o | + | > ss.o <- ss(o) # o집단의 평균인 검은색 선에서 개인 점수까지의 거리는 (오차는) 독립변인과 상관없이 랜덤하게 나타나는 것 |
| + | > ss.o # o집단의 것을 ss.o라고 부른다 | ||
| [1] 2179.19 | [1] 2179.19 | ||
| > | > | ||
| Line 163: | Line 182: | ||
| + col=rgb(.5, | + col=rgb(.5, | ||
| > abline(v=m.p, | > abline(v=m.p, | ||
| - | > ss.p | + | > ss.p <- ss(p) # p집단도 마찬가지이다. 이 집단 내의 sum of square값은 p 집단의 공통특징인 평균에서 개인점수가 랜덤하게 |
| + | > ss.p # 나타나는 것이고, 이것을 sum of square p라고 부른다 | ||
| [1] 2605.623 | [1] 2605.623 | ||
| > | > | ||
| + | > # 이 둘은 각 그룹의 평균을 중심으로 random 하게 나타나는 평균에서의 거리이다 (에러). | ||
| + | > # 따라서 우리는 이것을 sum of square within group이라고 부른다 | ||
| > ss.wit <- ss.o+ss.p | > ss.wit <- ss.o+ss.p | ||
| > ss.wit | > ss.wit | ||
| [1] 4784.813 | [1] 4784.813 | ||
| > | > | ||
| + | </ | ||
| + | {{.: | ||
| + | {{.: | ||
| + | |||
| + | |||
| + | < | ||
| > ss.bet | > ss.bet | ||
| [1] 598.7747 | [1] 598.7747 | ||
| Line 179: | Line 207: | ||
| > ss.tot | > ss.tot | ||
| [1] 5383.588 | [1] 5383.588 | ||
| - | > | + | > |
| + | </ | ||
| + | |||
| + | <code> | ||
| > df.tot <- length(op$values)-1 | > df.tot <- length(op$values)-1 | ||
| > df.bet <- nlevels(op$group) - 1 | > df.bet <- nlevels(op$group) - 1 | ||
| Line 190: | Line 221: | ||
| [1] 58 | [1] 58 | ||
| > | > | ||
| + | </ | ||
| + | |||
| + | < | ||
| > ms.tot <- ss.tot / df.tot | > ms.tot <- ss.tot / df.tot | ||
| > ms.bet <- ss.bet / df.bet | > ms.bet <- ss.bet / df.bet | ||
anova_note/output01.1765103951.txt.gz · Last modified: by hkimscil
