Article Text

Download PDFPDF

Coverage and factors associated with full immunisation among children aged 12–59 months in Bangladesh: insights from the nationwide cross-sectional demographic and health survey
  1. Abdur Razzaque Sarker1,2,
  2. Raisul Akram1,
  3. Nausad Ali1,
  4. Marufa Sultana1
  1. 1 International Centre for Diarrhoeal Disease Research Bangladesh, Dhaka, Bangladesh
  2. 2 Bangladesh Institute of Development Studies, Dhaka, Bangladesh
  1. Correspondence to Abdur Razzaque Sarker; razzaque{at}bids.org.bd

Abstract

Objective To estimate the coverage and factors associated with full immunisation coverage among children aged 12–59 months in Bangladesh.

Study design The study is cross sectional in design. Secondary dataset from Bangladesh Demographic and Health Survey was used for this analysis. Immunisation status was categorised as ‘fully immunised’ if the children had received all the eight recommended vaccine doses otherwise ‘partially/unimmunised’.

Settings Bangladesh.

Participant Children aged 12–59 months were the study participants. Participants were randomly selected through a two-stage stratified sampling design. A total of 6230 children were eligible for the analysis.

Results About 86% of the children (5356 out of 6230) were fully immunised. BCG has the highest coverage rate (97.1%) followed by oral polio vaccine 1 (97%) and pentavalent 1 (96.6%), where the coverage rate was the lowest for measles vaccine (88%). Coverage was higher in urban areas (88.5%) when compared with rural ones (85.1%). Full immunisation coverage was significantly higher among children who lived in the Rangpur division (adjusted OR (AOR)=3.46; 95% CI 2.45 to 4.88, p<0.001), were 48–59 months old (AOR=1.32; 95% CI 1.06 to 1.64, p=0.013), lived in a medium size family (AOR=1.56; 95% CI 1.32 to 1.86, p<0.001), had parents with a higher level of education (AOR=1.96; 95% CI 1.21 to 3.17, p=0.006 and AOR=1.55; 95% CI 1.05 to 2.29, p=0.026) and belonged to the richest families (AOR=2.2; 95% CI 1.5 to 3.21, p<0.001). The likelihood of being partially or unimmunised was higher among children who had the father as their sole healthcare decision-maker (AOR=0.69; 95% CI 0.51 to 0.92, p<0.012).

Conclusions There were significant variations of child immunisation coverage across socioeconomic and demographic factors. These findings will inform innovative approaches for immunisation programmes, and the introduction of relevant policies, including regular monitoring and evaluation of immunisation coverage—particularly for low-performing regions, so that the broader benefit of immunisation programmes can be achieved in all strata of the society.

  • health policy
  • public health
  • coverage
  • bangladesh

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Strengths and limitations of this study

  • Our study did a pioneer investigation of full immunisation coverage across various regions in Bangladesh using the latest nationwide demography and health survey dataset.

  • Multivariable logistic regression models were constructed to assess the potential factors associated with full immunisation coverage.

  • The study results can be generalised at the country level because of its large sample size and using the latest nationally representative demographic and health survey data of Bangladesh.

  • Due to the nature of the study, the possibility of recall bias with our results cannot be ignored.

Background

Universal immunisation programme of children against vaccine-preventable diseases (VPD) has been recognised as one of the most cost-effective programmes to diminish childhood mortalities and morbidities across the world.1 2 Every year, vaccination against VPDs prevents debilitating illness and disability, saving millions of young lives globally.3 Over the decades, remarkable improvements have been made towards the development of national immunisation programmes, with the Expanded Program on Immunization (EPI) being a major contribution to this success.4 EPI was formally established in 1974, with the support of WHO, with the goal of immunising every child against six VPD (ie, diphtheria, pertussis, tetanus (DPT), poliomyelitis, measles and tuberculosis) by 1990.5 The utmost priority was given to developing countries because of the higher prevalence and inadequate service delivery for immunisation within these settings.4 The government of Bangladesh had initiated the EPI with the support of Unicef and WHO through various outreach activities from 1979, with the overall objective to immunise all children by 1990 to prevent the VPDs and to eradicate poliomyelitis.5 6 Later, in 1995, national immunisation day was initiated by EPI to sustain the polio-free status and also to increase the immunisation coverage that was achieved through different operational activities.6 7 The implementation of the EPI has already been shown to be a great success globally. This can be seen through the significant improvement in child immunisation coverage and the eradication of poliomyelitis.8 In the line of success, Bangladesh has experienced impressive improvements in increasing immunisation coverage and a significant contribution to the reduction of childhood morbidity, mortality and also maintaining its polio-free status.9 10 Despite the success in reducing the child mortality rate by two-thirds since 1990, Bangladesh recorded 0.1 million child deaths in 201610 and was among the top 10 countries that had the highest childhood mortality globally.11 A recent study revealed that almost half of these deaths could be prevented through immunisation alone.3 12 Therefore, a greater focus on the completion of all recommended EPI vaccines is needed to further achieve gains in decreasing childhood morbidities and mortalities.7 Unquestionably, improving the utilisation of routine immunisation services, and easy access to all vaccines is the best option to improve the immunisation coverage.

Although a number of studies have reported findings on child immunisation, few of them have generated evidence about the sociodemographic factors associated with full vaccination among children aged 12–59 months.13–15 Furthermore, available studies have either focused on specified geographical settings, age groups or ethnic groups rather than the nationwide setting and have not identified the determinants of individual vaccines separately.14–21 This study thus sought to capture the full immunisation status using the latest country representative Demographic and Health Survey (DHS) dataset. The DHS provides reliable information on individual-level immunisation coverage, as well as a range of factors that might influence immunisation practices. As such, the objective of this study was to estimate the extent of full immunisation coverage and to investigate the determinants of full immunisation using the nationally representative data of Bangladesh. Analysis of a nationally representative dataset would thus allow for a generation of the current evidence of vaccination practices that could be useful for international comparison and will help in implementing the relevant immunisation policies and priorities for the betterment of the child’s health against VPD in Bangladesh.

Methods

Study design and population

This study was based on a secondary analysis of the latest Bangladesh Demographic and Health Survey (BDHS) 2014 dataset. The DHS is a nationally representative cross-sectional survey, which uses a two-stage stratified sampling design to cover the target population of the entire country. This survey used the sampling frame provided by the Bangladesh Bureau of Statistics. This survey was carried out for 6 months, from June 2014 to November 2014. The BDHS 2014 used three types of questionnaires: household questionnaire, woman’s questionnaire and community questionnaire. Through the women’s questionnaire, up-to-date information on sociodemographic, maternal and child health indications including individual-level vaccination coverage was collected. Participants of this study were children aged between 1 and 5 years. Childhood immunisation history was collected for all surviving children over the last 5 years. Immunisation data were collected mainly from the records on vaccine cards; if vaccine card was not available, mothers were asked to recall about the vaccination history of the respective child. Face-to-face interviews that were conducted with reproductive-aged women (15–49 years) were conducted for the collection of data using the structured questionnaire.22 Details about a sampling technique, survey design and quality control have been described elsewhere.23 All the DHS data were publicly accessible and were made available on request by the MEASURE DHS. Furthermore, approval was sought from and given by the MEASURE DHS programme office in the use of this dataset for this specific study. According to the DHS, written informed consent was obtained from all participants prior to the interviews.

Outcome variable

The outcome variable of the analysis was the children’s immunisation status and categorised as fully immunised and partially/unimmunised. Vaccination status was categorised as ‘fully immunised’ if the children had received doses of all the eight recommended vaccines: one dose of the vaccine against tuberculosis (BCG), three doses of pentavalent (DPT, Haemophilus influenzae type B and HepB), three doses of oral polio vaccine (OPV) and one dose of measles vaccines.23 Children aged 12–59 months of age were included in this analysis to capture the vaccination status of the children. Children younger than 12 months of age were excluded as they were not old enough to receive the full schedule of the routine vaccines. WHO recommended vaccination schedule for Bangladesh24 is presented in table 1. Immunisation histories for all vaccines were coded as dummy responses, with 1 for fully immunised and 0 otherwise. In a small number of cases, the immunisation cards were not available, and mothers responded with ‘don’t know’ while asked about their children’s vaccination status on certain vaccines. In that cases, the vaccination status of children was considered as ‘not fully immunised’ since such responses reflect a negative response regarding immunisation.14

Table 1

The Expanded Program on Immunizations schedule in Bangladesh

Explanatory variables

Explanatory variables were selected based on the published literature, prior knowledge and availability of variables in the BDHS 2014 dataset. Area of residence, administrative division, sex and age of children, birth order, family size, parental education, mass media exposure, household wealth index and the mother’s decision-making ability for both the children and their own healthcare matters were included as explanatory variables. In this analysis, the categorisation of continuous variables was done in light of previous literature where the child’s age was categorised into four groups at 12-month intervals, and maternal age was categorised into three groups (less than 20 years, 20–34 years and more than 34 years). Self-reported parental (both mother and father) educational attainment was used and categorised as ‘no education’, ‘primary’, ‘secondary’ and ‘higher’. No education refers to there not having any formal education, with primary was defined as completing grade 5, secondary as completing grade 10 and higher was defined as attaining more than grade 10. Family size was determined by the number of family members who lived together in a household, and defined by the BDHS. Family size was categorised on the basis of other published literature and categorised as ‘small’, ‘medium’ and ‘large’. A ‘small’ family size refers to the family consisting of less than four members, ‘medium’ as four to six members and ‘large’ if the family member exceeds six members. Mass media access was categorised as ‘yes’ if the family had access to televisions and radios and ‘no’ otherwise. We used the predetermined wealth index category provided in the dataset generated from selected household assets using principal component analysis and classified it into five groups, namely: ‘poorest’, ‘poorer’, ‘middle’, ‘rich’ and ‘richest’. Moreover, both the decision-making ability of the mothers for both their children and personal healthcare were categorised into four groups as ‘herself’, ‘jointly with husband’, ‘husband alone’ and ‘other’.

Statistical analysis

The original dataset comprised 8092 children, aged 0–59 months of age. In the case of our analysis, 1555 children were excluded as they were under the age of 12 months. Moreover, 307 cases of missing data were present, and these were cases where no information was provided (no responses, neither ‘yes’, ‘no’ nor ‘don’t know’) on vaccination information, and thus also excluded from the analysis. Finally, a total of 6230 children were eligible and included in the analysis. A proper sampling weight was used in this analysis to make the sample more representative of the population at the national level. Descriptive statistics, such as frequency with percentages, were executed to represent the background characteristics of study participants, and a proportion with 95% CI was used for presenting the coverage rate of fully vaccinated children. The association between each independent variable and the full vaccination uptake was investigated using a univariate analysis. Multivariable logistic regression that adjusted for all the selected independent variables was also constructed to determine the significant influencing factors for fully vaccination and results were presented in terms of adjusted OR(s) (AOR) with a 95% CI. Before the execution of a multivariate regression model, a bivariate analysis was conducted to trace out the significant factors for full immunisation. Variables significant at p ≤0.05 were included in the multivariate logistic regression analysis. Variance inflation factor (VIF) was calculated to detect any multicollinearity in the model. The low value of average VIF (3.07) confirms no notable multicollinearity among variables. Linktest command through Stata/SE V.13.0 indicated that the constructed model was well specified where the Hosmer-Lemeshow test statistics of goodness of fit indicates the acceptance of the model. Data cleaning, validation and all statistical analyses were performed using Stata/SE V.13.0 (Stata Corporation).

Patient and public involvement

No patients were involved in developing the research question, outcome measure and design of the study. We are unable to disseminate the results of the research directly to study participants.

Results

Descriptive statistics

Distribution of study participants across the sociodemographic variables is presented in table 2. The proportion of full vaccination was 86% (5356 out of 6230 children). Among the participants, 74.7% were from rural areas and 35.1% were from the Dhaka division. Participants were almost equally distributed by sex and age categories. Among the participants, 37.7% was the first-born child, and most of the children’s families (57.9%) consisted of four to six members. Majority of the mothers (80.4%) were aged between 20 and 34 years. 46.1% of the mothers and 29.8% of fathers had secondary-level education. Less than half of the mothers (40.9%) had access to mass media, for example, radio and television. In addition, a majority of the mothers reported that they had the capacity for decision-making with regards to both their own as well as their children’s healthcare.

Table 2

Baseline characteristics of study participants and their immunisation coverage in Bangladesh, 2014 (n=6230)

Vaccination coverage rate among study participants

Full immunisation coverage rate of study participants is also presented in table 2. Findings revealed that full immunisation coverage was slightly higher in urban areas (88.5%) as compared with rural areas (85.1%) and was found to be the highest in the Rangpur division (91.5%) followed by Khulna (89%), Dhaka (88.8%) and the lowest in Sylhet division (69.8%). The immunisation coverage across individual vaccines is shown in figure 1. Full immunisation coverage was similar across sex and age categories of the children. However, first-born child had slightly higher vaccine coverage rates (88.4%) than the other (87.7% for second and 81.6% for third or more birth order).

Figure 1

Coverage rate by individual vaccines (BCG, pentavalent, oral polio vaccine (OPV) and measles) across administrative divisions in Bangladesh (map obtained from open access source and edited according to study findings).

Full immunisation coverage was highest among the children whose parents had higher educational attainment (94.2% and 93.7% for higher educated mother and father, respectively). Vaccine coverage rate was also higher among the children of mothers who have access to mass media. Full vaccination coverage was found to increase in accordance with a higher economic status of the children’s family (73.9% for the poorest to 93.2% for the richest quintile). Full immunisation coverage rate was higher among the children whose mother had her own healthcare decision-making ability, including that of her children’s healthcare. The overall coverage rate of BCG vaccine was the highest (97.1%) followed by OPV-1 (97%), pentavalent 1 (96.6%), OPV-2 (95.8%) and pentavalent 2 (95.4%) where the coverage of measles was the lowest (88%) (figure 2). Furthermore, the coverage of each vaccine by urban and rural areas is also presented in figure 2.

Figure 2

Coverage rate by individual vaccines across place of residence of respondents (urban vs rural) in Bangladesh. OPV, oral polio vaccine.

Factors associated with full immunisation coverage

Table 3 shows factors associated with the full-immunisation coverage among children aged 12–59 months in terms of both unadjusted and AORs. The analysis shows that several sociodemographic factors, like administrative divisions, age of children, family size, parental education, wealth index and child healthcare decision-makers were associated with the full immunisation status in both models. From the univariate analysis, we found that the area of residence, childbirth order, access to mass media and mothers’ healthcare decision-maker were significant factors for full immunisation. Urban children were significantly 1.35 times more likely to be fully vaccinated than rural ones (95% CI 1.13 to 1.60, p=0.001), and the first-born child was 1.71 times more likely to be fully vaccinated than the child of third or more birth order (95% CI 1.45 to 2.03, p<0.001). We also observed that mothers of the children who had access to the mass media had higher odds than their counterparts (OR=2.01; 95% CI 1.72 to 2.36, p<0.001). Similarly, parental education and wealth status played a significant role for full immunisation coverage among children. Findings revealed that children of higher educated mother were significantly 4.93 times (95% CI 3.37 to 7.22, p<0.001) and children of higher educated father were 4.03 times (95% CI 2.98 to 5.44, p<0.001) more likely to be fully immunised than children of uneducated mothers and fathers. Findings also revealed that children from the richest households were significantly 4.27 times (95% CI 3.31 to 5.52, p<0.001) more likely to be fully immunised than the poorest.

Table 3

Unadjusted and adjusted effects of factors that are associated with full immunisation coverage

In the adjusted model, the likelihood of being fully vaccinated was significantly lower in the Sylhet division than all other regions. We found that the number of full immunised children was significantly higher in Rangpur division (AOR = 3.46; 95% CI 2.45 to 4.88, p<0.001) followed by the children from Dhaka (AOR = 2.59; 95% CI 2.05 to 3.28, p<0.001), Khulna (AOR = 2.33; 95% CI 1.62 to 3.33, p<0.001) and Rajshahi (AOR = 1.96; 95% CI 1.46 to 2.64, p<0.001) division. The odds of full immunisation differed across age groups of the children. For instance, children aged 48–59 months had the highest odds of being fully vaccinated (AOR=1.32; 95% CI 1.06 to 1.64, p=0.013) compared with children aged 12–23 months. Family size was also appeared to be a significant factor for full immunisation as medium households (four to six members) often fully immunised more than their counterpart (AOR=1.56; 95% CI 1.32 to 1.86, p<0.001). The findings also revealed that parental education was significantly associated with full immunisation. It was observed that the odds of being full immunisation were increased as the educational attainment of parents increased, and that children from higher educated mother (AOR=1.96; 95% CI 1.21 to 3.17, p=0.006) and higher educated father (AOR=1.55; 95% CI 1.05 to 2.29, p=0.026) were more likely to be full vaccinated compared with those with uneducated parents. Children from the richest households had the highest odds of being fully immunised (AOR=2.20; 95% CI 1.50 to 3.21, p<0.001) followed by the middle wealth quintile (AOR=1.78; 95% CI 1.37 to 2.30, p<0.001). The lowest odds were found among the children whose healthcare decision was made solely by their father (AOR=0.69; 95% CI 0.51 to 0.92, p=0.012) compared with children with healthcare decisions that were made by the mothers.

Discussion

The present study examined the extent of full immunisation coverage and determinants of full immunisation status. We found that the overall full immunisation coverage is impressive (86%) within South East Asia Region.25 The underlying success is due to the nature of a pluralistic health system of Bangladesh, where the public, private sector and non-governmental organisations actively participate to deliver healthcare services. In addition, the introduction of systematic outreach approaches helped to increase accessibility to routine immunisation where community health workers are directly involved in delivering vaccines.26 Further, community mobilisation and public awareness-related programmes have also played a significant role in generating demand for vaccine uptake.27

Despite the success, full immunisation coverage varies across the country. For instance, the full immunisation coverage is the highest in Rangpur while the lowest coverage is observed in the Sylhet region. Our study showed that various factors, such as parental education particularly mothers’ education, age of the children, family size, regional variation and wealth quintiles played a significant role for full immunisation coverage. The importance of parental education in improving child’s health is universally recognised. Similar to earlier studies on childhood vaccination, we found that parental education appeared to be a significant factor for childhood immunisation.14 28 This is due to their better knowledge and understanding of the recommended immunisation schedules than non-educated parents.29 Furthermore, educated parents are likely to be wealthier and have better access to health facilities and immunisation services.28 Like other studies in various settings, we observed that the children of educated mothers were often fully immunised when compared with non-educated mothers.14 28 30 It is well established that there is a positive relationship between the level of education and public health awareness. Therefore, community-based behaviour changes programmes, such as immunisation announcement through radio, television and using local drama and public announcement through miking (loud speaker), should be approached to target uneducated mothers in order to provide a better understanding of the beneficial role of immunisation so that they are encouraged to vaccinate their children in a timely fashion.

We observed that full immunisation coverage was higher among elder children (eg, 48–59 months of age) compared with the younger children (eg, 12–23 months of age). Policy-makers should provide more focus on continuous monitoring and promotional activities to increase the number of full immunisation coverage of all strata.14 Again, the size of the family appeared to be a significant factor of fully immunised children. Those who belonged to larger family size were more likely to be unimmunised. This is similar to the findings in other earlier studies conducted in Indonesia, Greece and Angola, where it was also reported that children form larger family size were less likely to be fully immunised.14 31 32 This might also reflect the socioeconomic status of the households as larger families consume more resources, which is troublesome for resource-poor families and would thus impact on their livelihoods. Therefore, policies should focus on awareness development, especially in targeting the larger family to increase the full immunisation coverage. Wealth status of the household plays a significant role in the immunising of their children. Our results are similar to many earlier studies which have shown that there is a significant positive relationship between wealth index and childhood immunisation completion rates; the chance of having fully immunised children increases in accordance with the wealth index of household.14 28 30 While the immunisation services in Bangladesh are completely free of charge, the indirect cost of vaccination, such as income loss and transportation cost, might be associated with the low demand for vaccination specially for poorest households.33 34 Children from urban areas have been reported to have better immunisation status compared with their rural counterparts. This result confirmed previous findings where similar association was also reported.14 35 This is due to better health services available in urban areas compared with the rural.36

Our study observed that child healthcare decision-making by mothers also contributed, playing a significant role in full immunisation coverage, with the child more likely to be fully immunised if the mother was the decision-maker. This might be due to more awareness of the mothers for their child’s health than the fathers. Therefore, the mothers’ autonomy is essential for their child’s immunisation. However, an earlier study showed that if the parents jointly decided on healthcare decision, their children were often fully immunised.14 Thus, there is a need to disseminate information about women autonomy in making healthcare decisions about their children to increase full immunisation coverage. Based on the findings, community-based behaviour change programme that targets parents might be helpful for developing awareness for their childhood immunisation. Like many low-income and middle-income countries, we also observed that geographical disparities may contribute to the full immunisation coverage in Bangladesh.14 37–39 Our study found that children who lived in Sylhet region were more often unimmunised when compared with other parts of the country. This is may be due to various supply-side and demand-side factors, such as the distance of health facilities and vaccination centres, fragile communication systems in some remote areas, a fear of the side effects of vaccines, religious conservatism, a low level of literacy and even lack of awareness about the benefit of vaccination for their children.39 Therefore, policy should target divisions with low immunisation coverage with an innovative immunisation approach that addresses both supply-side and demand-side barriers. For instance, households that are located away from or have difficulty in accessing immunisation services, especially in hilly areas, might benefit from outreach programme or mobile immunisation strategies. Furthermore, the use of mobile phones may be important vehicles for tracking and improving immunisation coverage in these rural, hard to reach areas of Bangladesh.40

Strengths and limitations

The study has several limitations. The study was based on secondary data where immunisation cards and mother’s report were considered as the source of information for their child’s immunisation status. Therefore, the potential of recall bias in our results cannot be ignored. Nonetheless, the mother’s report is considered to be a valid measure of coverage in the absence of a health card, especially in developing countries.41 The explanatory variables were selected based on previous studies and relied on the information available from the dataset. Therefore, there might be some other potential predictors that might be influenced by full immunisation which were not captured in this study. Despite the limitations, the study results can be generalised at a country level because the study used data from the latest nationally representative household survey. Thus, our findings are still significant and relevant in drawing attention to the health policy-makers in ensuring the benefit of vaccination for the betterment of child health.

Conclusion

Our results identified the presence of disparities for full immunisation coverage across regions and by types of vaccines in Bangladesh. Findings revealed that full immunisation status was significantly associated with regional variation, the age of the children, maternal age, parental education, family size, woman autonomy for child healthcare and wealth quintiles. Our study found that a large number of children from the Sylhet division, poor and larger family, and from lower parental education were not fully immunised. The study also identified that measles immunisation coverage was the lowest among the eight recommended vaccines. Therefore, to achieve maximum success and prevent children from VPD, it is mandatory to address the issues obstructing full vaccination. The findings of this study will serve to inform and support innovative approach for immunisation programmes, and the introduction of relevant policies, including the regular monitoring and evaluation of immunisation coverage particularly for low-performing regions of Bangladesh and targeting various vaccines in order to allow the broader benefit of the immunisation programme to be achieved in all strata of society.

Acknowledgments

ICDDRB is thankful to the Governments of Australia, Bangladesh, Canada, Sweden and the UK for providing core/unrestricted support. The authors would like to thank Health Economics and Financing Research Group of ICDDRB for their earlier comments in this research. The authors acknowledges with gratitude the commitment of Bangladesh Institute of Development Studies to its research efforts.

References

  1. 1.
  2. 2.
  3. 3.
  4. 4.
  5. 5.
  6. 6.
  7. 7.
  8. 8.
  9. 9.
  10. 10.
  11. 11.
  12. 12.
  13. 13.
  14. 14.
  15. 15.
  16. 16.
  17. 17.
  18. 18.
  19. 19.
  20. 20.
  21. 21.
  22. 22.
  23. 23.
  24. 24.
  25. 25.
  26. 26.
  27. 27.
  28. 28.
  29. 29.
  30. 30.
  31. 31.
  32. 32.
  33. 33.
  34. 34.
  35. 35.
  36. 36.
  37. 37.
  38. 38.
  39. 39.
  40. 40.
  41. 41.

Footnotes

  • Contributors ARS and MS participated in the design of the study. All authors performed the analysis and prepared the manuscript. ARS, RA and MS provided data analysis advice and revision of the manuscript. All authors read and approved the final manuscript.

  • Funding The authors have not declared a specific grant for this research from any funding agency in the public, commercial or not-for-profit sectors.

  • Competing interests None declared.

  • Ethics approval This study did not require ethical approval as it used unidentifiable secondary DHS dataset. According to the DHS, written informed consent was obtained from mothers/caretakers on behalf of the children enrolled in the survey. The DHS data are publicly accessible and were made available to us upon request by Measure DHS. No identifiable information was included in the dataset and no attempt was made to identify any individual interviewed in the survey.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Data sharing statement The electronic datasets analysed in this study are available for legitimate research purposes from the Measure DHS website (https://dhsprogram.com/data/).

  • Map disclaimer The depiction of boundaries on the map(s) in this article do not imply the expression of any opinion whatsoever on the part of BMJ (or any member of its group) concerning the legal status of any country, territory, jurisdiction or area or of its authorities. The map(s) are provided without any warranty of any kind, either express or implied.