Multivariate Poisson intensity estimation via low-rank tensor decomposition

In this work, we introduce new matrix- and tensor-based methodologies for estimating multivariate intensity functions of spatial point processes. By modeling intensity functions as infinite-rank tensors within function spaces, we develop new algorithms to reveal optimal bias-variance trade-off for infinite-rank tensor estimation. Our methods dramatically enhance estimation accuracy while simultaneously reducing computational complexity. To our knowledge, this work marks the first application of matrix and tensor techinques to spatial point processes. Extensive numerical experiments further demonstrate that our techniques consistently outperform current state-of-the-art methods.