R/custom.probdeg.R

In metaSeq: Meta-analysis of RNA-Seq count data in multiple studies

custom.probdeg <-
function (Mg, Dg, Mn, Dn, prec = 2) 
{
    tot <- length(Mn)
    gens <- names(Mg)
    Mruido <- round(Mn, prec)
    Druido <- round(Dn, prec)
    Mgen <- round(Mg, prec)
    Dgen <- round(Dg, prec)
    MDgen <- na.omit(cbind(Mgen, Dgen))
    MDunic <- unique(MDgen)
    Nres <- apply(MDunic, 1, NOISeq:::n.menor, S1 = Mruido, S2 = Druido)
    lugares <- apply(MDgen, 1, NOISeq:::busca, S = MDunic)
    Nconj <- Nres[lugares]
    names(Nconj) <- names(lugares)
    laprob <- Nconj/tot
    laprob <- laprob[gens]
    names(laprob) <- gens
    Nconj <- Nconj[gens]
    names(Nconj) <- gens
    laprob <- list(prob = laprob, numDE = Nconj, numNOISE = tot)
    laprob
}