Adaptive Combination of P-values (ADA) Method for Rare-Variant Association Testing

##########################################################################################
# This program was based on Cheung et al.'s code. Thank them for sharing their R code.
# Please cite:
# Lin W-Y, Lou XY, Gao G, Liu N (2014). Rare Variant Association Testing by Adaptive Combination of P-values. PLoS ONE, 9: e85728. [PMID: 24454922].
# and also the paper by Cheung et al. (2012) (listed below). Thank you.
# Any question, please contact: Wan-Yu Lin, linwy@ntu.edu.tw, Institute of Epidemiology and Preventive Medicine, National Taiwan University

#########################################################################
# Sigma-P Method for Rare Variant Analysis.
# Developer: Yee Him Cheung, Columbia University
# Email: patrick@ee.columbia.edu
# Copyright 2012
#
# Reference:
# Y.H. Cheung, G. Wang, S.M. Leal, S. Wang, "A Fast and Noise-Resilient
# Approach to Detect Rare-Variant Associations with Deep Sequencing Data
# for Complex Disorders", Genetic Epidemiology, Sep. 2012.
#########################################################################


The R code to implement the ADA method       Example file: ADAfile

In R, the code to implement this function is:


source("ADA.r")
ADATest('ADAfile.txt', 0.05, 1000, 1, 'additive', TRUE)        # wait for ~2 seconds


where 'ADAfile.txt' is the data file, in which the first column is the disease status (1: case; 0: control), and column 2 to the last column are the numbers of minor alleles.

The following input elements of this function are:
mafThr = 0.05,           (Exclude SNPs with combined MAF > 0.05)
num_perm = 1000,    (the number of permutations)
midp = TRUE,             (mid P-values according to the Fisher's exact test)
mode = 'additive',      (mode of inheritance = "additive")
twoSided = TRUE        (two-sided test)

Output: The output is the P-value of the ADA test.  [Lin W-Y, Lou XY, Gao G, Liu N (2014). Rare Variant Association Testing by Adaptive Combination of P-values. PLoS ONE, 9: e85728.]


About the simulation details in the ADA paper:
Simulation program to generate sequence-based data with dichotomous traits




Thanks for your interest.


Return to Wan-Yu Lin's homepage