coin-toss-in-stan.R

# load packages
library(rstan)
library(bayesplot)

# make a stan model
fancy_stan_model <- "
generated quantities {
    int<lower=0, upper=1> y_sim = bernoulli_rng(0.5); 
}
"

# compile the Stan model
compiled_fancy_stan_model <- stan_model(model_code = fancy_stan_model)

# fit the model
fit <- sampling(compiled_fancy_stan_model, 
                iter = 2000, 
                warmup = 500, 
                chains = 4, 
                seed = 123, 
                algorithm = "Fixed_param")

# summarize the results
posterior <- extract(fit)$y_sim
table(posterior)