Discrete Dynamics in Nature and Society

Volume 2014 (2014), Article ID 380392, 6 pages

http://dx.doi.org/10.1155/2014/380392

Research Article

Spectral Methods in Spatial Statistics

Kun Chen,¹ Lianmin Zhang,² and Maolin Pan³

¹School of Statistics, Southwestern University of Finance and Economics, Chengdu 611130, China

²School of Management and Engineering, Nanjing University, Nanjing 210093, China

³Department of Mathematics, Nanjing University, Nanjing 210093, China

Received 2 April 2014; Accepted 27 May 2014; Published 15 June 2014

Academic Editor: Xiaolin Xu

Copyright © 2014 Kun Chen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

When the spatial location area increases becoming extremely large, it is very difficult, if not possible, to evaluate the covariance matrix determined by the set of location distance even for gridded stationary Gaussian process. To alleviate the numerical challenges, we construct a nonparametric estimator called periodogram of spatial version to represent the sample property in frequency domain, because periodogram requires less computational operation by fast Fourier transform algorithm. Under some regularity conditions on the process, we investigate the asymptotic unbiasedness property of periodogram as estimator of the spectral density function and achieve the convergence rate.

1. Introduction

Recently, spatial statistics have attracted a number of researches from statisticians, geostatisticians, engineers, and econometricians [1–4]. The covariance clustering phenomena of neighborhood observations have always been the fundamental feature of spatial statistics, which means that closer observations have larger correlation and the correlation decays to zero for far separated ones. Therefore, the location information of spatial sited data plays a key role in constructing the probabilistic or statistical covariance function [5]. Using maximum likelihood estimation for the regularly spaced Gaussian data, the calculation of the covariance matrix inverse or determinant needs operations with the sample size . As becomes large, it is difficult to implement statistical inference due to computational limitations. Fortunately, spectral method is a powerful tool to analyze this complex spatial structure using fewer calculations with operations [6, 7].

This spectral method is widely applied in time series analysis [8], but little research has been done in spatial analysis. The spectral analysis, also known as harmonic analysis, is mainly based on the study of the spectral density function, which represents the underlying process property analog to covariance function. By spectral representation theorem, a positive definite covariance function in spatial domain is the Fourier transform of spectral density function in frequency domain. Operations performed on one domain have the corresponding operations on the other domain. In spatial domain, we usually use moment estimator for statistical inference. Under some regular assumptions, like ergodic or stationary, the nonparametric moment method provides asymptotically unbiased estimator. Because of the relationship between spatial and frequency domain, we can use the Fourier transform of moment estimator to make inference. The most important estimator in this sense is periodogram, which is asymptotically unbiased estimator of the spectral density function in time series analysis. Because the periodogram is the square of discrete Fourier transform (DFT), which decomposes the original sequence into components on different frequencies, in this case, the fast Fourier transform (FFT) provides a method to compute the DFT taking operations.

However, it is unknown whether the periodogram estimators have asymptotic unbiased property for spatial process. In this paper, based on the spatial version periodogram, we prove that the result still holds in spatial analysis and obtain that the convergence rate is a function of the sampling domain.

The periodogram method to study properties of regular stationary time series can be traced to Whittle [9], Ronsenblatt [10], and Priestly [11]. In some cases, the autocorrelation of time series decays very slowly, such that a pair of observations, when even separated far away, are still highly correlated. To overcome this difficulty, Heyde and Gay [12] introduced a smoothed version periodogram to estimate time series with long-range dependence. Geweke and Porter-Hudak [13] employed the log-periodogram as the response variable to construct a regression model for long memory time series.

The paper is organized as follows. Section 2 presents the basic definitions of periodogram and spectral density of spatial version. In Section 3, we give the main theorems on the asymptotic properties of periodogram. Concluding remarks and further research will be given in Section 4. Throughout the paper, the following notation will be adopted: denotes a term (a random variable) that is bounded (in probability); denotes a term (a random variable) that converges to zero (in probability); and denotes the Euclidean space with dimension .

2. Spectral Domain

We consider a weakly stationary Gaussian process with mean zero in which data are drawn from the regularly spaced grid on -dimensional Euclidean space, (see Figure 1). Here, we also call a random field on . A random field is weakly stationary if its mean is constant, and it has finite covariance function , such that , for any single location belonging to domain . Let be an orthogonal process on a regular gridded lattice, which is mean zero and has independent increments

for , with variance , where is the spectral density function. By spectral representation theorem [14], the real-valued weakly stationary random field with mean 0 has representation

Then the covariance function can be represented as

Figure 1: Simulation of regular spaced lattice data.

If the covariance function satisfies , the spectral density is the Fourier transform of the covariance function

If we further assume an additional assumption as in Fuentes [7],

(a),

where is the Euclidean norm given by , then the spectral density of has uniformly bounded first-order derivative. To simplify our discussion, we illustrate the idea by random fields on a two-dimensional grid; that is, , in what follows.

We estimate the spectral density function on the lattice by periodogram of spatial version

where is the total sample size; ; and is the spatial frequencies, computed at and for

In time series study, under some regular conditions on the stationary process, the periodogram is asymptotically unbiased estimator of spectral density; that is, as sample size , where the bias is of order . However, the periodogram is not a consistent estimator of the spectral density [8].

3. Main Results

Analogous to one-dimensional equally spaced time series observations, define the DFT on two-dimensional grid as

where and denote the real and imaginary parts, respectively, and

Since it is known that

the expectation of periodogram is equal to the variance of real and imaginary parts of DFT. To achieve this aim, we first show some asymptotic properties of Dirichlet kernel in time series. By simple calculation, we have

Also, Anderson [15] shows that

Also, the Feer kernel defined as holds an equation

Using the above formula (10), we obtain the next theorem on asymptotic properties of real and imaginary parts of DFT.

Theorem 1. Under Assumption (a), the variance value of real and imaginary parts of DFT is given by

where

where .

Proof. Consider

Applying the similar calculation procedure to , we can get the result that .

By the conclusion of Theorem 1 and formula (9), it follows that

Therefore, the bias of periodogram to spectral density is

The next theorem evaluates the bias in (17), such that we can prove the asymptotic unbiasedness property of periodogram for regular lattice random fields.

Theorem 2. Under Assumption , for , we have

Proof. By formula (17), it follows that

where and denote the first- and second-order derivatives of at and , respectively, where . By the fact that and by integration range , we have

Also,

because the integrand is odd, and the remaining terms hold as

because , for .

The idea of Theorem 2 can be directly extended to random fields on -dimensional regular gridded lattice. In conclusion, the periodogram of spatial version is asymptotically unbiased estimator of the spectral density function, in which the convergence rate is the minimum of sampling dimension, that is, .

4. Conclusions and Future Research

In this paper, we consider the nonparametric method to estimate the spectral density function. The spectral density function represents the population property and describes the behavior of underlying process. By constructing a moment estimator periodogram in frequency domain, which is a Fourier transform of covariance function in space domain, the spectral density can be estimated to be asymptotically unbiased. The spectral method outperforms traditional spatial method, because it provides significant computational advantages. In addition, we obtain the convergence rate as the minimum of spatial sampling dimension in the lattice. However, the consistency and asymptotic normality has not been established in this paper, which is important in statistical inference, hypothesis testing, and interval forecasts. On the other hand, by using the Whittle likelihood [9], we can approximate the maximum likelihood method [16] in spatial domain by Whittle estimator in frequency domain and thus estimate the specific parameter in the spectral density.

However, our asymptotic results are based on increasing domain asymptotics, allowing the fact that, as the sample size increases, the distance increases becoming sufficiently large, and thus the observations are nearly independent. However, if one is concerned to predict the value of a point in a bounded area, it is necessary to consider fixed domain asymptotics rather than increasing domain asymptotics. Specifically, even for regularly sited data, results from maximum likelihood estimate only exist in very limited cases. The difficulty under the fixed domain asymptotics is that there is at least one function of parameters that cannot be consistently estimated even for regularly gridded Gaussian process. This is indeed the case for Gaussian random fields with isotropic exponential covariance function , where is the sill parameter and is the range parameter. Ying [17] proved that when , the maximum likelihood estimator of is strongly consistent, but not possible for either or . Therefore, it needs more researches on the fixed domain asymptotics in spatial statistics.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors thank the guest editor and two anonymous referees for their helpful comments and suggestions. This work was partially supported by the National Natural Science Foundation of China (Grant no. 71301023), the Humanities and Social Sciences Foundation of Ministry of Education of China (Grant no. 13YJC630197), and the Basic Research Foundation (Natural Science) of Jiangsu Province (Grant no. BK20130582).

References

L. Anselin, Spatial Econometrics: Methods and Models, Kluwer Academic Publishers, Dordrecht, The Netherlands, 1988.

N. A. C. Cressie, Statistics for Spatial Data, Revised Edition, John Wiley & Sons, New York, NY, USA, 1993.

P. J. Diggle, J. A. Tawn, and R. A. Moyeed, “Model-based geostatistics,” Journal of the Royal Statistical Society C:. Applied Statistics, vol. 47, no. 3, pp. 299–350, 1998.

B. D. Ripley, Spatial Statistics, John Wiley & Sons, New York, NY, USA, 1981.

J. Besag, “Spatial interaction and the statistical analysis of lattice systems,” Journal of the Royal Statistical Society B: Methodological, vol. 36, pp. 192–236, 1974.

P. Whittle, “On stationary processes in the plane,” Biometrika, vol. 41, pp. 434–449, 1954.

M. Fuentes, “Approximate likelihood for large irregularly spaced spatial data,” Journal of the American Statistical Association, vol. 102, no. 477, pp. 321–331, 2007.

D. R. Brillinger, Time Series Data Analysis and Theory, Holden-Day, San Francisco, Calif, USA, 2nd edition, 1981.

P. Whittle, “Estimation and information in stationary time series,” Arkiv för Matematik, vol. 2, pp. 423–434, 1953.

M. R. Rosenblatt, Stationary Sequences and Random Fields, Birkhäuser, Boston, Mass, USA, 1985.

M. B. Priestley, Spectral Analysis and Time Series, Academic Press, 1981.

C. C. Heyde and R. Gay, “Smoothed periodogram asymptotics and estimation for processes and fields with possible long-range dependence,” Stochastic Processes and Their Applications, vol. 45, no. 1, pp. 169–182, 1993.

J. Geweke and S. Porter-Hudak, “The estimation and application of long memory time series models,” Journal of Time Series Analysis C, vol. 4, no. 4, pp. 221–238, 1983.

A. M. Yaglom, Correlation Theory of Stationary and Related Random Functions, Springer, New York, NY, USA, 1987.

T. W. Anderson, The Statistical Analysis of Time Series, John Wiley & Sons, New York, NY, USA, 1971.

K. V. Mardia and R. J. Marshall, “Maximum likelihood estimation of models for residual covariance in spatial regression,” Biometrika, vol. 71, no. 1, pp. 135–146, 1984.

Z. Ying, “Asymptotic properties of a maximum likelihood estimator with data from a Gaussian process,” Journal of Multivariate Analysis, vol. 36, no. 2, pp. 280–296, 1991.