On the Wick theorem for mixtures of centered Gaussian distributions

We consider centered conditionally Gaussian d-dimensional vectors X with random covariance matrix Ξ having an arbitrary probability distribution law on the set of nonnegative definite symmetric d × d matrices M _d ⁺. The paper deals with the evaluation problem of mean values \( E\left {\prod\nolimits_{i = 1}^{2n} {\left( {{c_i},X} \right)} } \right] \) for c _i ∈ ? ^d , i = 1, …, 2n, extending the Wick theorem for a wide class of non-Gaussian distributions. We discuss in more detail the cases where the probability law ?(Ξ) is infinitely divisible, the Wishart distribution, or the inverse Wishart distribution. An example with Ξ \( = \sum\nolimits_{j = 1}^m {{Z_j}{\sum_j}} \), where random variables Z _j , j = 1, …, m, are nonnegative, and Σ _j ∈ M _d ⁺, j = 1, …, m, are fixed, includes recent results from Vignat and Bhatnagar, 2008.