642

Deconvolution with unknown error distribution

Abstract

We consider the problem of estimating a density fXf_{X} using a sample Y1,...,YnY_{1},...,Y_{n} from fY=fXfϵf_{Y}=f_{X}*f_{\epsilon}, where fϵf_{\epsilon} is an unknown density function. We assume that an additional sample ϵ1,...,ϵm\epsilon_{1},...,\epsilon_{m} from fϵf_{\epsilon} is given. Estimators of fXf_{X} and its derivatives are constructed using nonparametric estimators of fYf_{Y} and fϵf_{\epsilon} and applying a spectral cut-off in the Fourier domain. In this paper the rate of convergence of the estimator is derived in the case of a known and an unknown density fϵf_{\epsilon} assuming that fXf_{X} belongs to a Sobolev space HpH_{p} and that the Fourier transform of fϵf_{\epsilon} descents polynomial, exponential or in some general form. It is shown that the proposed estimator is asymptotically optimal in a minimax sense if the density fϵf_{\epsilon} is known or has a Fourier transform with polynomial descent. Monte Carlo simulations demonstrate the reasonable performance of the estimator given an estimated error density fϵf_{\epsilon} compared with the estimator in the case when fϵf_{\epsilon} is known.

View on arXiv
Comments on this paper