/introduction_to_bayes_test_script.R
Last active Feb 12, 2017
Test script for "Introduction to Bayesian data analysis with 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) |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment