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

I had the same issue.
I fixed it with changing slope bn to HalfNormal and increasing sd

bn = pm.Normal(‘bn’, mu=0, sd=1, shape=2)

You can change sd to 10 first then if it still not working - move to HalfNormal