88
Views
0
CrossRef citations to date
0
Altmetric
Research Article

Computationally efficient techniques for spatial regression with differential regularization

ORCID Icon, ORCID Icon, ORCID Icon, & ORCID Icon
Pages 1971-1991 | Received 22 Nov 2022, Accepted 04 Jul 2023, Published online: 03 Aug 2023
 

Abstract

We investigate some computational aspects of an innovative class of PDE-regularized statistical models: Spatial Regression with Partial Differential Equation regularization (SR-PDE). These physics-informed regression methods can account for the physics of the underlying phenomena and handle data observed over spatial domains with nontrivial shapes, such as domains with concavities and holes or curved domains. The computational bottleneck in SR-PDE estimation is the solution of a computationally demanding linear system involving a low-rank but dense block. We address this aspect by innovatively using Sherman–Morrison–Woodbury identity. We also investigate the efficient selection of the smoothing parameter in SR-PDE estimates. Specifically, we propose ad hoc optimization methods to perform Generalized Cross-Validation, coupling suitable reformulation of key matrices, e.g. those based on Sherman–Morrison–Woodbury formula, with stochastic trace estimation, to approximate the equivalent degrees of freedom of the problem. These solutions permit high computational efficiency also in the context of massive data.

2020 AMS Subject Classification::

Acknowledgments

This work is based on preliminary explorations by Clelia Bambini and Luca Giussani. The authors also thank Michelle Carey for discussions on these topics.

Disclosure statement

No potential conflict of interest was reported by the author(s).

Notes

Additional information

Funding

C.d.F. and L.F. have been partially funded by the Italian Research Center on High-Performance Computing, Big Data and Quantum Computing (ICSC), European Union - Next Generation EU. The present research is within the framework of MUR grant Dipartimento di Eccellenza 2023–2027, Dipartimento di Matematica, Politenico di Milano.

Reprints and Corporate Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

To request a reprint or corporate permissions for this article, please click on the relevant link below:

Academic Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

Obtain permissions instantly via Rightslink by clicking on the button below:

If you are unable to obtain permissions via Rightslink, please complete and submit this Permissions form. For more information, please visit our Permissions help page.