add example Gibbs sampler code

Just finished my Gibbs sampling code, so here it is. Seems to match `lm()` for the BayesBridge data apart from a few coefficients, and quick plots show those variables are very collinear. In any case, subsequent runs are consistent, so at least it's converging to something.