Most results in nonparametric regression theory are developed only for the
case of additive noise. In such a setting many smoothing techniques including
wavelet thresholding methods have been developed and shown to be highly
adaptive. In this paper we consider nonparametric regression in exponential
families with the main focus on the natural exponential families with a
quadratic variance function, which include, for example, Poisson regression,
binomial regression and gamma regression. We propose a unified approach of
using a mean-matching variance stabilizing transformation to turn the
relatively complicated problem of nonparametric regression in exponential
families into a standard homoscedastic Gaussian regression problem. Then in
principle any good nonparametric Gaussian regression procedure can be applied
to the transformed data. To illustrate our general methodology, in this paper
we use wavelet block thresholding to construct the final estimators of the
regression function. The procedures are easily implementable. Both theoretical
and numerical properties of the estimators are investigated. The estimators are
shown to enjoy a high degree of adaptivity and spatial adaptivity with
near-optimal asymptotic performance over a wide range of Besov spaces. The
estimators also perform well numerically.Comment: Published in at http://dx.doi.org/10.1214/09-AOS762 the Annals of
Statistics (http://www.imstat.org/aos/) by the Institute of Mathematical
Statistics (http://www.imstat.org