R/ur.mat.R

ur.mat <-
  function(lambda,delta,y1,y2){
    
    logdpois1 <- apply(as.matrix(lambda),1,calc.ldpois,y=y1)
    logdpois.1<-apply(logdpois1,2,sum)
    
    logdpois2 <- apply(as.matrix(lambda/delta$nodes),1,calc.ldpois,y=y2)
    logdpois.2<-apply(logdpois2,2,sum)
    ur.mat<-exp(logdpois.1)*exp(logdpois.2)*(delta$weights)
    ur.mat
  }
MalathiSIDona/pathDESeq documentation built on May 8, 2019, 3:37 p.m.