penalized_nested_Gibbs function

Penalized nested Gibbs sampler for Bayesian linear regression