use log1p() where possible

log1p(-x) should be more accurate and faster than log(1-x)
restrict MCMC to sets that have observed elements