https://barnesanalytics.com/bayesian-regression-with-pymc3-in-python/