R/minP.MiRKAT.R

In JiyuanHu/massMap: massMap: a two-stage microbial association mapping framework with advanced FDR control

minP.MiRKAT <-
function(Qtest,Q0){
  if(is.vector(Q0)){
    res = mean(Q0>= Qtest)
  }else{
    true.false = t(apply(Q0,1,function(x)x>=Qtest))
    pMiRKATk = colMeans(true.false)
    res =min(pMiRKATk)
  }
  res
}