{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,5,29]],"date-time":"2026-05-29T12:15:14Z","timestamp":1780056914088,"version":"3.54.0"},"reference-count":57,"publisher":"Frontiers Media SA","license":[{"start":{"date-parts":[[2026,1,22]],"date-time":"2026-01-22T00:00:00Z","timestamp":1769040000000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"content-domain":{"domain":["frontiersin.org"],"crossmark-restriction":true},"short-container-title":["Front. Appl. Math. Stat."],"abstract":"<jats:sec>\n                    <jats:title>Introduction<\/jats:title>\n                    <jats:p>We extend intrinsic local polynomial regression (ILPR) to multivariate covariates on arbitrary Riemannian manifolds, providing a general framework to study the action of isometric maps on manifold-valued regression estimators.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Methods<\/jats:title>\n                    <jats:p>We derive a closed-form expression for the ILPR estimator on manifolds endowed with an Euclidean pullback metric (EPM), a broad class that encompasses infinitely many Riemannian metrics. Pseudocode is provided along with an explicit upper bound on computational complexity. For the multivariate local linear case, we characterize the asymptotic properties of the ILPR\u2013EPM estimator, including a closed-form expression for its bias, establishing statistical consistency under EPM.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Results<\/jats:title>\n                    <jats:p>Simulation studies on the manifold of symmetric positive definite (SPD) matrices\u2014benchmarking diffusion tensor imaging (DTI) data and high-dimensional covariance models\u2014show that the Log\u2013Cholesky metric achieves a favorable trade-off between estimation accuracy and computational efficiency. Using real neurophysiological data from the Harmonized-MultiNational quantitative electroencephalography (HarMNqEEG) dataset (N = 1,965), we apply ILPR to symmetric positive definite cross-spectral covariance matrices summarized via a 360-region cortical parcellation.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Discussion<\/jats:title>\n                    <jats:p>We show that Log\u2013Euclidean and Log\u2013Cholesky estimators yield statistically and neurophysiologically indistinguishable population-level results, while the Log\u2013Cholesky formulation substantially reduces computational cost. These results establish ILPR\u2013Log\u2013Cholesky as a statistically sound and computationally efficient intrinsic regression method for large-scale, high-dimensional manifold-valued data.<\/jats:p>\n                  <\/jats:sec>","DOI":"10.3389\/fams.2025.1690557","type":"journal-article","created":{"date-parts":[[2026,1,22]],"date-time":"2026-01-22T10:11:59Z","timestamp":1769076719000},"update-policy":"https:\/\/doi.org\/10.3389\/crossmark-policy","source":"Crossref","is-referenced-by-count":2,"title":["Multivariate intrinsic local polynomial regression on isometric Riemannian manifolds. Applications to symmetric positive data"],"prefix":"10.3389","volume":"11","author":[{"given":"Ronald Garc\u00eda","family":"Reyes","sequence":"first","affiliation":[{"name":"Clinical Hospital of Chengdu Brain Science Institute, University of Electronic Science and Technology of China","place":["Chengdu, China"]},{"name":"Neuroinformatics","place":["Cuban Neurosciences Center, Cuba"]},{"name":"China-Cuba Belt and Road Joint Laboratory on Neurotechnology and Brain\u2013Apparatus Communication, University of Electronic Science and Technology of China","place":["Chengdu, China"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Ying","family":"Wang","sequence":"additional","affiliation":[{"name":"Clinical Hospital of Chengdu Brain Science Institute, University of Electronic Science and Technology of China","place":["Chengdu, China"]},{"name":"China-Cuba Belt and Road Joint Laboratory on Neurotechnology and Brain\u2013Apparatus Communication, University of Electronic Science and Technology of China","place":["Chengdu, China"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Min","family":"Li","sequence":"additional","affiliation":[{"name":"Clinical Hospital of Chengdu Brain Science Institute, University of Electronic Science and Technology of China","place":["Chengdu, China"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Giani Ega\u00f1a","family":"Fern\u00e1ndez","sequence":"additional","affiliation":[{"name":"Faculty of Mathematics and Computer Sciences, University of Havana","place":["Havana, Cuba"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Ludovico","family":"Minati","sequence":"additional","affiliation":[{"name":"Clinical Hospital of Chengdu Brain Science Institute, University of Electronic Science and Technology of China","place":["Chengdu, China"]},{"name":"Center for Mind\/Brain Science (CIMeC), University of Trento","place":["Trento, Italy"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Pedro Antonio","family":"Valdes Sosa","sequence":"additional","affiliation":[{"name":"Clinical Hospital of Chengdu Brain Science Institute, University of Electronic Science and Technology of China","place":["Chengdu, China"]},{"name":"Neuroinformatics","place":["Cuban Neurosciences Center, Cuba"]},{"name":"China-Cuba Belt and Road Joint Laboratory on Neurotechnology and Brain\u2013Apparatus Communication, University of Electronic Science and Technology of China","place":["Chengdu, China"]}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"1965","published-online":{"date-parts":[[2026,1,22]]},"reference":[{"key":"B1","doi-asserted-by":"crossref","DOI":"10.1201\/9780203748725","volume-title":"Local Polynomial Modelling and Its Applications","author":"Fan","year":"2018"},{"key":"B2","doi-asserted-by":"publisher","first-page":"1346","DOI":"10.1214\/aos\/1176325632","article-title":"Multivariate locally weighted least squares regression","volume":"22","author":"Ruppert","year":"1994","journal-title":"Ann Stat"},{"key":"B3","doi-asserted-by":"crossref","DOI":"10.1201\/b14876","volume-title":"Kernel Smoothing","author":"Wand","year":"1994"},{"key":"B4","article-title":"New directions in shape analysis","volume-title":"The Art of Statistical Science, Vol. 115","author":"Kent","year":"1992"},{"key":"B5","doi-asserted-by":"publisher","first-page":"119190","DOI":"10.1016\/j.neuroimage.2022.119190","article-title":"Harmonized-multinational qEEG norms (HarMNqEEG)","volume":"256","author":"Li","year":"2022","journal-title":"Neuroimage"},{"key":"B6","doi-asserted-by":"publisher","first-page":"e26698","DOI":"10.1002\/hbm.26698","article-title":"EEG functional connectivity as a Riemannian mediator: an application to malnutrition and cognition","volume":"45","author":"Lopez Naranjo","year":"2024","journal-title":"Hum Brain Mapp"},{"key":"B7","doi-asserted-by":"crossref","first-page":"589","DOI":"10.1016\/B978-0-12-820480-1.00184-4","article-title":"Mapping effective connectivity","volume-title":"Encyclopediaof the Human Brain","author":"Reyes","year":"2025"},{"key":"B8","doi-asserted-by":"publisher","DOI":"10.1101\/2025.02.21.639413","article-title":"Lifespan mapping of EEG source spectral dynamics with \u03be- \u03b1 NET","author":"Reyes","year":"2025","journal-title":"bioRxiv"},{"key":"B9","doi-asserted-by":"publisher","DOI":"10.48550\/arXiv.2503.02553","article-title":"Exponential Speedup of the janashia-lagvilava matrix spectral factorization algorithm","author":"Wang","year":"2025","journal-title":"arXiv [preprint]"},{"key":"B10","doi-asserted-by":"publisher","first-page":"231","DOI":"10.1016\/j.neuroimage.2013.02.008","article-title":"A note on the phase locking value and its properties","volume":"74","author":"Aydore","year":"2013","journal-title":"Neuroimage"},{"key":"B11","doi-asserted-by":"publisher","first-page":"1713","DOI":"10.1109\/TPAMI.2008.75","article-title":"Pedestrian detection via classification on Riemannian manifolds","volume":"30","author":"Tuzel","year":"2008","journal-title":"IEEE Trans Pattern Anal Mach Intell"},{"key":"B12","doi-asserted-by":"publisher","first-page":"41","DOI":"10.1007\/s11263-005-3222-z","article-title":"A Riemannian framework for tensor computing","volume":"66","author":"Pennec","year":"2006","journal-title":"Int J Comput Vis"},{"key":"B13","doi-asserted-by":"publisher","first-page":"3334","DOI":"10.1214\/17-AOS1660","article-title":"Principal component analysis for functional data on Riemannian manifolds and spheres","volume":"46","author":"Dai","year":"2018","journal-title":"Ann Stat"},{"key":"B14","doi-asserted-by":"publisher","first-page":"1102","DOI":"10.1214\/09-AOAS249","article-title":"Non-Euclidean statistics for covariance matrices, with applications to diffusion tensor imaging","volume":"3","author":"Dryden","year":"2009","journal-title":"Ann Appl Stat"},{"key":"B15","doi-asserted-by":"publisher","DOI":"10.48550\/arXiv.2101.11872","article-title":"Robust extrinsic regression analysis for manifold valued data","author":"Lee","year":"2021","journal-title":"arXiv [preprint]"},{"key":"B16","doi-asserted-by":"publisher","first-page":"887","DOI":"10.1214\/18-BA1135","article-title":"Extrinsic Gaussian processes for regression and classification on manifolds","volume":"14","author":"Lin","year":"2019","journal-title":"Bayesian Anal"},{"key":"B17","doi-asserted-by":"publisher","first-page":"1261","DOI":"10.1080\/01621459.2016.1208615","article-title":"Extrinsic local regression on manifold-valued data","volume":"112","author":"Lin","year":"2017","journal-title":"J Am Stat Assoc"},{"key":"B18","doi-asserted-by":"publisher","first-page":"3533","DOI":"10.1214\/18-AOS1787","article-title":"Intrinsic Riemannian functional data analysis","volume":"47","author":"Lin","year":"2019","journal-title":"Ann Stat"},{"key":"B19","doi-asserted-by":"publisher","first-page":"691","DOI":"10.1214\/17-AOS1624","article-title":"Fr\u00e9chet regression for random objects with Euclidean predictors","volume":"47","author":"Petersen","year":"2019","journal-title":"Ann Stat"},{"key":"B20","doi-asserted-by":"publisher","first-page":"117464","DOI":"10.1016\/j.neuroimage.2020.117464","article-title":"Re-visiting Riemannian geometry of symmetric positive definite matrices for the analysis of functional connectivity","volume":"225","author":"You","year":"2021","journal-title":"Neuroimage"},{"key":"B21","doi-asserted-by":"publisher","first-page":"697","DOI":"10.1111\/j.1467-9868.2011.01022.x","article-title":"Local polynomial regression for symmetric positive definite matrices","volume":"74","author":"Yuan","year":"2012","journal-title":"J R Stat Soc B (Stat Methodol)"},{"key":"B22","first-page":"115","article-title":"New directions in shape analysis","volume-title":"The Art of Statistical Science","author":"Kent","year":"1992"},{"key":"B23","doi-asserted-by":"publisher","first-page":"291","DOI":"10.1080\/10618600.2017.1390468","article-title":"A new unified approach for the simulation of a wide class of directional distributions","volume":"27","author":"Kent","year":"2018","journal-title":"J Comput Graph Stat"},{"key":"B24","doi-asserted-by":"publisher","first-page":"2705","DOI":"10.1109\/CVPR.2014.352","article-title":"Multivariate General Linear Models (MGLM) on riemannian manifolds with applications to statistical analysis of diffusion weighted images","author":"Kim","year":"2014","journal-title":"Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition"},{"key":"B25","doi-asserted-by":"publisher","first-page":"361","DOI":"10.1093\/biomet\/asac055","article-title":"Additive models for symmetric positive-definite matrices and Lie groups","volume":"110","author":"Lin","year":"2022","journal-title":"Biometrika"},{"key":"B26","doi-asserted-by":"crossref","DOI":"10.1007\/978-3-319-91755-9","volume-title":"Introduction to Riemannian Manifolds, Vol. 176","author":"Lee","year":"2018"},{"key":"B27","doi-asserted-by":"publisher","first-page":"411","DOI":"10.1002\/mrm.20965","article-title":"Log-Euclidean metrics for fast and simple calculus on diffusion tensors","volume":"56","author":"Arsigny","year":"2006","journal-title":"Magn Reson Med"},{"key":"B28","doi-asserted-by":"crossref","first-page":"629","DOI":"10.1007\/978-3-642-15995-4_78","article-title":"Riemannian geometry applied to BCI classification","volume-title":"International Conference on Latent Variable Analysis and Signal Separation","author":"Barachant","year":"2010"},{"key":"B29","doi-asserted-by":"publisher","first-page":"1353","DOI":"10.1137\/18M1221084","article-title":"Riemannian geometry of symmetric positive definite matrices via Cholesky decomposition","volume":"40","author":"Lin","year":"2019","journal-title":"SIAM J Matrix Anal Appl"},{"key":"B30","first-page":"721","article-title":"Assessing LBICO, a bicoherence-derived multi-frequency index measured from EEG signals","volume":"16","author":"C\u00e1ceres","year":"2024","journal-title":"Rev Cubana Inf M\u00e9d"},{"key":"B31","doi-asserted-by":"publisher","first-page":"359","DOI":"10.1016\/0304-4076(91)90107-O","article-title":"Asymptotic optimality of generalized CL, cross-validation, and generalized cross-validation in regression with heteroskedastic errors","volume":"47","author":"Andrews","year":"1991","journal-title":"J Econom"},{"key":"B32","doi-asserted-by":"publisher","DOI":"10.48550\/arXiv.2111.02990","article-title":"The geometry of mixed-Euclidean metrics on symmetric positive definite matrices","author":"Thanwerdas","year":"2021","journal-title":"arXiv [preprint]"},{"key":"B33","doi-asserted-by":"publisher","first-page":"20","DOI":"10.2307\/1969989","article-title":"The imbedding problem for Riemannian manifolds","volume":"63","author":"Nash","year":"1956","journal-title":"Ann Math"},{"key":"B34","doi-asserted-by":"publisher","first-page":"1851","DOI":"10.1137\/22M1471729","article-title":"Theoretically and computationally convenient geometries on full-rank correlation matrices","volume":"43","author":"Thanwerdas","year":"2022","journal-title":"SIAM J Matrix Anal Appl"},{"key":"B35","doi-asserted-by":"publisher","first-page":"85","DOI":"10.1016\/S0377-0427(00)00393-9","article-title":"The ubiquitous Kronecker product","volume":"123","author":"Van Loan","year":"2000","journal-title":"J Comput Appl Math"},{"key":"B36","doi-asserted-by":"publisher","first-page":"165","DOI":"10.1016\/0024-3795(91)90332-Q","article-title":"Block Kronecker products and the VECB operator","volume":"149","author":"Koning","year":"1991","journal-title":"Linear Algebra Appl"},{"key":"B37","doi-asserted-by":"publisher","DOI":"10.48550\/arXiv.2011.01833","article-title":"Higher order differential analysis with vectorized derivatives","author":"Chac\u00f3n","year":"2020","journal-title":"arXiv [preprint]"},{"key":"B38","doi-asserted-by":"crossref","first-page":"337","DOI":"10.15625\/0866-7136\/14196","article-title":"Taylor expansion for matrix functions of vector variable using the Kronecker product","volume":"41","author":"Van Khang","year":"2019","journal-title":"Vietnam J Mech"},{"key":"B39","doi-asserted-by":"publisher","first-page":"474","DOI":"10.1016\/0022-2496(85)90006-9","article-title":"Matrix differential calculus with applications to simple, Hadamard, and Kronecker products","volume":"29","author":"Magnus","year":"1985","journal-title":"J Math Psychol"},{"key":"B40","doi-asserted-by":"publisher","first-page":"2200","DOI":"10.1016\/j.jmva.2010.05.005","article-title":"On the concept of matrix derivative","volume":"101","author":"Magnus","year":"2010","journal-title":"J Multivar Anal"},{"key":"B41","doi-asserted-by":"publisher","first-page":"981","DOI":"10.1016\/j.mechmachtheory.2010.03.005","article-title":"Consistent definition of partial derivatives of matrix functions in dynamics of mechanical systems","volume":"45","author":"Van Khang","year":"2010","journal-title":"Mech Mach Theory"},{"key":"B42","doi-asserted-by":"publisher","first-page":"172","DOI":"10.1016\/j.neucom.2012.12.039","article-title":"Classification of covariance matrices using a Riemannian-based kernel for BCI applications","volume":"112","author":"Barachant","year":"2013","journal-title":"Neurocomputing"},{"key":"B43","unstructured":"Reyes\n              RG\n            \n          \n          Multivariate Intrinsic Local Polynomial Regression on Isometric Riemannian Manifolds\n          \n          2023"},{"key":"B44","first-page":"1996","article-title":"Riemannian normality and random rotation-invariant statistics for high angular resolution diffusion imaging","volume":"34","author":"Schwartzman","year":"2006","journal-title":"Ann Stat"},{"key":"B45","doi-asserted-by":"publisher","first-page":"1085","DOI":"10.1198\/016214507000000581","article-title":"Statistical analysis of diffusion tensors in diffusion-weighted magnetic resonance image data (with discussion)","volume":"102","author":"Zhu","year":"2007","journal-title":"J Am Statist Ass"},{"key":"B46","doi-asserted-by":"publisher","first-page":"53","DOI":"10.1016\/S0720-048X(02)00328-5","article-title":"diffusion tensor imaging: recent advance and new techniques for diffusion tensor visualization","volume":"46","author":"Masutani","year":"2003","journal-title":"Eur J Radiol"},{"key":"B47","doi-asserted-by":"publisher","first-page":"93","DOI":"10.1016\/S1361-8415(02)00053-1","article-title":"Processing and visualization for diffusion tensor MRI","volume":"6","author":"Westin","year":"2002","journal-title":"Med Image Anal"},{"key":"B48","doi-asserted-by":"publisher","DOI":"10.48550\/arXiv.2203.09253","article-title":"Visualizing Riemannian data with Rie-SNE","author":"Bergsson","year":"2022","journal-title":"arXiv [preprint]"},{"key":"B49","doi-asserted-by":"publisher","first-page":"995","DOI":"10.1109\/TMI.2004.831793","article-title":"Principal geodesic analysis for the study of nonlinear statistics of shape","volume":"23","author":"Fletcher","year":"2004","journal-title":"IEEE Trans Med Imaging"},{"key":"B50","doi-asserted-by":"publisher","first-page":"250","DOI":"10.1016\/j.sigpro.2005.12.018","article-title":"Riemannian geometry for the statistical analysis of diffusion tensor data","volume":"87","author":"Fletcher","year":"2007","journal-title":"Signal Process"},{"key":"B51","doi-asserted-by":"publisher","first-page":"278","DOI":"10.1002\/cmr.a.20094","article-title":"Physical foundations, models, and methods of diffusion magnetic resonance imaging of the brain: a review","volume":"30","author":"Minati","year":"2007","journal-title":"Concepts Magn Reson A"},{"key":"B52","doi-asserted-by":"publisher","first-page":"171","DOI":"10.1038\/nature18933","article-title":"A multi-modal parcellation of human cerebral cortex","volume":"536","author":"Glasser","year":"2016","journal-title":"Nature"},{"key":"B53","doi-asserted-by":"publisher","first-page":"2787","DOI":"10.1090\/proc\/12953","article-title":"Positive definite matrices and the S-divergence","volume":"144","author":"Sra","year":"2016","journal-title":"Proc Am Math Soc"},{"key":"B54","doi-asserted-by":"crossref","first-page":"2052","DOI":"10.1109\/TCST.2012.2233476","article-title":"Direct causality detection via the transfer entropy approach","volume":"21","author":"Duan","year":"2013","journal-title":"IEEE Trans Control Syst Technol"},{"key":"B55","first-page":"2022","article-title":"Modeling neural activity in neurodegenerative diseases through a neural field model with variable density of neurons","author":"Reyes","year":"2022","journal-title":"bioRxiv"},{"key":"B56","doi-asserted-by":"crossref","first-page":"119299","DOI":"10.1016\/j.neuroimage.2022.119299","article-title":"A Riemannian approach to predicting brain function from the structural connectome","volume":"257","author":"Benkarim","year":"2022","journal-title":"NeuroImage"},{"key":"B57","first-page":"826","article-title":"Mean template for tensor-based morphometry using deformation tensors","volume-title":"International Conference on Medical Image Computing and Computer-Assisted Intervention","author":"Lepor\u00e9","year":"2007"}],"container-title":["Frontiers in Applied Mathematics and Statistics"],"original-title":[],"link":[{"URL":"https:\/\/www.frontiersin.org\/articles\/10.3389\/fams.2025.1690557\/full","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,1,22]],"date-time":"2026-01-22T10:12:05Z","timestamp":1769076725000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.frontiersin.org\/articles\/10.3389\/fams.2025.1690557\/full"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2026,1,22]]},"references-count":57,"alternative-id":["10.3389\/fams.2025.1690557"],"URL":"https:\/\/doi.org\/10.3389\/fams.2025.1690557","relation":{},"ISSN":["2297-4687"],"issn-type":[{"value":"2297-4687","type":"electronic"}],"subject":[],"published":{"date-parts":[[2026,1,22]]},"article-number":"1690557"}}