Geodesic slice sampling on the sphere

Fuente: arXiv
Gespeichert in:
Bibliographische Detailangaben
Hauptverfasser: Habeck, Michael, Hasenpflug, Mareike, Kodgirwar, Shantanu, Rudolf, Daniel
Format: Preprint
Veröffentlicht: 2023
Schlagworte:
Online-Zugang:
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
_version_ 1866918379022974976
author Habeck, Michael
Hasenpflug, Mareike
Kodgirwar, Shantanu
Rudolf, Daniel
author_facet Habeck, Michael
Hasenpflug, Mareike
Kodgirwar, Shantanu
Rudolf, Daniel
contents Probability measures on the sphere form an important class of statistical models and are used, for example, in modeling directional data or shapes. Due to their widespread use, but also as an algorithmic building block, efficient sampling of distributions on the sphere is highly desirable. We propose a shrinkage based and an idealized geodesic slice sampling Markov chain, designed to generate approximate samples from distributions on the sphere. In particular, the shrinkage-based version of the algorithm can be implemented such that it runs efficiently and has no tuning parameters. We verify reversibility and prove that under weak regularity conditions geodesic slice sampling is uniformly ergodic. Numerical experiments show that the proposed slice samplers achieve excellent mixing on challenging targets including distributions arising in rigid-registration problems and mixtures of von Mises-Fisher distributions. In these settings our approach outperforms standard samplers such as random-walk Metropolis-Hastings and Hamiltonian Monte Carlo.
format Preprint
id arxiv_https___arxiv_org_abs_2301_08056
institution arXiv
publishDate 2023
record_format arxiv
spellingShingle Geodesic slice sampling on the sphere
Habeck, Michael
Hasenpflug, Mareike
Kodgirwar, Shantanu
Rudolf, Daniel
Methodology
Probability
Statistics Theory
60J22, 58C35, 65C40
Probability measures on the sphere form an important class of statistical models and are used, for example, in modeling directional data or shapes. Due to their widespread use, but also as an algorithmic building block, efficient sampling of distributions on the sphere is highly desirable. We propose a shrinkage based and an idealized geodesic slice sampling Markov chain, designed to generate approximate samples from distributions on the sphere. In particular, the shrinkage-based version of the algorithm can be implemented such that it runs efficiently and has no tuning parameters. We verify reversibility and prove that under weak regularity conditions geodesic slice sampling is uniformly ergodic. Numerical experiments show that the proposed slice samplers achieve excellent mixing on challenging targets including distributions arising in rigid-registration problems and mixtures of von Mises-Fisher distributions. In these settings our approach outperforms standard samplers such as random-walk Metropolis-Hastings and Hamiltonian Monte Carlo.
title Geodesic slice sampling on the sphere
topic Methodology
Probability
Statistics Theory
60J22, 58C35, 65C40
url https://arxiv.org/abs/2301.08056