R/jtFunctions.R

Defines functions jt jtPTest

# 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
}

Try the gMWT package in your browser

Any scripts or data that you put into this service are public.

gMWT documentation built on March 26, 2020, 5:49 p.m.