An identity for the Wishart distribution with applications

AbstractLet Sp×p have a Wishart distribution with unknown matrix Σ and k degrees of freedom. For a matrix T(S) and a scalar h(S), an identity is obtained for E∑tr[h(S)T∑−1]. Two applications are given. The first provides product moments and related formulae for the Wishart distribution. Higher moments involving S can be generated recursively. The second application concerns good estimators of ∑ and ∑−1. In particular, identities for several risk functions are obtained, and estimators of ∑ (∑−1) are described which dominate aS(bS−1), a ≤ 1k (b ≤ k − p − 1). Haff [(1977) J. Multivar. Anal. 7 374–385; (1979) Ann. Statist. 7 No. 5; (1980) Ann. Statist. 8 used special cases of the identity to find unbiased risk estimators. These are unobtainable in closed form for certain natural loss functions. In this paper, we treat these case as well. The dominance results provide a unified theory for the estimation of ∑ and ∑−1