R/getPvalBH.R

In isglobal-brge/CNVassoc: Association analysis of CNVs and imputed SNPs incorporating uncertainty

getPvalBH <- function(x) {
        pvals <- as.numeric(unlist(x))
        p.adjust.M <- p.adjust.methods
        p.adj <- vapply(p.adjust.M, function(meth) p.adjust(pvals, 
            meth), rep(0, length(pvals)))
        dd <- data.frame(pval = pvals, pval.BH = p.adj[, 
            5])
        o <- order(dd[, 2])
        ddEnd <- data.frame(region = o, dd[o, 
            ])
        dimnames(ddEnd)[[1]] <- c(seq_len(nrow(ddEnd)))
        ddEnd
}