Spark Online Training by Edureka

LD22: LD statistics for two diallelic markers

Description Usage Arguments Value References Note Author(s) See Also Examples

Description

LD statistics for two SNPs.

It is possible to perform permutation test of r^2 by re-ordering the genotype through R's sample function, obtaining the haplotype frequencies by gc.em or genecounting, supplying the estimated haplotype frequencies to the current function and record x2, and comparing the observed x2 and that from the replicates.

Usage

1
LD22(h,n)

Arguments

h

a vector of haplotype frequencies

n

number of haplotypes

Value

The returned value is a list containing:

h

the original haplotype frequency vector

n

the number of haplotypes

D

the linkage disequilibrium parameter

VarD

the variance of D

Dmax

the maximum of D

VarDmax

the variance of Dmax

Dprime

the scaled disequilibrium parameter

VarDprime

the variance of Dprime

x2

the Chi-squared statistic

lor

the log(OR) statistic

vlor

the var[log(OR)] statistic

References

Zabetian CP, Buxbaum SG, Elston RC, Kohnke MD, Anderson GM, Gelernter J, Cubells JF. The structure of linkage disequilibrium at the DBH locus strongly influences the magnitude of association between diallelic markers and plasma dopamine beta-hydroxylase activity Am J Hum Genet 72: 1389-1400

Zapata C, Alvarez G, Carollo C (1997) Approximate variance of the standardized measure of gametic disequilibrium D'. Am. J. Hum. Genet. 61:771-774

Note

extracted from 2ld.c

Author(s)

Jing Hua Zhao

See Also

LDkl

Examples

1
2
3
4
5
6
7
## Not run: 
h <- c(0.442356,0.291532,0.245794,0.020319)
n <- 481*2
t <- LD22(h,n)
t

## End(Not run)


Search within the gap 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.