I'm just currently getting into the world of Bayesian machine learning, with a lot of frequentest stats background, and I frequently find myself limited in my modeling by my lack of knowledge about what distributions to use and when. For example, I'll have continuous data that is strongly right skewed and strictly greater than 0 by nature of some physical process. Which distribution should I use to model this?
I know I could probably use Gamma, or inverse Gaussian, chi square, or maybe a few others, but I'm not sure the implications of picking one of the above, or another distribution altogether. I know that Gamma and inverse Gaussian are in the exponential family and so have conjugate priors and are a little easier to work with, but in the age of MCMC, I'm looking for the distribution that has the best fit overall, so are there any rules of thumb for picking between these distributions? Are there certain behaviors of those distributions that are worth watching out for? Should I be taking a different approach altogether and trying to transform the data to something normal by taking logs or square roots?