Expected mean squares

In statistics, expected mean squares (EMS) are the expected values of certain statistics arising in partitions of sums of squares in the analysis of variance (ANOVA). They can be used for ascertaining which statistic should appear in the denominator in an F-test for testing a null hypothesis that a particular effect is absent.

Definition

When the total corrected sum of squares in an ANOVA is partitioned into several components, each attributed to the effect of a particular predictor variable, each of the sums of squares in that partition is a random variable that has an expected value. That expected value divided by the corresponding number of degrees of freedom is the expected mean square for that predictor variable.

Example

The following example is from Longitudinal Data Analysis by Donald Hedeker and Robert D. Gibbons.

Each of s treatments (one of which may be a placebo) is administered to a sample of (capital) N randomly chosen patients, on whom certain measurements [math]\displaystyle{ Y_{hij} }[/math] are observed at each of (lower-case) n specified times, for [math]\displaystyle{ h=1,\ldots,s, \quad i=1,\ldots,N_h }[/math] (thus the numbers of patients receiving different treatments may differ), and [math]\displaystyle{ j=1,\ldots, n. }[/math] We assume the sets of patients receiving different treatments are disjoint, so patients are nested within treatments and not crossed with treatments. We have

[math]\displaystyle{ Y_{hij} = \mu + \gamma_h + \tau_j + (\gamma\tau)_{hj} + \pi_{i(h)} + \varepsilon_{hij} }[/math]

where

• [math]\displaystyle{ \mu }[/math] = grand mean, (fixed)
• [math]\displaystyle{ \gamma_h }[/math] = effect of treatment [math]\displaystyle{ h }[/math], (fixed)
• [math]\displaystyle{ \tau_j }[/math] = effect of time [math]\displaystyle{ j }[/math], (fixed)
• [math]\displaystyle{ (\gamma\tau)_{hj} }[/math] = interaction effect of treatment [math]\displaystyle{ h }[/math] and time [math]\displaystyle{ j }[/math], (fixed)
• [math]\displaystyle{ \pi_{i(h)} }[/math] = individual difference effect for patient [math]\displaystyle{ i }[/math] nested within treatment [math]\displaystyle{ h }[/math], (random)
• [math]\displaystyle{ \varepsilon_{hij} }[/math] = error for patient [math]\displaystyle{ i }[/math] in treatment [math]\displaystyle{ h }[/math] at time [math]\displaystyle{ j }[/math]. (random)
• [math]\displaystyle{ \sigma_\pi^2 }[/math] = variance of the random effect of patients nested within treatments,
• [math]\displaystyle{ \sigma_\varepsilon }[/math] = error variance.

The total corrected sum of squares is

[math]\displaystyle{ \sum_{hij} (Y_{hij} - \overline Y)^2 \quad\text{where } \overline Y = \frac 1 n \sum_{hij} Y_{hij}. }[/math]

The ANOVA table below partitions the sum of squares (where [math]\displaystyle{ N = \sum_h N_h }[/math]):

source of variability	degrees of freedom	sum of squares	mean square	expected mean square
treatment	[math]\displaystyle{ s - 1 }[/math]	[math]\displaystyle{ \text{SS}_\text{Tr} = n\sum_{h=1}^s N_h(\overline Y_{h\cdot\cdot} - \overline Y_{\cdot\cdot\cdot})^2 }[/math]	[math]\displaystyle{ \dfrac{\text{SS}_\text{Tr}}{s-1} }[/math]	[math]\displaystyle{ \sigma_\varepsilon^2 + n \sigma_\pi^2 + D_\text{Tr} }[/math]
time	[math]\displaystyle{ n - 1 }[/math]	[math]\displaystyle{ \text{SS}_\text{T} = N \sum_{j=1}^n (\overline Y_{\cdot\cdot j} - \overline Y_{\cdot\cdot\cdot})^2 }[/math]	[math]\displaystyle{ \dfrac{\text{SS}_\text{T}}{n-1} }[/math]	[math]\displaystyle{ \sigma_\varepsilon^2 + D_\text{T} }[/math]
treatment × time	[math]\displaystyle{ (s-1)(n-1) }[/math]	[math]\displaystyle{ \text{SS}_\text{Tr T} = \sum_{h=1}^s \sum_{j=1}^n N_h (\overline Y_{h\cdot j} - \overline Y_{h\cdot\cdot} - \overline Y_{\cdot\cdot j} + \overline Y_{\cdot\cdot\cdot})^2 }[/math]	[math]\displaystyle{ \dfrac{\text{SS}_\text{Tr T}}{(n-1)(s-1)} }[/math]	[math]\displaystyle{ \sigma_\varepsilon^2 + D_\text{Tr T} }[/math]
patients within treatments	[math]\displaystyle{ N - s }[/math]	[math]\displaystyle{ \text{SS}_{\text{S}(\text{Tr})} = n\sum_{h=1}^s \sum_{i=1}^{N_h} (\overline Y_{hi\cdot} - \overline Y_{h\cdot\cdot})^2 }[/math]	[math]\displaystyle{ \dfrac{\text{SS}_{\text{S}(\text{Tr})}}{N-s} }[/math]	[math]\displaystyle{ \sigma_\varepsilon^2 + n\sigma_\pi^2 }[/math]
error	[math]\displaystyle{ (N - s)(n - 1) }[/math]	[math]\displaystyle{ \text{SS}_\text{E} = \sum_{h=1}^s \sum_{i=1}^{N_h} \sum_{j=1}^n (Y_{hij} - \overline Y_{h\cdot j} - \overline Y_{hi\cdot} + \overline Y_{h\cdot\cdot})^2 }[/math]	[math]\displaystyle{ \dfrac{\text{SS}_\text{E}}{(N-s)(n-1)} }[/math]	[math]\displaystyle{ \sigma_\varepsilon^2 }[/math]

Use in F-tests

A null hypothesis of interest is that there is no difference between effects of different treatments—thus no difference among treatment means. This may be expressed by saying [math]\displaystyle{ D_\text{Tr}=0, }[/math] (with the notation as used in the table above). Under this null hypothesis, the expected mean square for effects of treatments is [math]\displaystyle{ \sigma_\varepsilon^2 + n \sigma_\pi^2. }[/math]

The numerator in the F-statistic for testing this hypothesis is the mean square due to differences among treatments, i.e. it is [math]\displaystyle{ \left. \text{SS}_\text{Tr} \right/(s-1). }[/math] The denominator, however, is not [math]\displaystyle{ \left. \text{SS}_\text{E}\right/ \big( (N-s)(n-1) \big). }[/math] The reason is that the random variable below, although under the null hypothesis it has an F-distribution, is not observable—it is not a statistic—because its value depends on the unobservable parameters [math]\displaystyle{ \sigma_\pi^2 }[/math] and [math]\displaystyle{ \sigma_\varepsilon^2. }[/math]

[math]\displaystyle{ \frac{\left.\frac{\text{SS}_\text{Tr}}{\sigma_\varepsilon^2 + n\sigma_\pi^2} \right/(s-1)}{ \left. \frac{\text{SS}_\text{E}}{\sigma_\varepsilon^2} \right/ \big( (N-s)(n-1) \big)} \ne \frac{\text{SS}_\text{Tr}/(s-1)}{\text{SS}_\text{E}/\big((N-s)(n-1)\big)} }[/math]

Instead, one uses as the test statistic the following random variable that is not defined in terms of [math]\displaystyle{ \text{SS}_\text{E} }[/math]:

[math]\displaystyle{ F = \frac{\left.\frac{\text{SS}_\text{Tr}}{\sigma_\varepsilon^2 + n\sigma_\pi^2} \right/(s-1)}{ \left. \frac{\text{SS}_{\text{S}(\text{Tr})}}{\sigma_\varepsilon^2+ n\sigma_\pi^2} \right/ (N-s)} = \frac{\left. \text{SS}_\text{Tr} \right/(s-1)}{ \left. \text{SS}_\text{S(Tr)} \right/ (N-s)} }[/math]