Local polynomial regression with selection biased data

Let Y and X be real- and Rd-valued random variables. We consider the estimation of the nonparametric regression function m(x) = E(Y/X = x) when s greater than or equal to 1 independent selection-biased samples of (Y,X) are observed. This sampling scheme, which arises naturally in biological and epidemiological studies and many other fields, includes stratified samples, length-biased samples and other weighted distributions. A class of local polynomial estimators of m(x) is derived by smoothing Vardi's nonparametric maximum likelihood estimator of the underlying distribution function. Large sample properties, such as asymptotic distributions and asymptotic mean squared risks, are derived explicitly Unlike local polynomial regression with i.i.d. direct samples, we show here that kernel choices are important and optimal kernel functions may be asymmetric and discontinuous when the weight functions of the biased samples have jumps. A cross-validation criterion is proposed for the selection of data-driven bandwidths. Through a simple comparison, we show that our estimators are superior to other intuitive estimators of m(x).