Dear list,
when running a mood.test() (part of package "stats") on slightly longer
vectors (than the example from the help-page) we get the error-message
shown below : once both vectors tested are of length 50 this error oocurs.
Note, that this problem didn't occur with R-2.7.x (or even older versions).
> x <- rnorm(50,10,5)
> y <- rnorm(50,2,5)
>
> mood.test( x[1:45], y[1:45]) # works OK
Mood two-sample test of scale
data: x[1:45] and y[1:45]
Z = -0.1556, p-value = 0.8763
alternative hypothesis: two.sided
> mood.test( x, y) # works ONLY on R-2.7.x but NOT on
R-2.8.x
Mood two-sample test of scale
data: x and y
Z = NA, p-value = NA
alternative hypothesis: two.sided
Warning message:
In m * n * (N + 1L) * (N + 2L) * (N - 2L) :
NAs produced by integer overflow
>
> # for completness
> sessionInfo()
R version 2.8.1 (2008-12-22)
i386-pc-mingw32
locale:
LC_COLLATE=French_France.1252;LC_CTYPE=French_France.1252;LC_MONETARY=French_France.1252;LC_NUMERIC=C;LC_TIME=French_France.1252
attached base packages:
[1] stats graphics grDevices datasets tcltk utils
methods base
other attached packages:
[1] svSocket_0.9-5 svIO_0.9-5 R2HTML_1.59 svMisc_0.9-5
svIDE_0.9-5
loaded via a namespace (and not attached):
[1] tools_2.8.1
>
Thank's in advance,
Wolfgang and David
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Wolfgang Raffelsberger, PhD
Laboratoire de BioInformatique et Génomique Intégratives
CNRS UMR7104, IGBMC
1 rue Laurent Fries, 67404 Illkirch Strasbourg, France
Tel (+33) 388 65 3300 Fax (+33) 388 65 3276
wolfgang.raffelsberger (at) igbmc.fr
______________________________________________
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.
