ValueError: Bad initial energy: inf. The model might be misspecified

[junpenglao]

@erim

(‘outcome_likelihood’, array(-inf))

This is what causing the issue - you should not use epsilon = pm.Normal('epsilon', mu=0., sd= 100) as sigma for Normal distribution should be larger than 0. Try epsilon = pm.HalfNormal('epsilon', sd= 100)