Contour plot of the Flexible Dirichlet distribution in S^2
Contour plot of the Flexible Dirichlet distribution in S^2.
fd.contour(alpha, prob, tau, n = 100, x = NULL)
alpha |
A vector of the non-negative alpha parameters. |
prob |
A vector of the clusters' probabilities. It must sum to one. |
tau |
The non-negative scalar tau parameter. |
n |
The number of grid points to consider over which the density is calculated. |
x |
This is either NULL (no data) or contains a 3 column matrix with compositional data. |
The user can plot only the contour lines of a Dirichlet with a given vector of parameters, or can also add the relevant data should they wish to.
A ternary diagram with the points and the Flexible Dirichlet contour lines.
Michail Tsagris.
R implementation and documentation: Michail Tsagris mtsagris@uoc.gr.
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. (2017). A structured Dirichlet mixture model for compositional data: inferential and applicative issues. Statistics and Computing, 27, 963–983.
## Not run: x <- rfd(n = 50, a = c(12, 11, 10), p = c(0.25, 0.25, 0.5), tau = 4 ) mod <- fd.est(x, ini.iter = 10, final.iter = 20) fd.contour(mod$alpha, mod$prob, mod$tau, x = x) ## End(Not run)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.