R/mayers_hypertest_correct.R

In ZW-xjtlu/exomePeak2Test: Easy test of the exomePeak2 performance.

#assay_dds <- matrix(rpois(100, lambda = 10),nrow = 10,ncol = 10)
#
mayers_hypertest_correct <- function(IP_vec, input_vec, correct_IP, correct_input){
  total_white_balls <- sum(IP_vec)
  total_black_balls <- sum(input_vec)
  balls_draw <- IP_vec + input_vec
  white_balls_drawn <- IP_vec
  pvalue <- phyper(q = white_balls_drawn,
                   m = total_white_balls,
                   n = total_black_balls,
                   k = balls_draw, lower.tail = F)
  p.adjust(pvalue, method = "BH")
}