Higher order moments of multivariate normal distribution using matrix derivatives

Abstract

A general formula for the central moments of multivariate normal distribution is derived by differentiating its characteristic function using matrix derivatives. An explicit expression for the moments is obtained. Two applications of these results are given. The sixth order moments are arranged in a square matrix using the properties of commutation matrices and vec operators