Re multinational regression, you should take a look at Multinomial hierarchical regression with multiple observations per group ("Bad energy issue")
Re pymc3 style of trace plot, there is an issue: https://github.com/arviz-devs/arviz/issues/94 but nobody is actively working on it yet…