Estimation of unimodal densities without smoothness assumptions

The Grenander estimator of a decreasing density, which is defined as the derivative of the concave envelope of the empirical c.d.f., is known to be a very good estimator of an unknown decreasing density on the half-line R + when this density is not assumed to be smooth. It is indeed the maximum likelihood estimator and one can get precise upper bounds for its risk when the loss is measured by the L 1 -distance between densities. Moreover, if one restricts oneself to the compact subsets of decreasing densities bounded by H with support on [0, L] the risk of this estimator is within a fixed factor of the minimax risk. The same is true if one deals with the maximum likelihood estimator for unimodal densities with known mode. When the mode is unknown, the maximum likelihood estimator does not exist any more. We shall provide a general purpose estimator (together with a computational algorithm) for estimating nonsmooth unimodal densities. Its risk is the same as the risk of the Grenander estimator based on the knowledge of the true mode plus some lower order term. It can also cope with small departures from unimodality.

[1]  U. Grenander On the theory of mortality measurement , 1956 .

[2]  Tim Robertson,et al.  On Estimating a Density which is Measurable with Respect to a $\sigma$-Lattice , 1967 .

[3]  Prakasa Rao Estimation of a unimodal density , 1969 .

[4]  A note on estimating a unimodal density , 1969 .

[5]  Edward J. Wegman,et al.  Maximum Likelihood Estimation of a Unimodal Density Function , 1970 .

[6]  E. Wegman Maximum likelihood estimation of a unimodal density. II , 1970 .

[7]  R. Reiss On the Measurability and Consistency of Maximum Likelihood Estimates for Unimodal Densities , 1973 .

[8]  J. Wolfowitz,et al.  Asymptotically minimax estimation of concave and convex distribution functions , 1976 .

[9]  R. Reiss On minimum distance estimators for unimodal densities , 1976 .

[10]  R. Hasminskii Lower Bound for the Risks of Nonparametric Estimates of the Mode , 1979 .

[11]  Prakasa Rao Nonparametric functional estimation , 1983 .

[12]  P. Groeneboom Estimating a monotone density , 1984 .

[13]  S. Lo Estimation of a Unimodal Distribution Function , 1986 .

[14]  L. Birge On the Risk of Histograms for Estimating Decreasing Densities , 1987 .

[15]  Atsuyuki Kogure,et al.  Asymptotically Optimal Cells for a Historgram , 1987 .

[16]  L. Devroye A Course in Density Estimation , 1987 .

[17]  L. Birge Estimating a Density under Order Restrictions: Nonasymptotic Minimax Risk , 1987 .

[18]  F. T. Wright,et al.  Order restricted statistical inference , 1988 .

[19]  L. Birge,et al.  The Grenader Estimator: A Nonasymptotic Approach , 1989 .

[20]  Enno Mammen,et al.  Estimating a Smooth Monotone Regression Function , 1991 .

[21]  G. Kerkyacharian,et al.  Density estimation in Besov spaces , 1992 .

[22]  I. Johnstone,et al.  Ideal spatial adaptation by wavelet shrinkage , 1994 .