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

Description Usage Arguments Value Author(s) References Examples

View source: R/LRmoments.R

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 [email protected]

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)

euroMix documentation built on May 31, 2017, 4:51 a.m.