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!