Bad initial energy: inf. The model might be misspecified

[junpenglao]

In that case try: trace = pm.sample(1000, tune=1000, init='adapt_diag')

The jitter sometimes returns invalid starting value.