Fitting a Flexible Dirichlet distribution
Fitting a Flexible Dirichlet distribution.
fd.est(x, ini.iter = 50, final.iter = 100)
x |
A matrix or a dataframe containing the compositional data. |
ini.iter |
Number of iterations for the initial SEM step. Default value is 50. |
final.iter |
Number of iterations for the final EM step. Default value is 100. |
For more information see the references.
A list including:
alpha |
Estimated values of the parameter vector alpha. |
prob |
Estimated values of the parameter vector p. |
tau |
Estimated value of the parameter tau. |
loglik |
The Log-likelihood value. |
Michail Tsagris ported from the R package FlexDir. 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 ) ela <- fd.est(x, ini.iter = 10, final.iter = 20) ela ## End(Not run)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.