Dataset Information

Principal interactions analysis for repeated measures data: application to gene-gene and gene-environment interactions.

ABSTRACT: Many existing cohorts with longitudinal data on environmental exposures, occupational history, lifestyle/ behavioral characteristics, and health outcomes have collected genetic data in recent years. In this paper, we consider the problem of modeling gene-gene and gene-environment interactions with repeated measures data on a quantitative trait. We review possibilities of using classical models proposed by Tukey (1949) and Mandel (1961) using the cell means of a two-way classification array for such data. Although these models are effective for detecting interactions in the presence of main effects, they fail miserably if the interaction structure is misspecified. We explore a more robust class of interaction models that are based on a singular value decomposition of the cell-means residual matrix after fitting the additive main effect terms. This class of additive main effects and multiplicative interaction models (Gollob, 1968) provide useful summaries for subject-specific and time-varying effects as represented in terms of their contribution to the leading eigenvalues of the interaction matrix. It also makes the interaction structure more amenable to geometric representation. We call this analysis 'principal interactions analysis'. While the paper primarily focuses on a cell-mean-based analysis of repeated measures outcome, we also introduce resampling-based methods that appropriately recognize the unbalanced and longitudinal nature of the data instead of reducing the response to cell means. We illustrate the proposed methods by using data from the Normative Aging Study, a longitudinal cohort study of Boston area veterans since 1963. We carry out simulation studies under an array of classical interaction models and common epistasis models to illustrate the properties of the principal interactions analysis procedure in comparison with the classical alternatives.

SUBMITTER: Mukherjee B

PROVIDER: S-EPMC4046647 | biostudies-literature | 2012 Sep

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Principal interactions analysis for repeated measures data: application to gene-gene and gene-environment interactions.

Mukherjee Bhramar B Ko Yi-An YA Vanderweele Tyler T Roy Anindya A Park Sung Kyun SK Chen Jinbo J

Statistics in medicine 20120313 22

Many existing cohorts with longitudinal data on environmental exposures, occupational history, lifestyle/ behavioral characteristics, and health outcomes have collected genetic data in recent years. In this paper, we consider the problem of modeling gene-gene and gene-environment interactions with repeated measures data on a quantitative trait. We review possibilities of using classical models proposed by Tukey (1949) and Mandel (1961) using the cell means of a two-way classification array for s ...[more]

PMID: 22415818

Similar Datasets

Project description:Voxel-based principal-component analysis allows for an identification of patterns of glucose metabolism and amyloid deposition related to the conversion from mild cognitive impairment (MCI) to Alzheimer disease (AD). The present study aimed to validate these AD conversion-related patterns (ADCRPs) against neuropathologic findings. Methods: We included patients from the Alzheimer's Disease Neuroimaging Initiative who underwent autopsy and for whom 18F-FDG PET (30 AD, 6 MCI, 2 cognitively normal) and amyloid-β (Aβ) PET (17 AD, 3 MCI, 2 cognitively normal) were available. Pattern expression scores (PESs) of the 18F-FDG- and Aβ-ADCRP were compared with Braak tangle stage and Thal amyloid phase, respectively. Mean 18F-FDG uptake and mean 18F-AV-45 SUV ratio (SUVr) in regions of hypometabolism and elevated amyloid load typical of AD, respectively, were used as volume-of-interest-based PET measures. The diagnostic performance for identifying none-to-low vs. intermediate-to-high AD neuropathologic change (ADNC) was assessed for all biomarkers. Results: We observed significant associations between PES of 18F-FDG-ADCRP and Braak stage (ρ > 0.48, P < 0.005) and between PES of Aβ-ADCRP and Thal phase (ρ > 0.66, P < 0.001). PES of 18F-FDG-ADCRP, PES of Aβ-ADCRP, and their combination identified intermediate-to-high ADNC with an area under the receiver-operating-characteristic curve (AUC) of 0.80, 0.95, and 0.98 (n = 22), respectively. Mean 18F-FDG uptake and mean 18F-AV-45 SUVr in AD-typical regions were also significantly associated with Braak stage (|ρ| > 0.45, P < 0.01) and Thal phase (ρ > 0.55, P < 0.01), respectively. Volume-of-interest-based PET measures discriminated between ADNC stages with an AUC of 0.79, 0.88, and 0.90 for mean 18F-FDG uptake, mean 18F-AV-45 SUVr, and their combination (n = 22), respectively. Contemplating all subjects with available 18F-FDG PET and neuropathology information (n = 38), PES of 18F-FDG-ADCRP was a significant predictor of intermediate-to-high ADNC (AUC = 0.72), whereas mean 18F-FDG uptake was not (AUC = 0.66), although the difference between methods was not significant. Conclusion: PES of 18F-FDG-ADCRP, a measure of neurodegeneration, shows close correspondence with the extent of tau pathology, as assessed by Braak tangle stage. PES of Aβ-ADCRP is a valid biomarker of underlying amyloid pathology, as demonstrated by its strong correlation with Thal phase. The combination of ADCRPs performed better than 18F-FDG-ADCRP alone, although there was only negligible improvement compared with Aβ-ADCRP.

Project description:Longitudinal intervention studies with repeated measurements over time are an important type of experimental design in biomedical research. Due to the advent of “omics”-sciences (genomics, transcriptomics, proteomics, metabolomics), longitudinal studies generate increasingly multivariate outcome data. Analysis of such data must take both the longitudinal intervention structure and multivariate nature of the data into account. The ASCA+-framework combines general linear models with principal component analysis and can be used to separate and visualize the multivariate effect of different experimental factors. However, this methodology has not yet been developed for the more complex designs often found in longitudinal intervention studies, which may be unbalanced, involve randomized interventions, and have substantial missing data. Here we describe a new methodology, repeated measures ASCA+ (RM-ASCA+), and show how it can be used to model metabolic changes over time, and compare metabolic changes between groups, in both randomized and non-randomized intervention studies. Tools for both visualization and model validation are discussed. This approach can facilitate easier interpretation of data from longitudinal clinical trials with multivariate outcomes. Author summary Clinical trials are increasingly generating large amounts of complex biological data. Examples can include measuring metabolism or gene expression in tissue or blood sampled repeatedly over the course of a treatment. In such cases, one might wish to compare changes in not one, but hundreds, or thousands of variables simultaneously. In order to effectively analyze such data, both the study design and the multivariate nature of the data should be considered during data analysis. ANOVA simultaneous component analysis+ (ASCA+) is a statistical method which combines general linear models with principal component analysis, and provides a way to separate and visualize the effects of different factors on complex biological data. In this work, we describe how repeated measures linear mixed models, a class of models commonly used when analyzing changes over time and treatment effects in longitudinal studies, can be used together with ASCA+ for analyzing clinical trials in a novel method called repeated measures-ASCA+ (RM-ASCA+).

Dataset Information

Principal interactions analysis for repeated measures data: application to gene-gene and gene-environment interactions.

Publications

Principal interactions analysis for repeated measures data: application to gene-gene and gene-environment interactions.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets