7.4.2. Assuming the observations are normal, do the processes have the same variance?

7. Product and Process Comparisons
7.4. Comparisons based on data from more than two processes

7.4.2. Assuming the observations are normal, do the processes have the same variance?

Before comparing means, test whether the variances are equal

Techniques for comparing means of normal populations generally assume the populations have the same variance. Before using these ANOVA techniques, it is advisable to test whether this assumption of homogeneity of variance is reasonable. The following procedure is widely used for this purpose.

Bartlett's Test for Homogeneity of Variances

Null hypothesis

Bartlett's test is a commonly used test for equal variances. Let's examine the null and alternative hypotheses.

H0: sigma1**2 - sigma2**2 = ... = sigmak**2

against

Test statistic

Assume we have samples of size n_i from the i-th population, i = 1, 2, . . . , k, and the usual variance estimates from each sample:

where

s(i)**2 = SUM[j=1 to n(i)](x(ij) - xbar(i.))**2/(n(i) - 1)

Now introduce the following notation: _j = n_j - 1 (the _j are the degrees of freedom) and

s**2 = SUM[i=1 to k][nu(i)*s(i)**2]/nu

The Bartlett's test statistic M is defined by

M = nu*LOG(s**2) - SUM[i=1 to k][nu(i)*LOG(s(i)**2)]

Distribution of the test statistic

When none of the degrees of freedom is small, Bartlett showed that M is distributed approximately as Chi-Square(k-1)

. The chi-square approximation is generally acceptable if all the n_i are at least 5.

Bias correction

This is a slightly biased test, according to Bartlett. It can be improved by dividing M by the factor C = 1 +(1/(3*(k-1))*{[SUM[i=1 to k][1/v(i)]] - 1/nu}

Instead of M, it is suggested to use M/C for the test statistic.

Bartlett's test is not robust

This test is not robust, it is very sensitive to departures from normality.

An alternative description of Bartlett's test, which also describes how Dataplot implements the test, appears in Chapter 1.

Gear Data Example (from Chapter 1):

An illustrative example of Bartlett's test

Gear diameter measurements were made on 10 batches of product. The complete set of measurements appears in Chapter 1. Bartlett's test was applied to this dataset leading to a rejection of the assumption of equal batch variances at the .05 critical value level. applied to this dataset

The Levene Test for Homogeneity of Variances

The Levene test for equality of variances

Levene's test offers a more robust alternative to Bartlett's procedure. That means it will be less likely to reject a true hypothesis of equality of variances just because the distributions of the sampled populations are not normal. When non-normality is suspected, Levene's procedure is a better choice than Bartlett's.

Levene's test and its implementation in DATAPLOT were described in Chapter 1. This description also includes an example where the test is applied to the gear data. Levene's test does not reject the assumption of equality of batch variances for these data. This differs from the conclusion drawn from Bartlett's test and is a better answer if, indeed, the batch population distributions are non-normal.