MLEp: Maximum Likelihood Estimate of psi

In PEkit: Partition Exchangeability Toolkit

Description

Numerically searches for the MLE of ψ given an abundance vector with a binary search algorithm.

Usage

MLEp(abund)

Arguments

abund

An abundance vector.

Details

Numerically searches for the MLE of ψ as the root of equation

K=∑_{i=1}^nψ/(ψ+i-1),

where K is the observed number of different species in the sample. The right side of the equation is monotonically increasing when ψ>0, so a binary search is used to find the root. An accepted ψ sets value of the right side of the equation within R's smallest possible value of the actual value of K.

Value

The MLE of ψ.

References

W.J. Ewens, The sampling theory of selectively neutral alleles, Theoretical Population Biology, Volume 3, Issue 1, 1972, Pages 87-112, ISSN 0040-5809, <doi: 10.1016/0040-5809(72)90035-4>.

Examples

##Find the MLE of psi of the vector (1,2,2).
##The frequencies of the frequencies of the data vector are given as input:
MLEp(abundance(c(1,2,2)))

##Find the MLE of psi of a sample from the Poisson-Dirichlet distribution:
set.seed(1000)
x<-rPD(n=10000, psi=100)
MLEp(abundance(x))

PEkit documentation built on Nov. 22, 2021, 9:08 a.m.