Description Usage Arguments Details Value Author(s) References See Also Examples

Fitting a Dirichlet distribution via Newton-Rapshon.

1 |

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

`type` |
Type 1 uses a vectorised version of the Newton-Raphson (Minka, 2012). In high dimensions this is to be preferred. If the data are too concentrated, regardless of the dimensions, this is also to be preferrred. Type 2 uses the regular Newton-Raphson, with matrix multiplications. In small dimensions this can be considerably faster. |

`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:

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

`param` |
The estimated parameters. |

Michail Tsagris and Manos Papadakis

R implementation and documentation: Michail Tsagris <[email protected]> and Manos Papadakis <[email protected]>

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 |

Rfast documentation built on Nov. 18, 2017, 4 a.m.

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.