R/account.for.prior.R

In saturnin: Spanning Trees Used for Network Inference

Defines functions account.for.prior

Documented in account.for.prior

account.for.prior <-
function(prob,q0){
  p0 <- 2/nrow(prob)
  ratio <- (p0*(1-q0))/(q0*(1-p0))
  prob.q0 <- apply(prob,
                    1:2,
                    function(y) (1 + ratio*((1-y)/y))^{-1})
  return(prob.q0)
}

Try the saturnin package in your browser