Beispiel #1
0
 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);
 }
Beispiel #2
0
 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);
 }