rPD: Random sampling from the Poisson-Dirichlet Distribution
In ville-kinnula/SPEC: Supervised Classification under Partition Exchangeability

Description Usage Arguments Details Value References Examples

View source: R/PD_distribution.R

rPD samples randomly from the PD distribution with a given ψ by simulating the Hoppe urn model.

1	rPD(n, psi)

`n`	number of observations.
`psi`	dispersal parameter.

rPD samples random values with a given ψ from the Poisson-Dirichlet distribution by simulating the Hoppe urn model.

rPD returns a vector with a sample of size n from the Hoppe urn model with parameter ψ, along with its table of frequencies. given parameter ψ.

Hoppe, F.M. The sampling theory of neutral alleles and an urn model in population genetics. J. Math. Biology 25, 123–159 (1987). https://doi.org/10.1007/BF00276386.

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

## Get random sample from the PD distribution with different psi,
## and estimate the psi of the samples:
s1<-rPD(1000, 10)
s2<- rPD(1000, 50)
print(c(MLEp(table(table(s1)))$psi, MLEp(table(table(s2)))$psi))