# \donttest{
if (requireNamespace("rstanarm", quietly = TRUE)) {
set.seed(123)
years <- 2000:2010
direct <- data.frame(
Year = years,
A = 100 + cumsum(rnorm(11)),
B = 120 + cumsum(rnorm(11))
)
production <- data.frame(
Year = years,
A = 102 + cumsum(rnorm(11)),
B = 118 + cumsum(rnorm(11))
)
results <- fit_bayesian_glm_sectors(direct, production,
chains = 2, iter = 1000)
coefs <- extract_sector_coefficients(results$results)
print(coefs)
}
# }
Run the code above in your browser using DataLab