A Tutorial on Multivariate -Statistics and their Computation

This document aims to provide an accessible tutorial on the unbiased estimation of multivariate cumulants, using -statistics. We offer an explicit and general formula for multivariate -statistics of arbitrary order. We also prove that the -statistics are unbiased, using M\"obius inversion and rudimentary combinatorics. Many detailed examples are considered throughout the paper. We conclude with a discussion of -statistics computation, including the challenge of time complexity, and we examine a couple of possible avenues to improve the efficiency of this computation. The purpose of this document is threefold: to provide a clear introduction to -statistics without relying on specialized tools like the umbral calculus; to construct an explicit formula for -statistics that might facilitate future approximations and faster algorithms; and to serve as a companion paper to our Python library PyMoments, which implements this formula.
View on arXiv