4 files

Generalized Low-Rank Plus Sparse Tensor Estimation by Fast Riemannian Optimization

posted on 13.05.2022, 18:40 by Jian-Feng Cai, Jingyang Li, Dong Xia

We investigate a generalized framework to estimate a latent low-rank plus sparse tensor, where the low-rank tensor often captures the multi-way principal components and the sparse tensor accounts for potential model mis-specifications or heterogeneous signals that are unexplainable by the low-rank part. The framework flexibly covers both linear and generalized linear models, and can easily handle continuous or categorical variables. We propose a fast algorithm by integrating the Riemannian gradient descent and a novel gradient pruning procedure. Under suitable conditions, the algorithm converges linearly and can simultaneously estimate both the low-rank and sparse tensors. The statistical error bounds of final estimates are established in terms of the gradient of loss function. The error bounds are generally sharp under specific statistical models, for example, the sub-Gaussian robust PCA and Bernoulli tensor model. Moreover, our method achieves nontrivial error bounds for heavy-tailed tensor PCA whenever the noise has a finite 2+ε moment. We apply our method to analyze the international trade flow dataset and the statistician hypergraph coauthorship network, both yielding new and interesting findings. Supplementary materials for this article are available online.


Jian-Feng Cai’s research was partially supported by Hong Kong RGC grant GRF 16310620 and GRF 16309219. Dong Xia’s research was partially supported by Hong Kong RGC grant ECS 26302019 and GRF 16303320, 16300121.