R/prior_unif_dirichlet.R

In saturnin: Spanning Trees Used for Network Inference

Defines functions prior_unif_dirichlet

Documented in prior_unif_dirichlet

prior_unif_dirichlet <-
function(p,r,Neq = 0.5*r^2){
  a = r*diag(numeric(r)+1)
  b = matrix(1,r,r)
  return((Neq/(r^2))*array(c(a,rep(b,p)),c(r,r,p,p)))
}

Try the saturnin package in your browser