Sample from the posterior predictive of a 1-d binary population (very quickly)

mcmc_binary_superfast(obs, N, nsamples, a, b)

Arguments

obs

numeric vector representing the observed sample.

N

size of the population.

nsamples

number of samples to draw

a

first hyperparameter for beta prior

b

second hyperparameter for beta prior