Journal of Probability and Statistics

Volume 2009 (2009), Article ID 847830, 11 pages

http://dx.doi.org/10.1155/2009/847830

## A Note on the Properties of Generalised Separable Spatial Autoregressive Process

^{1}Department of Mathematics, Faculty of Science, University Putra Malaysia, 43400 Serdang, Malaysia^{2}Applied & Computational Statistics Laboratory, Institute for Mathematical Research, University Putra Malaysia, 43400 Serdang, Malaysia^{3}School of Mathematics and Statistics, The University of Sydney, NSW 2006, Australia

Received 27 February 2009; Accepted 20 July 2009

Academic Editor: Murray Clayton

Copyright © 2009 Mahendran Shitan and Shelton Peiris. 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

Spatial modelling has its applications in many fields like geology, agriculture, meteorology, geography, and so forth. In time series a class of models known as Generalised Autoregressive (GAR) has been introduced by Peiris (2003) that includes an index parameter . It has been shown that the inclusion of this additional parameter aids in modelling and forecasting many real data sets. This paper studies the properties of a new class of spatial autoregressive process of order 1 with an index. We will call this a *Generalised Separable Spatial Autoregressive* (GENSSAR) Model. The spectral density function (SDF), the autocovariance function (ACVF), and the autocorrelation function (ACF) are derived. The theoretical ACF and SDF plots are presented as three-dimensional figures.

#### 1. Introduction

Spatial modelling has its applications in many fields like geology, agriculture, meteorology, geography, and so forth. Spatial data can be classified as geostatistical data, lattice data, or point patterns. These differences are due to whether the spatial data has been observed on a continuous domain or at discrete locations. In point pattern analysis the domain is random and interest focuses on the location of events.

In this paper we concentrate on lattice data observed on a regular grid. Many models have been suggested in modelling spatial dependence like the Simultaneous Autoregression (SAR) [1], Conditional Autoregression (CAR) [2, 3], Moving Average (MA) [4], and Unilateral models [5].

For a two-dimensional *stationary process* we have the following definitions. Let be a sequence of spatial observations on a two-dimensional *regular* lattice. The mean function is E (a constant). The autocovariance function is = Cov and the autocorrelation function is given as .

Now, there exists a class of models that are known as *separable* models which have the property of a reflection symmetric correlation structure (i.e., ). The linear by linear process is defined as a stationary process where the autocovariance generating function of is defined as proportional to the product of two one-dimensional processes, and (see [6]) and the relationship may be represented as . As a consequence, its correlation structure can be expressed as a product of correlations (i.e., ). Basawa et al. [7] have considered separable models on a -dimensional lattice and have shown that the correlation structure is , where is the lag vector .

On the other hand, in the area of time series a class of models known as generalised autoregressive (GAR) models has been introduced by Peiris [8] by including an additional index parameter, . This is a natural extension of the standard AR model. It has been shown in Peiris [8] and Peiris et al*.* [9] that the additional index parameter plays an important role in modelling and forecasting real data sets. Shitan and Peiris [10] have also studied the estimation problem of the GAR model with a simulation study.

In this paper we will consider a special type of spatial model called a *Generalised Separable Spatial Autoregressive* (GENSSAR) Model. Some of its properties are discussed in Section 2. Finally in Section 3, some conclusions are drawn.

*The GENSSAR Model*

Let be a sequence of spatial observations on a two-dimensional *regular* lattice that satisfies
where is the usual backward shift operator acting in the th direction, is the backward shift operator acting in the th direction, and is a two-dimensional white noise process with mean zero and variance.

The term can be factored out as and hence (1.1) can be written as

The inclusion of the extra index parameter generalises the standard separable spatial model. Hence, we call the model defined in (1.1) as the Generalised Separable Spatial Autoregressive model or GENSSAR(1,1) model.

The following section reports some of its properties in detail.

#### 2. Some Properties of GENSSAR(1,1)

The solution of this process in (1.1) is given in Proposition 1.

Proposition 1. *For a process defined in (1.1), the solution is
**
where is the gamma function.*

*Proof. *From (1.2), we have Using binomial expansion, it follows that
Note that
Similarly
Substituting (2.3) and (2.4) into (2.2) and upon simplification completes the proof.

Proposition 2. *For a process defined in (1.1), the spectral density is given as
*

*Proof. *The proof is established by simplifying the following expression:

The following proposition provides an expression for the autocovariance function of the GENSSAR(1,1) process.

Proposition 3. *For a process defined in (1.1) the autocovariance of the process is given as
**
where is the hypergeometric function.*

*Proof. *We establish the previous proposition by integrating the spectral density as given in Proposition 2
Now by making use of the identity (see [8]),
we obtain
which completes the proof.

Corollary 4. *For a process defined in (1.1) the variance of the process is given as
*

*Proof. *This result is directly from Proposition 3 by letting .

The autocorrelation function (ACF) of the model in (1.1) is given as

*Remark 5. *Note when , we have the standard separable spatial model. Substituting in Proposition 3, we obtain

Using the following identity (see Abramowitz and Stegun [11, Page 556, Identity No. 15.1.8]): (2.13) reduces to Hence, Proposition 3 reduces to the autocovariance function of the standard separable spatial model when .

In Table 1, we have tabulated (to three decimal places) the ACF, computed by using (2.12) with , , . This is the standard separable model. Clearly, we can see that the numerical values computed by using (2.12) agree with the ACF of the standard separable model which is . Hence, this verifies (2.12).

In Table 2, we have tabulated (to three decimal places) the ACF, of the GENSSAR model () computed by using (2.12) with , , . While Table 3 shows the ACF values of the GENSSAR model () computed by (2.12) with , , .

Figures 1, 2, and 3 show the ACF for the three models considered in this paper.

From the tables and figures we can clearly see that the behaviour of ACF depends on the index parameter . When , the ACF decays slower than that of the standard separable model. On the other hand, when the ACF decays faster than the standard model. Hence, the GENSSAR model can be used to model many types of autocorrelation structure.

We also considered a further illustrative example when and were not equal to each other. That is, we chose the parameter values to be , , . In Table 4, we have tabulated (to three decimal places) the ACF, of the GENSSAR model () computed by using (2.12) with , , , and Figure 4 shows a plot of the ACF. Clearly we can see that the decay in the autocorrelation is more rapid along the axis as compared to the axis. Hence, we can model data whose autocorrelations decay at different rates in different directions.

For the models considered in this paper, the two-dimensional spectral densities for various parameter values are shown in Figures 5, 6, 7, 8, 9, 10, and 11.

#### 3. Conclusion

The objective of this research is to introduce a new class of models called *GENSSAR* models by including an additional index parameter and to establish some of its properties. We have established the autocovariance function. The GENSSAR(1,1) model is a more general model than the standard separable spatial AR(1,1) process. Due to the generality of this model, it is a useful model.

The authors are working on the other aspects of this model with applications and will be reported in a future paper.

#### Acknowledgments

The authors are very grateful to the reviewer and the editor for their valuable comments and suggestions to improve the quality of this paper. They express their thanks to the Department of Mathematics and the Institute of Mathematical Research, University Putra Malaysia for their support. They also wish to thank the School of Mathematics and Statistics, The University of Sydney, for their support during the first author's visit in 2008.

#### References

- P. Whittle, “On stationary processes in the plane,”
*Biometrika*, vol. 41, pp. 434–449, 1954. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. S. Bartlett, “Physical nearest neighbour models and non-linear time series,”
*Journal of Applied Probability*, vol. 8, pp. 222–232, 1971. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. E. Besag, “Spatial interaction and the statistical analysis of lattice systems,”
*Journal of the Royal Statistical Society B*, vol. 36, pp. 192–236, 1974. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. P. Haining, “The moving average model for spatial interaction,”
*Transactions of the Institute of British Geographers*, vol. 3, no. 2, pp. 202–225, 1978. View at Publisher · View at Google Scholar - S. Basu and G. C. Reinsel, “Properties of the spatial unilateral first order ARMA model,”
*Advances in Applied Probability*, vol. 25, no. 3, pp. 631–648, 1993. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. J. Martin, “A subclass of lattice processes applied to a problem in planar sampling,”
*Biometrika*, vol. 66, no. 2, pp. 209–217, 1979. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - I. V. Basawa, P. J. Brockwell, and V. M. Mandrekar, “Inference for spatial time series, computer science and statistics,” in
*Proceedings of the 22nd Symposium on the Interface*, Springer, New York, NY, USA, 1991. - M. S. Peiris, “Improving the quality of forecasting using generalized AR models: an application to statistical quality control,”
*Statistical Methods*, vol. 5, no. 2, pp. 156–171, 2003. View at Google Scholar · View at MathSciNet - S. Peiris, D. Allen, and A. Thavaneswaran, “An introduction to generalized moving average models and applications,”
*Journal of Applied Statistical Science*, vol. 13, no. 3, pp. 251–267, 2004. View at Google Scholar · View at MathSciNet - M. Shitan and S. Peiris, “Generalized autoregressive (GAR) model: a comparison of maximum likelihood and whittle estimation procedures using a simulation study,”
*Communications in Statistics Simulation and Computation*, vol. 37, no. 3, pp. 560–570, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - M. Abramowitz and I. Stegun,
*Handbook of Mathematical Functions*, Dover, New York, NY, USA, 1964.