R/jtFunctions.R

In gMWT: Generalized Mann-Whitney Type Tests

# Version: 30-11-2012, Daniel Fischer

jt <- function(X,g){
  goi <- sort(unique(g))
  diffProbs <- getComb(goi,"pairs",order=T)

  result <- 0
  for(sumRun in 1:dim(diffProbs)[1])
  {
    result <- result + sum(g==diffProbs[sumRun,1]) * sum(g==diffProbs[sumRun,2]) * getP.Cnaive(X[g==diffProbs[sumRun,1]],X[g==diffProbs[sumRun,2]])
  }
  result
}

jtPTest <- function(X,g,nper){
    permTestValues <- c(rep(0,nper))
    for(i in 1:nper){
	permGroups <- sample(g)
	permTestValues[i] <- jt(X,permGroups)
    }
    permTestValues
}