f1dr | 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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.