GENERATING FUNCTIONS AND SHORT RECURSIONS, WITH APPLICATIONS TO THE MOMENTS OF QUADRATIC FORMS IN NONCENTRAL NORMAL VECTORS
推广了Hillier等人(2009)的结果，证明了一类更广泛的统计对象存在短递归，并应用于非中心正态向量二次型的矩、乘积矩和幂的期望，进一步通过Horner方法得到超短递归，显著提高计算效率。
Recursive relations for objects of statistical interest have long been important for computation, and they remain so even with hugely improved computing power. Such recursions are frequently derived by exploiting relations between generating functions. For example, the top-order zonal polynomials that occur in much distribution theory under normality can be recursively related to other (easily computed) symmetric functions (power-sum and elementary symmetric functions; Ruben, 1962, Annals of Mathematical Statistics 33, 542–570; Hillier, Kan, and Wang, 2009, Econometric Theory 25, 211–242). Typically, in a recursion of this type the k th object of interest, d k , say, is expressed in terms of all lower order d j ’s. In Hillier et al. (2009) we pointed out that, in the case of top-order zonal polynomials and other invariant polynomials of multiple matrix argument, a fixed length recursion can be deduced. We refer to this as a short recursion. The present paper shows that the main results in Hillier et al. (2009) can be generalized and that short recursions can be obtained for a much larger class of objects/generating functions. As applications, we show that short recursions can be obtained for various problems involving quadratic forms in noncentral normal vectors, including moments, product moments, and expectations of ratios of powers of quadratic forms. For this class of problems, we also show that the length of the recursion can be further reduced by an application of a generalization of Horner’s method (cf. Brown, 1986, SIAM Journal on Scientific and Statistical Computing 7, 689–695), producing a super-short recursion that is significantly more efficient than even the short recursion.