recordHeights: Records the peak heights of the alleles present in the analyzed stains

Description

The recordHeights function records the peak heights of the alleles present in the analyzed stains. 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

recordHeights(x,y=NULL,geno,tabcsv,byloc=FALSE)

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.

byloc

logical indicating whether data should be displayed per locus (TRUE) or overall loci (FALSE, default)

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

Gill P, Puch-Solis R, Curran J. The low-template-DNA (stochastic) threshold-Its determination relative to risk analysis for national DNA databases. Forensic Science International: Genetics, 2009, 3, 104-111

Examples

Run this code

#load the exemple data
data(dropdata)
tabcsv<-dropdata$tabcsv
genot<-dropdata$genot
#individuals' labels are 1 and 2
#peak heights of heterozygote genotypes with non shared alleles
recordHeights(1,2,geno=genot,tabcsv=tabcsv)