[Rd] wilcox.exact function (PR#8856)
Peter Dalgaard
p.dalgaard at biostat.ku.dk
Fri May 12 16:48:23 CEST 2006
pat.hodgson at sympatico.ca writes:
> Full_Name: Patrick Hodgson
> Version: 2.0
> OS: solaris 2.9
> Submission from: (NULL) (65.94.128.161)
>
>
> The value reported for the parameter W in the function wilcox.exact appears to
> be incorrect. I have checked the reference in the help file for this function
> (Myles & Hollander 1973, as well as 2nd ed. 1999 by same authors) and it is
> clear that W is the sum of the ranks of the data set with the smaller number of
> members.
You should not use the bug repository to report on a (non-core)
package. However, it is the same W used by wilcox.test and
help(wilcox.test) has the explanation.
> Here is sample input/output:
> > x<-c(1,2,3,4,5,6,7,8,9,10)
> > y<-c(11,12,13,14,15)
> > library(exactRankTests)
> > wil<-wilcox.exact(x,y)
> > wil
>
> Exact Wilcoxon rank sum test
>
> data: x and y
> W = 0, p-value = 0.000666
> alternative hypothesis: true mu is not equal to 0
>
> In this case W should be the sum of the ranks of y (coincidentally also the
> values of the y data)= 11+12+13+14+15 = 65.
> W of zero is ridiculous.
>
> Similarly, in this case W should be 2+4+6+8 = 20, not 10 as reported.
> > x<-c(1,3,5,7,9)
> > y<-c(2,4,6,8)
> > wil<-wilcox.exact(x,y)
> > wil
>
> Exact Wilcoxon rank sum test
>
> data: x and y
> W = 10, p-value = 1
> alternative hypothesis: true mu is not equal to 0
>
> ______________________________________________
> R-devel at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-devel
>
--
O__ ---- Peter Dalgaard Øster Farimagsgade 5, Entr.B
c/ /'_ --- Dept. of Biostatistics PO Box 2099, 1014 Cph. K
(*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918
~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk) FAX: (+45) 35327907
More information about the R-devel
mailing list