dlevene {HardyWeinberg}

Calculate Levene's exact density for k alleles

Description

Function dlevene calculates Levene's exact density for a diploid system with k alleles.

Usage

dlevene(N)

Arguments

N

A lower triangular matrix with genotype counts

Details

The supplied matrix of genotype counts should be triangular, with the homozygote counts on the diagonal, and all heterozygote counts below the diagonal.

Value

a single real number

Author(s)

Jan Graffelman ([e-mail])

References

Levene, H. (1949) On a matching problem arising in genetics. Annals of Mathematical Statistics, 20, pp. 91-94.

See Also

HWExact

Examples

x <- c(AA=12,AB=19,AC=13,BB=7,BC=5,CC=0)
x <- toTriangular(x)
prob <- dlevene(x)
print(prob)