About this Journal Submit a Manuscript Table of Contents
Advances in Decision Sciences
Volume 2011 (2011), Article ID 748580, 8 pages
http://dx.doi.org/10.1155/2011/748580
Research Article

A Nonparametric Two-Sample Wald Test of Equality of Variances

1Centre for Computer-Assisted Research Computation and Its Applications, School of Mathematical and Physical Sciences, University of Newcastle, Newcastle, NSW 2308, Australia
2Centre for Statistical and Survey Methodology, School of Mathematics and Applied Statistics, University of Wollongong, Wollongong, NSW 2522, Australia
3School of Mathematical and Physical Sciences, University of Newcastle, Newcastle, NSW 2308, Australia

Received 31 August 2011; Accepted 17 December 2011

Academic Editor: YanXia Lin

Copyright © 2011 David Allingham and J. C. W. Rayner. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

We develop a test for equality of variances given two independent random samples of observations. The test can be expected to perform well when both sample sizes are at least moderate and the sample variances are asymptotically equivalent to the maximum likelihood estimators of the population variances. The test is motivated by and is here assessed for the case when both populations sampled are assumed to be normal. Popular choices of test would be the two-sample test if normality can be assumed and Levene’s test if this assumption is dubious. Another competitor is the Wald test for the difference in the population variances. We give a nonparametric analogue of this test and call it the test. In an indicative empirical study when both populations are normal, we find that when both sample sizes are at least 25 the test is nearly as robust as Levene’s test and nearly as powerful as the test.

1. Introduction: Testing Equality of Variances for Two Independent Samples

In the two-sample problem, we are given two independent random samples and . The location problem attracts most attention. Assuming that the samples are from normal populations, the pooled -test is used to test equality of means assuming equal variances and Welch’s test can be used when equality of variances is suspect but normality is not. When normality is in doubt, the Wilcoxon test is often used.

The corresponding dispersion problem is of interest to confirm the validity of, for example, the pooled -test, and when dispersion differences are of direct interest. For example, testing for reduced variability is of interest in confirming natural selection (see, e.g., [1, Section  5.5]) and if some processes are in control. In exploratory data analysis, it is sensible to assess if one population is more variable than another. If it is, the cause may be that one population is bimodal and the other is not; the consequences of this in both the scenario and the model can then be explored in depth.

The study here introduces a new test of equality of variances. The asymptotic null distribution of the test statistic is, but, depending on the populations sampled from, this may or may not be satisfactory for small to moderate sample sizes.

To assess this, and other aspects of the proposed test, an indicative empirical study is undertaken. We give comparisons when both populations sampled are normal, when both samples have the same sample size, and for 5% level tests. First, we derive a finite sample correction to the critical values based on the asymptotic null distribution when sampling from normal populations. Different corrections would be needed when sampling from different distributions. Next, we show that in moderate samples the new test is nearly as powerful as the test when normality may be assumed and, finally, that it is nearly as robust as Levene’s test when normality is in doubt. Moore and McCabe [2, page  519] claim that the “ test and other procedures for inference about variances are so lacking in robustness as to be of little use in practice.” The new test gives a counterexample to that proposition.

We acknowledge that the new test is the most effective for at least moderate sample sizes. In the normal case, each random sample should have at least 25 observations, which is what we would expect of a serious study aiming at reasonable power that cannot be hoped for with samples of size 10 or so. See Section 4.

We are aware of more expansive comparative studies such as [3, 4]. Our goal here is not to emulate these studies but to show that the new test is competitive in terms of test size, robustness, and power.

In Section 2, the new test is introduced. Sections 3, 4, and 5 give the results of an empirical investigation when the populations sampled are assumed to be normal. In Section 3, we investigate test size, showing that the asymptotic critical values should only be used for moderate to large sample sizes. For smaller sample sizes, a finite sample correction to the asymptotic 5%   critical value is given. This results in actual test sizes between 4.6% and 5.3%.

In Section 4, we show that when normality holds, the new test is not as powerful as the Levene test for smaller sample sizes but overtakes it for moderate samples of about 25. The new test is always inferior to the optimal test. However, the test has power that approaches that of the test, being at least 95% that of the test throughout most of the parameter space in samples of at least 80.

In Section 5, we show that if we sample from -distributions with varying degrees of freedom, the test is highly nonrobust for small degrees of freedom, as is well known for fat-tailed distributions. When both populations sampled are gamma distributions, chosen to model skewness differences from normality, or -distributions, chosen to model kurtosis differences, the new test performs far better than the test, and is competitive with the Levene test.

That the new test gives a good compromise between power and robustness and is valid when normality does not hold are strong reasons for preferring the new test for sample sizes that are at least moderate, and normality is dubious.

2. Competitor Tests for Equality of Variance

Initially we assume that we have two independent random samples from normal populations, for and 2. We wish to test against the alternative , with the population means being unknown nuisance parameters. If , in which , , then the are the unbiased sample variances, and the so-called test is based on their quotient, , say. It is well known and will be confirmed yet again in Section 5 that the null distribution of , , is sensitive to departures from normality. If is the cumulative distribution function of the distribution, and if is such that , then the test rejects at the level when and when .

Common practice when normality is in doubt is to use Levene’s test or a nonparametric test such as Mood’s test. Levene’s test is based on the ANOVA test applied to the residuals. There are different versions of Levene’s test using different definitions of residual. The two most common versions use residuals based on the group means, , and the group medians, , in which is the median of the th sample. The latter is called the Brown-Forsythe test. Again it is well known that these tests are robust in that when the population variances are equal but the populations themselves are not normal, they achieve levels close to nominal. However, this happens at the expense of some power. As the empirical study in this paper is intended to be indicative rather than exhaustive, we will henceforth make comparisons only with the Levene test, based on a test statistic we denote by .

We now construct a new test that we call the test. For univariate parameters , a Wald test statistic of against the alternative is based on , the maximum likelihood estimator of , usually via the test statistic , where is the asymptotic variance of evaluated when . Under the null hypothesis, this test statistic has an asymptotic distribution. As well as being equivalent to the likelihood ratio test, the test is also a Wald test for testing against .

Rayner [5] derived the Wald test for testing against . The test statistic is

Being a Wald test, the asymptotic distribution of is , while its exact distribution is not immediately obvious. However, is a 1-1 function of , so the two tests are equivalent. Since the exact distribution of is available, the test is the obvious test to use.

In , the variances var() are estimated optimally using the Rao-Blackwell theorem. This depends very strongly on the assumption of normality. If normality is in doubt, then we can estimate var() using results given, for example, in [6]. For a random sample, with population and sample central moments and , , [6] gives that

Applying [6, 10.5] to the numerator of , may be estimated to by , or, equivalently, by , where is the unbiased sample variance. It follows that, to order , may be estimated by . We thus propose a robust alternative to , given by

in which are the fourth central sample moments for the th sample, . We call the test based on the test. As the sample sizes increase, the distributions of the standardised sample variances approach standard normality, the denominator in will approximate , and will have asymptotic distribution . Thus, if is the point for which the distribution has weight in the right hand tail, then the test rejects at approximately the level when .

We emphasise that although the motivation for the derivation of is under the assumption of sampling from normal populations, it is a valid test statistic for testing equality of variances no matter what the populations sampled.

If the sample variances are equal to or asymptotically equivalent to the maximum likelihood estimators of the population variances, as is the case when sampling from normal populations, then the test is a Wald test for equality of variances in the sense described above. Since it does not depend on any distributional assumptions about the data, it can be thought of as a nonparametric Wald test. As such, it can be expected to have good properties in large samples.

We note that all the above test statistics are invariant under transformations , for constants , , and and for and .

The next three sections report an empirical study when the distributions sampled are assumed to be normal. As this is an indicative study, we fix the samples sizes to be equal, , say, and the significance level to be 5% throughout.

3. Test Size under Normality

Under the null hypothesis, the distribution of is known exactly, the distribution of is known approximately, and, as mentioned above, the distribution of is known asymptotically. In analysing data, these distributions are used to determine values and critical values. We now investigate their use in determining test size, the probability of rejecting the null hypothesis when it is true.

Two empirical assessments of test size will now be undertaken. Since the test statistics are scale invariant, it is sufficient under the null hypothesis to take both population variances to be one.

In the first assessment, we assume normality. For various values of the common sample size , we estimate the 5% critical points for each test by generating 100,000 pairs of random samples of size , calculating the test statistics, ordering them, and hence identifying the 0.95th percentile. The estimated critical points of approach the 5% critical point 3.841. These estimated critical points will subsequently be used in the power study to give tests with test the size of exactly 5%.

To see the extent of the error caused by using the asymptotic critical point 3.841, Figure 1 gives the proportion of rejections in 100,000 pairs of random samples for sample sizes up to 100. For , the proportion of rejections is nearly 20% and although for this has dropped to nearly 7%, most users would hope for observed test sizes closer to 5%.

748580.fig.001
Figure 1: Proportion of rejections of the test using the 5% critical point 3.841 for sample sizes up to 100.

The application of the test is improved by the use a Bartlett-type correction. This was found by plotting the estimated 5% critical points against and using standard curve fitting techniques to find that the exact critical points were well approximated between and 100 by . For larger , it is sufficient to use the asymptotic 5% value 3.84146, the error being at most 0.8%. We checked the exact probabilities of rejection under the null hypothesis when applying the test with critical value , and all were between 4.6% and 5.3%.

For levels other than 5%, and when sample sizes are unequal, further empirical work needs to be done to find critical values. However, as this study was intended to be indicative, we leave extensive tabulation of critical values for another time.

4. Power under Normality

For the , Levene, and tests, we estimate the power as the proportion of rejections from 100,000 pairs of random samples of size when the first sample is from an population and the second is from an population with . To compare like with like, we use estimated critical values that give exact 5% level tests. It is apparent that for sample sizes less than about 20 the Levene test is more powerful than the test and that for a sample size between approximately 20 and 30 the test takes over from the Levene test; thereafter the test is always more powerful than the Levene test. This is shown in Figure 2.

748580.fig.002
Figure 2: Power of the 5% level test (solid line) and test (dashed line) for various sample sizes.

Both the Levene and tests are always less powerful than the test. This is explored in Figure 3 that compares the Levene test to the test in the left hand panel and the test to the test in the right hand panel. What is given is a contour plot of the regions in which the ratios of the power of the stated test to the test are less than 95%, between 95% and 99.99%, and greater than 99.99%. Generally, for any given and , it is clear that the power of the Levene test is at most that of the test. For example, it appears that for approximately the power of the test is always at least 95% of that of the test, whereas there is a considerable region where the power of the Levene test are less than 95% of that of the test.

fig3
Figure 3: Contour plots of the test (a) and test (b) relative to the test showing regions in which the power ratios are less than 95%, between 95% and 99.99%, and greater than 99.99%.

5. Robustness

Even if the test has good power, the test is of little value unless it is robust in the sense that when the distributions from which we are sampling are not from the nominated population (here, the normal distribution), the values are reasonably accurate. It is thus of interest to estimate the proportion of rejections when the null hypothesis is true and both the populations from which we sample are not normal. We have looked at variable kurtosis and skewness. First, variable kurtosis was considered via -distributions with varying degrees of freedom, say. Second, variable skewness was considered through gamma distributions with probability density functions for . Thus, the scale parameter is set at one and the scale parameter was varied. As and increase, the distributions become increasingly more normal, the distribution sampled will be close enough to normal that we could expect the proportion of rejections to be close to the nominal.

In Figure 4, we plot the proportion of rejections for the Levene, , and tests when sampling from distributions, for , and from gamma distributions with scale parameter 1 and shape parameter , for . We show curves for each test with common sample sizes and 20. The critical values used for the test are the from Section 3; the critical values used for the Levene test are those estimated by simulation to give exactly the nominal level when normality holds. Thus, this comparison favours the Levene test.

fig4
Figure 4: Test sizes for the (dots), (dashes), and (solid line) tests for sample sizes of 10 and 20 for (a), (b) -distributions with varying degrees of freedom, , from 1 to 20, and (c), (d) gamma distributions with scale parameter 1 and varying shape parameter, , from 0.1 to 30.

For these distributions, samples are increasingly nonnormal as and decrease. It is apparent that the test performs increasingly poorly as these parameters decrease.

When sampling from the -distribution, the Levene test generally has exact level closer to the nominal level than the test except for very nonnormal samples (small ). However, the level of the test is almost always reasonable, and while for very small , the level is not as close to the exact level as perhaps we may prefer, the same is the case for the Levene test.

When sampling from the gamma distribution when , the test outperforms the Levene test at small values of and is slightly inferior for larger values of , although both tests have exact test sizes acceptably close to the nominal test size. When , the test is uniformly closer to the nominal test size than the Levene test. Although we only display results for and , sample sizes of more than 20 yield very similar conclusions.

6. Conclusion

First, we reflect on testing for equality of variances when it is assumed that the populations sampled are normal. The test is both the likelihood ratio test and a Wald test, and is the appropriate test to apply. When normality does not hold, the test is no longer an asymptotically optimal test, and its well-known nonrobustness means that tests such as the Levene are more appropriate for small to moderate sample sizes. However, for sample sizes of about 25 or more, the test is more powerful than the Levene, and with the small sample corrected critical values, it holds its nominal significance level well. For these sample sizes, it can be preferred to the Levene test.

Second, consider testing for equality of variances when both samples are drawn from the same population. If that population is nominated, then the test may be applied after determining critical values or using values calculated by Monte Carlo methods. When the sample variances are asymptotically equivalent to the maximum likelihood estimators of the population variances (as, e.g., is the case when sampling from normal populations but not Poisson populations), the test is a nonparametric Wald test and hence will have good power in sufficiently large samples. If the population is not specified, the test can be confidently applied when the sample sizes are large, using the asymptotic null distribution to calculate values or critical values.

References

  1. B. F. J. Manly, The Statistics of Natural Selection on Animal Populations, Chapman and Hall, London, UK, 1985.
  2. D. S. Moore and G. P. McCabe, Introduction to the Practice of Statistics, W. H. Freeman, New York, NY, USA, 5th edition, 2006.
  3. W. J. Conover, M. E. Johnson, and M. M. Johnson, “A comparative study of tests of homogeneity of variances, with applications to the outer continental shelf bidding data,” Technometrics, vol. 23, no. 4, pp. 351–361, 1981.
  4. D. D. Boos and C. Brownie, “Bootstrap methods for testing homogeneity of variances,” Technometrics, vol. 31, no. 1, pp. 69–82, 1989. View at Publisher · View at Google Scholar · View at Zentralblatt MATH
  5. J. C. Rayner, “The asymptotically optimal tests,” Journal of the Royal Statistical Society Series D, vol. 46, no. 3, pp. 337–346, 1997.
  6. A. Stuart and J. K. Ord, Kendall's Advanced Theory of Statistics, vol. 1, Edward Arnold, London, UK, 6th edition, 1994.