LRmoments: Calculates expectation, standard deviation and skewness of LR...

In euroMix: Calculations for DNA mixtures

Description

Exact numerical calculation

Usage

LRmoments(p = c(0.5, 0.5), kappaP = c(0, 1, 0),kappaD = c(1, 0, 0),log10=FALSE)

Arguments

p	Allele frequencies
kappaP	Probabilities of 0,1 and 2 alleles IBD corresponding to pedigree for HP
kappaD	Probabilities of 0,1 and 2 alleles IBD corresponding to pedigree for HD
log10	If TRUE, LR is log10 transformed.

Value

moments	expectation, standard deviation and skewness of LR under HP and HD
LRtable	Distribution of LR under HP and HD

Author(s)

Thore Egeland [e-mail]

References

Slooten and Egeland (2013, submitted)

Examples

LRmoments(kappaP=c(0,1,0)) #Motivating example
LRmoments(kappaP=c(0,0.25,0.75)) #skew(LR(HP))<0
#Appendix of Slooten and Egeland (2013, submitted)
## Not run: data(db)
p=db[db$Marker=="VWA",]$Freq
LRmoments(p=p,kappaP=c(0,1,0))
## End(Not run)

euroMix documentation built on May 2, 2019, 6:27 p.m.