% nbpreview --theme material notebook.ipynb
      ╭───────────────────────────────────────────────────────────────────────╮
[22]:with pm.Model() as model:                                             │
      │     alpha = pm.Normal("alpha", mu=0, sd=10)                           │
      │     beta = pm.Normal("beta", mu=0, sd=1)                              │
      │     epsilon = pm.HalfCauchy("epsilon", beta=5)                        │
      │                                                                       │
      │     mu = pm.Deterministic("mu", var=alpha + beta * x)                 │
      │     y_pred = pm.Normal("y_pred", mu=mu, sd=epsilon, observed=y)       │
      │ y_pred                                                                │
      ╰───────────────────────────────────────────────────────────────────────╯

[22]:  y_pred∼Normal(𝑚𝑢=mu, 𝑠𝑖𝑔𝑚𝑎=f(epsilon))