Hi, guys

I am trying to make a model for epidemic growth these days and i got some problem…

at first, my model was based on simple logistic growth model. and i found a reference using bayesian gamma GLM so i tried that way but actually confused about model…

through that way, i made model setting all the priors with simple non-informative ones just for trial.but it’s not working with ‘bad initial energy error’

with pm.Model() as model:

```
b0 = pm.Normal('b0',mu=0, sigma=10)
b1 = pm.HalfNormal('b1', sigma=10)
b2 = - pm.HalfNormal('b2', sigma=10)
sig = 2.5
theta = (b0 + b1*data['c']+b2*data['c^2'])
mu_ = 1/theta
y = pm.Gamma('y', mu = mu_, sigma = sig, observed=data['dc'].values)
```

i am a newbie of bayesian,PYMC so any kind of help would be really appreciated. thanks!