# NOT RUN {
library(dplyr)
library(yardstick)
# Compare model predictions to observation
predictions <- focal_vs_comp_ex %>%
mutate(growth_hat = predict(comp_bayes_lm_ex, focal_vs_comp_ex))
predictions %>%
rmse(truth = growth, estimate = growth_hat) %>%
pull(.estimate)
# Plot posterior parameters
comp_bayes_lm_ex %>%
autoplot()
# }
Run the code above in your browser using DataLab