Estimate Double Reduction in F1 Populations

Estimates double reduction in F1 populations by maximum likelihood.

f1dr(nvec, G1, G2)

Arguments

nvec

A vector containing the observed genotype counts, where nvec[[i]] is the number of individuals with genotype i-1. This should be of length ploidy+1.

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.

Value

A list with some or all of the following elements:

alpha

A vector of numerics of length floor(ploidy / 4), the estimated double reduction rate.

llike

The final log-likelihood.

See also

zygdist() for calculating the probability of offpring genotypes given parental genotypes and the double reduction rate.

Author

David Gerard

Examples

set.seed(1)
size <- 100
qvec <- zygdist(alpha = 0.1, G1 = 2, G2 = 2, ploidy = 4)
nvec <- c(stats::rmultinom(n = 1, size = size, prob = qvec))
f1dr(nvec = nvec, G1 = 2, G2 = 2)
#> $alpha
#> [1] 0.1352307
#> 
#> $llike
#> [1] -8.972345
#>