HWprob: Hardy-Weinberg probabilities

Description Usage Arguments Value Examples

View source: R/utils.R

Description

Hardy-Weinberg probabilities

Usage

1
HWprob(allele1, allele2, afreq, f = 0)

Arguments

allele1, allele2

Vectors of equal length, containing alleles in the form of indices of afreq

afreq

A numeric vector with allele frequencies

f

A single number in [0, 1]; the inbreeding coefficient

Value

A numeric vector of the same length as allele1 and allele2

Examples

1
2
3
p = 0.1; q = 1-p
hw = HWprob(c(1,1,2), c(1,2,2), c(p, q))
stopifnot(all.equal(hw, c(p^2, 2*p*q, q^2)))

pedprobr documentation built on Nov. 13, 2020, 5:06 p.m.