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

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

mcmc_multinomial_superfast(obs, N, nsamples, alpha)

Arguments

obs

numeric vector representing the observed sample.

N

size of the population.

nsamples

number of samples to draw

alpha

hyperparameter for dirichlet prior; can be either a single number or a numeric vector