Researcher profile

Philipp Wacker

Philipp Wacker contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

6 published item(s)

preprint2022arXiv

MAP estimators for nonparametric Bayesian inverse problems in Banach spaces

In order to rigorously define maximum-a-posteriori estimators for nonparametric Bayesian inverse problems for general Banach space valued parameters, we derive and prove certain previously postulated but unproven bounds on small ball probabilities. This allows us to prove existence of MAP estimators in the Banach space setting under very mild assumptions on the loglikelihood. As a similar statement so far (as far as the author is aware) only existed in the Hilbert space setting, this closes an important gap in the literature.

preprint2022arXiv

Nested sampling for physical scientists

We review Skilling's nested sampling (NS) algorithm for Bayesian inference and more broadly multi-dimensional integration. After recapitulating the principles of NS, we survey developments in implementing efficient NS algorithms in practice in high-dimensions, including methods for sampling from the so-called constrained prior. We outline the ways in which NS may be applied and describe the application of NS in three scientific fields in which the algorithm has proved to be useful: cosmology, gravitational-wave astronomy, and materials science. We close by making recommendations for best practice when using NS and by summarizing potential limitations and optimizations of NS.

preprint2022arXiv

The lion in the attic -- A resolution of the Borel--Kolmogorov paradox

The Borel--Kolmogorov paradox of conditioning with respect to events of prior probability zero has fascinated students and researchers since its discovery more than 100 years ago. Classical conditioning is only valid with respect to events of positive probability. If we ignore this constraint and condition on such sets, for example events of type $\{Y=y\}$ for a continuously distributed random variable $Y$, almost any probability measure can be chosen as the conditional measure on such sets. There have been numerous descriptions and explanations of the paradox' appearance in the setting of conditioning on a subset of probability zero. However, most treatments don't supply explicit instructions on how to avoid it. We propose to close this gap by defining a version of conditional measure which utilizes the Hausdorff measure. This makes the choice canonical in the sense that it only depends on the geometry of the space, thus removing any ambiguity. We describe the set of possible measures arising in the context of the Borel--Kolmogorov paradox and classify those coinciding with the canonical measure. The objective of this manuscript is to provide a manual for singular conditional probability: We give an explicit explanation in which settings ambiguity arises (and where not) and how to get rid of this ambiguity once and for all by a canonical choice.

preprint2020arXiv

On the Convergence of the Laplace Approximation and Noise-Level-Robustness of Laplace-based Monte Carlo Methods for Bayesian Inverse Problems

The Bayesian approach to inverse problems provides a rigorous framework for the incorporation and quantification of uncertainties in measurements, parameters and models. We are interested in designing numerical methods which are robust w.r.t. the size of the observational noise, i.e., methods which behave well in case of concentrated posterior measures. The concentration of the posterior is a highly desirable situation in practice, since it relates to informative or large data. However, it can pose a computational challenge for numerical methods based on the prior or reference measure. We propose to employ the Laplace approximation of the posterior as the base measure for numerical integration in this context. The Laplace approximation is a Gaussian measure centered at the maximum a-posteriori estimate and with covariance matrix depending on the logposterior density. We discuss convergence results of the Laplace approximation in terms of the Hellinger distance and analyze the efficiency of Monte Carlo methods based on it. In particular, we show that Laplace-based importance sampling and Laplace-based quasi-Monte-Carlo methods are robust w.r.t. the concentration of the posterior for large classes of posterior distributions and integrands whereas prior-based importance sampling and plain quasi-Monte Carlo are not. Numerical experiments are presented to illustrate the theoretical findings.

preprint2020arXiv

Pointwise defined version of conditional expectation with respect to a random variable

It is often of interest to condition on a singular event given by a random variable, e.g. $\{Y=y\}$ for a continuous random variable $Y$. Conditional measures with respect to this event are usually derived as a special case of the conditional expectation with respect to the random variables generating sigma algebra. The existence of the latter is usually proven via a non-constructive measure-theoretic argument which yields an only almost-everywhere defined quantity. In particular, the quantity $\mathbb E[f|Y]$ is initially only defined almost everywhere and conditioning on $Y=y$ corresponds to evaluating $\mathbb E[f|Y=y] = \mathbb E[f|Y]{Y=y}$, which is not meaningful because of $\mathbb E[f|Y]$ not being well-defined on such singular sets. This problem is not addressed by the introduction of regular conditional distributions, either. On the other hand it can be shown that the naively computed conditional density $f_{Z|Y=y}(z)$ (which is given by the ratio of joint and marginal densities) is a version of the conditional distribution, i.e. $\mathbb E[\{Z\in B\}|Y=y] = \int_B f_{Z|Y=y}(z) dz$ and this density can indeed be evaluated pointwise in $y$. This mismatch between mathematical theory (which generates an object which cannot produce what we need from it) and practical computation via the conditional density is an unfortunate fact. Furthermore, the classical approach does not allow a pointwise definition of conditional expectations of the form $\mathbb E[f|Y=y]$, only of conditional distributions $\mathbb E[\{Z\in B\}|Y=y]$. We propose a (as far as the author is aware) little known approach to obtaining a pointwise defined version of conditional expectation by use of the Lebesgue-Besicovich lemma without the need of additional topological arguments which are necessary in the usual derivation.

preprint2016arXiv

Probabilistic Estimates of the Maximum Norm of Random Neumann Fourier Series

We study the maximum norm behavior of $L^2$-normalized random Fourier cosine series with a prescribed large wave number. Precise bounds of this type are an important technical tool in estimates for spinodal decomposition, the celebrated phase separation phenomenon in metal alloys. We derive rigorous asymptotic results as the wave number converges to infinity, and shed light on the behavior of the maximum norm for medium range wave numbers through numerical simulations. Finally, we develop a simplified model for describing the magnitude of extremal values of random Neumann Fourier series. The model describes key features of the development of maxima and can be used to predict them. This is achieved by decoupling magnitude and sign distribution, where the latter plays an important role for the study of the size of the maximum norm. Since we are considering series with Neumann boundary conditions, particular care has to be placed on understanding the behavior of the random sums at the boundary.