# NOT RUN {
library("MASS")
n <- 2
params <- params_lm(coefs = MASS::mvrnorm(n, mu = c(.5,.6),
Sigma = matrix(c(.05, .01, .01, .05), nrow = 2)),
sigma <- rgamma(n, shape = .5, rate = 4))
print(params)
# }
Run the code above in your browser using DataLab