Objectives. We analyzed the likelihood of rural children (aged 6-24 months) being stunted according to whether they were enrolled in Mutuelles, a community-based health-financing program providing health insurance to rural populations and granting them access to health care, including nutrition services. Methods. We retrieved health facility data from the District Health System Strengthening Tool and calculated the percentage of rural health centers that provided nutrition-related services required by Mutuelles’ minimum service package. We used data from the 2010 Rwanda Demographic and Health Survey and performed multilevel logistic analysis to control for clustering effects and sociodemographic characteristics. The final sample was 1061 children. Results. Among 384 rural health centers, more than 90% conducted nutrition-related campaigns and malnutrition screening for children. Regardless of poverty status, the risk of being stunted was significantly lower (odds ratio = 0.60; 95% credible interval = 0.41, 0.83) for Mutuelles enrollees. This finding was robust to various model specifications (adjusted for Mutuelles enrollment, poverty status, other variables) or estimation methods (fixed and random effects). Conclusions. This study provides evidence of the effectiveness of Mutuelles in improving child nutrition status and supported the hypothesis about the role of Mutuelles in expanding medical and nutritional care coverage for children.
Our target population was rural children in Rwanda, who had a much higher prevalence of stunting than urban children in 2010 (47% vs 27%).25 We focused on rural children aged 6 to 24 months for the following reasons. First, it has been suggested that stunting is hard to reverse for children older than 2 years, and the first 2 years of life is the “window of opportunity” for interventions against stunting.2,31–34 Second, about 93% of rural children younger than 6 months had exclusive breastfeeding, and the prevalence of stunting for this age group was much lower (18%) than for the children aged 6 to 24 months (42%).25 We used data from the RDHS 2010 to conduct an individual-level analysis. The RDHS is a nationally representative, population-based survey conducted every 5 years to measure indicators of population health and nutrition, with a special emphasis on mothers and on children younger than 5 years. The RDHS also collects information on households’ and mothers’ sociodemographic characteristics, health insurance status, and utilization of health care. It has been widely used to provide national and regional evidence to policymakers in Rwanda.27,35,36 The RDHS 2010 collected individual information through a 2-stage sampling process. Villages, or primary sampling units, were selected at the first stage. Households in the selected villages were chosen at the second stage.37 We used the DHHST to study the availability of nutritional services included in the MSP across rural health centers in Rwanda. The DHHST is an ongoing Web-based database system built by the Rwandan Ministry of Health in 2009 for monitoring and strategic planning on strengthening health systems.36 The DHHST requires rural health centers to report to the database on an annual basis. The data provide information on (1) medical services (including nutritional care), (2) capacity building (infrastructure and staffing), and (3) revenues and expenditures. The RDHS 2010 measured the heights of children younger than 5 years from a randomly selected 50% subsample of households. Among the 3474 rural children younger than 5 years with a valid height measure, 1087 were aged 6 to 24 months. To identify the link between children’s Mutuelles membership and their nutritional status, we included only children who were either enrolled in the Mutuelles program or had no insurance. We excluded 23 rural children who reported other types of insurance. The final sample size was 1061 children: 838 enrolled in Mutuelles and 223 uninsured. In the DHHST data, of 389 rural health centers in Rwanda, 384 reported service provision in 2010. For individual-level analysis, we constructed a dichotomous variable to indicate whether a child in our study population was stunted, defined as height-for-age 2 standard deviations below the median of the international reference population recommended by the World Health Organization in 2006.38 Facility-level analysis identified a list of services that were included in the MSP of the Mutuelles and available in rural health centers. We treated each variable as binary and assigned a value of 0 or 1. The DHHST surveys included 30 questions on nutrition services in 2010: 10 promotional services, 7 preventive services, and 13 curative services. Online Panel B (available as a supplement to the online version of this article at http://www.ajph.org) lists these nutrition services and how they were delivered by community health workers. For individual-level analysis, we constructed a dichotomous variable to indicate whether a child was enrolled in Mutuelles or uninsured. According to the Mutuelles legislation enacted in 2008,35 enrollees were entitled to access the nutritional care (listed in online Panel B) when the services were available in rural health centers. Sociodemographic variables included a child’s gender, maternal characteristics (mother’s age, education, and height), and household wealth status. We constructed 2 dummy variables that accounted for a mother’s completion of primary school and age older than 30 years. Previous studies found that a mother’s height was significantly associated with her child’s stunting status.39 Following Özaltin et al.,39 we constructed a categorical variable to indicate the mother’s height (< 150.0 cm, 150.0–154.9 cm, 155.0–159.9 cm, or ≥ 160.0 cm). The RDHS 2010 had a wealth quintile variable that summarized a household’s assets (e.g., motorcycle), housing construction (e.g., floor), water source, and sanitation. Because about 45% of Rwanda’s rural population lived below the national poverty line (defined as US $0.45 per day per adult) in 2010,24 we regrouped the households into 2 groups: below the poverty line (households in the lowest and next-to-lowest wealth quintile) and above the poverty line (all others). We constructed a dummy variable to indicate whether a child was from a household living below the poverty line. Note that information about a household’s water and sanitation was reflected in its wealth quintile, and we did not construct separate variables to represent these factors. To determine whether Mutuelles enrollment was associated with the nutritional status of children living below the poverty line versus those living above the poverty line, we constructed an interaction variable between Mutuelles enrollment and poverty indicators. To control for district-level heterogeneity of health systems at the individual-level analysis, we used the reported number of community health workers in the DHSST to construct a variable indicating the number of district-level community health workers per capita in the 27 rural districts. As described in online Panel B, community health workers play a key role in implementing nutritional services in rural areas. We used a multilevel logistic regression model in individual-level analysis to estimate the association between children’s Mutuelles status and their likelihood of being stunted, controlling for child’s gender, poverty status, maternal factors (age, education, and height), and district- and village-level clustering effects. The multilevel model enabled us to control for clustering effects and correct for standard errors in the higher levels so as to obtain coefficients with more efficiency.40 Because of the process by which the RDHS sampled villages at the first stage, children from the same village could be more similar to each other than to children from other villages. As a result of decentralized health care delivery at the district level in Rwanda, children in the same district were more likely to share similarities to each other than to children of other districts.41 Ignoring these clustering effects could lead to underestimation of standard errors and assumption of statistical significance where it does not exist.42 To control for clustering effects, we included random effects at the village and the district level. More details on modeling are presented in the online Supplementary Methods (available as a supplement to the online version of this article at http://www.ajph.org). We applied a standard 2-step procedure recommended for discrete-outcomes multilevel models and started with a first-order marginal quasi-likelihood approach to get crude estimates, followed by a Bayesian Markov chain Monte Carlo approach in the second stage to improve the approximations.40,43 Details can be found in the online Supplementary Methods (see also Figures A and B, available as a supplement to the online version of this article at http://www.ajph.org). We report odd ratios and credible intervals. The credible intervals derived from the Markov chain Monte Carlo method indicate that, with 100 000 simulations, the true estimate will lie in the credible interval with a probability of 95% (Table A, available as a supplement to the online version of this article at http://www.ajph.org).43,44 We used the Stata version 14 command runmlwin (StataCorp LP, College Station, TX) to perform the multilevel statistical analysis.45–47 Because the inclusion of different covariates and the use of different estimation methods may alter the association between the exposure and the outcome variable, we examined the sensitivity of the results to model specifications. We estimated the main association between Mutuelles status and stunting (with or without covariates) using the logistic model with multilevel random effects (model 1 to model 3 in Table B, available as a supplement to the online version of this article at http://www.ajph.org). We then estimated the association between Mutuelles status and stunting among children in households below versus above the poverty line by adding the interaction variable between Mutuelles and poverty indicators (model 4 in Table B). We also tested the sensitivity of results to estimation methods by conducting an analysis with a fixed-effects model including the 27 district indicators in regression analysis. A fixed-effects model generates less biased estimates, whereas a random-effects model generates more efficient estimates.48
N/A