A mixed mover-stayer model for spatiotemporal two-state processes

Studies of recurring infection or chronic disease often collect longitudinal data on the disease status of subjects. Two-state transitional models are useful for analysis in such studies where, at any point in time, an individual may be said to occupy either a diseased or disease-free state and interest centers on the transition process between states. Here, two additional features are present. The data are spatially arranged and it is important to account for spatial correlation in the transitional processes corresponding to different subjects. In addition there are subgroups of individuals with different mechanisms of transitions. These subgroups are not known a priori and hence group membership must be estimated. Covariates modulating transitions are included in a logistic additive framework. Inference for the resulting mixture spatial Markov regression model is not straightforward. We develop here a Monte Carlo expectation maximization algorithm for maximum likelihood estimation and a Markov chain Monte Carlo sampling scheme for summarizing the posterior distribution in a Bayesian analysis. The methodology is applied to a study of recurrent weevil infestation in British Columbia forests.