ALL Metrics
-
Views
-
Downloads
Get PDF
Get XML
Cite
Export
Track
Research Article
Revised

Spatio-temporal characteristics of Tuberculosis in Ghana

[version 3; peer review: 1 approved, 1 approved with reservations]
PUBLISHED 11 Mar 2024
Author details Author details
OPEN PEER REVIEW
REVIEWER STATUS

Abstract

Background

Tuberculosis (TB) continues to be a prominent contributor to global mortality, standing as the second most fatal infectious disease and holding the seventh position among the top ten causes of death in Ghana. There is insufficient literature regarding the utilization of Bayesian hierarchical models, specifically within the framework of Integrated Nested Laplace Approximation (INLA), for examining the spatial and spatio-temporal dynamics of tuberculosis risk in Ghana. This study addresses this gap by determining TB hotspots regions in Ghana using the Bayesian modeling framework within the INLA.

Methods

TB data were sourced from the Ghana Health Service and National Tuberculosis Programme for the 10 administrative regions of Ghana, from 2008 to 2017. The relative risk of TB for each region and year was estimated utilizing Bayesian spatial and spatio-temporal modeling frameworks. Baseline predictors of TB risk were also considered. Maps for TB risks were created to visualized regions with TB hotspots. Model fitting and parameter estimation were conducted using R version 4.3.2.

Results

Among the baseline predictors, factors such as TB cure rate, TB success rate, knowledge about TB, HIV prevalence, percentage of literacy, and high income were found to be most significant in influencing the TB risk across the ten regions in Ghana. We noted an increased risk of TB infection in the Northern zone and the Eastern and Greater Accra regions in the Southern zone. Spatio-temporal distribution of TB infection risk was predominantly concentrated in the Southern zone. Clustering of TB risk was observed among neighboring regions.

Conclusion

To achieve a significant reduction in TB cases, it is essential to allocate resources to TB hotspots regions and also implement measures to control significant predictors of TB infection risk.

Keywords

Bayesian spatial and space-time models, Tuberculosis relative risk, baseline predictors and TB hot-spots.

Revised Amendments from Version 2

Based on the first reviewer’s comments, the following revision were made on version 1: authors proof-read the entire manuscript to address grammatical and spelling mistakes (comment 1), to reflect text/sentences suggestions (comments 3, 7,8 and 9), formatted some references (comment 4).  Thought provoking comments/questions raised by the second reviewer compelled authors to revise the entire document (version 2). Example, in version 2, authors have explicitly stated the research gap the study seeks to address (abstract and introduction section), revised the methodology for easy understanding and highlighted the robustness of the methodology to issues of parameter convergence, some criticisms, and benefits of the methods. In the data description section, we have also discussed the potential limitations or biases in this data collection process.  Authors have thoroughly discussed the study findings comparing the findings to existing literature. The revision also included policy direction and provided recommendations for future research (conclusion section). Each section was almost entirely revised, guided by the reviewer’s general comments 1.1 and 1.2.  In the course of revision, we have to re-run the models to clarify some unusual values of parameter estimates as pointed out by the second reviewer. Example, “dots” are added in Figures 1-3 (now Figures 2-4) to indicate values for years and regions. The results without covariate adjustments were redundant, hence were removed. The estimates (in the first row of Table 2, now Table 1) were unusual and were revised after re-analyses. Figure 7 (now Figure 6) was then produced after getting the results in Table 1.  We also re-run the spatio-temporal models to revise Table 5 (now Table 2), Table 6 (now in Table 3) and Figures were also revised accordingly. Statistical inferences in version 2 remain the same despite this extensive revision.

See the authors' detailed response to the review by Dr. Victor Abiola Adepoju
See the authors' detailed response to the review by Dioggban Jakperik

Background

Tuberculosis (TB) is an infectious disease caused by the bacillus Mycobacterium tuberculosis.1 It is a condition that transcends age, gender, race, and the health status of individuals worldwide.1 In 2019, around 10 million people tested positive for TB, resulting in an estimated 1.2 million deaths among those without HIV infection.1 Among adults aged 15 and above, males comprised 56% of global infections, while females accounted for 32% in the same age category. Children represented 12% of total cases, with approximately 8.2% of reported cases involving individuals with HIV.1 The distribution of TB infections in 2019 varied across regions, with Africa contributing to 25% of cases, South-East Asia and the Western Pacific accounting for 44% and 18%, respectively. In the same year, the Eastern Mediterranean, America, and Europe reported 8.2%, 2.9%, and 2.5% of the total cases, respectively.1

Nevertheless, the World Health Organization (WHO) has advocated for the development of tuberculosis (TB) vaccines.1 In addition to vaccines, there are available treatments for TB patients, and reports suggest that approximately 85% of patients can be successfully cured with a six-month drug regimen.1 According to available data, treatment drugs have played a crucial role in preventing over 60 million deaths from 2000 to 2020.1 Collectively, there was a 9% reduction in cases between 2015 and 2019, and a further 2.3% reduction between 2018 and 2019.1 Europe has successfully achieved a 19% reduction in cases and a 31% reduction in deaths between 2015 and 2019. In the same timeframe, Africa has seen a 16% reduction in cases and a 19% reduction in deaths. More comprehensive statistics on TB can be found in references.16

While global tuberculosis (TB) reports show a declining trend in cases and deaths, the World Health Organization’s (WHO) targets for reduction from 2015 to 2020 have not been met.1 Except for Europe, all other continents have fallen short of achieving the prescribed reduction levels.1 Consequently, TB remains among the top 10 causes of death, particularly in Africa and Ghana, where it outranks HIV/AIDS.1

Ghana has been affected by the respiratory disease and currently has challenges in eradicating TB. In response, the country implemented policies such as Directly Observed Therapy (DOT) and National Tuberculosis Programmes (NTPs) in 1994 to detect and treat TB.1,7,8 The initiation of the NTP resulted in 100% DOTs coverage by 2005, leading to an increased detection of TB cases for treatment each subsequent year.9,10 For example, the number of TB cases detected rose from 7,425 in 1996 to 15,286 in 2009.9,10 A more in-depth exploration of TB statistics in Ghana can be found in references.1,713 The incidence of tuberculosis (TB) has demonstrated a consistent decline, with figures decreasing from 45,000 cases in 2010 to 44,000 in 2017. Annually, approximately 44,000 individuals develop TB in Ghana, with an estimated 6,600 being children. The number of TB patients undergoing treatment ranged between 14,607 and 15,389 from 2010 to 2017, highlighting a pressing demand for increased resources dedicated to TB patient care. In 2013, the national TB prevalence was estimated at 290 cases per 100,000 population, indicating a burden fourfold higher than the World Health Organization’s (WHO) estimates for the same year (71 cases per 100,000 population). TB-related mortality is estimated at 10,000 deaths annually. Tuberculosis is the seventh leading cause of mortality in Ghana, contributing to 4.9% of all deaths and older males, aged above 45 years, bear the greatest burden of the disease. Also, TB ranks 7th among the leading causes of death for females and 8th for males, and it stands as the 4th primary cause of mortality within the category of communicable, maternal, neonatal, and nutritional diseases. The TB epidemic exhibits a generalized pattern, albeit with geographical variations in case notification, often associated with disparities in access to healthcare facilities.

Despite a decline in tuberculosis (TB) cases and deaths attributed to the implementation of mitigation and treatment strategies, TB remains a life-threatening disease, placing a substantial burden on the health infrastructure in Ghana. Consequently, TB has become a focal point of research across diverse fields. Notably, researchers from various backgrounds, such as authors,7,8,14 have investigated the dynamics of TB indicators and risk factors associated with TB in Ghana. Osei and colleagues,7,8 studied the trends of TB detection and treatment outcomes, employing logistic regression to analyze the relationship between patient and disease characteristics. In a separate study, Aryee and colleagues14 examined the dynamics of TB using Autoregressive Moving Average (ARIMA) methods, analyzing TB data recorded by Korle Bu Teaching Hospital from 2008 to 2017. On the other hand, Iddrisu and colleagues2 investigated the temporal and geographical patterns of TB prevalence in Ghana between 2015 and 2018.

The existing body of literature on the utilization of the convolution model or the Besag, York, Mollié (BYM) model within the Bayesian modeling framework, particularly in the application of spatial and spatio-temporal models for estimating and mapping disease risk in Ghana, is limited. This paper addresses this gap by identifying regions with hotspots and predictors of TB infection risk in Ghana using these methods. We primarily focused on studying the spatial and spatio-temporal patterns of TB relative risk using spatial and spatio-temporal models within the Bayesian modeling framework.5,1521 These models enable the incorporation of neighboring regions’ information to yield accurate estimates of disease risk within and between regions. The key hypotheses guiding this study are as follows: (1) TB infection risk exhibits clustering among specific regions in Ghana, (2) there is variation in TB infection risk among different regions in Ghana, and (3) TB preventive measures, education on TB, and economic variables serve as significant predictors of TB infection risk in Ghana. Parameter estimates are derived from the marginal posterior distribution using the Integrated Nested Laplace Approach (INLA) implemented through R software version 4.3.2.2023

Methods

This section provides a description of the data used in this study and specifies the Bayesian spatial and spatio-temporal models that were employed to estimate the relative risk of TB across the 10 old administrative regions of Ghana. The conditional auto-regressive (CAR) model and the random effect components of the Besag, York and Mollié (BYM) model for respectively identifying clustering and heterogeneity of disease risk in space and time are described.2,3,19,2430 The spatial model used is the BYM Model24,31 and the spatio-temporal models are based on three modeling frameworks20,21,32 developed by Knorr-Held and Ra’ser,33 Bernardinelli and colleagues34 and Waller and colleagues.35 The best fitting spatio-temporal model was selected and only results from this model were discussed.

TB cases

For this study, we utilized tuberculosis (TB) detection data sourced from the Ghana Health Service and National Tuberculosis Programme.9,13 The dataset comprises information on TB cases spanning the years 2008 to 2017, encompassing the 10 old administrative regions of Ghana. These regions consist of Ashanti, Brong Ahafo, Central, Eastern, Greater Accra, Northern, Upper East, Upper West, Volta, and Western, as illustrated in Figure 1. Although new regions were established after 2016, the TB data were collected and organized based on the original 10 administrative regions. One significant limitation is the reliance on existing data, which could contain errors if healthcare workers input information incorrectly. Additionally, using secondary data may compromise external validity if biases exist in the population under study. Another potential source of bias is non-random participant selection, which could make the estimates unreliable if the sample doesn’t accurately represent the entire population.

6ccd56a8-2b56-4f5c-a1e5-cd91dd262082_figure1.gif

Figure 1. Map representing 10 administrative regions in Ghana.

Outcome variable

The dependent variable in this study is the number of tuberculosis (TB) cases collected for each of the 10 administrative regions in Ghana, spanning the years 2008 to 2017. Figure 2 illustrates the TB trends in these regions over the specified time frame. Generally, there is an overall decline in TB cases observed across all regions of Ghana from 2008 to 2016, with the exception of the Brong Ahafo Region where TB cases increased. Notably, between 2016 and 2017, there was a significant increase in TB cases in the Northern and Upper East regions. Conversely, the Ashanti Region witnessed a decrease from 50 per 100,000 population in 2016 to 45 per 100,000 population in 2017. Specifically, TB cases in the Northern Region surged from 24 per 100,000 population in 2016 to 52 per 100,000 population in 2017. Similarly, the Upper East Region experienced an increase from 53 per 100,000 population in 2016 to 63 per 100,000 population in 2017. However, the changes in TB cases in the remaining regions remained relatively stable during this period.

6ccd56a8-2b56-4f5c-a1e5-cd91dd262082_figure2.gif

Figure 2. Trend of TB cases detection per 100,000 by region for 10 years from 2008 to 2017 in Ghana.

Figure 3 depicts the trajectory of the total number of tuberculosis (TB) cases for each region from 2008 to 2017. Notably, it illustrates that the Greater Accra Region consistently recorded the highest number of cases during this period, while the Northern Region registered the lowest. The Volta and Western Regions ranked second and third, respectively, with slightly lower case counts than the Greater Accra Region. Furthermore, Figure 4 displays the trend of total TB cases for each year or period. The data reveals that the peak in total TB cases occurred in 2011, whereas the lowest count was recorded in 2016. There is a discernible slow decline in TB cases from 2011 to 2017, marked by slight increases in 2012 and 2017.

6ccd56a8-2b56-4f5c-a1e5-cd91dd262082_figure3.gif

Figure 3. Trend of the total TB cases detection per 100,000 for each region from 2008 to 2017 in Ghana.

6ccd56a8-2b56-4f5c-a1e5-cd91dd262082_figure4.gif

Figure 4. Trend of the total TB cases detection per 100,000 population for each year.

Moreover, the variability of tuberculosis (TB) cases from 2008 to 2017 is illustrated through box-and-whisker plots in Figure 5. The overlapping box-and-whisker plots in Figure 5A suggest a lack of variability in TB cases across the years. These plots indicate a skew towards larger numbers from 2008 to 2014, shifting towards smaller numbers from 2015 to 2016. Particularly small numbers were noted in 2013, 2014, 2015, and 2017. Figure 5B presents the variability across regions. The non-overlapping box-and-whisker plots signify variability in TB cases among the regions. Additionally, it is observed that TB cases in most regions are skewed towards larger numbers, except for the Northern and Upper East regions.

6ccd56a8-2b56-4f5c-a1e5-cd91dd262082_figure5.gif

Figure 5. Box-and-whisker plot of TB cases by year (Figure A) and by region (Figure B).

Baseline predictors

This study considered various baseline regional characteristics that may influence the risk of tuberculosis (TB) infection. The set of baseline predictors encompasses factors such as the doctor-to-population ratio, nurse-to-population ratio, HIV prevalence, tuberculosis cure rate, tuberculosis success rate, wealth quantiles, and the proportions of men and women who are employed, unemployed, educated, and uneducated. The study also considered variables such as the proportions of individuals who have heard about TB, possess knowledge that TB is airborne, are aware that TB can be cured, and believe that TB status should be kept confidential. In the data analyses, all these baseline variables are thoroughly examined to identify significant predictors of TB cases.

Statistical methods

In this section, we discussed and defined the spatial and spatio-temporal models employed for estimating tuberculosis (TB) risk in Ghana. Considering that these models are within the Bayesian modeling framework, we initially provided a concise overview of Bayesian principles.

Bayesian modeling framework

For spatial data collected at a single time point, where yi,i=1,,n represents the observed TB cases for region i, consider a Poisson random variable with a probability mass function defined as Pryi|θ, where θ=θ1θn represents a vector indicating the TB risk for each region. The likelihood function for the Poisson variable yi is defined as Pry|θ=i=1nPryi|θ, assuming that the sample values of the vector y=y1yn given the parameter estimates θ are independent.19 In a Bayesian modeling framework, prior distributions of the unknown parameters θ in the likelihood function are required. The prior distribution reflects the existing knowledge about the parameters θ before observing the data yi.19 Within the Bayesian framework, all parameters are stochastic and are assigned appropriate distributions known as prior distributions.19

The Bayesian modeling framework integrates the likelihood function for the data with prior distributions for parameter estimates, giving rise to a probability distribution known as the posterior distribution.2,3,5,1921 This distribution, denoted as Prθ|y, signifies the probability distribution of the parameters given the observed data. It is proportional to the product of the likelihood function and the prior Prθdistributions, as defined by:

(1)
Prθy=PryθPrθpLyθPrθdθ

In equation (1), the denominated is often referred to as the normalizing constant and has been shown to integrate to 1. The implication is that the equation (1) can be alternatively expressed as Prθ|y Pry|θPrθ. The Bayesian modeling framework assumes that the parameter estimates are stochastic and hence assumes a distribution for each parameter estimate, including hyperparameters.3638 The main criticism of such an assumption is that it depends on subjective prior assumptions, which might result in extreme or confusing conclusions. However, these priors also represent an advantage since they provide a tool for statisticians to be transparent about their assumptions.3638 Overweighting or informative priors have the potential to mislead statistical inferences about the data being studied since parameter estimates would be much influenced by such priors rather than the data under study. In Bayesian modeling, conjugate or non-informative priors are recommended in order to allow the data to have more influence on the parameter estimates.3638

Using traceplots or density plots of parameter estimates,39 Markov Chain Monte Carlo (MCMC) methods determine when a model’s parameters have converged, greatly facilitating Bayesian inference for complex models. Nonetheless, these methods are computationally intensive and encounter various convergence issues. Integrated nested Laplace approximation (INLA) offers a computationally efficient alternative to MCMC for approximate Bayesian inference in latent Gaussian models.40 These models encompass a broad range, from generalized linear mixed models to spatial and spatio-temporal models.20,40 INLA employs a combination of analytical approximations and numerical algorithms for sparse matrices to approximate posterior distributions using closed-form expressions. This enables faster inference, avoid issues of sample convergence and mixing, thereby facilitating the analysis of large datasets and the exploration of alternative models. In this study, we obtained parameter estimates using INLA.40 INLA, facilitated by the R-INLA package, has been widely embraced by various researchers.20,21,41,42 Also, this package can directly be installed in R software by using the command install.packages(“INLA”,repos=c (getOption(“repos”),INLA=“https://inla.r-inla-download.org/R/stable”), dep=TRUE).

Spatial model: Besag, York and Mollié (BYM) Model

This research utilizes the Besag, York, and Mollié (BYM) Model as its spatial modeling approach.3,1921,24,31,43 The model consists of two main components: the conditional autoregressive (CAR) model, which captures and conveys information about the clustering of disease risk among neighboring regions, and the random effect component, which identifies heterogeneity or variability in disease risk among regions. The BYM model is defined asηi=μi+ϕi+ui, where μirepresents the mean, ui embodies the random effect structure, and ϕisignifies the CAR structure. It is assumed that the random effect component, ui, follows a Gaussian distribution.

A compelling feature of the Conditional Autoregressive (CAR) model is its capability to create a framework that facilitates the sharing of strength among neighboring regions. As a result, adjacent regions are inclined to display similar risks, whereas those that are more distant are more likely to exhibit variability in terms of the risk associated with the studied disease.44 In this investigation, we employed distances or boundary information between regions to define neighborhood properties within the CAR model.18,19,27

Let φ=i=12n represent the regions, and Ni=jφ:ij identifies regions that share boundaries with region i. Additionally, let ϕi, where iφ, have a stochastic distribution. It is assumed that iΩ follows a Gaussian distribution, as outlined by:

(2)
ϕi|ϕjiNijWijϕjτϕ2,

where Wij represents the spatial dependence parameter determining the impact of each observation on the Conditional Autoregressive (CAR) structure ϕi,where ϕj encompasses all observations except ϕi, and τϕ2 signifies the precision of ϕi. The spatial dependence parameter Wiiis non-zero if jS. To prevent autocorrelation, we designate Wij= 0. This is to avoid conducting regression involving any observation against itself.2426,35,44 The model specification (2) indicates that ϕi solely relies on a set of neighboring ϕj when the region j is part of the neighborhood set Ni of ϕi. For an in-depth explanation of the CAR structure, refer to.2426,35,44

Given that yi follows a Poisson distribution with parameter Eiexpηi, and μi=Eiexpηi, the Besag, York, and Mollié (BYM) model with covariates can be described as ηi=Xβ+ϕi+ui. Here, Ei=rNi represents the expected number of TB cases in region i, with r=t=1nyii=1nNi, where Ni denotes the at-risk population in region i. Furthermore, X denotes a design matrix of covariates, and β represents a vector of parameter estimates associated with these covariates.

Given θi=expηi, it implies that logθi=ηi, where ηi represents the linear predictor. Consequently, the relative risk for each region i can be expressed as θi=expXβ+ϕi+ui, and the log-log link function can be defined as logμi=logEi+expXβ+ϕi+ui.

To estimate parameters from the Bayesian posterior distribution, it is necessary to have the likelihood function of the data, along with prior distributions for the Conditional Autoregressive (CAR) structure, random effect, and the regression coefficients β. The prior distributions for the random effects and regression coefficients are assumed to conform to a Gaussian distribution. Consequently, the resulting posterior distribution can be specified as follows:

Prθβuϕτβ2τu2τϕ2|yEPryE|θβuϕτβ2τu2τϕ2PrβPruPrϕ.

The hyper-prior precision parameters τu2,τϕ2 and τβ2 follow the Gamma distributions: τu2 Gamma 0.5,0.005,τϕ2Gamma0.5,0.005 and τβ2Gamma0.5,0.01 repectively. The vector of regression coefficients, β, is governed by a Gaussian distribution specified as βN0τβ2. The estimated parameters, τu2 and τϕ2, represent the precision variance parameters for u and ϕ respectively. They are utilized to gauge the level of variability among regions and the clustering of risk among neighboring regions.19 Additionally, we empirically estimated the spatial structure effect as τϕ2=i=1nϕiϕ¯2n1, where ϕ¯ is the average of ϕ. This empirical estimation was then compared to the posterior marginal variance for the unstructured effect, provided by τu2 , provided proportion spatial= τϕ2τϕ2+τu2. We utilized Prθi>1|yi to identify regions with a high probability of demonstrating a risk greater than 1.

Spatio-temporal models

The TB case data utilized in this study were gathered over a period, and consequently, relying solely on spatial models may not suffice to capture and elucidate the spatio-temporal pattern of the relative risk of TB. This limitation arises because spatial models are confined to identifying heterogeneity and clustering of disease risk at a specific time point. In this section, we introduce three spatio-temporal models20 designed to investigate the spatio-temporal pattern of disease risk. These models vary in their space-time interaction structures.20,21,33,35,45

Let yij denotes observed TB cases in region i, year j. The distribution of TB cases is defined as yijPoissonEijexpηij, where the unknown relative risk θij for region i at year j can be expressed as θij=expηij and Eij is the expected number of TB cases in region i at year j. We defined the crude rate of TB cases for region i at year j as rijs=yijsNis and the expected TB cases in region i at year j is defined as Eij=ritsNij=yijsNijsNij, where Nij denotes the observed population, yijs is the TB cases in the standard population. This means that the overall crude rate of TB is given by r=injJyijsNij3 and the overall number of expected TB cases is defined by E=iNjTrijsNij=iNjJyijsNijsNij.3,2

The first spatio-temporal model34 presented in this study defines the linear predictor ηit as

(3)
ηij=μ+ϕi+ui+ϱ+δi×j,

where ϕi+ui component follow the BYM24 model’s specification, ϱ×j is the global linear time trend, j×δi is the interaction term between space and time defining the difference between ϱ and the area-specific time trend.20,46 The variable j represents a vector of temporal weights, and the intercept μ quantifies the average TB rate across all 10 regions. The differential trend δi signifies the interaction between time and space. The relative risk is expressed as θij=expηij, and the logarithm of the risk is logθij=ηij.This implies that θij=expηij=expμ+ϕi+ui+ϱ+δi×j. Consequently, the Poisson mean is μij=Eijexpμ+ϕi+ui+ϱ+δi×j, and the logarithm of the mean is defined as logμij=logEij+μ+ϕi+ui+ϱ+δi×j. These formulations suggest that each spatial unit possesses its own time trend with a spatial intercept μ+ϕi+ui and a slope ϱ+δi. This model assumes a linear time trend in each spatial unit. The parameters to be estimated are φ=ϱϕuδ, and the hyper-parameters are denoted as ψ=τϕτuτδ.

When accounting for TB risk factors Xi, model (3) can be expressed as model (4). The parameter estimates vector is now denoted as φ=βϱϕuδ, and the hyper-parameters are represented by ψ=τϕτuτβτδ

(4)
ηij=α+p=1PβpXi+ϕi+ui+ϱ+δi×j

It is established that if δi<0, the region-specific trend is less pronounced than the mean trend. Conversely, δi>0 indicates that the region-specific trend is steeper than the mean trend.20,46 The parameter δi is assumed as δiNormal0τδ.

The second spatio-temporal model35 defines the linear predictor as

(5)
ηij=α+ϕi+ui+ϑj+ωj

where the combination of terms ϕi+ui represents the BYM model, and the structures ϑj and ωj signify the temporally structured and random effects, respectively. This model assumes a non-parametric time trend. By incorporating covariates, the model (5) transforms into:

(6)
ηij=α+βiXi+ϕi+ui+ϑj+ωj

The parameters to be estimated encompass φ=αβϕuϑω, with corresponding hyper-parameters denoted as ψ=τϕτuτϑτω. The temporally structured effect, ϑj, is captured and modeled using a random walk through a neighboring structure.20 Specifically, for j=1, ϑj|ϑjNϑt+1τϑ; for j=2,,J1, ϑj|ϑjNϑj1+ϑj+12τθ2; and for j=J, ϑj|ϑjNϑj1τϑ. Finally, ωjis specified using a Gaussian exchangeable prior: ωjNormal0τω.

The third spatio-temporal model (7) proposed by33,47 is an extension of model (6) to allow for a space and time interaction πij to explain the difference in time trend of TB cases for various regions.

The third spatio-temporal model (7), proposed by Knorr-Held,33,47 extends model (6) by introducing a space and time interaction πij. This interaction term is included to account for variations in the time trend of TB cases across different regions.

(7)
ηij=μ+ϕi+ui+ϑj+ωj+πij

In this model, we estimate φ=μϕuϑωπ and ψ=τϕτuτϑτωτπ, where πij represents the interaction between ϕi and ui. It assumes no interaction between ϕi and ϑt; hence, πijN0τπ. By incorporating covariates into the model (7), we obtain model (8) and subsequently need to estimate θ=μβϕuϑωπ and ψ=τϕτuτϑτωτπ

(8)
ηij=μ+βiXi+ϕi+ui+ϑj+ωj+πij

Concerning the interaction term πij, if we assume the presence of spatial or temporal structure, then πijN0τδ.33 In this study, all precision parameters are assumed to follow the gamma distribution.20

Results

In this section, we applied the spatial and spatio-temporal models to the TB data. These models were implemented in R version 4.3.2 using the Integrated Nested Laplace Approximation (INLA) method.20,42 A comparison of the performance of the spatio-temporal models was conducted, and the best-performing model was selected based on the Deviance Information Criterion (DIC) developed by Refs. 20, 42, 48. The model with the lowest DIC was considered the most effective. We then presented and discussed the results obtained from both the BYM model and the best-fitting spatio-temporal model.

Results from the BYM model

The results from the BYM model with the posterior estimates are shown in Table 1. Among the baseline predictors, the significant predictors for TB infection risk in Ghana that yield accurate models include: Tuberculosis cure rate β1, Tuberculosis success rate β2, proportion of people with knowledge about TB β3, proportion of those who know that TB is airborne β4, HIV prevalence β5, proportion in high income group β6 and literacy β7.

Table 1. Summary statistics: posterior mean, standard deviation (Sd) and 95% credible interval for the fixed and random effects of the BYM model.

EstimateSd25%50%95%
Fixed effects
μ9.0855.543.706.11526.537
β11.0811.0271.0241.0811.141
β20.8551.0400.7890.8550.927
β31.0461.0171.0101.0461.083
β40.8971.0550.8060.8971.000
β50.4501.2970.2660.4500.762
β60.9461.0170.9140.9460.979
β71.1161.0451.0221.1161.219
Random effects
τu25.8020.623.8320.4180.01
τϕ1834.881810.43121.581299.346656.65

Figure 6A displays the map depicting the posterior mean for the region-specific relative risk θi of TB infection, while Figure 6B presents a map indicating the excess risk of TB. It can be observed Figure 6A that high risk of TB infection is distributed towards the Northern part of Ghana ranging between 1-1.2 for Upper East, Upper West and Northern regions. High risk of TB infection can be observed in Eastern and Greater Accra regions in the Southern part of Ghana. It can be observed in Figure 6B that Greater Accra region is the region with the highest probability of risk (0.8-1) exceeding the national risk 1 followed by Upper West, Northern and Eastern regions with risk between 0.5-0.8 and lowest risks (0.2-0.5) are observed in Upper East, Brong-Ahafo, Western, Volta, Ashanti, and Central regions.

6ccd56a8-2b56-4f5c-a1e5-cd91dd262082_figure6.gif

Figure 6. BYM model: Maps showing the distribution of the region-specific TB risk (Figure A) and distribution of region-specific exceedance probability (Figure B).

Once more, there is slight variability in TB infection risk was noted across various regions, alongside significant clustering or similarities of TB infection risk among neighboring regions shown in Figure 6.

The results in Table 1 confirm the similarity or clustering of risk between or among neighboring regions. This is indicated by the low variability captured by the precision of the spatial structure τϕ. The precision τu of the unstructured component of the BYM model indicates that risk is slightly heterogeneous among regions. The exponent of the posterior mean μ (overall mean effect) indicates that there is approximately 9-fold reduction in TB infection rate across the 10 regions in Ghana.

From Table 1, It can be observed that TB cure rate β1 increases the risk of TB infection by approximately 8%. This observation implies that as more cases are detected, more cases are cured and hence TB cases will in general decrease over time. This explains why TB success rate β2 leads to 14% reduced risk of TB infection. The results also revealed that knowledge about TB β3 significantly increases TB infection risk by approximately 5%. This behavior is expected because, as people become aware of TB, preventive measures are taken. The proportion of high-income group β6 is associated with 5% reduction in TB infection while proportion of literacy β7 is associated with 12% increase in TB cases. High income increases the use of health facilities and testing for TB, thus, leading to a reduction of TB cases. HIV prevalence β5 in the region leads to a 55% reduction in TB cases.

Best fitting spatio-temporal model for TB infection risk

In this section, we fitted the three spatio-temporal models, adjusting for covariates effects on TB risk. Let us call the model (4); Model I, model (6); Model II and model (8); Model III. The best model is then selected and results from such a model are reported. Table 2 presents the DIC, mean deviance D¯ and effective number of parameters pD components for the three space-time models. These indicators of a model’s performance suggest that Model I proved to be the best fitting model among the three candidates’ spatio-temporal models. Hence, we presented only the results from this model.

Table 2. Indicators of space-time models performance.

ModelD¯pDDIC
Model I518.617.81536.42
Model II547.110.56557.62
Model III546.811.00557.76

Results from spatio-temporal model: Model I

The results in Table 3 showed that TB infection risk decreases across the 10 regions of Ghana. However, there is no significant reduction in TB infection risk over the study period. The precision parameter τu indicates very low variability in the risk of TB detection among the regions and much clustering of risk between neighboring regions exhibited by high precision parameter τϕ values for the spatial structure. High precision characterized by τδ indicates lower variability associated with δi. This further indicates that there is no significant interaction between space and time as well as global trend ϱ and area-specific trend δi. Hence, the area-specific trend δi is less remarked than the mean trend.

Table 3. Summary statistics: posterior mean, standard deviation (sd) and 95% credible interval for the fixed and random effects of Model I.

Estimatesd25%50%95% CI
Fixed effects
μ0.1771.0750.0200.1811.400
t1.0061.0070.9921.0051.019
β10.9201.0120.8990.9200.943
β21.1141.0181.0761.1141.155
β30.9781.0080.9630.9780.993
β41.2451.0231.1891.2451.302
β52.2731.1201.8112.2772.829
β61.0381.0081.0211.0381.052
β70.8581.0190.8280.8580.891
Random effects
τu835.681220.692.33337.684258.23
τϕ1272.831543.4624.97738.525456.29
τδ521.70400.33114.09412.561581.45

The results also revealed that TB success rate significantly increases TB cases by 11%. Also, knowledge about TB significantly reduces TB cases by approximately 2%, while increasing TB cure rate, significantly reduces detection by 8%. Awareness that TB is airborne increases TB detection by approximately 25%. That is, more people are willing to participate in TB testing to know their status leading to more case detection. We also observed that HIV prevalence and high income significantly increases TB detection by 27% and approximately 4%, respectively. Literacy significantly reduces the risk of TB detection by approximately 14%.

The TB infection risk θi representing the spatial component of the spatio-temporal model is display in Figure 7A and Figure 7B is the map of exceedance risk. Both Figure 7A and Figure 7B showed that TB infection risk is clustered among regions in the Southern part of Ghana with Brong-Ahafo, Volta, Ashanti, Eastern, Greater Accra regions showing high risk (1-1.1) and the rest of the region show low risk between 0.9-1. These observations account for the low variability captured by both the unstructured and structured components of the area-specific trend. Furthermore, the precision τϑ of the random effect of the temporal structure ϑj suggests a clustering of the relative risk of TB over time. Maps showing spatio-temporal pattern of TB infection risk for 2008 and 2010 are respectively presented in Figure 8A and Figure 8B for illustration. Elevated TB infection was found in the Southern part of the country across all the years. However, the infection risk was consistently below 1 in all the years.

6ccd56a8-2b56-4f5c-a1e5-cd91dd262082_figure7.gif

Figure 7. Spatio-temporal model: regions-specific risk (Figure A) and exceedance probability (Figure B).

6ccd56a8-2b56-4f5c-a1e5-cd91dd262082_figure8.gif

Figure 8. Spatio-temporal pattern of TB infection risk for 2008 (Figure A) and 2010 (Figure B).

Discussion

Accurate assessment and mapping of a disease risks are crucial for the efficient allocation of resources and targeted interventions to address disease hotspots, where there is a high risk of the disease. This research provides valuable insights into tuberculosis (TB) infection hotspots across various regions in Ghana by employing maps that highlight areas with elevated TB infection risk. The study investigates the spatial and spatio-temporal patterns of TB prevalence using a Bayesian modeling framework. Within this framework, we explore the clustering and heterogeneity of TB relative risk, utilizing a conditional auto-regressive (CAR) spatial structure and an unstructured component for random effects24 respectively. The data utilized in our analysis are sourced from the Ghana Health Service and the National Tuberculosis Program. To fit the models, we employed the Integrated Nested Laplace Approximation through the R software.20,21,40

In this study, the model utilized to explore the spatial pattern of TB risk was based on the convolution model, specifically the Besag-York-Mollié (BYM) model.24 The findings indicate a heightened risk of tuberculosis (TB) infection in the Northern part of Ghana, specifically in the Upper East, Upper West, and Northern regions. Similarly, the Southern part, particularly the Eastern and Greater Accra regions, exhibits a high risk of TB infection. Notably, the Greater Accra region stands out with the highest probability of risk, surpassing the national risk level of 1, a trend consistent with the observation that more TB cases are reported in the Greater Accra regions, as noted by various authors.49 Additionally, the study reveals significant clustering or similarities in the risk of tuberculosis (TB) infection among adjacent regions in Ghana. Moreover, there is an observed decline in the rate of TB infections across the ten regions in Ghana, as observed in Refs. 13, 50. The clustering of disease risk between regions may stem from shared geographic factors or possibly be linked to characteristics such as demographics, clinical aspects, or epidemiological factors.51

The study findings indicate that the TB cure rate is associated with an increase in the risk of TB infection. This suggests that as more cases are detected and cured, the overall number of TB cases tends to decrease over time. Conversely, a higher TB success rate is linked to a reduction in the risk of TB infection, highlighting the effectiveness of successful treatment in lowering the prevalence of TB. Furthermore, increased knowledge about TB is shown to significantly raise TB infection risk; as more people are willing and available to be tested for TB, more TB infected cases are detected. On the socioeconomic front, a higher proportion of high-income individuals is associated with a reduction in TB infection, as it correlates with increased use of health facilities and testing for TB. However, a higher proportion of literacy is linked to an increase in TB cases, suggesting a complex relationship between literacy levels and TB infection. Notably, HIV prevalence is found to be associated with a substantial reduction in TB cases, underscoring the impact of HIV prevalence on TB infection rates. The reduced risk of TB infection may be due to the fact that regions with high prevelance of HIV might have received attention in terms of HIV counseling, screening for TB and treatment. For instance, interventions such as early HIV counselling and screening for TB patients and early diagnosis and initiation of antiretroviral therapy (ART) to coinfected individuals have all been shown to be effective in preventing TB disease.52,53

We also studied the spatio-temporal distribution of TB infection risk using the spatio-temporal modeling frameworks proposed by Bernardinelli and colleagues,46 Waller and colleagues35 and Knorr-Held and Raßer.54 To determine the most suitable spatio-temporal model, we fitted these models and employed their individual Deviance Information Criterion (DIC) for selection. Among these models, the space-time model with the interaction term proposed by Bernardinelli and colleagues46 demonstrated the optimal fit to the TB data, as evidenced by its lowest AIC value. The findings indicated a decreased risk of TB infection across the country. Nonetheless, there was no significant reduction in TB infection risk over the study period, and there remained substantial clustering of risk among neighboring regions. High TB infection risk was identified consistently in the Southern part of the country throughout the years, although the infection risk consistently stayed below 1. Within the Southern part of Ghana, TB infection risk showed clustering among regions, with Brong-Ahafo, Volta, Ashanti, Eastern, and Greater Accra regions exhibiting elevated risk levels.

The data utilized in this study possesses a spatial dimension, encompassing reported TB cases across various districts, towns, and villages within the region. As a result, there are often shared factors55 that contribute to the prevalence of TB. This highlights the importance of identifying regions with high TB incidence, and facilitating the allocation of essential resources to address variables recognized as potential contributors to TB infection in those specific areas. The significance of regional-level modeling in the context of TB decision-making in Ghana is considerable. By employing regional-level models, it becomes possible to comprehend the localized dynamics and specific factors contributing to TB prevalence within distinct regions of the country. This approach facilitates more targeted and effective decision-making in resource allocation, intervention strategies, and the deployment of healthcare measures. It provides a deeper understanding of the unique challenges and variations in TB transmission, ensuring that interventions are tailored to address the specific needs of different regions in Ghana.

There is a notable absence of comprehensive literature addressing the modeling and mapping of tuberculosis (TB) risk in Ghana, particularly in terms of statistical methodology and practical application. The current body of literature has predominantly focused on spatial and seasonal modeling and mapping of TB risk,49,56 with no exploration of spatio-temporal risk across different years by various authors. While some studies have conducted modeling and mapping of TB risk without considering seasonal or spatio-temporal effects,14 others have incorporated the seasonal pattern of TB risk.49,56 Notably, there is a lack of studies investigating the yearly or monthly spatio-temporal effects on TB infection risk in Ghana. In general, the literature lacks an in-depth exploration of statistical methodologies and applications for modeling and mapping TB infection risk in Ghana, especially in the domain of spatio-temporal modeling using advanced statistical methods. This study aims to address this gap and contribute to the existing literature on the modeling and mapping of TB risk in Ghana. It serves as a pioneering effort in utilizing these statistical methodologies to examine spatial and spatio-temporal aspects of TB risk in the country. Emphasizing the use of such statistical methodologies in studying disease risk in Ghana, the study underscores their ability to provide precise estimates of disease risk in both space and time.

Conclusion

The identification of TB hot-spots in the Northern part of Ghana and the Greater Accra and Eastern regions in the Southern part of Ghana suggests a need for targeted interventions in these areas. Factors such as TB cure rate, TB success rate, knowledge about TB, income status, and HIV prevalence were identified as significant predictors of TB infection risk. This suggests that TB control measures should focus on treating more TB infected patients, creating TB awareness among the general public, making healthcare facilities available and accessible for both the poor and rich, and ensuring people living with HIV are identified early for treatment.

Similar levels of TB infection risk were observed among neighboring regions in Ghana, suggesting that these areas share common geographical traits influencing TB risk. This highlights the necessity for further research to identify such factors, enabling targeted interventions to enhance TB control measures.

Our study underscores the critical importance of precise risk assessment and mapping for effective resource allocation and targeted interventions. This research contributes valuable insights into disease hotspots across various regions in Ghana, utilizing advanced mapping techniques. The convolution model, incorporating both structured and unstructured components, facilitated spatial modeling with Bayesian methodology, revealing substantial variability and clustering of risk among neighboring regions. The findings unveiled distinct patterns of risk distribution, with certain regions exhibiting heightened risk over specific years.

In conclusion, our study enhances the understanding of TB dynamics in Ghana by employing sophisticated modeling techniques and innovative spatial and spatio-temporal analyses. These insights provide a foundation for informed decision-making, resource allocation, and targeted interventions, contributing to the ongoing efforts against TB and improving public health outcomes. Regions sharing similar risk patterns, possibly due to common local climatic conditions, can guide targeted resource allocation, given the regional basis of resource allocation in Ghana. The utilization of advanced statistical methods, encompassing both spatial and spatio-temporal analyses, is encouraged to further enrich the existing statistical literature on modeling and mapping TB risk in Ghana. We also recommend application of these methodologies to study spatial and spatio-temporal patterns of TB infection risk using TB data on all the 16 administrative regions in Ghana.

Data availability

Source data The data used in this study can be found in the following links: https://open.africa/dataset/4176f749-cfa8-4e32-9418-86cef78f9db6/resource/0bcf9b54-3e35-4543-95cd-fd4de953edff/download/factsfigures_2018.pdf, https://www.who.int/teams/global-tuberculosis-programme/data https://www.stoptb.org/static_pages/GHA_Dashboard.html.

Author contribution

Conceptualization, Software. Formal analysis: AKI Methodology and Investigation: AKI, FKB Data curation, Writing (Original draft preparation), and Writing (Review and editing): AKI and EAA Validation.

Comments on this article Comments (0)

Version 3
VERSION 3 PUBLISHED 16 Feb 2022
Comment
Author details Author details
Competing interests
Grant information
Copyright
Download
 
Export To
metrics
Views Downloads
F1000Research - -
PubMed Central
Data from PMC are received and updated monthly.
- -
Citations
CITE
how to cite this article
Iddrisu AK, A. Amikiya E and Kwame Bukari F. Spatio-temporal characteristics of Tuberculosis in Ghana [version 3; peer review: 1 approved, 1 approved with reservations] F1000Research 2024, 11:200 (https://doi.org/10.12688/f1000research.109053.3)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
track
receive updates on this article
Track an article to receive email alerts on any updates to this article.

Open Peer Review

Current Reviewer Status: ?
Key to Reviewer Statuses VIEW
ApprovedThe paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approvedFundamental flaws in the paper seriously undermine the findings and conclusions
Version 2
VERSION 2
PUBLISHED 06 Jan 2023
Revised
Views
17
Cite
Reviewer Report 08 Feb 2024
Dr. Victor Abiola Adepoju, Department of HIV and Infectious Diseases, Jhpiego (An Affiliate of Johns Hopkins University), Abuja, Nigeria 
Approved with Reservations
VIEWS 17
Abstract
Clarity and Focus: The abstract states, "The spatio-temporal characteristics of TB in Ghana using Bayesian spatial and spatio-temporal regression models was analysed." However, it doesn't detail the novel aspects of these methods or their specific application to the ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Adepoju DVA. Reviewer Report For: Spatio-temporal characteristics of Tuberculosis in Ghana [version 3; peer review: 1 approved, 1 approved with reservations]. F1000Research 2024, 11:200 (https://doi.org/10.5256/f1000research.142548.r237979)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
  • Author Response 13 Apr 2024
    Abdul-Karim Iddrisu, Mathematics and Staistics, University of Energy and Natural Resources, Sunyani, +233, Ghana
    13 Apr 2024
    Author Response
    We would like to thank the reviewer for the many useful and insightful comments. We have now revised the manuscript substantially. We trust that the revised manuscript now meets the ... Continue reading
COMMENTS ON THIS REPORT
  • Author Response 13 Apr 2024
    Abdul-Karim Iddrisu, Mathematics and Staistics, University of Energy and Natural Resources, Sunyani, +233, Ghana
    13 Apr 2024
    Author Response
    We would like to thank the reviewer for the many useful and insightful comments. We have now revised the manuscript substantially. We trust that the revised manuscript now meets the ... Continue reading
Version 1
VERSION 1
PUBLISHED 16 Feb 2022
Views
17
Cite
Reviewer Report 21 Dec 2022
Dioggban Jakperik, Department of Statistics, School of Mathematical Sciences, C. K. Tedam University of Technology and Applied Sciences, Navrongo, Ghana 
Approved
VIEWS 17
The study modelled TB incidence in Ghana using Bayesian Spatio-Temporal modelling. This covered the old ten (10) administrative regions of Ghana. Regions with high incidents, regarded as hotspots have been identified. The relevant TB and patient characteristics have been considered. ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Jakperik D. Reviewer Report For: Spatio-temporal characteristics of Tuberculosis in Ghana [version 3; peer review: 1 approved, 1 approved with reservations]. F1000Research 2024, 11:200 (https://doi.org/10.5256/f1000research.120511.r156972)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
  • Author Response 09 Jan 2023
    Abdul-Karim Iddrisu, Mathematics and Staistics, University of Energy and Natural Resources, Sunyani, +233, Ghana
    09 Jan 2023
    Author Response
    Comment 1: "Grammatical errors need to be addressed."
    Response 1We thank the reviewer for the comment. In response, we have proof read the manuscript for spellings and grammatical errors
    ... Continue reading
COMMENTS ON THIS REPORT
  • Author Response 09 Jan 2023
    Abdul-Karim Iddrisu, Mathematics and Staistics, University of Energy and Natural Resources, Sunyani, +233, Ghana
    09 Jan 2023
    Author Response
    Comment 1: "Grammatical errors need to be addressed."
    Response 1We thank the reviewer for the comment. In response, we have proof read the manuscript for spellings and grammatical errors
    ... Continue reading

Comments on this article Comments (0)

Version 3
VERSION 3 PUBLISHED 16 Feb 2022
Comment
Alongside their report, reviewers assign a status to the article:
Approved - the paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations - A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approved - fundamental flaws in the paper seriously undermine the findings and conclusions
Sign In
If you've forgotten your password, please enter your email address below and we'll send you instructions on how to reset your password.

The email address should be the one you originally registered with F1000.

Email address not valid, please try again

You registered with F1000 via Google, so we cannot reset your password.

To sign in, please click here.

If you still need help with your Google account password, please click here.

You registered with F1000 via Facebook, so we cannot reset your password.

To sign in, please click here.

If you still need help with your Facebook account password, please click here.

Code not correct, please try again
Email us for further assistance.
Server error, please try again.