40
13

Mode-constrained estimation of a log-concave density

Abstract

We study nonparametric maximum likelihood estimation of a log-concave density function ff with a known mode mm. We develop asymptotic theory for the mode-constrained estimator, including, consistency, global rates of convergence, and local rates of convergence both in neighborhoods of the specified mode and away from the known mode. In local neighborhoods of the specified mode the constrained and unconstrained estimators differ, but away from the specified mode the constrained and unconstrained estimators are asymptotically equivalent. Software to compute the mode-constrained estimator is available in the R package \verb+logcondens.mode+. In a companion paper we use the mode-constrained MLE to develop a likelihood ratio test of the null hypothesis that the mode of ff equals a specified valued mm versus the alternative hypothesis that the mode of ff differs from the specified valued mm. We show that under the null hypothesis (and strict curvature of logf\log f at the mode) the natural likelihood ratio statistic is asymptotically pivotal: that is, it converges in distribution to a limiting distribution which is free of nuisance parameters, thus playing the role of the χ12\chi_1^2 distribution in classical parametric statistical problems. The test can be inverted to form confidence intervals.

View on arXiv
Comments on this paper