English

Optimal Randomized Approximations for Matrix based Renyi's Entropy

Machine Learning 2022-05-17 v1 Machine Learning

Abstract

The Matrix-based Renyi's entropy enables us to directly measure information quantities from given data without the costly probability density estimation of underlying distributions, thus has been widely adopted in numerous statistical learning and inference tasks. However, exactly calculating this new information quantity requires access to the eigenspectrum of a semi-positive definite (SPD) matrix AA which grows linearly with the number of samples nn, resulting in a O(n3)O(n^3) time complexity that is prohibitive for large-scale applications. To address this issue, this paper takes advantage of stochastic trace approximations for matrix-based Renyi's entropy with arbitrary αR+\alpha \in R^+ orders, lowering the complexity by converting the entropy approximation to a matrix-vector multiplication problem. Specifically, we develop random approximations for integer order α\alpha cases and polynomial series approximations (Taylor and Chebyshev) for non-integer α\alpha cases, leading to a O(n2sm)O(n^2sm) overall time complexity, where s,mns,m \ll n denote the number of vector queries and the polynomial order respectively. We theoretically establish statistical guarantees for all approximation algorithms and give explicit order of s and m with respect to the approximation error ε\varepsilon, showing optimal convergence rate for both parameters up to a logarithmic factor. Large-scale simulations and real-world applications validate the effectiveness of the developed approximations, demonstrating remarkable speedup with negligible loss in accuracy.

Keywords

Cite

@article{arxiv.2205.07426,
  title  = {Optimal Randomized Approximations for Matrix based Renyi's Entropy},
  author = {Yuxin Dong and Tieliang Gong and Shujian Yu and Chen Li},
  journal= {arXiv preprint arXiv:2205.07426},
  year   = {2022}
}
R2 v1 2026-06-24T11:18:03.411Z