Source: Ann. Statist.
Volume 38, Number 3
Using the asymptotical minimax framework, we examine convergence rates equivalency between a continuous functional deconvolution model and its real-life discrete counterpart over a wide range of Besov balls and for the L2-risk. For this purpose, all possible models are divided into three groups. For the models in the first group, which we call uniform, the convergence rates in the discrete and the continuous models coincide no matter what the sampling scheme is chosen, and hence the replacement of the discrete model by its continuous counterpart is legitimate. For the models in the second group, to which we refer as regular, one can point out the best sampling strategy in the discrete model, but not every sampling scheme leads to the same convergence rates; there are at least two sampling schemes which deliver different convergence rates in the discrete model (i.e., at least one of the discrete models leads to convergence rates that are different from the convergence rates in the continuous model). The third group consists of models for which, in general, it is impossible to devise the best sampling strategy; we call these models irregular.
We formulate the conditions when each of these situations takes place. In the regular case, we not only point out the number and the selection of sampling points which deliver the fastest convergence rates in the discrete model but also investigate when, in the case of an arbitrary sampling scheme, the convergence rates in the continuous model coincide or do not coincide with the convergence rates in the discrete model. We also study what happens if one chooses a uniform, or a more general pseudo-uniform, sampling scheme which can be viewed as an intuitive replacement of the continuous model. Finally, as a representative of the irregular case, we study functional deconvolution with a boxcar-like blurring function since this model has a number of important applications. All theoretical results presented in the paper are illustrated by numerous examples; many of which are motivated directly by a multitude of inverse problems in mathematical physics where one needs to recover initial or boundary conditions on the basis of observations from a noisy solution of a partial differential equation. The theoretical performance of the suggested estimator in the multichannel deconvolution model with a boxcar-like blurring function is also supplemented by a limited simulation study and compared to an estimator available in the current literature. The paper concludes that in both regular and irregular cases one should be extremely careful when replacing a discrete functional deconvolution model by its continuous counterpart.
Abramovich, F. and Silverman, B. W. (1998). Wavelet decomposition approaches to statistical inverse problems. Biometrika 85 115–129.
Antoniadis, A., Bigot, J. and Sapatinas, T. (2001). Wavelet estimators in nonparametric regression: A comparative simulation study. Journal of Statistical Software 6 Article 6.
Bender, C. M. and Orzag, S. A. (1978). Advanced Mathematical Methods for Scientists and Engineers. McGraw-Hill, New York.
Mathematical Reviews (MathSciNet): MR538168
Brown, L. D., Cai, T., Low, M. G. and Zhang, C.-H. (2002). Asymptotic equivalence theory for nonparametric regression with random design. Ann. Statist. 30 688–707.
Brown, L. D. and Low, M. G. (1996). Asymptotic equivalence of nonparametric regression and white noise. Ann. Statist. 24 2384–2398.
Casey, S. D. and Walnut, D. F. (1994). Systems of convolution equations, deconvolution, Shannon sampling, and the wavelet and Gabor transforms. SIAM Rev. 36 537–577.
Cavalier, L. and Raimondo, M. (2007). Wavelet deconvolution with noisy eigenvalues. IEEE Trans. Signal Process. 55 2414–2424.
Chesneau, C. (2008). Wavelet estimation via block thresholding: A minimax study under Lp-risk. Statist. Sinica 18 1007–1024.
Cirelson, B. S., Ibragimov, I. A. and Sudakov, V. N. (1976). Norm of Gaussian sample function. In Proceedings of the 3rd Japan-U.S.S.R. Symposium on Probability Theory. Lecture Notes in Math. 550 20–41. Springer, Berlin.
Mathematical Reviews (MathSciNet): MR458556
De Canditiis, D. and Pensky, M. (2004). Discussion on the meeting on “Statistical Approaches to Inverse Problems.” J. R. Stat. Soc. Ser. B Stat. Methodol. 66 638–640.
De Canditiis, D. and Pensky, M. (2006). Simultaneous wavelet deconvolution in periodic setting. Scand. J. Statist. 33 293–306.
Donoho, D. L. (1995). Nonlinear solution of linear inverse problems by wavelet-vaguelette decomposition. Appl. Comput. Harmon. Anal. 2 101–126.
Donoho, D. L. and Raimondo, M. (2004). Translation invariant deconvolution in a periodic setting. Int. J. Wavelets Multiresolut. Inf. Process. 14 415–432.
Golubev, G. (2004). The principle of penalized empirical risk in severely ill-posed problems. Probab. Theory Related Fields 130 18–38.
Golubev, G. K. and Khasminskii, R. Z. (1999). A statistical approach to some inverse problems for partial differential equations. Probl. Inf. Transm. 35 136–149.
Gradshtein, I. S. and Ryzhik, I. M. (1980). Tables of Integrals, Series, and Products. Academic Press, New York.
Härdle, W., Kerkyacharian, G., Picard, D. and Tsybakov, A. (1998). Wavelets, Approximation, and Statistical Applications. Lecture Notes in Statistics 129. Springer, New York.
Harsdorf, S. and Reuter, R. (2000). Stable deconvolution of noisy lidar signals. In Proceedings of EARSeL-SIG-Workshop LIDAR 16–17. Dresden/FRG.
Johnstone, I. M. (2002). Function Estimation in Gaussian Noise: Sequence Models. Unpublished Monograph. Available at http://www-stat.stanford.edu/~imj/.
Johnstone, I. M., Kerkyacharian, G., Picard, D. and Raimondo, M. (2004). Wavelet deconvolution in a periodic setting (with discussion). J. R. Stat. Soc. Ser. B Stat. Methodol. 66 547–573.
Johnstone, I. M. and Raimondo, M. (2004). Periodic boxcar deconvolution and Diophantine approximation. Ann. Statist. 32 1781–1804.
Kalifa, J. and Mallat, S. (2003). Thresholding estimators for linear inverse problems and deconvolutions. Ann. Statist. 31 58–109.
Kerkyacharian, G., Picard, D. and Raimondo, M. (2007). Adaptive boxcar deconvolution on full Lebesgue measure sets. Statist. Sinica 7 317–340.
Kolaczyk, E. D. (1994). Wavelet methods for the inversion of certain homogeneous linear operators in the presence of noisy data. Ph.D. dissertation, Dept. Statistics, Stanford Univ.
Lang, S. (1966). Introduction to Diophantine Approximations. Springer, New York.
Lattes, R. and Lions, J. L. (1967). Methode de Quasi-Reversibilite et Applications. Travoux et Recherche Mathematiques 15. Dunod, Paris.
Mathematical Reviews (MathSciNet): MR232549
Meyer, Y. (1992). Wavelets and Operators. Cambridge Univ. Press, Cambridge.
Müller, H. G. and Stadmüller, U. (1987). Variable bandwidth kernel estimators of regression curves. Ann. Statist. 15 182–201.
Neelamani, R., Choi, H. and Baraniuk, R. (2004). Forward: Fourier-wavelet regularized deconvolution for ill-conditioned systems. IEEE Trans. Signal Process. 52 418–433.
Park, Y. J., Dho, S. W. and Kong, H. J. (1997). Deconvolution of long-pulse lidar signals with matrix formulation. Appl. Optics 36 5158–5161.
Pensky, M. and Sapatinas, T. (2009a). Functional deconvolution in a periodic case: Uniform case. Ann. Statist. 37 73–104.
Pensky, M. and Sapatinas, T. (2009b). Diophantine approximation and the problem of estimation of the initial speed of a wave on a finite interval in a stochastic setting. Technical Report TR-07-2009, Dept. Mathematics and Statistics, Univ. Cyprus.
Petsa, A. and Sapatinas, T. (2009). Minimax convergence rates under the Lp-risk in the functional deconvolution model. Statist. Probab. Lett. 79 1568–1576. [Erratum: Statist. Probab. Lett. 79 1890 (2009).]
Reiss, M. (2008). Asymptotic equivalence for nonparametric regression with multivariate and random design. Ann. Statist. 36 1957–1982.
Schmidt, W. (1980). Diophantine Approximation. Lecture Notes in Math. 785. Springer, Berlin.
Mathematical Reviews (MathSciNet): MR568710
Willer, T. (2005). Deconvolution in white noise with a random blurring function. Preprint. Avaiable at arXiv:math/0505142v1.