implements a simple true-genotype-independent genotyping error model
ge_model_TGIE(L, epsilon = 0.01)
Arguments
- L
required locus specific information
- epsilon
the rate at which genotypes are incorrectly observed.
Examples
L <- list(
freqs = c(A = 0.6, B = 0.4),
geno_freqs = c(`A-A` = 0.36, `A-B` = 0.48, `B-B` = 0.16)
)
ge_model_TGIE(L, epsilon = 0.02)
#> A-A A-B B-B
#> A-A 0.980000000 0.01500000 0.005000000
#> A-B 0.013846154 0.98000000 0.006153846
#> B-B 0.008571429 0.01142857 0.980000000