library(dplyr)
library(bayesewi)
# This block loads in the data, and fits EWI models to each
da_data = read.csv("/users/eric.ward/downloads/da.examples.csv")
da_data = select(da_data, example, dec.date, DA.ppm) %>%
rename(y = DA.ppm, x = dec.date) %>%
filter(y != 0)
da_data$y = log(da_data$y)
models = list()
for(i in unique(da_data$example)) {
df = filter(da_data, example==i) %>%
select(-example)
df$x = round(365*(df$x - min(df$x)))
models[[i]] = fit_ewi(df, ewi_model="ar", iter = 2000, chains=3)
}
saveRDS(models, file="fitted_ar_models.rds")
pdf("DA_ewi_ar.pdf")
for(i in unique(da_data$example)) {
plot_estimates(models[[i]], alpha = 0.05, link_space=TRUE)
}
dev.off()
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.