Semiparametric Uncertainty Quantification via Isotonized Posterior for Deconvolutions

TL;DR

This paper introduces a fast, nonparametric Bayesian method for uncertainty quantification in deconvolution problems that avoids nuisance parameter estimation by projecting the posterior and isotonic adjustments, ensuring valid credible sets.

Contribution

It proposes a novel, computationally efficient Bayesian approach using posterior projection and isotonization to provide valid uncertainty quantification without nuisance parameter estimation in deconvolution.

Findings

Posterior credible sets achieve asymptotic frequentist coverage.

Method is computationally fast and robust across various noise distributions.

Avoids the need for estimating nuisance parameters in uncertainty quantification.

Abstract

We address the problem of uncertainty quantification for the deconvolution model \(Z = X + Y\), where \(X\) and \(Y\) are nonnegative random variables and the goal is to estimate the signal's distribution of \(X \sim F_0\) supported on~\([0,\infty)\), from observations where the noise distribution is known. Existing frequentist methods often produce confidence intervals for $F_0(x)$ that depend on unknown nuisance parameters, such as the density of \(X\) and its derivative, which are difficult to estimate in practice. This paper introduces a novel and computationally efficient nonparametric Bayesian approach, based on projecting the posterior, to overcome this limitation. Our method leverages the solution \(p\) to a specific Volterra integral equation as in \cite{74}, which relates the cumulative distribution function (CDF) of the signal, \(F_0\), to the distribution of the observables. We place a Dirichlet Process prior directly on the distribution of the observed data $Z$, yielding a simple, conjugate posterior. To ensure the resulting estimates for \(F_0\) are valid CDFs, we isotonize posterior draws taking the Greatest Convex Majorant of the primitive of the posterior draws and defining what we term the Isotonic Inverse Posterior. We show that this framework yields posterior credible sets for \(F_0\) that are not only computationally fast to generate but also possess asymptotically correct frequentist coverage after a straightforward recalibration technique for the so-called Bayes Chernoff distribution introduced in \cite{54}. Our approach thus does not require the estimation of nuisance parameters to deliver uncertainty quantification for the parameter of interest $F_0(x)$. The practical effectiveness and robustness of the method are demonstrated through a simulation study with various noise distributions for $Y$.