Two way ANOVA - overview
This page offers structured overviews of one or more selected methods. Add additional methods for comparisons (max. of 3) by clicking on the dropdown button in the right-hand column. To practice with a specific method click the button at the bottom row of the table
Two way ANOVA | Friedman test | Goodness of fit test | One sample Wilcoxon signed-rank test | Paired sample $t$ test | $z$ test for a single proportion |
|
---|---|---|---|---|---|---|
Independent/grouping variables | Independent/grouping variable | Independent variable | Independent variable | Independent variable | Independent variable | |
Two categorical, the first with $I$ independent groups and the second with $J$ independent groups ($I \geqslant 2$, $J \geqslant 2$) | One within subject factor ($\geq 2$ related groups) | None | None | 2 paired groups | None | |
Dependent variable | Dependent variable | Dependent variable | Dependent variable | Dependent variable | Dependent variable | |
One quantitative of interval or ratio level | One of ordinal level | One categorical with $J$ independent groups ($J \geqslant 2$) | One of ordinal level | One quantitative of interval or ratio level | One categorical with 2 independent groups | |
Null hypothesis | Null hypothesis | Null hypothesis | Null hypothesis | Null hypothesis | Null hypothesis | |
ANOVA $F$ tests:
| H0: the population scores in any of the related groups are not systematically higher or lower than the population scores in any of the other related groups
Usually the related groups are the different measurement points. Several different formulations of the null hypothesis can be found in the literature, and we do not agree with all of them. Make sure you (also) learn the one that is given in your text book or by your teacher. |
| H0: $m = m_0$
Here $m$ is the population median, and $m_0$ is the population median according to the null hypothesis. | H0: $\mu = \mu_0$
Here $\mu$ is the population mean of the difference scores, and $\mu_0$ is the population mean of the difference scores according to the null hypothesis, which is usually 0. A difference score is the difference between the first score of a pair and the second score of a pair. | H0: $\pi = \pi_0$
Here $\pi$ is the population proportion of 'successes', and $\pi_0$ is the population proportion of successes according to the null hypothesis. | |
Alternative hypothesis | Alternative hypothesis | Alternative hypothesis | Alternative hypothesis | Alternative hypothesis | Alternative hypothesis | |
ANOVA $F$ tests:
| H1: the population scores in some of the related groups are systematically higher or lower than the population scores in other related groups |
| H1 two sided: $m \neq m_0$ H1 right sided: $m > m_0$ H1 left sided: $m < m_0$ | H1 two sided: $\mu \neq \mu_0$ H1 right sided: $\mu > \mu_0$ H1 left sided: $\mu < \mu_0$ | H1 two sided: $\pi \neq \pi_0$ H1 right sided: $\pi > \pi_0$ H1 left sided: $\pi < \pi_0$ | |
Assumptions | Assumptions | Assumptions | Assumptions | Assumptions | Assumptions | |
|
|
|
|
|
| |
Test statistic | Test statistic | Test statistic | Test statistic | Test statistic | Test statistic | |
For main and interaction effects together (model):
| $Q = \dfrac{12}{N \times k(k + 1)} \sum R^2_i - 3 \times N(k + 1)$
Here $N$ is the number of 'blocks' (usually the subjects - so if you have 4 repeated measurements for 60 subjects, $N$ equals 60), $k$ is the number of related groups (usually the number of repeated measurements), and $R_i$ is the sum of ranks in group $i$. Remember that multiplication precedes addition, so first compute $\frac{12}{N \times k(k + 1)} \times \sum R^2_i$ and then subtract $3 \times N(k + 1)$. Note: if ties are present in the data, the formula for $Q$ is more complicated. | $X^2 = \sum{\frac{(\mbox{observed cell count} - \mbox{expected cell count})^2}{\mbox{expected cell count}}}$
Here the expected cell count for one cell = $N \times \pi_j$, the observed cell count is the observed sample count in that same cell, and the sum is over all $J$ cells. | Two different types of test statistics can be used, but both will result in the same test outcome. We will denote the first option the $W_1$ statistic (also known as the $T$ statistic), and the second option the $W_2$ statistic.
In order to compute each of the test statistics, follow the steps below:
| $t = \dfrac{\bar{y} - \mu_0}{s / \sqrt{N}}$
Here $\bar{y}$ is the sample mean of the difference scores, $\mu_0$ is the population mean of the difference scores according to the null hypothesis, $s$ is the sample standard deviation of the difference scores, and $N$ is the sample size (number of difference scores). The denominator $s / \sqrt{N}$ is the standard error of the sampling distribution of $\bar{y}$. The $t$ value indicates how many standard errors $\bar{y}$ is removed from $\mu_0$. | $z = \dfrac{p - \pi_0}{\sqrt{\dfrac{\pi_0(1 - \pi_0)}{N}}}$
Here $p$ is the sample proportion of successes: $\dfrac{X}{N}$, $N$ is the sample size, and $\pi_0$ is the population proportion of successes according to the null hypothesis. | |
Pooled standard deviation | n.a. | n.a. | n.a. | n.a. | n.a. | |
$ \begin{aligned} s_p &= \sqrt{\dfrac{\sum\nolimits_{subjects} (\mbox{subject's score} - \mbox{its group mean})^2}{N - (I \times J)}}\\ &= \sqrt{\dfrac{\mbox{sum of squares error}}{\mbox{degrees of freedom error}}}\\ &= \sqrt{\mbox{mean square error}} \end{aligned} $ | - | - | - | - | - | |
Sampling distribution of $F$ if H0 were true | Sampling distribution of $Q$ if H0 were true | Sampling distribution of $X^2$ if H0 were true | Sampling distribution of $W_1$ and of $W_2$ if H0 were true | Sampling distribution of $t$ if H0 were true | Sampling distribution of $z$ if H0 were true | |
For main and interaction effects together (model):
| If the number of blocks $N$ is large, approximately the chi-squared distribution with $k - 1$ degrees of freedom.
For small samples, the exact distribution of $Q$ should be used. | Approximately the chi-squared distribution with $J - 1$ degrees of freedom | Sampling distribution of $W_1$:
If $N_r$ is large, $W_1$ is approximately normally distributed with mean $\mu_{W_1}$ and standard deviation $\sigma_{W_1}$ if the null hypothesis were true. Here $$\mu_{W_1} = \frac{N_r(N_r + 1)}{4}$$ $$\sigma_{W_1} = \sqrt{\frac{N_r(N_r + 1)(2N_r + 1)}{24}}$$ Hence, if $N_r$ is large, the standardized test statistic $$z = \frac{W_1 - \mu_{W_1}}{\sigma_{W_1}}$$ follows approximately the standard normal distribution if the null hypothesis were true. Sampling distribution of $W_2$: If $N_r$ is large, $W_2$ is approximately normally distributed with mean $0$ and standard deviation $\sigma_{W_2}$ if the null hypothesis were true. Here $$\sigma_{W_2} = \sqrt{\frac{N_r(N_r + 1)(2N_r + 1)}{6}}$$ Hence, if $N_r$ is large, the standardized test statistic $$z = \frac{W_2}{\sigma_{W_2}}$$ follows approximately the standard normal distribution if the null hypothesis were true. If $N_r$ is small, the exact distribution of $W_1$ or $W_2$ should be used. Note: if ties are present in the data, the formula for the standard deviations $\sigma_{W_1}$ and $\sigma_{W_2}$ is more complicated. | $t$ distribution with $N - 1$ degrees of freedom | Approximately the standard normal distribution | |
Significant? | Significant? | Significant? | Significant? | Significant? | Significant? | |
| If the number of blocks $N$ is large, the table with critical $X^2$ values can be used. If we denote $X^2 = Q$:
|
| For large samples, the table for standard normal probabilities can be used: Two sided:
| Two sided:
| Two sided:
| |
n.a. | n.a. | n.a. | n.a. | $C\%$ confidence interval for $\mu$ | Approximate $C\%$ confidence interval for $\pi$ | |
- | - | - | - | $\bar{y} \pm t^* \times \dfrac{s}{\sqrt{N}}$
where the critical value $t^*$ is the value under the $t_{N-1}$ distribution with the area $C / 100$ between $-t^*$ and $t^*$ (e.g. $t^*$ = 2.086 for a 95% confidence interval when df = 20). The confidence interval for $\mu$ can also be used as significance test. | Regular (large sample):
| |
Effect size | n.a. | n.a. | n.a. | Effect size | n.a. | |
| - | - | - | Cohen's $d$: Standardized difference between the sample mean of the difference scores and $\mu_0$: $$d = \frac{\bar{y} - \mu_0}{s}$$ Cohen's $d$ indicates how many standard deviations $s$ the sample mean of the difference scores $\bar{y}$ is removed from $\mu_0.$ | - | |
n.a. | n.a. | n.a. | n.a. | Visual representation | n.a. | |
- | - | - | - | - | ||
ANOVA table | n.a. | n.a. | n.a. | n.a. | n.a. | |
- | - | - | - | - | ||
Equivalent to | n.a. | n.a. | n.a. | Equivalent to | Equivalent to | |
OLS regression with two categorical independent variables and the interaction term, transformed into $(I - 1)$ + $(J - 1)$ + $(I - 1) \times (J - 1)$ code variables. | - | - | - |
|
| |
Example context | Example context | Example context | Example context | Example context | Example context | |
Is the average mental health score different between people from a low, moderate, and high economic class? And is the average mental health score different between men and women? And is there an interaction effect between economic class and gender? | Is there a difference in depression level between measurement point 1 (pre-intervention), measurement point 2 (1 week post-intervention), and measurement point 3 (6 weeks post-intervention)? | Is the proportion of people with a low, moderate, and high social economic status in the population different from $\pi_{low} = 0.2,$ $\pi_{moderate} = 0.6,$ and $\pi_{high} = 0.2$? | Is the median mental health score of office workers different from $m_0 = 50$? | Is the average difference between the mental health scores before and after an intervention different from $\mu_0 = 0$? | Is the proportion of smokers amongst office workers different from $\pi_0 = 0.2$? Use the normal approximation for the sampling distribution of the test statistic. | |
SPSS | SPSS | SPSS | SPSS | SPSS | SPSS | |
Analyze > General Linear Model > Univariate...
| Analyze > Nonparametric Tests > Legacy Dialogs > K Related Samples...
| Analyze > Nonparametric Tests > Legacy Dialogs > Chi-square...
| Specify the measurement level of your variable on the Variable View tab, in the column named Measure. Then go to:
Analyze > Nonparametric Tests > One Sample...
| Analyze > Compare Means > Paired-Samples T Test...
| Analyze > Nonparametric Tests > Legacy Dialogs > Binomial...
| |
Jamovi | Jamovi | Jamovi | Jamovi | Jamovi | Jamovi | |
ANOVA > ANOVA
| ANOVA > Repeated Measures ANOVA - Friedman
| Frequencies > N Outcomes - $\chi^2$ Goodness of fit
| T-Tests > One Sample T-Test
| T-Tests > Paired Samples T-Test
| Frequencies > 2 Outcomes - Binomial test
| |
Practice questions | Practice questions | Practice questions | Practice questions | Practice questions | Practice questions | |