fit.dirichlet: Estimate Dirichlet Parameters
In dirichlet: Dirichlet and Generalized Dirichlet Distribution functions

Description Usage Arguments Value Author(s) See Also Examples

Estimate Dirichlet and Generalized Dirichlet distribution parameters.

1 2	fit.dirichlet(X, type) fit.genDirichlet(X, type)

`X`	Matrix with columns corresponding to the entries of `a` and rows corresponding to the number of simulations in the hypercube.
`type`	`mm` for Method of Moments estimation or `ml` for Maximum Likelihood estimation.

A list giving the k value or vector of k values for the Generalized Dirichlet and vector of p estimates. The normal Dirichlet parameterization can be obtained by k*p. The fit.dirichlet function returns two estimates of k. One based on the most likely parameter in the distribution, the other based on a weighted mean of the k estimates for each parameter.

Rob Carnell

rdirichlet, rGenDirichlet

  Z <- rdirichlet(100, c(4,3,2))
  fit.dirichlet(Z)
  fit.dirichlet(Z, "ml")
  
  Z2 <- rGenDirichlet(100, c(0.7, 0.2, 0.1), c(10, 20, 20))
  fit.genDirichlet(Z2)