Goodness-of-fit test comparing data following a uniform distribution to a Hardy-Weinberg distribution.

GOF_uniform determines the GOF of the heterozygous samples to the HWE model.

GOF_uniform(ref_counts, var_counts, allelefreq, SE, inbr = 0)

Arguments

ref_counts

Numeric list. Reference counts.

var_counts

Numeric list. Variant counts.

allelefreq

Number. Allele frequency.

SE

Number. Sequencing error rate.

inbr

Number. Degree of inbreeding (default = 0).

Value

Geometric mean of the ratios of the putative heterozygous samples as a measure of the GOF.

Examples

GOF_uniform(c(5, 8, 10, 3, 5, 6, 23), c(8, 8, 6, 4, 4, 10, 0), 0.5, 0.002, 0.12)
#> [1] 0.3810761
GOF_uniform(c(5, 0, 0, 3, 5, 1, 23), c(1, 8, 6, 2, 0, 10, 0), 0.5, 0.002)
#> [1] 0.9689088