allelic.exact.test: Fast Unbiased Exact Allelic Test

Description Usage Arguments Value See Also Examples

View source: R/allelic.exact.test.R

Description

This is the implementation in R+C of a new association test described in "A fast, unbiased and exact allelic exact test for case-control association studies" (Submitted). It appears that in most cases the classical chi-square test used for testing for allelic association on genotype data is biased. Our test is unbiased, exact but fast throught careful optimization.

Usage

1
 allelic.exact.test(d0,d1,d2,h0,h1,h2)

Arguments

Takes the 2x3 contingency table on which to compute the test

aa aA AA
[case (diseased)] d0 d1 d2
[control(healthy)] h0 h1 h2
d0

nb of first homozygous among cases

d1

nb of heterozygous among cases

d2

nb of second homozygous among cases

h0

nb of first homozygous among controls

h1

nb of heterozygous among controls

h2

nb of second homozygous among controls

Value

return the p-value of the test, or -1 if the sum of all cells in table is greater than TABLE\_OF\_LOG\_FACTORIALS\_SIZE, a C symbol defined in src/newallelic.c

See Also

chisq.test, fisher.test

Examples

1
 allelic.exact.test(160,80,60,160,160,30)

Example output

[1] 0.5074332

allelic documentation built on May 1, 2019, 10:31 p.m.