R/sampleSig2.R

Defines functions sampleSig2

Documented in sampleSig2

sampleSig2 <-
function(y, Px, v0, gamma0) {
  out = rinvgamma(1, shape=v0/2 + length(y)/2, scale = (gamma0 + t(y) %*% Px %*% y)/2)
  return(out)
}

Try the ARTIVA package in your browser

Any scripts or data that you put into this service are public.

ARTIVA documentation built on May 1, 2019, 6:31 p.m.