Hi Kelvin,
Have you fit the model in the unconstrainted domain setting? What are the results there?
Does setting the positive only constraint cause issues with sampling?
Are you using controls in your model as well? If so, do you have any transformations applied to them before model fit?