Next Article in Journal
Renewable k-Out-of-n System with the Component-Wise Strategy of Preventive System Maintenance
Previous Article in Journal
Adaptive Trajectory Tracking Algorithm for the Aerospace Vehicle Based on Improved T-MPSP
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Expected Values of Scalar-Valued Functions of a Complex Wishart Matrix

by
Daya K. Nagar
1,*,†,
Alejandro Roldán-Correa
1,† and
Saralees Nadarajah
2,†
1
Instituto de Matemáticas, Universidad de Antioquia, Calle 67, No. 53-108, Medellín 050010, Colombia
2
Department of Mathematics, University of Manchester, Manchester M13 9PL, UK
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2023, 11(9), 2162; https://doi.org/10.3390/math11092162
Submission received: 28 February 2023 / Revised: 27 April 2023 / Accepted: 1 May 2023 / Published: 4 May 2023

Abstract

:
The complex Wishart distribution has ample applications in science and engineering. In this paper, we give explicit expressions for E ( tr ( W h ) ) g ( tr ( W j ) ) i and E ( tr ( W h ) ) g ( tr ( W j ) ) i , respectively, for particular values of g, h, i, j, g + h + i + j 5 , where W follows a complex Wishart distribution. For specific values of g, h, i, j, we first write ( tr ( W h ) ) g ( tr ( W j ) ) i and ( tr ( W h ) ) g ( tr ( W j ) ) i in terms of zonal polynomials and then by using results on integration evaluate resulting expressions. Several expected values of matrix-valued functions of a complex Wishart matrix have also been derived.

1. Introduction

The complex Wishart distribution has applications in various fields, including statistics, finance, physics, and engineering. Properties such as marginal and conditional distributions of sub-matrices, distribution of diagonal elements, and distribution of the determinant of the complex Wishart matrix W, including, of course, its moments, are therefore of broad interest. The motivation of the present paper is therefore to derive expressions for expected values of functions of W to further enrich the already existing literature on the complex Wishart distribution. For some early work on the complex Wishart distribution the reader is referred to Goodman [1], James [2], Khatri [3,4], Krishnaiah [5], Srivastava [6], Shaman [7], and Tan [8]. Some related material is considered in Alfano et al. [9], Carmeli [10], Cunden, Dahlqvist, and O’Connell [11], Deng et al. [12], Ermolova and Tirkkonen [13], Gomez et al. [14], Kumar [15], Nielsen, Skriver, and Conradsen [16], Shakil and Ahsanullah [17], Tague and Caldwell [18], and Tralli and Conti [19]. Systematic treatment of the complex Wishart distribution is available in Andersen et al. [20].
If W C W m ( n , Σ ) , n > m 1 , Σ H = Σ > 0 , then its p.d.f. is
etr Σ 1 W det ( W ) n m Γ ˜ m n det ( Σ ) n , W H = W > 0 ,
where det ( A ) = determinant of A, A H denotes the conjugate transpose of A, A = A H > 0 means that A is Hermitian positive definite, and the complex multivariate gamma function Γ ˜ m ( a ) is defined by
Γ ˜ m ( a ) = π m ( m 1 ) / 2 i = 1 m Γ a i + 1 , R e ( a ) > m 1 .
The complex Wishart distribution can be derived as the joint distribution of sample variances and covariances from a complex multivariate normal population (Goodman [1]). The parameter n need not be an integer, but, when n is not an integer, W can no longer be interpreted as a matrix of sample variances and covariances.
The complex Wishart distribution frequently arises in multivariate analysis as distributions of complex random matrices (for example, see Gupta and Nagar [21]) and hence plays a pivotal role in various branches of science and engineering. For properties and different variations of the Wishart distribution the reader is referred to James [2], Nagar, Gupta and Sánchez [22], Latec and Massam [23], Nagar, Roldán-Correa, and Gupta [24], Di Nardo [25], Dharmawansa and McKay [26], Hillier and Can [27], and Tralli and Conti [19].
Several authors have derived expected values of functions of a complex Wishart matrix. In her doctoral thesis, Grace Wahba [28] gave expressions for E W 1 . Shaman [7] derived several expected values including Cov vec W 1 . For any constant square matrix A of order m, Tague and Caldwell [18] derived E ( W A W ) , E ( W 1 A W 1 ) , and many other useful results. Sultan and Tracy [29] gave expressions for E ( W W ) , E ( W W W ) , and E ( W W W W ) . Maiwald and Kraus [30] gave approximations for E tr W 1 H 1 tr W 1 H 2 , E tr W 1 H 1 W 1 H 2 , E tr W 1 H 1 W 1 H 2 tr W 1 H 3 , and E tr W 1 H 1 W 1 H 2 tr W 1 H 3 W 1 H 4 . Hélène Massam and her co-workers (see Graczyk, Letac, and Massam [31]) gave a number of results which include E tr W H 1 tr W H 2 , E tr W H 1 W H 2 , E tr W 1 H 1 tr W 1 H 2 , and E tr W 1 H 1 W 1 H 2 . Several results on expected values of functions of a complex Wishart matrix including E ( W ) , E W 1 , E W Σ 1 W , E W 1 Σ W 1 , E W Σ 1 W Σ 1 W , and E W 1 Σ W 1 Σ W 1 are available in Nagar and Gupta [32]. A recurrence relation for E tr ( W k ) , where k is a non-negative integer, and several special cases of E tr ( W k ) , W C W m n , I m are available in Pielaszkiewicz, von Rosen, and Singull [33]. A collection of moments of the Wishart distribution is given in Holgersson and Pielaszkiewicz [34].
In this article, we compute expected values of functions of complex Wishart and inverted complex Wishart matrices. By definition
E ( tr ( W h ) ) g ( tr ( W j ) ) i = W H = W > 0 etr ( Σ 1 W ) det ( W ) n m det ( Σ ) n Γ ˜ m ( n ) ( tr ( W h ) ) g ( tr ( W j ) ) i d W ,
and
E ( tr ( W h ) ) g ( tr ( W j ) ) i = W H = W > 0 etr Σ 1 W det ( W ) n m det ( Σ ) n Γ ˜ m ( n ) ( tr ( W h ) ) g ( tr ( W j ) ) i d W ,
where g, h, i, and j are non-negative integers, and W C W m ( n , Σ ) .
In Section 3 and Section 4, we give explicit expressions for E ( tr ( W h ) ) g ( tr ( W j ) ) i and E ( tr ( W h ) ) g ( tr ( W j ) ) i , respectively, for specific values of g , h , i , and j. More precisely, we give expressions for E tr ( W 2 ) , E tr ( W 3 ) , E tr ( W 2 ) tr ( W ) , E [ ( tr ( W 2 ) ) 2 ] , E [ tr ( W 2 ) ( tr ( W ) ) 2 ] , E tr ( W 3 ) tr ( W ) , E tr ( W 4 ) ) 2 , E tr ( W 4 ) , E tr ( W 2 ) ( tr ( W ) ) 3 , E ( tr ( W 2 ) ) 2 tr ( W ) , E tr ( W 2 ) ( tr ( W ) ) 3 , E tr ( W 3 ) ( tr ( W ) ) 2 , E tr ( W 3 ) tr ( W 2 ) , E tr ( W 4 ) tr ( W ) , E tr ( W 5 ) , E tr ( W 2 ) , E tr ( W 3 ) , E tr ( W 2 ) tr ( W 1 ) , E ( tr ( W 2 ) ) 2 , E tr ( W 2 ) ( tr ( W 1 ) ) 2 , E tr ( W 3 ) tr ( W 1 ) , E tr ( W 4 ) , E tr ( W 2 ) ( tr ( W 1 ) ) 3 , E ( tr ( W 2 ) ) 2 tr ( W 1 ) , E tr ( W 3 ) ( tr ( W 1 ) ) 2 , E tr ( W 3 ) tr ( W 2 ) , E tr ( W 4 ) tr ( W 1 ) , and E tr ( W 5 ) . Several special cases of these results for Σ = I m are given in Section 5. Further, these special cases have been used to evaluate expected values of matrix-valued functions of a complex Wishart matrix and inverted complex Wishart matrices. Finally, Section 6 closes the paper with concluding remarks.

2. Some Known Definitions and Results

The Pochhammer symbol, denoted by ( a ) n , is defined as ( a ) n = a ( a + 1 ) ( a + n 1 ) = ( a ) n 1 ( a + n 1 ) for n = 1 , 2 , , and ( a ) 0 = 1 . The ordered partition κ of k is defined by κ = k 1 , , k m , k 1 + + k m = k , k 1 k m 0 . The complex generalized hypergeometric coefficient [ a ] κ , for an ordered partition κ of k, is given as
[ a ] κ = i = 1 ( κ ) a i + 1 k i ,
where ( κ ) is the number of non-zero k i s. Using (5), the computation of [ a ] κ can be carried out for ordered partitions of k. These coefficients are in Table 1 for k 5 .
For an ordered partition κ of k, we define Γ ˜ m ( a , κ ) and Γ ˜ m ( a , κ ) as
Γ ˜ m ( a , κ ) = [ a ] κ Γ ˜ m ( a ) , Γ ˜ m ( a , 0 ) = Γ ˜ m ( a )
and
Γ ˜ m ( a , κ ) = ( 1 ) r Γ ˜ m ( a ) [ a + m ] κ , Re ( a ) > k 1 + m 1 ,
respectively. The coefficients [ a + m ] κ for k 5 are in Table 2. We denote by C ˜ κ ( X ) the zonal polynomial (James [2]) of an m × m complex symmetric matrix X corresponding to the ordered partition κ . For small values of k, explicit formulas for C ˜ κ ( X ) derived in Khatri [35] are listed in Appendix A.
Next, we give two integrals involving zonal polynomials. These results will be used to evaluate integrals in Section 3 and Section 4.
Lemma 1.
Let Z be a Hermitian positive definite matrix of order m and let T be an m × m Hermitian matrix. Then, for R e ( a ) > m 1 , we have
R H = R > 0 etr ( Z R ) det ( R ) a m C ˜ κ ( T R ) d R = Γ ˜ m ( a , κ ) det ( Z ) a C ˜ κ ( T Z 1 ) .
Lemma 2.
Let Z be a Hermitian positive definite matrix of order m and let T be an m × m Hermitian matrix. Then, for R e ( a ) > k 1 + m 1 , we have
R H = R > 0 etr ( Z R ) det ( R ) a m C ˜ κ ( T R 1 ) d R = Γ ˜ m ( a , κ ) det ( Z ) a C ˜ κ ( T Z ) = ( 1 ) k Γ ˜ m ( a ) [ a + m ] κ det ( Z ) a C ˜ κ ( T Z ) .
Lemmas 1 and 2 are in James [2] and Khatri [4].

3. Expressions for E tr W h g tr W j i

Consider the expected value of tr W 2 as
E tr ( W 2 ) = W H = W > 0 etr Σ 1 W det ( W ) n m det ( Σ ) n Γ ˜ m ( n ) ( tr ( W 2 ) ) d W .
Writing tr W 2 in terms of zonal polynomials (see Appendix A) and integrating the resulting expression by using (6), we obtain
E tr ( W 2 ) = W H = W > 0 etr Σ 1 W det ( W ) n m det ( Σ ) n Γ ˜ m ( n ) C ˜ ( 2 ) ( W ) d W W H = W > 0 etr Σ 1 W det ( W ) n m det ( Σ ) n Γ ˜ m ( n ) C ˜ 1 2 ( W ) d W = [ n ] ( 2 ) C ˜ ( 2 ) ( Σ ) [ n ] 1 2 C ˜ 1 2 ( Σ ) .
Now, substituting for [ n ] ( 2 ) , [ n ] 1 2 , C ˜ ( 2 ) ( Σ ) and C ˜ 1 2 ( Σ ) from Table 1 and Appendix A in the above expression, we have
E tr ( W 2 ) = n tr ( Σ ) 2 + n 2 tr ( Σ 2 ) .
Writing tr W 3 in terms of zonal polynomials by using results given in Appendix A and integrating the resulting expression by applying (6), we obtain
E tr ( W 3 ) = W H = W > 0 etr Σ 1 W det ( W ) n m det ( Σ ) n Γ ˜ m ( n ) ( tr ( W 3 ) ) d W = [ n ] ( 3 ) C ˜ ( 3 ) ( Σ ) 1 2 [ n ] ( 2 , 1 ) C ˜ ( 2 , 1 ) ( Σ ) + [ n ] 1 3 C ˜ 1 3 ( Σ ) .
Now, substituting for [ n ] ( 3 ) , [ n ] ( 2 , 1 ) [ n ] 1 3 , C ˜ ( 3 ) ( Σ ) , C ˜ ( 2 , 1 ) ( Σ ) and C ˜ 1 3 ( Σ ) above, we have
E tr ( W 3 ) = n tr ( Σ ) 3 + 3 n 2 tr ( Σ ) tr ( Σ 2 ) + n ( n 2 + 1 ) tr ( Σ 3 ) .
Following the procedure described above, we evaluate E tr W h g tr W j i explicitly for selected values of g, h, i, j given by
E tr ( W 2 ) tr ( W ) = n 2 tr ( Σ ) 3 + n ( n 2 + 2 ) tr ( Σ ) tr ( Σ 2 ) + 2 n 2 tr ( Σ 3 ) ,
E ( tr ( W 2 ) ) 2 = n 2 tr ( Σ ) 4 + 2 n ( n 2 + 2 ) tr ( Σ ) 2 tr ( Σ 2 ) + n 2 ( n 2 + 2 ) ( tr ( Σ 2 ) ) 2 + 8 n 2 tr ( Σ ) tr ( Σ 3 ) + 2 n ( 2 n 2 + 1 ) tr ( Σ 4 ) ,
E tr ( W 2 ) ( tr ( W ) ) 2 = n 3 tr ( Σ ) 4 + n 2 ( n 2 + 5 ) tr ( Σ ) 2 tr ( Σ 2 ) + n ( n 2 + 2 ) ( tr ( Σ 2 ) ) 2 + 4 n ( n 2 + 1 ) tr ( Σ ) tr ( Σ 3 ) + 6 n 2 tr ( Σ 4 ) ,
E tr ( W 3 ) tr ( W ) = n 2 tr ( Σ ) 4 + 3 n ( n 2 + 1 ) ( tr ( Σ ) ) 2 tr ( Σ 2 ) + 3 n 2 ( tr ( Σ 2 ) ) 2 + n 2 ( n 2 + 7 ) tr ( Σ ) tr ( Σ 3 ) + 3 n ( n 2 + 1 ) tr ( Σ 4 ) ,
E tr ( W 4 ) = n tr ( Σ ) 4 + 6 n 2 tr ( Σ ) 2 tr ( Σ 2 ) + n ( 2 n 2 + 1 ) ( tr ( Σ 2 ) ) 2 + 4 n ( n 2 + 1 ) tr ( Σ ) tr ( Σ 3 ) + n 2 ( n 2 + 5 ) tr ( Σ 4 ) ,
E tr ( W ) 3 tr ( W 2 ) = n 4 tr ( Σ ) 5 + n 3 ( 9 + n 2 ) tr ( Σ ) 3 tr ( Σ 2 ) + 3 n 2 ( 4 + n 2 ) tr ( Σ ) ( tr ( Σ 2 ) ) 2 + 2 n 2 ( 7 + 3 n 2 ) ( tr ( Σ ) ) 2 tr ( Σ 3 ) + 4 n ( 3 + 2 n 2 ) tr ( Σ 2 ) tr ( Σ 3 ) + 6 n ( 2 + 3 n 2 ) tr ( Σ ) tr ( Σ 4 ) + 24 n 2 tr ( Σ 5 ) ,
E ( tr ( W 2 ) ) 2 tr ( W ) = n 3 tr ( Σ ) 5 + 2 n 2 ( 4 + n 2 ) ( tr ( Σ ) ) 3 tr ( Σ 2 ) + n ( 8 + 6 n 2 + n 4 ) tr ( Σ ) ( tr ( Σ 2 ) ) 2 + 4 n ( 2 + 3 n 2 ) ( tr ( Σ ) ) 2 tr ( Σ 3 ) + 4 n 2 ( 4 + n 2 ) tr ( Σ 2 ) tr ( Σ 3 ) + 2 n 2 ( 13 + 2 n 2 ) tr ( Σ ) tr ( Σ 4 ) + 8 ( n + 2 n 3 ) tr ( Σ 5 ) ,
E tr ( W 3 ) ( tr ( W ) ) 2 = n 3 tr ( Σ ) 5 + n 2 ( 7 + 3 n 2 ) ( tr ( Σ ) ) 3 tr ( Σ 2 ) + 3 n ( 2 + 3 n 2 ) tr ( Σ ) ( tr ( Σ 2 ) ) 2 + n ( 6 + 13 n 2 + n 4 ) ( tr ( Σ ) ) 2 tr ( Σ 3 ) + n 2 ( 19 + n 2 ) tr ( Σ 2 ) tr ( Σ 3 ) + 6 n 2 ( 4 + n 2 ) tr ( Σ ) tr ( Σ 4 ) + 12 n ( 1 + n 2 ) tr ( Σ 5 ) ,
E tr ( W 3 ) tr ( W 2 ) = n 2 tr ( Σ ) 5 + 2 n ( 3 + 2 n 2 ) tr ( Σ ) 3 tr ( Σ 2 ) + 3 n 2 ( 4 + n 2 ) tr ( Σ ) ( tr ( Σ 2 ) ) 2 + n 2 ( 19 + n 2 ) ( tr ( Σ ) ) 2 tr ( Σ 3 ) + n ( 6 + 13 n 2 + n 4 ) tr ( Σ 2 ) tr ( Σ 3 ) + 6 n ( 2 + 3 n 2 ) tr ( Σ ) tr ( Σ 4 ) + 6 n 2 ( 3 + n 2 ) tr ( Σ 5 ) ,
E tr ( W 4 ) tr ( W ) = n 2 tr ( Σ ) 5 + ( 4 n + 6 n 3 ) tr ( Σ ) 3 tr ( Σ 2 ) + n 2 ( 13 + 2 n 2 ) tr ( Σ ) ( tr ( Σ 2 ) ) 2 + 4 n 2 ( 4 + n 2 ) tr ( Σ ) 2 tr ( Σ 3 ) + 4 n ( 2 + 3 n 2 ) tr ( Σ 2 ) tr ( Σ 3 ) + n ( 12 + 17 n 2 + n 4 ) tr ( Σ ) tr ( Σ 4 ) + 4 n 2 ( 5 + n 2 ) tr ( Σ 5 ) ,
E tr ( W 5 ) = n tr ( Σ ) 5 + 10 n 2 tr ( Σ ) 3 tr ( Σ 2 ) + 5 n ( 1 + 2 n 2 ) tr ( Σ ) ( tr ( Σ 2 ) ) 2 + 10 n ( 1 + n 2 ) tr ( Σ ) 2 tr ( Σ 3 ) + 5 n 2 ( 3 + n 2 ) tr ( Σ 2 ) tr ( Σ 3 ) + 5 n 2 ( 5 + n 2 ) tr ( Σ ) tr ( Σ 4 ) + n ( 8 + 15 n 2 + n 4 ) tr ( Σ 5 ) .

4. Expressions for E tr W h g tr W j i

The expected value of tr W 2 can be derived as
E tr ( W 2 ) = W H = W > 0 etr Σ 1 W det ( W ) n m det ( Σ ) n Γ ˜ m ( n ) ( tr ( W 2 ) ) d W .
Expressing tr X 2 in terms of zonal polynomials (see Appendix A) and integrating the resulting expression with the aid of (7), we can write
E tr ( W 2 ) = C ˜ ( 2 ) Σ 1 [ n + m ] ( 2 ) C ˜ 1 2 Σ 1 [ n + m ] 1 2 , n > m + 1 .
Now, substituting for [ n + m ] ( 2 ) , [ n + m ] 1 2 from Table 2 and for C ˜ ( 2 ) Σ 1 and C ˜ 1 2 Σ 1 from Appendix A, we have
E tr ( W 2 ) = tr Σ 1 2 + ( n m ) tr Σ 2 ( n m 1 ) ( n m ) ( n m + 1 ) , n > m + 1 .
Writing tr W 3 in terms of zonal polynomials (see Appendix A) and applying (7) to integrate the resulting expression, we arrive at
E tr ( W 3 ) = W H = W > 0 etr Σ 1 W det ( W ) n m det ( Σ ) n Γ ˜ m ( n ) ( tr ( W 3 ) ) d X = ( 1 ) 3 C ˜ ( 3 ) Σ 1 [ n + m ] ( 3 ) 1 2 C ˜ ( 2 , 1 ) Σ 1 [ n + m ] ( 2 , 1 ) + C ˜ 1 3 Σ 1 [ n + m ] 1 3 , n > m + 2 .
Now, substituting for [ n + m ] ( 3 ) , [ n + m ] ( 2 , 1 ) , [ n + m ] 1 3 , C ˜ ( 3 ) Σ 1 , C ˜ ( 2 , 1 ) Σ 1 and C ˜ 1 3 Σ 1 above, we have
E tr ( W 3 ) = 2 tr Σ 1 3 + 3 ( n m ) tr Σ 1 tr Σ 2 + ( n m ) 2 tr Σ 3 ( n m 2 ) ( n m 1 ) ( n m ) ( n m + 1 ) ( n m + 2 ) , n > m + 2 .
Similarly, following the procedure described above, we obtain
E tr ( W 2 ) tr ( W 1 ) = tr Σ 1 3 + 2 ( n m ) 1 + n m tr Σ 1 tr Σ 2 + 2 tr Σ 3 ( n m 2 ) ( n m 1 ) ( n m + 1 ) ( n m + 2 ) , n > m + 2 ,
E ( tr ( W 2 ) ) 2 = 1 ( n m ) ( n m 3 ) 7 [ { ( n m ) 2 + 6 } ( tr ( Σ 1 ) ) 4 + 2 ( n m ) { ( n m ) 2 + 6 } ( tr ( Σ 1 ) ) 2 tr ( Σ 2 ) + { ( n m ) 4 6 ( n m ) 2 + 18 } ( tr ( Σ 2 ) ) 2 + 8 { 2 ( n m ) 2 3 } tr ( Σ 1 ) tr ( Σ 3 ) + ( n m ) { 4 ( n m ) 2 6 } tr ( Σ 4 ) ] , n > m + 3 ,
E tr ( W 2 ) ( tr ( W 1 ) ) 2 = 1 ( n m 3 ) 7 [ { ( n m ) 2 4 } ( tr ( Σ 1 ) ) 4 + ( n m ) { ( n m ) 2 + 1 } ( tr ( Σ 1 ) ) 2 tr ( Σ 2 ) + { ( n m ) 2 + 6 } ( tr ( Σ 2 ) ) 2 + 4 { ( n m ) 2 + 1 } tr ( Σ 1 ) tr ( Σ 3 ) + 10 ( n m ) tr ( Σ 4 ) ] , n > m + 3 ,
E tr ( W 3 ) tr ( W 1 ) = 1 ( n m ) ( n m 3 ) 7 [ { 2 ( n m ) 2 3 } ( tr ( Σ 1 ) ) 4 + 3 ( n m ) { ( n m ) 2 + 1 } ( tr ( Σ 1 ) ) 2 tr ( Σ 2 ) + { 6 ( n m ) 2 9 } ( tr ( Σ 2 ) ) 2 + { ( n m ) 4 + 3 ( n m ) 2 + 12 } tr ( Σ 1 ) tr ( Σ 3 ) + 3 ( n m ) { 1 + ( n m ) 2 } tr ( Σ 4 ) ] , n > m + 3 ,
E tr ( W 4 ) = 1 ( n m 3 ) 7 [ 5 ( tr ( Σ 1 ) ) 4 + 10 ( n m ) ( tr ( Σ 1 ) ) 2 tr ( Σ 2 ) + { 2 ( n m ) 2 3 } ( tr ( Σ 2 ) ) 2 + 4 { ( n m ) 2 + 1 } tr ( Σ 1 ) tr ( Σ 3 ) + ( n m ) { ( n m ) 2 + 1 } tr ( Σ 4 ) ] , n > m + 3 ,
E tr ( W 2 ) ( tr ( W 1 ) ) 3 = 1 ( n m ) ( n m 4 ) 9 [ { ( n m ) 4 14 ( n m ) 2 + 24 } ( tr ( Σ 1 ) ) 5 + ( n m ) { 4 + ( n m ) 2 } { ( n m ) 2 9 } ( tr ( Σ 1 ) ) 3 tr ( Σ 2 ) + 3 { ( n m ) 4 + 24 } tr ( Σ 1 ) ( tr ( Σ 2 ) ) 2 + 2 { ( n m ) 2 6 } { 3 ( n m ) 2 + 8 } ( tr ( Σ 1 ) ) 2 tr ( Σ 3 ) + 10 ( n m ) { ( n m ) 2 + 12 } tr ( Σ 2 ) tr ( Σ 3 ) + 30 ( n m ) { ( n m ) 2 2 } tr ( Σ 1 ) tr ( Σ 4 ) + 84 ( n m ) 2 tr ( Σ 5 ) ] , n > m + 4 ,
E ( tr ( W 2 ) ) 2 tr ( W 1 ) = 1 ( n m ) ( n m 4 ) 9 [ ( n m ) { ( n m ) 2 2 } ( tr ( Σ 1 ) ) 5 + 2 { 24 + ( n m ) 4 } ( tr ( Σ 1 ) ) 3 tr ( Σ 2 ) + ( n m ) { ( n m ) 4 10 ( n m ) 2 + 114 } tr ( Σ 1 ) ( tr ( Σ 2 ) ) 2 + 20 ( n m ) { ( n m ) 2 2 } ( tr ( Σ 1 ) ) 2 tr ( Σ 3 ) + 4 { ( n m ) 4 + 24 } tr ( Σ 2 ) tr ( Σ 3 ) + 2 { 2 ( n m ) 4 + 25 ( n m ) 2 72 } tr ( Σ 1 ) tr ( Σ 4 ) + 24 ( n m ) { ( n m ) 2 2 } tr ( Σ 5 ) ] , n > m + 4 ,
E tr ( W 3 ) ( tr ( W 1 ) ) 2 = 1 ( n m ) ( n m 4 ) 9 [ 2 ( n m ) { ( n m ) 2 9 } ( tr ( Σ 1 ) ) 5 + { ( n m ) 2 6 } { 3 ( n m ) 2 + 8 } ( tr ( Σ 1 ) ) 3 tr ( Σ 2 ) + 15 ( n m ) { ( n m ) 2 2 } tr ( Σ 1 ) ( tr ( Σ 2 ) ) 2 + ( n m ) { 24 + ( n m ) 4 } ( tr ( Σ 1 ) ) 2 tr ( Σ 3 ) + { ( n m ) 4 + 60 ( n m ) 2 96 } tr ( Σ 2 ) tr ( Σ 3 ) + 6 { ( n m ) 4 + 24 } tr ( Σ 1 ) tr ( Σ 4 ) + 6 ( n m ) { 8 + 3 ( n m ) 2 } tr ( Σ 5 ) ] , n > m + 4 ,
E tr ( W 3 ) tr ( W 2 ) = 1 ( n m ) ( n m 4 ) 9 [ 2 { ( n m ) 2 + 12 } ( tr ( Σ 1 ) ) 5 + 5 ( n m ) { ( n m ) 2 + 12 } ( tr ( Σ 1 ) ) 3 tr ( Σ 2 ) + 3 { ( n m ) 4 + 24 } tr ( Σ 1 ) ( tr ( Σ 2 ) ) 2 + { ( n m ) 4 + 60 ( n m ) 2 96 } ( tr ( Σ 1 ) ) 2 tr ( Σ 3 ) + ( n m ) { ( n m ) 4 + 24 } tr ( Σ 2 ) tr ( Σ 3 ) + 30 ( n m ) { ( n m ) 2 2 } tr ( Σ 1 ) tr ( Σ 4 ) + 6 ( n m ) 2 { ( n m ) 2 2 } tr ( Σ 5 ) ] , n > m + 4 ,
E tr ( W 4 ) tr ( W 1 ) = 1 ( n m ) ( n m 4 ) 9 [ { 5 ( n m ) 2 24 } ( tr ( Σ 1 ) ) 5 + 10 ( n m ) { ( n m ) 2 2 } ( tr ( Σ 1 ) ) 3 tr ( Σ 2 ) + { 2 ( n m ) 4 + 25 ( n m ) 2 72 } tr ( Σ 1 ) ( tr ( Σ 2 ) ) 2 + 4 { ( n m ) 4 + 24 } ( tr ( Σ 1 ) ) 2 tr ( Σ 3 ) + 20 ( n m ) { ( n m ) 2 2 } tr ( Σ 2 ) tr ( Σ 3 ) + ( n m ) { ( n m ) 4 + 5 ( n m ) 2 + 84 } tr ( Σ 1 ) tr ( Σ 4 ) + 4 ( n m ) 2 { ( n m ) 2 + 5 } tr ( Σ 5 ) ] , n > m + 4 ,
E tr ( W 5 ) = 1 ( n m 4 ) 9 [ 14 ( tr ( Σ 1 ) ) 5 + 35 ( n m ) ( tr ( Σ 1 ) ) 3 tr ( Σ 2 ) + 15 { ( n m ) 2 2 } tr ( Σ 1 ) ( tr ( Σ 2 ) ) 2 + 5 { 3 ( n m ) 2 + 8 } ( tr ( Σ 1 ) ) 2 tr ( Σ 3 ) + 5 ( n m ) { ( n m ) 2 2 } tr ( Σ 2 ) tr ( Σ 3 ) + 5 ( n m ) { ( n m ) 2 + 5 } tr ( Σ 1 ) tr ( Σ 4 ) + ( n m ) 2 { ( n m ) 2 + 5 } tr ( Σ 5 ) ] , n > m + 4 .

5. Special Cases

By substituting Σ = I m into results given in Section 3 and Section 4, several special cases can be derived for W C W m n , I m . Thus, substituting appropriately, we obtain
E tr ( W 2 ) = n m ( m + n ) ,
E tr ( W 3 ) = n m ( m 2 + 3 n m + n 2 + 1 ) ,
E tr ( W 2 ) tr ( W ) = n m ( m + n ) ( m n + 2 ) ,
E ( tr ( W 2 ) ) 2 = n m [ n m 3 + 2 ( n 2 + 2 ) m 2 + n ( n 2 + 10 ) m + 2 ( 2 n 2 + 1 ) ] ,
E tr ( W 2 ) ( tr ( W ) ) 2 = n m ( m + n ) ( 2 + m n ) ( 3 + m n ) ,
E tr ( W 3 ) tr ( W ) = n m ( 3 + m n ) ( 1 + m 2 + 3 m n + n 2 ) ,
E tr ( W 4 ) = n m ( m + n ) ( 5 + m 2 + 5 m n + n 2 ) ,
E ( tr ( W ) ) 3 tr ( W 2 ) = n m ( m + n ) ( 2 + m n ) ( 3 + m n ) ( 4 + m n ) ,
E ( tr ( W 2 ) ) 2 tr ( W ) = n m ( 4 + m n ) [ ( 4 + m n ) ( m + n ) 2 + 2 ( 1 + m n ) ] ,
E tr ( W 3 ) tr ( W ) 2 = n m ( 3 + m n ) ( 4 + m n ) ( 1 + m 2 + 3 m n + n 2 ) ,
E tr ( W 3 ) tr ( W 2 ) = n m ( m + n ) ( 6 + m n ) ( 1 + m 2 + 3 m n + n 2 ) + 6 ( 2 + m n ) ,
E tr ( W 4 ) tr ( W ) = n m ( m + n ) ( 4 + m n ) ( 5 + m 2 + 5 m n + n 2 ) ,
E tr ( W 5 ) = n m ( 5 + m 2 + 5 m n + n 2 ) 2 + 5 ( n m ) 2 7 n 2 m 2 17 ,
E tr ( W 2 ) = m n ( n m 1 ) ( n m ) ( n m + 1 ) , n > m + 1 ,
E tr ( W 3 ) = m n ( m + n ) ( n m 2 ) 5 , n > m + 2 ,
E tr ( W 2 ) tr ( W 1 ) = m n ( m n + 2 m 2 ) ( n m 2 ) 5 , n > m + 2 ,
E ( tr ( W 2 ) ) 2 = m n [ m 3 n ( n 2 2 ) ( 2 m 2 m n 4 ) + 2 ] ( n m ) ( n m 3 ) 7 , n > m + 3 ,
E tr ( W 2 ) ( tr ( W 1 ) ) 2 = m n [ m 2 { ( n m ) 2 9 } + 5 m n + 10 ] ( n m 3 ) 7 , n > m + 3 ,
E tr ( W 3 ) tr ( W 1 ) = m n [ ( n 2 m 2 ) ( 3 + m n m 2 ) 3 ( m 2 1 ) ] ( n m ) ( n m 3 ) 7 , n > m + 3 ,
E tr ( W 4 ) = m n ( m 2 + 3 m n + n 2 + 1 ) ( n m 3 ) 7 , n > m + 3 ,
E tr ( W 2 ) ( tr ( W 1 ) ) 3 = m n ( n m ) ( n m 4 ) 9 [ m 3 ( n m ) 4 + 9 m 2 ( n m ) 3 + m ( 40 14 m 2 ) ( n m ) 2 + ( 84 60 m 2 ) ( n m ) + 24 m ( m 2 1 ) ] , n > m + 4 ,
E ( tr ( W 2 ) ) 2 tr ( W 1 ) = m n ( n m ) ( n m 4 ) 9 [ m 2 ( n m ) 4 + ( m 3 + 8 m ) ( n m ) 3 + ( 2 m 2 + 24 ) ( n m ) 2 + ( 26 m 2 m 3 ) ( n m ) + 48 ( m 2 1 ) ] , n > m + 4 ,
E ( tr ( W 3 ) ( tr ( W 1 ) ) 2 = m n ( n m ) ( n m 4 ) 9 [ m 2 ( n m ) 4 + ( 2 m 3 + 7 m ) ( n m ) 3 + ( 8 m 2 + 18 ) ( n m ) 2 + ( 42 m 18 m 3 ) ( n m ) 48 ( m 2 1 ) ] , n > m + 4 ,
E ( tr ( W 3 ) tr ( W 2 ) = m n ( m + n ) ( n m ) ( n m 4 ) 9 [ m ( n m ) 3 + ( m 2 + 6 ) ( n m ) 2 + 12 ( m n 1 ) ] , n > m + 4 ,
E ( tr ( W 4 ) tr ( W 1 ) = m n ( n m ) ( n m 4 ) 9 [ m ( n m ) 4 + ( 5 m 2 + 4 ) ( n m ) 3 + ( 5 m 3 + 21 m ) ( n m ) 2 + 20 n 20 m 3 + 24 m ] , n > m + 4 ,
E tr ( W 5 ) = m n ( m + n ) ( m 2 + 5 m n + n 2 + 5 ) ( n m 4 ) 9 , n > m + 4 .
Further, using the unitary invariance (Khatri, Khattree, and Gupta [36], Gupta, Nagar, and Vélez-Carvajal [37], Nagar and Gupta [32], and Shaman [7]) of W C W m n , I m , we can easily derive
E ( W 2 ) = n ( m + n ) I m ,
E ( W 3 ) = n ( m 2 + 3 n m + n 2 + 1 ) I m ,
E ( W 4 ) = n ( m + n ) ( 5 + m 2 + 5 m n + n 2 ) I m ,
E ( W 5 ) = n m 4 + 10 n m 3 + 5 ( 3 + 4 n 2 ) m 2 + 10 m n ( 4 + n 2 ) + 8 + 15 n 2 + n 4 I m ,
E ( W 2 ) = n ( n m 1 ) ( n m ) ( n m + 1 ) I m , n > m + 1 ,
E ( W 3 ) = n ( m + n ) ( n m 2 ) ( n m 1 ) ( n m ) ( n m + 1 ) ( n m + 2 ) I m , n > m + 2 ,
E ( W 4 ) = n ( m 2 + 3 m n + n 2 + 1 ) ( n m 3 ) 7 I m , n > m + 3 ,
E ( W 5 ) = n ( m + n ) ( m 2 + 5 m n + n 2 + 5 ) ( n m 4 ) 9 I m , n > m + 4 .

6. Discussion and Conclusions

In Section 3 and Section 4, we evaluated as many as twenty-six expected values of the type E [ ( tr ( W h ) ) g ( tr ( W j ) ) i ] and E [ ( tr ( W h ) ) g ( tr ( W j ) ) i ] , W C W m n , Σ . In Section 5, these expected values were simplified for W C W m n , I m . Further, using the unitary invariance of W C W m n , I m , we also computed a number of expected values of matrix-valued functions of W.
For k 5 , explicit formulas for C ˜ κ ( X ) available in Khatri [35] are reproduced in Appendix A. By using these expressions for zonal polynomials, through a linear equation system, we expressed traces such as ( tr ( X h ) ) g ( tr ( X j ) ) i , g + h + i + j 5 in terms of zonal polynomials. These results are also summarized in Appendix A. Further, using these expressions, for specific values of g, h, i, j, and Lemmas 1 and 2 adequately, we computed several expected values given in Section 3 and Section 4. Clearly, the method applied in this paper is simple, straightforward, gives explicit expressions for expected values of functions of W and W 1 , and does not require any advanced mathematical tool. The present method requires results on zonal polynomials and two lemmas. To find moments of higher order, for example E tr W 6 , we will require zonal polynomials for k = 6 .

Author Contributions

Conceptualization, D.K.N., A.R.-C. and S.N.; methodology, D.K.N., A.R.-C. and S.N.; writing—original draft preparation, D.K.N., A.R.-C. and S.N.; writing—review and editing, D.K.N., A.R.-C. and S.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

For small values of k, explicit formulas for C ˜ κ ( X ) are available in Khatri [35] as
C ˜ ( 1 ) ( X ) = tr ( X ) ,
C ˜ ( 2 ) ( X ) = 1 2 ( tr ( X ) ) 2 + tr ( X 2 ) ,
C ˜ ( 1 2 ) ( X ) = 1 2 ( tr ( X ) ) 2 tr ( X 2 ) ,
C ˜ ( 3 ) ( X ) = 1 6 ( tr ( X ) ) 3 + 3 tr ( X ) tr ( X 2 ) + 2 tr ( X 3 ) ,
C ˜ ( 2 , 1 ) ( X ) = 2 3 ( tr ( X ) ) 3 tr ( X 3 ) ,
C ˜ ( 1 3 ) ( X ) = 1 6 ( tr ( X ) ) 3 3 tr ( X ) tr ( X 2 ) + 2 tr ( X 3 ) .
C ˜ ( 4 ) ( X ) = 1 24 ( tr ( X ) ) 4 + 6 ( tr ( X ) ) 2 tr ( X 2 ) + 3 ( tr ( X 2 ) ) 2 + 8 tr ( X ) tr ( X 3 ) + 6 tr ( X 4 ) ,
C ˜ ( 3 , 1 ) ( X ) = 3 24 3 ( tr ( X ) ) 4 + 6 ( tr ( X ) ) 2 tr ( X 2 ) 3 ( tr ( X 2 ) ) 2 6 tr ( X 4 ) ,
C ˜ ( 2 2 ) ( X ) = 2 24 2 ( tr ( X ) ) 4 + 6 ( tr ( X 2 ) ) 2 8 tr ( X ) tr ( X 3 ) .
C ˜ ( 2 , 1 2 ) ( X ) = 3 24 3 ( tr ( X ) ) 4 6 ( tr ( X ) ) 2 tr ( X 2 ) 3 ( tr ( X 2 ) ) 2 + 6 tr ( X 4 ) ,
C ˜ ( 1 4 ) ( X ) = 1 24 ( tr ( X ) ) 4 6 ( tr ( X ) ) 2 tr ( X 2 ) + 3 ( tr ( X 2 ) ) 2 + 8 tr ( X ) tr t ( X 3 ) 6 tr t ( X 4 ) ,
C ˜ ( 5 ) ( X ) = 1 120 ( tr ( X ) ) 5 + 10 ( tr ( X ) ) 3 tr ( X 2 ) + 15 tr ( X ) ( tr ( X 2 ) ) 2 + 20 ( tr ( X ) ) 2 tr ( X 3 ) + 20 tr ( X 2 ) tr ( X 3 ) + 30 ( tr ( X ) ) tr ( X 4 ) + 24 tr ( X 5 ) ,
C ˜ ( 4 , 1 ) ( X ) = 4 120 4 ( tr ( X ) ) 5 + 20 ( tr ( X ) ) 3 tr ( X 2 ) + 20 ( tr ( X ) ) 2 tr ( X 3 ) 20 tr ( X 2 ) tr ( X 3 ) 24 tr ( X 5 ) ,
C ˜ ( 3 , 2 ) ( X ) = 5 120 [ 5 ( tr ( X ) ) 5 + 10 ( tr ( X ) ) 3 tr ( X 2 ) + 15 tr ( X ) ( tr ( X 2 ) ) 2 20 ( tr ( X ) ) 2 tr ( X 3 ) + 20 tr ( X 2 ) tr ( X 3 ) 30 tr ( X ) tr ( X 4 ) ] ,
C ˜ ( 3 , 1 2 ) ( X ) = 6 120 6 ( tr ( X ) ) 5 30 tr ( X ) ( tr ( X 2 ) ) 2 + 24 tr ( X 5 ) ,
C ˜ ( 2 2 , 1 ) ( X ) = 5 120 [ 5 ( tr ( X ) ) 5 10 ( tr ( X ) ) 3 tr ( X 2 ) + 15 tr ( X ) ( tr ( X 2 ) ) 2 20 ( tr ( X ) ) 2 tr ( X 3 ) 20 tr ( X 2 ) tr ( X 3 ) + 30 tr ( X ) tr ( X 4 ) ] ,
C ˜ ( 2 , 1 3 ) ( X ) = 4 120 [ 4 ( tr ( X ) ) 5 20 ( tr ( X ) ) 3 tr ( X 2 ) + 20 ( tr ( X ) ) 2 tr ( X 3 ) + 20 tr t ( X 2 ) tr ( X 3 ) 24 tr ( X 5 ) ] ,
C ˜ ( 1 5 ) ( X ) = 1 120 [ ( tr ( X ) ) 5 10 ( tr ( X ) ) 3 tr ( X 2 ) + 15 tr ( X ) ( tr ( X 2 ) ) 2 + 20 ( tr ( X ) ) 2 tr ( X 3 ) 20 tr ( X 2 ) tr ( X 3 ) 30 tr ( X ) tr ( X 4 ) + 24 tr ( X 5 ) ] .
By using the above expressions for zonal polynomials it is straightforward, through a linear equation system, to express traces such as tr X h g tr X j i in terms of zonal polynomials. Thus, after some algebra, we obtain
tr ( X 2 ) = C ˜ ( 2 ) ( X ) C ˜ 1 2 ( X ) ,
tr ( X 3 ) = C ˜ ( 3 ) ( X ) 1 2 C ˜ ( 2 , 1 ) ( X ) + C ˜ 1 3 ( X ) ,
tr ( X ) tr ( X 2 ) = C ˜ ( 3 ) ( X ) C ˜ 1 3 ( X ) ,
( tr ( X ) ) 2 tr ( X 2 ) = C ˜ ( 4 ) ( X ) + 1 3 C ˜ ( 3 , 1 ) ( X ) 1 3 C ˜ 2 , 1 2 ( X ) C ˜ 1 4 ( X ) ,
( tr ( X 2 ) ) 2 = C ˜ ( 4 ) ( X ) 1 3 C ˜ ( 3 , 1 ) ( X ) + C ˜ 2 2 ( X ) 1 3 C ˜ 2 , 1 2 ( X ) + C ˜ 1 4 ( X ) ,
tr ( X ) tr ( X 3 ) = C ˜ ( 4 ) ( X ) 1 2 C ˜ 2 2 ( X ) + C ˜ 1 4 ( X ) ,
tr ( X 4 ) = C ˜ ( 4 ) ( X ) 1 3 C ˜ ( 3 , 1 ) ( X ) + 1 3 C ˜ 2 , 1 2 ( X ) C ˜ 1 4 ( X ) ,
( tr ( X ) ) 3 tr ( X 2 ) = C ˜ ( 5 ) ( X ) + 1 2 C ˜ ( 4 , 1 ) ( X ) + 1 5 C ˜ ( 3 , 2 ) ( X ) 1 5 C ˜ ( 2 2 , 1 ) ( X ) 1 2 C ˜ ( 2 , 1 3 ) ( X ) C ˜ ( 1 5 ) ( X ) ,
tr ( X ) ( tr ( X 2 ) ) 2 = C ˜ ( 5 ) ( X ) + 1 5 C ˜ ( 3 , 2 ) ( X ) 1 3 C ˜ 3 , 1 2 ( X ) + 1 5 C ˜ 2 2 , 1 ( X ) + C ˜ 1 5 ( X ) ,
( tr ( X ) ) 2 tr ( X 3 ) = C ˜ ( 5 ) ( X ) + 1 4 C ˜ ( 4 , 1 ) ( X ) 1 5 C ˜ ( 3 , 2 ) ( X ) 1 5 C ˜ ( 2 2 , 1 ) ( X ) + 1 4 C ˜ ( 2 , 1 3 ) ( X ) + C ˜ ( 1 5 ) ( X ) ,
tr ( X 2 ) tr ( X 3 ) = C ˜ ( 5 ) ( X ) 1 4 C ˜ ( 4 , 1 ) ( X ) + 1 5 C ˜ ( 3 , 2 ) ( X ) 1 5 C ˜ ( 2 2 , 1 ) ( X ) + 1 4 C ˜ ( 2 , 1 3 ) ( X ) C ˜ ( 1 5 ) ( X ) ,
tr ( X ) tr ( X 4 ) = C ˜ ( 5 ) ( X ) 1 5 C ˜ ( 3 , 2 ) ( X ) + 1 5 C ˜ 2 2 , 1 ( X ) C ˜ 1 5 ( X ) ,
tr ( X 5 ) = C ˜ ( 5 ) ( X ) 1 4 C ˜ ( 4 , 1 ) ( X ) + 1 6 C ˜ 3 , 1 2 ( X ) 1 4 C ˜ 2 , 1 3 ( X ) + C ˜ 1 5 ( X ) .

References

  1. Goodman, N.R. Statistical analysis based on a certain multivariate complex Gaussian distribution (An introduction). Ann. Math. Stat. 1963, 34, 152–177. [Google Scholar] [CrossRef]
  2. James, A.T. Distributions of matrix variate and latent roots derived from normal samples. Ann. Math. Stat. 1964, 35, 475–501. [Google Scholar] [CrossRef]
  3. Khatri, C.G. Classical statistical analysis based on a certain multivariate complex Gaussian distribution. Ann. Math. Stat. 1965, 36, 98–114. [Google Scholar] [CrossRef]
  4. Khatri, C.G. On certain distribution problems based on positive definite quadratic functions in normal vectors. Ann. Math. Stat. 1966, 37, 468–479. [Google Scholar] [CrossRef]
  5. Krishnaiah, P.R. Some recent developments on complex multivariate distributions. J. Multivar. Anal. 1976, 6, 1–30. [Google Scholar] [CrossRef]
  6. Srivastava, M.S. On the complex Wishart distribution. Ann. Math. Stat. 1965, 36, 313–315. [Google Scholar] [CrossRef]
  7. Shaman, P. The inverted complex Wishart distribution and its application to spectral estimation. J. Multivar. Anal. 1980, 10, 51–59. [Google Scholar] [CrossRef]
  8. Tan, W.Y. Some distribution theory associated with complex Gaussian distribution. Tamkang J. 1968, 7, 263–302. [Google Scholar]
  9. Alfano, G.; Tulino, A.; Lozano, A.; Verdú, S. Random Matrix Transforms and Applications via Non-Asymptotic Eigenanalysis; International Zurich Seminar on Communications: Zurich, Switzerland, 2006; pp. 18–21. [Google Scholar]
  10. Carmeli, M. Statistical Theory and Random Matrices in Physics; Marcel Dekker: New York, NY, USA, 1983. [Google Scholar]
  11. Cunden, F.D.; Dahlqvist, A.; O’Connell, N. Integer moments of complex Wishart matrices and Hurwitz numbers. Ann. Inst. Henri Poincaré D 2021, 8, 243–268. [Google Scholar] [CrossRef]
  12. Deng, X.; López-Martínez, C.; Chen, J.; Han, P. Statistical Modeling of Polarimetric SAR Data: A Survey and Challenges. Remote Sens. 2008, 9, 348. [Google Scholar] [CrossRef]
  13. Ermolova, N.Y.; Tirkkonen, O. Distribution of diagonal elements of a general central complex Wishart matrix. IEEE Commun. Lett. 2012, 16, 1373–1376. [Google Scholar] [CrossRef]
  14. Gomez, L.; Alvarez, L.; Mazorra, L.; Frery, A.C. Classification of complex Wishart matrices with a diffusion-reaction system guided by stochastic distances. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 2015, 373, 20150118. [Google Scholar] [CrossRef] [PubMed]
  15. Kumar, S. Eigenvalue statistics for the sum of two complex Wishart matrices. EPL 2014, 107, 60002. [Google Scholar] [CrossRef]
  16. Nielsen, A.A.; Skriver, H.; Conradsen, K. Complex Wishart Distribution Based Analysis of Polarimetric Synthetic Aperture Radar Data. In Proceedings of the 2007 International Workshop on the Analysis of Multi-Temporal Remote Sensing Images, Leuven, Belgium, 18–20 July 2007; pp. 1–6. [Google Scholar]
  17. Shakil, M.; Ahsanullah, M. Some inferences on the distribution of the Demmel condition number of complex Wishart matrices. Spec. Matrices 2017, 5, 127–138. [Google Scholar] [CrossRef]
  18. Tague, J.A.; Caldwell, C.I. Expectations of useful complex Wishart forms. Multidimens. Syst. Signal Process. 1994, 5, 263–279. [Google Scholar] [CrossRef]
  19. Tralli, V.; Conti, A. Noncentral complex Wishart matrices: Moments and correlation of minors. Random Matrices Theory Appl. 2022, 11, 2250009. [Google Scholar] [CrossRef]
  20. Andersen, H.H.; Højbjerre, M.; Sørensen, D.; Eriksen, P.S. Linear and Graphical Models; For the Multivariate Complex Normal Distribution, Lecture Notes in Statistics, 101; Springer: New York, NY, USA, 1995. [Google Scholar]
  21. Gupta, A.K.; Nagar, D.K. Matrix Variate Distributions; Chapman and Hall/CRC: Boca Raton, FL, USA, 2000. [Google Scholar]
  22. Nagar, D.K.; Gupta, A.K.; Sánchez, L.E. A class of integral identities with Hermitian matrix argument. Proc. Am. Math. Soc. 2006, 134, 3329–3341. [Google Scholar] [CrossRef]
  23. Letac, G.; Maasam, H. The noncentral Wishart as an exponential family, and its moments. J. Multivar. Anal. 2008, 99, 1393–1417. [Google Scholar] [CrossRef]
  24. Nagar, D.K.; Roldán-Correa, A.; Gupta, A.K. Extended matrix variate gamma and beta functions. J. Multivar. Anal. 2013, 122, 53–69. [Google Scholar] [CrossRef]
  25. Di Nardo, E. On a symbolic representation of non-central Wishart random matrices with applications. J. Multivar. Anal. 2014, 125, 121–135. [Google Scholar] [CrossRef]
  26. Dharmawansa, P.; McKay, M.R. Extreme eigenvalue distributions of some complex correlated non-central Wishart and gamma-Wishart random matrices. J. Multivar. Anal. 2011, 102, 847–868. [Google Scholar] [CrossRef]
  27. Hillier, G.; Kan, R. Properties of the inverse of a noncentral Wishart matrix. Econom. Theory 2022, 38, 1092–1116. [Google Scholar] [CrossRef]
  28. Wahba, G. Cross-Spectral Distribution Theory for Mixed Spectra and Estimation of Prediction Filter Coefficients. Ph.D. Thesis, Stanford University, Stanford, CA, USA, 1966. [Google Scholar]
  29. Sultan, S.A.; Tracy, D.S. Moments of the complex Wishart distribution. Int. J. Math. Stat. Sci. 1999, 8, 51–62. [Google Scholar]
  30. Maiwald, D.; Kraus, D. Calculation of moments of complex Wishart and complex inverse Wishart distributed matrices. IEE Proc.-Radar Sonar Navig. 2000, 147, 162–168. [Google Scholar] [CrossRef]
  31. Graczyk, P.; Letac, G.; Massam, H. The complex Wishart distribution and the symmetric group. Ann. Stat. 2003, 31, 287–309. [Google Scholar] [CrossRef]
  32. Nagar, D.K.; Gupta, A.K. Expectations of functions of complex Wishart matrix. Acta Appl. Math. 2011, 113, 265–288. [Google Scholar] [CrossRef]
  33. Pielaszkiewicz, J.; von Rosen, D.; Singull, M. On E i = 0 k Tr W m i , where W W p ( I , n ) . Commun. Stat.—Theory Methods 2017, 46, 2990–3005. [Google Scholar] [CrossRef]
  34. Holgersson, T.; Pielaszkiewicz, J. A Collection of Moments of the Wishart Distribution. In Recent Developments in Multivariate and Random Matrix Analysis; Holgersson, T., Singull, M., Eds.; Springer: Cham, Switzerland, 2020; pp. 147–162. [Google Scholar]
  35. Khatri, C.G. On the moments of traces of two matrices in three situations for complex multivariate normal populations. Sankhyā Indian J. Stat. Ser. A 1970, 32, 65–80. [Google Scholar]
  36. Khatri, C.G.; Khattree, R.; Gupta, R.D. On a class of orthogonal invariant and residual independent matrix distributions. Sankhyā Indian J. Stat. Ser. B 1991, 53, 1–10. [Google Scholar]
  37. Gupta, A.K.; Nagar, D.K.; Vélez-Carvajal, A.M. Unitary invariant and residual independent matrix distributions. Comput. Appl. Math. 2009, 28, 63–86. [Google Scholar] [CrossRef]
Table 1. Values of [ a ] κ .
Table 1. Values of [ a ] κ .
[ a ] κ i = 1 ( κ ) a i + 1 k i
[ a ] ( 1 ) a
[ a ] ( 2 ) a ( a + 1 )
[ a ] 1 2 a ( a 1 )
[ a ] ( 3 ) a ( a + 1 ) ( a + 2 )
[ a ] ( 2 , 1 ) a ( a + 1 ) ( a 1 )
[ a ] 1 3 a ( a 1 ) ( a 2 )
[ a ] ( 4 ) a ( a + 1 ) ( a + 2 ) ( a + 3 )
[ a ] ( 3 , 1 ) a ( a + 1 ) ( a + 2 ) ( a 1 )
[ a ] 2 2 a 2 ( a + 1 ) ( a 1 )
[ a ] 2 , 1 2 a ( a + 1 ) ( a 1 ) ( a 2 )
[ a ] 1 4 a ( a 1 ) ( a 2 ) ( a 3 )
[ a ] ( 5 ) a ( a + 1 ) ( a + 2 ) ( a + 3 ) ( a + 4 )
[ a ] ( 4 , 1 ) a ( a + 1 ) ( a + 2 ) ( a + 3 ) ( a 1 )
[ a ] ( 3 , 2 ) a 2 ( a + 1 ) ( a + 2 ) ( a 1 )
[ a ] 3 , 1 2 a ( a + 1 ) ( a + 2 ) ( a 1 ) ( a 2 )
[ a ] 2 2 , 1 a 2 ( a + 1 ) ( a 1 ) ( a 2 )
[ a ] 2 , 1 3 a ( a + 1 ) ( a 1 ) ( a 2 ) ( a 3 )
[ a ] 1 5 a ( a 1 ) ( a 2 ) ( a 3 ) ( a 4 )
Table 2. Values of [ a + m ] κ .
Table 2. Values of [ a + m ] κ .
[ a + m ] κ i = 1 ( κ ) a + m i + 1 k i
[ a + m ] ( 1 ) ( a m )
[ a + m ] ( 2 ) ( a m 1 ) ( a m )
[ a + m ] 1 2 ( a m ) ( a m + 1 )
[ a + m ] ( 3 ) ( a m 2 ) ( a m 1 ) ( a m )
[ a + m ] ( 2 , 1 ) ( a m 1 ) ( a m ) ( a m + 1 )
[ a + m ] 1 3 ( a m ) ( a m + 1 ) ( a m + 2 )
[ a + m ] ( 4 ) ( a m 3 ) ( a m 2 ) ( a m 1 ) ( a m )
[ a + m ] ( 3 , 1 ) ( a m 2 ) ( a m 1 ) ( a m ) ( a m + 1 )
[ a + m ] 2 2 ( a m 1 ) ( a m ) 2 ( a m + 1 )
[ a + m ] 2 , 1 2 ( a m 1 ) ( a m ) ( a m + 1 ) ( a m + 2 )
[ a + m ] 1 4 ( a m ) ( a m + 1 ) ( a m + 2 ) ( a m + 3 )
[ a + m ] ( 5 ) ( a m 4 ) ( a m 3 ) ( a m 2 ) ( a m 1 ) ( a m )
[ a + m ] ( 4 , 1 ) ( a m 3 ) ( a m 2 ) ( a m 1 ) ( a m ) ( a m + 1 )
[ a + m ] ( 3 , 2 ) ( a m 2 ) ( a m 1 ) ( a m ) 2 ( a m + 1 )
[ a + m ] 3 , 1 2 ( a m 2 ) ( a m 1 ) ( a m ) ( a m + 1 ) ( a m + 2 )
[ a + m ] 2 2 , 1 ( a m 1 ) ( a m ) 2 ( a m + 1 ) ( a m + 2 )
[ a + m ] 2 , 1 3 ( a m 1 ) ( a m ) ( a m + 1 ) ( a m + 2 ) ( a m + 3 )
[ a + m ] 1 5 ( a m ) ( a m + 1 ) ( a m + 2 ) ( a m + 3 ) ( a m + 4 )
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Nagar, D.K.; Roldán-Correa, A.; Nadarajah, S. Expected Values of Scalar-Valued Functions of a Complex Wishart Matrix. Mathematics 2023, 11, 2162. https://doi.org/10.3390/math11092162

AMA Style

Nagar DK, Roldán-Correa A, Nadarajah S. Expected Values of Scalar-Valued Functions of a Complex Wishart Matrix. Mathematics. 2023; 11(9):2162. https://doi.org/10.3390/math11092162

Chicago/Turabian Style

Nagar, Daya K., Alejandro Roldán-Correa, and Saralees Nadarajah. 2023. "Expected Values of Scalar-Valued Functions of a Complex Wishart Matrix" Mathematics 11, no. 9: 2162. https://doi.org/10.3390/math11092162

APA Style

Nagar, D. K., Roldán-Correa, A., & Nadarajah, S. (2023). Expected Values of Scalar-Valued Functions of a Complex Wishart Matrix. Mathematics, 11(9), 2162. https://doi.org/10.3390/math11092162

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop