# Generate data
[Link](123)
x <- rnorm(50, mean = 5, sd = 2)
# Known variance
sigma2 <- 4 # (since sd = 2)
# Prior
mu0 <- 5
tau2 <- 1
# Sample stats
n <- length(x)
xbar <- mean(x)
# Bayes estimate of mu (SELF)
mu_bayes <- ((n / sigma2) * xbar + (1 / tau2) * mu0) / ((n /
sigma2) + (1 / tau2))
mu_bayes