ok. thank you. I have now implemented the method of non center parametrization as suggested. I am just curious, Why is it so that doing
…
pm.Normal(“test”, mu=1e-10, sd=1e-20)
…
gives a “derivative is zero” error for the sampling, while shifting a the distribution by using the below works? As far as I can tell, the resulting distributions are equal.
....
test = pm.Normal("test", mu=0, sd=1)
shift = test*1e-20 + 1e-10
shifted = pm.Deterministic("shifted", shift)
....