Monday 7 October 2013

Bartlett's Test: Equality of variances

There are statistical tests that involve two or more datasets, and that rely on the assumption that all the datasets have the same variance. For example, the commonly used version of the Unpaired Student's t-test assumes that the two datasets have the same variance. 

Another widely used test that assumes that three or more datasets have the same variance is the one-way and two-way ANOVA test.

Bartlett's test can be used to determine whether the datasets (or groups) have the same variance, and relies on the assumption that all the datasets come from a Normal distribution.

Suppose we have $k$ groups/datasets.

We first calculate the variance for each of the $k$ groups as follows ($i=1,2,..,k$):
$\Large \sigma_i^2=\frac{1}{N_i-1}\sum_{n=1}^{N_i}(x_i[n]-\mu_i)^2$
where $N_i$ is the number of samples in the $ith$ group, and $\mu_i$ is the mean of the $ith$ group.
$\Large \mu_i=\frac{1}{N_i}\sum_{n=1}^{N_i}x_i[n]$

Next we calculate the pooled variance $\sigma_p^2$
$\Large \sigma_p^2=\frac{1}{N-k}(\sum_{i=1}^k(N_i-1)\sigma_i^2)$
where $N$ is the total number of samples for all $k$ groups, i.e.

Once we have the variances for the $k$ groups and the pooled variance, we are in a position to calculate the Bartlett statistic $X^2$, which is given by

$\Large X^2=\frac{(N-k)ln(\sigma_p^2)-\sum_{i=1}^k(N_i-1)ln(\sigma_i^2)}{1+\frac{1}{3(k-1)}(\sum_{i=1}^k(\frac{1}{N_i-1})-\frac{1}{N-k})}$

Now the statistic $X^2$ (approximately) has a chi-squared distribution $\chi_{k-1}^2$.

Assuming that the Null Hypothesis is that the variances for all $k$ groups are the same, we reject this (i.e. the test result is significant, hence the assumption that the variances are all the same is rejected) if the calculated $X^2$ exceeds $\chi_{(k-1),\alpha}^2$, the critical value of the chi-squared distribution.

This will be the case when $X^2$ falls within the right tail region of the chi-squared distribution.

The parameter $\alpha$ is 1 minus the significance level. For a 5% significance level, $\alpha=1-0.05=0.95$ - note that this is a one-tailed test, as the variance is always greater than zero.

To find the value $\chi_{(k-1),\alpha}^2$, you need to calculate the inverse CDF of the chi-squared distribution for given $\alpha$ and $k$. This can be done using various statistical packages, or using tables, or (yes you guessed correctly..) using my iOS App SciStatCalc.

No comments:

Post a Comment

Logistic Regression Calculator and ROC Curve Plotter

This blog post implements a Logistic Regression calculator for a binary output. Consider a binary outcome response variable \(Y\...