Nonasymptotic bounds for the quadratic risk of the Grenander estimator

There is an enormous literature on the so-called Grenander estimator, which is merely the nonparametric maximum likelihood estimator of a nonincreasing probability density on [0, 1] (see, for instance, Grenander (1981)), but unfortunately, there is no nonasymptotic (i.e. for arbitrary finite sample size n) explicit upper bound for the quadratic risk of the Grenander estimator readily applicable in practice by statisticians. In this paper, we establish, for the first time, a simple explicit upper bound 2n−1/2 for the latter quadratic risk. It turns out to be a straightforward consequence of an inequality valid with probability one and bounding from above the integrated squared error of the Grenander estimator by the Kolmogorov–Smirnov statistic.

Harmonic moments of inhomogeneous branching processes

We study the harmonic moments of Galton-Watson processes that are possibly inhomogeneous and have positive values. Good estimates of these are needed to compute unbiased estimators for noncanonical branching Markov processes, which occur, for instance, in the modelling of the polymerase chain reaction. By convexity, the ratio of the harmonic mean to the mean is at most 1. We prove that, for every square-integrable branching mechanism, this ratio lies between 1-A/k and 1-A/k for every initial population of size k>A. The positive constants A and Aͤ are such that A≥Aͤ, are explicit, and depend only on the generation-by-generation branching mechanisms. In particular, we do not use the distribution of the limit of the classical martingale associated with the Galton-Watson process. Thus, emphasis is put on nonasymptotic bounds and on the dependence of the harmonic mean upon the size of the initial population. In the Bernoulli case, which is relevant for the modelling of the polymerase chain reaction, we prove essentially optimal bounds that are valid for every initial population size k≥1. Finally, in the general case and for sufficiently large initial populations, similar techniques yield sharp estimates of the harmonic moments of higher degree.

Harmonic moments of inhomogeneous branching processes

We study the harmonic moments of Galton-Watson processes that are possibly inhomogeneous and have positive values. Good estimates of these are needed to compute unbiased estimators for noncanonical branching Markov processes, which occur, for instance, in the modelling of the polymerase chain reaction. By convexity, the ratio of the harmonic mean to the mean is at most 1. We prove that, for every square-integrable branching mechanism, this ratio lies between 1-A/k and 1-A/k for every initial population of size k>A. The positive constants A and Aͤ are such that A≥Aͤ, are explicit, and depend only on the generation-by-generation branching mechanisms. In particular, we do not use the distribution of the limit of the classical martingale associated with the Galton-Watson process. Thus, emphasis is put on nonasymptotic bounds and on the dependence of the harmonic mean upon the size of the initial population. In the Bernoulli case, which is relevant for the modelling of the polymerase chain reaction, we prove essentially optimal bounds that are valid for every initial population size k≥1. Finally, in the general case and for sufficiently large initial populations, similar techniques yield sharp estimates of the harmonic moments of higher degree.