I just stumbled upon this series (link) of tutorials on statistical modelling. What I find very interesting is the parallel of frequentist and bayesian methods (with STAN code included). Here you can find how to run tests such as ANOVA and the like in bayesian way. I think it wold be cool to convert these tutorials to pyMC3! I may start working on that, as I am trying to figure out how to run a multi factor ANOVA myself.
The INLA tutorial on this is excellent - I found it a few days ago as well.
I wonder who is behind these tutorials.