Related papers: Simulation of hitting times for Bessel processes w…
Diffusion processes arise in many fields, and so simulating the path of a diffusion is an important problem. It is usually necessary to make some sort of approximation via model-discretization, but a recently introduced class of algorithms,…
We describe a new algorithm for the numerical simulation of quantum spin and boson systems. The method is based on the Trotter decomposition in imaginary time and a decoupling by auxiliary Ising spins. It can be applied, in principle, to…
While an important topic in practice, the estimation of the number of non-noise components in blind source separation has received little attention in the literature. Recently, two bootstrap-based techniques for estimating the dimension…
Computer experiments are becoming increasingly important in scientific investigations. In the presence of uncertainty, analysts employ probabilistic sensitivity methods to identify the key-drivers of change in the quantities of interest.…
We study statistical inference on unit roots and cointegration for time series in a Hilbert space. We develop statistical inference on the number of common stochastic trends embedded in the time series, i.e., the dimension of the…
We present the mathematical basis of a new approach to the analysis of temporal coding. The foundation of the approach is the construction of several families of novel distances (metrics) between neuronal impulse trains. In contrast to most…
Simulating samples from arbitrary probability distributions is a major research program of statistical computing. Recent work has shown promise in an old idea, that sampling from a discrete distribution can be accomplished by perturbing and…
Numerical methods for the simulation of cavitation processes have been developed for more than 50 years. The rich variety of physical phenomena triggered by the collapse of a bubble has several applications in medicine and environmental…
We consider the one-dimensional squared Bessel process given by the stochastic differential equation (SDE) \begin{align*} dX_t = 1\,dt + 2\sqrt{X_t}\,dW_t, \quad X_0=x_0, \quad t\in[0,1], \end{align*} and study strong (pathwise)…
The Poisson process of order $i$ is a weighted sum of independent Poisson processes and is used to model the flow of clients in different services. In the paper below we study some extensions of this process, for different forms of the…
We propose a novel mesh refinement scheme based on signal processing for boundary integral simulations of inviscid droplet dynamics with axial symmetry. A key idea is to directly access the Fourier coefficients of a principal curvature as a…
This article deals with the numerical resolution of backward stochastic differential equations. Firstly, we consider a rather general case where the filtration is generated by a Brownian motion and a Poisson random measure. We provide a…
We consider discrete time dynamical systems and show the link between Hitting Time Statistics (the distribution of the first time points land in asymptotically small sets) and Extreme Value Theory (distribution properties of the partial…
In the usual Bayesian setting, a full probabilistic model is required to link the data and parameters, and the form of this model and the inference and prediction mechanisms are specified via de Finetti's representation. In general, such a…
Approximating integrals is a fundamental task in probability theory and statistical inference, and their applied fields of signal processing, and Bayesian learning, as soon as expectations over probability distributions must be computed…
We provide asymptotic results and develop high frequency statistical procedures for time-changed L\'evy processes sampled at random instants. The sampling times are given by first hitting times of symmetric barriers whose distance with…
We present a method to compute the magnetic moment of a bulk, finite-size, three-dimensional, anisotropic superconductor. Our numerically implemented perturbative procedure is based on a solution of the nonlinear Maxwell- London equations,…
The Bessel process models the local eigenvalue statistics near $0$ of certain large positive definite matrices. In this work, we consider the probability \begin{align*} \mathbb{P}\Big( \mbox{there are no points in the Bessel process on }…
In this paper, we construct the Bessel line ensemble, a countable collection of continuous random curves. This line ensemble is stationary under horizontal shifts with the Bessel point process as its one-time marginal. Its finite…
We consider systems of multiple Brownian particles in one dimension that repel mutually via a logarithmic potential on the real line, more specifically the Dyson model. These systems are characterized by a parameter that controls the…