FD.generate {FlexDir} | R Documentation |
The Flexible Dirichlet Random Generation
Description
Random generation from the Flexible Dirichlet distribution with parameters a
, p
and t
.
Usage
FD.generate(n, a, p, t)
Arguments
n |
number of points on the simplex to be generated. |
a |
vector of the non-negative alpha parameters. |
p |
vector of the clusters' probabilities. It must sum to one. |
t |
non-negative scalar tau parameter. |
Details
Vectors a
and p
must be of the same length.
The Flexible Dirichlet distribution derives from the normalization of a basis of positive dependent random variables obtained by starting from a basis of independent equally scaled gamma random variables, and randomly allocating to the i
-th element a further independent gamma random variable.
References
Ongaro, A. and Migliorati, S. (2013) A generalization of the Dirichlet distribution. Journal of Multivariate Analysis, 114, 412–426.
Migliorati, S., Ongaro, A. and Monti, G. S. (2016) A structured Dirichlet mixture model for compositional data: inferential and applicative issues. Statistics and Computing, 1–21.
See Also
FD.estimation
, FD.density
, FD.theorcontours
, FD.subcomposition
, FD.amalgamation
Examples
n <- 100
alpha <- c(12,7,15)
prob <- c(0.3,0.4,0.3)
tau <- 8
data <- FD.generate(n,alpha,prob,tau)
data