Tutorials › ANOVA › Choosing Between Classical and Welch One-Way Tests

Heteroscedastic ANOVA · Tutorial 559 of 595

Choosing Between Classical and Welch One-Way Tests

Use the study design and groupwise variability to choose a defensible one-way test, then interpret its result without relying on an automatic variance-pretest rule.

Intermediate 12 min read

What You'll Learn

  • State the shared null hypothesis and assumptions for classical and Welch one-way tests
  • Compare group standard deviations and sample sizes when assessing variance heterogeneity
  • Explain why balanced allocation can reduce, but does not eliminate, the effects of unequal variances
  • Choose a test using design knowledge and the inferential goal rather than a variance pretest alone
  • Work through examples where the balance and variability patterns lead to different choices
  • Describe how to report the choice and its limitations

Why the Choice Depends on More Than a Variance Test

A one-way analysis of variance asks whether the population means of several groups are equal. The classical test estimates one common within-group variance by pooling the group variances. Welch’s test instead allows the group variances to differ and adjusts the test statistic and its reference degrees of freedom. The practical choice depends on the plausibility of a common variance, the pattern of group sample sizes, and the study’s inferential goal.

A variance pretest may seem like a simple decision rule: use the classical test if the pretest is not significant, and Welch’s test otherwise. That rule is unreliable. A nonsignificant pretest does not establish equal variances, particularly when samples are small; and choosing the final procedure based on a preliminary test can affect the properties of the overall testing process. The choice should instead be guided by the design, groupwise variability summaries, and reasonable knowledge of how the data were generated.

The Two Tests and Their Assumptions

Definition: For \(k\) independent groups, let \(n_i\), \(\bar{x}_i\), and \(s_i\) be the sample size, sample mean, and sample standard deviation in group \(i\), for \(i=1,\ldots,k\). The one-way equal-means null hypothesis is \(H_0:\mu_1=\cdots=\mu_k\); the alternative is that at least one population mean differs.

Both tests require independent observations within and across groups. For the classical one-way \(F\) test to have its exact reference distribution, observations within each group are normally distributed and all groups share a common population variance, \(\sigma^2\). Welch’s test does not require a common variance. It uses an approximate \(F\) reference distribution, so the within-group distributions should still be reasonably well behaved, especially in small samples: severe skewness or influential outliers can undermine either method.

The classical test pools the within-group variation. Define \(N=\sum_i n_i\), the grand mean \(\bar{x}_{\cdot\cdot}=\sum_i n_i\bar{x}_i/N\), the between-group sum of squares \(SSB=\sum_i n_i(\bar{x}_i-\bar{x}_{\cdot\cdot})^2\), and the within-group sum of squares \(SSW=\sum_i(n_i-1)s_i^2\). Its statistic is \(F=MSB/MSW\), where \(MSB=SSB/(k-1)\) and \(MSW=SSW/(N-k)\); under \(H_0\) and the classical assumptions, the reference degrees of freedom are \(k-1\) and \(N-k\).

Welch’s test gives more weight to group means estimated with greater precision. Its weights are \(w_i=n_i/s_i^2\), with \(W=\sum_i w_i\) and weighted mean \(\bar{x}_w=\sum_iw_i\bar{x}_i/W\). Define \(A=\sum_i(1-w_i/W)^2/(n_i-1)\). The Welch statistic and its approximate denominator degrees of freedom are:

$$ F_W = \frac{\displaystyle \sum_{i=1}^{k} w_i(\bar{x}_i-\bar{x}_w)^2/(k-1)} {\displaystyle 1+\frac{2(k-2)}{k^2-1}A}, \qquad df_1=k-1,\qquad df_2=\frac{k^2-1}{3A}. $$

The Welch \(p\)-value is calculated from an \(F\) distribution with these approximate degrees of freedom. Both procedures test equality of population means; choosing Welch does not change the scientific question. It changes how the uncertainty in group means is estimated and calibrated.

How Balance Changes the Assessment

Groupwise variability summaries should be read alongside sample sizes. A useful first check is to compare the \(s_i\) values (or their squares), not just a single pooled summary. The range of standard deviations can flag potential heterogeneity, but there is no universal ratio that automatically determines the test. Consider whether the differences are substantial in context and whether they are supported by the design or prior knowledge.

Balance matters because the classical test pools variances while the precision of each group mean depends on that group’s own variance and sample size. A short expectation calculation explains why equal sample sizes can make the classical test less sensitive to unequal variances, without making the equal-variance assumption true.

Justification: Suppose group \(i\) has population variance \(\sigma_i^2\). Under independent sampling, \(E(MSW)=\sum_i(n_i-1)\sigma_i^2/(N-k)\). If every group has the same size \(m\), this becomes \(\sum_i\sigma_i^2/k\), the average group variance. Also, the variance of a difference between two sample means is \(\operatorname{Var}(\bar{X}_i-\bar{X}_j)=\sigma_i^2/m+\sigma_j^2/m\). Thus equal allocation avoids pairing a much larger group with a much smaller group in the pooled estimate and mean comparison. It does not make the variances equal or guarantee the exact classical \(F\) calibration under heterogeneity.

With unequal allocation, the relationship between sample size and variance is particularly important. A small group with a large variance, or a large group with a small variance, can make the pooled estimate a poor match to the uncertainty in the group means. The earlier material on unequal allocation explains why these patterns can affect the classical test’s error rate. Welch is generally a safer choice when meaningful heterogeneity is plausible, particularly alongside unequal sample sizes.

Worked Examples: Applying the Choice

Worked Example: Balanced Groups With Noticeably Different Spreads

Suppose three hypothetical independent groups have \(n_1=n_2=n_3=10\), means \(10,12,14\) units, and standard deviations \(2,5,8\) units. Assume the study design supports independence and that plots show no severe skewness or influential outliers. The equal sample sizes are helpful, but the standard deviations range from 2 to 8 units, so a common variance is not especially convincing. We calculate both tests to see how the choice plays out.

For the classical test, \(\bar{x}_{\cdot\cdot}=(10+12+14)/3=12\) units because the group sizes are equal. Thus \(SSB=10[(10-12)^2+(12-12)^2+(14-12)^2]=80\), and \(MSB=80/2=40\). The pooled within-group sum of squares is \(SSW=9(2^2+5^2+8^2)=9(93)=837\); hence \(MSW=837/27=31.00\). The statistic is \(F=40/31=1.290\), with \(df=(2,27)\), giving \(p\approx0.292\).

For Welch’s test, the weights are \(w_1=10/4=2.500\), \(w_2=10/25=0.4000\), and \(w_3=10/64=0.1563\), so \(W=3.0563\). The weighted mean is \(\bar{x}_w=[2.5(10)+0.4(12)+0.1563(14)]/3.0563\approx10.47\) units. Substitution gives \(\sum_iw_i(\bar{x}_i-\bar{x}_w)^2\approx3.44\). The correction term is \(A\approx(1-2.5/3.0563)^2/9+(1-0.4/3.0563)^2/9+(1-0.1563/3.0563)^2/9\approx0.1876\). Therefore \(F_W=(3.44/2)/[1+0.25(0.1876)]\approx1.640\), with \(df_1=2\) and \(df_2=8/[3(0.1876)]\approx14.21\). The approximate \(p\)-value is \(0.229\).

Neither test provides strong evidence against equal means at \(\alpha=0.05\). That does not settle which method is preferable: the noticeable spread differences make Welch a defensible primary choice, despite balanced allocation. Balance reduces one source of concern; it does not erase the evidence about variability.

Worked Example: Unequal Allocation and Unequal Spreads

Consider hypothetical independent groups with sample sizes \(8,20,40\), means \(20,24,22\) units, and standard deviations \(2,5,8\) units. Suppose the observations within groups are reasonably well behaved, but the summaries indicate markedly different spreads. The allocation is also unbalanced, and the smallest group has the smallest variance while the largest group has the largest variance. These features argue against relying on a pooled variance.

For comparison, the grand mean is \(\bar{x}_{\cdot\cdot}=[8(20)+20(24)+40(22)]/68=1520/68\approx22.35\) units. The between-group sum of squares is \(SSB=8(20-22.35)^2+20(24-22.35)^2+40(22-22.35)^2\approx103.53\). Thus \(MSB\approx103.53/2=51.76\). The pooled within-group sum of squares is \(SSW=7(2^2)+19(5^2)+39(8^2)=28+475+2496=2999\), giving \(MSW=2999/65\approx46.14\). The classical statistic is \(F=51.76/46.14\approx1.122\), with \(df=(2,65)\) and \(p\approx0.332\).

For Welch’s test, the weights are \(w_1=8/4=2.000\), \(w_2=20/25=0.8000\), and \(w_3=40/64=0.6250\), giving \(W=3.425\). The weighted mean is \(\bar{x}_w=[2(20)+0.8(24)+0.625(22)]/3.425\approx21.30\) units. The weighted squared-deviation sum is \(2(20-21.30)^2+0.8(24-21.30)^2+0.625(22-21.30)^2\approx9.52\). Also, \(A\approx(1-2/3.425)^2/7+(1-0.8/3.425)^2/19+(1-0.625/3.425)^2/39\approx0.07278\). So \(F_W=(9.52/2)/[1+0.25(0.07278)]\approx4.674\), with \(df_1=2\) and \(df_2=8/[3(0.07278)]\approx36.64\). The approximate \(p\)-value is \(0.0155\).

At the 0.05 level, Welch’s test provides evidence that not all three population means are equal, whereas the classical test does not. This example illustrates why the procedures can lead to different conclusions when spreads differ and allocation is unbalanced. The choice should still be made from design and variability considerations—not by selecting whichever test gives the smaller \(p\)-value.

Worked Example: Similar Spreads and Roughly Comparable Group Sizes

Suppose three hypothetical independent groups have sample sizes \(12,15,18\), means \(30,32,31\) units, and standard deviations \(4,4.5,4.2\) units. Assume the design and diagnostic plots support independence and show no severe departures from within-group normality. The standard deviations are fairly similar, the sample sizes are reasonably comparable, and there is no substantive reason to expect different population variances. In this setting, the classical test is a reasonable primary analysis.

The grand mean is \(\bar{x}_{\cdot\cdot}=[12(30)+15(32)+18(31)]/45=31.07\) units. Then \(SSB=12(30-31.07)^2+15(32-31.07)^2+18(31-31.07)^2\approx26.80\), so \(MSB=26.80/2=13.40\).

The within-group sum of squares is \(SSW=11(4^2)+14(4.5^2)+17(4.2^2)=176+283.5+299.88=759.38\). Therefore \(MSW=759.38/(45-3)=18.08\), and \(F=13.40/18.08\approx0.741\), with \(df=(2,42)\) and \(p\approx0.482\). At \(\alpha=0.05\), the result does not provide evidence against equal population means. Welch would also be a valid alternative, but the available information gives little practical reason to prefer it here.

A Practical Decision and Reporting Process

1
Check the design and distributions.
Confirm that observations are independent and inspect group distributions for severe skewness or influential outliers. A variance-robust test does not fix dependence or problematic outliers.
2
Summarize each group’s spread and size.
Report \(n_i\) and \(s_i\) (or \(s_i^2\)) for every group. Judge the size of the differences in context, and note whether allocation is balanced or uneven.
3
Choose a method before interpreting the result.
Use the classical test when a common variance is plausible and the allocation and diagnostics are not concerning. Prefer Welch when meaningful variance differences are plausible, especially with unequal group sizes.
4
State the choice and its basis.
Report the test, statistic, degrees of freedom, \(p\)-value, and the design or variability considerations supporting the method. If useful, a secondary analysis can assess sensitivity, but should not replace a clear primary rationale.

Common Mistake: Treating a Variance Pretest as an Automatic Switch

The mistaken reasoning is: “The variance test had \(p>0.05\), so the population variances are equal; therefore, the classical test is required.” A nonsignificant result means only that the pretest did not provide sufficient evidence against its own null hypothesis. It does not prove equal variances, and a pretest can have little power to detect important differences when samples are small. The opposite automatic rule is also problematic: a significant variance pretest does not by itself describe whether the differences are consequential for the planned mean comparison.

Instead, evaluate groupwise spreads together with sample sizes and knowledge of the design. If heterogeneity is credible, Welch is often a sensible choice because it does not impose a common-variance model. If homogeneity is plausible and the design supports pooling, the classical test remains appropriate. Report the reasoning; do not imply that a pretest established an assumption as true.

Key takeaway: Choose between classical and Welch one-way tests using the variance pattern, allocation, design knowledge, and inferential goal. Balanced samples can make the classical test more robust to unequal variances, but they do not guarantee equal variances or exact calibration.

Check Your Understanding

Use the assumptions and decision principles above to answer these questions.

  1. What assumptions does the classical one-way \(F\) test require for its exact reference distribution, and which of these does Welch relax?
  2. Why is a nonsignificant variance pretest not proof that population variances are equal?
  3. Two studies have the same groupwise standard deviations. Why might the study with strongly unequal group sizes raise greater concern about using a pooled-variance test?
  4. What does balanced allocation help with when variances differ, and what does it not guarantee?
  5. In a report using Welch’s test, what information should accompany the test statistic so readers can interpret the result?
Clinical Biostats