PGLS using BRMS model and uncertainty in data

I have a set of infraspecific regression analyses of 10 measurements, for 50+ species. I want to test whether the regression coefficients differ statistically between species using a PGLS, and as the regression coefficients differ in their confidence intervals, I would like to use a Bayesian approach. How might I go about doing this?