should σ ∼Exponential (1) not translate to sigma <- rexp(n, **1**) instead of : https://github.com/corriebar/Statistical-Rethinking/blob/15f24e45214376f1eee223449514032a6b5c9a1c/Chapter_4/chapter4_Ex.Rmd#L129
