Code covered by the BSD License  

Highlights from
FKtest

4.0

4.0 | 1 rating Rate this file 12 Downloads (last 30 days) File Size: 3.13 KB File ID: #25040

FKtest

by

 

15 Aug 2009 (Updated )

Fligner-Killeen test for homogeneity of variances.

| Watch this File

File Information
Description

There are several tests for homogeneity of variances. Fligner and Killeen (1976) suggest ranking (i) the |x_ij - mean| as normal type scores (a) by (0.5 + i/2(N+1)). Where i is the group and j the observation. The replacement of the mean by the median, a modification, is an attempt to improve the robustness of the test. From these scores is formulated a statistic based on a Chi-squared or a F distribution.

According to Conover et al. (1981), the Fligner-Killeen test by median is one of the best tests to use on the basis of robustness and power. Also, by several simulations, the Type I error rate and power is slightly larger when F approximation was used than when the Chi-squared approximation was used. Some of these tests are very sensitive to outliers, but Fligner test is not. Fligner test is the most robust against departures from normality.

Syntax: function FKtest(X,o)
      
Inputs:
X - Data matrix; Data=column 1, Group=column 2
o - By median=1 (default), mean=2

Outputs:
Complete analysis of the homogeneity of variances test by a Chi-squared and F approximation

Required Products Statistics Toolbox
MATLAB release MATLAB 7 (R14)
Tags for This File   Please login to tag files.
Please login to add a comment or rating.
Comments and Ratings (2)
16 Jul 2012 Alexander

I just verified this using fligner.test in R (used the data provided by the author in the commments).

One issue I had was the line of code:

A = [a' X(:,2)];

Had to change this to...

A = [a X(:,2)];

To avoid an error with concatenation

31 Aug 2011 James

Antonio: have you validated your code by comparing the results between yours and those given by R fligner.test?

Updates
17 Aug 2009

It was added an appropriate format to cite this file.

19 Aug 2009

Text was improved.

Contact us