TY - CONF

T1 - A spatio-temporal model based on the SVD to analyze large spatio-temporal datasets

AU - Onorati, Rossella

PY - 2009

Y1 - 2009

N2 - A common problem in the analysis of space-time data is to compress a large dataset in order to extract the underlying trends. Empirical orthogonal function (EOF) analysis is a useful tool for examining both the temporal and the spatial variation in atmospherical and physical process and a convenient method of performing this is the Singular Value Decomposition (SVD). Many spatio-temporal models for measurements Z(s; t) at location s at time t, can be written as a sum of a systematic component and a residual component: Z = M+E, where Z, M and E are all T x N matrices. Our approach permits modeling ofincomplete data matrices using an "EM-like" iterative algorithm for the SVD. We model the trend, M, by linear combinations of smooth temporal basis functions derived from left (temporal) singular vectors of the SDV of Z with dimension of the model chosen by cross-validation. We further decompose by SVD the spatio-temporal residual matrix E computed as residuals from regressions at each site (column) of the observations on the smoothed temporal basis functions. Finally we fit an autoregressive model to the columns (time series) of residuals from the SVD of E. Our aim is to illustrate a simple model to characterize trends and model the variability in large spatio-temporal data matrices. The methodology is demostrated with a spatiotemporal dataset.

AB - A common problem in the analysis of space-time data is to compress a large dataset in order to extract the underlying trends. Empirical orthogonal function (EOF) analysis is a useful tool for examining both the temporal and the spatial variation in atmospherical and physical process and a convenient method of performing this is the Singular Value Decomposition (SVD). Many spatio-temporal models for measurements Z(s; t) at location s at time t, can be written as a sum of a systematic component and a residual component: Z = M+E, where Z, M and E are all T x N matrices. Our approach permits modeling ofincomplete data matrices using an "EM-like" iterative algorithm for the SVD. We model the trend, M, by linear combinations of smooth temporal basis functions derived from left (temporal) singular vectors of the SDV of Z with dimension of the model chosen by cross-validation. We further decompose by SVD the spatio-temporal residual matrix E computed as residuals from regressions at each site (column) of the observations on the smoothed temporal basis functions. Finally we fit an autoregressive model to the columns (time series) of residuals from the SVD of E. Our aim is to illustrate a simple model to characterize trends and model the variability in large spatio-temporal data matrices. The methodology is demostrated with a spatiotemporal dataset.

KW - SVD

KW - Spatio-temporal processes

KW - SVD

KW - Spatio-temporal processes

UR - http://hdl.handle.net/10447/52889

UR - http://www2.stat.unibo.it/ties2009/doc/bookTIES2009Bologna.pdf

M3 - Other

ER -