vartestn

Multiple-sample tests for equal variances

Syntax

vartestn(x)

vartestn(x,Name,Value)

vartestn(x,group)

vartestn(x,group,Name,Value)

p = vartestn(___)

[p,stats] = vartestn(___)

Description

vartestn(x) returns a summary table of statistics and a box plot for a Bartlett test of the null hypothesis that the columns of data vector x come from normal distributions with the same variance. The alternative hypothesis is that not all columns of data have the same variance.

vartestn(x,Name,Value) returns a summary table of statistics and a box plot for a test of unequal variances with additional options specified by one or more name-value pair arguments. For example, you can specify a different type of hypothesis test or change the display settings for the test results.

vartestn(x,group) returns a summary table of statistics and a box plot for a Bartlett test of the null hypothesis that the data in each categorical group comes from normal distributions with the same variance. The alternative hypothesis is that not all groups have the same variance.

vartestn(x,group,Name,Value) returns a summary table of statistics and a box plot for a test of unequal variances with additional options specified by one or more name-value pair arguments. For example, you can specify a different type of hypothesis test or change the display settings for the test results.

p = vartestn(___) also returns the p-value of the test, p, using any of the input arguments in the previous syntaxes.

[p,stats] = vartestn(___) also returns the structure stats containing information about the test statistic.

Examples

Test Data for Equal Variances

Load the sample data.

load examgrades

Test the null hypothesis that the variances are equal across the five columns of data in the students’ exam grades matrix, grades.

vartestn(grades)

ans = 7.9086e-08

The low $$p$$-value, p = 0, indicates that vartestn rejects the null hypothesis that the variances are equal across all five columns, in favor of the alternative hypothesis that at least one column has a different variance.

Test Grouped Data for Equal Variances

Load the sample data.

load carsmall

Test the null hypothesis that the variances in miles per gallon (MPG) are equal across different model years.

vartestn(MPG,Model_Year)

ans = 0.8327

The high $$p$$-value, p = 0.83269, indicates that vartestn does not reject the null hypothesis that the variances in miles per gallon (MPG) are equal across different model years.

Test for Equal Variances Using Levene’s Test

Load the sample data.

load carsmall

Use Levene’s test to test the null hypothesis that the variances in miles per gallon (MPG) are equal across different model years.

p = vartestn(MPG,Model_Year,'TestType','LeveneAbsolute')

p = 0.6320

The high $$p$$-value, p = 0.63195, indicates that vartestn does not reject the null hypothesis that the variances in miles per gallon (MPG) are equal across different model years.

Test for Equal Variances Using the Brown-Forsythe Test

Load the sample data.

load examgrades

Test the null hypothesis that the variances are equal across the five columns of data in the students’ exam grades matrix, grades, using the Brown-Forsythe test. Suppress the display of the summary table of statistics and the box plot.

[p,stats] = vartestn(grades,'TestType','BrownForsythe','Display','off')

p = 1.3121e-06

stats = struct with fields:
    fstat: 8.4160
       df: [4 595]

The small $$p$$-value, p = 1.3121e-06, indicates that vartestn rejects the null hypothesis that the variances are equal across all five columns, in favor of the alternative hypothesis that at least one column has a different variance.

Input Arguments

x — Sample data

matrix | column vector

Sample data, specified as a matrix or column vector. If a grouping variable group is specified, then x must be a column vector. If a grouping variable is not specified, x must be a matrix. In either case, vartestn treats NaN values as missing values and ignores them.

Data Types: single | double

group — Grouping variable

categorical array | logical or numeric vector | character array | string array | cell array of character vectors

Grouping variable, specified as a categorical array, logical or numeric vector, character array, string array, or cell array of character vectors with one row for each element of x. Each unique value in a grouping variable defines a group. vartestn treats NaN values as missing values and ignores them.

For example, if Gender is a cell array of character vectors with values 'Male' and 'Female', you can use Gender as a grouping variable to test your data by gender.

Example: Gender

Data Types: categorical | single | double | logical | string | cell | char

Name-Value Pair Arguments

Specify optional comma-separated pairs of Name,Value arguments. Name is the argument name and Value is the corresponding value. Name must appear inside quotes. You can specify several name and value pair arguments in any order as Name1,Value1,...,NameN,ValueN.

Example: 'TestType','BrownForsythe','Display','off' specifies a Brown-Forsythe test and omits the plot of the results.

'Display' — Display settings for test results

'on' (default) | 'off'

Display settings for test results, specified as the comma-separated pair consisting of 'Display' and one of the following.

'on'	Display a box plot and table of summary statistics.
'off'	Do not display a box plot and table of summary statistics.

Example: 'display','off'

'TestType' — Type of hypothesis test

'Bartlett' (default) | 'LeveneQuadratic' | 'LeveneAbsolute' | 'BrownForsythe' | 'OBrien'

Type of hypothesis test to perform, specified as the comma-separated pair consisting of 'TestType' and one of the following.

'Bartlett'	Bartlett’s test.
'LeveneQuadratic'	Levene’s test computed by performing ANOVA on the squared deviations of the data values from their group means.
'LeveneAbsolute'	Levene’s test computed by performing ANOVA on the absolute deviations of the data values from their group means.
'BrownForsythe'	Brown-Forsythe test computed by performing ANOVA on the absolute deviations of the data values from the group medians.
'OBrien'	O’Brien’s modification of Levene’s test with W = 0.5.

Example: 'TestType','OBrien'

Output Arguments

p — p-value

scalar value in the range [0,1]

p-value of the test, returned as a scalar value in the range [0,1]. p is the probability of observing a test statistic as extreme as, or more extreme than, the observed value under the null hypothesis. Small values of p cast doubt on the validity of the null hypothesis.

stats — Test statistics

structure

Test statistics for the hypothesis test, returned as a structure containing:

• chistat: Value of the test statistic.

• df: Degrees of freedom of the test.

More About

Bartlett’s Test

Bartlett’s test is used to test whether multiple data samples have equal variances, against the alternative that at least two of the data samples do not have equal variances.

The test statistic is

$$T=\frac{\left(N-k\right)\ln s_p{}^2-{\displaystyle \sum _{i=1}^k\left(N_i-1\right)\ln s_i{}^2}}{1+\left(1/\left(3\left(k-1\right)\right)\right)\left(\left({\displaystyle \sum _{i=1}^{k}1/\left(N_i-1\right)}\right)-1/\left(N-k\right)\right)},$$

where $$s_i{}^2$$ is the variance of the ith group, N is the total sample size, N_i is the sample size of the ith group, k is the number of groups, and $$s_p{}^2$$ is the pooled variance. The pooled variance is defined as

$$s_p{}^2={\displaystyle \sum _{i=1}^k\left(N_i-1\right)s_i{}^2/\left(N-k\right).}$$

The test statistic has a chi-square distribution with k – 1 degrees of freedom under the null hypothesis.

Bartlett’s test is sensitive to departures from normality. If your data comes from a nonnormal distribution, Levene’s test could provide a more accurate result.

Levene, Brown-Forsythe, and O’Brien Tests

The Levene, Brown-Forsythe, and O’Brien tests are used to test if multiple data samples have equal variances, against the alternative that at least two of the data samples do not have equal variances.

The test statistic is

$$W=\frac{\left(N-k\right){\displaystyle \sum _{i=1}^k{N}_i{\left({\overline{Z}}_{i.}-{\overline{Z}}_{\mathrm{..}}\right)}^2}}{\left(k-1\right){\displaystyle \sum _{i=1}^k{\displaystyle \sum _{j=1}^{N_i}{\left(Z_{ij}-{\overline{Z}}_{i.}\right)}^2}}},$$

where N_i is the sample size of the ith group, and k is the number of groups. Depending on the type of test specified with the TestType name-value pair arguments, Z_ij can have one of four definitions:

• If you specify LeveneAbsolute, vartestn uses $$Z_{ij}=\left|Y_{ij}-{\overline{Y}}_{i.}\right|$$, where $${\overline{Y}}_{i.}$$ is the mean of the ith subgroup.

• If you specify LeveneQuadratic, vartestn uses $$Z_{ij}{}^2={\left(Y_{ij}-{\overline{Y}}_{i.}\right)}^2$$, where $${\overline{Y}}_{i.}$$ is the mean of the ith subgroup.

• If you specify BrownForsythe, vartestn uses $$Z_{ij}=\left|Y_{ij}-{\tilde{Y}}_{i.}\right|$$, where $${\tilde{Y}}_{i.}$$ is the median of the ith subgroup.

• If you specify OBrien, vartestn uses

  $$Z_{ij}=\frac{\left(0.5+n_i-2\right)n_i{\left(y_{ij}-{\overline{y}}_i\right)}^2-0.5\left(n_i-1\right){\sigma }_i{}^2}{\left(n_i-1\right)\left(n_i-2\right)},$$

  where n_i is the size of the ith group, σ_i² is its sample variance.

In all cases, the test statistic has an F-distribution with k – 1 numerator degrees of freedom, and N – k denominator degrees of freedom.

The Levene, Brown-Forsythe, and O’Brien tests are less sensitive to departures from normality than Bartlett’s test, so they are useful alternatives if you suspect the samples come from nonnormal distributions.