Increasing Power for Tests of Genetic Association in the Presence of Phenotype and/or Genotype Error by Use of Double-Sampling

Derek Gordon, Rockefeller University
Yaning Yang, Rockefeller University
Chad Haynes, Rockefeller University
Stephen J. Finch, Stony Brook University
Nancy R. Mendell, Stony Brook University
Abraham M. Brown, Burke Medical Research Institute
Vahram Haroutunian, Mount Sinai School of Medicine

Abstract

Phenotype and/or genotype misclassification can: significantly increase type II error probabilities for genetic case/control association, causing decrease in statistical power; and produce inaccurate estimates of population frequency parameters. We present a method, the likelihood ratio test allowing for errors (LRTae) that incorporates double-sample information for phenotypes and/or genotypes on a sub-sample of cases/controls. Population frequency parameters and misclassification probabilities are determined using a double-sample procedure as implemented in the Expectation-Maximization (EM) method. We perform null simulations assuming a SNP marker or a 4-allele (multi-allele) marker locus. To compare our method with the standard method that makes no adjustment for errors (LRTstd), we perform power simulations using a 2^k factorial design with high and low settings of: case/control samples, phenotype/genotype costs, double-sampled phenotypes/genotypes costs, phenotype/genotype error, and proportions of double-sampled individuals. All power simulations are performed fixing equal costs for the LRTstd and LRTae methods. We also consider case/control ApoE genotype data for an actual Alzheimer's study.

The LRTae method maintains correct type I error proportions for all null simulations and all significance level thresholds (10%, 5%, 1%). LRTae average estimates of population frequencies and misclassification probabilities are equal to the true values, with variances of 10e-7 to 10e-8. For power simulations, the median power difference LRTae-LRTstd at the 5% significance level is 0.06 for multi-allele data and 0.01 for SNP data. For the ApoE data example, the LRTae and LRTstd p-values are 5.8 x 10e-5 and 1.6 x 10e-3, respectively. The increase in significance is due to adjustment in the LRTae for misclassification of the most commonly reported risk allele. We have developed freely available software that performs our LRTae statistic.

Submitted: August 6, 2004 · Accepted: September 28, 2004 · Published: October 6, 2004

Recommended Citation

Gordon, Derek; Yang, Yaning; Haynes, Chad; Finch, Stephen J.; Mendell, Nancy R.; Brown, Abraham M.; and Haroutunian, Vahram (2004) "Increasing Power for Tests of Genetic Association in the Presence of Phenotype and/or Genotype Error by Use of Double-Sampling ," Statistical Applications in Genetics and Molecular Biology: Vol. 3 : Iss. 1, Article 26.
Available at: http://www.bepress.com/sagmb/vol3/iss1/art26

Readers' Reactions

Fei Ji, Yaning Yang, Chad Haynes, Stephen J. Finch, and Derek Gordon, Computing Asymptotic Power and Sample Size for Case-Control Genetic Association Studies in the Presence of Phenotype and/or Genotype Misclassification Errors (January 2006)

 
 
 
 

ISSN: 1544-6115 ©1999-2008 The Berkeley Electronic Press™ All rights reserved.

To submit, subscribe, recommend this journal to your library, or sign up for email alerts, please visit: http://www.bepress.com/sagmb