Fitting a Dirichlet distribution via Newton-Rapshon.

1 | ```
diri.nr2(x, tol = 1e-07)
``` |

`x` |
A matrix containing the compositional data. Zeros are not allowed. |

`tol` |
The tolerance level idicating no further increase in the log-likelihood. |

Maximum likelihood estimation of the parameters of a Dirichlet distribution is performed via Newton-Raphson. Initial values suggested by Minka (2012) are used.

A list including:

`iters` |
The number of iterations required. |

`loglik` |
The value of the log-likelihood. |

`param` |
The estimated parameters. |

Michail Tsagris and Manos Papadakis

R implementation and documentation: Michail Tsagris <mtsagris@yahoo.gr> and Manos Papadakis <papadakm95@gmail.com>

Minka Thomas (2012). Estimating a Dirichlet distribution. Technical report.

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

1 2 3 |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.