373 Autologistic Regression Models for Spatio-Temporal Binary Data
constants in spatial point processes, which Zheng and Zhu (2008) used for computing the
MCMLE. Friel et al. (2009) proposed a fast computation method for the estimation of the
normalizing constant based on a reduced dependence approximation of the likelihood
function. Later, we describe statistical inference based on MPLE, MCMLE, and Bayesian
hierarchical modeling.
17.2.2.1 Maximum Pseudo-Likelihood Estimation
Maximum pseudo-likelihood, rst introduced by Besag (1975) for autologistic models, is a
popular approach to the statistical inference for autologistic regression models. The MPLE
is the value of θ that maximizes the product of the full conditional distributions,
θ
˜
= arg max
θ
L
PL
(Y; θ),
where the pseudo-likelihood function for a spatio-temporal autologistic model is
L
PL
(Y; θ) = p(y
i,t
|y
i
,t
: (i
, t
) = (i, t))
i,t
exp{
k
p
=0
θ
k
x
k,i,t
y
i,t
+
j∈N
i
θ
p+1
y
i,t
y
j,t
+ θ
p+2
y
i,t
(y
i,t−1
+ y
i,t+1
)}
=
p
.
i,t
1 + exp{
k=0
θ
k
x
k,i,t
+
j∈N
i
θ
p+1
y
j,t
+ θ
p+2
(y
i,t−1
+ y
i,t+1
)}
(17.18)
Although the pseudo-likelihood function (17.18) is not the true likelihood except in the
trivial case of spatio-temporal independence, it can be shown that MPLEs are consistent
and asymptotically normal under suitable regularity conditions (Guyon, 1995).
To maximize the pseudo-likelihood function and obtain the MPLE of θ, it is straightfor-
ward to apply the standard logistic regression that assumes independence, which can be
implemented by, for example, proc logistic in SAS or the function glm in R. The corre-
sponding standard errors and approximate condence intervals can be obtained by a para-
metric bootstrap. Specically, in the parametric bootstrap, M resamples of spatio-temporal
binary responses are drawn according to the spatio-temporal autologistic regression model
using Gibbs sampling or perfect sampling. For each resample, an MPLE is computed and
the M resampled MPLEs are used to obtain an estimate of the variance of the MPLE based
on the original data. In particular, perfect sampling uses coupling and upon coalescence of
the coupled Markov chains, the resulting Monte Carlo samples are guaranteed to be from
the target distribution (e.g., Propp and Wilson, 1996; Møller, 1999).
17.2.2.2 Monte Carlo Maximum Likelihood Estimation
The maximum pseudo-likelihood approach is computationally efcient, but is statisti-
cally less efcient than maximum likelihood (Gumpertz et al., 1997; Wu and Huffer, 1997;
Zheng and Zhu, 2008). An alternative approach is Monte Carlo maximum likelihood
(MCML), where the normalizing constant is approximated using MCMC and thus direct
maximization of likelihood function can be obtained.