drop_model: Defines a (semi-continuous) drop model

In simDNAmixtures: Simulate Forensic DNA Mixtures

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: locus_names. Character vector. size_regression. Function, see read_size_regression.

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. \Sexpr[results=rd]{tools:::Rd_expr_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)

simDNAmixtures documentation built on April 15, 2025, 1:11 a.m.