Spatial StatisticsPub Date : 2025-06-01Epub Date: 2025-04-23DOI: 10.1016/j.spasta.2025.100901
Duncan Lee
{"title":"Computationally efficient spatio-temporal disease mapping for big data","authors":"Duncan Lee","doi":"10.1016/j.spasta.2025.100901","DOIUrl":"10.1016/j.spasta.2025.100901","url":null,"abstract":"<div><div>Disease mapping models estimate the spatio-temporal variation in population-level disease risks or rates across a set of <span><math><mi>K</mi></math></span> areal units for <span><math><mi>N</mi></math></span> time periods, aiming to identify temporal trends and spatial hotspots. Highly parameterised Bayesian hierarchical models with over <span><math><mrow><mi>K</mi><mi>N</mi></mrow></math></span> random effects are commonly used to estimate this spatio-temporal variation, which are assigned autoregressive and conditional autoregressive prior distributions. These models work well when there are tens of thousands of data points, but are likely to be computationally burdensome when this rises to hundreds of thousands or above. This paper proposes a computationally efficient alternative, which can fit a range of spatio-temporal disease trends almost as well as existing highly parameterised models but only takes around 5% to 40% of the time to implement. It achieves this by modelling the average spatial and temporal trends in the data with autoregressive type random effects, which are augmented by an observation-driven process using functions of earlier data as additional covariates in the model. The efficacy of this methodology is tested by simulation, before being applied to the motivating study that estimates the spatio-temporal trends in asthma, cancer, coronary heart and chronic obstructive pulmonary disease prevalences for <span><math><mrow><mi>K</mi><mo>=</mo><mn>32</mn><mo>,</mo><mn>751</mn></mrow></math></span> small areas over <span><math><mrow><mi>N</mi><mo>=</mo><mn>13</mn></mrow></math></span> years in England.</div></div>","PeriodicalId":48771,"journal":{"name":"Spatial Statistics","volume":"67 ","pages":"Article 100901"},"PeriodicalIF":2.1,"publicationDate":"2025-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143873638","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Spatial StatisticsPub Date : 2025-06-01Epub Date: 2025-04-24DOI: 10.1016/j.spasta.2025.100904
Nicolas Coloma, William Kleiber
{"title":"Random elastic space–time (REST) prediction","authors":"Nicolas Coloma, William Kleiber","doi":"10.1016/j.spasta.2025.100904","DOIUrl":"10.1016/j.spasta.2025.100904","url":null,"abstract":"<div><div>Statistical modeling and interpolation of space–time processes has gained increasing relevance over the last few years. However, real world data often exhibit characteristics that challenge conventional methods such as nonstationarity and temporal misalignment. For example, high frequency solar irradiance data are typically observed at fine temporal scales, but at sparse spatial sampling, so space–time interpolation is necessary to support solar energy studies. The nonstationarity and phase misalignment of such data challenges extant approaches. We propose random elastic space–time (REST) prediction, a novel method that addresses temporally-varying phase misalignment by combining elastic alignment and conventional kriging techniques. Moreover, uncertainty in both amplitude and phase alignment can be readily quantified in a conditional simulation framework, whereas conventional space–time methods only address amplitude uncertainty. We illustrate our approach on a challenging solar irradiance dataset, where our method demonstrates superior predictive distributions compared to existing geostatistical and functional data analytic techniques.</div></div>","PeriodicalId":48771,"journal":{"name":"Spatial Statistics","volume":"67 ","pages":"Article 100904"},"PeriodicalIF":2.1,"publicationDate":"2025-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143886951","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Spatial StatisticsPub Date : 2025-06-01Epub Date: 2025-03-27DOI: 10.1016/j.spasta.2025.100899
Murat Okatan
{"title":"Fast computation of the statistical significance test for spatio-temporal receptive field estimates obtained using spike-triggered averaging of binary pseudo-random sequences","authors":"Murat Okatan","doi":"10.1016/j.spasta.2025.100899","DOIUrl":"10.1016/j.spasta.2025.100899","url":null,"abstract":"<div><h3>Background</h3><div>Spatio-temporal receptive fields (STRFs) of visual neurons are often estimated using spike-triggered averaging (STA) with binary pseudo-random stimulus sequences. An exact analytical test—called the STA-BPRS test—has been developed to determine the statistical significance of each pixel in the STRF estimate. However, computing this test can take minutes to days, or even longer, for certain neurons.</div></div><div><h3>New method</h3><div>Here, the STA-BPRS test is accelerated by approximating the null distribution of STRF pixel estimates with a Normal distribution. This methodological refinement significantly reduces computation time, making large-scale data analysis feasible.</div></div><div><h3>Results</h3><div>The approximate test is systematically validated on real mouse retinal ganglion cell data and synthetic spike train data, demonstrating that it yields identical significance thresholds to the exact test. For neurons where, exact computation would be prohibitively long (e.g., hundreds of years), the approximate test completes in seconds or minutes.</div></div><div><h3>Comparison with existing methods</h3><div>Few approaches address pixel-by-pixel significance in STA-based STRF estimates. While subspace methods like spike-triggered covariance exist for STRF estimation, they typically do not provide direct voxel-wise or pixel-wise p-values. The proposed method specifically accelerates an exact distribution-based test.</div></div><div><h3>Conclusions and impact</h3><div>The proposed Normal approximation drastically reduces computation time, enabling high-throughput analysis of STRF mapping from spike data. This advancement may foster broader adoption of precise statistical tests of STRFs in large-scale, high-density electrophysiological recordings. Moreover, fast detection of significant STRF features could facilitate closed-loop experiments where stimuli dynamically adapt to changing STRF structures.</div></div>","PeriodicalId":48771,"journal":{"name":"Spatial Statistics","volume":"67 ","pages":"Article 100899"},"PeriodicalIF":2.1,"publicationDate":"2025-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143767964","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Spatial StatisticsPub Date : 2025-06-01Epub Date: 2025-02-17DOI: 10.1016/j.spasta.2025.100886
Antonella Congedi, Sandra De Iaco, Donato Posa
{"title":"A term structure geostatistical model with correlated residuals: A comparative analysis","authors":"Antonella Congedi, Sandra De Iaco, Donato Posa","doi":"10.1016/j.spasta.2025.100886","DOIUrl":"10.1016/j.spasta.2025.100886","url":null,"abstract":"<div><div>The growth of financial markets and the emerging derivative instruments require the development of advanced techniques for forecasting the term structure of interest rates. In this context, two significant dimensions, i.e. maturity and time, need to be jointly considered in the modeling procedure. In the literature, the Nelson–Siegel model is commonly used to explain the dependence of the interest rates on maturity and time. However, it cannot be excluded that the residuals obtained from Nelson–Siegel estimates are still correlated. At this purpose, a geostatistical approach is adopted and an innovative modeling solution is provided. Indeed, differently from the existing contributions, this paper proposes a dynamic model for predicting the term structure of spot interest rates, where the joint evolution with respect to time and maturity is considered for both the deterministic and the stochastic parts of the model. The relevance as well as the potentiality of the geostatistical modeling techniques extended to treat observations not strictly referred to a geographic system, has been properly underlined. For comparative reasons, different hypotheses on the random field, utilized to describe the interest rates and its trend component, are also assumed and a comparison among predictive performance of alternative models is discussed.</div></div>","PeriodicalId":48771,"journal":{"name":"Spatial Statistics","volume":"67 ","pages":"Article 100886"},"PeriodicalIF":2.1,"publicationDate":"2025-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143519048","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Spatial StatisticsPub Date : 2025-06-01Epub Date: 2025-04-05DOI: 10.1016/j.spasta.2025.100898
Jakub J. Pypkowski, Adam M. Sykulski, James S. Martin
{"title":"Isotropy testing in spatial point patterns: nonparametric versus parametric replication under misspecification","authors":"Jakub J. Pypkowski, Adam M. Sykulski, James S. Martin","doi":"10.1016/j.spasta.2025.100898","DOIUrl":"10.1016/j.spasta.2025.100898","url":null,"abstract":"<div><div>Several hypothesis testing methods have been proposed to validate the assumption of isotropy in spatial point patterns. A majority of these methods are characterised by an unknown distribution of the test statistic under the null hypothesis of isotropy. Parametric approaches to approximating the distribution involve simulation of patterns from a user-specified isotropic model. Alternatively, nonparametric replicates of the test statistic under isotropy can be used to waive the need for specifying a model. In this paper, we first present a general framework which allows for the integration of a selected nonparametric replication method into isotropy testing. We then conduct a large simulation study comprising application-like scenarios to assess the performance of tests with different parametric and nonparametric replication methods. In particular, we explore distortions in test size and power caused by model misspecification, and demonstrate the advantages of nonparametric replication in such scenarios.</div></div>","PeriodicalId":48771,"journal":{"name":"Spatial Statistics","volume":"67 ","pages":"Article 100898"},"PeriodicalIF":2.1,"publicationDate":"2025-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143799115","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Spatial StatisticsPub Date : 2025-06-01Epub Date: 2025-04-28DOI: 10.1016/j.spasta.2025.100903
Gianfranco Piras , Mauricio Sarrias
{"title":"A J-test for spatial autoregressive binary models","authors":"Gianfranco Piras , Mauricio Sarrias","doi":"10.1016/j.spasta.2025.100903","DOIUrl":"10.1016/j.spasta.2025.100903","url":null,"abstract":"<div><div>Spatial autoregressive binary models are well established in spatial statistics and econometric literature. Recently, different estimation methods have been proposed that account for logistic as well as probit regressions. In spatial models the choice of the spatial weighting matrix is crucial to reflect the amount of correlation in the data. This article proposes a simple <span><math><mi>J</mi></math></span>-test procedure for spatial autoregressive binary model. Since the <span><math><mi>J</mi></math></span>-test is a non-nested test, it can be used, among other things, to test the specification of the spatial weighting matrix. The <span><math><mi>J</mi></math></span>-test is based on augmenting the null model with the predictor from the alternative model(s). After defining these predictors, we develop the theory and derive the steps for the <span><math><mi>J</mi></math></span>-test. We also evaluate the finite sample properties in the context of a Monte Carlo experiment. An empirical application on firms’ decisions to reopen in the aftermath of Hurricane Katrina for New Orleans is also presented.</div></div>","PeriodicalId":48771,"journal":{"name":"Spatial Statistics","volume":"67 ","pages":"Article 100903"},"PeriodicalIF":2.1,"publicationDate":"2025-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143894998","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Spatial StatisticsPub Date : 2025-06-01Epub Date: 2025-02-27DOI: 10.1016/j.spasta.2025.100890
Guilherme Ludwig , Yuan Wang , Tingjin Chu , Haonan Wang , Jun Zhu
{"title":"Bayesian analysis and variable selection for spatial count data with an application to Rio de Janeiro gun violence","authors":"Guilherme Ludwig , Yuan Wang , Tingjin Chu , Haonan Wang , Jun Zhu","doi":"10.1016/j.spasta.2025.100890","DOIUrl":"10.1016/j.spasta.2025.100890","url":null,"abstract":"<div><div>Statistical analysis has been successfully applied to crime data for identification of crime hot spots and prediction of future crimes. In this paper, our main objective is to identify key factors for gun violence in Rio de Janeiro and study the relationship between these key factors and the number of reported events. We use a Bayesian hierarchical stochastic Poisson regression model for spatial counts, which enables us to address the over-dispersed count data and to handle the spatial correlation. Moreover, we propose a variable selection method for key factor identification based on the spike-and-slab prior distribution for the regression coefficients. A new Gibbs sampler is developed for sampling from the posterior distributions with the help of augmentation of Pólya-Gamma auxiliary variables. Simulation studies are used to demonstrate the performance of our proposed approach. Our analysis of the gun violence data in Rio de Janeiro reveals the relationship between violence events and socio-demographic covariates as well as an interpretable spatial random effect that accounts for unmeasured covariate information.</div></div>","PeriodicalId":48771,"journal":{"name":"Spatial Statistics","volume":"67 ","pages":"Article 100890"},"PeriodicalIF":2.1,"publicationDate":"2025-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143519049","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Spatial StatisticsPub Date : 2025-06-01Epub Date: 2025-03-12DOI: 10.1016/j.spasta.2025.100892
Miao Long, Zhimeng Sun
{"title":"Bayesian adaptive Lasso estimation for partially linear hierarchical spatial autoregressive model","authors":"Miao Long, Zhimeng Sun","doi":"10.1016/j.spasta.2025.100892","DOIUrl":"10.1016/j.spasta.2025.100892","url":null,"abstract":"<div><div>This paper presents a Bayesian adaptive Lasso estimation approach for partially linear hierarchical spatial autoregressive models. Despite advancements in spatial modeling, two key gaps remain: the lack of non-linear components in hierarchical spatial autoregressive models to capture complex spatial relationships, and the insufficient application of dimensionality reduction techniques to address high-dimensionality and overfitting. This paper addresses these issues by combining partially linear models with spatial autoregressive structures and incorporating dimensionality reduction techniques to enhance model efficiency and mitigate overfitting. The hierarchical structure facilitates multi-level modeling, accommodating complex data relationships. The Bayesian adaptive Lasso technique ensures effective variable selection and regularization, improving model interpretability and performance. Simulations and real data applications demonstrate the proposed method’s excellent performance. This work offers valuable insights for researchers and practitioners in dealing with spatially correlated data in various fields.</div></div>","PeriodicalId":48771,"journal":{"name":"Spatial Statistics","volume":"67 ","pages":"Article 100892"},"PeriodicalIF":2.1,"publicationDate":"2025-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143627809","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Similarity and geographically weighted regression considering spatial scales of features space","authors":"Shifeng Yu , Xiaoyu Hu , Yehua Sheng , Chenmeng Zhao","doi":"10.1016/j.spasta.2025.100897","DOIUrl":"10.1016/j.spasta.2025.100897","url":null,"abstract":"<div><div>Unlike a geographically weighted regression (GWR), a similarity and geographically weighted regression (SGWR) calculates weights using attribute and geographic similarities, thereby effectively improving the accuracy of the model. Nevertheless, SGWR does not set an attribute similarity bandwidth. This leads to its inability to measure the scale of variation of spatial processes in feature space. In addition, owing to the solving method used, SGWR can get stuck in local optima. To address these issues, this study proposed an improved similarity and geographically weighted regression model (SGWR-GD) that adds bandwidth to the attribute similarity kernel function. This parameter gives SGWR-GD the ability to measure the scale of change of spatial processes in the attribute dimension and thus enhances the flexibility of modelling. When solving the model, SGWR-GD first calculated the gradient of the model's Modified Akaike Information Criterion (AICc) with respect to the two bandwidths and the impact ratio. Subsequently, the optimal global solution of the model was obtained based on a gradient descent algorithm with box constraints. SGWR-GD and SGWR were applied to five different datasets and the accuracies of their fitting results were compared. SGWR-GD significantly improved the accuracy of the model compared to SGWR. In addition, the SGWR-GD stably determined the global optimal solution for each parameter. Simultaneously, the distribution of local residuals was also more stable.</div></div>","PeriodicalId":48771,"journal":{"name":"Spatial Statistics","volume":"67 ","pages":"Article 100897"},"PeriodicalIF":2.1,"publicationDate":"2025-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143760349","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Spatial StatisticsPub Date : 2025-06-01Epub Date: 2025-04-05DOI: 10.1016/j.spasta.2025.100896
Pilar García-Soidán , Tomás R. Cotos-Yáñez
{"title":"Nonparametric approaches for direct approximation of the spatial quantiles","authors":"Pilar García-Soidán , Tomás R. Cotos-Yáñez","doi":"10.1016/j.spasta.2025.100896","DOIUrl":"10.1016/j.spasta.2025.100896","url":null,"abstract":"<div><div>The estimation of the spatial quantiles provides information on the thresholds of a spatial variable. This methodology is particularly appealing for its application to data of pollutants, so as to assess their level of risk. A spatial quantile can be approximated through different mechanisms, proposed in the statistics literature, although these approaches suffer from several drawbacks, regarding their lack of optimality or the fact of not leading to direct approximations. Thus, the current work introduces alternative procedures, which try to overcome the aforementioned issues by employing order statistics, similarly as done for independent data. With this aim, the available observations are appropriately transformed to yield a sample of the process at each target site, so that the data obtained are then ordered and used to derive the spatial quantile at the corresponding location. The new methodology can be directly applied to data from processes that are either stationary or that deviate from this condition for a non-constant trend and, additionally, it can be even extended to heteroscedastic data. Simulation studies under different scenarios have been accomplished, whose results show the adequate performance of the proposed estimators. A further step of this research is the application of the new approaches to data of nitrogen dioxide concentrations, to exemplify the potential of the quantile estimates to check the thresholds of a pollutant at a specific moment, as well as their evolution over time.</div></div>","PeriodicalId":48771,"journal":{"name":"Spatial Statistics","volume":"67 ","pages":"Article 100896"},"PeriodicalIF":2.1,"publicationDate":"2025-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143783992","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}