rdirichlet: Draw a single sample from a Dirichlet
In EpiNow2: Estimate and Forecast Real-Time Infection Dynamics
rdirichlet R Documentation
Draw a single sample from a Dirichlet
Description
Base R does not provide an rdirichlet(). We use the
gamma-normalisation method also used by the Stan model:
draw an independent Gamma(alpha_i, 1) per bin and rescale by
the segment sum. Bins with alpha == 0 stay at zero so
structural zeros (e.g. the t = 0 generation-time bin) are
preserved.
Usage
rdirichlet(alpha)
Arguments
alpha
A non-negative numeric vector of concentration
parameters.
Value
A numeric vector the same length as alpha, summing
to 1 over the positive-alpha entries.
References
Stan discourse, "Ragged array of simplexes",
https://discourse.mc-stan.org/t/ragged-array-of-simplexes/1382/21.
EpiNow2 documentation built on June 17, 2026, 1:07 a.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
| rdirichlet | R Documentation |
Draw a single sample from a Dirichlet
Description
Base R does not provide an rdirichlet(). We use the
gamma-normalisation method also used by the Stan model:
draw an independent Gamma(alpha_i, 1) per bin and rescale by
the segment sum. Bins with alpha == 0 stay at zero so
structural zeros (e.g. the t = 0 generation-time bin) are
preserved.
Usage
rdirichlet(alpha)
Arguments
alpha |
A non-negative numeric vector of concentration parameters. |
Value
A numeric vector the same length as alpha, summing
to 1 over the positive-alpha entries.
References
Stan discourse, "Ragged array of simplexes", https://discourse.mc-stan.org/t/ragged-array-of-simplexes/1382/21.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.