Joint spatial Bayesian modeling for studies combining longitudinal and cross-sectional data

Design for intervention studies may combine longitudinal data collected from sampled locations over several survey rounds and cross-sectional data from other locations in the study area. In this case, modeling the impact of the intervention requires an approach that can accommodate both types of data, accounting for the dependence between individuals followed up over time. Inadequate modeling can mask intervention effects, with serious implications for policy making. In this paper we use data from a large-scale larviciding intervention for malaria control implemented in Dar es Salaam, United Republic of Tanzania, collected over a period of almost 5 years. We apply a longitudinal Bayesian spatial model to the Dar es Salaam data, combining follow-up and cross-sectional data, treating the correlation in longitudinal observations separately, and controlling for potential confounders. An innovative feature of this modeling is the use of Ornstein–Uhlenbeck process to model random time effects. We contrast the results with other Bayesian modeling formulations, including cross-sectional approaches that consider individual-level random effects to account for subjects followed up in two or more surveys. The longitudinal modeling approach indicates that the intervention significantly reduced the prevalence of malaria infection in Dar es Salaam by 20% whereas the joint model did not suggest significance within the results. Our results suggest that the longitudinal model is to be preferred when longitudinal information is available at the individual level.

[1]  L. Pritchett,et al.  Estimating Wealth Effects Without Expenditure Data—Or Tears: An Application To Educational Enrollments In States Of India* , 2001, Demography.

[2]  H. Rue,et al.  Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations , 2009 .

[3]  P. Guttorp.  Stochastic modeling of scientific data , 1995 .

[4]  H. Rue,et al.  An explicit link between Gaussian fields and Gaussian Markov random fields: the stochastic partial differential equation approach , 2011 .

[5]  G. Killeen,et al.  A tool box for operational mosquito larval control: preliminary results and early lessons from the Urban Malaria Control Programme in Dar es Salaam, Tanzania , 2008, Malaria Journal.

[6]  Gerry F. Killeen,et al.  The Importance of Drains for the Larval Development of Lymphatic Filariasis and Malaria Vectors in Dar es Salaam, United Republic of Tanzania , 2010, PLoS neglected tropical diseases.

[7]  S. Ahmed.  Markov Chain Monte Carlo: Stochastic Simulation for Bayesian Inference (2nd ed.), by Dani Gamerman and Hedibert F. Lopes , 2008 .

[8]  Stefan Dongus,et al.  Participatory mapping of target areas to enable operational larval source management to suppress malaria vector mosquitoes in Dar es Salaam, Tanzania , 2007, International journal of health geographics.

[9]  J. Besag,et al.  Bayesian image restoration, with two applications in spatial statistics , 1991 .

[10]  G. Killeen,et al.  Microbial Larvicide Application by a Large-Scale, Community-Based Program Reduces Malaria Infection Prevalence in Urban Dar Es Salaam, Tanzania , 2009, PloS one.

[11]  Gianluca Baio,et al.  Spatial and spatio-temporal models with R-INLA. , 2013, Spatial and spatio-temporal epidemiology.

[12]  Leonhard Held,et al.  Spatio‐temporal disease mapping using INLA , 2011 .

[13]  M. Castro,et al.  Impact of Community-Based Larviciding on the Prevalence of Malaria Infection in Dar es Salaam, Tanzania , 2013, PloS one.

[14]  Roderick J. A. Little,et al.  Statistical Analysis with Missing Data: Little/Statistical Analysis with Missing Data , 2002 .

[15]  Neal Alexander,et al.  Bayesian Disease Mapping: Hierarchical Modeling in Spatial Epidemiology , 2011 .

[16]  L Knorr-Held,et al.  Bayesian modelling of inseparable space-time variation in disease risk. , 2000, Statistics in medicine.