rpowell at horizonenv.com
2006-Jun-14 13:41 UTC
[Rd] Bug in nortest cvm.test package (PR#8980)
I believe there to be a bug in the cvm.test module of the nortest package authored by Juergen Gross. I do not know how to contact the author directly. I've been running some normality tests using the nortest package. For some of my datasets the Cramer-von Mises normality test generates an extremely high probability (e.g., 1.637e+31) and indicates normality when the other tests do not. I don't believe the generated probability should exceed a value of "1". Below are the results of the tests and then, below that, I've provided the column of the datafile that results in the high value. Sincerely, Bob Powell Robert M. Powell Senior Project Scientist/Project Mgr. Horizon Environmental Corporation Suite 1700, Fisher Building 3011 W. Grand Blvd. Detroit, MI 48202 Tel (313) 874-0124 Fax (313) 874-0135 Cell (248) 977-0317 http://www.horizonenv.com > shapiro.test(Mn.Total..ug.kg) Shapiro-Wilk normality test data: Mn.Total..ug.kg W = 0.5269, p-value < 2.2e-16 > shapiro.test(ln.Mn.Total..ug/kg) Error in inherits(x, "factor") : object "ln.Mn.Total..ug" not found > shapiro.test(ln.Mn.Total..ug.kg) Shapiro-Wilk normality test data: ln.Mn.Total..ug.kg W = 0.943, p-value = 1.301e-06 > sf.test(Mn.Total..ug.kg) Shapiro-Francia normality test data: Mn.Total..ug.kg W = 0.5173, p-value < 2.2e-16 > sf.test(ln.Mn.Total..ug.kg) Shapiro-Francia normality test data: ln.Mn.Total..ug.kg W = 0.9403, p-value = 3.461e-06 > cvm.test(Mn.Total..ug.kg) Cramer-von Mises normality test data: Mn.Total..ug.kg W = 4.024, p-value = 1.637e+31 > cvm.test(ln.Mn.Total..ug.kg) Cramer-von Mises normality test data: ln.Mn.Total..ug.kg W = 0.356, p-value = 7.647e-05 > lillie.test(Mn.Total..ug.kg) Lilliefors (Kolmogorov-Smirnov) normality test data: Mn.Total..ug.kg D = 0.2325, p-value < 2.2e-16 > lillie.test(ln.Mn.Total..ug.kg) Lilliefors (Kolmogorov-Smirnov) normality test data: ln.Mn.Total..ug.kg D = 0.0869, p-value = 0.002011 OK, next the data column: Mn-Total, ug/kg 460000 400000 510000 600000 770000 180000 210000 1000000 2600000 490000 260000 760000 840000 400000 430000 230000 660000 1200000 370000 230000 290000 320000 400000 660000 620000 1100000 350000 450000 880000 960000 570000 870000 820000 970000 1100000 450000 730000 390000 640000 380000 340000 1400000 870000 260000 430000 290000 290000 550000 540000 240000 470000 650000 390000 380000 410000 220000 400000 380000 500000 1000000 520000 690000 500000 520000 260000 630000 630000 410000 290000 1300000 860000 600000 450000 330000 390000 580000 270000 460000 360000 300000 3600000 540000 370000 460000 150000 560000 630000 810000 1300000 470000 620000 520000 540000 440000 1900000 760000 750000 440000 550000 3300000 780000 3600000 1400000 840000 450000 970000 610000 490000 280000 960000 350000 1100000 770000 510000 460000 590000 760000 610000 550000 460000 470000 650000 370000 170000 630000 530000 720000 850000 570000 860000 6000000 1700000 1400000 470000 1200000 360000 370000 1100000 310000 430000 480000 670000 350000 530000 380000 490000 490000 380000 630000 370000 710000 1500000 380000 550000 560000 560000 350000 340000 250000 390000 670000 720000 440000 380000 240000 340000 350000 310000 250000 320000 360000 360000 490000 490000 490000 710000 670000 1200000 530000 500000 490000
rpowell at horizonenv.com wrote:> I believe there to be a bug in the cvm.test module of the nortest > package authored by Juergen Gross. I do not know how to contact the > author directly.What about typing library(help=nortest) It reveals the address at once. Please send the report to the package maintainer. Please do not send bug reports on contributed packages to R-bugs! Uwe Ligges> I've been running some normality tests using the nortest package. For > some of my datasets the Cramer-von Mises normality test generates an > extremely high probability (e.g., 1.637e+31) and indicates normality > when the other tests do not. I don't believe the generated > probability should exceed a value of "1". > > Below are the results of the tests and then, below that, I've > provided the column of the datafile that results in the high value. > > Sincerely, > Bob Powell > > Robert M. Powell > Senior Project Scientist/Project Mgr. > Horizon Environmental Corporation > Suite 1700, Fisher Building > 3011 W. Grand Blvd. > Detroit, MI 48202 > Tel (313) 874-0124 > Fax (313) 874-0135 > Cell (248) 977-0317 > http://www.horizonenv.com >[SNIP]