gibbs sampler seems to do ok, but something weird with NormalVar object