DM.MoM {HMP} | R Documentation |
Method-of-Moments (MoM) Estimators of the Dirichlet-Multinomial Parameters
Description
Method-of-Moments (MoM) estimators of the Dirichlet-multinomial parameters: taxa proportions and overdispersion.
Usage
DM.MoM(data)
Arguments
data |
A matrix of taxonomic counts(columns) for each sample(rows). |
Details
Given a set of taxa-count vectors , the methods of moments (MoM) estimator of the set of parameters
and
is given as follows (Mosimann, 1962; Tvedebrink, 2010):
and
where ,
and
,
and
with
.
Value
A list providing the MoM estimator for overdispersion, the MoM estimator of the RAD-probability mean vector, and the corresponding loglikelihood value for the given data set and estimated parameters.
References
Mosimann, J. E. (1962). On the compound multinomial distribution, the multivariate -distribution, and correlations among proportions. Biometrika 49, 65-82.
Tvedebrink, T. (2010). Overdispersion in allelic counts and theta-correction in forensic genetics. Theor Popul Biol 78, 200-210.
Examples
data(throat)
fit.throat <- DM.MoM(throat)
fit.throat