Abstract:
|
Joint analysis of multiple phenotypes in genome-wide association studies (GWASs) is important for detecting pleiotropic effects. For secondary phenotypes in case-control studies, the over-sampling of the cases relative to the primary (disease) phenotype can induce bias in genetic association analysis with secondary phenotypes. To account for this ascertainment bias, we propose using weighted estimating equations for fitting the scaled mean model. In contrast to earlier work by Schifano et al. (2013) where the weight is the inverse of fixed, known sampling fractions/probabilities of the cases and controls, the proposed weight is a population-based measure, defined as the inverse of the conditional probability of case-control status. The proposed method is robust to departures from normality of multiple continuous phenotypes and the misspecification of within-individual correlation across the multiple phenotypes. We show through simulation that the proposed method out-performs other methods and other weighting strategies. We illustrate our approach using a lung cancer case-control study.
|