plot(
plot(),
plot(),
...
)
vline!([2]; color=["blue"], label="Real coeff", linestyle=[:dash])
https://github.com/sswatson/cheatsheets/blob/master/plotsjl-cheatsheet.pdf
surface(-5.0:0.1:5.0, 0.1:0.1:5.0,(x, y) -> pdf(prior, x, y))
xlabel!("mean")
ylabel!("variance")
heatmap(-3:0.1:3,-3:0.1:3, (x, y) -> pdf(Normal([1, x]'beta - y, sigma), 0.0))