Likelihood Estimation and Inference for the Autologistic Model

The autologistic model is commonly used to model spatial binary data on the lattice. However, if the lattice size is too large, then exact calculation of its normalizing constant poses a major difficulty. Various different methods for estimation of model parameters, such as pseudo-likelihood, have been proposed to overcome this problem. This article presents a method to estimate the normalizing constant in an efficient manner. In particular, this allows tasks such as maximum likelihood estimation and inference for model parameters. We also consider the true likelihood approximated by the product of likelihoods for which the normalizing constant can be found by an analytic computational method by wrapping the lattice on the cylinder. This gives a simulation-free method of inference. We compare estimates of model parameters based on our new methods with the commonly used pseudo-likelihood approach. Although we have not considered Bayesian inferences here, the method can be straightforwardly extended to find posterior distributions. We apply our methods to the well-known endive data and to simulated data and find that our methods give substantially increased accuracy of estimation of model parameters.

[1]  David Bruce Wilson,et al.  Exact sampling with coupled Markov chains and applications to statistical mechanics , 1996, Random Struct. Algorithms.

[2]  C. Geyer,et al.  Constrained Monte Carlo Maximum Likelihood for Dependent Data , 1992 .

[3]  Xiao-Li Meng,et al.  Simulating Normalizing Constants: From Importance Sampling to Bridge Sampling to Path Sampling , 1998 .

[4]  J. Heikkinen,et al.  Fully Bayesian Approach to Image Restoration with an Application in Biogeography , 1994 .

[5]  F. Huang,et al.  Improvements of the Maximum Pseudo-Likelihood Estimators in Various Spatial Statistical Models , 1999 .

[6]  Yosihiko Ogata,et al.  Comparison of Two Methods for Calculating the Partition Functions of Various Spatial Statistical Models , 2001 .

[7]  M. Gu,et al.  Maximum likelihood estimation for spatial models by Markov chain Monte Carlo stochastic approximation , 2001 .

[8]  J. Besag Statistical Analysis of Non-Lattice Data , 1975 .

[9]  B. Ripley Statistical inference for spatial processes , 1990 .

[10]  J. Besag Spatial Interaction and the Statistical Analysis of Lattice Systems , 1974 .

[11]  Anthony N. Pettitt,et al.  Efficient calculation of the normalizing constant of the autologistic and related models on the cylinder and lattice , 2003 .

[12]  Anthony N. Pettitt,et al.  Binary probability maps using a hidden conditional autoregressive Gaussian process with an application to Finnish common toad data , 2000 .

[13]  Hulin Wu,et al.  Modelling the distribution of plant species using the autologistic regression model , 1997, Environmental and Ecological Statistics.

[14]  J. Besag On the Statistical Analysis of Dirty Pictures , 1986 .

[15]  J. Besag Nearest‐Neighbour Systems and the Auto‐Logistic Model for Binary Data , 1972 .

[16]  S. T. Buckland,et al.  An autologistic model for the spatial distribution of wildlife , 1996 .

[17]  Y. Ogata A Monte Carlo method for high dimensional integration , 1989 .

[18]  J. Propp,et al.  Exact sampling with coupled Markov chains and applications to statistical mechanics , 1996 .

[19]  J.-P. Dubois,et al.  Assessing the risk of soil contamination in the Swiss Jura using indicator geostatistics , 1997, Environmental and Ecological Statistics.