Random generation of Dirichlet-Multinomial samples.
Usage
Dirichlet.multinomial(Nrs, shape)
Arguments
- Nrs
A vector specifying the number of reads or sequence depth for each sample.
- shape
A vector of Dirichlet parameters for each taxa.
Details
Adapted from HMP::Dirichlet.multinomial(), since HMP was removed from CRAN