c:ms:2023:w07_anova_note
Differences
This shows you the differences between two versions of the page.
| Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
| c:ms:2023:w07_anova_note [2023/04/18 13:30] – hkimscil | c:ms:2023:w07_anova_note [2023/06/04 12:38] (current) – [Post hoc test] hkimscil | ||
|---|---|---|---|
| Line 125: | Line 125: | ||
| </ | </ | ||
| - | ====== Output | + | |
| + | ===== ANOVA in R: Output ===== | ||
| < | < | ||
| > # | > # | ||
| Line 433: | Line 434: | ||
| d.bc | d.bc | ||
| - | # se | + | # mse (ms within) |
| - | # from the a.res.sum output | + | # a.res.sum == summary(aov(values ~ group, data=comb3)) |
| a.res.sum | a.res.sum | ||
| # mse = 50 | # mse = 50 | ||
| mse <- 50 | mse <- 50 | ||
| + | # 혹은 fansy way from comb3 data.frame | ||
| + | # 15 는 각 그룹의 df | ||
| + | sse.ch <- sum(tapply(comb3$values, | ||
| + | sse.ch | ||
| + | mse.ch <- sse.ch/45 | ||
| + | mse.ch | ||
| se <- sqrt(mse/ | se <- sqrt(mse/ | ||
| Line 467: | Line 474: | ||
| </ | </ | ||
| + | < | ||
| + | plot(TukeyHSD(a.res, | ||
| + | </ | ||
| + | |||
| + | < | ||
| + | pairwise.t.test(comb3$values, | ||
| + | </ | ||
| ===== post hoc test: output ===== | ===== post hoc test: output ===== | ||
| - | < | + | < |
| > m.a | > m.a | ||
| [1] 26 | [1] 26 | ||
| Line 544: | Line 558: | ||
| </ | </ | ||
| + | {{: | ||
| ====== R square or Eta square ====== | ====== R square or Eta square ====== | ||
| SS toal | SS toal | ||
c/ms/2023/w07_anova_note.1681824629.txt.gz · Last modified: by hkimscil
