f1dr {hwep} | R Documentation |
Estimates double reduction in F1 populations by maximum likelihood.
f1dr(nvec, G1, G2)
nvec |
A vector containing the observed genotype counts,
where |
G1 |
The dosage of parent 1. Should be an integer between |
G2 |
The dosage of parent 2. Should be an integer between |
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.
David Gerard
zygdist()
for calculating the probability of
offpring genotypes given parental genotypes and the double reduction
rate.
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)