Simplified version for debugging.
FIXED non-convergence issue by recalculating "prev.posterior" conditioned on current value of "sd.i" at EACH gibbs loop.
Reverted to old test run.