dPD: The Poisson-Dirichlet distribution

In PEkit: Partition Exchangeability Toolkit

Description

Distribution function for the Poisson-Dirichlet distribution.

Usage

dPD(abund, psi)

Arguments

abund	An abundance vector.
psi	Dispersal parameter ψ. Accepted input values are positive real numbers, "a" for absolute value ψ=1 by default, or "r" for relative value ψ=n, where n is the size of the input sample.

Details

Given an abundance vector abunds, calculates the probability of a data vector x given by the Poisson-Dirichlet distribution. The higher the dispersal parameter ψ, the higher the amount of distinct observed species. In terms of the paintbox process, a high ψ increases the size of the continuous part p_0 of the process, while a low ψ will increase the size of the discrete parts p_{\neq 0}.

Value

The probability of the Poisson-Dirichlet distribution for the input abundance vector, e.g. an exchangeable random partition, and a dispersal parameter ψ.

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

## Get a random sample from the Poisson Dirichlet distribution, and
## find the probability of such a sample with psi=5:
set.seed(111)
s <- rPD(n=100,psi=5)
a=abundance(s)
dPD(a, psi=5)

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