key: cord-0025136-pryzftjs
authors: Negesse, Yilkal; Mankelkl, Gosa; Setegn, Melsew; Fetene, Gossa
title: Multilevel analysis of factors associated with HIV among women of reproductive age (15–49 years old) in Ethiopia: Bayesian approach
date: 2021-12-22
journal: Womens Health (Lond)
DOI: 10.1177/17455065211067638
sha: e22f8aee2346639256fb2748f7900e9b831c2096
doc_id: 25136
cord_uid: pryzftjs

BACKGROUND: Human immunodeficiency virus remains the leading cause of morbidity and mortality throughout the world. Sub-Saharan Africa regions are the most affected regions and accounted for 67% of HIV infections worldwide, and 72% of the world’s AIDS-related deaths. OBJECTIVE: To estimate the prevalence of HIV and identify factors associated with it among women of reproductive age in Ethiopia. METHODS: This study was conducted based on the 2016 Ethiopian Demographic and Health Surveys data. The data were weighted using sampling weight for probability sampling and non-response to restore the representativeness of the data and get valid statistical estimates. Then, a total of 14,161 weighted sample women were used to investigate the study. Finally, a multilevel analysis was done based on the Bayesian approach to identify factors associated with HIV among women of reproductive age in Ethiopia. RESULTS: This study showed the prevalence of HIV among reproductive age group women was 0.85%. Being rural resident (adjusted odds ratio = 0.20; 95% CrI = 0.1–0.4), secondary education level (adjusted odds ratio = 0.20; 95% CrI = 0.1–0.4), rich wealth status (adjusted odds ratio = 4; 95% CrI = 3–6), married women but living separately (adjusted odds ratio = 2.3; 95% CrI = 1.2–4.5), long distance from the health facility (adjusted odds ratio = 0.4; 95% CrI = 0.3–0.5), and exposure to media (adjusted odds ratio = 2.9; 95% CrI = 1.8–4.7) were significantly associated with HIV. CONCLUSION: Being rural residents, women whose marital status is separated, wealthy, travel a long distance to get health facility, and are exposed to media are risky to be infected by HIV. Whereas being a rural resident and educated are preventive factors for HIV. Therefore, the government of Ethiopia and the ministry of health should consider those factors when they design HIV prevention and control strategies.

The human immunodeficiency virus (HIV) is a type of virus that affects the immune system of the human cell and results in a decline in CD4 + cell count and immune function, which can lead to life-threatening opportunistic infec-of which, 36.2 million adults and 1.8 million were children of age 0-14 years. 4 And also in this region, women of childbearing age comprise 61% of people living with HIV, accounting for over 12 million. 4 Sub-Saharan Africa is home to only 12% of the global population, however, accounts for 71% of the global burden of HIV infection. 5 Of the estimated 6000 new infections that occur globally each day, two out of three are in sub-Saharan Africa with young women continuing to bear a disproportionate burden. 5 The HIV epidemic has been showing remarkable variations across population subgroups, regions, and countries, at the subnational level between provinces. 6 According to the Ethiopian HIV/AIDS Prevention and Control Office (HAPCO) estimates, 0.9% of the adult population was HIV positive. Nevertheless, this distribution was highly heterogeneous, with the lowest prevalence (0.1%) in Somali and the highest in the Gambela Region (4.8%). 7 HIV/AIDS is a leading cause of death and disease burden, especially in sub-Saharan Africa. 2, 3 The HIV epidemic has had an adverse effect on efforts to reduce maternal mortality rates. 8 The government of Ethiopia has committed to ending AIDS as a public health threat by 2030. Therefore, to prevent and control the pandemic of HIV, it is very important to know factors associated with HIV precisely using an appropriate method of statistical analysis. Traditional multiple regression techniques treat the units of analysis as independent observations. 9 Therefore, standard errors for the coefficients of higherlevel predictor variables are highly affected as a result of ignoring grouping. Therefore, standard errors of regression coefficients will be underestimated and leading to an overstatement of statistical significance. Finally, the result is misleading. 9 But in Ethiopia, there are no studies conducted at the national level that shows predictors associated with HIV by considering the clustering effect using the Bayesian analysis approach among women of reproductive age. The Bayesian analysis approach is one of the data analysis approaches independent of the frequentist analysis approach. In the Bayesian analysis approach, the estimates of the parameters are estimated from the posterior distribution which is the combination of the prior information and the likelihood of the data. 10 Therefore, this study aimed to estimate the prevalence of HIV and to identify its associated factors among women of reproductive age in Ethiopia based on the 2016 Ethiopian Demographic and Health Survey (EDHS) data.

Secondary data analysis was performed based on the 2016 EDHSs data. The EDHSs are nationally representative cross-sectional surveys performed in nine regions and two city administrations in the country every 5 years. In each of the surveys, stratified two-stage sampling of clusters was carried out. Stratification was achieved by separating each region into urban and rural areas. For 2016 EDHS data, a total of 645 enumeration areas (EAs) were randomly selected proportional to the EA size in the first stage. At the second stage, on average 27-32 households per EA were selected. [11] [12] [13] The data were accessed from the Measure Demographic and Health Surveys (DHS) website (http://www.dhsprogram.com) after permission was granted through an online request by explaining the objective of this study. All reproductive age group women in the selected EAs and households during the surveys were included in this study. After excluding mothers whose test did not provide a clear negative or positive result, a total of 14,161 weighted sample women were used for this study. The comprehensive procedure for sampling and sample size determination technique, in general, about the survey, is described in the complete EDHS report. [11] [12] [13] Variables Interviewers collected finger-prick blood specimens from women age 15-49 years who consented to HIV testing. 14 The protocol for blood specimen collection and analysis was based on the anonymously linked protocol developed for the DHS Program. 14 Blood samples were dried overnight and packaged for storage the following morning. Samples were periodically collected from the field and transported to the laboratory at the Ethiopian Public Health Institute (EPHI) in Addis Ababa. 14 Upon arrival at EPHI, each blood sample was tested to know their HIV status. 14 The response variable is dichotomous with possible values Y i = Yes (coded as "Yes = 1") if the ith woman is positive and Y i = No (coded as "No = 0"), if the woman is negative for HIV. The independent variables were classified as community and individual-level variables. Place of residence (recoded as rural and urban) and region of the study participants were considered as community-level variables. Whereas marital status (categorized as married, never married, separated, divorced, and widowed), educational level (recoded as no formal education, primary, secondary, and higher), occupation (coded as had and had no work), sex of household head, desire for children (coded as yes and no), heard about HIV (coded as yes and no) media exposure, and wealth status were considered as individuallevel variables. Women who were either reading a newspaper, use the Internet, listening to the radio, or watched television at least once a week were considered as exposed to media, otherwise, they are considered as not exposed to media. 14 In 2016 EDHSs data, wealth status was measured using different assets of the households. In the data set, the categories for the wealth index were presented as poorest, poorer, middle, richer, and richest. But, in this study, the categories are recategorized as "poor," "middle," and" "rich." Poorest and poorer are merged and named as "poor" and richest and richer are merged and named as "rich." But the category middle is taken as it is.

Before any statistical analysis, the data were weighted using sampling weight for probability sampling and nonresponse to restore the representativeness of the survey and get reliable statistical estimates.

Multilevel analysis. To see the association between HIV status and independent variables, we applied multilevel analysis using Bayesian statistics. In this study, the hierarchy of data was two levels. Women in the household were considered as level one units and EAs were considered as level two units. Woman in the household was nested within EAs. The dependent variable was represented by

the category is dichotomous. Therefore, Generalized Linear Mixed Model (GLMM) was fitted based on the Bayesian approach. To estimate the parameters of the variable and the extent of random variations between clusters, we used the Brms-R package. It uses No-U-Turn Sampler (NUTS), an extension to HMC uses a recursive algorithm to build a set of likely candidate points that spans a wide swath of the target distribution, stopping automatically when it starts to double back and retrace its steps. 15 These features allow it to converge to high-dimensional target distributions much more quickly than simpler methods such as the random walk Metropolis or Gibbs sampler. 15 Therefore, this package helped us to fit multilevel models based on the Bayesian approach. 16 We used flat prior with beta distribution (1, 1) to determine the regression coefficients and gamma distribution (0.001, 0.001) to estimate the variance. And also we used iteration = 10,000, number of iterations that was discarded = 1000, number of cores = 2, chains = 2 adapt delta (controls divergent transition) = 0.95, and initials (the starting values of the iterations) = 0.

Four multilevel binary logistic regression models were fitted using compound symmetry covariance structure. The first model (a model without covariate) was fitted to determine the extent of cluster variation in HIV. The second model was fitted with individual-level factors. The third model was fitted with community variables. The fourth model was fitted with both individual-and communitylevel factors. Accordingly, all models include a random intercept (random intercept only) at the cluster level to capture heterogeneity among clusters. Then, the models were compared based on their Widely Applicable Information Criteria (WAIC) and the model with a small WAIC value was selected. Finally, all interpretations and inferences were made based on this model. To consider the variation of HIV between the EAs, we used the intraclass correlation coefficient (ICC) value greater than 5%. We used the 95% posterior credible interval test of significance and the interval that does not contain 1 is considered statistically significant. Since the results obtained from the posterior distribution are not reliable until the chain has reached its stationary distribution, 17 we used the most popular and straightforward convergence assessment methods. The methods we used to monitor the convergence of the algorithm were chains of the time series plots mixed well, smooth density plot, Rhat = 1, and Bulk_ESS and Tail_ ESS greater than 1000.

The EDHS was conducted based on the permission of the government, and informed consent was taken and participants' confidentiality was assured during that time. For this study, we accessed the data set based upon request (www.dhsprogram.com online), and there was no ethical approval required. Moreover, there are no names of individuals or household addresses in the data files.

In this study, a total of 14,159 women of reproductive age were included from all parts of the country. From this, 76.2% of the respondents were rural residents. The majority of the study participants were from Oromia regional state followed by Amhara regional state with 22%. Sixtyfive percent of the women were married and 26% of them were never married. More than three-fourth of the households were led by males ( Table 1 ).

As shown in Table 2 , the prevalence of HIV among women of reproductive age was 0.85%. Among urban residents, 2.4% were positive for HIV. And also 0.4% of rural residents were positive for HIV which is lower than urban residents. Based on region, in Gambela regional state, the highest prevalence was observed (3.9%) followed by Addis Ababa city administration (3.6%).

The variance between clusters (EA) being positive for HIV among reproductive age group women is 2.4. Therefore, we reject the null hypothesis which states that the variation across the cluster was zero. So that the variation of to be HIV positive among women of reproductive age in Ethiopia between clusters was none zero. The results also displayed that the ICC was 39%, meaning that roughly 39% of the variability in the prevalence of HIV among women of reproductive age was attributable to the clusters.

Among the four fitted models, the model fitted with individual-and community-level factors has the smallest WAIC (2147) Therefore, this model is the best-fitted model, and every interpretation and report were made based on this model.

Among the factors included in the fitted model (a model fitted by both individual-and community-level factors), residence, educational level, wealth status, occupation, marital status, distance to the health facility, and exposure to media was significantly associated with HIV in Ethiopia ( Table 2 ). The odds of being positive for HIV among rural resident women was 80% (adjusted odds ratio (AOR) = 0.2; 95% CrI = 0.1-0.4) less likely as compared to urban resident women. Based on marital status, women who were separated are 2.3 times (AOR = 2.3; 95% CrI = 1.2-4.5) more likely to be positive for HIV as compared to married women. Similarly, women who were rich in wealth were four times (AOR = 4; 95% CrI = 3-6) more likely to be infected by HIV as compared to poor women (Table 3 ).

The prevalence of HIV among reproductive age group women in Ethiopia was 0.85%. This value is lower as compared to the pooled estimates of the Sub-Saharan Africa regions. 18 This variation might be due to the introduction of the Health Extension Program (HEP), implementation of different programs, including education and stigmareduction programs, behavior change initiatives, improving access to health care, and expansion of HIV testing for all age groups across the territory. 19 The country also has committed to reducing new HIV infections by 50% by 2020 and to ending AIDS as a public health threat by 2030. This is reflected in the Country's Health Sector Transformation Plan II 2015-2020, where one of the major indicators is a reduction of HIV incidence rate. 20 The multilevel binary logistic regression analysis revealed that place of residence, educational level, occupation, wealth status, marital status, distance to the health facility, and exposure to media have an association with HIV among reproductive age group women in Ethiopia. The risk of being positive for HIV among rural women is less likely as compared to women who were urban residents. This finding is supported by the study conducted in Malawi. 21 This could be, in rural areas, there might be no commercial sex workers who are at risk for HIV. Similarly, the odds of being positive among women whose educational level was secondary were higher as compared to women who had no formal education. The odds of having HIV among women who had no work were more likely than women who had work. This result is supported by a study conducted in Sub-Saharan Africa. 22 This might be due to the richest women tend to be more educated, more mobile, and more likely to live in urban areas, where HIV is more prevalent. The odds of having HIV among separated women were 2.3 times more likely as compared to women who were married. This finding is consistent with a study conducted in South Africa. 23 This could be due to the fact that separated women are exposed to low social cohesion, and the resulting destabilization of sexual relationships may increase the risk of HIV. Women who had travel long distances to reach health facilities were less likely to have HIV as compared to women who live near health facilities. This could be due to the fact that women who lived far from health facilities are mostly rural residents, are less likely to attend nightclubs, and have no commercial sex workers around their living environment, which are the commonest predisposing event to gain HIV.

The odds of having HIV among women who had not been exposed to media were three times more likely than those women who were exposed to media. This study is consistent with a study conducted in Bangladesh, 24 India, 25 and Iran. 26 This might be due to women who had not been exposed to media did not augment their HIV-related knowledge. In addition, they did not enhance their knowledge concerning HIV prevention modalities.

To get fine estimates of the parameters, fitting a multilevel model by considering all the national regional states of Ethiopia was the strength of the study. As a limitation, we 

From the factors included in the selected model, place of residence, educational level, marital status, wealth status, occupation, distance to the health facility, and exposure to media have an association with HIV. Therefore, Ethiopian government, the ministry of health, and any other responsible body consider these factors when they design HIV prevention and control strategies.

participation, voluntary participation was ensured during the interviews. 14 The detail about ethical approval and consent to participate is available in the EDHS's report. 14

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

Yilkal Negesse https://orcid.org/0000-0003-4751-3615

Gossa Manikelkil https://orcid.org/0000-0002-7834-8838

Supplemental material for this article is available online.

Radiation oncology in Africa: improving access to cancer care on the African continent

Global HIV/AIDS response: epidemic update and health sector progress towards universal access: progress report

Estimates of global, regional, and national incidence, prevalence, and mortality of HIV, 1980-2015: the Global Burden of Disease Study

Maintaining HIV care during the COVID-19 pandemic

HIV prevalence and trends in sub-Saharan Africa: no decline and large subregional differences

The evolving epidemiology of HIV/AIDS

HIV/AIDS in Ethiopia: data from the 2011 Ethiopia demographic and health survey

Maternal health and HIV

Multilevel analysis of matching behavior

Multilevel analysis of determinants of anemia prevalence among children aged 6-59 months in Ethiopia: classical and Bayesian approaches

Addis Ababa: Central Statistical Authority and ORC Macro, 2001. 12. CSA and ICF. Ethiopia Demographic and Health Survey

Central Statistical Authority and ORC Macro. Ethiopia Demographic and Health Survey

Demographic and Health Survey 2016: key indicators report

The No-U-Turn sampler: adaptively setting path lengths in Hamiltonian Monte Carlo

Advanced Bayesian multilevel modeling with the R package brms

Eliciting expert knowledge for Bayesian logistic regression in species habitat modelling

Recent HIV prevalence trends among pregnant women and all women in sub-Saharan Africa: implications for HIV estimates

Revised classification system for HIV infection and expanded surveillance case definition for AIDS among adolescents and adults

Spatial analysis of factors associated with HIV infection in Malawi: indicators for effective prevention

HIV infection does not disproportionately affect the poorer in sub-Saharan Africa

Does marital status matter in an HIV hyperendemic country? Findings from the 2012 South African National HIV Prevalence, Incidence and Behaviour Survey

Educating women for HIV prevention: does exposure to mass media make them more knowledgeable

Access to media and HIV knowledge in India

Effects of media, interpersonal communication and religious attitudes on HIV-related stigma in Tehran, Iran

The authors acknowledge the Demographic and Health Surveys (DHS) center allowing us to access the data set.

Y.N., G.M., M.S., and G.F. were participated in this paper from the inception to design, acquisition of data, data cleaning and coding, data analysis, and interpretation, and drafting and finally revising of the manuscript. Then, the final draft of the manuscript was prepared by Y.N. Finally, all authors read and approved the revised manuscript.

The data were accessed from the Measure DHS website (http:// www.dhsprogram.com) after permission was granted through an online request by explaining the objective of this study. The data sets analyzed during this study are available from the corresponding author upon reasonable request.

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

Ethical clearance for the EDHS data was given by the Ethiopian Health and Nutrition Research Institute (EHNRI) Review Board; the National Research Ethics Review Committee (NRERC) at the Ministry of Science and Technology. 14 After the respondents were informed about the survey and consent was taken for their