Exploring the spatial and temporal dynamics of the health index of Bangladesh using Bayesian framework

Abstract

Background:

Health is vital for survival and societal well-being, with global goals focusing on improved access and achieving progress in public health. In Bangladesh, investigating health determinants and disparities is crucial for effective interventions.

Objective:

This study investigates spatial and temporal variations in the influence of socioeconomic and health-related factors on the health index at a detailed district level in Bangladesh, aiming to discern disparities in health indicators among districts. Based on 21 years of secondary data across 64 districts in Bangladesh, this research deploys four Bayesian spatial-temporal models to analyze the health index and relevant variables.

Design and methods:

The study identifies the optimal model by evaluating linear, analysis of variance, and autoregressive models of orders one and two. Comparative analysis utilizing the WAIC and DIC aids model selection. The MCMC method enhances precision in determining health index determinants, providing valuable insights for policymakers and researchers.

Results:

Our findings reveal that the AR (2) model demonstrated superior performance among all the spatiotemporal models. The health index demonstrated a consistent upward trend overall, but regional variations in residents’ health indices were noted. Positive spatial autocorrelation was identified in the health index, emphasizing interconnected health outcomes. Spatiotemporal models underscored the significance of the income index, child mortality, life expectancy, and the proportion of the active population as crucial factors shaping the health index in Bangladesh.

Conclusions:

This study provides actionable insights for targeted health interventions in districts with specific needs, emphasizing focused policies on poverty, education, and healthcare to enhance overall well-being in Bangladesh.

Keywords

health index spatial-temporal approach Bayesian inference MCMC life expectancy child mortality Bangladesh

Introduction

The role of human capital in the development process involves prioritizing the preservation, growth, and enhancement of the well-being of human populations and is now widely acknowledged as a pivotal strategy for fostering sustainable development, particularly in developing nations like Bangladesh.¹ Access to representative and comparable health data within a nation is essential for formulating evidence-based health policies, planning educational and awareness initiatives, and conducting public health research. Governments globally allocate substantial budgets to healthcare expenditures, recognizing its indispensability for enhancing human productivity and sustaining a healthy nation. Since gaining independence in 1971, Bangladesh has demonstrated a steadfast commitment to progress, achieving notable success in economic growth, poverty reduction, healthcare, and education. This remarkable development trajectory, often termed the “Bangladesh Miracle,” has transformed the nation from a struggling, war-torn country into one of the world’s most promising emerging economies.² However, the government must address the concept of balanced regional growth, a constitutional requirement.³ Disparities in regional development are common in nations like Bangladesh, reflecting uneven growth. Previous plans and initiatives for the country’s development lacked a regional approach, resulting in variations in socioeconomic development across regions. These differences manifest in literacy rates, health indices, income levels, child health, life expectancy, and other critical indicators. In navigating the challenges of regional disparities, the government must adopt strategic and targeted interventions to ensure balanced and inclusive development across the nation. Smart and region-specific policies and initiatives are essential to address the existing imbalances and pave the way for holistic and sustainable progress.

In a session of the United Nations General Assembly, nations pledged to prioritize health investments in the resolution titled “Transforming our world: the 2030 Agenda for Sustainable Development.” This commitment involved a dedication to achieving universal health coverage and reducing health inequalities for people of all ages.^4,5 Bangladesh actively implemented the “Health, Population, and Nutrition Sector Program (HPNSP)” in 2022, showcasing the nation’s commitment to advancing healthcare, population management, and nutrition initiatives.⁵ Various macro-level studies have identified a range of factors linked to the overall health status, including the level of socioeconomic development, educational attainment, environmental factors, healthcare expenditures, urbanization, and lifestyle choices.^6,7 Various factors have been utilized in health studies to elucidate individuals’ health status. Numerous studies, including references,^6–10 underscore the significant role of income as a determinant of health. Education is closely linked to better healthcare outcomes. An educated population is more likely to have access to healthcare services, leading to improved public health. Education also emerges as a crucial factor influencing health outcomes in various ways.^1,6,9,11,12 Individuals with higher levels of education tend to secure better employment opportunities and higher incomes and engage in fewer risky behaviors. Investments in education enhance productivity and income levels, according to research, which helps the economy flourish.¹³ Urbanization is identified as another determinant of health^7,9 presenting both positive and negative impacts. Positive aspects, such as increased access to medical facilities and information, and negative impacts, like pollution, which can affect overall health.¹⁰ In the contemporary landscape, life expectancy stands out as a key determinant for gaging population health. According to a study,⁹ there are regional implications for life expectancy because of uneven health facilities. Life expectancy serves as a statistical measure representing the average time an individual or creature is expected to live from birth to their current age. Again, according to WHO, a more active population is required for a healthier world, which is also an important determinant of the health index.¹⁴

This study pioneers the application of spatial patterns, autocorrelation, and spatio-temporal models in assessing the distribution of health indices in Bangladesh from 2001 to 2021. Employing a Bayesian Hierarchical Spatial-temporal model, the research identifies the most suitable spatial model, offering valuable insights for decision-makers.¹⁵ This innovative approach, a pioneering effort in Bangladesh, integrates spatial and temporal dimensions, providing a comprehensive understanding of health facility distribution. The study’s findings aim to inform targeted interventions and serve as a foundational resource for future spatio-temporal research endeavors.

Methods and materials

Data source

The Global Data Lab’s Area Database (https://globaldatalab.org/areadata/) was considered the main source of data that is accessed and free to use.¹⁶ The GDL is an independent data and research center at the Nijmegen School of Management (https://www.ru.nl/nsm/) of Radboud University (https://www.ru.nl/). Bangladesh Bureau of Statistics (BBS) is also used as a secondary source for our data.¹⁷ The data was collected for the time duration of 2001–2021. Sixty four districts of Bangladesh (the second level of administration) are considered as the spatial domain in this study; therefore, all the data is collected for each district and each predefined period (yearly). In cases where district-level data were missing, information from neighboring districts was used to impute the missing values and ensure a complete subnational dataset.

Study variables

The health index refers to a composite measure that assesses various health-related factors, including income, healthcare accessibility, quality of life, and health outcomes in different regions in Bangladesh, which is the main response variable of this study. In this research, we have considered the health index according to the district for the predefined period as our response variable. To support our research objectives, we included several covariates drawn from existing literature and their relevance to health outcomes in Bangladesh.

These include the literacy rate, which influences health awareness, preventive behavior, and the ability to access and interpret health information^11,14; the income index, which reflects economic capacity to afford healthcare services and maintain a healthy lifestyle^8,10; and life expectancy and child mortality, which are both widely recognized as direct measures of population health status.^6,9,14,18 The proportion of the active population (working-age individuals) is included as it affects productivity, dependency ratios, and resource distribution, while population density captures the effects of congestion, resource competition, and disease transmission.^7,19,20 Urbanization reflects access to infrastructure, sanitation, and health services, but also potential exposure to pollution and urban stressors.^21,22 These covariates provide a broader context for understanding the factors influencing the health index.^6-13

Additionally, the study’s spatial domain encompassed 64 distinct district names, while the temporal domain spanned 21 consecutive years across these districts. Consequently, each variable within the dataset comprised a total of 1344 observations (64 districts multiplied by 21 years), reflecting the comprehensive nature of the study across both spatial and temporal dimensions.

Spatio-temporal modeling

Examining, evaluating, and modeling spatiotemporal data characteristics that indicate the geometric or geographic properties of the data along with temporal characteristics is referred to as spatio-temporal analysis. The focus of spatio-temporal analysis is on variation in the variable in combinations of place and time units over the study area and study period.

Spatio-temporal pattern analysis

An integral aspect of spatio-temporal analysis involves examining the patterns of a response variable over space and time. Spatial autocorrelation is employed as a method to analyze these patterns. Here, autocorrelation in spatio-temporal data has been analyzed using Moran’s I (Global and Local Moran’s I) and Geary’s C²³ to observe the spatial pattern of the response variable (health index).

Global Moran’s I

Global Moran’s I is a widely used measure of overall spatial autocorrelation, defined as

l_{i} = \frac{n}{\sum_{i = 1}^{n} {(y_{i} - \bar{y})}^{2}} \times \frac{\sum_{i = 1}^{n} \sum_{j = 1}^{n} w_{i j} (y_{i} - \bar{y}) (y_{j} - \bar{y})}{\sum_{i = 1}^{n} \sum_{j = 1}^{n} w_{i j}}

Where $y_{i}$ is the health index, $\bar{y}$ is the arithmetic mean across the entire study area (district), $w_{i j}$ is the spatial weight between i and j, and n is the total number of features. A significantly positive Moran’s I indicates positive spatial autocorrelation, while a negative value implies spatial dispersion.

Local Moran’s I (LISA)

To detect spatial heterogeneity, Local Indicators of Spatial Association (LISA) were computed. The local Moran’s I for district $i$ is given by

l_{i} = \frac{n}{\sum_{i = 1}^{n} {(y_{i} - \bar{y})}^{2}} \times \frac{\sum_{i = 1}^{n} \sum_{j = 1}^{n} w_{i j} (y_{i} - \bar{y}) (y_{j} - \bar{y})}{\sum_{i = 1}^{n} \sum_{j = 1}^{n} w_{i j}}

This statistic identifies localized clusters, including high-high (HH) and low-low (LL) groupings, as well as spatial outliers (high-low or low-high), which are essential for targeted policy decisions.

Geary’s C

In addition to Moran’s I, Geary’s C was used to capture local dissimilarities and spatial contrast:

C = \frac{(n - 1) \sum_{i = 1}^{n} \sum_{j = 1}^{n} w_{i j} {(y_{i} - y_{j})}^{2}}{2 W \sum_{i = 1}^{n} {(y_{i} - \bar{y})}^{2}} W

is the sum of all spatial weights. While both Moran’s I and Geary’s C assess spatial autocorrelation, Geary’s C emphasizes local dissimilarity. A value of $C < 1$ suggests positive spatial autocorrelation, while $C > 1$ implies negative spatial autocorrelation.

Although in theory Geary’s coefficient and Moran’s index are comparable, Geary’s coefficient is based on a sample, and Moran’s index is based on the population.^21,22 In this study, spatio-temporal autocorrelation checking has been done as an extension of Moran’s I and Geary’s C.²⁴

Bayesian hierarchical model

In this context, the response variable, denoted as the health index $y_{i t}$ where i = 1. . ., 64 and t = 1. . ., 21 was hypothesized to adhere to a Gaussian distribution, $y_{i t} ~ N (μ_{i t}, {σ^{2}}_{i t})$ . And the probability density function denoted as

f (y_{i t}, μ_{i t}, {σ^{2}}_{i t}) = \frac{1}{\sqrt{2 π {σ^{2}}_{i t}}} e^{(- \frac{1}{2 {σ^{2}}_{i t}} {(y_{i t} - μ_{i t})}^{2})}, μ_{i t} ϵ ℝ, y_{i t} ϵ ℝ

Where, $y_{i t}$ is health index, $µ_{i t}$ is the mean and ${σ^{2}}_{i t}$ is the variance of defined response variable at space-time combination. Let’s say that $x_{i t}$ represents the p-dimensional spatio-temporal covariates at the space-time combination. Then the model can be written as

y ~ x_{i t}^{'} β + ψ_{i t} + ϵ_{i t} i = 1, \dots, n, t = 1, \dots, T

Where $ϵ_{i t} ~ N (0, ϑ^{2})$ (distributed independently), $ψ_{i t}$ are spatio temporal random effects, and $ϑ^{2}$ is the error variance.

For independent error general linear regression model matrix of unknown regression coefficients is denoted as $β$ and the random effects $ψ_{i t}$ are assumed to have the following diverse range of models based on various assumptions regarding their architecture. To extend the process and parameter models to both space and time dimensions, several models have been used in research.²⁵ Knorr-Held also proposed an analysis of the variance type model with or without the inclusion of the space-time interaction.²⁶ The general spatial-temporal models were considered as

ψ_{i t} = {\begin{matrix} β_{1} + ϕ_{i} + (β_{2} + δ_{i}) \frac{t - \bar{t}}{T}, L i n e a r m o d e l o f t r e n d \\ ϕ_{i} + δ_{i} + γ_{i t}, A N O V A m o d e l \\ ϕ_{i t} + δ_{t}, A u t o R e g r e s s i v e m o d e l \end{matrix}

Here, $β_{1}$ and $β_{2}$ are overall intercept and slope (trend) parameters which were given a flat prior distribution in the Bayesian model. The parameters $ϕ_{i}$ and $δ_{i}$ are incremental intercept and slope parameters for the $i^{t h}$ region, $ϕ = ϕ_{1 \dots,} ϕ_{64}$ and $δ = δ_{1 \dots,} δ_{64}$ are assigned the NCAR ( $ϕ | ρ_{i n t}, {τ^{2}}_{i n t}, W)$ and NCAR ( $δ | ρ_{s l o}, {τ^{2}}_{s l o}, W)$ distribution.²⁶ Here $ρ_{i n t}, ρ_{s l o}, {τ^{2}}_{s l o}, {τ^{2}}_{s l o}$ are autoregression and variance parameters for the intercept ( $ϕ_{i}$ ) and slope ( $δ_{i}$ ) processes. In the ANOVA model, the three sets of parameters $ϕ_{i}$ , $δ_{i}$ , and $γ_{i t}$ are all random effects, also following NCAR prior distributions.²⁶ In the model using NCAR priors, the W represents the spatial neighborhood matrix. This matrix is a critical component of Conditional Autoregressive (CAR) and its extension, the Non-Centered Autoregressive (NCAR) priors. In this study, the utilization of the “bmstdr” package and the Bcartime function, among other tools within the R programing language, is integral to performing geographical analysis and model fitting.²⁷ Additional spatial and temporal domains are incorporated, employing MCMC methods to estimate priors for posterior estimation. The “CARBayes” package defaults to an inverse gamma prior (shape = 1, scale = 0.01) for the unknown variance parameter in the Gaussian model. The “ggplot2” package in RStudio was used to display the indices, and shapefiles are read using the “readOGR” function. To select the appropriate model, Bayesian model selection criteria such as the Log Marginal Predictive Likelihood (LMPL), Deviance Information Criterion (DIC)²⁸ and Watanabe-Akaike information criterion (WAIC)²⁹ were used in this study. RMSE and MAE were used to evaluate the accuracy of the best model.

Results

Exploratory data analysis

Exploratory Data Analysis (EDA) is a crucial phase in the data analysis process where the primary goal is to gain insights and a deeper understanding of the dataset. It involves the use of statistical graphics, plots, and summary statistics to explore the main characteristics, patterns, and relationships visually and quantitatively within the data.

Distribution of response variable

It is assumed that $Y_{i t}$ denotes the annual health index of people of $i^{t h}$ district and $t^{t h}$ year. Here i denotes for district and i = 1, 2, . . ., 64 districts in Bangladesh and t = 1, 2, . . ., 21 for 21 years from 2001 to 2021. The following figure of spatial distribution and time series representation shows the variance in the health index in Bangladesh over time in 64 districts.

From the above graph, it is seen that over the past 21 years, the health index has been rising in Bangladesh (Figure 1(a)). In the most recent year of 2019, Bangladesh had a maximum health index of 81.41. Figure 1(a) demonstrates that the annual health index steadily rose with a decline in the year 2009, reaching its peak in 2019, but after that year it started to decline due to the global COVID-19 pandemic scenario. After the coronavirus storm, the annual health index is once again rising. Again Figure 1(b), which is a spatial representation of the health index represents that, compared to other areas, the capital (Dhaka), which includes Bangladesh’s southeast and northeast, had relatively good health improvement. Additionally, compared to other areas of Bangladesh, the coastline region has fewer health facilities.

Figure 1.

Annual Health Index: (a) time series plot and (b) spatial distributions.

Correlation matrix of variables

Utilizing pairwise correlations, Figure 2 illustrates Bangladesh’s health index connections with relevant covariables. Notably, strong positive correlations exist between the health index and the life expectancy of people. Income, education, and active population also strongly align with the health index (0.67, 0.77, 0.73), supporting the growth of high-quality healthcare facilities in financially stable, educated, and active communities.

Figure 2.

Correlation matrix of the selected variables.

Urbanization and population density have a positive impact on health, specifically in urban areas with higher indices, as indicated by correlation coefficients of 0.14 and 0.17, respectively. However, the overall correlation is less, influenced by regional disparities. The inverse correlation between urbanization and child mortality underscores effective medical services. Notably, no discernible link between the health index and Dhaka from a distance indicates minimal geographical impact. This study offers compelling evidence that financial standing, educational attainment, and an active population foster health consciousness and ensure widespread access to healthcare in Bangladesh.

Spatial analysis of influencing factors

Figure 3 underscores pronounced regional differences in Bangladesh, highlighting Dhaka’s high population density compared to sparser mountainous areas and distant northeastern districts.

Figure 3.

District-wise variation of significant covariates.

Living conditions vary widely, with substandard levels prevalent in many areas, while higher income indices are concentrated in Dhaka and the northwest. The literacy rate distribution mirrors these trends, with Dhaka and the eastern regions leading in terms of active population, while the western districts lag.

Child mortality, as depicted in Figure 3, has varied significantly over the past 21 years, with Chuadanga, Kushtia, and Jhenaidah reporting lower rates than Feni, Noakhali, and Lakshmipur, correlating negatively with health indices. Dhaka City and the southeast exhibit higher life expectancy, positively linked to health indices, whereas Sunamganj and Habiganj show lower life expectancy. Urbanization is concentrated in Dhaka, Chittagong, and select districts, emphasizing Bangladesh’s uneven development. Figure 3 visually portrays Dhaka as the most populous district, with the mountainous region having a lower population density. The northeastern districts remain distant from the capital, encapsulating Bangladesh’s geographical diversity.

Spatio-temporal modeling

Global Moran’s I and Geary’s C for spatial pattern analysis

To examine spatial patterns using autocorrelation with Global Moran’s I let set the hypothesis as,

H₀ = $\begin{array}{l} Spatial interdependence is absent among \\ neighboring districts concerning the health index \end{array}$

H₁ = $\begin{array}{l} Spatial interdependence is present among \\ neighboring districts concerning the health index \end{array}$

The hypothesis mentioned above is considered while analyzing autocorrelation with Global Moran’s I. The statistics for the global Moran’s I and Geary’s C of the health index are shown in Table 1 below.

Table 1.

Moran’s I and Geary’s C statistics under randomization for health index.

Statistics	Global Moran’s I	Geary’s C
Indicator statistic	0.724	0.542
Expectation	−0.0158	1.000
Variance	0.0068	0.0083
Standard deviate	8.982	5.0174
p-Value	$2.2 \times 10^{- 10}$	$2.619 \times 10^{- 7}$

A good spatial autocorrelation between regional domains was found using Moran I’s, which produced a score of 0.724. The obtained p-value of $2.2 \times 1 0^{- 10}$ , which is less than 0.05, suggests that the null hypothesis should be rejected using MC simulation of the 19999 Global Moran, that is there is spatial dependence between neighboring districts of Bangladesh due to health index. In both cases, the null hypothesis should not be taken into consideration.

Again, Figure 4(a)’s density plot of Monte Carlo permutation results illustrates the observed Moran’s I value at 0.724, indicating positive spatial autocorrelation. The Geary’s C statistic, computed at 0.542 (within [0,1)), further confirms positive autocorrelation between districts' health indices.

Figure 4.

(a) Density plot of Global Moran’s I and (b) graphical representation of spatial – autocorrelation.

Both analyses unveil clusters in health indices. The “Moran scatter plot” (Figure 4(b)) reinforces spatial autocorrelation, depicting data points close to one another. The significantly positive Moran’s I suggest spatial clustering, prompting further exploration for potential hotspots. This insight guides subsequent model refinement in the MCMC setting.

Local Moran’s I

The spatial weights of items form the foundation for the Local Moran’s I statistic. Figure 5’s cluster map illustrates districts with notable geographic homogeneity, deemed significant at a 90% confidence level. This determination takes into consideration the effects of both spatial lag and the spatial weights assigned to adjacent districts. North-east and southeast regions of Bangladesh without Brahamanbaria, Comilla, and Chandpur districts have a high Low-Low health facilities cluster. Again, the south-west part of Bangladesh including four northeast districts Panchagarh, Thakurgoan, Dinajpur, and Nilphamari shows a dark red color, indicating that districts have a high health index with a high health index in neighborhood districts. At the same time, the brown color indicates that four districts are in the High-Low cluster and the dark blue color represents five districts are in the Low-Heigh cluster.

Figure 5.

Spatial clustering (local Moran’s I) of health index.

Bayesian spatiotemporal modeling

Exploratory spatial data analysis (ESDA) faces challenges in capturing the temporal dynamics of geographical features, despite its ability to depict complex spatial patterns.³⁰ For a comprehensive understanding of both spatial and temporal effects, Bayesian spatio-temporal modeling proves effective. It is advisable to initiate modeling with an Independent Error Bayesian model before incorporating spatial and temporal dimensions. This step allows for a clear appreciation of the enhancements introduced when integrating spatial and temporal effects into the model. The linear predictor $x_{i t}^{’} β$ can be specified in the model as,

Educational Index + Income Index + Child Mortality + Active Population + Life Expectancy + \log (Population Density) + Urbanization + Distance

Population density is in log scale and distance from Dhaka divided by 10 for modeling the health index. After 20,000 burn-in iterations, all the models were run for 120,000 iterations. To reduce autocorrelation in the MCMC setting, the data were kept following a thinning of 10 iterations. The estimated results from the fitted bmstdr model, encompassing confidence intervals, are presented in Table 2 for all five models.

Table 2.

Estimated parameters with their credible interval for all proposed models.

	Independent error Bayesian GLM	Bayesian spatio-temporal model
	Independent error Bayesian GLM	Linear model	ANOVA model	AR (1) model	AR (2) model
	Estimate	Estimate	Estimate	Estimate	Estimate
Model Parameters	Credible interval	Credible interval	Credible interval	Credible interval	Credible interval
Intercept	−15.953 (−19.628, −12.249)	−16.387 (−20.027, −12.733)	−15.374 (−19.182, −11.523)	−15.438 (−19.220, −11.715)	−15.442 (−19.200, −11.716)
Educational index	−0.052 (−0.076, −0.027)	−0.056 (−0.080, −0.030)	−0.073 (−0.100, −0.047)	−0.063 (−0.089, −0.037)	−0.063 (−0.090, −0.036)
Income index	0.123 (−0.011, 0.057)	0.126 (0.008, 0.162)	0.253 (0.016, 0.391)	0.138 (0.002, 0.075)	0.138 (0.035, 0.175)
Child mortality	−0.018 (−0.030, −0.007)	−0.018 (−0.030, −0.007)	−0.022 (−0.034, −0.010)	−0.021 (−0.033, −0.010)	−0.021 (−0.033, −0.009)
Active population	0.100 (0.066, 0.135)	0.108 (0.072, 0.143)	0.125 (0.088, 0.162)	0.112 (0.075, 0.149)	0.112 (0.074, 0.149)
Life expectancy	1.265 (1.218, 1.313)	1.267 (1.220, 1.313)	1.228 (1.179, 1.278)	1.245 (1.196, 1.295)	1.246 (1.196, 1.295)
Log (Pop_ Density)	−0.001 (−0.130, 0.128)	−0.01 (−0.139, 0.117)	0.027 (−0.124, 0.178)	0.016 (−0.126, 0.153)	0.014 (−0.129, 0.158)
Urbanization	0.005 (−0.002, 0.013)	0.003 (−0.004, 0.011)	−0.005 (−0.013, 0.004)	−0.001 (−0.009, 0.008)	−0.001 (−0.010, 0.008)
Distance	0.004 (−0.004, 0.013)	0.005 (−0.003, 0.014)	0.003 (−0.006, 0.012)	0.003 (−0.006, 0.012)	0.003 −15.442
${τ^{2}}_{S}$	-	-	0.131 (0.074, 0.224)	-	-
${τ^{2}}_{T}$	-	-	0.196 (0.102, 0.364)	0.311 (0.130, 0.710)	0.304 (0.129, 0.769)
$ν^{2}$	1.625 (1.507, 1.753)	1.591 (1.473, 1.715)	1.551 (1.433, 1.676)	1.44 (1.227, 1.593)	1.435 (1.185, 1.592)
$ρ_{S}$	-	-	0.578 (0.163, 0.938)	0.838 (0.305, 0.975)	0.828 (0.253, 0.975)
$ρ_{T}$	-	-	0.370 (0.025, 0.821)	0.145 (0.006, 0.414)	-
$α$	-	−0.305 (−0.550, −0.056)	-	-	-
${τ^{2}}_{i n t}$	-	0.135 (0.075, 0.230)	-	-	-
${τ^{2}}_{s l o}$	-	0.382 (0.142, 0.864)	-	-	-
$ρ_{i n t}$	-	0.566 (0.154, 0.930)	-	-	-
$ρ_{s l o}$	-	0.477 (0.051, 0.921)	-	-	-
$ρ 1_{T}$	-	-	-	-	0.022 (−0.385, 0.464)
$ρ 2_{T}$	-	-	-	-	0.023 (−0.466, 0.458)

Table 2 incorporates various spatiotemporal indicators for the five described models. The findings consistently highlight the substantial influence of factors on the health index across all models. The parameter estimates in Table 2 show that the proportion of the active population, the life expectancy and the income are significant predictors while considering health development. Again, child mortality is not significant in the model, which is like pairwise correlation. In Independent Error Bayesian GLM, there is a high impact of life expectancy on health development, which can be seen from the above explanations. Autoregression variance $ν^{2}$ is 1.625, as the spatial and temporal variance is not considered yet. Here, for the linear model $ρ_{i n t}$ = 0.566 indicates autoregression for the intercept is greater than the slop $ρ_{s l o}$ = 0.477 and ${τ^{2}}_{i n t}$ = 0.135 means the variance for the intercept is lower than the variance for the slop, ${τ^{2}}_{s l o}$ = 0.382. The ANOVA model in Table 2 displays a spatial correlation value of $ρ_{S}$ is 0.578, indicating significant spatial influence, and it is higher than the value of $ρ_{T} = 0.370$ which indicates that the health index is strongly dependent on spatial variation along with temporal effect. The values of ${τ^{2}}_{S}$ and ${τ^{2}}_{T}$ for the ANOVA model in Table 2 are non-zero, which supports the necessity of taking the spatial and temporal domain into account in the model as well as the existence of geographical and temporal variation in the health index. Table 2 lists the parameter estimates for the selected AR(1) model, and according to estimates, the temporal autocorrelation parameter $ρ_{T}$ is roughly almost six times smaller than the spatial correlation value $ρ_{S .}$ Therefore, it may be inferred from the spatial and temporal autocorrelation values that there is a greater impact on space than time. The following table shows that the income index, percentage of active population income, and life expectancy of people significantly affect the health index. And the child mortality index is still insignificant for this model. The above Table 2 also presents the spatial correlation value $ρ_{S}$ = 0.994, which indicates the existence of high spatial correlation. Both the AR parameters show a negative lagged correlation, which can be interpreted from the confidence interval of the estimated value.

Selecting a well-fitted model is of utmost importance to accurately capture the data patterns. Table 3 and Figure 6 serve as representations of the model selection and validation criteria, facilitating the identification of the most suitable model to capture and express the underlying patterns within the data effectively.

Table 3.

Assessment criteria for the chosen models’ performance.

Models	Model choice criteria						Model validation criteria
Models	DIC	p. d	WAIC	p. w	LMPL	loglikelihood	RMSE	MAE	CRPS	CVG
IEBGLM	4480.339	9.999	4488.121	17.395	−2244.06	−2230.17	1.102	0.734	0.804	97.76
Linear	4487.867	47.451	4496.772	54.082	−2248.53	−2196.48	1.245	0.802	0.886	97.76
AR (1)	4440.753	136.84	4397.204	152.812	−2233.91	−2083.54	1.251	0.808	0.709	97.015
AR (2)	4432.09	139.72	4336.386	174.175	−2235.16	−2076.33	1.247	0.704	0.616	97.015
ANOVA	4457.907	52.65	4466.956	59.182	−2233.61	−2176.3	1.232	0.813	0.785	97.015

Values in bold indicate statistical significance at the 5% level.

Figure 6.

Comparison plot for performance measures of the selected models.

Table 3 consolidates the performance measures of the chosen models, utilizing both model selection and validation criteria. Additionally, Figure 6 illustrates model comparison plots based on the MCMC sample, further aiding in the decision-making process to identify the most suitable model. Table 3 and Figure 6 provide insights into model performance and validation criteria, highlighting that the AR(2) model surpasses others by yielding the lowest WAIC and DIC values. Moreover, AR(2) demonstrates superior log marginal predictive likelihood (LMPL) and log-likelihood values compared to alternative spatiotemporal models in the context of our health-related data. The model’s excellence is further underscored by lower MAE and CRPS values in model validation criteria. Notably, the spatio-temporal linear model exhibits the highest convergence rate. Therefore, it is established that the autoregressive model of order two is the best model for the dataset of health index by the model performance measures.

Figure 7 represents the coverage percentage of MCMC estimation. To assess the best-fitted model, 10,000 samples from the MCMC are employed to compute residuals and their standard deviation for spatial model fitting. Figure 8 presents a residual plot for each district and a plot depicting the standard error of the residuals. The standard error plot indicates minimal variation in residual values. Notably, the plots in Figure 8 do not unveil any significant spatial patterns warranting further investigation. Instead, they depict spatially aggregated residuals and the standard deviation of residuals derived from the spatio-temporal AR(2) model.

Figure 7.

Coverage percentage using MCMC Bayesian estimation.

Figure 8.

(a) Aggregated spatial residuals and (b) standard error of the residuals.

The temporally fitted plot spanning the entire 21-year period indicates a balanced distribution of fitted values. There is neither a significant overestimation nor underestimation, suggesting the AR (2) model’s accuracy within the spatial-temporal framework. The proximity of the lower and upper lines to the fitted values, as depicted in Figure 9, further reinforces the model’s reliability in capturing the temporal dynamics of the data.

Figure 9.

Time-series representation of observed and fitted health index values along with 95% confidence intervals.

Discussion

This study investigates the spatiotemporal effects of various factors on the health index in Bangladesh. Regional disparities are a common phenomenon, seen in both developed and developing nations, with more pronounced variations in the latter, exemplified by Bangladesh. In our study,^31,32 we employed spatial autocorrelation to find clusters and Bayesian spatial-temporal models to assess influential factors of health index in our country. The study aims to uncover spatial variations in health facilities at regional levels in the country. The results reveal that, at the district level in Bangladesh, the health index has shown growth over the years (2001–2021). However, notable inequities persist, indicating varying rates among districts. Consistent with findings in other studies^6–8,33 income per capita demonstrates a positive impact on health. Higher income levels contribute to increased access to high-quality goods and services, improved housing, and enhanced medical care services, collectively influencing overall health status.^7,34 Contrary to the conventional positive or negative associations found in many studies, our analysis reveals a positive correlation between income and the health index in our country, acknowledging regional disparities. Regarding spatio-temporal effects, the findings on education differ from those of previous studies. Typically, higher levels of education are associated with increased sensitivity and awareness about health, prompting individuals to take proactive measures to enhance both the quantity and quality of their health.⁷ In our analysis, we observed no significant effect of education on the health index, accounting for both spatial and temporal variations. This finding contrasts with the expectations based on general trends, emphasizing the nuanced and context-specific nature of the relationship between education and health in our study. This may be due to multicollinearity among predictors such as income and life expectancy, which are often strongly correlated with education. In Bayesian models, multicollinearity can inflate posterior variances and obscure individual effects, even with weakly informative priors.^35,36 While our model accounted for spatial and temporal dependencies, further sensitivity analysis, such as testing alternative prior settings or excluding correlated covariates, could help clarify the true role of education. A study conducted by Das indicated the potential to improve the health index in areas with low levels of education through systematic awareness-raising programs. The implementation of such programs could lead to increased service coverage in those regions.³⁷ However, it’s crucial to note that certain variables, such as urbanization, population density, and distance, exhibited insignificance in their impact on health development in Bangladesh. This outcome can be attributed to the economic landscape of Bangladesh, which is characteristic of a developing country. In particular, the uneven distribution of urbanization across different districts within Bangladesh likely contributed to the insignificance of the urbanization variable in the context of health development. It is noteworthy that, in some studies, the argument is made that individuals in urban areas possess greater access to medical care services and health information, thereby suggesting a positive impact of urbanization on the health index.^10,38 In our study, we also found a positive impact of life expectancy and the proportion of the active population on the health index. Simultaneously, the observed decreasing trend in child mortality is indicative of a positive effect on the health index. Consistent with prior research, our study aligns with the recommendation to mitigate child mortality by ensuring the effective utilization of health services.^39,40 Our study underscores regional disparities in health indicators, emphasizing the need for targeted community-level interventions aimed at poverty reduction, universal education, and addressing geographical disparities in healthcare access.^32,40 Specifically, our findings suggest that efforts should be differentially focused on enhancing child health care utilization, enough health facilities, and improving socioeconomic status, especially in disadvantaged regions.

Policy implications, practical applications, and limitations

To enhance the applicability of our findings in real-world policymaking, especially in the context of Bangladesh, it is essential to translate statistical insights into concrete actions. To make these findings more actionable, practical interventions such as deploying community health clinics, community health workers, and telemedicine initiatives could be prioritized in low-performing districts like Kurigram, Bandarban, and Sunamganj, where access remains limited.⁴¹ Insights from successful regional models, including Nepal’s Female Community Health Volunteers (FCHVs) and India’s Accredited Social Health Activists (ASHAs), show how community-based health outreach can improve maternal and child health.^42,43 Such models can be adapted in Bangladesh through collaboration with organizations like BRAC or Gonoshasthaya Kendra. Better policy targeting can also be ensured by district-level health monitoring and resource allocation that is directed by the spatial clusters found in this study. This study offers a basis for well-informed, fair, and context-sensitive health policymaking in Bangladesh by coordinating evidence with workable strategy.^13,44 Despite the strengths of this study, some limitations should be acknowledged. The reliance on secondary data over 21 years may involve inconsistencies in data collection methods, missing values, and potential measurement errors. Although imputation from neighboring districts was applied where necessary, these factors may still influence the robustness of the results.

Conclusion

The health index serves as a comprehensive measure of overall health, combining various indicators. This study investigates the spatiotemporal effects on the health index in 64 districts in Bangladesh, utilizing Bayesian spatial-temporal models. The AR (2) model emerges as the optimal choice, capturing spatial and temporal dynamics effectively. Again, the spatial pattern of the health index shows four clusters of the 64 districts, which shows how similar different districts are within the clusters. Findings indicate both spatial and temporal influences on the health index, highlighting significant impacts of life expectancy, income, proportion of active population, and child mortality on the health index. The findings of this study offer policymakers a concise overview of various districts in Bangladesh, encompassing crucial health indicators. This information empowers them to formulate essential policies and implement targeted development strategies based on a comprehensive understanding of the health landscape across different regions. The study underscores the need for targeted community-level programs, particularly in underprivileged districts, emphasizing a reevaluation of government strategies. Policymakers should prioritize initiatives based on research outcomes. To achieve universal health coverage and overcome barriers to healthcare access, the study recommends decentralizing healthcare services and initiating health promotion initiatives in regions with fewer facilities. The study provides a holistic understanding of regional development disparities, offering valuable insights for future research on health indices in Bangladesh.

Supplemental Material

sj-docx-1-phj-10.1177_22799036251351156 – Supplemental material for Exploring the spatial and temporal dynamics of the health index of Bangladesh using Bayesian framework

Supplemental material, sj-docx-1-phj-10.1177_22799036251351156 for Exploring the spatial and temporal dynamics of the health index of Bangladesh using Bayesian framework by Afroza Sultana, Akher Ali, Sifat Ar Salan and Mohammad Alamgir Kabir in Journal of Public Health Research

Footnotes

Acknowledgements

The authors are thankful to the academic editor, associate editor and three anonymous reviewers for their valuable comments and suggestions that helped to enhance the quality of the manuscript.

Abbreviations

WAIC: Watanabe-Akaike Information Criterion, DIC: Deviance Information Criterion, MCMC: Markov Chain Monte Carlo, AR(2): Autoregressive of order two, ANOVA: Analysis of Variance, CAR: Conditional Autoregressive, LMPL: Log Marginal Predictive Likelihood, RMSE: Route Mean Square Error, MAE: Mean Absolute Error, GLM: Generalized Linear Model, IEBGLM: Independent Error Bayesian GLM, AR(1): Autoregressive of order one, p.d: Penalty of deviance information criteria, p.w: Watanabe Akaike Information Criteria, CRPS: Continuous Ranked Probability Score, CVG: Convergence.

ORCID iD

Afroza Sultana

Author Contributions

Afroza Sultana: Conceptualization, Data curation, Formal analysis, Methodology, Software, Visualization, Investigation, Writing – original draft. Md. Akher Ali: Data curation, Visualization, Writing original draft. Md. Sifat Ar Salan: Conceptualization, Formal analysis, Methodology, Software, Supervision, Investigation, Writing original draft, Validation, Writing, Reviewing and editing. Mohammad Alamgir Kabir: Conceptualization, Methodology, Supervision, Validation, Writing, Reviewing and editing. All authors approved the final version of the manuscript.

Funding

The author(s) received no financial support for the research, authorship, and/or publication of this article

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article

Code sharing statement

The R-code used in this study will be made available to the first author upon reasonable request

Supplemental material

Supplemental material for this article is available online.

References

von Schirnding

. The World Summit on Sustainable Development: reaffirming the centrality of health. Global Health 2005; 1(1): 8.

Akther

Bangladesh’s journey towards SDG 4: a review of progresses. Int J Multidiscip Inf Res Rev 2022; 2(1): 43–51.

Khan

II . Education and development in Bangladesh: a study from spatial perspective. J Bangladesh Inst Plann 2010; 3(1): 1–16.

Jha

16. United Nations General Assembly (UNGA). Yearb Int Environ Law 2016; 27: 506–514.

Rahman

MM.

Achieving Sustainable Development Goals in Bangladesh: an organizational analysis. SSRN Electron J 2020. DOI: 10.2139/ssrn.3779081.

Fayissa

Gutema

Estimating a health production function for Sub-Saharan Africa (SSA). Appl Econ 2005; 37(2): 155–164.

Lei

Liu

, et al. Quantitative study on socioeconomic determinants of life expectancy in Beijing, China. J Évid Med 2009; 2(2): 92–98.

Bayati

Akbarian

Kavosi

, et al. Socioeconomic Determinants of Health in Western Pacific Region: a panel data analysis. Soc Welf Q 2013; 12(47): 111–130.

Kabir

Determinants of life expectancy in developing countries. J Dev Areas 2008; 41(2): 185–204.

10.

Thornton

Estimating a health production function for the US: some new evidence. Appl Econ 2002; 34(1): 59–62.

11.

Auster

Leveson

Sarachek

. The production of health, an exploratory study. In: Essays in the economics of health and medical care. 1972, pp.135–158. http://www.nber.org/books/fuch72-1 (accessed 13 July 2025).

12.

Determinants of health outcomes in industrialised countries: a pooled, cross-country, time-series analysis. OECD Econ Stud 2000; 30: 53–78.

13.

World Bank. Poverty assessment for Bangladesh, https://datatopics.worldbank.org/education/country/bangladesh (2021, accessed 17 October 2023).

14.

WHO. Global action plan on physical activity 2018-2030: more active people for a healthier world. World Health Organization, 2019.

15.

Smith

MJD

Goodchild Michael

Longley

Geospatial analysis: a comprehensive guide to principles, techniques and software tools. Troubador Publishing Ltd, 2007.

16.

GDL. Global Data Lab, https://globaldatalab.org//ddw/ (2021, accessed 17 October 2023).

17.

BBS. Statistical Year book of Bangladesh. Bangladesh Bureau of Statistics, 2021. http://www.bbs.gov.bd (accessed 17 October 2023).

18.

Kabir

Islam

Ahmed

, et al. Factors influencing infant and child mortality in Bangladesh. J Med Sci 2001; 1(5): 292–295.

19.

Bairagi

Sutradhar

Alam

Levels, trends and determinants of child mortality in Matlab, Bangladesh, 1966-1994. Asia-Pacific Population Journal 1999; 14(2): 1–12.

20.

Hobcraft

McDonald

Rutstein

SO.

Socio-economic factors in infant and child mortality: a cross-national comparison. Popul Stud 1984; 38(2): 193–223.

21.

Chen

New Approaches for calculating Moran’s index of Spatial Autocorrelation. PLoS One 2013; 8(7): e68336.

22.

Chen

Spatial autocorrelation approaches to testing residuals from least squares regression. PLoS One 2016; 11(1): e0146865.

23.

Geary

RC.

The contiguity ratio and statistical mapping. Inc Stat 1954; 5(3): 115.

24.

Shen

Spatio-temporal autocorrelation measures for nonstationary series: a new temporally detrended spatio-temporal Moran’s index. Phys Lett 2016; 380(1-2): 106–116.

25.

Sahu

. Hierarchical Bayesian models for space–time air pollution data, Vol. 30. Elsevier, 2012, pp.477–495. https://doi.org/10.1016/B978-0-444-53858-1.00016-8

26.

Knorr-Held

Rasser

Bayesian detection of clusters and discontinuities in disease maps. Biometrics 2000; 56(1): 13–21.

27.

Sahu

. bmstdr: Bayesian Modeling of Spatio-Temporal Data with R. Technical Report, Submitted, 2021. https://www.soton.ac.uk/sks/bmbook/bmstdr-vignette.html

28.

Spiegelhalter

Best

Carlin

, et al. Bayesian measures of model complexity and fit. J R Stat Soc Series B Stat Methodol 2002; 64(4): 583–639.

29.

Watanabe

Asymptotic equivalence of Bayes cross validation and widely applicable information criterion in singular learning theory. J Mach Learn Res 2010; 11: 3571–3595.

30.

Wang

Lam

NSN

. Extending Getis–Ord statistics to account for local space–time autocorrelation in spatial panel data. Prof Geogr 2020; 72(3): 411–420.

31.

Feng

Shamim

, et al. High-resolution mapping of reproductive tract infections among women of childbearing age in Bangladesh: a spatial-temporal analysis of the demographic and health survey. BMC Public Health 2021; 21(1): 342.

32.

Raheem

Khan

Hossain

MS.

Regional disparities in maternal and child health indicators: cluster analysis of districts in Bangladesh. PLoS One 2019; 14(2): e0210697.

33.

Schoder

Zweifel

Flat-of-the-curve medicine: a new perspective on the production of health. Health Econ Rev 2011; 1(1): 1–10.

34.

Nandi

Hossain

Roy

, et al. An investigation of the relation between life expectancy & socioeconomic variables using path analysis for Sustainable Development Goals (SDG) in Bangladesh. PLoS One 2023; 18(2): e0275431.

35.

Congdon

Applied bayesian modelling. 2nd ed. John Wiley & Sons, 2014.

36.

Ghosh

Natarajan

Stroud

TWF

, et al. Generalized linear models for small-area estimation. J Am Stat Assoc 1998; 93(441): 273–282.

37.

Das

Mia

Hanifi

SMA

, et al. Health literacy in a community with low levels of education: findings from Chakaria, a rural area of Bangladesh. BMC Public Health 2017; 17(1): 203.

38.

Gulcan

A study on life expectancy in Turkey. EMU J Pharm Sci 2020; 2(3): 86–96.

39.

Akter

Factors influencing health service utilization among mothers for under-five children: a cross-sectional study in Khulna district of Bangladesh. PLoS One 2022; 17(9): e0274449.

40.

Kabir

MR.

Adopting Andersen’s behavior model to identify factors influencing maternal healthcare service utilization in Bangladesh. PLoS One 2021; 16(11): e0260502.

41.

Moazzem

Shibly

ASA

. Power sector in the 8th five year plan: reflection on its strategy and initiatives (138; CPD Working Paper), 2021.

42.

Agarwal

Curtis

Angeles

, et al. The impact of India’s accredited social health activist (ASHA) program on the utilization of maternity services: a nationally representative longitudinal modelling study. Hum Resour Health 2019; 17(1): 68.

43.

Khatri

Mishra

Khanal

Female Community Health Volunteers in community-based health programs of Nepal: Future perspective. Front Public Health 2017; 5: 181. DOI: 10.3389/fpubh.2017.00181.

44.

Sheikh

MD.

Challanges and prospects of community-based health insurance in Bangladesh. Doctoral Thesis. University of Strathclyde, 2024.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

0.03 MB