English

Estimating Higher-Order Moments Using Symmetric Tensor Decomposition

Numerical Analysis 2020-10-06 v2 Numerical Analysis

Abstract

We consider the problem of decomposing higher-order moment tensors, i.e., the sum of symmetric outer products of data vectors. Such a decomposition can be used to estimate the means in a Gaussian mixture model and for other applications in machine learning. The ddth-order empirical moment tensor of a set of pp observations of nn variables is a symmetric dd-way tensor. Our goal is to find a low-rank tensor approximation comprising rpr \ll p symmetric outer products. The challenge is that forming the empirical moment tensors costs O(pnd)O(pn^d) operations and O(nd)O(n^d) storage, which may be prohibitively expensive; additionally, the algorithm to compute the low-rank approximation costs O(nd)O(n^d) per iteration. Our contribution is avoiding formation of the moment tensor, computing the low-rank tensor approximation of the moment tensor implicitly using O(pnr)O(pnr) operations per iteration and no extra memory. This advance opens the door to more applications of higher-order moments since they can now be efficiently computed. We present numerical evidence of the computational savings and show an example of estimating the means for higher-order moments.

Keywords

Cite

@article{arxiv.1911.03813,
  title  = {Estimating Higher-Order Moments Using Symmetric Tensor Decomposition},
  author = {Samantha Sherman and Tamara G. Kolda},
  journal= {arXiv preprint arXiv:1911.03813},
  year   = {2020}
}
R2 v1 2026-06-23T12:10:29.682Z