library(rjags)
read_rdump <- rstan::read_rdump
set.seed(1234)
line_data <- read_rdump("line_data.R")
line_model <- jags.model("line.jag", data=line_data, n.adapt = 0, n.chains = 2)
line_mcmc_list <-
coda.samples(line_model, c("beta", "tau"), n.iter = 100, thin = 1)
save(line_mcmc_list, file="line_mcmc_list.Rdata")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.