The hospitalisation risk of chronic circulatory and respiratory diseases associated with coal mining in the general population in Queensland, Australia.

Queensland is the main coal mining state in Australia where populations in coal mining areas have been historically exposed to coal mining emissions. Although a higher risk of chronic circulatory and respiratory diseases has been associated with coal mining globally, few studies have investigated these associations in the Queensland general population. This study estimates the association of coal production with hospitalisations for chronic circulatory and respiratory diseases in Queensland considering spatial and temporal variations during 1997-2014. An ecological analysis used a Bayesian hierarchical spatiotemporal model to estimate the association of coal production with standardised rates of each, chronic circulatory and respiratory diseases, adjusting for sociodemographic factors and considering the spatial structure of Queensland's statistical areas (SA2) in the 18-year period. Two specifications; with and without a space-time interaction effect were compared using the integrated nested Laplace approximation -INLA approach. The posterior mean of the best fit model was used to map the spatial, temporal and spatiotemporal trends of risk. The analysis considered 2,831,121 hospitalisation records. Coal mining was associated with a 4 % (2.4-5.5) higher risk of hospitalisation for chronic respiratory diseases in the model with a space-time interaction effect which had the best fit. An emerging higher risk of either chronic circulatory and respiratory diseases was identified in eastern areas and some coal-mining areas in central and southeast Queensland. There were important disparities in the spatiotemporal trend of risk between coal -and non-coal mining areas for each, chronic circulatory and respiratory diseases. Coal mining is associated with an increased risk of chronic respiratory diseases in the Queensland general population. Bayesian spatiotemporal analyses are robust methods to identify environmental determinants of morbidity in exposed populations. This methodology helps identifying at-risk populations which can be useful to support decision-making in health. Future research is required to investigate the causality links between coal mining and these diseases.

Exploring Dengue Dynamics: A Multi-Scale Analysis of Spatio-Temporal Trends in Ibagué, Colombia.

Our study examines how dengue fever incidence is associated with spatial (demographic and socioeconomic) alongside temporal (environmental) factors at multiple scales in the city of Ibagué, located in the Andean region of Colombia. We used the dengue incidence in Ibagué from 2013 to 2018 to examine the associations with climate, socioeconomic, and demographic factors from the national census and satellite imagery at four levels of local spatial aggregation. We used geographically weighted regression (GWR) to identify the relevant socioeconomic and demographic predictors, and we then integrated them with environmental variables into hierarchical models using integrated nested Laplace approximation (INLA) to analyze the spatio-temporal interactions. Our findings show a significant effect of spatial variables across the different levels of aggregation, including human population density, gas and sewage connection, percentage of woman and children, and percentage of population with a higher education degree. Lagged temporal variables displayed consistent patterns across all levels of spatial aggregation, with higher temperatures and lower precipitation at short lags showing an increase in the relative risk (RR). A comparative evaluation of the models at different levels of aggregation revealed that, while higher aggregation levels often yield a better overall model fit, finer levels offer more detailed insights into the localized impacts of socioeconomic and demographic variables on dengue incidence. Our results underscore the importance of considering macro and micro-level factors in epidemiological modeling, and they highlight the potential for targeted public health interventions based on localized risk factor analyses. Notably, the intermediate levels emerged as the most informative, thereby balancing spatial heterogeneity and case distribution density, as well as providing a robust framework for understanding the spatial determinants of dengue.

Spatiotemporal joint analysis of PM2.5 and Ozone in California with INLA approach.

The substantial threat of concurrent air pollutants to public health is increasingly severe under climate change. To identify the common drivers and extent of spatiotemporal similarity of PM2.5 and ozone (O3), this paper proposed a log Gaussian-Gumbel Bayesian hierarchical model allowing for sharing a stochastic partial differential equation and autoregressive model of order one (SPDE-AR(1)) spatiotemporal interaction structure. The proposed model, implemented by the approach of integrated nested Laplace approximation (INLA), outperforms in terms of estimation accuracy and prediction capacity for its increased parsimony and reduced uncertainty, especially for the shared O3 sub-model. Besides the consistently significant influence of temperature (positive), extreme drought (positive), fire burnt area (positive), gross domestic product (GDP) per capita (positive), and wind speed (negative) on both PM2.5 and O3, surface pressure and precipitation demonstrate positive associations with PM2.5 and O3, respectively. While population density relates to neither. In addition, our results demonstrate similar spatiotemporal interactions between PM2.5 and O3, indicating that the spatial and temporal variations of these pollutants show relatively considerable consistency in California. Finally, with the aid of the excursion function, we see that the areas around the intersection of San Luis Obispo and Santa Barbara counties are likely to exceed the unhealthy O3 level for USG simultaneously with other areas throughout the year. Our findings provide new insights for regional and seasonal strategies in the co-control of PM2.5 and O3. Our methodology is expected to be utilized when interest lies in multiple interrelated processes in the fields of environment and epidemiology.

Spatiotemporal Changes in the Slavery-Inequality Relationship: The Diffusion of the Legacy of Slavery.

Despite the persistence of relationships between historical racist violence and contemporary Black-White inequality, research indicates, in broad strokes, that the slavery-inequality relationship in the United States has changed over time. Identifying the timing of such change across states can offer insights into the underlying processes that generate Black-White inequality. In this study, we use integrated nested Laplace approximation models to simultaneously account for spatial and temporal features of panel data for Southern counties during the period spanning 1900 to 2018, in combination with data on the concentration of enslaved people from the 1860 census. Results provide the first evidence on the timing of changes in the slavery-economic inequality relationship and how changes differ across states. We find a region-wide decline in the magnitude of the slavery-inequality relationship by 1930, with declines traversing the South in a northeasterly-to-southwesterly pattern over the study period. Different paces in declines in the relationship across states suggest the expansion of institutionalized racism first in places with the longest-standing overt systems of slavery. Results provide guidance for further identifying intervening mechanisms-most centrally, the maturity of racial hierarchies and the associated diffusion of racial oppression across institutions, and how they affect the legacy of slavery in the United States.

Spatiotemporal hierarchical Bayesian analysis to identify factors associated with COVID-19 in suburban areas in Colombia.

Introduction: The pandemic had a profound impact on the provision of health services in Cúcuta, Colombia where the neighbourhood-level risk of Covid-19 has not been investigated. Identifying the sociodemographic and environmental risk factors of Covid-19 in large cities is key to better estimate its morbidity risk and support health strategies targeting specific suburban areas. This study aims to identify the risk factors associated with the risk of Covid-19 in Cúcuta considering inter -spatial and temporal variations of the disease in the city's neighbourhoods between 2020 and 2022. Methods: Age-adjusted rate of Covid-19 were calculated in each Cúcuta neighbourhood and each quarter between 2020 and 2022. A hierarchical spatial Bayesian model was used to estimate the risk of Covid-19 adjusting for socioenvironmental factors per neighbourhood across the study period. Two spatiotemporal specifications were compared (a nonparametric temporal trend; with and without space-time interaction). The posterior mean of the spatial and spatiotemporal effects was used to map the Covid-19 risk. Results: There were 65,949 Covid-19 cases in the study period with a varying standardized Covid-19 rate that peaked in October-December 2020 and April-June 2021. Both models identified an association of the poverty and stringency indexes, education level and PM10 with Covid-19 although the best fit model with a space-time interaction estimated a strong association with the number of high-traffic roads only. The highest risk of Covid-19 was found in neighbourhoods in west, central, and east Cúcuta. Conclusions: The number of high-traffic roads is the most important risk factor of Covid-19 infection in Cucuta. This indicator of mobility and connectivity overrules other socioenvironmental factors when Bayesian models include a space-time interaction. Bayesian spatial models are important tools to identify significant determinants of Covid-19 and identifying at-risk neighbourhoods in large cities. Further research is needed to establish causal links between these factors and Covid-19.

Modelling spatiotemporal variation in under-five malaria risk in Ghana in 2016-2021.

BACKGROUND: Ghana is among the top 10 highest malaria burden countries, with about 20,000 children dying annually, 25% of which were under five years. This study aimed to produce interactive web-based disease spatial maps and identify the high-burden malaria districts in Ghana. METHODS: The study used 2016-2021 data extracted from the routine health service nationally representative and comprehensive District Health Information Management System II (DHIMS2) implemented by the Ghana Health Service. Bayesian geospatial modelling and interactive web-based spatial disease mapping methods were employed to quantify spatial variations and clustering in malaria risk across 260 districts. For each district, the study simultaneously mapped the observed malaria counts, district name, standardized incidence rate, and predicted relative risk and their associated standard errors using interactive web-based visualization methods. RESULTS: A total of 32,659,240 malaria cases were reported among children < 5 years from 2016 to 2021. For every 10% increase in the number of children, malaria risk increased by 0.039 (log-mean 0.95, 95% credible interval = - 13.82-15.73) and for every 10% increase in the number of males, malaria risk decreased by 0.075, albeit not statistically significant (log-mean - 1.82, 95% credible interval = - 16.59-12.95). The study found substantial spatial and temporal differences in malaria risk across the 260 districts. The predicted national relative risk was 1.25 (95% credible interval = 1.23, 1.27). The malaria risk is relatively the same over the entire year. However, a slightly higher relative risk was recorded in 2019 while in 2021, residing in Keta, Abuakwa South, Jomoro, Ahafo Ano South East, Tain, Nanumba North, and Tatale Sanguli districts was associated with the highest malaria risk ranging from a relative risk of 3.00 to 4.83. The district-level spatial patterns of malaria risks changed over time. CONCLUSION: This study identified high malaria risk districts in Ghana where urgent and targeted control efforts are required. Noticeable changes were also observed in malaria risk for certain districts over some periods in the study. The findings provide an effective, actionable tool to arm policymakers and programme managers in their efforts to reduce malaria risk and its associated morbidity and mortality in line with the Sustainable Development Goals (SDG) 3.2 for limited public health resource settings, where universal intervention across all districts is practically impossible.

Assessing the survival time of women with breast cancer in Northwestern Ethiopia: using the Bayesian approach.

BACKGROUND: Despite the significant weight of difficulty, Ethiopia's survival rate and mortality predictors have not yet been identified. Finding out what influences outpatient breast cancer patients' survival time was the major goal of this study. METHODS: A retrospective study was conducted on outpatients with breast cancer. In order to accomplish the goal, 382 outpatients with breast cancer were included in the study using information obtained from the medical records of patients registered at the University of Gondar referral hospital in Gondar, Ethiopia, between May 15, 2016, and May 15, 2020. In order to compare survival functions, Kaplan-Meier plots and the log-rank test were used. The Cox-PH model and Bayesian parametric survival models were then used to examine the survival time of breast cancer outpatients. The use of integrated layered Laplace approximation techniques has been made. RESULTS: The study included 382 outpatients with breast cancer in total, and 148 (38.7%) patients died. 42 months was the estimated median patient survival time. The Bayesian Weibull accelerated failure time model was determined to be suitable using model selection criteria. Stage, grade 2, 3, and 4, co-morbid, histological type, FIGO stage, chemotherapy, metastatic number 1, 2, and >=3, and tumour size all have a sizable impact on the survival time of outpatients with breast cancer, according to the results of this model. The breast cancer outpatient survival time was correctly predicted by the Bayesian Weibull accelerated failure time model. CONCLUSIONS: Compared to high- and middle-income countries, the overall survival rate was lower. Notable variables influencing the length of survival following a breast cancer diagnosis were weight loss, invasive medullar histology, comorbid disease, a large tumour size, an increase in metastases, an increase in the International Federation of Gynaecologists and Obstetricians stage, an increase in grade, lymphatic vascular space invasion, positive regional nodes, and late stages of cancer. The authors advise that it is preferable to increase the number of early screening programmes and treatment centres for breast cancer and to work with the public media to raise knowledge of the disease's prevention, screening, and treatment choices.

Joint quantile disease mapping with application to malaria and G6PD deficiency.

Statistical analysis based on quantile methods is more comprehensive, flexible and less sensitive to outliers when compared to mean methods. Joint disease mapping is useful for inferring correlation between different diseases. Most studies investigate this link through multiple correlated mean regressions. We propose a joint quantile regression framework for multiple diseases where different quantile levels can be considered. We are motivated by the theorized link between the presence of malaria and the gene deficiency G6PD, where medical scientists have anecdotally discovered a possible link between high levels of G6PD and lower than expected levels of malaria initially pointing towards the occurrence of G6PD inhibiting the occurrence of malaria. Thus, the need for flexible joint quantile regression in a disease mapping framework arises. Our model can be used for linear and nonlinear effects of covariates by stochastic splines since we define it as a latent Gaussian model. We perform Bayesian inference using the R integrated nested Laplace approximation, suitable even for large datasets. Finally, we illustrate the model's applicability by considering data from 21 countries, although better data are needed to prove a significant relationship. The proposed methodology offers a framework for future studies of interrelated disease phenomena.

A Gibbs-INLA algorithm for multidimensional graded response model analysis.

In this paper, we propose a novel Gibbs-INLA algorithm for the Bayesian inference of graded response models with ordinal response based on multidimensional item response theory. With the combination of the Gibbs sampling and the integrated nested Laplace approximation (INLA), the new framework avoids the cumbersome tuning which is inevitable in classical Markov chain Monte Carlo (MCMC) algorithm, and has low computing memory, high computational efficiency with much fewer iterations, and still achieve higher estimation accuracy. Therefore, it has the ability to handle large amount of multidimensional response data with different item responses. Simulation studies are conducted to compare with the Metroplis-Hastings Robbins-Monro (MH-RM) algorithm and an application to the study of the IPIP-NEO personality inventory data is given to assess the performance of the new algorithm. Extensions of the proposed algorithm for application on more complicated models and different data types are also discussed.

Forecasting freshwater cyanobacterial harmful algal blooms for Sentinel-3 satellite resolved U.S. lakes and reservoirs.

This forecasting approach may be useful for water managers and associated public health managers to predict near-term future high-risk cyanobacterial harmful algal blooms (cyanoHAB) occurrence. Freshwater cyanoHABs may grow to excessive concentrations and cause human, animal, and environmental health concerns in lakes and reservoirs. Knowledge of the timing and location of cyanoHAB events is important for water quality management of recreational and drinking water systems. No quantitative tool exists to forecast cyanoHABs across broad geographic scales and at regular intervals. Publicly available satellite monitoring has proven effective in detecting cyanobacteria biomass near-real time within the United States. Weekly cyanobacteria abundance was quantified from the Ocean and Land Colour Instrument (OLCI) onboard the Sentinel-3 satellite as the response variable. An Integrated Nested Laplace Approximation (INLA) hierarchical Bayesian spatiotemporal model was applied to forecast World Health Organization (WHO) recreation Alert Level 1 exceedance >12 µg L-1 chlorophyll-a with cyanobacteria dominance for 2192 satellite resolved lakes in the United States across nine climate zones. The INLA model was compared against support vector classifier and random forest machine learning models; and Dense Neural Network, Long Short-Term Memory (LSTM), Recurrent Neural Network (RNN), and Gneural Network (GNU) neural network models. Predictors were limited to data sources relevant to cyanobacterial growth, readily available on a weekly basis, and at the national scale for operational forecasting. Relevant predictors included water surface temperature, precipitation, and lake geomorphology. Overall, the INLA model outperformed the machine learning and neural network models with prediction accuracy of 90% with 88% sensitivity, 91% specificity, and 49% precision as demonstrated by training the model with data from 2017 through 2020 and independently assessing predictions with data from the 2021 calendar year. The probability of true positive responses was greater than false positive responses and the probability of true negative responses was less than false negative responses. This indicated the model correctly assigned lower probabilities of events when they didn't exceed the WHO Alert Level 1 threshold and assigned higher probabilities when events did exceed the threshold. The INLA model was robust to missing data and unbalanced sampling between waterbodies.

A scalable approach for short-term disease forecasting in high spatial resolution areal data.

Short-term disease forecasting at specific discrete spatial resolutions has become a high-impact decision-support tool in health planning. However, when the number of areas is very large obtaining predictions can be computationally intensive or even unfeasible using standard spatiotemporal models. The purpose of this paper is to provide a method for short-term predictions in high-dimensional areal data based on a newly proposed "divide-and-conquer" approach. We assess the predictive performance of this method and other classical spatiotemporal models in a validation study that uses cancer mortality data for the 7907 municipalities of continental Spain. The new proposal outperforms traditional models in terms of mean absolute error, root mean square error, and interval score when forecasting cancer mortality 1, 2, and 3 years ahead. Models are implemented in a fully Bayesian framework using the well-known integrated nested Laplace estimation technique.

Are ride-hailing services safer than taxis? A multivariate spatial approach with accommodation of exposure uncertainty.

Despite many research efforts on ride-hailing services and taxis, limited studies have compared the safety performance of the two modes. A major challenge is the need for reliable mode-specific exposure data to model their safety outcomes. Moreover, crash frequencies of the two modes by injury severities tend to be spatially and inherently correlated. To fully address these issues, this study proposes a novel multivariate conditional autoregressive model considering measurement errors in mode-specific exposures (MVCARME). More specially, a classical measurement error structure accommodates the uncertainty of estimated mode-specific exposures, and a multivariate spatial specification is adopted to capture potential spatial and inherent correlations. The model estimation is accelerated by an integrated nest Laplace approximation method. The census tracts in the city of Chicago are set as the spatial analysis unit. The mode-specific exposures (vehicle-mile-traveled) in each census tract are estimated by trip assignments using ride-hailing and taxi trip data in 2019. The modeling results indicate that both ride-hailing crashes and taxi crashes are positively associated with transportation factors (e.g., vehicle-mile-traveled, mode-specific vehicle-mile-traveled, and traffic signal numbers), land use factors (i.e., number of educational and alcohol-related sites), and demographic factors (e.g., median household income, transit ratio, and walk ratio). By comparison, the proposed model outperforms the others (i.e., negative binomial models and multivariate conditional autoregressive model) by yielding the lowest deviance information criterion (DIC), Watanabe-Akaike information criterion (WAIC), mean absolute error (MAE), and root-mean-square error (RMSE). According to the results of t-tests, ride-hailing services are found to be prone to a higher risk of minor injury crashes compared with taxis, despite no significant difference between the risks of severe injury crashes. Methodologically, this study adds a robust safety evaluation approach for comparing crash risks of different modes to the literature. At the same time, practically, it provides researchers, practitioners, and policy-makers insights into the safety management of various mobility alternatives.

A joint Bayesian framework for missing data and measurement error using integrated nested Laplace approximations.

Measurement error (ME) and missing values in covariates are often unavoidable in disciplines that deal with data, and both problems have separately received considerable attention during the past decades. However, while most researchers are familiar with methods for treating missing data, accounting for ME in covariates of regression models is less common. In addition, ME and missing data are typically treated as two separate problems, despite practical and theoretical similarities. Here, we exploit the fact that missing data in a continuous covariate is an extreme case of classical ME, allowing us to use existing methodology that accounts for ME via a Bayesian framework that employs integrated nested Laplace approximations (INLA) and thus to simultaneously account for both ME and missing data in the same covariate. As a useful by-product, we present an approach to handle missing data in INLA since this corresponds to the special case when no ME is present. In addition, we show how to account for Berkson ME in the same framework. In its broadest generality, the proposed joint Bayesian framework can thus account for Berkson ME, classical ME, and missing data, or any combination of these in the same or different continuous covariates of the family of regression models that are feasible with INLA. The approach is exemplified using both simulated and real data. We provide extensive and fully reproducible Supporting Information with thoroughly documented examples using R-INLA and inlabru.

A Bayesian multivariate spatial approach for illness-death survival models.

Illness-death models are a class of stochastic models inside the multi-state framework. In those models, individuals are allowed to move over time between different states related to illness and death. They are of special interest when working with non-terminal diseases, as they not only consider the competing risk of death but also allow us to study the progression from illness to death. The intensity of each transition can be modelled including both fixed and random effects of covariates. In particular, spatially structured random effects or their multivariate versions can be used to assess spatial differences between regions and among transitions. We propose a Bayesian methodological framework based on an illness-death model with a multivariate Leroux prior for the random effects. We apply this model to a cohort study regarding progression after an osteoporotic hip fracture in elderly patients. From this spatial illness-death model, we assess the geographical variation in risks, cumulative incidences and transition probabilities related to recurrent hip fracture and death. Bayesian inference is done via the integrated nested Laplace approximation.

Mapping the risk of respiratory infections using suburban district areas in a large city in Colombia.

BACKGROUND: Acute respiratory infections (ARI) in Cúcuta -Colombia, have a comparatively high burden of disease associated with high public health costs. However, little is known about the epidemiology of these diseases in the city and its distribution within suburban areas. This study addresses this gap by estimating and mapping the risk of ARI in Cúcuta and identifying the most relevant risk factors. METHODS: A spatial epidemiological analysis was designed to investigate the association of sociodemographic and environmental risk factors with the rate of ambulatory consultations of ARI in urban sections of Cúcuta, 2018. The ARI rate was calculated using a method for spatial estimation of disease rates. A Bayesian spatial model was implemented using the Integrated Nested Laplace Approximation approach and the Besag-York-Mollié specification. The risk of ARI per urban section and the hotspots of higher risk were also estimated and mapped. RESULTS: A higher risk of IRA was found in central, south, north and west areas of Cúcuta after adjusting for sociodemographic and environmental factors, and taking into consideration the spatial distribution of the city's urban sections. An increase of one unit in the percentage of population younger than 15 years; the Index of Multidimensional Poverty and the rate of ARI in the migrant population was associated with a 1.08 (1.06-1.1); 1.04 (1.01-1.08) and 1.25 (1.22-1.27) increase of the ARI rate, respectively. Twenty-four urban sections were identified as hotspots of risk in central, south, north and west areas in Cucuta. CONCLUSION: Sociodemographic factors and their spatial patterns are determinants of acute respiratory infections in Cúcuta. Bayesian spatial hierarchical models can be used to estimate and map the risk of these infections in suburban areas of large cities in Colombia. The methods of this study can be used globally to identify suburban areas and or specific communities at risk to support the implementation of prevention strategies and decision-making in the public and private health sectors.

Assessing survival time of outpatients with cervical cancer: at the university of Gondar referral hospital using the Bayesian approach.

BACKGROUND: Cervical cancer is the 4th most common cancer in women worldwide. as well as the 4th most common cause of cancer-related death. The main objective of this study was to identify factors that affect the survival time of outpatients with cervical cancer. METHODS: A retrolective study including outpatients with cervical cancer was carried out in a hospital. To achieve the aim, 322 outpatients with cervical cancer were included in the study based on the data taken from the medical records of patients enrolled from May 15, 2018, to May 15, 2022, at the University of Gondar referral hospital, Gondar, Ethiopia. The Kaplan-Meier plots and log-rank test were used for the comparison of survival functions; the Cox-PH model and Bayesian parametric survival models were used to analyze the survival times of outpatients with cervical cancer. Integrated nested Laplace approximation methods have been applied. RESULTS: Out of a total of 322 patients, 118 (36.6%) died as outpatients. The estimated median survival time for patients was 42 months. Using model selection criteria, the Bayesian log-normal accelerated failure time model was found to be appropriate. According to the results of this model, oral contraceptive use, HIV, stage, grade, co-morbid disease, history of abortion, weight, histology type, FIGO stage, radiation, chemotherapy, LVSI, metastatic number, regional nodes examined, and tumor size all have a significant impact on the survival time of outpatients with cervical cancer. The Bayesian log-normal accelerated failure time model accurately predicted the survival time of cervical cancer outpatients. CONCLUSIONS: The findings of this study suggested that reductions in weight, treatment, the presence of comorbid disease, the presence of HIV, squamous cell histology type, having a history of abortion, oral contraceptive use, a large tumor size, an increase in the International Federation of Gynecologists and Obstetricians stage, an increase in metastasis number, an increase in grade, positive regional nodes, lymphatic vascular space invasion, and late stages of cancer all shortened the survival time of cervical cancer outpatients.

Bayesian spatio-temporal modelling and mapping of malaria and anaemia among children between 0 and 59 months in Nigeria.

BACKGROUND/M&M: A vital aspect of disease management and policy making lies in the understanding of the universal distribution of diseases. Nevertheless, due to differences all-over host groups and space-time outbreak activities, data are subject to intricacies. Herein, Bayesian spatio-temporal models were proposed to model and map malaria and anaemia risk ratio in space and time as well as to ascertain risk factors related to these diseases and the most endemic states in Nigeria. Parameter estimation was performed by employing the R-integrated nested Laplace approximation (INLA) package and Deviance Information Criteria were applied to select the best model. RESULTS: In malaria, model 7 which basically suggests that previous trend of an event cannot account for future trend i.e., Interaction with one random time effect (random walk) has the least deviance. On the other hand, model 6 assumes that previous event can be used to predict future event i.e., (Interaction with one random time effect (ar1)) gave the least deviance in anaemia. DISCUSSION: For malaria and anaemia, models 7 and 6 were selected to model and map these diseases in Nigeria, because these models have the capacity to receive strength from adjacent states, in a manner that neighbouring states have the same risk. Changes in risk and clustering with a high record of these diseases among states in Nigeria was observed. However, despite these changes, the total risk of malaria and anaemia for 2010 and 2015 was unaffected. CONCLUSION: Notwithstanding the methods applied, this study will be valuable to the advancement of a spatio-temporal approach for analyzing malaria and anaemia risk in Nigeria.

Determining factors associated with cholera disease in Ethiopia using Bayesian hierarchical modeling.

BACKGROUND: Cholera is a diarrheal disease caused by infection of the intestine with the gram-negative bacteria Vibrio cholera. It is caused by the ingestion of food or water and infected all age groups. This study aimed at identifying risk factors associated with cholera disease in Ethiopia using the Bayesian hierarchical model. METHODS: The study was conducted in Ethiopia across regions and this study used secondary data obtained from the Ethiopian public health institute. Latent Gaussian models were used in this study; which is a group of models that contains most statistical models used in practice. The posterior marginal distribution of the Latent Gaussian models with different priors is determined by R-Integrated Nested Laplace Approximation. RESULTS: There were 2790 cholera patients in Ethiopia across the regions. There were 81.61% of patients are survived from cholera outbreak disease and the rest 18.39% have died. There was 39% variation across the region in Ethiopia. Latent Gaussian models including random and fixed effects with standard priors were the best model to fit the data based on deviance. The odds of surviving from cholera outbreak disease for inpatient status are 0.609 times less than the outpatient status. CONCLUSIONS: The authors conclude that the fitted latent Gaussian models indicate the predictor variables; admission status, aged between 15 and 44, another sick person in a family, dehydration status, oral rehydration salt, intravenous, and antibiotics were significantly associated with cholera outbreak disease.

Spatiotemporal Variation and Predictors of Unsuppressed Viral Load among HIV-Positive Men and Women in Rural and Peri-Urban KwaZulu-Natal, South Africa.

Unsuppressed HIV viral load is an important marker of sustained HIV transmission. We investigated the prevalence, predictors, and high-risk areas of unsuppressed HIV viral load among HIV-positive men and women. Unsuppressed HIV viral load was defined as viral load of ≥400 copies/mL. Data from the HIV Incidence District Surveillance System (HIPSS), a longitudinal study undertaken between June 2014 to June 2016 among men and women aged 15−49 years in rural and peri-urban KwaZulu-Natal, South Africa, were analysed. A Bayesian geoadditive regression model which includes a spatial effect for a small enumeration area was applied using an integrated nested Laplace approximation (INLA) function while accounting for unobserved factors, non-linear effects of selected continuous variables, and spatial autocorrelation. The prevalence of unsuppressed HIV viral load was 46.1% [95% CI: 44.3−47.8]. Predictors of unsuppressed HIV viral load were incomplete high school education, being away from home for more than a month, alcohol consumption, no prior knowledge of HIV status, not ever tested for HIV, not on antiretroviral therapy (ART), on tuberculosis (TB) medication, having two or more sexual partners in the last 12 months, and having a CD4 cell count of <350 cells/µL. A positive non-linear effect of age, household size, and the number of lifetime HIV tests was identified. The higher-risk pattern of unsuppressed HIV viral load occurred in the northwest and northeast of the study area. Identifying predictors of unsuppressed viral load in a localized geographic area and information from spatial risk maps are important for targeted prevention and treatment programs to reduce the transmission of HIV.

Past, present, and future trends of overweight and obesity in Belgium using Bayesian age-period-cohort models.

BACKGROUND: Overweight and obesity are one of the most significant risk factors of the twenty-first century related to an increased risk in the occurrence of non-communicable diseases and associated increased healthcare costs. To estimate the future impact of overweight, the current study aimed to project the prevalence of overweight and obesity to the year 2030 in Belgium using a Bayesian age-period-cohort (APC) model, supporting policy planning. METHODS: Height and weight of 58,369 adults aged 18+ years, collected in six consecutive cross-sectional health interview surveys between 1997 and 2018, were evaluated. Criteria used for overweight and obesity were defined as body mass index (BMI) ≥ 25, and BMI ≥ 30. Past trends and projections were estimated with a Bayesian hierarchical APC model. RESULTS: The prevalence of overweight and obesity has increased between 1997 and 2018 in both men and women, whereby the highest prevalence was observed in the middle-aged group. It is likely that a further increase in the prevalence of obesity will be seen by 2030 with a probability of 84.1% for an increase in cases among men and 56.0% for an increase in cases among women. For overweight, it is likely to see an increase in cases in women (57.4%), while a steady state in cases among men is likely. A prevalence of 52.3% [21.2%; 83.2%] for overweight, and 27.6% [9.9%; 57.4%] for obesity will likely be achieved in 2030 among men. Among women, a prevalence of 49,1% [7,3%; 90,9%] for overweight, and 17,2% [2,5%; 61,8%] for obesity is most likely. CONCLUSIONS: Our projections show that the WHO target to halt obesity by 2025 will most likely not be achieved. There is an urgent necessity for policy makers to implement effective prevent policies and other strategies in people who are at risk for developing overweight and/or obesity.