BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Robert E. Gaunt

Robert E. Gaunt contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.PRmath.STStatistics Theorymath.CAComputationMethodology

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
6topics
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

14 published item(s)

preprint2026arXiv

Integrals of products of four modified Bessel functions

We evaluate definite integrals involving the product of four modified Bessel functions of the first and second kind and a power function. We provide general formulas expressed in terms of the Meijer $G$-function and generalized hypergeometric and Lauricella $F_C$ functions, and study a number of special cases in which the integrals can be evaluated in terms of simpler special functions or indeed take an elementary form. As a consequence, we deduce some new formulas for definite integrals of products of four Airy functions.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Stein's method for the matrix normal distribution

This work presents the first systematic development of Stein's method for matrix distributions. We establish the basic essential ingredients of Stein's method for matrix normal approximation: we derive a generator-based Stein identity from a matrix Ornstein--Uhlenbeck diffusion with two-sided scales, provide an explicit semigroup representation for the solution of the Stein equation, and obtain regularity estimates for the solution. The new methodology is illustrated with three statistical applications, these being smooth Wasserstein distance bounds to quantify the matrix central limit theorem, a Wasserstein distance bound for the matrix normal approximation of the centered matrix $T$ distribution, and the derivation of Stein's method-of-moments estimators for scale parameters of the matrix normal distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Bounding Kolmogorov distances through Wasserstein and related integral probability metrics

We establish general upper bounds on the Kolmogorov distance between two probability distributions in terms of the distance between these distributions as measured with respect to the Wasserstein or smooth Wasserstein metrics. These bounds generalise existing results from the literature. To illustrate the broad applicability of our general bounds, we apply them to extract Kolmogorov distance bounds from multivariate normal, beta and variance-gamma approximations that have been established in the Stein's method literature.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Bounds for an integral involving the modified Lommel function of the first kind

Simple upper and lower bounds are established for the integral $\int_0^x\mathrm{e}^{-βu}u^νt_{μ,ν}(u)\,\mathrm{d}u$, where $x>0$, $0<β<1$, $μ+ν>-2$, $μ-ν\geq-3$, and $t_{μ,ν}(x)$ is the modified Lommel function of the first kind. Our bounds complement and improve on existing bounds for this integral, by either being sharper or increasing the range of validity. Our bounds also generalise recent bounds for an integral involving the modified Struve function of the first kind, and in some cases more a direct approach lead to sharper bounds when our general bounds are specialised to the modified Struve case.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Bounds for the chi-square approximation of Friedman&#39;s statistic by Stein&#39;s method

Friedman&#39;s chi-square test is a non-parametric statistical test for $r\geq2$ treatments across $n\ge1$ trials to assess the null hypothesis that there is no treatment effect. We use Stein&#39;s method with an exchangeable pair coupling to derive an explicit bound on the distance between the distribution of Friedman&#39;s statistic and its limiting chi-square distribution, measured using smooth test functions. Our bound is of the optimal order $n^{-1}$, and also has an optimal dependence on the parameter $r$, in that the bound tends to zero if and only if $r/n\rightarrow0$. From this bound, we deduce a Kolmogorov distance bound that decays to zero under the weaker condition $r^{1/2}/n\rightarrow0$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On algebraic Stein operators for Gaussian polynomials

The first essential ingredient to build up Stein&#39;s method for a continuous target distribution is to identify a so-called \textit{Stein operator}, namely a linear differential operator with polynomial coefficients. In this paper, we introduce the notion of \textit{algebraic} Stein operators (see Definition \ref{def:algebraic-Stein-Operator}), and provide a novel algebraic method to find \emph{all} the algebraic Stein operators up to a given order and polynomial degree for a target random variable of the form $Y=h(X)$, where $X=(X_1,\dots, X_d)$ has i.i.d$.$ standard Gaussian components and $h\in \mathbb{K}[X]$ is a polynomial with coefficients in the ring $\mathbb{K}$. Our approach links the existence of an algebraic Stein operator with \textit{null controllability} of a certain linear discrete system. A \texttt{MATLAB} code checks the null controllability up to a given finite time $T$ (the order of the differential operator), and provides all \textit{null control} sequences (polynomial coefficients of the differential operator) up to a given maximum degree $m$. This is the first paper that connects Stein&#39;s method with computational algebra to find Stein operators for highly complex probability distributions, such as $H_{20}(X_1)$, where $H_p$ is the $p$-th Hermite polynomial. Some examples of Stein operators for $H_p(X_1)$, $p=3,4,5,6$, are gathered in the Appendix and many other examples are given in the Supplementary Information.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Stein factors for variance-gamma approximation in the Wasserstein and Kolmogorov distances

We obtain new bounds for the solution of the variance-gamma (VG) Stein equation that are of the correct form for approximations in terms of the Wasserstein and Kolmorogorov metrics. These bounds hold for all parameters values of the four parameter VG class. As an application we obtain explicit Wasserstein and Kolmogorov distance error bounds in a six moment theorem for VG approximation of double Wiener-Itô integrals.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Stein&#39;s Method Meets Computational Statistics: A Review of Some Recent Developments

Stein&#39;s method compares probability distributions through the study of a class of linear operators called Stein operators. While mainly studied in probability and used to underpin theoretical statistics, Stein&#39;s method has led to significant advances in computational statistics in recent years. The goal of this survey is to bring together some of these recent developments and, in doing so, to stimulate further research into the successful field of Stein&#39;s method and statistics. The topics we discuss include tools to benchmark and compare sampling methods such as approximate Markov chain Monte Carlo, deterministic alternatives to sampling methods, control variate techniques, parameter estimation and goodness-of-fit testing.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The basic distributional theory for the product of zero mean correlated normal random variables

The product of two zero mean correlated normal random variables, and more generally the sum of independent copies of such random variables, has received much attention in the statistics literature and appears in many application areas. However, many important distributional properties are yet to be recorded. This review paper fills this gap by providing the basic distributional theory for the sum of independent copies of the product of two zero mean correlated normal random variables. Properties covered include probability and cumulative distribution functions, generating functions, moments and cumulants, mode and median, Stein characterisations, representations in terms of other random variables, and a list of related distributions. We also review how the product of two zero mean correlated normal random variables arises naturally as a limiting distribution, with an example given for the distributional approximation of double Wiener-Itô integrals.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Bounds for modified Lommel functions of the first kind and their ratios

The modified Lommel function $t_{μ,ν}(x)$ is an important special function, but to date there has been little progress on the problem of obtaining functional inequalities for $t_{μ,ν}(x)$. In this paper, we advance the literature substantially by obtaining a simple two-sided inequality for the ratio $t_{μ,ν}(x)/t_{μ-1,ν-1}(x)$ in terms of the ratio $I_ν(x)/I_{ν-1}(x)$ of modified Bessel functions of the first kind, thereby allowing one to exploit the extensive literature on bounds for this ratio. We apply this result to obtain two-sided inequalities for the condition numbers $xt_{μ,ν}&#39;(x)/t_{μ,ν}(x)$, the ratio $t_{μ,ν}(x)/t_{μ,ν}(y)$ and the modified Lommel function $t_{μ,ν}(x)$ itself that are given in terms of $xI_ν&#39;(x)/I_ν(x)$, $I_ν(x)/I_ν(y)$ and $I_ν(x)$, respectively. The bounds obtained in this paper are quite accurate and often tight in certain limits. As an important special case we deduce bounds for modified Struve functions of the first kind and their ratios, some of which are new, whilst others extend the range of validity of some results given in the recent literature.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On bounds for the mode and median of the generalized hyperbolic and related distributions

Except for certain parameter values, a closed form formula for the mode of the generalized hyperbolic (GH) distribution is not available. In this paper, we exploit results from the literature on modified Bessel functions and their ratios to obtain simple but tight two-sided inequalities for the mode of the GH distribution for general parameter values. As a special case, we deduce tight two-sided inequalities for the mode of the variance-gamma (VG) distribution, and through a similar approach we also obtain tight two-sided inequalities for the mode of the McKay Type I distribution. The analogous problem for the median is more challenging, but we conjecture some monotonicity results for the median of the VG and McKay Type I distributions, from we which we conjecture some tight two-sided inequalities for their medians. Numerical experiments support these conjectures and also lead us to a conjectured tight lower bound for the median of the GH distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Stein&#39;s method for functions of multivariate normal random variables

By the continuous mapping theorem, if a sequence of $d$-dimensional random vectors $(\mathbf{W}_n)_{n\geq1}$ converges in distribution to a multivariate normal random variable $Σ^{1/2}\mathbf{Z}$, then the sequence of random variables $(g(\mathbf{W}_n))_{n\geq1}$ converges in distribution to $g(Σ^{1/2}\mathbf{Z})$ if $g:\mathbb{R}^d\rightarrow\mathbb{R}$ is continuous. In this paper, we develop Stein&#39;s method for the problem of deriving explicit bounds on the distance between $g(\mathbf{W}_n)$ and $g(Σ^{1/2}\mathbf{Z})$ with respect to smooth probability metrics. We obtain several bounds for the case that the $j$-component of $\mathbf{W}_n$ is given by $W_{n,j}=\frac{1}{\sqrt{n}}\sum_{i=1}^nX_{ij}$, where the $X_{ij}$ are independent. In particular, provided $g$ satisfies certain differentiability and growth rate conditions, we obtain an order $n^{-(p-1)/2}$ bound, for smooth test functions, if the first $p$ moments of the $X_{ij}$ agree with those of the normal distribution. If $p$ is an even integer and $g$ is an even function, this convergence rate can be improved further to order $n^{-p/2}$. These convergence rates are shown to be of optimal order. We apply our general bounds to some examples, which include the distributional approximation of asymptotically chi-square distributed statistics; the approximation of expectations of smooth functions of binomial and Poisson random variables; rates of convergence in the delta method; and a quantitative variance-gamma approximation of the $D_2^*$ statistic for alignment-free sequence comparison in the case of binary sequences.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Multivariate normal approximation of the maximum likelihood estimator via the delta method

We use the delta method and Stein&#39;s method to derive, under regularity conditions, explicit upper bounds for the distributional distance between the distribution of the maximum likelihood estimator (MLE) of a $d$-dimensional parameter and its asymptotic multivariate normal distribution. Our bounds apply in situations in which the MLE can be written as a function of a sum of i.i.d. $t$-dimensional random vectors. We apply our general bound to establish a bound for the multivariate normal approximation of the MLE of the normal distribution with unknown mean and variance.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Wasserstein and Kolmogorov error bounds for variance-gamma approximation via Stein&#39;s method I

The variance-gamma (VG) distributions form a four parameter family that includes as special and limiting cases the normal, gamma and Laplace distributions. Some of the numerous applications include financial modelling and approximation on Wiener space. Recently, Stein&#39;s method has been extended to the VG distribution. However, technical difficulties have meant that bounds for distributional approximations have only been given for smooth test functions (typically requiring at least two derivatives for the test function). In this paper, which deals with symmetric variance-gamma (SVG) distributions, and a companion paper \cite{gaunt vgii}, which deals with the whole family of VG distributions, we address this issue. In this paper, we obtain new bounds for the derivatives of the solution of the SVG Stein equation, which allow for approximations to be made in the Kolmogorov and Wasserstein metrics, and also introduce a distributional transformation that is natural in the context of SVG approximation. We apply this theory to obtain Wasserstein or Kolmogorov error bounds for SVG approximation in four settings: comparison of VG and SVG distributions, SVG approximation of functionals of isonormal Gaussian processes, SVG approximation of a statistic for binary sequence comparison, and Laplace approximation of a random sum of independent mean zero random variables.

0 citations0 reviewsNo code linkedOpen access