Bayesian Regression

Simulate data

Number of data points:

Model definition

\[model(x, \beta) = \beta_1 + \beta_2 x\]

Inference

Prior

Likelihood

\[\log p(y \mid \beta, x) = \sum_i \log \mathrm{Normal}(y_i \mid \beta_1 + \beta_2 x_i, \beta_3)\]

Sampling from Posterior

Sampling...

Posterior prediction

The idea is that we propagate the uncertainty represented by the parameter posterior p(β | data) through the model. The prediction interval also takes the measurement uncertainty into account.