global_stutter_model {simDNAmixtures}

Global stutter model where the expected stutter rate is constant across alleles and loci

Description

Global stutter model where the expected stutter rate is constant across alleles and loci

Usage

global_stutter_model(back_stutter_rate, forward_stutter_rate, size_regression)

Arguments

back_stutter_rate

Numeric. (Optional)

forward_stutter_rate

Numeric. (Optional)

size_regression

Function, see read_size_regression.

Details

When a pg_model is constructed (see gamma_model), a stutter model can optionally be applied. In the global stutter model, the expected stutter rate is constant across all loci and for all parent alleles.

Value

Object of class stutter_model to be used by e.g. gamma_model.

See Also

allele_specific_stutter_model for a stutter model where the expected stutter rate depends on the allele and locus.

Examples

# the stutter model needs a size regression to determine fragment length
# of stutter products
size_regression <- read_size_regression(system.file("extdata",
"GlobalFiler_SizeRegression.csv",package = "simDNAmixtures"))

# define a stutter model with an expected back stutter rate of 10%
stutter_model <- global_stutter_model(back_stutter_rate = 0.1,
                                     size_regression = size_regression)

stutter_model