Stacked Grenander estimator of a discrete distribution

In this paper we combine isotonic regression with least-squares cross-validation based model-mix method to estimate a discrete distribution with an infinite support. The new estimator is proved to be strongly consistent with $\sqrt{n}$-rate of convergence, it behaves well for moderate sized data sets and provides a trade-off between goodness-of-fit and shape constraints.