euroMix (version 1.1.1)

LRmoments: Calculates expectation, standard deviation and skewness of LR under HP and HD

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

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)