R/each.Fisher.ignore.test.R

each.Fisher.ignore.test <-
function (p) 
{
    loc <- setdiff(1:length(p), which(is.na(p)))
    if (length(loc) >= 1) {
        p <- p[loc]
        Xsq <- -2 * sum(log(p))
        pval <- 1 - pchisq(Xsq, df = 2 * length(p))
        return(pval)
    }
    else {
        return(NA)
    }
}

Try the metaSeq package in your browser

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

metaSeq documentation built on Nov. 8, 2020, 7:26 p.m.