A Geometric Perspective on the Power of Principal Component Association Tests in Multiple Phenotype Studies. Issue 527 (3rd July 2019)

• Record Type: Journal Article
• Authors: Liu, Zhonghua; Lin, Xihong
• Is Part Of: Journal of the American Statistical Association. Volume 114:Issue 527(2019)
• Page Start: 975
• Page End: 990
• Publication Date: 2019-07-03
• Subjects: Dimension reduction -- Genome-wide Association Studies (GWAS) -- Omnibus test -- Power analysis -- Principal angle -- Summary statistics; Statistics -- Periodicals; Statistics -- Periodicals; Statistiques -- Périodiques; États-Unis -- Statistiques -- Périodiques; 519.5
• DOI: 10.1080/01621459.2018.1513363
• Languages: English
• ISSNs: 0162-1459

Abstract

ABSTRACT: Joint analysis of multiple phenotypes can increase statistical power in genetic association studies. Principal component analysis, as a popular dimension reduction method, especially when the number of phenotypes is high dimensional, has been proposed to analyze multiple correlated phenotypes. It has been empirically observed that the first PC, which summarizes the largest amount of variance, can be less powerful than higher-order PCs and other commonly used methods in detecting genetic association signals. In this article, we investigate the properties of PCA-based multiple phenotype analysis from a geometric perspective by introducing a novel concept called principal angle. A particular PC is powerful if its principal angle is 0 ° and is powerless if its principal angle is 90 ° . Without prior knowledge about the true principal angle, each PC can be powerless. We propose linear, nonlinear, and data-adaptive omnibus tests by combining PCs. We demonstrate that the Wald test is a special quadratic PC-based test. We show that the omnibus PC test is robust and powerful in a wide range of scenarios. We study the properties of the proposed methods using power analysis and eigen-analysis. The subtle differences and close connections between these combined PC methods are illustrated graphically in terms of their rejection boundaries. Our proposed tests have convex acceptance regions and hence are admissible. The p -values for the proposed tests can be efficiently calculated analytically and the proposed tests have been implemented in a publicly available R package MPAT . We conduct simulation studies in both low- and high-dimensional settings with various signal vectors and correlation structures. We apply the proposed tests to the joint analysis of metabolic syndrome-related phenotypes with datasets collected from four international consortia to demonstrate the effectiveness of the proposed combined PC testing procedures. Supplementary materials for this article, including a standardized description of the materials available for reproducing the work, are available as an online supplement.