R/gibbs_betabin.R

In bayesball/TeachBayes: Teaching Bayesian Inference

Defines functions gibbs_betabin

Documented in gibbs_betabin

gibbs_betabin <- function(n, a, b, p = 0.5, iter = 1000){
  x <- matrix(0, iter, 2)
  for(k in 1:iter){
    y <- rbinom(1, size = n, prob = p )
    p <- rbeta(1, y + a, n - y + b )
    x[k, ] <- c(y, p)
  }
  x
}