Multinomial Dirichlet Sampler
MultDirichletResample(Y, alpha = rep(0.65, nrow(Y)), n_samples = 2000)
| Y | count data (D x N) (e.g., taxa x samples) |
|---|---|
| alpha | D-vector of Dirichlet Prior parameters |
| n_samples | number of posterior samples to collect |
Posterior samples as array of dimensions: D x N x n_samples