drop_model {simDNAmixtures} | R Documentation |
Defines a (semi-continuous) drop model
Description
Defines a (semi-continuous) drop model
Usage
drop_model(dropout_probabilities, drop_in_rate = 0, freqs, model_settings)
Arguments
dropout_probabilities |
Numeric vector with values between 0 and 1. Dropout probabilities for each contributor. |
drop_in_rate |
Numeric vector of length one. Expected number of drop-ins per locus. Default is 0. |
freqs |
Optionally a list with allele frequencies (needed when drop_in_rate > 0). See read_allele_freqs. |
model_settings |
List. Possible parameters:
|
Details
Define the classes semi-continuous drop-model. The model may then be used to sample DNA profiles using the sample_mixture_from_genotypes function. Alternatively, to sample many models and profiles in one go with parameters according to a specified distribution, the sample_mixtures function can be used.
Value
Object of class pg_model
.
References
Slooten, K. (2017). Accurate assessment of the weight of evidence for DNA mixtures by integrating the likelihood ratio. Forensic Science International: Genetics, 27, 1-16. doi:10.1016/j.fsigen.2016.11.001
See Also
sample_mixture_from_genotypes, sample_mixtures, gamma_model, log_normal_model.
Examples
gf <- gf_configuration()
freqs <- read_allele_freqs(system.file("extdata","FBI_extended_Cauc_022024.csv",
package = "simDNAmixtures"))
settings <- list(locus_names = gf$autosomal_markers, size_regression = gf$size_regression)
model <- drop_model(dropout_probabilities = c(0.1),
drop_in_rate = 1e-3,
freqs = freqs, model_settings = settings)
g <- sample_contributor_genotypes(contributors = "U1", freqs = freqs,
loci = settings$locus_names)
# genotype
g
# sample with dropout
sample_mixture_from_genotypes(g, model)