Skip to main content
Advertisement
  • Loading metrics

Joint modelling of systolic and diastolic blood pressure and its associated factors among women in Ghana: Multivariate response multilevel modelling methods

  • Justice Moses K. Aheto ,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    justiceaheto@yahoo.com, jmkaheto@ug.edu.gh

    Affiliations Department of Biostatistics, School of Public Health, College of Health Sciences, University of Ghana, Accra, Ghana, WorldPop, School of Geography and Environmental Science, University of Southampton, Southampton, United Kingdom, College of Public Health, University of South Florida, Tampa, Florida, United States of America

  • Tracy Gates,

    Roles Data curation, Writing – original draft, Writing – review & editing

    Affiliation College of Public Health, University of South Florida, Tampa, Florida, United States of America

  • Rahmatu Babah,

    Roles Data curation, Writing – original draft, Writing – review & editing

    Affiliation Department of Biostatistics, School of Public Health, College of Health Sciences, University of Ghana, Accra, Ghana

  • Wisdom Takramah

    Roles Writing – original draft, Writing – review & editing

    Affiliation Department of Biostatistics, School of Public Health, College of Health Sciences, University of Ghana, Accra, Ghana

Abstract

Elevated blood pressure is the leading cause of cardiovascular diseases related mortality and a major contributor to non-communicable diseases globally, especially in sub-Saharan Africa where about 74.7 million people live with hypertension. In Ghana, hypertension is epidemic with prevalence of over 30% and experiencing continuing burden with its associated morbidity and mortality. Using the 2014 Ghana Demographic and Health Survey, we analyzed data on 4744 women aged 15–49 years residing in 3722 households. We employed univariate and multivariate response multilevel linear regression models to analyze predictors of systolic blood pressure (SBP) and diastolic blood pressure (DBP). Geospatial maps were produced to show the regional distribution of hypertension prevalence in Ghana. Stata version 17 and R version 4.2.1 were used to analyze the data. Of the 4744 woman, 337 (7.1%) and 484 (10.2%) were found to be hypertensive on SBP and DBP, respectively. A combined prevalence of 12.3% was found. Older ages 25–34 (OR 2.45, 95%CI: 1.27, 3.63), 35–44 (OR 8.72, 95%CI: 7.43, 10.01), 45–49 (OR 15.85, 95%CI: 14.07, 17.64), being obese (OR 5.10, 95%CI: 3.62, 6.58), and having no education (OR -2.05, 95%CI: -3.40, -0.71) were associated with SBP. For DBP, we found the associated factors to be older ages 25–34 (OR 3.29, 95%CI: 2.50, 4.08), 35–44 (OR 6.78, 95%CI: 5.91, 7.64), 45–49 (OR 10.05, 95%CI: 8.85, 11.25), being obese (OR 4.20, 95%CI: 3.21, 5.19), and having no education (OR -1.23, 95%CI: -2.14, -0.33). Substantial residual household level differences in SBP (15%) and DBP (14%) were observed. We found strong residual correlation of SBP and DBP on individual women (r = 0.73) and household-level (r = 0.81). The geospatial maps showed substantial regional differences in the observed and reported hypertension prevalence. Interventions should be targeted at the identified high-risk groups like older age groups and those who are obese, and the high-risk regions.

Introduction

The leading cause of cardiovascular disease and premature mortality worldwide, hypertension is the “silent killer” affecting more than a billion people, the majority of whom live in low- and middle-income countries (LMICs) [1, 2]. Hypertension occurs when blood pressure is too high, thus exerting too much force on artery walls leading to a reduction in the critical flow of blood and oxygen throughout the body [3]. Hypertension is a significant risk factor for cardiovascular diseases, including stroke, myocardial infarction, heart failure, and atherosclerosis, which account for most of the non-communicable disease deaths globally [1, 2, 46]. Blood pressure consists of two measurements; systolic blood pressure which measures the pressure in the blood vessels when the heart contracts or beats while diastolic pressure measures the pressure in the vessels when the heart rests between beats [3]. Hypertension is diagnosed when a person’s systolic blood pressure is ≥140mmHg, or diastolic blood pressure ≥90mmHg [3].

Termed the “silent killer”, hypertension is often asymptomatic, resulting in many people being unaware of their condition [1, 3, 7]. Therefore, hypertension diagnosis from blood pressure measurements is critical for people to know their status, receive necessary treatment, and adapt important lifestyle changes. Diagnosis must be supplemented with education and accessibility of treatment. In 2019, it was estimated that less than half of hypertension cases in sub-Saharan Africa were diagnosed, only 22% of men and 45% of women with hypertension were treated, and only 9% of men and 13% of women had their condition under control [4]. These critical gaps in diagnosis and effective treatment must be addressed to reduce the millions of preventable deaths and years of life lost annually.

Because of its global prevalence, detrimental outcomes, and preventability, the reduction of hypertension has been prioritized around the world leading to significant improvements in most countries [2, 3]. However, little improvements have been made in sub-Saharan Africa despite Africa being home to the highest hypertension prevalence rates globally [24, 7]. A significant contributor to the limited progress in addressing hypertension in sub-Saharan Africa is the reality of competing interests for limited resources posed by highly visible communicable diseases such as malaria, HIV/AIDS, and tuberculosis [8]. Studies in Ghana, an Anglophone country in sub-Saharan Africa, have found increasing prevalence of hypertension and associated cardiovascular morbidity and mortality despite the initiation of numerous preventative efforts including the National Health Insurance Scheme, national policy for the Prevention and Control of Non-Communicable Diseases, community-based hypertension improvement program, and May Measurement Months [1, 2, 9]. According to a meta-analysis conducted in 2020, the pooled prevalence of hypertension in Ghana is greater than 30% (30.1% for females and 34.0% for males) [8]. Additionally, it is estimated that the diagnosis rate of hypertension in Ghana is between 16.4% and 54.1%, with only 1.7% to 12.7% of hypertension cases under control [10].

Because of the high prevalence rates, and low diagnosis and control rates, it is critical to understand the numerous factors influencing hypertension, and their relationship to each other, in order to enhance prevention efforts by addressing specific concerns in targeted locations. Hypertension is induced through a conglomeration of effects within a complex network of predisposing, enabling, and reinforcing factors. Many studies have found that predisposing factors contributing to hypertension include genetics, sex, race, age, marital status, type and place of residence, weight, pregnancy, health insurance, employment, educational attainment, and lack of understanding of hypertension [13, 5, 7, 8, 1113]. Enabling factors include alcohol consumption, sedentary lifestyle, tobacco usage, diet, and socioeconomic status [13, 5, 7, 8, 11]. Factors reinforcing the presence of hypertension in Ghana include the limited availability and accessibility of healthcare, shortcomings of the medical and insurance system, cultural perceptions, favorable to being overweight, and cost of healthy lifestyles (including diet and physical activity) [1, 2, 5, 8, 11, 12]. Understanding the influence of these demographic and socio-ecological factors on hypertension is critical to addressing them, reducing prevalence, and improving treatment outcomes.

Enhancing the understanding of local risk factors for hypertension is critical to addressing the deadly hypertension epidemic in Ghana. Previous studies have analyzed risk factors of hypertension in Ghana using univariate models [1, 5, 12, 14]; however, no research has yet analyzed these factors using multivariate response multilevel regression modeling approaches to effectively account for the effects of clustering and multilevel interactions. To support targeted and optimal hypertension policy and intervention strategies by policymakers and program managers, there is a need to employ a more complex model that considers maternal, household, and community factors simultaneously, while investigating their interactions across multiple levels.

To this end, we attempt to develop a sophisticated multivariate response multilevel model to jointly analyze systolic blood pressure (SBP) and diastolic blood pressure (DBP) among women in Ghana to understand the complexity of the data to be able to target interventions to improve prevention, diagnosis, and treatment efforts. The multivariate response multilevel modelling is critical in this study because we are interested in modelling both SBP and DBP outcomes measured on same woman nested within a household and we need to draw conclusions about the degree to which the residual correlations depend on these women and their households. We are also interested in quantifying unobserved household level differences in SBP and DBP outcomes, and to investigate specific effect of a given risk factor across these outcomes simultaneously. Also, the multivariate response multilevel modelling approach will permit us to conduct a single test of a joint effect of a risk factor on these two outcomes without multiplying the effect of Type I error which will lead to higher accuracy and reliability in estimates relative to modelling these outcomes separately, especially when these outcomes are at least moderately correlated based on the empirical data [15, 16].

Materials and methods

Data source and study population

The data for the current study was extracted from the Ghana demographic and health survey (GDHS) database. The Ghana demographic and health surveys are nationally representative cross-sectional population-based surveys that are carried out every five years by Ghana Statistical Service (GSS) in collaboration with United States Agency for International Development (USAID) and the demographic and health survey (DHS) program. The Ghana Statistical Service has been mandated to implement demographic and health surveys by collecting, analyzing, and providing accurate information on nutrition, malaria, maternal and child health, fertility, gender, and family planning.

The study design used in demographic and health surveys involves a two-stage sampling technique. The enumeration areas (EA) were selected at random using the probability proportional to size technique at the first stage. The second stage involves the selection of households from the EAs using a systematic sampling technique. The three main DHS questionnaires comprise the household questionnaire, women’s questionnaire, and man’s questionnaire. The sampling frame used in the GDHS was obtained from the 2010 Ghana Population and Housing Census (PHC) which were commissioned by Ghana Statistical Service. The total number of primary sampling unit or enumeration areas (EAs) selected for the 2014 GDHS was 427. A total of 9396 women aged 15–49 years and 4388 men aged 15–59 drawn from 11,385 occupied households were interviewed [17]. The current study extracted biometrics data from the women’s questionnaire where systolic blood pressure and diastolic blood pressure were measured on the women. Detailed description of the methods and procedures utilized in the survey are available elsewhere [17]. In the present study, we extracted data on 4744 women residing in 3727 households.

Outcome variable

The current study considered systolic blood pressure and diastolic blood pressure as the two outcome variables. These continuous outcome variables are measured in units of millimeters of mercury (mmHg).

Covariates

The selection of the covariates for the current study was based on the variables used in similar studies and on their availability in the GDHS database [5, 14, 18]. Thus, the current study considered factors such as current age, vegetable consumption, household size, cooking fuel, water source, toilet type, education, wealth, wall type, and obesity.

Statistical analysis

Data cleaning, recording and validation were done in STATA version 17. We examined any possible presence of outliers in the quantitative variables using boxplot. We also performed checks on the categorical variables to ensure that the categories were mutually exclusive and exhaustive using a simple frequency distribution. Descriptive statistics were performed to describe the characteristics of the study sample.

Previous studies only applied multilevel model to explore any possible individual-level and household-level variations in these outcomes of interest separately [14, 19, 20]. The inferential statistic in the current study focuses on joint modeling of SBP and DBP and their associated factors within a multilevel framework. Clinically, measurement of high blood pressure/hypertension is based on the combination of systolic blood pressure and diastolic blood pressure. Even though we could specify two separate regression models to estimate the correlates of systolic blood pressure and diastolic blood pressure as in previous studies [14, 19, 20], a reviewed literature indicated the computational advantage of joint modeling of outcomes. Our data emanated from a complex survey where individuals were nested in households and given that we are modelling SBP and DBP jointly, we developed and applied the multivariate response multilevel models [15, 16] as the most appropriate statistical method to analyze the data to answer the study objectives. Thus, we jointly model the SBP and DBP of individual women nested within households while simultaneously adjusting for potential risk factors.

Model description

Let yij, i = 1,…,N, j = 1,…,nij, be the response for the ith subject nested in jth household. Let Y1ij represents systolic blood pressure for the ith subject nested in jth household (SBP), Y2ij represents diastolic blood pressure (DBP) for the ith subject nested in jth household. Thus, the equations specified to jointly model Y1ij and Y2ij continuous data are defined by [11, 12]: (1) where Xij denotes the covariate vectors of fixed effects for the outcome variables, β1 and β2 denote vector of regression coefficients for SBP and DBP respectively, h1j and h2j denotes the random effects at the household level for SBP and DBP respectively, ε1ij and ε2ij represent the error terms or residuals at the individual level for SBP and DBP respectively.

The joint distribution for the random effects that combines both response trajectories is given by: (2) where G denotes the structure of the variance-covariance matrix for the random effect and is defined as: (3)

The residual correlation [15, 16, 21, 22] between systolic blood pressure (SBP) and diastolic blood pressure (DBP) at the household-level (AEH) and women-level (AEW) were estimated using the covariance matrix of the random effects obtained from the multivariate response multilevel regression models.

For the households, we have: (4) and for the individual women, we have: (5)

The intraclass correlation coefficient (ICC) which coincide with the variance partition coefficient (VPC) under the random intercept model provides information on the proportion of total variance in the outcome variable that could be attributable to the household level [21, 23, 24]. Thus, the VPC was estimated for systolic blood pressure (SBP) and diastolic blood pressure (DBP) at the household level using the equations below: (6) (7) where and denote the variance of the random effects at the household-level h1j and h2j respectively, and denote the variance of the residuals at the individual level ε1ij and ε2ij respectively.

The joint and separate models were evaluated using Akaike Information criterion (AIC), Bayesian Information criterion (BIC) and likelihood ratio tests (LRT). The model with the lowest AIC and BIC score was selected as the most preferred model. Also, a significant test on LRT was used to select the preferred model. We examined the residuals for multivariate normality assumption using quantile-quantile plots. The runmlwin [25] program that was designed to run MLwiN multilevel models from within Stata was used to fit both the univariate response multilevel and the multivariate response multilevel regression models in Stata version 17.

Geospatial mapping of the observed and reported hypertension prevalence by regions in Ghana

We downloaded and extracted regional shapefiles for Ghana using Geographic Information Systems (GIS) and geospatial approaches. The shapefiles for the 10 regions that were in existence at the time of the survey were obtained from the GADM website freely available at https://gadm.org/download_country40.html and terms of use/license available at https://gadm.org/license.html. This will allow geospatial mapping of the observed hypertension prevalence based on SBP, DBP, and those ever told to have had hypertension. The packages rgdal, tmap, and leaflet in R version 4.2.1 were used to prepare the data and for the geospatial mapping.

Ethics approval and consent to participate

The protocol for the 2014 GDHS was reviewed and approved by the Ghana Health Service Ethical Review Committee and the Institutional Review Board of ICF International. An informed consent statement was read to the respondent, who accepted or declined to participate. The study obtained consent to participate from all participants, and a parent or guardian provided consent prior to participation by a child or adolescent [17]. Though this study did not require ethical approval because the GDHS dataset is publicly available and this study was granted permission to use the data by the MEASURE DHS Program, it is worth mentioning that ethical clearance procedures were strictly adhered to before, during and after the GDHS to ensure protection of human rights as required by the US Department of Health and Human Services. Detailed information on demographic and health surveys (DHS) data and ethical standards are publicly available at http://goo.gl/ny8T6X.

Results

Summary of background characteristics

We analyzed data on a total of 4744 women residing in 3727 households out of which 7.1% and 10.2% were found to be hypertensive on SBP and DBP, respectively with a combined (i.e., SBP and/or DBP) prevalence of 12.3%, and 6.9% of the women said they were ever told to have had hypertension. Majority were aged 15–24 years (35.6%) while 12.8% of the participants were obese. A total of 2688 (56.7%) had secondary education or higher and majority of them belonged to poorer/poorest households while 3362 (70.9%) reside in households with more than 3 household members. About 14% of them did not eat vegetables in the past 7 days prior to the survey (Table 1).

thumbnail
Table 1. Background characteristics of the study participants.

https://doi.org/10.1371/journal.pgph.0001613.t001

Factors associated with SBP and DBP in the univariate response multilevel models

In the univariate models, we found age, obesity, and education to be associated with both SBP and DBP. Type of cooking fuel was only associated with DBP while toilet type was only associated with SBP. Older ages and being obese were associated with increased SBP and DBP while having no education was associated with decreased SBP and DBP. Use of primitive fuel for cooking was associated with decreased DBP and having unimproved toilet in households was associated with increased SBP. We found substantial residual household level differences in SBP and DBP of 15% and 14%, respectively (Table 2).

thumbnail
Table 2. Parameter estimates for separate univariate response multilevel models of SBP and DBP.

https://doi.org/10.1371/journal.pgph.0001613.t002

Factors associated with SBP and DBP in the multivariate response multilevel models

To begin with, we estimated the empirical correlation between SBP and DBP, and determined its statistical significance. We found a strong positive and statistically significant correlation of 0.76 (p<0.0001) between SBP and DBP. Furthermore, we fitted a single level (i.e., women level) multivariate response model and compared this to the multivariate response multilevel model (i.e., household and women level). We observed that the multivariate response multilevel model provided a better fit to the data compared to the single level multivariate response model (χ2 = 36.6, p<0.0001), suggesting substantial unobserved household level variations in SBP and DBP based on the multivariate response multilevel regression model after adjusting for the same set of risk factors considered in both models.

The results from the multivariate response multilevel regression model are presented in Table 3. Generally, the results are similar to those presented in Table 1 except the quantification of the residual correlation between SBP and DBP on individual women and their households and testing the impact of a given risk factors across SBP and DBP which is not possible under the separate model results presented in Table 1. We found that older ages 25–34 (OR 2.45, 95%CI: 1.27, 3.63), 35–44 (OR 8.72, 95%CI: 7.43, 10.01), and 45–49 (OR 15.85, 95%CI: 14.07, 17.64), being obese (OR 5.10, 95%CI: 3.62, 6.58), having no education (OR -2.05, 95%CI: -3.40, -0.71), and unimproved toilet facility (OR 1.75, 95%CI: 0.56, 2.93) were associated with SBP. Furthermore, older ages 25–34 (OR 3.29, 95%CI: 2.50, 4.08), 35–44 (OR 6.78, 95%CI: 5.91, 7.64), 45–49 (OR 10.05, 95%CI: 8.85, 11.25), being obese (OR 4.20, 95%CI: 3.21, 5.19), use of primitive cooking fuel (OR -2.06, 95%CI: -3.40, -0.72), and having no education (OR -1.23, 95%CI: -2.14, -0.33) were associated with DBP (Table 3).

thumbnail
Table 3. Parameter estimates for multivariate response multilevel models of SBP and DBP.

https://doi.org/10.1371/journal.pgph.0001613.t003

Substantial residual household level differences in SBP and DBP were observed where 15% and 14% of the differences observed in SBP and DBP respectively could be attributable to unobserved household level differences after adjusting for the risk factors considered in our models. Furthermore, the residual correlation on individual women and household-level effects for SBP and DBP were 0.73 and 0.81 respectively, suggesting strong residual correlation of these outcomes on women and their households (Table 4). The residual correlation measures whether after accounting for the risk factors in our model, the factors yet to be identified at the women, household or community levels that predict SBP and DBP are the same. Thus, the higher the residual correlation, the more similar these as-yet unidentified factors that could influence both SBP and DBP will be and vice-versa.

thumbnail
Table 4. Analysis of variance and residual correlation estimates for the multivariate response multilevel models of SBP and DBP.

https://doi.org/10.1371/journal.pgph.0001613.t004

In Table 5, we tested the hypothesis that the effect of a risk factor on SBP is the same as on DBP. Our tests showed that the effects of women’s age, and toilet type on SBP were significantly different from that of DBP, but no such differences exist for obesity, cooking fuel and women education.

thumbnail
Table 5. Test for equality of risk factor effect across SBP and DBP.

https://doi.org/10.1371/journal.pgph.0001613.t005

Geospatial mapping of observed and reported hypertension prevalence

We found remarkable differences in the prevalence of hypertension among women across the regions of Ghana. Based on SBP, the hypertension prevalence was highest in Greater Accra (15.7%) followed by Ashanti (11.6%) and Eastern (11.6%) regions with Upper West recording the lowest prevalence of 5.3%. For DBP, the prevalence was highest in Greater Accra (17.2%) followed by Western (13.8%) with the lowest prevalence of 5.2% recorded in Upper West region. For the combined prevalence based on both SBP and DBP, Greater Accra recorded the highest (16.8%) followed by Western (12.8%) with Upper West recording the lowest (5.1%). For those who were ever told to have had hypertension, majority of them were from Greater Accra (21.8%) followed by Volta (12.3%) with minority coming from the Upper West region (4.0%) (Figs 14).

thumbnail
Fig 1. Regional distribution of observed hypertension prevalence based on systolic blood pressure measurements among women in 2014 in Ghana.

Source: This map was produced by the authors. The link to download the shapefile is available at https://gadm.org/download_country40.html.

https://doi.org/10.1371/journal.pgph.0001613.g001

thumbnail
Fig 2. Regional distribution of observed hypertension prevalence based on diastolic blood pressure measurements among women in 2014 in Ghana.

Source: This map was produced by the authors. The link to download the shapefile is available at https://gadm.org/download_country40.html.

https://doi.org/10.1371/journal.pgph.0001613.g002

thumbnail
Fig 3. Regional distribution of observed hypertension prevalence based on both systolic and diastolic blood pressure measurements among women in 2014 in Ghana.

Source: This map was produced by the authors. The link to download the shapefile is available at https://gadm.org/download_country40.html.

https://doi.org/10.1371/journal.pgph.0001613.g003

thumbnail
Fig 4. Regional distribution of reported hypertension prevalence among women who were ever told they had hypertension in Ghana.

Source: This map was produced by the authors. The link to download the shapefile is available at https://gadm.org/download_country40.html.

https://doi.org/10.1371/journal.pgph.0001613.g004

Discussion

In this study, we set out to jointly analyze SBP and DBP under multivariate response multilevel modelling methods while simultaneously adjusting for critical risk factors, the first of its kind in this area. Furthermore, we produced geospatial maps to show differences in regional spatial distribution of observed hypertension prevalence based on SBP, DBP, both SBP and DBP (i.e., combined), and reported hypertension. The hypertension prevalence based on SBP and DBP were 7.1% and 10.2% respectively, while those who reported ever been told to have hypertension, and hypertension based on the combined prevalence of SBP and/or DBP were respectively 12.3%, and 6.9%. Substantial regional differences in the hypertension prevalence were found for SBP, DBP, the combined prevalence and the reported. The national hypertension prevalence found in this study ranged from 5.1% in the Upper West region to 16.8% in the Greater Accra region. The prevalence based on SBP ranged from 5.3% in the Upper West region to 15.7% in the Greater Accra region. For DBP, the prevalence ranged from 5.2% in the Upper West region to 17.2% in the Greater Accra region and for those who reported that they were ever told to have had hypertension, the prevalence ranged from 4% in the Upper West region to 21.8% in the Greater Accra region. Thus, Upper West and Greater Accra regions consistently recorded the lowest and the highest hypertension prevalence respectively, on all the four outcome measures. The national prevalence based on SBP and DBP, and a combined (i.e., SBP and/or DBP) prevalence among this group of women aged 15–49 in Ghana was lower than the combined prevalence of 20% reported for women aged 15–49 years in Lesotho [18], 41.4% in Guatemala [26], and 33.8% in Ghana [12].

Of particular interest in this study is the quantification of the residual correlation between SBP and DBP on the individual women and their households in which they reside. This residual correlation measures whether after accounting for the risk factors in our model, the factors yet to be identified to predict SBP and DBP are similar or the same. This is critical because one cannot assume that two different outcomes (SBP and DBP in our case) collected on the same woman at the same time residing in a household are not correlated or similar on the same woman and her household. Ignoring of this similarity in these outcomes as in the fitting of separate multilevel models for these outcomes could lead to imprecision in model parameters [15, 16, 21]. This is one of the innovations introduced in this study. We found strong residual correlation between SBP and DBP on the individual women and their households, an indication that the as-yet unidentified factors could be predictive of both SBP and DBP. Thus, an intervention being rollout that aimed at identifying other factors to further understand hypertension prevalence based on SBP could similarly address the problem of hypertension based on DBP without the need to rollout separate interventions for SBP and DBP. Also, we are interested in the quantification of the unobserved household level differences in SBP and DBP outcomes which represents variations in household level SBP and DBP outcomes that cannot be explained by the available risk factors considered in our models. It is a common knowledge that households constitute critical predictors of socioeconomic differences in health and overall well-being of individuals because they influence their opportunities and govern exposure to various opportunities, risks, and resources over their life course [14, 27, 28]. Significant unobserved household-level differences in SBP and DBP among the women were found, suggesting that there are unanswered questions about why a woman residing in a certain household is likely to have increased SBP and DBP compared to a woman residing in another household or vice versa. This finding is consistent with a previous study that found substantial residual household level variations in hypertension prevalence [14].

Critical risk factors identified to be associated with both SBP and DBP were age, obesity, and education. Specifically, older ages of women and being obese were predictive of increased SBP and DBP, and women with no education had decreased SBP and DBP. Also, having unimproved toilet type in households was associated with increased SBP while using primitive fuel for cooking in households was associated with decreased DBP. Our finding that older ages was associated with increased SBP and DBP is not surprising because generally the older the age, the higher the chance of one developing hypertension and this is consistent with previous studies that found that high blood pressure increase with an increasing age [12, 13, 18, 2931]. Consistent with previous studies that reported increase in hypertension prevalence among those who are obese [12, 13, 30, 31], being obese was associated with increased SBP and DBP in this study. Obesity and age were reported to be the most powerful factors associated with hypertension in previous studies [12, 18, 29, 30]. A previous study [31] reported that increase in years of education reduced hypertension prevalence which is not consistent with our finding that women with no education had decreased SBP and DBP. This inconsistency could be attributable to the fact that the study populations are not that similar and that our modelling procedures are more sophisticated with continuous outcomes compared to dichotomous outcome for hypertension used in the previous study [31]. The finding that using primitive fuel for cooking in households was associated with decreased DBP is inconsistent with previous studies which found that biomass fuels for cooking was associated with higher DBP, SBP, and higher risk of hypertension [32, 33]. This could be due to differences in study settings and statistical modelling approaches. Also, there is lack of available literature on the mechanism through which unimproved toilet facility was associated with increased SBP. Further studies are therefore recommended to identify the possible pathways. The findings also revealed that the effect of women’s age, and toilet type differ significantly across SBP and DBP outcomes.

We addressed a critical public health problem facing women of reproductive age, especially in the presence of high maternal mortality in the country as previous studies have reported a link between hypertension versus maternal mortality and morbidity [18]. To improve the general cardiovascular health outcome of women, our modelling approach and the identified risk factors should be considered as part of all overall strategy aimed at addressing hypertension and the general cardiovascular health outcomes in this population. The association of both outcomes (i.e., SBP and DBP) with obesity and type of toilet suggest the need for modification of lifestyle in the form of dietary intake and dietary knowledge and education supported with strong public health nutrition messages, and general improvement in socioeconomic conditions of households. Older women should be targeted with interventions aimed at controlling their high blood pressure while at the same time targeting younger women with interventions aimed at preventing the risk of developing hypertension in their old ages.

Strengths and limitations of the study

The study is a population-level based with national coverage, making our findings to be readily generalized to the population of women within reproductive age in Ghana and other similar countries. Our modelling approach is very novel and captures residual similarity between SBP and SBP on individual women and their households while simultaneously quantifying unobserved household level differences in these outcomes after adjusting for selected risk factors, leading to more precise estimates in our parameters. In addition, the modelling approach helps avoid multiplication of Type I errors by testing the effect of a given risk factor across both outcomes simultaneously. Despite these strengths, our finding should be interpreted with caution. We utilized data which is cross-sectional and could not establish causative effects. As with all studies, we could not account for all potential risk factors because the DHS data did not collect data on additional factors such as data on policies implemented in Ghana to address the problem of hypertension, including other individual, household and community level factors.

Conclusion

The current study employed novel statistical methods, first of its kind in this area to capture residual similarity between SBP and DBP on women and their households, and residual household level variations in these outcomes, and identified critical risk factors. We addressed the statistical problem of ignoring correlation between SBP and DBP outcomes on the same woman and the household she lives which helps avoid imprecision in our model parameters, leading to better statistical inference and sound policy decisions. Substantial residual correlation between SBP and DBP on women and their households were found, while significant unobserved household level differences in SBP and DBP were observed. Women’s age, obesity, education, toilet type and type of cooking fuel were found to be predictive of SBP and/or DBP. Our findings that substantial residual correlation between SBP and DBP on women and their households, and significant unobserved household level differences in SBP and DBP outcomes warrant further research to identify as-yet unidentified predictors of SBP and DBP outcomes in this population of women. There is also the need to conduct qualitative studies to establish why women from certain households are more likely to have increased SBP and DBP compared to women from other households and vice versa. We encourage researchers working on hypertension predictors and the role of the higher-level hierarchy effects (e.g., household and community group levels) to jointly model SBP and DBP outcomes under multivariate response multilevel modelling framework as utilized in this study.

Acknowledgments

Special thanks go to DHS program for providing access to the datasets used in the current study.

References

  1. 1. Tannor EK, Nyarko OO, Adu-Boakye Y, Owusu Konadu S, Opoku G, Ankobea-Kokroe F, et al. Burden of hypertension in Ghana—Analysis of awareness and screening campaign in the Ashanti Region of Ghana. JRSM Cardiovascular Disease. 2022;11:20480040221075521. Epub 20220128. pmid:35111304; PubMed Central PMCID: PMC8801650.
  2. 2. Mills KT, Stefanescu A, He J. The global epidemiology of hypertension. Nature Reviews Nephrology. 2020;16(4):223–37. pmid:32024986
  3. 3. World Health Organization. Hypertension. WHO. 25 August 2021. [Accessed on 28th September 2022 at https://www.who.int/news-room/fact-sheets/detail/hypertension].
  4. 4. Zhou B, Carrillo-Larco RM, Danaei G, Riley LM, Paciorek CJ, Stevens GA, et al. Worldwide trends in hypertension prevalence and progress in treatment and control from 1990 to 2019: a pooled analysis of 1201 population-representative studies with 104 million participants. The Lancet. 2021;398(10304):957–80. pmid:34450083
  5. 5. Konlan KD, Lee H, Lee M, Kim Y, Lee H, Abdulai JA. Risk factors associated with the incidence and prevalence of hypertension in Ghana: an integrated review (2016–2021). International Journal of Environmental Health Research. 2022:1–16. pmid:35578383
  6. 6. World Health Organization. Noncommunicable diseases. WHO. 16 September 2022. [Accessed on 8th October 2022 at https://www.who.int/news-room/fact-sheets/detail/noncommunicable-diseases].
  7. 7. Dosoo DK, Nyame S, Enuameh Y, Ayetey H, Danwonno H, Twumasi M, et al. Prevalence of Hypertension in the Middle Belt of Ghana: A Community-Based Screening Study. International Journal of Hypertension. 2019;2019:1–7. pmid:31687204
  8. 8. Atibila F, Hoor Gt, Donkoh ET, Wahab AI, Kok G. Prevalence of hypertension in Ghanaian society: a systematic review, meta-analysis, and GRADE assessment. Systematic Reviews. 2021;10(1). pmid:34364395
  9. 9. Laar AK, Adler AJ, Kotoh AM, Legido-Quigley H, Lange IL, Perel P, et al. Health system challenges to hypertension and related non-communicable diseases prevention and treatment: perspectives from Ghanaian stakeholders. BMC Health Services Research. 2019;19(1). pmid:31615529
  10. 10. Laar AK, Adler AJ, Prieto-Merino D, Der RMM, Mangortey D, Dirks R, et al. An integrated community-based hypertension improvement program improves hypertension awareness in Ghana. Journal of Global Health Science. 2021;3(2).
  11. 11. Oyekale AS. Effect of Obesity and Other Risk Factors on Hypertension among Women of Reproductive Age in Ghana: An Instrumental Variable Probit Model. International Journal of Environmental Research and Public Health. 2019;16(23). pmid:31779087
  12. 12. Acheampong K, Nyamari JM, Ganu D, Appiah S, Pan X, Kaminga A, et al. Predictors of Hypertension among Adult Female Population in Kpone-Katamanso District, Ghana. International Journal of Hypertension. 2019;2019:1876060. pmid:31308975
  13. 13. Kohler IV, Sudharsanan N, Bandawe C, Kohler H-P. Aging and hypertension among the global poor—Panel data evidence from Malawi. PLOS Global Public Health. 2022;2(6):e0000600. pmid:36962748
  14. 14. Aheto JMK, Dagne GA. Multilevel modeling, prevalence, and predictors of hypertension in Ghana: Evidence from Wave 2 of the World Health Organization’s Study on global AGEing and adult health. Health Science Reports. 2021;4(4):e453. pmid:34938897
  15. 15. Snijders TAB, Bosker RJ. Multilevel analysis: an introduction to basic and advanced multilevel modeling. 2nd ed. London: SAGE; 2012.
  16. 16. Thum YM. Hierarchical Linear Models for Multivariate Outcomes. Journal of Educational and Behavioral Statistics; 1997. p. 77–108.
  17. 17. Ghana Statistical Service GSS, Ghana Health Service GHS, International ICF. Ghana Demographic and Health Survey 2014. Rockville, Maryland, USA: GSS, GHS, and ICF International, 2015.
  18. 18. Lebuso M, De Wet- Billings N. The prevalence and socio-demographic correlates of hypertension among women (15–49 years) in Lesotho: a descriptive analysis. BMC Public Health. 2022;22(1):562. pmid:35317759
  19. 19. Mt Shakeri, Khajavian N, Esmaily H, Yusefi R, Navipour E. Multilevel regression modeling for risk factors associated with hypertension. Navid No. 2019;22(69):25–33.
  20. 20. Le C, Jun D, Yichun L, Zhankun S, Keying Z. Multilevel analysis of the determinants of pre-hypertension and hypertension in rural southwest China. Public Health Rep. 2011;126(3):420–7. pmid:21553671; PubMed Central PMCID: PMC3072865.
  21. 21. Aheto MJK. Modelling Malnutrition among Under-Five- Year-Old Children in Ghana. Lancaster University, United Kingdom: Lancaster University; 2016. https://eprints.lancs.ac.uk/id/eprint/82143
  22. 22. Gao F, Philip Miller J, Xiong C, Luo J, Beiser JA, Chen L, et al. Estimating correlation between multivariate longitudinal data in the presence of heterogeneity. BMC Medical Research Methodology. 2017;17(1):124. pmid:28818061
  23. 23. Aheto JMK, Keegan TJ, Taylor BM, Diggle PJ. Childhood Malnutrition and Its Determinants among Under-Five Children in Ghana. Paediatric and Perinatal Epidemiology. 2015;29(6):552–61. pmid:26332093
  24. 24. Steele F, Goldstein H, Rao C, Sinharay S. Multilevel Models in Psychometrics. Psychometrics. 2007;26:401–20. WOS:000311036800013.
  25. 25. Leckie G, Charlton C. runmlwin—A Program to Run the MLwiN Multilevel Modelling Software from within Stata. Journal of Statistical Software. 2013;52(11):1–40.
  26. 26. Pickens CM, Flores-Ayala R, Addo OY, Whitehead RD Jr, Palmieri M, Ramirez-Zea M, et al. Prevalence and Predictors of High Blood Pressure Among Women of Reproductive Age and Children Aged 10 to 14 Years in Guatemala. Preventing Chronic Disease. 2020;17:E66. pmid:32701434
  27. 27. Adekanmbi VT, Kayode GA, Uthman OA. Individual and contextual factors associated with childhood stunting in Nigeria: a multilevel analysis. Maternal and Child Nutrition. 2013;9(2):244–59. pmid:22004134.
  28. 28. Kandala NB, Madungu TP, Emina JB, Nzita KP, Cappuccio FP. Malnutrition among children under the age of five in the Democratic Republic of Congo (DRC): does geographic location matter? BMC Public Health. 2011;11:261. pmid:21518428; PubMed Central PMCID: PMC3111378.
  29. 29. Ibrahim NK, Hijazi NA, Al-Bar AA. Prevalence and Determinants of Prehypertension and Hypertension among Preparatory and Secondary School Teachers in Jeddah. J Egypt Public Health Assoc. 2008;83(3–4):183–203. pmid:19302774.
  30. 30. Ibrahim MM, Damasceno A. Hypertension in developing countries. Lancet. 2012;380(9841):611–9. pmid:22883510.
  31. 31. Ware LJ, Chidumwa G, Charlton K, Schutte AE, Kowal P. Predictors of hypertension awareness, treatment and control in South Africa: results from the WHO-SAGE population survey (Wave 2). J Hum Hypertens. 2019;33(2):157–66. Epub 2018/10/31. pmid:30382179.
  32. 32. Deng Y, Gao Q, Yang D, Hua H, Wang N, Ou F, et al. Association between biomass fuel use and risk of hypertension among Chinese older people: A cohort study. Environment International 2020, 138:105620. pmid:32179315
  33. 33. Arku RE, Ezzati M, Baumgartner J, Fink G, Zhou B, Hystad P, et al. Elevated blood pressure and household solid fuel use in premenopausal women: Analysis of 12 Demographic and Health Surveys (DHS) from 10 countries. Environ Res 2018, 160:499–505. pmid:29107891