Geographic weighted regression analysis of hot spots of anemia and its associated factors among children aged 6–59 months in Ethiopia: A geographic weighted regression analysis and multilevel robust Poisson regression analysis

Introduction Anemia among children aged 6–59 months remains a major public health problem in low-and high-income countries including Ethiopia. Anemia is associated with significant consequences on the health of children such as under-five morbidity and mortality, increased risk of infection, and poor academic performance. The prevalence of anemia in Ethiopia has varied across areas. Therefore, this study aimed to investigate the geographic weighted regression analysis of anemia and its associated factors among children aged 6–59 months in Ethiopia. Methods This study was based on the 2016 Ethiopian Demographic and Health Survey (EDHS) data. A total weighted sample of 8482 children aged 6–59 months was included. For the spatial analysis, Arc-GIS version 10.7 and SaTScan version 9.6 statistical software were used. Spatial regression was done to identify factors associated with the hotspots of anemia and model comparison was based on adjusted R2 and Corrected Akaike Information Criteria (AICc). For the associated factors, the multilevel robust Poisson regression was fitted since the prevalence of anemia was greater than 10%. Variables with a p-value < 0.2 in the bi-variable analysis were considered for the multivariable analysis. In the multivariable multilevel robust Poisson regression analysis, the adjusted prevalence ratio with the 95% confidence interval was reported to declare the statistical significance and strength of association. Results The prevalence of anemia among children aged 6–59 months was 57.56% (95%CI: 56.50%, 58.61%) with significant spatial variation across regions in Ethiopia. The significant hot spot areas of anemia among children aged 6–59 months were detected in the central, west, and east Afar, Somali, Dire Dawa, Harari, and northwest Gambella regions. Mothers who had anemia, a child aged 23–59 months, mothers aged 15–19 years, and coming from a household with a poorer or poorest household were significant predictors of the spatial variations of anemia among children aged 6–59 months. In the multilevel robust Poisson analysis, born to mothers aged 30–39 (APR = 0.84, 95% CI: 0.76, 0.92) and 40–49 years (APR = 0.73, 95% CI: 0.65, 0.83), mothers who didn’t have formal education (APR = 1.10, 95% CI: 1.00, 1.20), Children in the poorest household wealth index (APR = 1.17, 95% CI: 1.06, 1.29), being 4–6 (APR = 1.08, 95% CI: 1.02, 1.13) and above 6 order of birth (APR = 1.15, 95% CI: 1.07, 1.23), children born to anemic mothers (APR = 1.24, 95% CI: 1.19, 1.29), children aged 24–59 months (APR = 0.70, 95% CI: 0.68, 0.73), stunted children (APR = 1.09, 95% CI: 1.04, 1.13) and underweight children (APR = 1.07, 95% CI: 1.03, 1.13) were significantly associated with anemia among children aged 6–59 months. Conclusion and recommendation Anemia is still a public health problem for children in Ethiopia. Residing in a geographic area where a high proportion of children born to mothers aged 15–19 years, a child aged 6–23 months, coming from a household with poorer or poorest wealth index, and mothers with anemia increased the risk of experiencing anemia among children aged 6–59 months. Maternal education, maternal age, child age, household wealth, stunting, underweight, birth order, and maternal anemia were significant predictors of anemia among children. The detailed map of anemia hot spots among children aged 6–59 months and its predictors could assist program planners and decision-makers to design targeted public health interventions.


Conclusion and recommendation
Anemia is still a public health problem for children in Ethiopia. Residing in a geographic area where a high proportion of children born to mothers aged 15-19 years, a child aged 6-23 months, coming from a household with poorer or poorest wealth index, and mothers with anemia increased the risk of experiencing anemia among children aged 6-59 months. Maternal education, maternal age, child age, household wealth, stunting, underweight, birth order, and maternal anemia were significant predictors of anemia among children. The detailed map of anemia hot spots among children aged 6-59 months and its predictors could assist program planners and decision-makers to design targeted public health interventions.

Background
Anemia is the commonest public health problem in Low-and Middle-Income Countries (LMICs) [1][2][3]. Globally, an estimated 1.62 billion people are anemic, of these more than 43% occurred in LMICs particularly in Asia and Africa [4,5]. Under-five children are highly affected by anemia compared to the general population [6,7]. According to World Health Organization (WHO), an estimated 293 million under-five children are anemic globally with a prevalence of 47.4% [8]. In Ethiopia, the prevalence of anemia among under-five children was 57% [7,9]. Anemia during childhood has been linked to developmental delay, recurrent infections, reduced working capacity, poor school performance, and dilated cardiomyopathy [10,11].
Previous studies on anemia among under-five children revealed that residence, maternal educational status, taking drugs for intestinal parasites, child nutrition status (stunting, wasting, and underweight), maternal age, child age, household wealth index, maternal anemia status, child-size at birth, birth order, parity, type of water source, type of toilet facility, type of birth, sex of children, and media exposure were significantly associated with anemia [24][25][26][27].
According to studies reported on the prevalence and associated factors of anemia, the prevalence of anemia has significant variation across regions in Ethiopia [7,[28][29][30][31][32][33][34]. Therefore, this study aimed to investigate spatial regression analysis of anemia and its associated factors among children aged 6-59 months. The findings of this study will help policymakers, program planners, and other health care programs in guiding health programs and prioritize prevention and intervention programs. Besides, mapping hotspot areas of anemia will provide a deeper understanding of the impacts of already implemented interventions in each region of the country.

Study area, data source, and study period
This study was based on the 2016 Ethiopian Demographic and Health Survey (EDHS) data. The 2016 EDHS was the fourth DHS in Ethiopia, which was conducted every 5 years. The EDHS is mainly aimed to generate updated health and health-related indicators such as maternal mortality, child mortality, family planning, vaccination, and maternal health care service utilization. Ethiopia is administratively divided into nine geographical regions (Tigray, Afar, Amhara, Oromia, Somalia, Benishangul-Gumuz, Southern Nation Nationality, and People's Region (SNNPR), Gambella, and Harari) and two self-administrative cities (Addis Ababa and Dire Dawa). Each region is subdivided into zones, each zone into Weredas, and each Wereda into Kebeles (which is the lowest administrative unit). A multistage stratified sampling technique was applied to select the study designs. In the first Enumeration Areas (EAs) was randomly selected and in the second on average 28 households per clusters/EAs.

Sample and population
Hemoglobin testing was carried out among children aged 6-59 months in the selected households using HemoCue rapid testing methodology. For the test, a drop of capillary blood was taken from a child's fingertip or heel and was drawn into the micro cuvette which was then analyzed using the photometer that displays the hemoglobin concentration. Then, anemic status was determined based on the hemoglobin level. The EDHS has several datasets such as men (MR file), women (IR file), children (KR file), birth (BR file), and household (HR file) datasets. For this study, we used the Kids Record dataset (KR file), and a total weighted sample of 8482 children aged 6-59 months was included.

Study variables
The dependent variable was the anemia status of children aged 6-59 months, which was categorized into anemic (hemoglobin level < 11 g/dl) and not anemic (hemoglobin �11.0 g/dl). It was assessed based on the hemoglobin concentration in blood adjusted for altitude.
In EDHS, before determining a child is anemic or not, hemoglobin adjustment for altitude was done by subtracting or adding the adjusted Hgb value to each individual observed Hgb value.
The Hgb adjustment was made using the formula; The adjustment for altitude was done to take into account the reduction in oxygen saturation of the blood. The independent variables considered in this study were region, residence, maternal age, maternal educational status, household wealth status, media exposure, maternal anemia status, sex of children, type of birth, age of children, size of child at birth, water source, type of toilet facility, parity, birth order, taking drugs for intestinal parasites in the last 6 months, wasting status (Z-scores for Weight-for-Height (WHZ)), underweight status (Zscores for Weight-for-Age (WAZ)) and stunting status (Z-scores for Height-for-Age (HAZ)).

Data management and analysis
Factors associated with anemia. Data extraction, coding, and analysis were done using Stata version 14 and Arc-GIS version 10.6 statistical software. The weighted data were used for analysis to restore the representativeness of the data. Since the EDHS data has a hierarchical nature, the Intra-class Correlation Coefficient (ICC) was estimated to assess the clustering effect. The ICC indicated that there was a significant clustering effect (ICC>10%). This study was a cross-sectional study and the prevalence of anemia was greater than 10%, and if we reported the odds ratio it could overestimate the association between anemia and the independent variables. In such cases, the prevalence ratio is the best measure of association, and therefore, multilevel Poisson regression analysis with robust variance was fitted to identify predictors of anemia. Variables with a p-value<0.2 in the bi-variable multilevel Poisson regression analysis were considered for the multivariable analysis. Deviance was used to verify model fitness, and a model with the lowest deviance was considered the best-fit model. Finally, the Adjusted Prevalence Ratio (APR) with its 95% confidence interval (CI) was reported, and variables with p value<0.05 in the multivariable analysis were considered as significant predictors of anemia among under-five children.
Spatial analysis. The global spatial autocorrelation (Global Moran's I) was done to assess whether the spatial distribution of anemia among under-five children in Ethiopia was dispersed, clustered, or randomly distributed [35]. Global Moran's I is a spatial statistic used to measure spatial autocorrelation by taking the entire data set and produce a single output value that ranges from -1 to +1. Moran's, I value close to −1 indicates that anemia among under-five children is dispersed, whereas Moran's I close to +1 indicates anemia among under-five children is clustered and if Moran's I close to 0 revealed that anemia among under-five children is randomly distributed. A statistically significant Moran's I (p < 0.05) value showed that anemia among under-five children is non-random. The hotspot analysis was done using the Getis-OrdGi � statistics to explore how spatial autocorrelation varies over the study location by calculating GI � statistic for each area. Z-score is computed to determine the statistical significance of clustering, and the p-value is computed for the significance. Statistical output with high GI � indicates "hotspot" whereas low GI � means a "cold spot" [36].
Spatial regression analysis. The Ordinary Least Square (OLS) regression and Geographic Weighted Regression (GWR) statistical analysis were employed for exploring the spatial relationship between anemia among under-five children and the explanatory variables. The outcome variable for spatial regression analysis was the percentage of anemia among under-five children at the EA level. A neighborhood or bandwidth is the distance band or the number of neighbors used for each regression equation, it is the most important parameter for spatial regression as it controls the degree of smoothening in the model. The complexity of spatial regression model depends not only by the number of variables in the model but also the bandwidth. There are three choice of band width methods such as AICc, CV and bandwidth parameter. For this study we have used adaptive kernel whose bandwidth was found by minimizing the AICc value.
Ordinary Least Squares (OLS) regression. The spatial regression modeling was performed to identify predictors of the spatial heterogeneity of anemia among under-five children. OLS is a global statistical model for testing and explaining the relationship between the dependent and independent variables [37]. It uses a single equation to estimate the relationship between the dependent and independent variables and assumes stationarity or consistent relationship across the study area. The OLS was used as a diagnostic tool and for selecting the appropriate predictors (concerning their relationship with anemia) for the Geographic Weighted Regression (GWR) model [38].
The OLS can automatically check the multicollinearity between independent variables (redundancy among explanatory variables). The multicollinearity was assessed using the Variance Inflation Factor (VIF). If the VIF values are greater than 10 in the OLS model, it indicates the existence of multicollinearity among the explanatory variables and should apply to leave one out an approach based on the VIF values. Besides, the autocorrelation statistic was applied to detect whether there is spatial autocorrelation or clustering of the residuals which violates the assumptions of OLS. The spatial independence of the residuals was assessed with the global spatial autocorrelation coefficient Moran's I value. The Moran's I value ranges from +1 (positive autocorrelation) and -1 (negative autocorrelation).
Geographically Weighted Regression (GWR). A local spatial statistical technique that assumes the non-stationarity in relationships/ heterogeneity in the relationship between the dependent and explanatory variables across EAs [38][39][40]. The GWR analysis is considered when the Koenker statistics is significant (p-value<0.05), which means the relationships between the dependent and the independent variable change from location to location. In the GWR analysis, the coefficients of the explanatory variables take different values across the study area. Mapping the GWR coefficients associated with the explanatory variables, which are produced using the GWR, provides insight for targeted interventions. The corrected Akaike Information Criteria (AICc) and adjusted R-squared for model comparison of OLS (global model) and GWR (local) model. A model with the lowest AICc value and a higher adjusted Rsquared value was considered as the best-fitted model for the data.

Descriptive results
A total weighted sample of 8482 children aged 6 to 59 months was included. Of these, more than half (51.85%) of the children were males. The majority (43.87%) of the children were in the Oromia region and 7621 (89.85%) were from the rural areas. Nearly half of the mothers (48.76%) were aged 20-29 years, and 5684 (67.02%) of the mothers didn't attain formal education. About 1988 (23.44%) and 1142 (13.46%) of the mothers fall within the poorest and richest household index quintiles, respectively. Nearly one-third (30.1%) of the children's mothers were anemic. Regarding children's nutritional status, about 40.73%, 25.29%, and 9.38% of the children were stunted, underweight and wasted, respectively ( Table 1).

Prevalence of anemia among under-five children in Ethiopia
The prevalence of anemia among under-five children was 57.56% (95%CI: 56.50%, 58.61%). The highest prevalence of anemia among under-five children was observed in Somali

Factors associated with anemia among under-five children
Random effect analysis results. In the null model, the ICC value was 18.8% (95% CI: 15.99%, 21.98%), indicated that about 18.8% of the overall variability in anemia was explained by the between cluster variation while the remaining 81.20% was attributed to the individuallevel variation. Besides, the Likelihood Ratio (LR) test was (LR test vs. logistic model: X 2 (01) = 512.36, p< 0.0001), which showed that the mixed-effect models were the best-fitted model for this data compared to the standard model.
Fixed effect analysis results. In the multivariable mixed-effect Poisson regression with a robust variance; maternal age, maternal education, household wealth index, birth order, maternal anemia, age of children, stunting, and underweight were significantly associated with anemia among under-five children. The prevalence of anemia among children born to mothers aged 30-39 and 40-49 years was decreased by 16% (APR = 0.84, 95% CI: 0.76, 0.92) and 27% (APR = 0.73, 95% CI: 0.65, 0.83) compared to children born to mothers aged less than 20

Spatial distribution of anemia among under-five children.
The highest prevalence of anemia among children aged 6-59 months was observed in Afar, Somali, Tigray, Dire Dawa, and east Amhara regions (Fig 2). The spatial distribution of anemia among children aged 6-59 months showed significant spatial variation across the country with a global Moran's I value of 0.089 (p-value<0.01) (Fig 3). The statistically significant hotspot areas of anemia were identified in the central, west, and east Afar, Somali, Dire Dawa, Harari, and northwest Gambella regions. While significant cold spot areas were detected in the northwest SNNPRs, Addis Ababa, Benishangul-gumuz, central and southwest Amhara regions (Fig 4).
The global ordinary least square regression analysis results. The OLS model was calibrated to diagnose multicollinearity among the independent variables and the mean VIF was less than 10. In the OLS analysis, the model explained about 26% (adjusted R 2 = 0.22) of the variation in anemia among children aged 6-59 months with AICc = -169.06. The Joint F-statistics and Wald statistics were significant (p<0.05), which proves that the model was statistically significant. The spatial distribution of residuals was normally distributed as the Jarque-Bera statistics were non-significant (the residuals were normally distributed) (p = 0.18). The Koenker statistics were statistically significant, indicates that the relationship between the independent variables and the dependent variable was non-stationary or heterogeneous across the study areas. This indicates that GWR should be applied (since the Koenker statistics showed the non-stationarity in the relationship) as it assumes the spatial heterogeneity of the relationship between independent and dependent variables across space. The proportion of women aged 15-19 years, the proportion of women who were anemic, the proportion of women who attained formal education, and the proportion of children aged 6-23 months were significantly associated with the percentage of anemia among children aged 6-59 months in the OLS model ( Table 3).
The geographically weighted regression analysis result. The GWR analysis showed that there was a significant improvement over the global model (OLS). The AICc value decreased from -169.06 (for the OLS model) to -257.41 (for the GWR model). The difference was 88.35 implied that the GWR best explains the spatial heterogeneity of anemia among children aged 6-59 months. Besides, the adjusted R 2 was 0.38, the model's ability to explain anemia among children aged 6-59 months has been improved by using GWR since the adjusted R 2 was 0.38, indicates that GWR improved the model explaining the power of the OLS model by about 12% (Table 4). In the geographically weighted regression analysis, the proportion of women who had formal education, the proportion of women aged 15-19 years, the proportion of women with poverty, the proportion of being birth order above 6 th , proportion of women who had anemia, proportion of children aged 6-23 months, the proportion of stunted children and proportion of wasted children were considered as explanatory variables in the GWR model since it was significant in the multilevel robust Poisson regression analysis as well as have good R 2 in the exploratory analysis. The proportion of women who had anemia had a positive relationship with the proportion of anemia among children aged 6-59 months. As the proportion of women who had anemia increased, the percentage of anemia among under-five children increased in the entire Amhara, Benishangul-gumuz, Gambella, and SNNP regions. The geographic area with redcolored points indicates the highest coefficient of the proportion of maternal anemia (Fig 5). The proportion of mothers aged 15-49 years was significantly associated with the increased risk of anemia among children aged 6-59 months, with the highest effect of mothers age observed in southeast Amhara, east Tire, west Afar, Gambella, and southwest SNNP regions (Fig 6). The proportion of mothers in the poorest household wealth status showed strong and positively associated with increased risk of anemia among under-five children in Somali regions (Fig 7). The proportion of children aged 6-23 months had a significant positive association with anemia among children aged 6-59 months (Fig 8).

PLOS ONE
Geographic weighted regression analysis of hot spots of anemia and its associated factors

Discussion
Anemia among children aged 6-59 months remains a common public health problem in Ethiopia. In this study, the prevalence of anemia among children aged 6-59 months in Ethiopia was 57.56% (95%CI: 56.50%, 58.61%) ranged from 42.66% in the Amhara region to 83.24% in the Somali region. This was higher than a study reported in Ghana [41] and China [42]. Even though the combined strategies particularly iron supplementation and infectious disease management (such as malaria and helminths infections) are being introduced by the WHO to combat anemia, anemia remains a serious health care problem in Ethiopia. Besides, the spatial distribution of anemia among children aged 6-59 months in Ethiopia was non-random and the hotspot areas of anemia were identified in the central, west, and east Afar, Somali, Dire-Dawa, Harari, and northwest Gambella regions. Though numerous interventions have been implemented to reduce anemia among children such as nutritional interventions like food fortification and supplementation [43][44][45], iron-folate supplementation during pregnancy [46,47], appropriate breastfeeding practice [48], and improving personal hygiene [49], there was high spatial heterogeneity of anemia across areas. The potential reason may be due to the long-standing prevalence of severe malnutrition among under-five children, because of insufficient dietary intake of nutrients in Ethiopia [50,51]. Besides, Ethiopian children's are highly affected by infectious diseases such as malaria, hookworms, Schistosoma, and visceral leishmaniasis, due to their frequent exposure to poor sanitation and environmental conditions that favor the transmission and spread of parasites [52][53][54]. In the spatial regression analysis, maternal age, child age, poverty, and maternal anemia was significant predictors of hotspot areas of anemia among under-five children. There is a positive relationship between poverty and hotspots of anemia in Somali regions. Inadequate wealth in the community is associated with an increased risk of intestinal infections such as hookworm, amoebiasis, and ascariasis and these could increase the risk of anemia [55,56]. Besides, children in the poorest household are prone to undernutrition like lack of folate, iron vitamin B12, and vitamin A [57]. An increased proportion of mothers aged 15-19 years increases the odds of anemia among under-five children in southeast Amhara, east Tigray, west Afar, Gambella, and southwest SNNP regions. The possible explanation might be due to adolescent pregnant mothers are more likely to give low birth weight, preterm or small for gestational age babies [58,59], in turn, these might contribute to low hemoglobin levels in the blood. The proportion of children aged 6-23 months had a significant positive association with anemia among children aged 6-59 months, this could be due to children aged 6-23 months are at higher risk of exposure to infectious diseases like foreign body aspiration, pneumonia, diarrheal diseases, and intestinal infections because of exposure to complementary feeding, which increases the malabsorption of iron or folate and increases the risk of anemia [60]. In addition, an increased proportion of anemic mothers increases the odds of anemia among children. This could be due to both mothers and children mostly share a common home environment, which involves mutual exposure to a common set of physical, socioeconomic, and dietary conditions [61]. Also, maternal anemia might be associated with poor birth outcomes such as low birth weight and prematurity of the child, which might lead to limited fetal iron stores and the amount of iron secreted by the breast milk might be insufficient for the daily iron requirement of the child [62].
In the multivariable mixed-effect Poisson regression with a robust variance; maternal age, maternal education, household wealth index, birth order, maternal anemia, age of children, stunting, and underweight were significantly associated with anemia among under-five children.
A child born to a mother aged 30-39 years and 40-49 years had decreased prevalence of anemia compared to a child born to a mother aged less than 20 years. These findings were in line with previous studies reported in Ethiopia [63] and Bangladesh [64]. Anemia is a major public health problem worldwide primarily affects women particularly teenage pregnant girls. Adolescents are particularly susceptible because of their rapid growth and associated high iron requirements. So, pregnancy and lactation increase their nutritional demand, in turn, babies born to teenagers are at higher risk of anemia [65]. Children aged 6-23 months are at higher risk of anemia compared to children aged 24-59 months. This is consistent with study findings reported in Ethiopia [66], and Uganda [67]. It could be due to complementary feeding is initiated after 6 months of birth and during this, a child is exposed to contaminated food and malabsorption syndrome, this could increase nutritional deficiency anemia. In addition, an increased iron requirement due to rapid growth, low availability of foods rich in iron, and lack of diet variety. Iron intake is also likely to improve with age as a result of a more varied diet, including the introduction of meat and other iron-containing foods. The prevalence of anemia was higher among children born to anemic mothers compared to children born to non-anemic mothers. It is in line with study findings in Ghana [68], and Bangladesh [64]. This might be due to children born to anemic mothers are more likely to have nutritional deficiencies like folate, iron, vitamin B12, and vitamin A. Besides, anemic mothers might have underlying diseases such as malaria, HIV/AIDS, or other genetic diseases, and these could transmit to the newborn transplacental, which in turn, could result in anemia. Moreover, mothers and children have common nutritional sources and practices. Another most important significant predictors of anemia among children were maternal education. A child born to a mother who did not have formal education had a higher prevalence of anemia compared to a child born to a mother who attained secondary education or above. It is consistent with study findings reported in Korea [69]. The possible explanation might be due to children born to mothers who did not attain formal education are more likely to consume iron-rich food like meat and poultry compared to children of educated mothers [70]. In addition, educated mothers are more likely to utilize child health services, which can have a positive effect on their children's health outcomes, and improved mothers' level of education results in the corresponding improvement in child feeding practice [71]. Educated mothers exclusively breastfeed their children and initiate appropriate complementary feeding after six months of gestation [72].
High birth order was significantly associated with an increased prevalence of anemia among children aged 6-59 months. This is consistent with studies reported in Sub-Saharan Africa [73] and India [62]. This might be due to as the birth order increases there is maternal nutritional depletion and higher-order births are more likely to be low birth weight, this could increase the risk of anemia [74].
In this study, child nutritional status was significantly associated with anemia among children aged 6-59 months. Children who were stunted and/or underweight were more likely to be anemic than their counterparts. This could be due to anemia and malnutrition often share common causes, it is expected that multiple nutrition problems would co-occur in the same individuals. Low intake of iron-rich foods and diminished nutrient absorption caused by changes in the gastrointestinal epithelium in malnourished individuals contribute towards the development of anemia.

Strength and limitations of the study
This study has several strengths. It is based on the nationally representative DHS data which was weighted and a multilevel model was fitted that enables us to generalize these findings at the national level. Besides, spatial and geographic weighted regression analyses were conducted, the findings can assist the policymakers and program planners to design spatially targeted public health interventions to reduce the incidence of anemia. The findings of this study should be interpreted considering the following limitations. First, important variables such as underlying medical conditions like malaria, HIV/AIDS, visceral leishmaniasis, and hookworm infections, etc were not considered in the analysis as these variables were not collected in EDHS 2016. Also, we are unable to show the cause-effect relationship between the dependent and independent variables as the DHS data has cross-sectional nature. The geographic locations (GPs) of enumeration areas were displaced up to 2 kilometers in urban and 5 kilometers for most enumeration areas in rural and 10 kilometers for 1% of clusters in rural areas for the sake of privacy, could affect the estimated cluster effects in the spatial regression.

Conclusion
The prevalence of anemia among children aged 6-59 months in Ethiopia was high, and there was a significant spatial variation of anemia among children aged 6-59 months across regions in Ethiopia. Significant hotspot areas of anemia were identified in the central, west, and east Afar, Somali, Dire Dawa, Harari, and northwest Gambella regions. Being a resident in a geographic area with high community poverty, a high proportion of maternal anemia, a high proportion of children aged 6-23 months, and mothers aged 15-19 years increased the risk of experiencing anemia among children aged 6-59 months. Maternal age, maternal education, birth order, maternal anemia, child age, household wealth index, stunting, and underweight were significant predictors of anemia. Therefore, public health interventions targeting hotspot areas of anemia through empowering women with education help reduce the incidence of anemia among children aged 6-59 months.