thanks for the quick response. You’re right. I was thinking more from an API standpoint. There could be a builtin method to learn a parametric distribution from a set of samples, perhaps there is.
I found this recommendation from @junpenlao
I want to look further into the comparative efficiency vs fidelity of these methods, but at first glance, the method with pm.Potential looks like the best option.