Related papers: Peak Inference for Gaussian Random Fields on a Lat…
We study the peak height distribution of certain non-stationary Gaussian random fields. The explicit peak height distribution of smooth, non-stationary Gaussian processes in 1D with general covariance is derived. The formula is determined…
We present exact formulas for both the expected number and the height distribution of local maxima (peaks) in two distinct categories of smooth, non-centered Gaussian fields: (i) nonstationary Gaussian processes and (ii) stationary planar…
Let $\{f(t): t\in T\}$ be a smooth Gaussian random field over a parameter space $T$, where $T$ may be a subset of Euclidean space or, more generally, a Riemannian manifold. For any local maximum of $f(t)$ located at $t_0$ in the interior of…
We study the probability distribution of the maximum $M_S $ of a smooth stationary Gaussian field defined on a fractal subset $S$ of $\R^n$. Our main result is the equivalent of the asymptotic behavior of the tail of the distribution…
This paper proposes a new approach for Bayesian and maximum likelihood parameter estimation for stationary Gaussian processes observed on a large lattice with missing values. We propose an MCMC approach for Bayesian inference, and a Monte…
The explicit formulae for the height distribution and expected number of local maxima have been obtained for isotropic Gaussian random fields on certain low-dimensional Euclidean space or low-dimensional spheres.
Peak inference concerns the use of local maxima ("peaks") of a noisy random field to detect and localize regions where underlying signal is present. We propose a peak inference method that first subjects observed peaks to a significance…
Our focus is on the design and analysis of efficient Monte Carlo methods for computing tail probabilities for the suprema of Gaussian random fields, along with conditional expectations of functionals of the fields given the existence of…
Motivated by the problem of testing for the existence of a signal of known parametric structure and unknown ``location'' (as explained below) against a noisy background, we obtain for the maximum of a centered, smooth random field an…
We study the probability distribution $F(u)$ of the maximum of smooth Gaussian fields defined on compact subsets of $\R^d$ having some geometric regularity. Our main result is a general formula for the density of $F$. Even though this is an…
Random fields in nature often have, to a good approximation, Gaussian characteristics. We present the mathematical framework for a new and simple method for investigating the non-Gaussian contributions, based on counting the maxima and…
In this paper we develop a novel inferential approach based on geometric records for estimating the tail index of heavy-tailed distributions. We construct a maximum likelihood estimator for the Pareto model and establish its strong…
We introduce efficient Markov chain Monte Carlo methods for inference and model determination in multivariate and matrix-variate Gaussian graphical models. Our framework is based on the G-Wishart prior for the precision matrix associated…
Markov chain Monte Carlo methods are often deemed too computationally intensive to be of any practical use for big data applications, and in particular for inference on datasets containing a large number $n$ of individual data points, also…
Sampling from a lattice Gaussian distribution is emerging as an important problem in various areas such as coding and cryptography. The default sampling algorithm --- Klein's algorithm yields a distribution close to the lattice Gaussian…
We study power approximation formulas for peak detection using Gaussian random field theory. The approximation, based on the expected number of local maxima above the threshold $u$, $\mathbb{E}[M_u]$, is proved to work well under three…
Stochastic Differential Equations (SDEs) are used as statistical models in many disciplines. However, intractable likelihood functions for SDEs make inference challenging, and we need to resort to simulation-based techniques to estimate and…
Fractional Gaussian fields provide a rich class of spatial models and have a long history of applications in multiple branches of science. However, estimation and inference for fractional Gaussian fields present significant challenges. This…
In this paper we consider Bayesian estimation for the parameters of inverse Gaussian distribution. Our emphasis is on Markov Chain Monte Carlo methods. We provide complete implementation of the Gibbs sampler algorithm. Assuming an…
We study some aspects of a Monte Carlo method invented by Maggs and Rossetto for simulating systems of charged particles. It has the feature that the discretized electric field is updated locally when charges move. Results of simulations of…