What's the best statistic for a simple test of genetic association in a case-control study?
Genome-wide genetic association studies typically start with univariate statistical tests of each marker. In principle, this single-SNP scanning is statistically straightforward--the testing is done with standard methods (e.g. chi(2) tests, regression) that have been well studied for decades. However, a number of different tests and testing procedures can be used. In a case-control study, one can use a 1 df allele-based test, a 1 or 2 df genotype-based test, or a compound procedure that combines two or more of these statistics. Additionally, most of the tests can be performed with or without covariates included in the model. While there are a number of statistical papers that make power comparisons among subsets of these methods, none has comprehensively tackled the question of which of the methods in common use is best suited to univariate scanning in a genome-wide association study. In this paper, we consider a wide variety of realistic test procedures, and first compare the power of the different procedures to detect a single locus under different genetic models. We then address the question of whether or when it is a good idea to include covariates in the analysis. We conclude that the most commonly used approach to handle covariates--modeling covariate main effects but not interactions--is almost never a good idea. Finally, we consider the performance of the statistics in a genome scan context.
Kuo
Chia-Ling
CL
Department of Biostatistics, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, Pennsylvania, USA.
Feingold
Eleanor
E
D000483:Alleles; D016022:Case-Control Studies; D003627:Data Interpretation, Statistical; D056726:Genetic Association Studies; D020022:Genetic Predisposition to Disease; D055106:Genome-Wide Association Study; D005838:Genotype; D006579:Heterozygote; D006801:Humans; D008957:Models, Genetic; D015233:Models, Statistical; D012044:Regression Analysis