Compositional (version 1.0)

diri.est: Fiiting a Dirichlet distribution

Description

Estimation of the parameters of a fitted Dirichlet distribution.

Usage

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

Arguments

x

A matrix containing the 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". If you want to estimate the parameters via the entropy, the

Value

• 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".

Details

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

References

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.

See Also

diri.nr, diri.contour

Examples

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