Deconvolution of distribution functions without integral transforms

Abstract

We study the recovery of the distribution function FX of a random variable X that is subject to an independent additive random error . To be precise, it is assumed that the target variable X is available only in the form of a blurred surrogate Y = X + . The distribution function FY then corresponds to the convolution of FX and F, so that the reconstruction of FX is some kind of deconvolution problem. Those have a long history in mathematics and various approaches have been proposed in the past. Most of them use integral transforms or matrix algorithms. The present article avoids these tools and is entirely confined to the domain of distribution functions. Our main idea relies on a transformation of a first kind to a second kind integral equation. Thereof, starting with a right-lateral discrete target and error variable, a representation for FX in terms of available quantities is obtained, which facilitates the unbiased estimation through a Y-sample. It turns out that these results even extend to cases in which X is not discrete. Finally, in a general setup, our approach gives rise to an approximation for FX as a certain Neumann sum. The properties of this sum are briefly examined theoretically and visually. The paper is concluded with a short discussion of operator theoretical aspects and an outlook on further research. Various plots underline our results and illustrate the capabilities of our functions with regard to estimation.

0

Turn this paper into a full lesson

ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…