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

Description

Exact numerical calculation

Usage

1
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 Thore.Egeland@nmbu.com

References

Slooten and Egeland (2013, submitted)

Examples

1
2
3
4
5
6
7
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)