Non-parametric confidence bands in deconvolution density estimation
Uniform confidence bands for densities "f" via non-parametric kernel estimates were first constructed by Bickel and Rosenblatt. In this paper this is extended to confidence bands in the deconvolution problem "g"="f"*"ψ" for an ordinary smooth error density "ψ". Under certain regularity conditions, we obtain asymptotic uniform confidence bands based on the asymptotic distribution of the maximal deviation ("L"<sub>∞</sub>-distance) between a deconvolution kernel estimator <formula format="inline"><file name="rssb_599_mu1.gif" type="gif" /></formula> and "f". Further consistency of the simple non-parametric bootstrap is proved. For our theoretical developments the bias is simply corrected by choosing an undersmoothing bandwidth. For practical purposes we propose a new data-driven bandwidth selector that is based on heuristic arguments, which aims at minimizing the "L"<sub>∞</sub>-distance between <formula format="inline"><file name="rssb_599_mu2.gif" type="gif" /></formula> and "f". Although not constructed explicitly to undersmooth the estimator, a simulation study reveals that the bandwidth selector suggested performs well in finite samples, in terms of both area and coverage probability of the resulting confidence bands. Finally the methodology is applied to measurements of the metallicity of local F and G dwarf stars. Our results confirm the 'G dwarf problem', i.e. the lack of metal poor G dwarfs relative to predictions from 'closed box models' of stellar formation. Copyright 2007 Royal Statistical Society.
