-1-

Non-Central Multivariate Chi-Square and Gamma Distributions
Thomas Royen
TH Bingen, University of Applied Sciences
e-mail: thomas.royen@t-online.de
A ( p - 1) -variate integral representation is given for the cumulative distribution function of the general p - variate
non-central gamma distribution with a non-centrality matrix of any admissible rank. The real part of products of
well known analytical functions is integrated over arguments from ( ,  ). To facilitate the computation, these
formulas are given more detailed for p  3. These ( p - 1) -variate integrals are also derived for the diagonal of a
non-central complex Wishart matrix. Furthermore, some alternative formulas are given for the cases with an
associated “one-factorial” ( p p)- correlation matrix R , i.e. R differs from a suitable diagonal matrix only by a
matrix of rank 1, which holds in particular for all (3  3) - R with no vanishing correlation.

1. Introduction and Notation
The following notations are used: ( n ) stands for  n1 ... n p  n with n1 ,..., np 

0

. The spectral norm of any

square matrix B is denoted by ||B|| and |B| is its determinant. For a symmetrical matrix, A  0 means positive
definiteness of A, and I or I p is always an identity matrix. The symbol  also denotes a non-singular covariance matrix ( ik ) with 1  ( ik ) , and R  (rik ) is a non-singular correlation matrix with R1  (r ik ). The
abbreviation Lt stands for “Laplace transform” or “Laplace transformation”. Formulas from the NIST Handbook of
Mathematical Functions are cited by HMF and their numbers.
The p - variate non-central chi-square distribution with ν degrees of freedom, an “associated” covariance matrix

 and a “non-centrality” matrix  (the  p2 ( , , ) - distribution) is defined as the joint distribution of the diagonal elements of a non-central Wp ( , , ) - Wishart matrix. To avoid confusion with different multivariate distributions having univariate (non-central) chi-square marginal distributions, this distribution can also be called a (noncentral) “Wishart chi-square distribution”. Without loss of generality (w.l.o.g.) we can standardize  to a correlation matrix R. The ( p ) - matrix Z of the corresponding Wishart matrix ZZ T has  independent N p ( j , R) normally distributed columns and   12 MM T with the ( p ) - matrix M  ( 1 ,...,  ). Therefore, rank()

min( p, ). (In the literature, frequently 1 is called “non-centrality” matrix.) This distribution can be applied
to power calculations for many statistical tests with several correlated chi-square test statistics, (see e.g. section
6 in [13]), but also in wireless communication at least the trivariate case occurs. For p  3 the corresponding cumulative distribution function (cdf) is not easy to compute. The complexity of the  p2 ( , R, ) - distribution (or of
__________________________________________________________________________________________
Key words and phrases: non-central multivariate chi-square distribution, non-central multivariate gamma distribution
2010 Mathematics Subject Classifications: 62H10, 62E15

-2-

the slightly more general  p ( , R, ) - distribution defined in (2.2) below) depends - apart from the dimension

p - essentially on two ranks, the rank k of  and the minimal possible rank m of R  V , where V is a suitable
non-singular - not necessarily positive definite (pos. def.) - diagonal matrix. For many applications, only distributions with rank()  1 are required, e.g. with identical expectation vectors  j in the above matrix M . These
cases can be reduced to central  p21 ( , ) - distributions by theorem 4.1 in [13] or by a similar theorem in [8].
For p  2 see formula (3.24) and (3.25) in section 3. However, a bivariate correlation matrix is always “one-factorial”, i.e. m  1. The  p2 ( , R, ) -cdf with m  1 and V  0 was already derived in [13] from a suitable representation of the Lt of the corresponding probability density (pdf). For the corresponding  p ( , R, ) - cdf see
[14]. These distributions are included in section 3 and extended to cases with an indefinite V . In particular, all
R33 with no vanishing correlation are “one-factorial” in this extended sense or limit cases with | V | 0. In [2] a
special 32 ( , R, ) - pdf is found in formula (49) with rank()  1 and r13  0. This is a limit case of the distributions with a “one-factorial” R33 with v 22  0. A simpler form of the  p ( , R, ) - cdf, p  3, of such limit
cases is also derived in section 3.
The most general  p ( , R, ) - cdf is represented by a p - variate integral over the cube
section 4 and by a ( p  1) - variate one over

p 1

p

 ( ,  ] p in

in section 5. These latter integrals are written in a more

detailed form in formula (5.7) and (5.8) for p  2 and p  3 to facilitate the direct computation. At least for the
above mentioned cases with m  1 and an indefinite V33 these bivariate integrals are supposed to be more
favourable for a numerical evaluation than the corresponding formulas in section 3. A further advantage is the
visibility of the integrand if p  3. The extension with a non-central complex Wishart matrix is found in section 6.
More general integrals over

p

for convolutions of non-central gamma distributions were already given in [16]. In

particular, for p  1 an integral over (0,  ) is found there for the cdf of a pos. def. quadratic form of (non-central) Gaussian random variables. For p  2 e.g. a bivariate integral is obtained for the cdf of Jensen’s bivariate
chi-square distribution (see also [5] and [17]). A survey on formulas for the central  p2 ( , R) - distribution and
some similar distributions is found in the appendix of [3].

2. Some Preliminaries and Special Functions
Definition. Any non-singular covariance matrix  p p is called “positive m-factorial” or only “m-factorial” if m is
the lowest possible rank of A in a representation   V  AAT with any pos. def. diagonal matrix V , where

A may contain a mixture of real or pure imaginary columns. If m is the lowest possible rank in such a representation with a non-singular not necessarily pos. def. diagonal matrix V , then  is called “real m-factorial”.
If  p p is m-factorial with V  W 2 then W W  I p  BBT , where B  WA has p rows b j and m columns

b and - w.l.o.g. - bT b  0,   . Each  p p 0 is at most ( p-1)- factorial with W 2   I p , where   0
is the lowest eigenvalue of .

-3-

The  p2 ( , R, ) - pdf has the Lt
fˆ (t1 ,..., t p ; , R, ) |I p  2 RT | /2 etr(2T ( I p  2 RT ) 1 ),

(2.1)

T  diag (t1 ,..., t p )  0,  ,  12 MM T ,rank()  min( p, ).
The p - variate non-central standardized gamma - pdf (in the extended sense of Krishnamoorthy and Parthasarathy [6]) - the  p ( , R, ) - pdf - can be defined by the Lt

gˆ (t1 ,..., t p ; , R, ) |I p  RT | etr(T ( I p  RT ) 1 ).

(2.2)

For positive integers   2 this distribution differs from the  p2 ( , R, ) - distribution only by a scale factor 2,
but this distribution is also defined for some non-integer values 2 . The function g ( x1 ,..., x p ; , R, ) with the
Lt from (2.2) is not for all   0 a pdf, but additionally to 2 1,..., p  2 , p  3,rank()  2 , at least all
real values 2  p  1, p  2, are admissible with any   0, see e.g. [7] or [9] (and also [11] if rank()  1 ).
Recently, in [20] the central  p ( , R) - distribution was also shown to exist for all values 2  [( p  1) / 2]
without any restrictions for R. For an m-factorial R with WRW  I p  BBT the  p ( , R) - cdf is given by

E



p
j 1



G (w2j x j , 12 b j S2 b jT ) ,2 

or 2  m  1  [( p  1) / 2],

(2.3)

with the univariate non-central gamma cdf G ( x, y) from (2.8) and the expectation referring to the Wm (2 , I m ) Wishart matrix S 2 , see [13] or [15]. If B is real, then in (2.3) all non-integer values 2  m  1  0 are admissible too. Hence, the  p ( , R, ) - distribution also exists for all real values 2  m with such m-factorial R
and rank()  1. Formula (2.3) was also used in [19] to extend the proof of the Gaussian correlation conjecture
in [18]. If |I  RT |1 is infinitely divisible, then |I  RT | is the Lt of a  p ( , R) - pdf for all   0 and then
formula (2.2) is the Lt of a  p ( , R, ) - pdf for all values 2  1 if rank()  1. Sufficient and necessary
conditions for R entailing infinite divisibility are found in [1] and [4]. Special infinitely divisible  p ( , R) - distributions arise by a one-factorial R  W 2  aaT with a real column a.
If R is real one-factorial with a real or pure imaginary column a and possibly with one wk2  0, we have the
relations
WRW  I p  bbT with wj 2  1  a 2j ,b j  w j a j , :| WRW | 1  bT b,
(WRW )1  I p   1bbT .

(2.4)

We need the following functions: The univariate gamma density

g ( x)  (( ))1 x 1e x , ,x  0,
with shape parameter  and the corresponding cdf

(2.5)

-4x

G ( x)   g ( )d  P( , x)   ( , x) / ( )

(2.6)

0

with the notations from HMF 8.2.4 for the incomplete gamma function  ( , x), which can be extended to

, cut

along the negative real axis. Later on, we shall use the fact that x   ( , x) is a single valued function on

.

The non-central gamma density g ( x, y ) with non-centrality parameter y is given by

yn



g ( x, y)  e y n0 g n ( x)

n!

 e x y ( xy )( 1)/2 I 1 (2 xy )  e y g ( x)0 F1 ( ; xy)

(2.7)

with the modified Bessel function I 1 and the hypergeometric function 0 F1 . The corresponding non-central cdf
is


G ( x, y)  e y n0 G n ( x)

yn
n!

.

(2.8)

Besides,


G* ( x, y) : e yG ( x, y)   k 1 g k ( x)0 F1 (  k ; xy)

(2.9)

will be used.
For n 

we have the relations

G1/2n ( x, y)  12 (erf ( x  y )  erf ( x  y ))  e y k 1 g1/2k ( x)0 F1 ( 12  k; xy),

(2.10)

G1n ( x, y)  G1 ( x, y)  e y k 1 g1k ( x)0 F1 (1  k; xy),

(2.11)

n

n

Gn ( x, y)  1  Qn ( 2 y , 2 x )

with the Marcum Q - function Qn of order n. Furthermore, with

0, z  0 


G0 ( z ) :  12 , z  0 


1, z  0 

(2.12)

we have the integral representation, (see section 2 in [13]),

Gn ( x, y)  e x y ( xy )n/2 1

 y cos( n ) xy cos(( n1) )

0

y 2 xy cos( ) x

The functions G ( x, y) can be extended to x, y 

e2

xy cos( )

d  G0 ( x  y),n 

0

, y  0. (2.13)

too with the principal value of x if   .

Furthermore, with x  , y  , we need the functions

11y (G ( x)  y1 e x ( y 1)G ( xy )), y  1,  0 



n
1 (G ( x)  y1 e x ( y 1)G ( xy )),  1,  ,
(
x
,
y
)

G
(
x
)
y

 1 y  1

 n 0   n

 1


 xg ( x)  (1    x)G ( x), y  1. 
verified by Lt.

(2.14)

-5-

3. Distributions with “One-Factorial“ Correlation Matrices and the Trivariate Distributions
Some earlier results are compiled here and extended by theorem 1. The  p ( , R, ) - cdf is given for all nonsingular “real one-factorial” correlation matrices R from (2.4) and any admissible rank of   0 and additionally
for the limit cases with wk2  1  ak2  0 for one index k . In particular, all non-singular R33  (rik ) with

r12 r13r23  0 are real one-factorial or limit cases of such matrices. We have the following subclasses: pure imaginary or real column a with max a 2j  1 , real column a with one value ak2  1 and the limit cases with one
value ak  1 (replacing a by a avoids ak  1). An essential simplification is obtained for rank()  1 . The
cases with max a 2j  1 were already given in [13] and [14]. The totally different unified ( p  1) - variate integral
representation in section 5 is suitable for feasible computations for p  3 and might be preferred for a numerical
evaluation at least for cases with ak2  1 or for the here not included case with one vanishing correlation.
The key for the formulas in theorem 1 for the  p ( , R, ) - cdf are suitable representations of the Lt from (2.2)
with R as in (2.4), which were already given by the formulas (2.16) and (2.17) in [13] for the corresponding
 p2 ( , R, ) - distribution with max a 2j  1. With
R  W 2  aaT  WRW  I p  bbT ,W 2  diag (w12 ,..., w2p ), b j  w j a j ,

(3.1)

(one wk2  0 admissible),

D  R1R1, z j  (1  wj 2t j )1 , u j  wj 2t j z j  1  z j , Z  diag ( z1 ,..., z p ),U  I p  Z ,

(3.2)

T  diag (t1 ,..., t p )  0,
we have the following identities
|I p  RT | etr(T ( I p  RT )1 )) 

 bTUW WUb 
|Z | (1  bTUb) etr(W WU ) exp 

T
1

b
Ub



(3.3)

 aT ZDZa 
exp(aT Da)|Z | (1  bTUb) etr(W 1DW 1U ) exp 
.
T
 1  b Ub 

(3.4)

With rank()  1,   T  ( i j ),D  dd T  (di d j ),
these formulas become





 (bTUW  )2 
|Z | (1  b Ub) exp  j 1 w u j exp 

T
 1  b Ub 




T

p

2
j

2
j

(3.5)

and





 (aT Zd )2 
exp((a d ) )|Z | (1  b Ub) exp  j 1 d w u j exp 
,
T
 1  b Ub 
T

2



T



p

2
j

2
j

(3.6)

respectively. With  | WRW | 1  bT b we can also write

1  bTUb   (1   1bT Zb).

(3.7)

-6-

If there is a negative value wk2  1  ak2 , we obtain a negative  and the above formulas will be used later with

tk  (ak2 1)1 entailing zk  0 and 1   1bT Zb  0.
The identity (3.3) follows by straightforward calculation from
1


T ( I p  RT )1  T (Z 1  aaT T )1  TZ ( I p  aaT TZ )1  TZ  I p 
aaT TZ  
T
1  a TZa


1


TZ  I p 
W 1bbTW 1TZ  .
T
1  b Ub



For the equality of the representations in (3.3) and (3.4) only the equation
1

1

a Da  tr(W DW U ) 
T

aT ZDZa

  tr(W WU ) 

1  bT Ub

bT UW WUb

(3.8)

1  bTUb

has to be verified. The left hand side can be written as
bTW 1 DW 1b  tr (W 1 DW 1U ) 

bT ( I  U )W 1 DW 1 ( I  U )b
1  bTUb

.

After having inserted

W 1DW 1  (WRW )1W W (WRW )1  (I   1bbT )W W (I   1bbT )
the equality (3.8) can be verified by elementary calculations.
For the limit case  : 1  ak2  0 the limit
 t y 
lim(1   t )( n) exp  
exp(ty)
 0
 1 t  

of the Lt of  1 g n ( 1 x,  1 y) implies
0,x  y 
lim G  n ( 1 x,  1 y )  
.
  0
1,x  y 

For aT Da 

(3.9)

dii , j  i 
dij* ai a j ,dij*  

1i  j  p
2dij , j  i 



we need
(aT Da )n

n!



d (n1 ,..., n p ) j 1 a j j ,d (n1 ,..., n p )  
p

(2 n )

n



*n

1 d ij ,
ij
n
ij !
1i  j  p

(3.10)

where the summation is taken over the nij with n j  i j nij  i j n ji , j  1,..., p.
Theorem 1. The  p ( , R, ) - cdf with the non-singular correlation matrix R  W 2  aaT
( WRW  I p  bbT ),W  diag (w1 ,..., wp ) with possibly one imaginary value wk , is given by

G( x1 ,..., x p ; , R, ) 
exp(aT Da)



0

 


n 0



d (n1 ,..., n p ) j 1 a j j G n j (w2j x j , d jj wj 2  b2j y) g n ( y)dy,
p

(2 n )

max a 2j  1,D  (dij )  R 1R 1 ,

n

(3.11)

-7-

or by

exp(aT Da)etr(DW 2 ) 




0

n 0

 

  n





d (n1 ,..., n p ) j 1 a j j G* n j (w2j x j , d jj wj 2   1b2j y) g n ( y)dy,
p

(2 n )

n

(3.12)

max a 2j  1ormax a 2j  1,  | WRW | 1  bT b, G* n ( x, y)  e yG n ( x, y) from (2.9).

With rank()  1,   T  ( i j ),D  dd T  (di d j ) these formulas are simplified to
G( x1 ,..., x p ; , R, ) 

exp((aT d ) 2 ) 



0

 

p

0

j 1



( 2 n)! 
n  0 n! 



 

(2 n )

(a j d j )
nj!
j 1

p

   G (w x ,  w  b y  2b  w y
1
f ( )   B( 12 ,   12 )  (sin  )2( 1) ,max a
0

2
j

j

2
j

2
j

2
j

1/2

j

j

j



or to



nj

2
j


G  n j (w2j x j , d 2j wj 2  b 2j y ) g  n ( y )dy 


cos( )) f ( ) g ( y)ddy,

(3.13)
(3.14)

 1,



exp((aT d )2 )exp  j 1 d 2j wj 2 




0

n 0

 



p

( 2n)! 
n! 

  n



 

(2 n )

(a j d j )
nj!
j 1

p



nj


G*  n j (w2j x j , d 2j wj 2   1b 2j y ) g  n ( y )dy 


(3.15)

exp((aT d )2 )exp  j 1 d 2j wj 2 
 

p

0

j 1

  


0

p

G* (w2j x j , d 2j wj 2   1b2j y  2 1/2b j d j wj 1 y1/2 cos( )) f ( ) g ( y)ddy,

(3.16)

max a 2j  1 ormax a 2j  1.

If exactly one value ak2  1, then

G( x1 ,..., x p ; , R, ) 


nj
2
2
2
  (2 n ) d (n1 ,..., n p ) a j G n j (wj x j , d jj w j  b j y) g n ( y)dy,
j k


and, if rank()  1, this is simplified to
exp(aT Da)

xk

0

exp((a d ) )
T

2



xk

0

 G (w x , 
j k

2
j

j



n 0


2
j

( 2n)! 
n 0 n! 





p

 
(2 n )

j 1

(a j d j )
nj!

nj


  G n j (w2j x j , d 2j wj 2  b2j y ) g  n ( y )dy 
j k


w2j  b2j y  2b j j w j y1/2 cos( )) f ( ) g ( y)dyd ,

(3.17)

(3.18)
(3.19)

 (, y) | [0,  ], y  0, k2  y  2 k y1/2 cos( )  xk .

Remarks. With y  ( k cos( )  ( xk   k2 sin 2  )1/2 )2 and 1  arcsin(| k |1 x1/2
k ) we obtain



(, y) |0  y  y ( ),0     ,0 k2  xk
( , y) | y ( )  y  y ( ),  1     , k2  xk , k  0
(, y) | y ()  y  y ( ),0    1, k2  xk , k  0

(3.20)

-8-

For p  3 the coefficients d (n1 ,..., n p ) are

d (n1 , n2 ) 

0 n12

d (n1 , n2 , n3 ) 

 
 min( n ,n ) i  j
1

2

[ n1 / 2 ] [ n2 / 2 ] [ n3 / 2 ]

  
n 0 n  n
i j
22 0 33  0

11

n

( dij* ) ij
nij ! ,ni

 n12 even,nii  12 (ni  n12 )

(3.21)

n

( dij* ) ij
nij ! ,nij

 12 (ni  n j  nk )  (n ii  n jj  nkk ),i  j ,k  i, j.

(3.22)

Proof of Theorem 1. All the representations for the cdf in theorem 1 can be verified by the Lt of the correspondding densities. The Lt gˆ  n (t , y ) of w2 g n (w2 x, y) is
(1  w2t )( n) exp(w2t (1  w2t )1 y)  z neuy ,u  1  z,

which also holds for w2  0 if t  w2  0 and consequently z  0. Then, from the formulas (3.11), (3.12) we
obtain the Lt in (3.3) and (3.4), from the formulas (3.13), (3.14) the Lt in (3.5), and from (3.15), (3.16) the Lt in
(3.6). For the verification of the formulas (3.14), (3.16) the integral representation of the modified Bessel functions

I 1 from HMF10.32.1 is used for the integration over  within the corresponding integral for the Lt. For the
verification of the cases with one wk2  1  ak2  0, the variable tk of the Lt is supposed to be larger than  wk2 ,
which is equivalent to zk  0. If for the moment the integrals for the corresponding asserted representations of
the pdf are denoted by h( x1 ,..., x p ), we obtain with a sufficiently large   0 a Lt of exp(  xi h( x1 ,..., x p )
coinciding with the Lt from (2.2) with T   I p ,T  0, instead of T . Therefore, h( x1 ,..., x p ) has the Lt. in (2.2),
which proves the corresponding representation. Finally, the formulas for the limit cases with ak2  1 follow from
(3.9).
The relationship between a  p2 ( , R, ) - pdf with   12 MM T  12 x p 1 T ,M p  ( y1 ,..., y  ),x p 1 




j 1

y 2j  and a central  p21 ( , ) - pdf with
 R   T






T



1

is described in theorem 4.1 in [13] (for a similar theorem see [8]). For the corresponding densities, this theorem
gives the relation
f ( x1 ,..., x p ; , R, )  ( 12 g /2 ( 12 x p 1 ))1 f ( x1 ,..., x p 1; , ).

(3.23)

Here, a slightly more general version is shown for the relation between a  p1 ( , ) - pdf and a  p ( , 0 , ) -

    T
pdf with    0 T
 



T
 and   x p 1 .
1

Theorem 2. The conditional density ( g ( x p 1 ))1 g ( x1 ,..., x p 1;  , ) is equal to g ( x1 ,..., x p ;  , 0 , x p 1 T ).
Proof. With Tp  diag (t1 ,..., t p ) the Lt of g ( x p 1 ) g ( x1 ,..., x p ; , 0 , x p 1 T ) is given by

-9

|I p   0Tp |



 exp   x

p 1

( T Tp ( I p   0Tp ) 1  t p 1  g ( x p 1 )dx p 1 

0



|I p  0Tp |

1  t

p 1

  T Tp ( I p  0Tp )1 



,

which coincides with the Lt |I p 1  T | of g ( x1 ,..., x p1; , ) since
|I p 1  T | (1  t p 1 )|I p  0Tp   T Tp  t p 1 (1  t p 1 )1 T Tp |=(1  t p 1 )|I p  0Tp  (1  t p 1 )1 T Tp |
(1  t p 1 )|I p  0Tp ||I p  (1  t p 1 )1 ( I p  0Tp )1 T Tp | and |I p  (1  t p 1 )1 ( I p  0Tp )1 T Tp |
1  (1  t p 1 )1 T Tp ( I p  0Tp )1 .

With the correlation r in a (2  2) - correlation matrix we also use the notation  2 ( , r ,  ). With the central case
from formula (3.11) and theorem 2 with
1 r 
0  

r 1

and a one-factorial correlation matrix

R  diag ((1  12 )1/2 ,(1   22 )1/2 ,1)diag ((1  12 )1/2 ,(1   22 )1/2 ,1)  W 2  aaT ,WRW  I3  bbT ,
we obtain the 2 ( , r, x3 T ) - cdf

G( x1 , x2 ; , r , x3 )  ( g ( x3 ))
T

1



 2

G ((1   j2 ) 1 w2j x j , b2j y) w32 g ( w32 x3 , b32 y) g ( y)dy.
0  
j 1


(3.24)

In a similar way we can apply formula (3.12) for the central case if max a 2j  1.
Alternatively, we can use the absolutely convergent series for a trivariate density in formula (3.5b) in [12]. With

R1  (r ij ) and Q  (qij ),qij  r ij (r ii r jj )1/2 , this leads to
G( x1 , x2 ;  , r , x3 T ) |Q|  g ( x3 ) 


 (1) N

N 2

1

(r

2

33

g (r 33 x3 ) G ((1   j2 ) 1 r jj x j ) 
j 1

2



m1  m2  m3  M  N

q( , m1 , m2 , m3 , Q)r 33 g  M m3 (r 33 x3 ) G  M m j ((1   j2 ) 1 r jj x j ) ,

)

(3.25)

j 1

M  [ N / 2], N  N  2M ,q( , m1 , m2 , m3 , Q) 

1
( )

min( m j )

22 m N (  M  m)

m0

(2m   N )! j 1 (m j  m)!



3

3

 qij2 mk  N ,
k 1

i  j,k  i, j.
A further alternative is the integral over (0,  ) in (5.7), which also holds for rank()  2. If theorem 2 is applied
with p  3, again the formulas (3.14) and (3.16) are obtained.

4. A p-Variate Integral Representation for the General Non-Central Distribution Function
If a p - variate pdf or cdf with a known characteristic function is not explicitly available, it can be represented by a

p - variate integral over

p

by means of the Fourier inversion formula under rather general assumptions. In this

- 10 -

section an alternative p - variate integral representation for the general  p ( , , ) - cdf (and the pdf) is given,
which requires only integration over ( ,  ) p . A further advantage is the always possible reduction of this integral to a ( p  1) - variate one as shown in the following section. This facilitates the computation at least for

p  3 with a visible integrand. The simple idea is the representation of the required probability by an integral over
a p - variate Fourier series (or a p - variate Laurent series).
The Lt
|I p  T | etr(T ( I p  T )1 )  etr(1 )|I p  T | etr(( I p  T ) 1  1 )

(4.1)

with any non-singular covariance matrix  and T  diag (t1 ,.., t p )  0 is not for all   0 and all non-centrality
matrices   0 the Lt of a pdf. The following integral representation also holds for real measures whose density
has the above Lt. For admissible values  and rank() of a  p ( , , ) - pdf see the explanations behind
(2.2) and (2.3). For a suitable form of the Lt we use the matrices
B  (v)1  I p ,D  v111

(4.2)

with any scale factor v, leading to ||B|| < 1, e.g.
1
1
v  12 (min
 max
)

(4.3)

with the extremal eigenvalues of , which leads to ||B||  (max  min )(max  min )1. With
Z  diag ( z1 ,..., z p ),z j  (1  v1t j )1 ,

(4.4)

we obtain
|I p  T | =|I p  vv-1T | | vΣ||(v
 )-1  I p  I p  v-1T ||v||B  Z -1| |v||Z -1|| I p  BZ |

and
tr(( I p  vv1T )1 1 )  tr(( B  Z 1 )1 v111 )  tr(Z ( I p  BZ )1 D).

Hence, the Lt from (4.1) can be written as

 p

|v| etr(1 )|Z | |I p  BZ | etr(Z ( I p  BZ )1 D)    zj  b( z1 ,..., z p )
 j 1 
with a p - variate analytical function

(4.5)

b( y1 ,..., y p )    (n1 ,..., np ) y1n1  ...  y p p ,
n

(4.6)

which is abs. convergent if max|y j |||B||1.
Theorem 3. With the functions



from (2.14) the  p ( , , ) - cdf is given by

G ( x1 ,..., x p ;  , , )  |v| etr(  1 )(2 )  p


p

etr(Y 1 ( I p  BY 1 ) 1 D)
|I p  BY 1 |

p




j 1



(vx j , y j ) d j ,
i j

p

 ( ,  ] p ,v from (4.3), the matrices B and D from (4.2), Y  diag ( y1 ,..., y p ), y j  re

,r ||B||.

- 11 -

Proof. The asserted integral representation of the theorem is equivalent to a corresponding representation for

the  p ( , , ) - pdf with the derivatives
 (vx j , y j ). Comparing this representation with the Lt from (4.5)
x j
we obtain

     (n ,..., n )

(2 ) p

1

p

n j

p

y
j 1 j

  

p



vg  n j (vx j ) y j j d1...d p 
n

j 1

p

  (n ,..., n )
1

p

p
j 1

vg  n j (vx j ) ,

which is abs. convergent because of lim   n max
n

of this series is

 z    (n ,..., n )
p



j 1

j

1

p

p
j 1



p
j 1



vg n j (vx j ) | j 1 n j  n  0 for every   0. The Lt
p

n

z j j , coinciding with the Lt in (4.5). Integration over x1 ,..., x p

concludes the proof.
Remarks. More generally, the scale factor v  w2 can be replaced by a scale matrix W 2  diag ( w12 ,..., w2p ) 
0. Then, with

TW  W 2T ,W  W W ,W  W W , z j  (1  wj 2t j ) 1 ,  diag ( z1 ,..., z p ),B  W1  I p , and

D  W1W W1  W 111W 1

(4.7)

the Lt in (4.1) becomes

|W W | etr(1 )|Z | |I p  BZ | etr(Z ( I p  BZ )1 D) 



p
j 1

zj

   (n ,..., n )
1

p

p
j 1

n

z j j , (4.8)

where the series is abs. convergent if max z j  ||B||1 ||W -1-1W -1  I p ||1 , i.e. for all sufficiently large values

t j if ||B||  1. This can occur in particular with a “natural” scaling if
Q  (qij )  W 11W 1

(4.9)

is a correlation matrix with 1  ( ij ) and w2j   jj .
Then, with any scale matrix W and the matrices B,D from (4.7), the  p ( , , ) - cdf is given by

G( x1 ,..., x p ; , , ) 
|W W | etr( )(2 )


1

p


p

i j

y j  re

etr(Y 1 ( I p  BY 1 ) 1 D)
1



|I p  BY |

p




j 1



( w2j x j , y j )d j ,

(4.10)

,r ||B||.

The proof is the same as for theorem 3, but by the series expansion and the termwise Lt the correct Lt is
obtained first only for sufficiently large values t j with max z j ||B||1 , if ||B|| 

- 12 -

5. The (p-1)-Variate Integral Representation
In the integral representation from (4.10) we can perform the integration over one argument  j , w.l.o.g. j  p.
With

 Bpp
B  (W W )1  I p   T
 bp

bp 
 Dpp
1 1
1
1
 ,D  W   W   T
bpp 
 dp
i j

Ypp  diag ( y1 ,..., y p 1 ), y j  re

dp 
,
d pp 

(5.1)

, j  ( ,  ],

(5.2)

y0  y0 ( y1 ,..., y p 1 )  bTp (Ypp  Bpp ) 1 bp  bpp ,

(5.3)

q  q( y1 ,..., y p 1 )  (bTp (Ypp  Bpp ) 1 , 1) D(bTp (Ypp  Bpp ) 1, 1)T ,

(5.4)

L ( y1 ,..., y p1 ) 
etr(Ypp1 ( I p 1  BppYpp1 )1 Dpp )|I p 1  BppYpp1|  etr((Ypp  Bpp )1 Dpp )|Ypp  Bpp | |Ypp | ,

and the functions G* from (2.9) and

(5.5)

from (2.14) we obtain



Theorem 4. With the above notations the  p ( , , ) - cdf is given by

G( x1 ,..., x p ; , , ) 
etr(  1 )


|W W |



1
(2 )

p 1



L ( y1 ,..., y p 1 )
p1



(1  y0 )



p 1

G* ((1  y0 )w2p x p , (1  y0 )1 q)
j 1

q



( w2j x j , y j )d j ,r  ||B||.



Remark. lim (1  y0 ) G* (1  y0 )w2p x p , 1 y  ((  1))1 w2p x p ) 0 F1   w2p qx p ).
y0 1

0

For the proof we use two lemmas:
Lemma 1. With the notations (5.1) to (5.4) the equation
q
etr((Y  B) 1 D)|Y  B|  etr((Ypp  B pp ) 1 Dpp )|Ypp  B pp | exp( y p  y )(y p  y 0 ) 
0

holds for r  |y j |>||B||.
Proof. For

 App
A : Y  B   T
 bp

bp 

y p  bpp 

we obtain with the Schur complement
1
|A|= |App | (y p  bpp  bTp App
bp ) |Ypp  Bpp |(y p  y0 ),

(5.6)

- 13 1
1
 A1  ( y p  y0 )1 App
bp bTp App
A1   pp
1

( y p  y0 )1 bTp App


1
( y p  y0 )1 App
bp 

1
( y p  y0 )


and
1
1
1
1
1
tr( A1D)  tr( App
Dpp  ( y p  y0 )1 ( App
bpbTp App
Dpp  App
bp d Tp ))  ( y p  y0 ) 1 (d pp  bpT App
dp) 
1
tr( App
Dpp )  (y p  y0 )1 q,

which implies equation (5.6).
Lemma 2. Let be q 

any number, Sr   y  ||y| r and y0 

\ {1} any number with |y0 |r |y|,

then

1
2 i

S exp(( y  y )
0

r

1

 ( x, y)( y  y0 )

q)



y 1dy  (1 y0 ) G* ((1  y0 ) x,(1  y0 )1 q).

Proof. With the binomial series we get

1
2 i


r x

S

y 1

 ( x, y) ( y y ) n
0

dy  21 i

S 
r


m 0

g m ( x) y m




k 0

  y

( nk ) y0
( n)k! y

k

 n1

dy 



1
g nk ( x)(  n  k ) y0k /k !  g n ( x) exp( xy0 )  (1  y0 )  n (1  y0 ) g n ((1  y0 ) x).
( n) 
k 0

Multiplication by qn /n! and summation over n yields
(1  y0 ) x G* ((1  y0 ) x,(1  y0 )1 q)).

Integration over x provides the assertion.
Proof of Theorem 4. If y p in Y is replaced by a variable y with any value |y|, then the equation
|Y  B||Ypp  Bpp |( y  bpp  bTp (Ypp  Bpp )1 bp )  0

always has a unique solution y  y0  bTp (Ypp  Bpp )1 bp  bpp , where |y0 | r since ||Bpp || ||B|| r  | y j |,

1  j  p  1. Hence, with lemma 1 and lemma 2, theorem 4 is obtained by integration over  p in the integral
from (4.10).
For p  2 we obtain with

 2
 1
  1 2

 1 2  2
 0
2
2
 ,w j   j (1   ),B  
2
2 
 

 
i
1 1
1
1
 ,D  (dij )  W   W , y  re ,r  ||,
0 

the cdf G( x1 , x2 ; , , )  (1   2 ) etr(1 )  `

1





 exp(d11 y 1 ) *  (1   2 y 1 ) x2 d11  2 y 1  2d12   d 22 y 
G  2
,

2 1 
2
y  2
 (1   y )
  2 (1   )




0




 
x1
, y  d.
 2
2
  1 (1   )  

(5.7)

- 14 -

For p  3 we find with 1  ( ij ),w2j   jj , the correlation matrix Q  (qij )  W 11W 1 ,B  Q  I 3 ,D 
(dij )  W 1 1 1W 1 , y0 ( y1 , y2 ) 
1
y1 y2  q122



q232 y1  q132 y2  2q12 q13 q23
y1 y2  q

2
12

q( y1 , y2 )
, y j  rei , j  1,2,r ||B||,

1  y0 ( y1 , y2 )

(q13 y2  q12 q23 ,q23 y1  q12 q13 ,q122  y1 y2 ) D(q13 y2  q12 q23 ,q23 y1  q12 q13 ,q122  y1 y2 )T
y1 y2  q232 y1  q132 y2  2q12 q13q23  q122

 d 22 y1  d11 y2  2d12 q12 
q122 
exp 
 1 

y1 y2  q122

  y1 y2 

, and L ( y1 , y2 ) 



the cdf

G( x1 , x2 , x3 ; , , )  |Q| etr(1 ) 
1

 



q( y1 , y2 )  2
*
33

G
(1

y
(
y
,
y
))

x
,



0
1
2
3

1  y0 ( y1 , y2 )  j 1

 (1  y0 ( y1 , y2 ))


 

2 2 0 

L ( y1 , y2 )





jj


x j , y j  d1d2 . (5.8)


The values (1   2 y 1 )  and (1  q122 ( y1 y2 )1 ) are defined by the binomial series and the functions 



( x, y )

and (1  y0 ) G* ((1  y0 ) x,(1  y0 )1 q), x  0, are single-valued functions. If, e.g., one element qik  0,

i  k , then ||Q - I 3 ||  1. If ||Q - I3 ||  1, then also theorem 4 can be applied with W 2  vI p ,v from (4.3), B 
v11  I3 with ||B||<min(r ,1), and y0 ,q , L from (5.3), (5.4), (5.5).

6. The Extension to the Distribution of the Diagonal of a Non-Central Complex Wishart Matrix
Let X be a ( p  ) - matrix with  independent CN p (  j , ) - distributed circular complex Gaussian column
vectors x j . The diagonal of XX * ,( X *  X T ), has a pdf with the Lt from (4.8) with the notations from (4.7) and

 instead of  , but now with a pos. def. Hermitian covariance matrix  and a pos. semi-def. Hermitian “noncentrality” matrix   MM * M  (1 ,...,  ). For the corresponding cdf F ( x1 ,..., x p ; , , ) we obtain again
the integral representation from (4.10), but with  replaced by  . The ( p  1) - variate integral representation of
theorem 4 can also be retained for F ( x1 ,..., x p ; , , ) with the same proof if we replace within the notations
(5.1), (5.3), (5.4) the row vectors bTp and d Tp by b*p and d *p respectively, and q by
1
1
1
1
q( y1 ,..., y p 1 )  b*p App
Dpp App
bp  d *p App
bp  b*p App
d p  d pp

with the non-Hermitian ( p  1)  ( p  1) - matrix App  Ypp  Bpp .

- 15 -

References
[1] R. B. Bapat, Infinite divisibility of multivariate gamma distributions and M-matrices, Sankhyā 51 (1989),
73-78.
[2] P. Dharmawansa and M. R. McKay, Diagonal distributions of a complex Wishart matrix: A new trivariate
non-central chi-squared density, J. Multivariate Anal. 100 (2009), 561-580.
[3] T. Dickhaus and T. Royen, A survey on multivariate chi-square distributions and their application in testing
multiple hypotheses, Statistics 49 (2015), 427- 454.
[4] R. C. Griffiths, Characterization of infinitely divisible multivariate gamma distributions, J. Multivariate Anal. 15
(1984), 13-20.
[5] D. R. Jensen, A generalization of the multivariate Rayleigh distribution, Sankhyā Ser A 32 (1970), 193-208.
[6] A. S. Krishnamoorthy and M. Parthasarathy, A multivariate gamma type distribution,. Ann. Math. Stat. 22
(1951), 549-557.
[7] E. Mayerhofer, On the existence of non-central Wishart distributions, J. Multivariate Anal. 114 (2013), 448456.
[8] K. S. Miller and H. Sackrowitz, Relationships between biased and unbiased Rayleigh distributions, SIAM J.
Appl. Math. 15 (1967), 1490-1495.
[9] R. J. Muirhead, Aspects of Multivariate Statistical Theory, 2nd ed. (2005), John Wiley & Sons, New York.
[10] F. W. J. Olver, D. W. Lozier, R. F. Boisvert and C. W. Clark, NIST Handbook of Mathematical Functions,
Cambrigde University Press (2010), Cambridge, New York.
[11] D. St. P. Richards and S. D. Peddada, Proof of a conjecture of M. L. Eaton on the characteristic function of
the Wishart distribution, Ann. Probab. 19 (1991), 868-874.
[12] T. Royen, Expansions for the multivariate chi-square distribution, J. Multivariate Anal. 38 (1991), 213-232.
[13] T. Royen, On some central and non-central multivariate chi-square distributions, Statistica Sinica 5 (1995),
373-397.
[14] T. Royen, Multivariate gamma distributions (Update), in Encyclopedia of Statistical Sciences, Update
Volume 1 (1997), 419-425, John Wiley and Sons, Inc., New York.
[15] T. Royen, Integral representations and approximations for multivariate gamma distributions, Ann. Inst.
Statist. Math. 59 (2007), 499-513.
[16] T. Royen, Integral representations for convolutions of non-central multivariate gamma distributions, (2007),
arXiv:0704.0539 [].
[17] T. Royen, Some representations for convolutions of multivariate gamma distributions, Far East J. Theor.
Stat. 45 (2013), 21-49.
[18] T. Royen, A simple proof of the Gaussian correlation conjecture extended to some multivariate gamma distributions, Far East J. Theor. Stat. 48 (2014), 139-145.
[19] T. Royen, Some probability inequalities for multivariate gamma and normal distributions, Far East J. Theor.
Stat. 51 (2015), 17-36.
[20] T. Royen, A note on the existence of the multivariate gamma distribution, (2016), arXiv:1606.04747[math.PR]

