Yeah! it works, i know it because of the degrees of freedom. This is what i get now
> g <- interaction(ano, edadysexo, zona, estacion) > fligner.test(rojos ~ g) Fligner-Killeen test of homogeneity of variances data: rojos by g Fligner-Killeen:med chi-squared = 249.7591, df = 87, p-value < 2.2e-16 I didn't repair on d.f. I would pay more attention next time. Thanks a lot, ----- Mario Garrido Escudero PhD student Dpto. de Biología Animal, Ecología, Parasitología, Edafología y Qca. Agrícola Universidad de Salamanca -- View this message in context: http://r.789695.n4.nabble.com/different-results-from-fligner-test-tp4283312p4283753.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.