# TODO: Add comment
# 
# Author: E.Korsching 10.9.2009
###############################################################################



Wilcoxon.Mann.Whitney.U.test <- function(x, ncx=ncx, ncy=ncy){
    #Wilcoxon rank sum test: mu=0: H0: >alpha: mu1=mu2, Ha: <=alpha: mu1 unequal mu2
	# can produce NaN values

    #
    f1 <- function(x, ncx, ncy){
        u.p <- wilcox.test(x=x[1:ncx], y=x[(ncx+1):(ncx+ncy)],
                           alternative="two.sided", mu=0, paired=F, exact=T, correct=T)$p.value

        return(u.p)
    }
    #
    result1 <- data.frame(apply(x, 1, f1, ncx=ncx, ncy=ncy))
    names(result1) <- c("u.p")
    #
    return(result1)
}