Journal
JOURNAL OF THE AMERICAN STATISTICAL ASSOCIATION
Volume 114, Issue 526, Pages 759-772Publisher
AMER STATISTICAL ASSOC
DOI: 10.1080/01621459.2018.1442340
Keywords
Adaptive thresholding; Basis covariance; Centered log-ratio covariance; High dimensionality; Microbiome; Regularization
Categories
Funding
- NIH [CA127334, GM097505]
- NSFC [11671018, 71532001]
Ask authors/readers for more resources
High-dimensional compositional data arise naturally in many applications such as metagenomic data analysis. The observed data lie in a high-dimensional simplex, and conventional statistical methods often fail to produce sensible results due to the unit-sum constraint. In this article, we address the problem of covariance estimation for high-dimensional compositional data and introduce a composition-adjusted thresholding (COAT) method under the assumption that the basis covariance matrix is sparse. Our method is based on a decomposition relating the compositional covariance to the basis covariance, which is approximately identifiable as the dimensionality tends to infinity. The resulting procedure can be viewed as thresholding the sample centered log-ratio covariance matrix and hence is scalable for large covariance matrices. We rigorously characterize the identifiability of the covariance parameters, derive rates of convergence under the spectral norm, and provide theoretical guarantees on support recovery. Simulation studies demonstrate that the COAT estimator outperforms some existing optimization-based estimators. We apply the proposed method to the analysis of a microbiome dataset to understand the dependence structure among bacterial taxa in the human gut.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available