Taylor & Francis Group
Browse
uasa_a_1699422_sm5328.pdf (328.3 kB)

Optimal Estimation of Wasserstein Distance on a Tree With an Application to Microbiome Studies

Download (328.3 kB)
Version 2 2020-01-23, 14:41
Version 1 2019-12-07, 06:29
journal contribution
posted on 2020-01-23, 14:41 authored by Shulei Wang, T. Tony Cai, Hongzhe Li

The weighted UniFrac distance, a plug-in estimator of the Wasserstein distance of read counts on a tree, has been widely used to measure the microbial community difference in microbiome studies. Our investigation however shows that such a plug-in estimator, although intuitive and commonly used in practice, suffers from potential bias. Motivated by this finding, we study the problem of optimal estimation of the Wasserstein distance between two distributions on a tree from the sampled data in the high-dimensional setting. The minimax rate of convergence is established. To overcome the bias problem, we introduce a new estimator, referred to as the moment-screening estimator on a tree (MET), by using implicit best polynomial approximation that incorporates the tree structure. The new estimator is computationally efficient and is shown to be minimax rate-optimal. Numerical studies using both simulated and real biological datasets demonstrate the practical merits of MET, including reduced biases and statistically more significant differences in microbiome between the inactive Crohn’s disease patients and the normal controls. Supplementary materials for this article are available online.

Funding

This research was supported by NIH grants R01GM123056 and R01GM129781.

History

Usage metrics

    Journal of the American Statistical Association

    Licence

    Exports

    RefWorks
    BibTeX
    Ref. manager
    Endnote
    DataCite
    NLM
    DC