A Bayesian nonparametric approach to log-concave density estimation
Abstract
The estimation of a log-concave density on $\mathbb{R}$ is a canonical problem in the area of shape-constrained nonparametric inference. We present a Bayesian nonparametric approach to this problem based on an exponentiated Dirichlet process mixture prior and show that the posterior distribution converges to the log-concave truth at the (near-) minimax rate in Hellinger distance. Our proof proceeds by establishing a general contraction result based on the log-concave maximum likelihood estimator that prevents the need for further metric entropy calculations. We also present two computationally more feasible approximations and a more practical empirical Bayes approach, which are illustrated numerically via simulations.
- Publication:
-
arXiv e-prints
- Pub Date:
- March 2017
- DOI:
- 10.48550/arXiv.1703.09531
- arXiv:
- arXiv:1703.09531
- Bibcode:
- 2017arXiv170309531M
- Keywords:
-
- Mathematics - Statistics Theory;
- 62G07;
- 62G20
- E-Print:
- 39 pages, 17 figures. Simulation studies were significantly expanded and one more theorem has been added