TY - JOUR

T1 - Efficient calculation of the normalizing constant of the autologistic and related models on the cylinder and lattice.

AU - Pettitt, Anthony

AU - Friel, N.

AU - Reeves, R.

N1 - RAE_import_type : Journal article RAE_uoa_type : Statistics and Operational Research

PY - 2003/2

Y1 - 2003/2

N2 - Motivated by the autologistic model for the analysis of spatial binary data on the two-dimensional lattice, we develop efficient computational methods for calculating the normalizing constant for models for discrete data defined on the cylinder and lattice. Because the normalizing constant is generally unknown analytically, statisticians have developed various ad hoc methods to overcome this difficulty. Our aim is to provide computationally and statistically efficient methods for calculating the normalizing constant so that efficient likelihood-based statistical methods are then available for inference. We extend the so-called transition method to find a feasible computational method of obtaining the normalizing constant for the cylinder boundary condition. To extend the result to the free-boundary condition on the lattice we use an efficient path sampling Markov chain Monte Carlo scheme. The methods are generally applicable to association patterns other than spatial, such as clustered binary data, and to variables taking three or more values described by, for example, Potts models.

AB - Motivated by the autologistic model for the analysis of spatial binary data on the two-dimensional lattice, we develop efficient computational methods for calculating the normalizing constant for models for discrete data defined on the cylinder and lattice. Because the normalizing constant is generally unknown analytically, statisticians have developed various ad hoc methods to overcome this difficulty. Our aim is to provide computationally and statistically efficient methods for calculating the normalizing constant so that efficient likelihood-based statistical methods are then available for inference. We extend the so-called transition method to find a feasible computational method of obtaining the normalizing constant for the cylinder boundary condition. To extend the result to the free-boundary condition on the lattice we use an efficient path sampling Markov chain Monte Carlo scheme. The methods are generally applicable to association patterns other than spatial, such as clustered binary data, and to variables taking three or more values described by, for example, Potts models.

U2 - 10.1111/1467-9868.00383

DO - 10.1111/1467-9868.00383

M3 - Journal article

VL - 65

SP - 235

EP - 246

JO - Journal of the Royal Statistical Society: Series B (Statistical Methodology)

JF - Journal of the Royal Statistical Society: Series B (Statistical Methodology)

SN - 1369-7412

IS - 1

ER -