Check this tutorial on PyMC3: https://github.com/fonnesbeck/intro_stat_modeling_2017/blob/...