- Research article
- Open Access
- Open Peer Review
Influence of socioeconomic position and gender on current cigarette smoking among people living with HIV in sub-Saharan Africa: disentangling context from composition
BMC Public Healthvolume 16, Article number: 998 (2016)
Smoking is still gaining ground in Sub-Saharan Africa, especially among socially disadvantaged groups. People living with HIV represent a subgroup with a significantly elevated prevalence of cigarette smoking. The objective of the study was to examine the influence of individual-, neighbourhood- and country-level socioeconomic position on current cigarette smoking among people living with HIV in Sub-Saharan Africa.
We applied multivariable multilevel logistic regression analysis on Demographic and Health Survey data collected between 2003 and 2012 in sub-Saharan Africa. We identified 31,270 individual living with HIV (Level 1) nested within 7,054 neighbourhoods (Level 2) from 19 countries (Level 3).
After adjustment for individual-, neighbourhood- and country-level factors, respondents, the following significant independent risk factors for increasing odds of being a current cigarette smokers among people living with HIV: male gender (odds ratio [OR] = 62.49; 95 % credible interval [CrI] 45.93 to 78.28), from the poorer households (OR = 1.62, 95 % CrI 1.38 to 1.90); living in urban areas (OR = 1.24, 95 % CrI 1.09 to 1.41), from neighbourhoods with low poverty rate (OR = 1.25, 95 % CrI 1.09 to 1.43), illiteracy rate (OR = 1.28, 95 % CrI 1.14 to 1.42), low unemployment rate (OR = 1.11, 95 % crI 1.01 to 1.43); and from countries with low socio-economic deprivation (OR = 1.53, 95 CrI 1.08 to 1.96). About 3.4 % and 39.4 % variation in cigarette smoking behaviour among people living with HIV is conditioned by differences between neighbourhoods and countries.
Gender, education and socioeconomic context are independently associated with current cigarette smoking among people living with HIV in sub-Saharan Africa.
Tobacco use and HIV infection are two major causes of death globally that continue to grow . The intersection of these two epidemics represents an area of growing clinical and public health importance, especially in sub-Saharan Africa. Living with HIV is associated with a two-fold increase in the likelihood of smoking cigarettes [2–4]. Individual living with HIV are more prone to the adverse health effects of tobacco use than those without [5–7]. In addition to the increased risk of numerous AIDS- and non-AIDS related disease, cigarette smoking has an adverse impact upon the health-related quality of life of people living with HIV . Individual living with HIV in care loose more years of life to smoking than to HIV Infection .
Without objective and sound information about factors associated with cigarette smoking behaviours among people living with HIV, it will be difficult to plan effective care and strategies for smoking cessation. In the context of sub-Saharan Africa, much research has focused on individual compositional factors associated with cigarette smoking behaviours for both individuals living with HIV [10–13] and for general population [14–17]. It is well established that individual-level characteristics such as age, educational attainment, occupation and income are associated with smoking behaviours [10–17]. Indeed, we found no published studies that had examined contextual factors associated with current cigarette smoking among people living with HIV. This neglect is important given the central role of neighbourhoods in forming smoking habits [18–23], as they shape individual opportunities and expose residents to multiple risks and resources over the life course [24, 25]. The purpose of this study was to develop and test a model of factors associated with current cigarette smoking among people living with HIV that includes individual-level compositional socioeconomic characteristics along with contextual socioeconomic characteristics defined at the neighbourhood- and country level. We further, examined how much of the variation in cigarette smoking behaviour among people living with HIV is conditioned by differences between neighbourhoods and countries.
Study design and data
Data for this cross-sectional study were obtained from Demographic and Health Surveys (DHS), which are nationally representative household surveys conducted in low- and middle-income countries. This study used data from 19 recent DHS surveys conducted between 2003 and 2012 in sub-Saharan Africa available as of October 2014 and that included rapid HIV test results and questions on self-reported tobacco use. The DHS uses a multi-stage, stratified sampling design with households as the sampling unit . Within each sample household, all women and men meeting the eligibility criteria are interviewed. Because the surveys are not self-weighting, weights are calculated to account for unequal selection probabilities as well as for non-response. With weights applied, survey findings represent the full target populations. The DHS surveys include a household questionnaire, a women’s questionnaire, and in most countries, a men’s questionnaire. All three DHS questionnaires are implemented across countries with similar interviewer training, supervision, and implementation protocols.
For the HIV testing, blood spots were collected on filter paper from a finger prick and transported to a laboratory for testing. The laboratory protocol includes an initial enzyme-linked immunosorbent assay (ELISA) test, and then retesting of all positive tests and 5–10 % of the negative tests with a second ELISA. For those with discordant results on the two ELISA tests, a new ELISA or a Western Blot is performed . Participation in the test was voluntary and before collecting blood samples each selected participant was asked to provide informed consent to the testing . In order to ensure confidentiality, the HIV test results were anonymously linked to individual questionnaire information .
Respondents were explicitly asked “Do you currently smoke cigarettes?” Those who responded ‘yes’ to this question were defined as current cigarette smokers, whereas those who responded ‘no’ were defined as current non-smokers.
Individual level factors
The following individual-level factors were included in the models: sex of the respondent (male versus female), respondents’ age in completed years (18 to 24, 25 to 34, 35 to 44 or 45 or older), educational attainment (no education, primary, secondary or higher); marital status (never married versus ever married) and occupation (working or not working). DHS did not collect direct information on household income and expenditure. We used DHS wealth index as a proxy indicator for socioeconomic position. The methods used in calculating DHS wealth index have been described elsewhere [28, 29]. Briefly, an index of economic status for each household was constructed using principal components analysis based on the following household variables: number of rooms per house, ownership of car, motorcycle, bicycle, fridge, television and telephone as well as any kind of heating device. From these criteria the DHS wealth index tertiles (poor, middle, and rich) were calculated and used in the subsequent modelling.
We used the term neighbourhood to describe clustering within the same geographical living environment. Neighbourhoods were based on sharing a common primary sample unit within the DHS data. The sampling frame for identifying primary sample unit in the DHS is usually the most recent census. The unit of analysis was chosen for two reasons. First, primary sample unit is the most consistent measure of neighbourhood across all the surveys , and thus the most appropriate identifier of neighbourhood for this cross-region comparison. Second, for most of the DHS conducted, the sample size per cluster meet the optimum size with a tolerable precision loss .
The following neighbourhood-level factors were included in the models: place of residence (rural or urban area), neighbourhood poverty-, illiteracy- and unemployment rates. We categorized neighbourhood poverty-, illiteracy- and unemployment rates into two categories (low and high), to allow for non-linear effects and provide results that were more readily interpretable in the policy arena. Median values served as the reference group for comparison.
Country-level data were collected from the reports published by the United Nations Development Program . At country-level, we included percentage rural population and intensity of deprivation. Intensity of deprivation is average percentage of deprivation experienced by people in multidimensional poverty. Like wealth index, intensity of deprivation was computed using principal component based on data on household deprivations in education, health and living standards, however, at the country-level . The country-level variables were also categorized into two (low and high) levels.
In the descriptive statistics the distribution of respondents by key variables were expressed as percentages.
We used multivariable logistic multilevel regression models to analyse the association between individual compositional and contextual factors associated with current cigarette smoking among people living with HIV. We specified a 3-level model for binary response reporting current cigarette smoking or not-currently smoking, for people living HIV (at level 1), in a neighbourhood (at level 2) living in a country (at level 3) (see Fig. 1).
We constructed five models. The first model, an empty or unconditional model without any explanatory variables, was specified to decompose the amount of variance that existed between country and neighbourhood levels. The second model contained only individual-level factors, the third model contained only neighbourhood-level factors, and fourth model contained only country-level factors. Finally, the fifth model simultaneously controlled for individual-, neighbourhood- and country-level factors (Full Model).
Fixed effects (measures of association)
The results of fixed effects (measures of association) were reported as odds ratios (ORs) with their 95 % credible intervals (CrIs). Bayesian statistical inference provides probability distributions for measures of association (ORs), which can be summarized with 95 % credible intervals (95 % CrI), rather than 95 % confidence intervals (95 % CI). A 95 % credible interval can be interpreted as there being a 95 % probability that the parameter takes a value in the specified range.
Random effects (measures of variation)
The possible contextual effects were measured by the intraclass correlation (ICC) and median odds ratio (MOR). We measured similarity between respondents in the same neighbourhood and within the same country using ICC. The ICC represents the percentage of the total variance in the probability of reporting current cigarette smoking that is related to the neighbourhood- and country-level, i.e. measure of clustering of odds of reporting cigarette smoking in the same neighbourhood and country. The ICC was calculated by the linear threshold (latent variable method) . Following the ideas of Larsen et. al. on neighbourhood effects , we reported the random effects in terms of odds. The MOR measures the second or third level (neighbourhood or country) variance as odds ratio and estimates the probability of being a current cigarette smoker that can be attributed to neighbourhood and country context. MOR equal to one indicates no neighbourhood or country variance. Conversely, the higher the MOR, the more important are the contextual effects for understanding the probability of being a current cigarette smoker.
Model fit and specifications
We checked for multi-collinearity among explanatory variables examining the variance inflation factor (VIF) , all diagonal elements in the variance-covariance (τ) matrix for correlation between −1 and 1, and diagonal elements for any elements close to zero. None of the results of the tests provided reasons for concern. Thus, the models provide robust and valid results. The MLwinN software, version 2.31, was used for the analyses [36, 37]. Parameters were estimated using the Markov Chain Monte Carlo procedure . The Bayesian Deviance Information Criterion was used as a measure of how well our different models fitted the data. A lower value on Deviance Information Criterion indicates a better fit of the model .
The countries, year of data collection, and the surveys characteristics are listed in Table 1. The surveys were conducted between 2003 and 2012. The median number of neighbourhoods sampled was 569, ranging from 76 in Sao Tome and Principle to 9008 in Cote D’Ivoire. The median number of women and men living with HIV in the analysis was 370 and 199, respectively. The median prevalence of current cigarette smoking among women living with HIV was 0.6 %. The prevalence of current cigarette smoking among men living with HIV ranged from 10.3 % in Sao Tome and Principle to as much as 41.6 % in Lesotho. Table 2 presents the descriptive statistics for the final pooled sample. For this analysis, we analysed information on 31,270 people living with HIV (Level 1) nested within 7,052 neighbourhoods (Level 2) from 19 countries (Level 3) in sub-Saharan Africa. More than half of the respondents were female (58 %) and the majority of women were between the 18 and 34 years of age (60 %). Thirty-six percent of the respondents had no formal education. Most of the respondents were currently employed (70 %) and never married (81 %). Most of the respondents were living in rural.(60 %), high poverty rate. (68 %), high illiteracy rate. (61 %), and high unemployment rate. (55 %) neighbourhoods.
Measures of associations (fixed effects)
The results of different models are shown in Table 3. In the fully adjusted model controlling for the effects of individual-, neighbourhood- and societal-level factors, men living with HIV were more likely to be current cigarette smokers than women living with HIV (OR = 62.49, 95 % CrI 45.93 to 78.28). Respondents with primary education were significantly more likely to be current cigarette smokers than those with secondary or higher education (OR = 1.38; 95 % CrI 1.24 to 1.53). Respondents from the poorer households were 62 % more likely to be current cigarette smokers as those from the richer households (OR = 1.62, 95 % CrI 1.38 to 1.90). Respondents currently working were significantly more likely to be current smokers than those not working (OR = 1.27, 95 % CrI 1.11 to 1.46). Respondents from urban areas were more likely to be current smokers than those from rural areas (OR 1.24; 95 % CrI 1.09 to 1.43). Respondents from low poverty (OR = 1.25, 95 % CrI 1.09 to 1.43) and illiteracy (OR = 1.28, 95 % CrI 1.14 to 1.42) rates neighbourhoods were more likely to be current cigarette smokers than those from high rates neighbourhoods. Finally, respondents from low intensity of deprivation were significantly more likely to be current cigarette smokers than those from high intensity of deprivation (OR = 1.53, 95 CrI 1.08 to 1.96).
Measures of variations (random effects)
As shown in Table 3, in Model 1 (unconditional model), there was a significant variation odds of reporting current cigarette smoking across the countries (σ2 = 0.18, 95 % CrI 0.07 to 0.40) and across the neighbourhoods (σ2 = 1.96, 95 % CrI 1.70 to 2.23). According to the intra-country and intra-neighbourhood correlation coefficient, 3.4 % and 39.4 % of the variance in odds of reporting cigarette smoking could be attributed to the country- and neighbourhood-level factors, respectively. Results from the median odds ratio (MOR) also confirmed evidence of neighbourhood and societal contextual phenomena shaping individual cigarette smoking behaviour. From the full model (Model 5), it was estimated that if a respondent moved to another neighbourhood or another country with a higher probability of cigarette smoking, the median increase in their odds of being current cigarette smoker would be 1.14 (95 % CrI 1.08 to 1.27) and 1.54-fold (95 % CrI 1.31 to 1.90).
To our knowledge, the current study is the first multilevel examination of smoking behaviour among people living with HIV in sub-Saharan Africa using national representative data and a very large number of respondents (31,270). We found that at individual-level, respondents with primary education (vs those with secondary or higher), those currently working, and poorer households were significantly more likely to be current cigarette smokers. Gender appeared to be a very important risk factor, men living with HIV were significantly more likely to be current cigarette smokers than women living with HIV. The findings corroborate those of previous studies that examined the association between cigarette smoking and individual socioeconomic position both for people living with HIV [10–13] and the general population [14–17] in sub-Saharan Africa. More importantly, the findings uncover new evidence by demonstrating that neighbourhood- and country-level factors influence cigarette smoking behaviours above and beyond individual level factors.
Furthermore, our findings reveal a striking example of Simpson’s paradox [39, 40], a situation where the association between two variables (e.g. cigarette smoking and area poverty rate) are reversed when a third variable (e.g. individual wealth index) is considered. The results of including only neighbourhood-level factors showed that living in neighbourhood with low rates of socioeconomic position decreases the odds of being a current cigarette smoker. However, in the full model, when individual-, neighbourhood- and country-level factors were adjusted for simultaneously, there is a evidence that living in neighbourhood with low rates of socioeconomic position increases the odds of being a current cigarette smoker. The findings corroborate those of previous studies from high-income countries that had examined the effect of contextual effects on smoking in the general population [18–23].
We found evidence of geographical clustering in current cigarette smoking behaviours. About 3.4 % and 39.4 % of the variation in cigarette smoking behaviour among people living with HIV, is conditioned by differences between neighbourhoods and countries, respectively. If a respondent moved to another neighbourhood or another country with a higher probability of cigarette smoking, their odds of becoming a cigarette smoker may increase by about 14 % and 54 %, respectively. It is instinctual that people living with HIV from the same neighbourhood may be more similar to each other in relation to their current cigarette smoking behaviours than to others from other neighbourhoods , i.e. contextual phenomenon expresses itself as clustering of individual current cigarette smoking behaviours within neighbourhood. On these grounds, we might conclude that there is some evidence for a possible neighbourhood and country contextual phenomenon shaping a common individual cigarette smoking behaviours. These findings underscore the need to implement public health prevention strategies not only at the high-risk individual level, but also high-risk neighbourhoods, such as urban slum areas. In addition, there is a need for longitudinal studies to explain how deleterious behaviours are transmitted among individual residents, i.e. mechanisms that connect the people and the observed contextual factors. Similarly, further decomposition analyses could provide further evidence about important factors that could explain the disparities in cigarette smoking behaviour among high-risk individuals and high-risk places.
Our findings should be considered in light of the following limitations. Firstly, the cross-sectional nature of the data limits our ability to draw any causal inferences on the reported association reported. Secondly, we did not measure the length of time that participants had spent in their neighbourhoods and the extent of their exposure to the neighbourhood environment. We were, thus, unable to determine whether the associations of neighbourhood characteristics with current cigarette smoking were due to cumulated effects. Finally, one important limitation is that DHS surveys do not collect data on household income or expenditure, the traditional indicators used to measure wealth. The assets-based wealth index used here is only a proxy indicator of household economic status, and it does not always produce results similar to those obtained from direct measurements of income and expenditure where such data are available or can be collected reliably [28, 42]. In addition, we recognize potential data limitations that should be borne in mind while interpreting our findings. The observed magnitude of the association could have been under-estimated, since some African countries with huge populations and high HIV burden such as Nigeria could not be included in the model because of non-availability of comparable data. These countries were missing at not random, it is not clear why the DHS program did not include testing for HIV Infection in some countries.
Despite these limitations, the study strengths are significant. It is a large, population-based study with national coverage from 19 countries with high response rates. The DHS has some important advantages when compared with other surveys. They are often nationally representative, allowing for conclusions that cover the entire nation. In addition, variables in DHS were operationalized in the same way and making it possible for numerical values comparable across countries. There are advantages to studying factors associated with current cigarette smoking using a multilevel approach; we are able to provide more robust evidence about individual compositional and contextual measures of socioeconomic position associated with current cigarette smoking. The Bayesian approach we adopted has the additional advantage of being able to produce a far more robust estimate with better properties and yields unbiased estimates [43, 44]. Bryan and Jenkins state in their excellent discussion of this issue recommended : “a third strategy would be to move beyond the classical (‘frequentist’) statistical framework used by most applied social science researchers and to make greater use of Bayesian methods of estimation and inference, as there is some evidence that they perform better in the mall number of countries.”
In conclusion, individual compositional and contextual measures of socioeconomic position were independently associated with current cigarette smoking among people living with HIV in sub-Saharan Africa, which underscores the need to implement cigarette smoking prevention strategies not only at the individual level taking into account socioeconomic position, but also at the contextual levels.
Furber AS, Maheswaran R, Newell JN, Carroll C. Is smoking tobacco an independent risk factor for HIV infection and progression to AIDS? A systemic review. Sex Transm Infect. 2007;83(1):41–6.
Burkhalter JE, Springer CM, Chhabra R, Ostroff JS, Rapkin BD. Tobacco use and readiness to quit smoking in low-income HIV-infected persons. Nicotine & tobacco research : official journal of the Society for Research on Nicotine and Tobacco. 2005;7(4):511–22.
Gritz ER, Vidrine DJ, Lazev AB, Amick 3rd BC, Arduino RC. Smoking behavior in a low-income multiethnic HIV/AIDS population. Nicotine & tobacco research : official journal of the Society for Research on Nicotine and Tobacco. 2004;6(1):71–7.
Mamary EM, Bahrs D, Martinez S. Cigarette smoking and the desire to quit among individuals living with HIV. AIDS Patient Care STDs. 2002;16(1):39–42.
Niaura R, Shadel WG, Morrow K, Tashima K, Flanigan T, Abrams DB. Human immunodeficiency virus infection, AIDS, and smoking cessation: the time is now. Clinical infectious diseases : an official publication of the Infectious Diseases Society of America. 2000;31(3):808–12.
Calvo M, Laguno M, Martinez M, Martinez E. Effects of Tobacco Smoking on HIV-Infected Individuals. AIDS Rev. 2015;17(1):47–55.
Duval X, Baron G, Garelik D, Villes V, Dupre T, Leport C, Lert F, Peretti-Watel P, Ravaud P, Spire B. Living with HIV, antiretroviral treatment experience and tobacco smoking: results from a multisite cross-sectional study. Antivir Ther. 2008;13(3):389–97.
Turner J, Page-Shafer K, Chin DP, Osmond D, Mossar M, Markstein L, Huitsing J, Barnes S, Clemente V, Chesney M. Adverse impact of cigarette smoking on dimensions of health-related quality of life in persons with HIV infection. AIDS Patient Care STDs. 2001;15(12):615–24.
Helleberg M, Afzal S, Kronborg G, Larsen CS, Pedersen G, Pedersen C, Gerstoft J, Nordestgaard BG, Obel N. Mortality attributable to smoking among HIV-1-infected individuals: a nationwide, population-based cohort study. Clinical infectious diseases : an official publication of the Infectious Diseases Society of America. 2013;56(5):727–34.
Iliyasu Z, Gajida AU, Abubakar IS, Shittu O, Babashani M, Aliyu MH. Patterns and predictors of cigarette smoking among HIV-infected patients in northern Nigeria. Int J STD AIDS. 2012;23(12):849–52.
Jaquet A, Ekouevi DK, Aboubakrine M, Bashi J, Messou E, Maiga M, Traore HA, Zannou M, Guehi C, Ba-Gomis FO, et al. Tobacco use and its determinants in HIV-infected patients on antiretroviral therapy in West African countries. The international journal of tuberculosis and lung disease : the official journal of the International Union against Tuberculosis and Lung Disease. 2009;13(11):1433–9.
Louwagie GM, Ayo-Yusuf OA. Tobacco use patterns in tuberculosis patients with high rates of human immunodeficiency virus co-infection in South Africa. BMC Public Health. 2013;13:1031.
Waweru P, Anderson R, Steel H, Venter WD, Murdoch D, Feldman C. The prevalence of smoking and the knowledge of smoking hazards and smoking cessation strategies among HIV- positive patients in Johannesburg, South Africa. South African medical journal = Suid-Afrikaanse tydskrif vir geneeskunde. 2013;103(11):858–60.
Caleyachetty R, Tait CA, Kengne AP, Corvalan C, Uauy R, Echouffo-Tcheugui JB. Tobacco use in pregnant women: analysis of data from Demographic and Health Surveys from 54 low-income and middle-income countries. The Lancet Global health. 2014;2(9):e513–20.
Hosseinpoor AR, Parker LA, Tursan d’Espaignet E, Chatterji S. Social determinants of smoking in low- and middle-income countries: results from the World Health Survey. PLoS One. 2011;6(5), e20331.
Hosseinpoor AR, Parker LA, Tursan d’Espaignet E, Chatterji S. Socioeconomic inequality in smoking in low-income and middle-income countries: results from the World Health Survey. PLoS One. 2012;7(8), e42843.
Pampel F. Tobacco use in sub-Sahara Africa: estimates from the demographic health surveys. Soc Sci Med. 2008;66(8):1772–83.
Chaix B, Guilbert P, Chauvin P. A multilevel analysis of tobacco use and tobacco consumption levels in France: are there any combination risk groups? Eur J Public Health. 2004;14(2):186–90.
Chuang YC, Li YS, Wu YH, Chao HJ. A multilevel analysis of neighborhood and individual effects on individual smoking and drinking in Taiwan. BMC Public Health. 2007;7(147):151.
Duncan C, Jones K, Moon G. Smoking and deprivation: are there neighbourhood effects? Soc Sci Med. 1999;48(4):497–505.
Frohlich KL, Potvin L, Gauvin L, Chabot P. Youth smoking initiation: disentangling context from composition. Health Place. 2002;8(3):155–66.
Karvonen S, Sipila P, Martikainen P, Rahkonen O, Laaksonen M. Smoking in context - a multilevel approach to smoking among females in Helsinki. BMC Public Health. 2008;8:134.
Mathur C, Erickson DJ, Stigler MH, Forster JL, Finnegan Jr JR. Individual and neighborhood socioeconomic status effects on adolescent smoking: a multilevel cohort-sequential latent growth analysis. Am J Public Health. 2013;103(3):543–8.
Leventhal T, Brooks-Gunn J. The neighborhoods they live in: the effects of neighborhood residence on child and adolescent outcomes. Psychol Bull. 2000;126(2):309–37.
Sampson RJ. The neighborhood context of well-being. Perspect Biol Med. 2003;46(3 Suppl):S53–64.
Rutstein SO, Rojas G. Guide to DHS Statistics: Demographic and Health Surveys Methodology. Calverton: ORC Macro; 2006.
Macro ICF. HIV prevalence estimates from Demographic and Health Surveys. Calverton: ICF Macro International; 2010.
Montgomery MR, Gragnolati M, Burke KA, Paredes E. Measuring living standards with proxy variables. Demography. 2000;37(2):155–74.
Vyas S, Kumaranayake L. Constructing socio-economic status indices: how to use principal components analysis. Health Policy Plan. 2006;21(6):459–68.
Griffiths P, Madise N, Whitworth A, Matthews Z. A tale of two continents: a multilevel comparison of the determinants of child nutritional status from selected African and Indian regions. Health Place. 2004;10(2):183–99.
Kravdal O. A simulation-based assessment of the bias produced when using as contextual variables in multilevel models. Demographic-Research. 2006;15(1):1–20.
United Nations Development Programme. Human Development Report 2014. Sustaining Human Progress: Reducing Vulnerabilities and Building Resilience. NY, USA: United Nations Development Programme; 2014.
Snijders T, Bosker R. multilevel analysis – an introduction to basic and advanced multilevel modelling. Thousand Oaks: SAGE publications; 1999.
Larsen K, Merlo J. Appropriate assessment of neighborhood effects on individual health: integrating random and fixed effects in multilevel logistic regression. Am J Epidemiol. 2005;161(1):81–8.
Tu YK, Kellett M, Clerehugh V, Gilthorpe MS. Problems of correlations between explanatory variables in multiple regression analyses in the dental literature. Br Dent J. 2005;199(7):457–61.
Browne WJ. MCMC Estimation in MLwiN v2.31. University of Bristol: Centre for Multilevel Modelling; 2014.
Rasbash J, Charlton C, Browne WJ, Healy M, Cameron B. MLwiN Version 2.31. University of Bristol: Centre for Multilevel Modelling; 2014.
Spiegelhalter DJ, Best NG, Carlin BP, Van Der Linde A. Bayesian measures of model complexity and fit. J R Stat Soc Ser B Stat Methodol. 2002;64(4):583–639.
Hernan MA, Clayton D, Keiding N. The Simpson’s paradox unraveled. Int J Epidemiol. 2011;40(3):780–5.
Simpson EH. The interpretation of interaction in contingency tables. J R Stat Soc Ser B. 1951;13:238–41.
Merlo J, Chaix B, Yang M, Lynch J, Rastam L. A brief conceptual tutorial of multilevel analysis in social epidemiology: linking the statistical concept of clustering to the idea of contextual phenomenon. J Epidemiol Community Health. 2005;59(6):443–9.
Filmer D, Pritchett LH. Estimating wealth effects without expenditure data--or tears: an application to educational enrollments in states of India. Demography. 2001;38(1):115–32.
Browne WJ, Draper D. A comparison of Bayesian and likelihood-based methods for fitting multilevel models. Bayesian Anal. 2006;1:473–514.
Stegmueller D. How Many Countries for Multilevel Modeling? A Comparison of Frequentist and Bayesian Approaches. Am J Polit Sci. 2013;57(3):748–61.
Bryan ML, Jenkins SP. Multilevel Modelling of Country Effects: A Cautionary Tale. Eur Sociol Rev. 2015.
The authors are grateful to DHS Program for providing them with the DHS data.
This paper presents independent research funded by the Swedish Research Council for Health, Working Life and Welfare under its FORTE Marie Curie International Post Doc to OAU (Grant Reference no: 2012–0064).
Availability of data and materials
The data supporting this article is available at: http://dhsprogram.com/data/available-datasets.cfm
OAU, AK and TM were involved in the conception of the study. OAU carried out data extraction. OAU conducted statistical analysis under supervision of AK and TM. OAU drafted the paper with contributions from the co-authors. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
This study was based on an analysis of existing survey data with all identifier information removed. The survey was approved by the Ethics Committee of the ICF Macro at Fairfax, Virginia in the USA and by the National Ethics Committees in their respective countries. All study participants gave informed consent before participation and all information was collected confidentially.