  Westling and Carone (2020) proposed a framework for studying the large sample
distributional properties of generalized Grenander-type estimators, a versatile
class of nonparametric estimators of monotone functions. The limiting
distribution of those estimators is representable as the left derivative of the
greatest convex minorant of a Gaussian process whose covariance kernel can be
complicated and whose monomial mean can be of unknown order (when the degree of
flatness of the function of interest is unknown). The standard nonparametric
bootstrap is unable to consistently approximate the large sample distribution
of the generalized Grenander-type estimators even if the monomial order of the
mean is known, making statistical inference a challenging endeavour in
applications. To address this inferential problem, we present a
bootstrap-assisted inference procedure for generalized Grenander-type
estimators. The procedure relies on a carefully crafted, yet automatic,
transformation of the estimator. Moreover, our proposed method can be made
``flatness robust'' in the sense that it can be made adaptive to the (possibly
unknown) degree of flatness of the function of interest. The method requires
only the consistent estimation of a single scalar quantity, for which we
propose an automatic procedure based on numerical derivative estimation and the
generalized jackknife. Under random sampling, our inference method can be
implemented using a computationally attractive exchangeable bootstrap
procedure. We illustrate our methods with examples and we also provide a small
simulation study. The development of formal results is made possible by some
technical results that may be of independent interest.
