[R] error in mood.test

[Wolfgang Raffelsberger]

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