(1) Gibbs steps outlined in section 3 (2) begun implementation in diff-expr-r.r (3) 3 references in writeup.bib for simulating from gamma distributions