doing MCMC with P as free param