Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Trend and determinants of quality of family planning counseling in Ethiopia: Evidence from repeated PMA cross-sectional surveys, (2014–2019)

Abstract

Introduction

The modern contraceptive prevalence rate (mCPR) among married women has increased by nearly five-fold in Ethiopia from 8.1% in 2000 to 37% in 2019. Despite this increase, receipt of high quality contraceptive counselling, as measured by the percentage of contraceptive users who were told about other methods, counseled on side effects and counseled what to do in the event that they encountered side effects, has declined in recent years. The quality of family planning counseling service measured by using these three components, known as the Method Information Index(MII), is an index designed to measure quality and a key indicator of the FP2020 initiative. The effects of potential client and service provider-level factors on receipt of high quality counseling and its progress over time have not been well studied in Ethiopia.

Methods

We pooled data from seven Performance Monitoring for Action (PMA), formerly PMA2020, survey-rounds to examine the trend and effect of potential factors on receiving high quality of family planning counseling service in Ethiopia. Data from a total of 15,597 women aged 15 to 49 from seven survey-rounds were used in the analysis. To account for the study design and unequal probabilities of selection from target-populations for sampled women, design-based analysis was used to compute proportions. Multilevel ordinal regression model with enumeration area as a second level were employed to examine potential factors associated with quality of family planning counseling service.

Results

We found that the percentage of women who received high quality family planning counseling service declined from 39% (95%CI: 33%, 44%) in 2015 to 12% in 2019 (95% CI: 10%, 14%) nationally. Amhara region had the lowest percentage of women receiving high quality counseling at both the earliest(2014) and latest(2019) survey rounds(17% and 6%, respectively). Results show that lack of media exposure about family planning, having no formal education, using short-acting methods, and getting the service from pharmacy were the main factors associated with receiving low quality family planning counseling service.

Conclusions

Given the importance of continuous provision of information on the range of family planning methods, it is imperative to use media and particularly regional media which can effectively address the rural populations in local languages as an important vehicle of information on family planning. Interventions aimed at improving quality of family planning counseling need to be mindful of regional disparities in the severity of the problem to ensure equity in service access. To improve the coverage of high quality family planning counseling service, there is an urgent need to re-visit the format of family planning counselling services.

1 Background

The 2030 Sustainable Development Goals (SDGs) calls on countries “by 2030, to ensure universal access to sexual and reproductive health-care services, including for family planning, information and education, and the integration of reproductive health into national strategies and programmes” [1]. In 2019, among 1.9 billion women of reproductive age (15–49 years) worldwide, 1.1 billion had a need for family planning [2]. Of these 842 million were users of modern contraceptive methods and 190 million women who wanted to avoid pregnancy were not using any contraceptive method [1]. In sub-Saharan Africa, annual rate of changes (percentage points) in modern contraceptive use among women of reproductive age varied substantially from 0.77 to 3.17 [3].

Counseling plays a key role in enhancing family planning services, providing contraception related information and supporting family planning and fertility goals for women of childbearing age [4]. The Method Information Index (MII) is the current standard proposed by the Track 20 initiative to monitor quality of family planning counseling [5, 6]. The simple index uses three questions—whether the respondent received counseling on multiple contraceptive methods, whether she was told about side-effects, and if so, whether she was told what to do about side-effects—that are available in the majority of population-based surveys, including the Demographic and Health Survey (the DHS) and Performance Monitoring for Action (PMA) [7]. Adaptations have been proposed to the MII, including a fourth question assessing whether women were told they could switch to another method. This question was added in PMA 2019 Survey, however, and for trend purposes is not included in the current analysis.

Despite substantial global efforts to improve conceptualization, measurement, and delivery of high-quality family planning services, studies consistently find that receipt of high-quality services is low. In a study exploring the levels and trends of the MII across 25 countries, urban residence, education, household wealth, and method type were generally found to be associated with receiving higher quality care [8]. In all countries combined, the median MII values show that the contraceptive information received by women varied by method type and was highest for implant users and lowest among women relying on sterilization.

In Ethiopia specifically, despite increases in contraceptive coverage over time the percentage of contraceptive users who received all components of Method Information Index(MII) has not changed substantially [9, 10]. Studies on family planning counseling services in Ethiopia consistently show that high quality of family planning counseling service coverage is below 40% [1114].

Understanding which individual, social, and health system factors are associated with the receipt of high-quality family planning counseling is critical in order to identify inequities and improve overall service delivery. Multiple studies have assessed facility and provider characteristics that are associated with the receipt of quality family planning counseling services [1518], but relatively few studies have assessed the client-level and social factors that may be related to receipt of high quality counseling service.

A previous study in Ethiopia found that in addition to these factors, women who received family planning services from a public provider had higher odds of receiving all components of the MII. Little additional evidence exists, however, on the client-level factors that may be associated with receiving high quality services, resulting in a call for additional research into identifying factors associated with quality [19].

Studies that have explored client-level and social factors associated with higher quality services have generally been limited by use of one to two cross-sectional surveys which limit the ability to detect trends in receipt of high quality family planning counseling service over time. Increasingly, however, it is common to see repeated cross-sectional surveys carried out in multiple rounds over short periods of time, particularly surveys conducted by the Performance Monitoring for Action (PMA) project. Pooling such data improves the power of test due to an increased sample size and can be used to better investigate trends over-time. A study done by [12] used pooled data from five survey rounds in Ethiopia to compute the percentage of women who received high counseling service across survey years by different factors. However, the analysis on factors of quality of family planning counseling service was limited only to one survey round (i.e., 2018 survey). To our knowledge, there is no study that considered repeated cross-sectional surveys to identify factors related to quality of counseling service in Ethiopia.

In the present study, we use repeated cross-sectional data from seven-rounds of a national survey to: i) establish the trends in quality of family planning counseling among all women aged 15–49 years, ii) identify client-level and social factors associated with quality of counseling, and iii) identify regional disparities in quality of counseling. Findings from this study can help to inform planning interventions to improve the quality of family planning counseling services in Ethiopia.

2 Materials and methods

2.1 Study design and data source

The data for this study were obtained from seven rounds of the Performance Monitoring for Action (PMA) surveys, which were undertaken yearly since 2014 (formerly as Performance Monitoring and Accountability 2020 (PMA2020)). The PMA surveys are repeated cross-sectional surveys based on a multistage stratified cluster sampling design with urban-rural stratification for Amhara, Oromia, Tigray, and Southern Nations, Nationalities, and Peoples’ (SNNP) regions, and for the remaining regions, regions served as the strata, without additional urban/rural stratification. The primary sampling units or Enumeration Areas (EAs) are selected using probability proportional-to-size method for which the sample selection probability depends on the size of population. For the first round (referred to as PMA2014/Ethiopia) and latest round considered in this study (referred to as PMA2019/Ethiopia) of data collection the survey conducted in 200 and 265 enumeration areas, respectively, which were selected by the Central Statistics Agency (CSA) to be representative at the national and regional level. Prior to data collection, all households in each enumeration area were listed and mapped by the resident enumerators to create a frame for the second stage of the sampling process. This mapping and listing process was done in the first week of data collection in each enumeration area. Once listing was completed, 35 households were randomly selected by field supervisors using a phone-based random number-generating application. All occupants in selected households were enumerated and from this list, all eligible women between the ages of 15 to 49 were administered the female-household questionnaire after informed consent was obtained. In PMA surveys considered in this analysis, the response rate is high across different survey years which ranges from 95.6% to 99.2%. In all survey-rounds, data were collected using smartphones.

The details of the PMA2020 survey design and sampling techniques are described elsewhere [20], while updates made to the sampling design in 2019 are further described by Zimmerman and colleagues [21]. PMA uses standard data collection procedures which are consistent across-rounds, and consistent in content, allowing pooling the data and making comparisons over time. To assess trend over time and identify factors linked with the outcome variable, described below the pooled data from seven-rounds, spanning 2014 to 2019 were considered.

2.2 Sample

In this study, data from nationally representative cross-sectional female and household surveys collected on yearly bases from 2014 to 2019 among women aged 15–49 were utilized. A total of 15,597 women who were asked about the three questions, described below, were considered as our analytical sample.

2.3 Variables

The outcome variable in this study, quality of counseling service, was derived based on the responses to the three equally weighted questions used in the MII. Women were asked to report whether, at the visit where they first received their current contraceptive method, they were: i) told about other methods, ii) counseled on side effects, and iii) counseled on what to do if she experienced side effects. Originally, the data related to the variable of interest were collected using possible Yes or No questions to assess quality of counseling. Using these questions, the MII was assessed two separate ways; first, as a binary variable, indicating whether a woman received all components of the MII (“high-quality”) or not, and second, as an ordinal variable. We defined the ordinal variable with four categories indicating “high” quality (all three components received), “moderate” quality (any two of three), “low” quality (any one of three), and “no” counseling (no components received).

In order to assess the effect of different factors on receiving high quality family planning counseling service, the outcome variable for the ith women from jth enumeration area(EA), Yij, is dichotomized as follows:

This categorization was used to compute design-based prevalence of high quality family planning counseling service by different socio-demographic factors, described further below.

The analysis using the dichotomized outcome cannot provide information about different counseling levels simultaneously, thus further disaggregation was done based on the following categorization: (1)

After reviewing related literature, possible factors linked with the quality of family planning counseling were extracted from the survey data. Table 1 presents the list of factors we hypothesize are associated with the receipt of quality of family planning counseling service considered in this study. We define long-acting methods as female sterilization, implant and IUD and short-acting methods as the injectable, pill, emergency contraception, male condom, the Standard Days Method, and lactational amehorrhea method (LAM).

thumbnail
Table 1. Description of exploratory variables considered in the study.

https://doi.org/10.1371/journal.pone.0267944.t001

2.4 Statistical analysis

Nationally representative survey data-sets are often complex in nature for two reasons: i) the use of stratified multistage cluster sampling, and ii) unequal probabilities of selection from target-populations for sampled elements, often as a result of oversampling of key subgroups. Thus to account for the study design and unequal probabilities of selection from target-populations for sampled women, sampling weights were used during the estimation of proportions and their respective confidence interval [22].

The types of models used for data analysis depend on the nature and measurement scale of the outcome variable under consideration. The level of measurement of the variable is both binary and ordinal. When the response categories are ordered, the use of this ordering yields more parsimoniously parameterized models [23]. In such models, the ordinal nature of the response is taken into account by considering the cumulative probabilities. The data from different women are assumed to be independent, but due to the clustered nature of the data, observations from the same cluster/EA may not be independent.

Women living in the same EA may share similar health facilities, and the quality of counseling services available to them may be more alike than women living in another EA. Hence, our modeling approach should takes into account the correlated nature of the data by introducing a random effects term to account for a variety of situations, including EA heterogeneity, and unobserved covariates within that EA. As a result, in this study we employ a random effects model for ordinal outcomes.

Let μij be the probability of the ith subject from the jth EA being in the response category k, μij = P(yij = k). Further, we let cumulative probability of the response in category k or above represented by Πijk = P(yijk). The random intercept proportional odds model (POM) is given by: (2) where k is the level of the ordered category. The parameter β0k is the intercept for category k, usually is considered as nuisance parameter of little interest ([23]), Xij is a vector of exploratory variables described in Table 1 with associated regression coefficients β. The random intercept bj is assumed to vary randomly among EA’s according to a distribution. The variance of the random effect () represents how much variability there is between clusters with respect to the response variable.

Model 2 assumes an identical effect of the predictors for each cumulative probability. The model assumes that, conditional on the random effects, the coefficients that describe the relationship between, for example, the lowest versus all higher categories of the response variable are the same as those that describe the relationship between the next lowest category and all higher categories, etc. This is called the proportional odds assumption or the parallel regression assumption [23]. Specifically, the model implies that odds ratios for describing effects of explanatory variables on the response variable are the same for each of the possible ways of collapsing the response to a binary variable. The proportional odds (PO) assumption for each variable were tested using the Brant test ([24]), and held across all variables.

Further, to quantify the relative importance of clusters as a source of variation about receiving high counseling service, EA-level random effects variance was expressed in terms of Variance Partition Coefficients (VPC).

The proportion of total variation of receiving high family planning counseling attributable to EA random effect () is computed as follows [25, 26]: where is EA (level-2) random effect variance, and π2/3 = 3.29 is individual level variance is equal to the variance of a logistic distribution [27]. In Model 2, fixed effects parameter has a conditional interpretation. It refers to the consequence of changing the value of an explanatory variable, for which the fixed effect is the coefficient, for a given value of the random effect and the other fixed effects.

In all exploratory data analyses, sample weights calculated based on the multistage sampling design were considered. Data analyses were done using Stata 16.1 software [28], and graphs presented in this manuscript were generated using the R software [29] by ggplot2 package [30]. In this study, for statistical significance, the considered alpha(type-I error) is 0.05.

2.5 Ethics approval statement

The study protocol was approved by Institutional Review Boards of College of Health Science at Addis Ababa University, Ethiopia and at the Bloomberg School of Public Health at Johns Hopkins University in Baltimore, USA.

3 Results

3.1 Descriptive results

Table 2 presents client, service provider, and regional level characteristics of the sampled study population by year. Across all years, the majority of women lived in rural areas. Approximately 50% of respondents were age 25–34 years, and less than 20% of respondents were nulliparous. On average, around 50% of women got received their contraceptive methods from health centers. The use of long-acting contraception methods (female sterilization, implant and IUD) increased from 21.81% in 2014 to 36.94% in 2019. The percentage of respondents who reported hearing a message about family planning from the media decreased across survey years.

thumbnail
Table 2. Sample characteristics of PMA survey respondents across survey years.

https://doi.org/10.1371/journal.pone.0267944.t002

The overall weighted prevalence of receiving high quality family planning counseling service was 28.61% (with 95% CI: 27.68%—29.55%). Table A1 in S1 Appendix presents the percentage distribution of women who received each level of family planning counseling quality by different factors and over the survey-rounds. The percentage of women who received high quality family planning counseling service was 27.6% in 2014, increased to 38.6% in 2015, and thereafter declined to a low of 11.8% in 2019 (Table A1 in S1 Appendix, Fig 1). Receiving a method from a health center or health post, using a long-acting method, and having media exposure to family planning were associated with a higher prevalence of receiving high quality family planning counseling service across all survey rounds.

thumbnail
Fig 1. Percentage of women receiving different levels of family planning counseling service over survey years.

https://doi.org/10.1371/journal.pone.0267944.g001

Table 3 presents a summary of the percentage change in receiving high quality family planning counseling service. Compared with 2018, there was a significant decline (by 59.4%) in the coverage of providing high counseling service in 2019. Overall, the percentage of the women who received high quality family planning counseling service declined between 2014 and 2019 by 57.3 percentage points.

thumbnail
Table 3. Percentage change in receipt of high-quality family planning services across survey rounds.

https://doi.org/10.1371/journal.pone.0267944.t003

The percentage of high quality family planning counseling changed over time with different rate across regions of the country (Fig 2, Table A1 in S1 Appendix). The coverage in the receipt of high quality of counseling in Tigray, Addis Ababa and SSNP region was above the national level while, the coverage was found to be lower in Amhara region (Fig 2).

thumbnail
Fig 2. Percentage of women receiving high quality family planning counseling service over survey years by region.

https://doi.org/10.1371/journal.pone.0267944.g002

Fig 3 presents overall trends in the percentage of women who received high quality family planning counseling service across survey year by different factors. The results reveal that, irrespective of the factors examineds, the quality of counseling increased from the year 2014 to 2015 and start declining thereafter.

thumbnail
Fig 3. Trends in receipt of high quality family planning counseling by different factors.

Percentage of women who received high family planning counseling service according to: a) type of method used, b) exposure to family Planning via media, and c) source of contraceptive method.

https://doi.org/10.1371/journal.pone.0267944.g003

3.2 Modeling results

Table 4 presents the multilevel modeling results by considering one factor at a time (Model A) and after controlling for observed explanatory variables (Model B). Conditional on EA level random effects, the odds of getting counseling service varied across categories for parity, method type, education, media exposure, method source, region of residence, and survey year.

thumbnail
Table 4. Parameter estimates from the multilevel ordinal logistic regression model of receiving high quality family planning counseling service in Ethiopia.

https://doi.org/10.1371/journal.pone.0267944.t004

Women who reported being exposed to family planning information through the media, those who received their family planning method from health post or health center, women who attained secondary or higher education, and women using long-acting methods were significantly more likely to receive high counseling service on family planning (Table 4).

Among women with the same random effect value, relative to women who were using short-acting methods, women who were using long-acting methods had 1.90 times the estimated odds of getting “high” counseling service instead of “moderate” or “low” or “no” counseling and the estimated odds of getting “high” or “moderate” counseling service instead of “low” or “no” counseling service and the estimated odds of getting “high” or “moderate” or “low” counseling service instead of “no” counseling service. Among women with the same random effect, compared to women who got their method from a hospital, women who got their method from a health post had 2.18 times the estimated odds of getting “high” counseling service instead of “moderate” or “low” or “no” counseling service and the estimated odds of getting “high” or “moderate” counseling service instead of “low” or “no” counseling service and the estimated odds of getting “high” or “moderate” or “low” counseling service instead of “no” counseling service (Table 4).

For Model B (Table 4), the variance partition coefficient (VPC) is , which indicates that 17.63% if the variance in counseling service can be attributed to differences between EA’s.

4 Discussion

This study examined factors associated with receipt of high quality family planning counseling services and its trend in Ethiopia from 2014 to 2019. We found that the percentage of women who received high-quality services was low throughout the study period, reaching a high of 38.6% in 2015 and declining sharply to 11.8% in 2019. Significant variation in receipt of high quality counseling was found at the client-level, service delivery point-level, and by region of residence. Nationally, our finding on the trend analysis in the percentage of women receiving high quality family planning counseling service declined overtime. This findings point out some of the gains in improving service coverage were not complimented with improved service quality. In addition, the fact that exposure to family planning message from the media and health professionals has decreased over the years indicate perhaps less attention was given to maintaining and improving the existing public awareness on family planning methods. While the level of high-quality counseling did not exceed 40% at any point, there was a sharp decline in 2019, falling almost 60% from the year before. Among client-level factors, parity, education, media exposure, and type of contraceptive method were found to be statistically significantly associated with high quality family planning counseling services.

Except in Addis Ababa, high-quality family planning counseling peaked in 2015 in all other regions and declined in the following years. The observed relatively better family planning counseling service in 2015 may be due to the fact the Ethiopian health sector transformation plan started in this year [31]. As compared with other regions, high quality counseling service was low in Amhara region and sharply declined since 2017 (Fig 2). This aligns with Hrusa’s finding using PMA data [12].

At the client-level, Jain found that education and wealth were generally associated with higher quality services [8]. We found no differences by wealth, though we did see that women who attended secondary or higher education were significantly more likely to report receiving higher quality services than women without education. Similarly, women who heard about family planning in the media were also more likely to indicate they received higher quality counseling services. These results provide support for Jain’s hypothesis that women with more education may have greater method knowledge and greater ability to recall the information exchanged during their family planning visit. Engagement with family planning media may also lead to greater knowledge of different contraceptive methods and their side effects, allowing women to ask questions of providers that may otherwise be overlooked. Continued training to ensure that comprehensive contraceptive information is delivered in simple, understandable language may alleviate some educational disparities, as would investment in ongoing family planning media campaigns. Messaging that is developed in conjunction with regional media, and distributed in local languages, may support greater general knowledge of family planning and support women in asking questions and receiving high quality services. Media exposure enables a more informative client-provider interaction because it’s enhances client knowledge and empowers them to be more inquisitive in the family planning visit.

In line with findings in other countries, education, parity have impact on receiving high family planning counseling service [8, 15, 3234]. A study by [14] showed, 58.9% of the counseling sessions did not maintain the privacy of clients during the consultation and 74.2% of them were not told about the possible side effects related to the use of a method. A study done in the northern part of Ethiopia shows that only 22.4% of the clients chose contraceptive method by themselves and 74.1% of them with assistance of the provider [13]. Variation in access to and quality of family planning services by parity and method-type has been documented extensively across settings and has identified consistent barriers to high quality family planning services among unmarried adolescents, including provider bias and lack of training [14, 35, 36]. Policies and programs designed to expand and strengthen youth friendly services must emphasize that the delivery of comprehensive information is a critical component of high-quality youth friendly services, in addition to improving privacy, confidentiality, and reducing provider bias.

Finally, women who reported using a long-acting method were significantly more likely to report receiving high quality services than women using short-acting methods. This aligns with Jain’s findings that the MII was highest among implant users, followed by the IUD [8]. While it is positive that women who choose long-acting methods are made aware of potential side effects and how to address them, the majority of contraceptive users in Ethiopia use the injectable [9]. Even as the method mix changes in Ethiopia, with greater use of long-acting methods, it is important to ensure that women who continue to prefer short-acting methods are given comprehensive information.

At the service provider-level, method source was significantly associated with quality of family planning counseling service. Women who received care from pharmacies had significantly lower odds of receiveing high-quality family planning counseling than women who received care at either a health post or health center. In addition to differences in service modality and training of service providers, this likely reflects differences in characteristics of users and method choice; other studies have shown that women who use pharmacies tend to live in urban areas, be younger, and rely on short-term methods such as emergency contraceptives and condoms which are typically not preceded by sufficient counseling [37, 38].

We did not find any differences by urban or rural residence, which is contrary to the recent Hrusa study in Ethiopia that found rural women were more likely to receive high quality services [12]. This may be due to the fact that the paper did not account for clustering at the EA level, which may be particularly important among rural women who have fewer health facilities available to them than urban women. While Jain found that across 25 countries, urban women were more likely to receive high quality counseling, Ethiopia has invested heavily in its public health system, particularly in rural areas, to address disparities in access to and use of family planning services [31]. This, combined with the generally higher level of services available at public health facilities (i.e., health post and health centers), may explain why urban and rural differences are not pronounced.

The results presented in this study should be considered in light of some limitations. One limitation is the fact that the MII does not capture the full range of elements necessary to measure high-quality contraceptive counseling. Additionally, due to the nature of the household survey which asks women to recall information about the first time they received their current method, there may be a recall bias. Further, covariates considered in this study were limited to those available in the survey which lacks supply-side variables such as availability of methods at facilities, availability of trained staff, job aids, information leaflets for clients and other facility level variables. Since the analysis results in this study were derived from cross-sectional female and household surveys, and quality of counseling is heavily influenced by facility-level factors, future studies should consider using PMA service delivery point (SDP) surveys.

5 Conclusions and recommendation

Quality of family planning services is low in Ethiopia, with evidence of a sharp decline in recent years. Significant disparities by region, method source and by socio-demographic characteristics exist, particularly among uneducated women, media exposure about family planning and by choice of method. Given the importance of continuous provision of information on the range of family planning methods and their side effects, it is imperative to develop simple, culturally appropriate messages in local languages and collaborate with regional media as an important vehicle of information on family planning. Interventions aimed at improving quality of family planning counseling need to be mindful of the regional disparities and ensure equity in service quality. Further, the results of this study reveal that high-quality family planning counseling service tended to be more available at health posts and health centers compared to hospitals and is particularly poor at pharmacies. Thus, the Ministry of Health and partners should support training on contraceptive counseling for pharmacists and related auxiliary health staff to improve the information provided through the private sector.

Acknowledgments

The authors greatly appreciate the comments of Colin Baynes and Orvalho Augusto. The paper has been considerably strengthened due to their suggestions.

References

  1. 1. United Nations, Department of Economic and Social Affairs, Population Division (2019). Family Planning and the 2030 Agenda for Sustainable Development: Data Booklet.; 2019. Available from: http://creativecommons.org/licenses/by/3.0/igo/.
  2. 2. Gujo AB, Kare AP. Utilization of Long-Acting Reversible Contraceptives and Associated Factors Among Reproductive Age Women Attending Governmental Health Institutions for Family Planning Services in Wondo Genet District, Sidama, National Regional State, Southern Ethiopia. Health services research and managerial epidemiology. 2021;8. pmid:33796626
  3. 3. Ahmed S, Choi Y, Rimon JG, et al. Trends in contraceptive prevalence rates in sub-Saharan Africa since the 2012 London Summit on Family Planning: results from repeated cross-sectional surveys. Lancet Global Health. 2019;7(7):e904–e911. pmid:31109881
  4. 4. PATH. PATH and United Nations Population Fund. Meeting the Need: Strengthening Family Planning Programs; 2016.
  5. 5. Aparna J, Kumudha A, Elizabeth T, Arupendra M, Rajib A. Adding a Question About Method Switching to the Method Information Index Is a Better Predictor of Contraceptive Continuation. 2019;7(2):289–299.
  6. 6. Chang KT, Mukanu M, B B, Hameed W, Kalamar AM, Grépin KA, et al. Evaluating Quality of Contraceptive Counseling: An Analysis of the Method Information Index. Studies in family Planning. 2019;50:25–42. pmid:30666641
  7. 7. Jain AK, Townsend J, S R. Proposed metrics to measure quality: Overview; 2018. Available from: https://www.pmadata.org/countries/ethiopia.
  8. 8. Jain A. Examining progress and equity in information received by women using a modern method in 25 developing countries. International perspectives on sexual and reproductive health. 2016;42(3):131–140. pmid:28825904
  9. 9. Shiferaw S, Seme A, Zimmerman L. PMA Ethiopia Results Brief, 2019; 2019. Available from: https://www.pmadata.org/countries/ethiopia.
  10. 10. Shiferaw S, Seme S, Yirgu R, Ejigu B. Examining trends and associated factors of quality of counseling on family planning in Ethiopia: Further analysis of PMA data and a qualitative study; 2020.
  11. 11. Teshome A, Birara M, Romins S. Quality of family planning counseling among women attending prenatal care at a hospital in Addis Ababa, Ethiopia. International Journal of Gynecology and Obstetrics. 2017;137(2). pmid:28170078
  12. 12. Hrusa G, Spigt M, Dejene T, Shiferaw S. Quality of Family Planning Counseling in Ethiopia: Trends and determinants of information received by female modern contraceptive users, evidence from national survey data. PLOS ONE. 2020;15(2):e0228714. pmid:32040485
  13. 13. Fantahun M. Quality of family planning services in northwest Ethiopia. EthiopJHealth Dev. 2005;19(3):195–202.
  14. 14. Abdulreshid M, Dadi H. Assessment of Family Planning Counseling Provided for Postpartum Women and Associated Factors. International Journal of Reproductive Medicine. 2020;. pmid:32047803
  15. 15. Assaf S, Wang W, Mallick L. Quality of care in family planning services at health facilities in Senegal; 2015.
  16. 16. Johnson D, Ugaz J. Determinants of Quality of Family Planning Counseling among Private Health Facilities in Lagos. Studies in Family Planning. 2016;47(3):239–251. pmid:27482638
  17. 17. Okour AM, Saadeh RA, Zaqoul M. Evaluation of Family Planning Counselling in North Jordan. Sultan Qaboos University medical journal. 2017;17(4). pmid:29372086
  18. 18. Puri MC, Moroni M, Pearson E, Pradhan E, Shah IH. Investigating the quality of family planning counselling as part of routine antenatal care and its effect on intended postpartum contraceptive method choice among women in Nepal. BMC women’s health. 2020;20(1). pmid:32070339
  19. 19. Tessema GA, Gomersall JS, Mahmood MA, Laurence CO. Factors Determining Quality of Care in Family Planning Services in Africa: A Systematic Review of Mixed Evidence. Plos ONE. 2016;11(11). pmid:27812124
  20. 20. Zimmerman L, Olson H. PMA2020 Principal Investigators Group, Tsui A, Radloff S. PMA2020: Rapid Turn-Around Survey Data to Monitor Family Planning Service and Practice in Ten Countries. Stud Fam Plann. 2017;48(3):293–303. pmid:28885679
  21. 21. Zimmerman L, Desta S, Yihdego M, Ann Rogers A, Amogne A, et al. Protocol for PMA-Ethiopia: A new data source for cross-sectional and longitudinal data of reproductive, maternal, and newborn health [version 1; peer review: 2 approved]. Gates Open Res. 2020;4(126). pmid:33150302
  22. 22. Carle AC. Fitting multilevel models in complex survey data with design weights: recommendations. BMC Medical Research Methodology. 2009;9:49–62. pmid:19602263
  23. 23. Agresti A. Categorical Data Analysis. New York: John Wiley & Sons; 2012.
  24. 24. Brant R. Assessing Proportionality in the Proportional Odds Model for Ordinal Logistic Regression. Biometrics. 1990;46(4):1171–1178. pmid:2085632
  25. 25. Leckie G, French R. Module 11: Three-level multilevel models MLwiN practical. LEMMA VLE Module. 2013;11:1–56.
  26. 26. Leyland AH, Groenewegen PP. Multilevel Modelling for Public Health and Health Services Research. Gewerbestrasse 11, 6330 Cham, Switzerland: Springer, Cham; 2020.
  27. 27. Evans M, Hastings N, Peacock B. Statistical Distributions. New York: John Wiley & Sons; 1993.
  28. 28. StataCorp. Stata Statistical Software: Release 16.1; 2020. Available from: https://www.stata.com/.
  29. 29. R Core Team. R: A Language and Environment for Statistical Computing; 2020. Available from: https://www.R-project.org/.
  30. 30. Wickham H. ggplot2: Elegant Graphics for Data Analysis. Springer-Verlag New York; 2016. Available from: https://ggplot2.tidyverse.org.
  31. 31. MoH. Health Sector Transformation Plan; 2015. Available from: https://www.globalfinancingfacility.org/sites/gff_new/files/Ethiopia-health-system-transformation-plan.pdf.
  32. 32. Mensch B, Arends-Kuenning M, Jain A. The impact of the quality of family planning services on contraceptive use in Peru. Stud Fam Plann. 1996;27(2):59–75. pmid:8714304
  33. 33. Chakraborty NM, Chang K, Bellows B, Grépin KA, Hameed W, Kalamar A, et al. Association Between the Quality of Contraceptive Counseling and Method Continuation: Findings From a Prospective Cohort Study in Social Franchise Clinics in Pakistan and Uganda. Global Health: Science and Practice. 2019;7(1):87–102. pmid:30846566
  34. 34. Hossain S, Sripad P, Zieman B, et al. Measuring quality of care at the community level using the contraceptive method information index plus and client reported experience metrics in Bangladesh. Journal of Global Health. 2021;11. pmid:33763221
  35. 35. Brittain AW, Loyola Briceno AC, Pazol K, Zapata LB, Decker E, et al. Youth-Friendly Family Planning Services for Young People: A Systematic Review Update. American journal of preventive medicine. 2018;55(5):725–735. pmid:30342635
  36. 36. Tilahun M, Mengistie B, Egata G, et al. Health workers attitudes toward sexual and reproductive health services for unmarried adolescents in Ethiopia. Reprod Health. 2012;9(19). pmid:22943476
  37. 37. Radovich E, Dennis ML, Wong KLM, Ali M, Lynch CA, Cleland J, et al. Who Meets the Contraceptive Needs of Young Women in Sub-Saharan Africa? Journal of Adolescent Health. 2018;62(3):273–280. pmid:29249445
  38. 38. Corroon M, Kebede E, Spektor G, Speizer I. Key Role of Drug Shops and Pharmacies for Family Planning in Urban Nigeria and Kenya. 2016;4(4):594–609.