Saved in:
Bibliographic Details
Main Authors: Yao, Tsung-Hung, Wu, Zhenke, Bharath, Karthik, Baladandayuthapani, Veerabhadran
Format: Preprint
Published: 2024
Subjects:
Online Access:https://arxiv.org/abs/2401.11515
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866916706742435840
author Yao, Tsung-Hung
Wu, Zhenke
Bharath, Karthik
Baladandayuthapani, Veerabhadran
author_facet Yao, Tsung-Hung
Wu, Zhenke
Bharath, Karthik
Baladandayuthapani, Veerabhadran
contents Ultrametric matrices are a class of covariance matrices that arise in latent tree models. As a parameter space in a statistical model, the set of ultrametric matrices is neither convex nor a smooth manifold. Focus in the literature has hitherto been restricted to estimation through projections and relaxation-based techniques, and inferential methods are lacking. Motivated by this, we establish a bijection between the set of positive definite ultrametric matrices and the set of rooted, leaf-labeled trees equipped with the stratified geometry of the well-known phylogenetic treespace. Using the pullback geometry under the bijection and by adapting sampling algorithms in Bayesian phylogenetics, we develop algorithms to sample from the posterior distribution on the set of ultrametric matrices in a Bayesian latent tree model where the tree may be binary or multifurcating. We demonstrate the utility of the algorithms in simulation studies, and illustrate them on a pre-clinical cancer application to quantify uncertainty about treatment trees that identify treatments with high mechanism similarity that target correlated pathways.
format Preprint
id arxiv_https___arxiv_org_abs_2401_11515
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle Geometry-driven Bayesian Inference for Ultrametric Covariance Matrices
Yao, Tsung-Hung
Wu, Zhenke
Bharath, Karthik
Baladandayuthapani, Veerabhadran
Methodology
Ultrametric matrices are a class of covariance matrices that arise in latent tree models. As a parameter space in a statistical model, the set of ultrametric matrices is neither convex nor a smooth manifold. Focus in the literature has hitherto been restricted to estimation through projections and relaxation-based techniques, and inferential methods are lacking. Motivated by this, we establish a bijection between the set of positive definite ultrametric matrices and the set of rooted, leaf-labeled trees equipped with the stratified geometry of the well-known phylogenetic treespace. Using the pullback geometry under the bijection and by adapting sampling algorithms in Bayesian phylogenetics, we develop algorithms to sample from the posterior distribution on the set of ultrametric matrices in a Bayesian latent tree model where the tree may be binary or multifurcating. We demonstrate the utility of the algorithms in simulation studies, and illustrate them on a pre-clinical cancer application to quantify uncertainty about treatment trees that identify treatments with high mechanism similarity that target correlated pathways.
title Geometry-driven Bayesian Inference for Ultrametric Covariance Matrices
topic Methodology
url https://arxiv.org/abs/2401.11515