Compositional: diri.est – R documentation

Become an expert in R — Interactive courses, Cheat Sheets, certificates and more!

diri.est

Fitting a Dirichlet distribution

Estimation of the parameters of a fitted Dirichlet distribution.

diri.est(x, type = "mle")

`x`	A matrix containing compositional data.
`type`	If you want to estimate the parameters use type="mle". If you want to estimate the mean vector along with the precision parameter, the second parametrisation of the Dirichlet, use type="prec".

Maximum likelihood estimation of the parameters of a Dirichlet distribution is performed.

A list including:

`loglik`	The value of the log-likelihood.
`param`	The estimated parameters.
`phi`	The estimated precision parameter, if type = "prec".
`a`	The estimated mean vector, if type = "prec".
`runtime`	The run time of the maximisation procedure.

Michail Tsagris

R implementation and documentation: Michail Tsagris mtsagris@uoc.gr and Giorgos Athineou <gioathineou@gmail.com>.

Ng Kai Wang, Guo-Liang Tian and Man-Lai Tang (2011). Dirichlet and related distributions: Theory, methods and applications. John Wiley \& Sons.

Aitchison J. (1986). The statistical analysis of compositional data. Chapman \& Hall.

x <- rdiri( 100, c(5, 7, 1, 3, 10, 2, 4) )
diri.est(x)
diri.est(x, type = "prec")

Compositional Data Analysis

v4.6

GPL (>= 2)

Authors

Michail Tsagris [aut, cre], Giorgos Athineou [aut], Abdulaziz Alenazi [ctb]

Initial release

2021-04-27