I am modeling data that I think is gamma distributed so I write a PyMC3 model as below.
import pymc3 as pm import arviz as az with pm.Model() as model2: