Hierarchical Agent-Based Spatio-Temporal Dynamic Models for Discrete-Valued Data 355
spatial heterogeneity in the interaction probabilities p that correlate with the gradient of a
potential eld. To briey describe how this generalization could be implemented, consider
a potential eld α(X) on the grid that depends on a set of variables X inuencing the spread
of disease. One possible hypothesis for disease spread might be that it responds to changes
in the potential eld α; for example, an abrupt change in a landscape feature such as a
boundary between land and water would lead to an increase in the gradient perpendicu-
lar to the direction of the boundary. On a discrete spatial support, such as the one we are
considering here for the rabies example, the rst-order gradient could be summarized as a
set of velocities a
i
, one in each direction from the cell of interest i to its neighboring cells in
N
i
. These gradient values a
i
(or some function of them) can then serve as hyperparameters
in the Dirichlet model for the spatially varying interaction probabilities p
i
. Of course, the
Dirichlet is only one way to model p
i
; other stochastic models or deterministic functions
could be used to link p
i
to a
i
,or X.
Returning now to the model described in (16.1–16.5), we can express the posterior dis-
tribution for all unknowns as proportional to the conditionally factored joint distribution
such that
T
m
[p, a, φ, ψ|{y
i,t
, ∀i, t}] ∝ [y
i,t
|θ
i,t
(p, φ, ψ)][p|a][a][φ][ψ] . (16.6)
t=1
i=1
Also, recall that the presence probabilities θ
i,t
are a function of the other model parameters;
thus, inference can be obtained for them as derived quantities in the model. For exam-
ple, the posterior mean and standard deviation of the space–time series for θ
i,t
are shown
in Figure 16.2. The left panel of Figure 16.2 shows the posterior mean for θ
i,t
and gives
us a quantitative understanding of presence probability along the front of the epidemic
over time. Similarly, the posterior standard deviation for θ
i,t
shown in the right panel of
Figure 16.2 allows us to visualize the uncertainty pertaining to presence probability for
all grid cells and times. Such statistical products could be used for short-term forecasting;
for example, forecasting the presence for the next year. Such forecasts require a reasonable
timescale to accommodate the required computation time.
Overall, the relatively simple statistical ABM presented in (16.1–16.5) represents a funda-
mentally different approach to modeling dynamics (i.e., bottom-up rather than top-down)
that is quite general and capable of representing complicated dynamics based on only a
simple set of rules describing the behavior among agents. As a reminder, in this example
we let the grid cells act as agents, but there is no reason why individual raccoons could not
serve as agents being explicitly modeled if sufcient individual-level data were collected
to gain the desired statistical inference. Scale is a critical component of all models, but it
seems especially important in statistical ABMs because we seek to invert the models and
are thus limited by the scale on which the data were collected.
16.4 Hierarchical First-Order Emulators and ABMs
Parameter estimation, calibration, and validation can be difcult in deterministic and
stochastic ABMs given the computational cost of simulation and, in most cases, the
nonlinear relationships in the parameters. The use of statistical emulators (or surrogates)