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

Description Usage Arguments Value Author(s) References Examples

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@gmail.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)


Search within the euroMix package
Search all R packages, documentation and source code

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.