Test script for "Introduction to Bayesian data analysis with R"

introduction_to_bayes_test_script.R

# Prior to the tutorial make sure that the script below runs without error on your R installation.
# What you need is a working installation of JAGS: http://mcmc-jags.sourceforge.net/
# and the rjags R package that can be installed from R by running
# install.packages("rjags")

# Generating some fake data
set.seed(123)
y <- rbinom(30, size = 1, prob = 0.2015)

# Fitting a simple binomial model using JAGS
library(rjags)

model_string <- "model {
for(i in 1:length(y)) {
  y[i] ~ dbern(theta)
}

theta ~ dunif(0, 1)
}"

jags_model <- jags.model(textConnection(model_string), data = list(y = y))
mcmc_samples <- coda.samples(jags_model, variable.names = "theta", n.iter = 1000)
plot(mcmc_samples)
summary(mcmc_samples)