anova
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revisionNext revisionBoth sides next revision | ||
anova [2020/05/20 15:36] – [Example] hkimscil | anova [2020/05/20 15:51] – [Example 2] hkimscil | ||
---|---|---|---|
Line 580: | Line 580: | ||
F< | F< | ||
====== Example 2 ====== | ====== Example 2 ====== | ||
- | < | + | < |
+ | x1 <- c(15, 20, 14, 13, 18, 16, 13, 12, 18, 41) | ||
x2 <- c(21, 25, 29, 18, 26, 22, 26, 24, 28, 21) | x2 <- c(21, 25, 29, 18, 26, 22, 26, 24, 28, 21) | ||
x3 <- c(28, 30, 32, 28, 26, 30, 25, 36, 20, 15) | x3 <- c(28, 30, 32, 28, 26, 30, 25, 36, 20, 15) | ||
</ | </ | ||
- | < | + | < |
x1 x2 x3 | x1 x2 x3 | ||
1 15 21 28 | 1 15 21 28 | ||
Line 634: | Line 635: | ||
< | < | ||
- | > colnames(xs[1]) <- " | + | > colnames(xs) |
- | > colnames(xs[2]) <- "cond" | + | |
</ | </ | ||
- | < | + | < |
+ | # cf | ||
# lengthofelements <- length(x1) | # lengthofelements <- length(x1) | ||
# varofvariable <- var(x1)</ | # varofvariable <- var(x1)</ | ||
< | < | ||
- | df_x1 | + | df.total <- length(xs$wrong) - 1 |
- | df_x2 | + | ss.total <- var(xs$wrong)*df_tot |
- | df_x3 | + | var.total <- ss.total/ |
+ | var.total.r <- var(xs$wrong) | ||
- | ss_x1 | + | df.x1 <- length(x1)-1 |
- | ss_x2 | + | df.x2 <- length(x2)-1 |
- | ss_x3 | + | df.x3 <- length(x3)-1 |
+ | ss.x1 <- var(x1)*df.x1 | ||
+ | ss.x2 <- var(x2)*df.x2 | ||
+ | ss.x3 <- var(x3)*df.x3 | ||
- | df_bet | + | ss.within <- ss.x1 + ss.x2 + ss.x3 |
- | ss_bet | + | df.within <- df.x1 + df.x2 + df.x3 |
+ | ss.between <- ss.total - ss.within | ||
+ | df.between <- df.total - df.within | ||
- | df_tot | + | ms.between <- ss.between/ |
- | ss_tot | + | ms.within <- ss.within/ |
+ | f.value <- ms.between/ | ||
- | df_with | + | ss.between |
- | ss_with | + | df.between |
- | df_bet | + | ss.within |
- | ss_bet | + | df.within |
- | </ | + | ms.between |
+ | ms.within | ||
+ | |||
+ | f.value | ||
+ | |||
- | < | ||
- | df_tot <- length(xs$ind) - 1 | ||
- | ss_tot <- var(xs$values)*df_tot | ||
- | var_tot <- var(xs$values) | ||
- | df_x1 <- length(x1)-1 | ||
- | df_x2 <- length(x2)-1 | ||
- | df_x3 <- length(x3)-1 | ||
- | ss_x1 <- var(x1)*df_x1 | ||
- | ss_x2 <- var(x2)*df_x2 | ||
- | ss_x3 <- var(x3)*df_x3 | ||
</ | </ | ||
====== E.G. 1 (R) ====== | ====== E.G. 1 (R) ====== |
anova.txt · Last modified: 2022/09/30 09:02 by hkimscil