[Rd] Help: StatXact

[Torsten Hothorn]

Help needed!

Has anyone access to StatXact? I just hacked exact two-sided p-values for
rank tests (for package exactDistr, which will move to CRAN/contrib as
exactRankTests soon ;-) and would like to compare the results of my
implementation to that of StatXact. Could someone please calculate the
exact one-sided (both greater and less) and two-sided p-values?

# Data from the StatXact-4 manual, page 221
treat <- c(94, 108, 110, 90)
contr <- c(80, 94, 85, 90, 90, 90, 108, 94, 78, 105, 88)

My results: 

> wilcox.exact(x,y)

         Wilcoxon rank sum test 

data:  x and y 
W = 52, p-value = 0.8963 
alternative hypothesis: true mu is not equal to 0 

> wilcox.exact(x,y, alternative="l")

         Wilcoxon rank sum test 

data:  x and y 
W = 52, p-value = 0.5656 
alternative hypothesis: true mu is less than 0 

> wilcox.exact(x,y, alternative="g")

         Wilcoxon rank sum test 

data:  x and y 
W = 52, p-value = 0.4636 
alternative hypothesis: true mu is greater than 0 


Thanks!

Torsten

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-devel 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-devel-request@stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._