R help - Oct 2009 - wilcox.exact() problem

Dear R friends,


here I write again about the wilcox.exact() problem. I want to compare two sets
of categorical data, and in one case it says "negative length vectors not
allowed", and in the other one I get the error "cannot allocate vector
of length ...".

On http://rapidshare.com/files/298621893/wilcox.exact_trouble.Rdata.html you can
download the data that cause the trouble. (Attention: 10 Downloads only, so you
may have to send me an e-mail so I can mail you the Rdata file directly.)
> library(exactRankTests)
  Package ?exactRankTests? is no longer under development.
  Please consider using package ?coin? instead.
> wilcox.exact(group_1,group_2,exact=T)
Error in .Call("cpermdist2", ma = as.integer(m), mb = as.integer(col),
:
  negative length vectors are not allowed
> wilcox.exact(group_3,group_4,exact=T)
Error in .Call("cpermdist2", ma = as.integer(m), mb = as.integer(col),
:
  cannot allocate vector of length 1621367542

1621367542 is the same as 2^30.59, so the memory can't be allocated on my
32-bit computer.

Can somebody point out the problem to me?


Kind regards,


David
-- 
Jetzt kostenlos herunterladen: Internet Explorer 8 und Mozilla Firefox 3.5 -
sicherer, schneller und einfacher! http://portal.gmx.net/de/go/chbrowser

On Tue, 27 Oct 2009, David Croll wrote:
> Dear R friends,
>
> here I write again about the wilcox.exact() problem. I want to compare two
sets of categorical data, and in one case it says "negative length vectors
not allowed", and in the other one I get the error "cannot allocate
vector of length ...".
>
> On http://rapidshare.com/files/298621893/wilcox.exact_trouble.Rdata.html
you can download the data that cause the trouble. (Attention: 10 Downloads only,
so you may have to send me an e-mail so I can mail you the Rdata file directly.)
>
>> library(exactRankTests)
>  Package ?exactRankTests? is no longer under development.
>  Please consider using package ?coin? instead.
>> wilcox.exact(group_1,group_2,exact=T)
> Error in .Call("cpermdist2", ma = as.integer(m), mb =
as.integer(col),  :
>  negative length vectors are not allowed
>> wilcox.exact(group_3,group_4,exact=T)
> Error in .Call("cpermdist2", ma = as.integer(m), mb =
as.integer(col),  :
>  cannot allocate vector of length 1621367542
>
> 1621367542 is the same as 2^30.59, so the memory can't be allocated on
my 32-bit computer.
>
> Can somebody point out the problem to me?
The sample size is too large (~4000 observations per group) to solve this 
problem exactly. The error message could maybe be improved, but the 
message is clear: This is too large to deal with.

However, this is not a problem. With several thousand observations, 
standard normal approximations should work sufficiently well. And if you 
don't believe it, then you can look at approximate solutions that draw a 
sufficiently large number of permutations. Both is easily available when 
using wilcox_test() in "coin" as the startup message of
"exactRankTests"
suggests.

hth,
Z
>
> Kind regards,
>
>
> David
> -- 
> Jetzt kostenlos herunterladen: Internet Explorer 8 und Mozilla Firefox 3.5
-
> sicherer, schneller und einfacher! http://portal.gmx.net/de/go/chbrowser
>
> ______________________________________________
> R-help at 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.