with pm.Model() as model_5_7: a = pm.Normal(‘a’, mu=10, sd=100) bn = pm.Normal(‘bn’, mu=0, sd=1, shape=2) sigma = pm.Uniform(‘sigma’, lower=0, upper=1) mu = pm.Deterministic(‘mu’, a + bn[0] * dcc[‘neocortex.perc’] + bn[1] * dcc[‘log_mass’]) kcal = pm.Normal(‘kcal’, mu=mu, sd=sigma, observed=dcc…

I would say it is. You can also test if increasing the number of tuning steps helps (you are currently using 1000, you can increase to 3000 for instance).

There were 6 divergences after tuning. Increase `target_accept` or reparameterize

Questions

ayao December 21, 2020, 12:19am 6

thank you!~

Topic		Replies	Views
What to do with divergences in my case Questions	7	2495	September 21, 2017
Advice with hierarchical model Questions	9	2258	October 8, 2018
Overcoming divergence in MvNormal covariance model? Questions	1	481	January 28, 2019
Reducing divergences : Implementation challenges Questions	7	2213	August 3, 2021
There were 12 divergences after tuning. Increase `target_accept` or reparameterize Questions	3	2769	August 26, 2021

There were 6 divergences after tuning. Increase `target_accept` or reparameterize

Related topics