Abstract:In this era of big data, all scientific disciplines are evolving fast to cope up with the enormity of the available information. So is statistics, the queen of science. Big data are particularly relevant to spatio-temporal statistics, thanks to much-improved technology in satellite based remote sensing and Geographical Information Systems. However, none of the existing approaches seem to meet the simultaneous demand of reality emulation and cheap computation. In this article, with the Levy random fields as the starting point, e construct a new Bayesian nonparametric, nonstationary and nonseparable dynamic spatio- temporal model with the additional realistic property that the lagged spatio-temporal correlations converge to zero as the lag tends to infinity. Although our Bayesian model seems to be intricately structured and is variable-dimensional with respect to each time index, we are able to devise a fast and efficient parallel Markov Chain Monte Carlo (MCMC) algorithm for Bayesian inference. Our simulation experiment brings out quite encouraging performance from our Bayesian Levy-dynamic approach. We finally apply our Bayesian Levy-dynamic model and methods to a sea surface temperature dataset consisting of 139,300 data points in space and time. Although not big data in the true sense, this is a large and highly structured data by any standard. Even for this large and complex data, our parallel MCMC algorithm, implemented on 80 processors, generated 110,000 MCMC realizations from the Levy-dynamic posterior within a single day, and the resultant Bayesian posterior predictive analysis turned out to be encouraging. Thus, it is not unreasonable to expect that with significantly more computing resources, it is feasible to analyse terabytes of spatio-temporal data with our new model and methods.

Stationary process approximation for the analysis of large spatial datasets

Gaussian predictive process models for large spatial data sets

Improving the performance of predictive process modeling for large datasets

Nonstationary Spatial Process Models with Spatially Varying Covariance Kernels

Bayesian statistics in spatial epidemiology]

Hierarchical Nearest-Neighbor Gaussian Process Models for Large Geostatistical Datasets

Adaptive Bayesian nonstationary modeling for large spatial datasets using covariance approximations

A Scalable Partitioned Approach to Model Massive Nonstationary Non-Gaussian Spatial Datasets

High-Dimensional Bayesian Geostatistics

Bayesian Nonstationary Spatial Modeling for Very Large Datasets

High-dimensional Multivariate Geostatistics: A Bayesian Matrix-Normal Approach

Spatial Factor Modeling: A Bayesian Matrix-Normal Approach for Misaligned Data

Bayesian inference for high-dimensional nonstationary Gaussian processes

Full-scale approximations of spatio-temporal covariance models for large datasets

Spatial Process Approximations: Assessing Their Necessity

Sparse Approximate Inference for Spatio-Temporal Point Process Models

Bayesian Levy-Dynamic Spatio-Temporal Process: Towards Big Data Analysis

A full scale approximation of covariance functions for large spatial data sets

Multivariate Hawkes processes with spatial covariates for spatiotemporal event data analysis

A Case Study Competition Among Methods for Analyzing Large Spatial Data

A Variational Approach for Modeling High-dimensional Spatial Generalized Linear Mixed Models