powered by
Provide a vector of genotype counts and this will return a matrix of genotype log-likelihoods.
simgl(nvec, rd = 10, seq = 0.01, bias = 1, od = 0.01)
A matrix of genotype log-likelihoods. The rows index the individuals and the columns index the genotypes. This is natural log (base e).
A vector of counts. nvec[k] is the number of folks with a genotype of k-1.
nvec[k]
Read depth. Lower is more uncertain.
Sequencing error rate. Higher means more uncertain.
Bias. 1 means no bias.
Overdispersion. Typical value is like 0.01. Higher means more uncertain.
David Gerard
set.seed(1) simgl(nvec = c(1, 2, 1, 1, 3))
Run the code above in your browser using DataLab