On Missing Data and Genotyping Errors in Association Studies

Transcription

1 On Missing Data and Genotyping Errors in Association Studies Department of Biostatistics Johns Hopkins Bloomberg School of Public Health May 16, 2008

2 Specific Aims of our R01 1 Develop and evaluate new statistical methods to prioritize genes through proper ranking in genome-wide association (GWA) studies that address GxE interactions. 2 Develop and evaluate new statistical methods to localize causal genes as part of linkage and fine mapping studies while considering GxE interactions. 3 Develop and evaluate new statistical methods to identify higher order interactions between environmental variables and SNPs in candidate genes studies.

3 Specific Aims (cont.) 4 Adapt existing and develop new statistical methods to address imprecise and missing environmental and genetic measurements. 5 Develop and disseminate efficient algorithms for GxE analyses, and apply these methods in several ongoing genetic studies of complex diseases.

4 Missing Data There are mainly three types of missing / unobserved data in genetic association studies: 1 Missing observations in some environmental variables. 2 Missing data at SNPs selected for genotyping. 3 Genotypes of SNPs not selected.

5 Approaches for Missing Data The most common approach for dealing with missing data is to omit the observations that have missing records in the model s covariates. This approach can have several shortcomings, including: Loss of power. Bias in the parameter estimates. A good reference on this topic is Greenland and Finkle (1995). Some other used approaches are: To impute a value from the marginal distribution of the covariate. To create an extra level indicating missingness, if the covariate is a factor. These choices tend to be not so great either.

6 Approaches for Missing Data Multiple imputation can be used to draw valid statistical inference from data with missing values when the data are missing at random (Little and Rubin 1987, Schafer 1997). In essence, multiple imputation acknowledges the uncertainty due to missing data, instead of simply ignoring it: several complete data sets are generated, and the uncertainty in the model parameter estimates incorporates the standard errors of the parameter estimates as well as the variability between the parameter estimates from the replicate data sets. While the hypothesis of missing at random cannot formally be tested, it is a lot less stringent than the requirement of missing completely at random, which is the underlying assumption made when observations are omitted.

7 Missing Environmental Data Number of Pairs Odds Ratio Confidence Interval XPD Lys751Gln original data set ( ) multiple imputations ( ) XPD Gln751Gln original data set ( ) multiple imputations ( ) Positive Family History original data set ( ) multiple imputations ( )

8 Missing Environmental Data Family History not complete Family History complete AA AC CC na AA AC CC na raw numbers case control percentages case control Reference: Brewster AM et al (2006). Polymorphisms of the DNA Repair Genes XPD (Lys751Gln) and XRCC1 (Arg399Gln and Arg194Trp): Relationship to Breast Cancer Risk and Familial Predisposition to Breast Cancer. Breast Cancer Res Treat, 95(1):

9 Missing Environmental Data 4 XPD Lys751Gln odds ratio XPD Gln751Gln odds ratio Family History odds ratio combined original The missing data were imputed using decision trees. Reference: Dai J et al (2006). Imputation Methods to Improve Inference in SNP Association Studies. Genetic Epidemiology, 30(8):

10 Why Become a Biostatistician? Because people appreciate your help analyzing their data, and that means that people surely will like you.

11 Tree-based Imputation Classification trees are great for categorical data!

12 Dummy Levels for Missingness hazards ratio SNP 1 SNP 2 SNP 3 SNP 4 SNP 5 SNP 6 SNP 7 SNP 8 SNP 9 SNP 10 SNP 11 SNP 12 SNP 13 Unpublished data.

13 Incorporating Genotype Uncertainty The confidence in genotype calls can differ substantially between SNPs! 4 Concordant 2 AA AB BB Discordant called BB (AB) Sense Antisense

14 Missingness at Random? From the white paper, updates/brlmm algorithm.affx

15 Incorporating Genotype Uncertainty 100 [ truncated ] log(ratio) Ratio correlation(genotype,crlmm) / correlation(genotype,brlmm)

16 Incorporating Genotype Uncertainty Easy SNP Difficult SNP Median Rank [ 10,000 simulations ] True Genotype CRLMM (continuous) CRLMM (called) 75 BRLMM (all) BRLMM (selected) RR: θ RR: θ

17 Incorporating Genotype Uncertainty 100 [ truncated ] log(ratio) Ratio correlation(genotype,crlmm) / correlation(genotype,crlmm[called])

18 Incorporating Genotype Uncertainty - CNVs 5.0 A D B C E

19 Incorporating Genotype Uncertainty - CNVs Deletion Normal LOH Amplification A D B C E 2 1 Van ICE A D B E Van ICE Mb

20 A HapMap Sample Deletion Normal LOH Amplification 1 Van ICE Van ICE Mb

21 Many HapMap Samples

22 De Novo Deletion