I get "Bad initial energy" even on one of the examples

Try adding init='adapt_diag' when you call pm.sample.
More here: Frequently Asked Questions

1 Like