# NOT RUN {
n <- 2
coefs1 <- MASS::mvrnorm(n, mu = c(.5,.6),
Sigma = matrix(c(.05, .01, .01, .05), nrow = 2))
sigma1 <- rgamma(n, shape = .5, rate = 4)
params1 <- params_lm(coefs = coefs1, sigma = sigma1)
coefs2 <- MASS::mvrnorm(n, mu = c(.2,.9),
Sigma = matrix(c(.08, .02, .02, .08), nrow = 2))
sigma2 <- rgamma(n, shape = .9, rate = 4)
params2 <- params_lm(coefs = coefs2, sigma = sigma2)
params_list <- params_lm_list(params1, params2)
print(params_list)
# }
Run the code above in your browser using DataLab