recordDrop: Records the allelic dropout events matched with individual DNA proxies

Description

The recordDrop function records the dropout events from experimental data. The function aims to facilitate the manipulation of experimental data used for the estimation of the probability of allelic dropout (cf. the references sections).

Usage

recordDrop(x, y, geno, tabcsv,s=40)

Arguments

numeric label of the contributing individual, if the stain is a mixture, x should give the label of the first individual contributing to the mixture

numeric label of the second contributing individual, default is NULL. If the stain is a mixture, y should give the label of the second individual contribution to the mixture. This argument is skipped if the stain is not a mixture (default

geno

a matrix giving the genotypes of the individuals contributing to the analyzed data for each locus. An individual genotype is given in rows for each locus in column. A homozygous carrying allele 9 is coded '9/9', a heterozygous carrying alleles 8 and 9

tabcsv

a matrix giving the validation table of the analysed DNA stain. tabcsv must have a ``genemapper'' validation table structure, namely, information about the present alleles and the corresponding peak heights must be given.

numeric giving the detection threshold for alleles in Relative fluorescence units (RFU) , default is set to 40 RFUS. An observed allele with a peak height smaller (

Value

A list of length the number of analyzed loci, each component of the list is a matrix with the following information: - The names of expected alleles - The expected allele counts for the first contributor (when date is a mixture) - The expected allele counts for the second contributor (when date is a mixture) - The observed alleles - The observed peak heights - The dropout variable D, takes 1 if the allele has dropped out, 0 otherwise

References

Tvedebrink T, Eriksen PS, Mogensen HS, Morling N. Estimating the probability of allelic drop-out of STR alleles in forensic genetics. Forensic Science International: Genetics, 2009, 3(4), 222-226.

Examples

Run this code

#load the exemple data
data(dropdata)
tabcsv<-dropdata$tabcsv
genot<-dropdata$genot
#individuals' labels are 1 and 2
#record the dropout the surviving peak heights for heterozygotes whith non shared alleles
recordDrop(1,2,geno=genot,tabcsv=tabcsv,s=40)