Book Image

Bayesian Analysis with Python - Second Edition

By : Osvaldo Martin
4.5 (2)
Book Image

Bayesian Analysis with Python - Second Edition

4.5 (2)
By: Osvaldo Martin

Overview of this book

The second edition of Bayesian Analysis with Python is an introduction to the main concepts of applied Bayesian inference and its practical implementation in Python using PyMC3, a state-of-the-art probabilistic programming library, and ArviZ, a new library for exploratory analysis of Bayesian models. The main concepts of Bayesian statistics are covered using a practical and computational approach. Synthetic and real data sets are used to introduce several types of models, such as generalized linear models for regression and classification, mixture models, hierarchical models, and Gaussian processes, among others. By the end of the book, you will have a working knowledge of probabilistic modeling and you will be able to design and implement Bayesian models for your own data science problems. After reading the book you will be better prepared to delve into more advanced material or specialized statistical modeling if you need to.
Table of Contents (11 chapters)
9
Where To Go Next?

Exercises

  1. Check the following definition of a probabilistic model. Identify the likelihood, the prior, and the posterior:
  1. For the model in exercise 1, how many parameters have the posterior? In other words, how many dimensions does it have?
  2. Write down Bayes' theorem for the model in exercise 1.
  1. Check the following model. Identify the linear model and identify the likelihood. How many parameters does the posterior have?
  1. For the model in exercise 1, assume that you have a dataset with 57 data points coming from a Gaussian with a mean of 4 and a standard deviation of 0.5. Using PyMC3, compute:
    • The posterior distribution
    • The prior distribution
    • The posterior predictive distribution
    • The prior predictive distribution

Tip: Besides pm.sample(), PyMC3 has other functions to compute samples.

  1. Execute model_g using NUTS (the default sampler) and then...