R/GeneratePvalues.R

In globalGSA: Global Gene-Set Analysis for Association Studies.

Defines functions GeneratePvalues

Documented in GeneratePvalues

GeneratePvalues <-
function(data, B, addit=FALSE, covariable=NULL, family=binomial) {
  pvalors <- numeric(ncol(data)-1)
  pvalors <- sort(runPvalues(data, addit, covariable,family))
  pvalors <- rbind(pvalors, t(sapply(1:B, function(i) sort(runPermut(data, addit, covariable, family)))))
  rownames(pvalors) <- NULL  
  return(pvalors)
}

Try the globalGSA package in your browser