forward simulation seems to be working ok, but MCMC does not seem to converge with lbda=0.1, mu=0.9 and rad=2.0