Description Details Author(s) References Examples
Tests for marginal and pairwise SNP associations with binary phenotypes in case-control, case-population, and case-control-population studies.
Package: | CCpop |
Type: | Package |
Version: | 1.0 |
Date: | 2014-03-18 |
License: | GPL-2 |
Inputs to test functions are vectors (for marginal tests) and matrices (for pairwise tests) of genotypic counts and are given separately for cases, controls, and population cohorts (where relevant). The constrained maximum likelihood estimation based tests (Kaufman and Rosset, 2014) also require a value for the known phenotypic prevalence in the population (a crude estimate is usually sufficient). When applicable, tests assuming Hardy-Weinberg equilibrium (HWE) and linkage equilibrium (LE) will be considerably more powerful.
See marginal.assoc.tests for single SNP tests, and pairwise.assoc.tests for pairwise joint, pure-interaction, and conditional tests.
Shachar Kaufman <shachark@post.tau.ac.il>
Kaufman, S., & Rosset, S. (2014). Exploiting Population Samples To Enhance Genome-Wide Association Studies of Disease. Genetics, genetics-114.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 | ## An example marginal/pairwise association
# Controls
t0 = matrix(c(375, 240, 46,
640, 405, 62,
300, 169, 19), nrow = 3, byrow = TRUE)
# Cases
t1 = matrix(c(317, 162, 15,
459, 209, 22,
120, 76, 13), nrow = 3, byrow = TRUE)
# Independent population sample, marginalized for SNP1 and SNP2
tp1 = c(2410, 4253, 1945)
tp2 = c(4972, 3140, 496)
## The prevalence of the studied disease in the population
prevalence = 0.001
marginal.assoc.test.pop.hwe.kpy(t0 = rowSums(t0), t1 = rowSums(t1), tp = tp1, prevalence)
marginal.assoc.test.pop.hwe.kpy(t0 = colSums(t0), t1 = colSums(t1), tp = tp2, prevalence)
pairwise.assoc.test.pop.hwe.le.kpy(t0, t1, tp1, tp2, prevalence)
conditional.assoc.test.pure.pop.hwe.le.kpy(t0, t1, tp1, tp2, prevalence)
|
Loading required package: nloptr
$pen
[1] 0.0012653534 0.0009941329 0.0006774466
$statistic
[1] 54.86703
$p.value
[1] 1.21836e-12
$pen
[1] 0.0011115798 0.0008813761 0.0006267568
$statistic
[1] 26.73425
$p.value
[1] 1.565783e-06
$pen
[,1] [,2] [,3]
[1,] 0.0014033346 0.0011397123 0.0006703782
[2,] 0.0011427526 0.0008270668 0.0005537271
[3,] 0.0006722005 0.0006768101 0.0007364398
$statistic
[1] 89.73377
$p.value
[1] 5.266646e-16
$statistic
[1] 34.86674
$p.value
[1] 4.573408e-06
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.