Spatial StatisticsPub Date : 2024-03-05DOI: 10.1016/j.spasta.2024.100818
Carolina Euán , Ying Sun , Brian J. Reich
{"title":"Regime-based precipitation modeling: A spatio-temporal approach","authors":"Carolina Euán , Ying Sun , Brian J. Reich","doi":"10.1016/j.spasta.2024.100818","DOIUrl":"10.1016/j.spasta.2024.100818","url":null,"abstract":"<div><p>In this paper, we propose a new regime-based model to describe spatio-temporal dynamics of precipitation data. Precipitation is one of the most essential factors for multiple human-related activities such as agriculture production. Therefore, a detailed and accurate understanding of the rain for a given region is needed. Motivated by the different formations of precipitation systems (convective, frontal, and orographic), we proposed a hierarchical regime-based spatio-temporal model for precipitation data. We use information about the values of neighboring sites to identify such regimes, allowing spatial and temporal dependence to be different among regimes. Using the Bayesian approach with R INLA, we fit our model to the Guanajuato state (Mexico) precipitation data case study to understand the spatial and temporal dependencies of precipitation in this region. Our findings show the regime-based model’s versatility and compare it with the truncated Gaussian model.</p></div>","PeriodicalId":48771,"journal":{"name":"Spatial Statistics","volume":"60 ","pages":"Article 100818"},"PeriodicalIF":2.3,"publicationDate":"2024-03-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S2211675324000095/pdfft?md5=34516482aa33a4d0c7231ce4614fe6c6&pid=1-s2.0-S2211675324000095-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140074621","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 : 2024-03-01DOI: 10.1016/j.spasta.2024.100820
Mohammad Moradi , Jennifer Brown
{"title":"Mapping using an adaptive sampling design","authors":"Mohammad Moradi , Jennifer Brown","doi":"10.1016/j.spasta.2024.100820","DOIUrl":"10.1016/j.spasta.2024.100820","url":null,"abstract":"<div><p>Interpolation is commonly used in the construction of maps and images when there is limited information for some of the sites. The accuracy of interpolation methods depends, in part, on the location of the sample sites where more complete information has been gathered. An initial survey design where the sample sites are spaced so there is wide-spread coverage is desirable. However, when there is considerable variation in the variable of interest, other design features may be preferable. Here we introduce an adaptive design where in the first stage of site selection gives wide-spread coverage, and in subsequent stages additional sites are selected adjacent to areas of high variability.</p></div>","PeriodicalId":48771,"journal":{"name":"Spatial Statistics","volume":"60 ","pages":"Article 100820"},"PeriodicalIF":2.3,"publicationDate":"2024-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140074924","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 : 2024-02-05DOI: 10.1016/j.spasta.2024.100817
Mevin B. Hooten , Michael R. Schwob , Devin S. Johnson , Jacob S. Ivan
{"title":"Geostatistical capture–recapture models","authors":"Mevin B. Hooten , Michael R. Schwob , Devin S. Johnson , Jacob S. Ivan","doi":"10.1016/j.spasta.2024.100817","DOIUrl":"https://doi.org/10.1016/j.spasta.2024.100817","url":null,"abstract":"<div><p>Methods for population estimation and inference have evolved over the past decade to allow for the incorporation of spatial information when using capture–recapture study designs. Traditional approaches to specifying spatial capture–recapture (SCR) models often rely on an individual-based detection function that decays as a detection location is farther from an individual’s activity center. Traditional SCR models are intuitive because they incorporate mechanisms of animal space use based on their assumptions about activity centers. We modify the SCR model to accommodate a wide range of space use patterns, including for those individuals that may exhibit traditional elliptical utilization distributions. Our approach uses underlying Gaussian processes to characterize the space use of individuals. This allows us to account for multimodal and other complex space use patterns that may arise due to movement. We refer to this class of models as geostatistical capture–recapture (GCR) models. We adapt a recursive computing strategy to fit GCR models to data in stages, some of which can be parallelized. This technique facilitates implementation and leverages modern multicore and distributed computing environments. We demonstrate the application of GCR models by analyzing both simulated data and a data set involving capture histories of snowshoe hares in central Colorado, USA.</p></div>","PeriodicalId":48771,"journal":{"name":"Spatial Statistics","volume":"59 ","pages":"Article 100817"},"PeriodicalIF":2.3,"publicationDate":"2024-02-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S2211675324000083/pdfft?md5=09305eb130f1cfc623cdc920435039a4&pid=1-s2.0-S2211675324000083-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139699661","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 : 2024-02-05DOI: 10.1016/j.spasta.2024.100816
Qi Zhang , Alexandra M. Schmidt , Yogendra P. Chaubey
{"title":"Modeling left-censored skewed spatial processes: The case of arsenic drinking water contamination","authors":"Qi Zhang , Alexandra M. Schmidt , Yogendra P. Chaubey","doi":"10.1016/j.spasta.2024.100816","DOIUrl":"https://doi.org/10.1016/j.spasta.2024.100816","url":null,"abstract":"<div><p>Commonly, observations from environmental processes are spatially structured and present skewed distributions. Recently, different models have been proposed to model spatial processes in their original scale. This work was motivated by modeling the levels of arsenic groundwater concentration in Comilla, a district of Bangladesh. Some of the observations are left censored. We propose spatial gamma models and explore different parametrizations of the gamma distribution. The gamma model naturally accounts for the skewness present in the data and the fact that arsenic levels are positive. We compare our proposed approaches with two skewed models proposed in the literature. Inference is performed under the Bayesian paradigm and interpolation to unobserved locations of interest naturally accounts for the estimation of the parameters in the proposed model. For the arsenic dataset, one of our proposed gamma models performs best in comparison to previous spatial models for skewed data, in terms of scoring rules criteria. Moreover, under the skewed models, some of the lower limits of the 95% posterior predictive distributions provide negative values violating the assumption that observations are strictly positive. The gamma distribution provides a reasonable, and simpler, alternative to account for the skewness present in the data and provide forecasts that are within the valid values of the observations.</p></div>","PeriodicalId":48771,"journal":{"name":"Spatial Statistics","volume":"59 ","pages":"Article 100816"},"PeriodicalIF":2.3,"publicationDate":"2024-02-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S2211675324000071/pdfft?md5=4e6657431bd93da01f8723d3f2fdc303&pid=1-s2.0-S2211675324000071-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139714900","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}
{"title":"Robust interaction detector: A case of road life expectancy analysis","authors":"Zehua Zhang , Yongze Song , Lalinda Karunaratne , Peng Wu","doi":"10.1016/j.spasta.2024.100814","DOIUrl":"10.1016/j.spasta.2024.100814","url":null,"abstract":"<div><p>Spatial stratified heterogeneity, revealing the disparity mechanisms across spatial strata, can be effectively quantified using the geographical detector (GD). GD requires reasonable spatial discretization strategies to investigate the spatial association between the target variable and numerical independent variables. In previous studies, the Robust Geographical Detector (RGD) optimized spatial strata for examining the power of determinants (PD) of individual variables, which demonstrate more robust spatial discretization than other models. However, the GD's interaction detector that explores PD of the interaction of two variables still needs to be enhanced by the robust spatial discretization. This study develops a Robust Interaction Detector (RID), an improved interaction detector, using change detection algorithms for the robust spatial stratified heterogeneity analysis with multiple explanatory variables. RID is applied in a road life expectancy analysis in Western Australia. Results show that RID presents higher PD values than previous GD models, ensuring the growth of PD value with more spatial strata. The RID model indicates that the interactions between various transport variables and elevation are strongly associated with road life expectancy from the perspective of spatial patterns. The developed RID model provides significant potential for enhanced geospatial factor analysis across diverse fields.</p></div>","PeriodicalId":48771,"journal":{"name":"Spatial Statistics","volume":"59 ","pages":"Article 100814"},"PeriodicalIF":2.3,"publicationDate":"2024-01-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S2211675324000058/pdfft?md5=f61d206ff82268fb072a2711dc2fed1e&pid=1-s2.0-S2211675324000058-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139516662","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 : 2024-01-18DOI: 10.1016/j.spasta.2024.100812
Yuhan Ma , Kyuhee Shin , GyuWon Lee , Joon Jin Song
{"title":"Spatial classification in the presence of measurement error","authors":"Yuhan Ma , Kyuhee Shin , GyuWon Lee , Joon Jin Song","doi":"10.1016/j.spasta.2024.100812","DOIUrl":"10.1016/j.spasta.2024.100812","url":null,"abstract":"<div><p>In recent decades, spatial classification has received considerable attention in a wide array of disciplines. In practice, binary response variable is often subject to measurement error, misclassification. To account for the misclassified response in spatial classification, we proposed validation data-based adjustment methods that use interval validation data to rectify misclassified responses. Regression calibration and multiple imputation methods are utilized to correct the misclassified outcomes at the locations where the gold-standard device is not available. Generalized linear mixed model and indicator Kriging are applied for spatial classification at unsampled locations. Simulation studies are performed to compare the proposed methods with naive methods that ignore the misclassification. It was found that the proposed models significantly improve prediction accuracy. Additionally, the proposed models are applied for precipitation detection in South Korea.</p></div>","PeriodicalId":48771,"journal":{"name":"Spatial Statistics","volume":"59 ","pages":"Article 100812"},"PeriodicalIF":2.3,"publicationDate":"2024-01-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S2211675324000034/pdfft?md5=0a15ee5f09dc0fe93583f96e7eac46cf&pid=1-s2.0-S2211675324000034-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139501559","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 : 2024-01-17DOI: 10.1016/j.spasta.2024.100813
Kellie McClernon, Katherine Goode, Daniel Ries
{"title":"A comparison of model validation approaches for echo state networks using climate model replicates","authors":"Kellie McClernon, Katherine Goode, Daniel Ries","doi":"10.1016/j.spasta.2024.100813","DOIUrl":"10.1016/j.spasta.2024.100813","url":null,"abstract":"<div><p>As global temperatures continue to rise, climate mitigation strategies such as stratospheric aerosol injections (SAI) are increasingly discussed, but the downstream effects of these strategies are not well understood. As such, there is interest in developing statistical methods to quantify the evolution of climate variable relationships during the time period surrounding an SAI. Feature importance applied to echo state network (ESN) models has been proposed as a way to understand the effects of SAI using a data-driven model. This approach depends on the ESN fitting the data well. If not, the feature importance may place importance on features that are not representative of the underlying relationships. Typically, time series prediction models such as ESNs are assessed using out-of-sample performance metrics that divide the times series into separate training and testing sets. However, this model assessment approach is geared towards forecasting applications and not scenarios such as the motivating SAI example where the objective is using a data driven model to capture variable relationships. In this paper, we demonstrate a novel use of climate model replicates to investigate the applicability of the commonly used repeated hold-out model assessment approach for the SAI application. Simulations of an SAI are generated using a simplified climate model, and different initialization conditions are used to provide independent training and testing sets containing the same SAI event. The climate model replicates enable out-of-sample measures of model performance, which are compared to the single time series hold-out validation approach. For our case study, it is found that the repeated hold-out sample performance is comparable, but conservative, to the replicate out-of-sample performance when the training set contains enough time after the aerosol injection.</p></div>","PeriodicalId":48771,"journal":{"name":"Spatial Statistics","volume":"59 ","pages":"Article 100813"},"PeriodicalIF":2.3,"publicationDate":"2024-01-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S2211675324000046/pdfft?md5=a6ba0350eda3c86948baceceabdf144c&pid=1-s2.0-S2211675324000046-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139496081","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 : 2024-01-17DOI: 10.1016/j.spasta.2023.100799
Hiroshi Yamada
{"title":"Spatial Smoothing Using Graph Laplacian Penalized Filter","authors":"Hiroshi Yamada","doi":"10.1016/j.spasta.2023.100799","DOIUrl":"10.1016/j.spasta.2023.100799","url":null,"abstract":"<div><p>This paper considers a filter for smoothing spatial data. It can be used to smooth data on the vertices of arbitrary undirected graphs with arbitrary non-negative spatial weights. It consists of a quantity analogous to Geary’s <span><math><mi>c</mi></math></span>, which is one of the most prominent measures of spatial autocorrelation. In addition, the quantity can be represented by a matrix called the graph Laplacian in spectral graph theory. We show mathematically how spatial data becomes smoother as a parameter, called the smoothing parameter, increases from 0 and is fully smoothed as the parameter goes to infinity, except for the case where the spatial data is originally fully smoothed. We also illustrate the results numerically and apply the spatial filter to climatological/meteorological data. In addition, as supplementary investigations, we examine how the sum of squared residuals and the effective degrees of freedom vary with the smoothing parameter. Finally, we review two closely related literatures to the spatial filter. One is the intrinsic conditional autoregressive model and the other is the eigenvector spatial filter. We clarify how the spatial filter considered in this paper relates to them. We then mention future research.</p></div>","PeriodicalId":48771,"journal":{"name":"Spatial Statistics","volume":"60 ","pages":"Article 100799"},"PeriodicalIF":2.3,"publicationDate":"2024-01-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139496155","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 : 2024-01-17DOI: 10.1016/j.spasta.2024.100811
Daniela Cisneros , Jordan Richards , Ashok Dahal , Luigi Lombardo , Raphaël Huser
{"title":"Deep graphical regression for jointly moderate and extreme Australian wildfires","authors":"Daniela Cisneros , Jordan Richards , Ashok Dahal , Luigi Lombardo , Raphaël Huser","doi":"10.1016/j.spasta.2024.100811","DOIUrl":"https://doi.org/10.1016/j.spasta.2024.100811","url":null,"abstract":"<div><p>Recent wildfires in Australia have led to considerable economic loss and property destruction, and there is increasing concern that climate change may exacerbate their intensity, duration, and frequency. Hazard quantification for extreme wildfires is an important component of wildfire management, as it facilitates efficient resource distribution, adverse effect mitigation, and recovery efforts. However, although extreme wildfires are typically the most impactful, both small and moderate fires can still be devastating to local communities and ecosystems. Therefore, it is imperative to develop robust statistical methods to reliably model the full distribution of wildfire spread. We do so for a novel dataset of Australian wildfires from 1999 to 2019, and analyse monthly spread over areas approximately corresponding to Statistical Areas Level 1 and 2 (SA1/SA2) regions. Given the complex nature of wildfire ignition and spread, we exploit recent advances in statistical deep learning and extreme value theory to construct a parametric regression model using graph convolutional neural networks and the extended generalised Pareto distribution, which allows us to model wildfire spread observed on an irregular spatial domain. We highlight the efficacy of our newly proposed model and perform a wildfire hazard assessment for Australia and population-dense communities, namely Tasmania, Sydney, Melbourne, and Perth.</p></div>","PeriodicalId":48771,"journal":{"name":"Spatial Statistics","volume":"59 ","pages":"Article 100811"},"PeriodicalIF":2.3,"publicationDate":"2024-01-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S2211675324000022/pdfft?md5=ae5571a2ce9cc73b6ee09f62cde7a192&pid=1-s2.0-S2211675324000022-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139487270","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 : 2024-01-08DOI: 10.1016/j.spasta.2023.100809
Yuanfeng Wang, Yunquan Song
{"title":"Variable selection via penalized quasi-maximum likelihood method for spatial autoregressive model with missing response","authors":"Yuanfeng Wang, Yunquan Song","doi":"10.1016/j.spasta.2023.100809","DOIUrl":"https://doi.org/10.1016/j.spasta.2023.100809","url":null,"abstract":"<div><p>Spatial autoregressive model is widely concerned in the economic field, whereas when the data is missing, variable selection and parameter estimation of the model is quite challenging. Based on this, we discuss the variable selection in spatial autoregressive model with missing data. Under the condition that errors are independent and identically distributed, we have developed a penalized quasi-maximum likelihood method to achieve variable selection and parameter estimation simultaneously in the presence of missing responses. The method’s theoretical properties, including consistency and asymptotical normality, are established under certain assumptions. Meanwhile, an improved expectation–maximization algorithm is provided for optimizing the penalized quasi-maximum likelihood function. Simulations are conducted to examine the proposed method and assess the finite-sample performance. Additionally, we present a practical example to illustrate the method’s application.</p></div>","PeriodicalId":48771,"journal":{"name":"Spatial Statistics","volume":"59 ","pages":"Article 100809"},"PeriodicalIF":2.3,"publicationDate":"2024-01-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S2211675323000842/pdfft?md5=559f1254900f3364b51f6f266c8dfb77&pid=1-s2.0-S2211675323000842-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139433320","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}