Gibbs sampler corrected, it works OK now

Numerical instability when low number of samples are drawn from POI
has been dealt with. If really no samples are drawn from POI, it
end up with error().