R/sample_R.R

sample_R <-
function (e, priR) {
  S <- t(e) %*% e + priR[2] * diag(ncol(e))
  v <- priR[1]
  sig <- solve(rwish(solve(S), v))
  cov2cor(sig)
}
goldingn/BayesComm documentation built on May 17, 2019, 7:41 a.m.