Spatio-temporal clustering analysis using generalized lasso with an application to reveal the spread of Covid-19 cases in Japan

• Published in: Computational statistics 2024-05, Vol.39 (3), p.1513-1537
• Main Authors: Rahardiantoro, Septian, Sakamoto, Wataru
• Format: Article
• Language: English
• Subjects: Cluster analysis; Clustering; COVID-19; Economic Theory/Quantitative Economics/Mathematical Methods; Mathematical models; Mathematics and Statistics; Original Paper; Parameters; Probability and Statistics in Computer Science; Probability Theory and Stochastic Processes; Regularization; Statistics; Tuning
• ISSN: 0943-4062; 1613-9658
• DOI: 10.1007/s00180-023-01331-x

This study addressed the issue of determining multiple potential clusters with regularization approaches for the purpose of spatio-temporal clustering. The generalized lasso framework has flexibility to incorporate adjacencies between objects in the penalty matrix and to detect multiple clusters. A generalized lasso model with two L 1 penalties is proposed, which can be separated into two generalized lasso models: trend filtering of temporal effect and fused lasso of spatial effect for each time point. To select the tuning parameters, the approximate leave-one-out cross-validation (ALOCV) and generalized cross-validation (GCV) are considered. A simulation study is conducted to evaluate the proposed method compared to other approaches in different problems and structures of multiple clusters. The generalized lasso with ALOCV and GCV provided smaller MSE in estimating the temporal and spatial effect compared to unpenalized method, ridge, lasso, and generalized ridge. In temporal effects detection, the generalized lasso with ALOCV and GCV provided relatively smaller and more stable MSE than other methods, for different structure of true risk values. In spatial effects detection, the generalized lasso with ALOCV provided higher index of edges detection accuracy. The simulation also suggested using a common tuning parameter over all time points in spatial clustering. Finally, the proposed method was applied to the weekly Covid-19 data in Japan form March 21, 2020, to September 11, 2021, along with the interpretation of dynamic behavior of multiple clusters.