zygdist: Zygote dosage probabilities.

In dcgerard/hwep: Hardy-Weinberg Equilibrium in Polyploids

Zygote dosage probabilities.

Description

Calculates the distribution of an offspring dosages given parental dosages (G1 and G2), the ploidy of the species (ploidy), and the double reduction parameter (alpha).

Usage

zygdist(alpha, G1, G2, ploidy)

Arguments

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.
G1	The dosage of parent 1. Should be an integer between 0 and ploidy.
G2	The dosage of parent 2. Should be an integer between 0 and ploidy.
ploidy	The ploidy of the species. This should be an even positive integer.

Value

A vector of probabilities. The ith element is the probability that the offspring will have dosage i-1.

Author(s)

David Gerard

Examples

zygdist(alpha = c(0.5, 0.1), G1 = 4, G2 = 5, ploidy = 8)

dcgerard/hwep documentation built on Sept. 13, 2023, 8:49 a.m.