arXiv:1705.00163v1 [math.PR] 29 Apr 2017

A Proof of the Explicit Formula for
Product Moments of
Multivariate Gaussian Random Variables
Iickho Song∗
School of Electrical Engineering
Korea Advanced Institute of Science and Technology
291 Daehag Ro, Yuseong Gu, Daejeon 34141 Korea

Abstract
A detailed proof of a recent result on explicit formulae for the product moments
E {X1a1 X2a2 · · · Xnan } of multivariate Gaussian random variables is provided in
this note.
Keywords: jointly Gaussian random variables, Gaussian random vector, joint
moment, multivariate Gaussian random variables, product moment

1. Introduction
Recently, explicit formulae for the product moments of multivariate Gaussian
random variables have been derived in Song and Lee (2015). The formulae not
only are simpler and more concise than any of the other ones but also allow us
5

to promptly compute any term of the product moments.
Since the publication of the results, a huge number of requests has been
received for the details of the proof. In this note, we provide a detailed proof of
the main result of Song and Lee (2015).
∗ Corresponding author. Tel.: +82-42-350-3445. This study was supported by the National Research Foundation of Korea, under Grant NRF-2015R1A2A1A01005868, for which
the author wishes to express his thanks.
Email address: i.song@ieee.org (Iickho Song)

Preprint submitted to https://arxiv.org/

May 2, 2017

2. Proof of the Explicit Formula
10

2.1. The Explicit Formula
n

Let {Xi }i=1 be multivariate Gaussian random variables with mean vector
µ = (µ1 , µ2 , · · · , µn ) and n × n variance-covariance matrix Φ = [ϕij ], where
µi = E {Xi } denotes the expected value of Xi and ϕij = E {Xi Xj } − µi µj
denotes the covariance between Xi and Xj if i 6= j and variance of Xi if i = j.
15

One of the main results of Song and Lee (2015) is as follows:
Theorem 1. For the product moments



)
(n
n
n Y
n
Y
Y
X
Y
L
l
µj a,j 
ϕijij  
da,l 
Xiai =
E
i=1

(1)

j=1

i=1 j=i

l∈Sa

of multivariate Gaussian random variables, we have
n
Q

k=1

da,l =

n Q
n
Q

2Ml

lij !

i=1 j=i

for l ∈ Sa , where Ml =

n
P

ak !
!

n
Q

La,j !

j=1

!

(2)

lkk . 

k=1

20


Here, the elements of the two sets a = {a1 , a2 , · · · , an } and l = {lij }ni=1

n
j=i

are nonnegative integers,

La,k = ak − lkk −

n
X

ljk

(3)

j=1

with lji = lij , and Sa denotes the collection of all the sets l such that
{La,k ≥ 0}nk=1 .

(4)

2.2. Proof
Let us prove Theorem 1 by mathematical induction. First, it is clear from
Lemma 2.2, Theorem 1, and Theorem 2 of Song and Lee (2015) that (2) holds
25

for n = 1, 2, and 3, respectively.

2

Next, assume that (2) is true when n = m − 1. Then, since

am−1
am
am
E {X1a1 X2a2 · · · Xm
} = E X1a1 X2a2 · · · Xm−1
}
E {Xm

(5)

when ϕ1m = ϕ2m = · · · = ϕm−1,m = 0, we get [da2 ,l2 ]l1m 7→0,l2m 7→0,··· ,lm−1,m 7→0 =
da1 ,l1 dam ,lmm , or
[da2 ,l2 ]l1m 7→0,l2m 7→0,··· ,lm−1,m 7→0 =
m
Q

ak !

k=1
m−1
Q m−1
Q

2Ml2

!

lij ! lmm !

i=1 j=i

(6)

!

La1 ,j ! (am − 2lmm )!

j=1

using (2) with n = m − 1 and

da1 ,l =

m−1
Q

a1 !
2l l! (a1 − 2l)!

(7)
m−1

30

(Eq. (5), Song and Lee (2015)) for l = 0, 1, · · · , ⌊ a21 ⌋, where a1 = {ai }i=1 ,
a2 = a1 ∪ {am }, l1 = {lij } with j = i, i + 1, · · · , m − 1 for i = 1, 2, · · · , m − 1,
m

and l2 = l1 ∪{lim }i=1 . Here, the symbol 7→ denotes a substitution: For example,
α 7→ β means the substitution of α with β.
Next, employing the result
∂
E
∂ϕij

35

(

n
Y

Xkak

k=1

)

= ai aj E









a −1

Xiai −1 Xj j

n
Y

Xkak

k=1
k6=i,j






(8)





for (i, j) = (1, m), (2, m), · · · , (m − 1, m) in Price (1958), and then taking some
steps (similar to those from (16) to (20) in Song and Lee (2015)), we will get
[da2 ,l2 ]lim 7→lim +1 =

ai !am ! [da2 ,l2 ]lim 7→0,ai 7→ai −lim −1,am 7→am −lim −1
(ai − lim − 1)! (am − lim − 1)! (lim + 1)!

(9)

and then, by substituting lim + 1 with lim in (9), we have
da2 ,l2 =

ai !am ! [da2 ,l2 ]lim 7→0,ai 7→ai −lim ,am 7→am −lim

(10)

(ai − lim )! (am − lim )! (lim )!

for i = 1, 2, · · · , m − 1. Now, letting l2m = 0 in (10) with i = 1, we get
[da2 ,l2 ]l2m 7→0 =

a1 !am ! [da2 ,l2 ]l1m 7→0,l2m 7→0,a1 7→a1 −l1m ,am 7→am −l1m
(a1 − l1m )! (am − l1m )!l1m !
3

,

(11)

which can be used into (10) with i = 2 to produce
da2 ,l2 =

40

a2 !am !
a1 ! (am − l2m )!
×
(a2 − l2m )! (am − l2m )!l2m ! (a1 − l1m )! (am − l1m − l2m )!l1m !
× [da2 ,l2 ]l1m 7→0,l2m 7→0,a1 7→a1 −l1m ,a2 7→a2 −l2m ,am 7→am −l1m −l2m .
(12)

Subsequently, letting l3m = 0 in (12), we get
[da2 ,l2 ]l3m 7→0 =

a1 !a2 !am !
l1m !l2m ! (a1 − l1m )! (a2 − l2m )! (am − l1m − l2m )!

× [da2 ,l2 ]lkm 7→0 for k=1,2,3, a1 7→a1 −l1m ,a2 7→a2 −l2m ,am 7→am −l1m −l2m , (13)
which can be employed into (10) with i = 3 to produce
da2 ,l2

=

a3 !am !
(a3 − l3m )! (am − l3m )!l3m !
a1 !a2 ! (am − l3m )!
×
l1m !l2m ! (a1 − l1m )! (a2 − l2m )! (am − l1m − l2m − l3m )!
× [da2 ,l2 ]lkm 7→0,ak 7→ak −lkm for k=1,2,3, am 7→am −l1m −l2m −l3m . (14)

From (14), it can be formulated

da2 ,l2

=



i
Q

ak ! am !


i
i
P
Q
lkm !
(ak − lkm )!
am −
k=1



i
Q

lkm !

k=1



k=1

k=1

× [da2 ,l2 ]



lkm 7→0,ak 7→ak −lkm for k=1,2,3, am 7→am −

(15)

i
P

lkm

k=1

with i = 3. Next, letting l4m = 0 in (15), we will get [da2 ,l2 ]l4m 7→0 which can be
employed into (10) with i = 4 to produce (15) with i = 4. After proceeding the
45

above procedures until i = m − 1, we will eventually get
m−1 
Q
ak ! am !
k=1

 m−1

da2 ,l2 = m−1
m−1
P
Q
Q
lkm !
lkm !
(ak − lkm )!
am −
k=1

× [da2 ,l2 ]

k=1

k=1

lkm 7→0,ak 7→ak −lkm for k=1,2,··· ,m−1, am 7→am −

m−1
P

(16)
lkm

k=1

from [da2 ,l2 ]lm−1,m 7→0 (obtained by letting lm−1,m = 0 in (15) with i = m − 2)

4

and (10) with i = m − 1. Finally, by combining (6) into (16), we can get
m−1 
Q
ak ! am !
k=1

 m−1

da2 ,l2 = m−1
m−1
P
Q
Q
lkm !
lkm !
(ak − lkm )!
am −
k=1
k=1
k=1

m−1

m−1
P
Q
lkm !
(ak − lkm )!
am −
k=1
k=1
!
×
m−1
Q m−1
Q
Ml2
lkj ! lmm !
2
k=1 j=i

× m−1 
Q
j=1

1

[La1 ,j ]ak 7→ak −lkm

for k=1,2,··· ,m−1

1

×
m−1
P
lkm − 2lmm !
am −



!

k=1

m
Q

k=1

=
2

Ml2

m Q
m
Q

i=1 j=i

ak !
!

lij !

m
Q

La2 ,j !

j=1

!,

(17)

which implies that (2) holds when n = m also.

3. Conclusion
50

We have provided a detailed proof of an elegant explicit formulae for product moments E {X1a1 X2a2 · · · Xnan } of multivariate Gaussian variables derived
recently.

References
Price, R., 1958. A useful theorem for nonlinear devices having Gaussian inputs.
55

IRE Trans. Inform. Theory 4, 69–72.
Song. I., Lee, S., 2015, Explicit formulae for product moments of multivariate
Gaussian random variables, Stat., Prob. Lett. 100, 27-34.

5

