dgamete: Gamete dosage probability

In hwep: Hardy-Weinberg Equilibrium in Polyploids

Gamete dosage probability

Description

Estimates the probability of a gamete dosage given the parent dosage (G), the parent ploidy (ploidy), and the double reduction parameter (alpha). This is for biallelic loci.

Usage

dgamete(x, alpha, G, ploidy, log_p = FALSE)

Arguments

x	A vector of numerics in seq(0, ploidy/2). The dosage of the gametes.
alpha	A numeric vector containing the double reduction parameter(s). This should be a vector of length floor(ploidy/4) where alpha[i] is the probability of exactly i pairs of IBDR alleles being in the gamete. Note that sum(alpha) should be less than 1, as 1 - sum(alpha) is the probability of no double reduction.
G	The dosage of the parent. Should be an integer between 0 and ploidy.
ploidy	The ploidy of the species. This should be an even positive integer.
log_p	A logical. Should we return the log-probability (TRUE) or not (FALSE)? Defaults to FALSE.

Value

A vector of length length(x), containing the (log) probabilities of a gamete carrying a dosage of x from a parent of dosage G who has ploidy ploidy and a double reduction rate alpha.

Author(s)

David Gerard

Examples

dgamete(x = 0:2, alpha = 0, G = 2, ploidy = 4)

hwep documentation built on May 31, 2023, 9:06 p.m.