[R] Help with perm.test please

[Torsten Hothorn]

> Now I understand that you don't need to be a
> weatherman to see which way the wind is blowing here
> and clearly you would reject the null of equality.
> 
> I thought that perm.test creates a rx2 contingency
> table and then does a network algorithm ala stat exact
> (or something like that). In this example it doesn't
> seem that the contingency table would be that big. Any
> suggestions or comments ?

no. perm.test uses another algorithm which basically builds a matrix where
the number of columns is the sum of all observations in the x and y
sample. However, if you have enough memory you can set 

#define PERM_MAX_N 200

to whatever you want in exactRankTests/src/permdist.c l.44 (710 in your
example :-)

I will think of something more flexible ...

Torsten

> 
> Thanks,
> 
> Steve
> 
> But I get the following error:
> 
> > perm.test(x,y,exact=T)
> Error in cperm(eq, m, paired) : N > 200 in cpermdistr2
> 
> 
> 
> __________________________________________________
> 
> Yahoo! Finance - Get real-time stock quotes
> 
> -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
> r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
> Send "info", "help", or "[un]subscribe"
> (in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
> _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
> 

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._