[Rd] fligner.test (ctest) (PR#6739)
karel.zvara at mff.cuni.cz
karel.zvara at mff.cuni.cz
Mon Apr 5 14:13:00 CEST 2004
Full_Name: Karel Zvara
Version: 1.8.1
OS: MS Winows 2000
Submission from: (NULL) (195.113.30.163)
The test statistics of the fligner.test (ctest package) depends on the order of
cases:
> fligner.test(count~spray,data=InsectSprays)
Fligner-Killeen test for homogeneity of variances
data: count by spray
Fligner-Killeen:med chi-squared = 14.4828, df = 5, p-value =
0.01282
> fligner.test(count~spray,data=InsectSprays[order(InsectSprays$count),])
Fligner-Killeen test for homogeneity of variances
data: count by spray
Fligner-Killeen:med chi-squared = 5.4527, df = 5, p-value =
0.3632
> fligner.test(count~spray,data=InsectSprays[order(sample(72)),])
Fligner-Killeen test for homogeneity of variances
data: count by spray
Fligner-Killeen:med chi-squared = 2.7349, df = 5, p-value =
0.7408
The problem lies, I think, in the ordering of centered values x:
x <- unlist(tapply(x, g, function(u) u - median(u)))
More information about the R-devel
mailing list