























In this paper, we propose a novel Bayesian approach for nonparametric estimation in Wicksell's problem. This has important applications in astronomy for estimating the distribution of the positions of the stars in a galaxy given projected stellar positions and in materials science to determine the 3D microstructure of a material, using its 2D cross sections. We deviate from the classical Bayesian nonparametric approach, which would place a Dirichlet Process (DP) prior on the distribution function of the unobservables, by directly placing a DP prior on the distribution function of the observables. Our method offers computational simplicity due to the conjugacy of the posterior and allows for asymptotically efficient estimation by projecting the posterior onto the \( \mathbb{L}_2 \) subspace of increasing, right-continuous functions. Indeed, the resulting Isotonized Inverse Posterior (IIP) satisfies a Bernstein--von Mises (BvM) phenomenon with minimax asymptotic variance \( g_0(x)/2γ\), where \( γ> 1/2 \) reflects the degree of Hölder continuity of the true cdf at \( x \). Since the IIP gives automatic uncertainty quantification, it eliminates the need to estimate \( γ\). Our results provide the first semiparametric Bernstein--von Mises theorem for projection-based posteriors with a DP prior in inverse problems.
此内容由惯性聚合(RSS阅读器)自动聚合整理,仅供阅读参考。 原文来自 — 版权归原作者所有。