Abstract
The simultaneous spatiotemporal modeling of multiple related diseases strengthens inferences by borrowing information between related diseases. Numerous research contributions to spatiotemporal modeling approaches exhibit their strengths differently with increasing complexity. However, contributions that combine spatiotemporal approaches to modeling of multiple diseases simultaneously are not so common. We present a full Bayesian hierarchical spatio-temporal approach to the joint modeling of Human Immunodeficiency Virus and Tuberculosis incidences in Kenya. Using case notification data for the period 2012-2017, we estimated the model parameters and determined the joint spatial patterns and temporal variations. Our model included specific and shared spatial and temporal effects. The specific random effects allowed for departures from the shared patterns for the different diseases. The space-time interaction term characterized the underlying spatial patterns with every temporal fluctuation. We assumed the shared random effects to be the structured effects and the disease-specific random effects to be unstructured effects. We detected the spatial similarity in the distribution of Tuberculosis and Human Immunodeficiency Virus in approximately 29 counties around the western, central and southern regions of Kenya. The distribution of the shared relative risks had minimal difference with the Human Immunodeficiency Virus disease-specific relative risk whereas that of Tuberculosis presented many more counties as high-risk areas. The flexibility and informative outputs of Bayesian Hierarchical Models enabled us to identify the similarities and differences in the distribution of the relative risks associated with each disease. Estimating the Human Immunodeficiency Virus and Tuberculosis shared relative risks provide additional insights towards collaborative monitoring of the diseases and control efforts.