powered by
Estimates double reduction in F1 populations by maximum likelihood.
f1dr(nvec, G1, G2)
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.
floor(ploidy / 4)
llike
The final log-likelihood.
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.
nvec[[i]]
i-1
ploidy+1
The dosage of parent 1. Should be an integer between 0 and ploidy.
0
ploidy
The dosage of parent 2. Should be an integer between 0 and ploidy.
David Gerard
zygdist() for calculating the probability of offpring genotypes given parental genotypes and the double reduction rate.
zygdist()
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)
Run the code above in your browser using DataLab