A Proof of the Explicit Formula for Product Moments of Multivariate Gaussian Random Variables

Iickho Song

DOI: https://doi.org/10.48550/arXiv.1705.00163

2017-04-29

Abstract:A detailed proof of a recent result on explicit formulae for the product moments $E \left \{ X_1^{a_1} X_2^{a_2} \cdots X_n^{a_n}\right \}$ of multivariate Gaussian random variables is provided in this note.

Probability,Statistics Theory,Computation

What problem does this paper attempt to address?

An explicit Wishart moment formula for the product of two disjoint principal minors

Christian Genest,Frédéric Ouimet,Donald Richards

2024-09-23

Abstract:This paper provides the first explicit formula for the expectation of the product of two disjoint principal minors of a Wishart random matrix, solving a part of a broader problem put forth by Samuel S. Wilks in 1934 in the Annals of Mathematics. The proof makes crucial use of hypergeometric functions of matrix argument and their Laplace transforms. Additionally, a Wishart generalization of the Gaussian product inequality conjecture is formulated and a stronger quantitative version is proved to hold in the case of two minors.

Probability,Statistics Theory
A fast algorithm for computing product moments of multivariate normal random variables

Raymond Kan Jiening Pan a Rotman School of Management,University of Toronto,Toronto,Ontario,Canadab School of Finance,Nankai University,Tianjin,China

DOI: https://doi.org/10.1080/03610918.2024.2360136

2024-06-04

Communications in Statistics - Simulation and Computation

Abstract:We provide a simple identity that decomposes a product moment of multivariate normal random variables as a sum of various products of univariate moments of one of the random variables and multivariate moments of the other random variables. The new identity allows for much faster computation of the product moments of multivariate normal random variables than existing methods.

statistics & probability
A simple proof for the almost sure convergence of the largest singular value of a product of Gaussian matrices

Thiziri Nait Saada,Alireza Naderi

2024-09-30

Abstract:Let $m \geq 1$ and consider the product of $m$ independent $n \times n$ Gaussian matrices $\mathbf{W} = \mathbf{W}_1 \dots \mathbf{W}_m$, each $\mathbf{W}_{i}$ with i.i.d. normalised $\mathcal{N}(0, n^{-1/2})$ entries. It is shown in Penson et al. (2011) that the empirical distribution of the squared singular values of $\mathbf{W}$ converges to a deterministic distribution compactly supported on $[0, u_m]$, where $u_m = \frac{{(m+1)}^{m+1}}{m^m}$. This generalises the well-known case of $m=1$, corresponding to the Marchenko-Pastur distribution for square matrices. Moreover, for $m=1$, it was first shown by \cite{Geman} that the largest squared singular value almost surely converges to the right endpoint (aka ``soft edge'') of the support, i.e. $s_1^2(\mathbf{W}) \xrightarrow{a.s.} u_1$. Herein, we present a proof for the general case $s_1^2(\mathbf{W}) \xrightarrow{a.s.} u_m$ when $m\geq 1$. Although we do not claim novelty for our result, the proof is simple and does not require familiarity with modern techniques of free probability.

Probability
Local Times of N-parameter Gaussian Processes

ZY Lin,ZM Cheng

DOI: https://doi.org/10.1080/10020070512331341850

IF: 4.269

2005-01-01

Progress in Natural Science Materials International

Abstract:An N-parameter Gaussian stationary process X = X(t) : tR(+)(N) is introduced and the existence and joint continuity of its local times is presented. And the moments of local times are estimated. Furthermore moduli of continuity and large increment results for the local times are established.
Using Sums-of-Squares to Prove Gaussian Product Inequalities

Oliver Russell,Wei Sun

DOI: https://doi.org/10.48550/arXiv.2205.02127

2022-10-14

Abstract:The long-standing Gaussian product inequality (GPI) conjecture states that $E [\prod_{j=1}^{n}X_j^{2m_j}]\geq\prod_{j=1}^{n}E[X_j^{2m_j}]$ for any centered Gaussian random vector $(X_1,\dots,X_n)$ and $m_1,\dots,m_n\in\mathbb{N}$. In this paper, we describe a computational algorithm involving sums-of-squares representations of multivariate polynomials that can be used to resolve the GPI conjecture. To exhibit the power of this novel method, we apply it to prove two new GPIs: $E[X_1^{2m_1}X_2^{6}X_3^{4}]\ge E[X_1^{2m_1}]E[X_2^{6}]E[X_3^{4}]$ and $E[X_1^{2m_1}X_2^{2}X_3^{2}X_4^{2}]\ge E[X_1^{2m_1}]E[X_2^{2}]E[X_3^{2}]E[X_4^{2}]$.

Probability
Explicit expressions for joint moments of $n$-dimensional elliptical distributions

Baishuai Zuo,Chuancun Yin,Narayanaswamy Balakrishnan

DOI: https://doi.org/10.48550/arXiv.2007.09349

2020-08-02

Abstract:Inspired by Stein's lemma, we derive two expressions for the joint moments of elliptical distributions. We use two different methods to derive $E[X_{1}^{2}f(\mathbf{X})]$ for any measurable function $f$ satisfying some regularity conditions. Then, by applying this result, we obtain new formulae for expectations of product of normally distributed random variables, and also present simplified expressions of $E[X_{1}^{2}f(\mathbf{X})]$ for multivariate Student-$t$, logistic and Laplace distributions.

Statistics Theory,Risk Management
A Note on the Limiting Mean Distribution of Singular Values for Products of Two Wishart Random Matrices

Lun Zhang

DOI: https://doi.org/10.1063/1.4818978

IF: 1.469

2013-01-01

Journal of Mathematical Physics

Abstract:The product of M complex random Gaussian matrices of size N has recently been studied by Akemann, Kieburg, and Wei. They showed that, for fixed M and N, the joint probability distribution for the squared singular values of the product matrix forms a determinantal point process with a correlation kernel determined by certain biorthogonal polynomials that can be explicitly constructed. We find that, in the case M = 2, the relevant biorthogonal polynomials are actually special cases of multiple orthogonal polynomials associated with Macdonald functions (modified Bessel functions of the second kind) which was first introduced by Van Assche and Yakubovich. With known results on asymptotic zero distribution of these polynomials and general theory on multiple orthogonal polynomial ensembles, it is then easy to obtain an explicit expression for the distribution of squared singular values for the product of two complex random Gaussian matrices in the limit of large matrix dimensions.
A Central Limit Theorem with Explicit Lyapunov Exponent and Variance for Products of Random Non-invertible Matrices

Audrey Benson,Hunter Gould,Phanuel Mariano,Grace Newcombe,Joshua Vaidman

DOI: https://doi.org/10.1007/s10955-024-03335-3

2024-10-02

Journal of Statistical Physics

Abstract:The theory of products of random matrices and Lyapunov exponents have been widely studied and applied in the fields of biology, dynamical systems, economics, engineering and statistical physics. We consider the product of an i.i.d. sequence of random non-invertible matrices with real entries. Given some mild moment assumptions we prove an explicit formula for the Lyapunov exponent and prove a central limit theorem with an explicit formula for the variance in terms of the entries of the matrices. We also give examples where exact values for the Lyapunov exponent and variance are computed. An important example where non-invertible matrices are essential is the random Hill's equation, which has numerous physical applications, including the astrophysical orbit problem.

physics, mathematical
Symmetric functions and the explicit moment problem for abelian groups

Roger Van Peski

2024-02-26

Abstract:Recently, Sawin and Wood (arXiv:math/2210.06279) proved a formula for the distribution of a random abelian group $G$ in terms of its $H$-moments $\mathbb{E} [\#\operatorname{Sur}(G,H)]$. We show that properties of Macdonald polynomials yield an alternate proof.

Combinatorics,Number Theory,Probability
Limit theorems for a strongly irreducible product of independent random matrices under optimal moment assumptions

Axel Péneau

2024-08-21

Abstract:Let $ \nu $ be a probability distribution over the linear semi-group $ \mathrm{End}(E) $ for $ E $ a finite dimensional vector space over a locally compact field. We assume that $ \nu $ is proximal, strongly irreducible and that $ \nu^{*n}\{0\}=0 $ for all integers $ n\in\mathbb{N} $. We consider the random sequence $ \overline\gamma_n := \gamma_0 \cdots \gamma_{n-1} $ for $ (\gamma_k)_{k \ge 0} $ independents of distribution law $ \nu $. We define the logarithmic singular gap as $ \mathrm{sqz} = \log\left( \frac{\mu_1}{\mu_2} \right) $ , where $ \mu_1 $ and $ \mu_2 $ are the two largest singular values. We show that $ (\mathrm{sqz}(\overline{\gamma}_n))_{n\in\mathbb{N}} $ escapes to infinity linearly and satisfies exponential large deviations estimates below its escape rate. With the same assumptions, we also show that the image of a generic line by $ \overline{\gamma}_n $ as well as its eigenspace of maximal eigenvalue both converge to the same random line $l_\infty $ at an exponential speed.If we moreover assume that the push-forward distribution $N(\nu)$ is $ \mathrm{L}^p $ for $ N:g\mapsto\log\left(\|g\|\|g^{-1}\|\right) $ and for some $ p\ge 1 $, then we show that $ \log|w(l_\infty)| $ is $ \mathrm{L}^p $ for all unitary linear form $ w $ and the logarithm of each coefficient of $ \overline{\gamma}_n $ is almost surely equivalent to the logarithm of the norm. To prove these results, we do not rely on any classical results for random products of invertible matrices with $ \mathrm{L}^1 $ moment assumption. Instead we describe an effective way to group the i.i.d factors into i.i.d random words that are aligned in the Cartan projection. We moreover have an explicit control over the moments.

Probability
A central limit theorem with explicit Lyapunov exponent and variance for products of $2\times2$ random non-invertible matrices

Audrey Benson,Hunter Gould,Phanuel Mariano,Grace Newcombe,Joshua Vaidman

2024-06-15

Abstract:The theory of products of random matrices and Lyapunov exponents have been widely studied and applied in the fields of biology, dynamical systems, economics, engineering and statistical physics. We consider the product of an i.i.d. sequence of $2\times 2$ random non-invertible matrices with real entries. Given some mild moment assumptions we prove an explicit formula for the Lyapunov exponent and prove a central limit theorem with an explicit formula for the variance in terms of the entries of the matrices. We also give examples where exact values for the Lyapunov exponent and variance are computed. An important example where non-invertible matrices are essential is the random Hill's equation, which has numerous physical applications, including the astrophysical orbit problem.

Probability,Mathematical Physics
On the moments of the variance-gamma distribution

Robert E. Gaunt

2023-06-04

Abstract:We obtain new closed-form formulas for the moments and absolute moments of the variance-gamma distribution. We thus deduce new formulas for the moments and absolute moments of the product of two correlated zero mean normal random variables.

Probability,Statistics Theory
A recursion formula for mixed trace moments of isotropic Wishart matrices and the Gaussian unitary/orthogonal ensembles

Ben Deitmar

2024-04-10

Abstract:Exact recursion formulas for mixed moments of four fundamental random matrix ensembles are derived. The reason such recursive formulas are possible is closely related to properties of polygon gluings studied by Harer and Zagier as well as Akhmedov and Shakirov. The proofs of the formulas are direct and written in such a way that they do not rely on understanding of polygon gluings.

Probability,Statistics Theory
The Gaussian product inequality conjecture for multinomial covariances

Frédéric Ouimet

2023-08-24

Abstract:In this short note, we find an equivalent combinatorial condition only involving finite sums under which a centered Gaussian random vector with multinomial covariance matrix satisfies the Gaussian product inequality (GPI) conjecture. These covariance matrices are relevant since their off-diagonal elements are negative, which is the hardest case to cover for the GPI conjecture, as mentioned by [Russell, O., \& Sun, W. 2022. Some new {G}aussian product inequalities. {\em J. Math. Anal. Appl.}, {\bf 515}(2), Paper No. 126439, 21 pp.].

Probability,Statistics Theory
Moment formulae for general point processes

Laurent Decreusefond,Ian Flint

DOI: https://doi.org/10.1016/j.jfa.2014.04.014

2012-11-21

Abstract:The goal of this paper is to generalize most of the moment formulae obtained in [Pri11]. More precisely, we consider a general point process \mu, and show that the relevant quantities to our problem are the so-called Papangelou intensities. Then, we show some general formulae to recover the moment of order n of the stochastic integral of a random process. We will use these extended results to study a random transformation of the point process.

Probability
Central and noncentral moments of the multivariate hypergeometric distribution

Frédéric Ouimet

2024-04-14

Abstract:In this short note, explicit formulas are developed for the central and noncentral moments of the multivariate hypergeometric distribution. A numerical implementation is provided in Mathematica for fast evaluations. This work complements the paper by Ouimet (2021), where analogous formulas were derived and implemented in Mathematica for the multinomial distribution.

Probability,Statistics Theory,Computation
Product formulas for multiple stochastic integrals associated with Lévy processes

Paolo Di Tella,Christel Geiss,Alexander Steinicke

DOI: https://doi.org/10.1007/s13348-024-00456-6

2024-11-08

Collectanea mathematica

Abstract:In the present paper, we obtain an explicit product formula for products of multiple integrals w.r.t. a random measure associated with a Lévy process. As a building block, we use a representation formula for products of martingales from a compensated-covariation stable family. This enables us to consider Lévy processes with both jump and Gaussian part. It is well known that for multiple integrals w.r.t. the Brownian motion such product formulas exist without further integrability conditions on the kernels. However, if a jump part is present, this is, in general, false. Therefore, we provide here sufficient conditions on the kernels which allow us to establish product formulas. As an application, we obtain explicit expressions for the expectation of products of iterated integrals, as well as for the moments and the cumulants for stochastic integrals w.r.t. the random measure. Based on these expressions, we show a central limit theorem for the long time behaviour of a class of stochastic integrals. Finally, we provide methods to calculate the number of summands in the product formula.

mathematics, applied
Explicit Formulas for the Fourth Moment of Mixed Exponential Sums

Chen Hua,Chen Jianhua,Cai Guangxing,Ai Xiaochuan,Zhang Silan

DOI: https://doi.org/10.1016/j.jnt.2012.10.004

IF: 0.679

2013-01-01

Journal of Number Theory

Abstract:We use elementary methods to obtain an explicit formula for the fourth power mean∑m=1q∑χmodq|∑a=1q′χ(a)e(mak+naq)|4, where e(y)=e2πiy, χ is a Dirichlet character modulo q and ∑′a=1q denotes the summation over all a with (a,q)=1. This extends a result of Z.F. Xu et al. by avoiding the restriction (k,q)=1.
Exact Distribution for the Product of Two Correlated Gaussian Random Variables.

Guolong Cui,Xianxiang Yu,Salvatore Iommelli,Lingjiang Kong

DOI: https://doi.org/10.1109/lsp.2016.2614539

2016-01-01

IEEE Signal Processing Letters

Abstract:This letter considers the distribution of product for two correlated real Gaussian random variables with nonzero means and arbitrary variances, which arises widely in radar and communication societies. We determine the exact probability density function (PDF) in terms of an infinite sum of modified Bessel functions of second kind, which includes some existent results, i.e., zero-means and/or independent variables, as special cases. Then, we study the approximation error and convergence rate when finite summations are exploited in practice. Finally, we evaluate the PDF behaviors of the derived expression as well as the Monte Carlo simulations.
On Wilks' joint moment formulas for embedded principal minors of Wishart random matrices

C. Genest,F. Ouimet,D. Richards

DOI: https://doi.org/10.1002/sta4.706

2024-06-14

Stat

Abstract:In 1934, the American statistician Samuel S. Wilks derived remarkable formulas for the joint moments of embedded principal minors of sample covariance matrices in multivariate Gaussian populations, and he used them to compute the moments of sample statistics in various applications related to multivariate linear regression. These important but little‐known moment results were extended in 1963 by the Australian statistician A. Graham Constantine using Bartlett's decomposition. In this note, a new proof of Wilks' results is derived using the concept of iterated Schur complements, thereby bypassing Bartlett's decomposition. Furthermore, Wilks' open problem of evaluating joint moments of disjoint principal minors of Wishart random matrices is related to the Gaussian product inequality conjecture.

statistics & probability

A Proof of the Explicit Formula for Product Moments of Multivariate Gaussian Random Variables

An explicit Wishart moment formula for the product of two disjoint principal minors

A fast algorithm for computing product moments of multivariate normal random variables

A simple proof for the almost sure convergence of the largest singular value of a product of Gaussian matrices

Local Times of N-parameter Gaussian Processes

Using Sums-of-Squares to Prove Gaussian Product Inequalities

Explicit expressions for joint moments of $n$-dimensional elliptical distributions

A Note on the Limiting Mean Distribution of Singular Values for Products of Two Wishart Random Matrices

A Central Limit Theorem with Explicit Lyapunov Exponent and Variance for Products of Random Non-invertible Matrices

Symmetric functions and the explicit moment problem for abelian groups

Limit theorems for a strongly irreducible product of independent random matrices under optimal moment assumptions

A central limit theorem with explicit Lyapunov exponent and variance for products of $2\times2$ random non-invertible matrices

On the moments of the variance-gamma distribution

A recursion formula for mixed trace moments of isotropic Wishart matrices and the Gaussian unitary/orthogonal ensembles

The Gaussian product inequality conjecture for multinomial covariances

Moment formulae for general point processes

Central and noncentral moments of the multivariate hypergeometric distribution

Product formulas for multiple stochastic integrals associated with Lévy processes

Explicit Formulas for the Fourth Moment of Mixed Exponential Sums

Exact Distribution for the Product of Two Correlated Gaussian Random Variables.

On Wilks' joint moment formulas for embedded principal minors of Wishart random matrices