void LSB::draw_beta_given_gamma(){ const Selector & inc(mod_->inc()); Ominv = inc.select(pri_->siginv()); Spd ivar = Ominv + inc.select(suf()->xtx()); Vec b = inc.select(suf()->xty()) + Ominv * inc.select(pri_->mu()); b = rmvn_suf(ivar, b); mod_->set_included_coefficients(b); }
void LS::draw_beta(){ ivar = pri_->siginv() + suf_->xtx(); ivar_mu = pri_->siginv() * pri_->mu() + suf_->xty(); ivar_mu = rmvn_suf(ivar, ivar_mu); mod_->set_beta(ivar_mu); }