Skip to contents

Genotype frequencies of an F1 population under a generalized model.

Source: R/higher_freq.R
gf_freq.Rd

Genotype frequencies of an F1 population under a generalized model.

Usage

gf_freq(
  p1_g,
  p1_ploidy,
  p1_gamma = NULL,
  p1_alpha = NULL,
  p1_beta = NULL,
  p1_type = c("mix", "polysomic"),
  p1_add_dr = TRUE,
  p2_g,
  p2_ploidy,
  p2_gamma = NULL,
  p2_alpha = NULL,
  p2_beta = NULL,
  p2_type = c("mix", "polysomic"),
  p2_add_dr = TRUE,
  pi = 0,
  nudge = sqrt(.Machine$double.eps)
)

Arguments

p1_g, p1_ploidy, p1_gamma, p1_alpha, p1_beta, p1_type, p1_add_dr

The first parent's version of the parameters in gamfreq().

p2_g, p2_ploidy, p2_gamma, p2_alpha, p2_beta, p2_type, p2_add_dr

The second parent's version of the parameters in gamfreq().

pi

The proportion of outliers.

nudge

Zeros go to nudge

Value

A vector of genotype frequencies. Element i is the probability and offspring has genotype i - 1.

See also

gamfreq().

Author

David Gerard

Examples

q <- gf_freq(
  p1_g = 2,
  p1_ploidy = 4,
  p1_gamma = c(0.1, 0.9),
  p1_type = "mix",
  p2_g = 2,
  p2_ploidy = 6,
  p2_alpha = 0.1,
  p2_type = "polysomic",
  pi = 0.05)