I'm having trouble with the wilcox.test command in R. To demonstrate the anomalous behavior of wilcox.test, consider 

> wilcox.test(c(1.5,5.5), c(1:10000), exact = F)$p.value
[1] 0.01438390
> wilcox.test(c(1.5,5.5), c(1:10000), exact = T)$p.value
[1] 6.39808e-07 (this calculation takes noticeably longer). 
> wilcox.test(c(1.5,5.5), c(1:20000), exact = T)$p.value
(R closes/crashes)


I believe that wilcox.test(c(1.5,5.5), c(1:10000), exact = F)$p.value yields a bad result because of the normal approximation which R uses when exact = F. Any suggestions for how to compute wilcox.test(c(1.5,5.5), c(1:20000), exact = T)$p.value? 

--Greg

	[[alternative HTML version deleted]]


