Martingale Posterior Distributions for Log-concave Density Functions
Preprint
Abstract
The family of log-concave density functions contains various kinds of common probability distributions. Due to the shape restriction, it is possible to find the nonparametric estimate of the density, for example, the nonparametric maximum likelihood estimate (NPMLE). However, the associated uncertainty quantification of the NPMLE is less well developed. The current techniques for uncertainty quantification are Bayesian, using a Dirichlet process prior combined with the use of Markov chain Monte Carlo (MCMC) to sample from the posterior. In this paper, we start with the NPMLE and use a version of the martingale posterior distribution to establish uncertainty about the NPMLE. The algorithm can be implemented in parallel and hence is fast. We prove the convergence of the algorithm by constructing suitable submartingales. We also illustrate results with different models and settings and some real data, and compare our method with that within the literature.Keywords
All Related Versions
This publication has 4 references indexed in Scilit:
- A Bayesian nonparametric approach to log-concave density estimationBernoulli, 2020
- A fast algorithm for univariate log‐concave density estimationAustralian & New Zealand Journal of Statistics, 2018
- Maximum likelihood estimation of a log-concave density and its distribution function: Basic properties and uniform consistencyBernoulli, 2009
- Uniform Convergence in Probability and Stochastic EquicontinuityEconometrica, 1991