powered by
Prior: sigma2 ~ Inv-Gamma(a0, b0) (shape a0, scale b0) Conditional: sigma2 | rest ~ Inv-Gamma(a0 + n/2, b0 + 0.5 * RSS)
update_sigma2_ig(y, X, beta, x, a0, b0)