Associations of childhood, maternal and household dietary patterns with childhood stunting in Ethiopia: proposing an alternative and plausible dietary analysis method to dietary diversity scores

Background Identifying dietary patterns that consider the overall eating habits, rather than focusing on individual foods or simple counts of consumed foods, better helps to understand the combined effects of dietary components. Therefore, this study aimed to use dietary patterns, as an alternative method to dietary diversity scores (DDSs), and investigate their associations with childhood stunting in Ethiopia. Methods Mothers and their children aged under 5 years (n = 3788) were recruited using a two-stage random cluster sampling technique in two regions of Ethiopia. Socio-demographic, dietary and anthropometric data were collected. Dietary intake was assessed using standardized dietary diversity tools. Household, maternal and child DDSs were calculated and dietary patterns were identified by tetrachoric (factor) analysis. Multilevel linear and Poisson regression analyses were applied to assess the association of DDSs and dietary patterns with height-for-age z score (HAZ) and stunting, respectively. Results The overall prevalence of stunting among children under-five was 38.5% (n = 1459). We identified three dietary patterns each, for households (“fish, meat and miscellaneous”, “egg, meat, poultry and legume” and “dairy, vegetable and fruit”), mothers (“plant-based”, “egg, meat, poultry and legume” and “dairy, vegetable and fruit” and children (“grain based”, “egg, meat, poultry and legume” and “dairy, vegetable and fruit”). Children in the third tertile of the household “dairy, vegetable and fruit” pattern had a 0.16 (β = 0.16; 95% CI: 0.02, 0.30) increase in HAZ compared to those in the first tertile. A 0.22 (β = 0.22; 95% CI: 0.06, 0.39) and 0.19 (β = 0.19; 0.04, 0.33) increase in HAZ was found for those in the third tertiles of “dairy, vegetable and fruit” patterns of children 24–59 months and 6–59 months, respectively. Those children in the second (β = −0.17; 95% CI: -0.31, −0.04) and third (β = −0.16; 95% CI: -0.30, −0.02) tertiles of maternal “egg, meat, poultry and legume” pattern had a significantly lower HAZ compared to those in the first tertile. No significant associations between the household and child “egg, meat, poultry and legume” dietary patterns with HAZ and stunting were found. Statistically non-significant associations were found between household, maternal and child DDSs, and HAZ and stunting. Conclusion A higher adherence to a “dairy, vegetable and fruit” dietary pattern is associated with increased HAZ and reduced risk of stunting. Dietary pattern analysis methods, using routinely collected dietary data, can be an alternative approach to DDSs in low resource settings, to measure dietary quality and in determining associations of overall dietary intake with stunting. Electronic supplementary material The online version of this article (10.1186/s12937-018-0316-3) contains supplementary material, which is available to authorized users.

(Continued from previous page) Conclusion: A higher adherence to a "dairy, vegetable and fruit" dietary pattern is associated with increased HAZ and reduced risk of stunting. Dietary pattern analysis methods, using routinely collected dietary data, can be an alternative approach to DDSs in low resource settings, to measure dietary quality and in determining associations of overall dietary intake with stunting.
Keywords: Dietary data analysis, Dietary diversity score, Dietary pattern, Stunting

Background
Although the global burden of stunting decreased between 1990 and 2015 by more than 25%, it has continued to be a major nutrition-related risk factor causing 257 deaths per 100,000 globally [1]. In sub-Saharan Africa, 11.8 disabilityadjusted life years (DALYs) and 136,455 childhood deaths in 2015 were attributable to stunting [2]. In Ethiopia, 38% of children under-five years of age were stunted [3] and it was a risk factor for 960,742 DALYs and 11,065 deaths in 2015 [2]. Stunting also halts the development of societies by negatively affecting mental and physical health [4]. Suboptimal nutrition is a major contributor to stunting in developing countries [5], although there are other causal and contextual factors [4].
Although there are well-established methods, collation and analysis of dietary data have remained challenging in low-income countries (LICs) for various reasons, including high costs, lack of centralized platforms for dietary data, little investment in research, low capacity and technical complexity [6]. As a result, dietary assessment is mainly dependent on approaches which require low cost and provide low quality. Dietary diversity assessment has remained the most commonly used method of data collection, analysis and interpretation approach in LICs. Dietary diversity scores (DDSs) of households, women and children [7][8][9] are important tools and the most commonly used indicators of assessing the adequacy of nutrient intake. In many studies, it has been also demonstrated that DDSs were useful indicators of micronutrient status [10][11][12][13] and a higher DDS is associated with a lower risk of stunting [14][15][16]. The indicators are relatively simple and suitable for use in large surveys [3]. However, data collected for the purpose of DDS analysis are qualitative, and in most cases, they are dichotomized (yes/no) [8] restricting further analyses. Thus, the analyses depend on a simple count of food groups and do not consider the correlations of the food groups and their impact on nutritional (disease) outcomes. In addition, because the main purpose of DDS analysis is on the number, rather than the type of foods consumed, this may ignore the antagonist, additive and synergistic effect of food groups.
Currently the focus of nutritional epidemiology is to investigate the patterns of multiple food and nutrient intakes without ignoring the interactions. Methodological development over the last two decades enables us to explore the association between diet and disease outcomes through a systematic consideration of the correlation between the components of the overall diet, that is dietary patterns [17,18]. A study by Humphries et al. reported that total food expenditures (using food groups for child DDSs) did not significantly predict HAZ in Ethiopia. In this study, household food group expenditure index, determined by factor analysis of disaggregated food expenditure, was found to be a significant predictor for HAZ [19]. This leads to a premise that a mere aggregate availability and accessibility of the included food groups are not the determining factors for HAZ, rather the specific types of food groups available and their consumption pattern. Another study in the same cohort strengthens this conclusion [20]. Therefore, identifying dietary patterns that consider the overall eating habits, rather than focusing on individual foods (simple counting of consumed foods), better reflect the complexity of dietary intakes and help to understand the combined effect of diet components [21]. In this study, we aimed to identify household, maternal and child dietary patterns and investigate their associations with childhood stunting in Ethiopia using the same dietary data collected for determining DDS. In addition, the study compares the findings with the estimates of associations between DDSs and stunting. To the best of our knowledge, this is the first study investigating the aforementioned objectives.

Study area and participants
A cross-sectional study was conducted in the South Nations, Nationalities and People (SNNP) and Tigray (northern Ethiopia) regions between June and September 2014. The two regions are geographically located at opposite ends of Ethiopia, in the south and north, with differences in agroecology, subsistence farming being the most common occupation in both regions. The SNNP is a larger geographic area and has a greater population size compared to the Tigray region. This study was part of a larger project of the Alive and Thrive's (A&T) impact evaluation for community-based interventions. The major objectives of the evaluation included assessment of infant and young child feeding (IYCF) practices and stunting prevalence. The baseline and progress evaluation of the project were conducted between June and September 2010 and 2013, respectively.

Sample size calculation and sampling technique
The sample size was calculated based on the 2010 baseline and 2013 progress evaluation surveys' estimates. It took into account an intracluster correlation of 0.03-0.04. A total of 75 clusters (enumeration areas [EAs]), a one-sided test, a power of 80%, and a significance level of α = 0.05 were included in the calculations. Based on these considerations, a minimum sample size of 2950 was required (Table 1). However, in total, data were collected from 3788 children and their mothers.
A two-stage cluster sampling technique was used to select households with children under five. In the first stage (primary sampling unit), EAs were selected from 89 districts (the second smallest administrative units in Ethiopia). The EA is a geographical unit devised by the Central Statistical Authority (CSA) of Ethiopia, which consists of 150-200 households. This is the smallest cluster used in Demographic and Household Surveys and roughly coincides with the kebele (the smallest administrative units) boundaries (Fig. 1). A total of 75 EAs (26 from Tigray and 49 from SNNP), from 56 districts (19 from Tigray and 37 from SNNP) were selected using probability proportional to size (PPS) sampling in relation to the population of the EAs.
In the second stage, children between 0 to 59.9 months (n = 3788) were selected. A complete household listing with the number of children residing in each household, in each selected cluster, was developed in collaboration with the local health and administrative offices. This list included identification of all eligible candidates for the survey (mothers of those children under 60 months of age). From this list, three sampling frames were developed: children aged 0-5.9 months, 6-23.9 months, and 24-59.9 months. From each sampling frame, study subjects were selected using a systematic random sampling (SRS). Households selected to participate in one age category were not included in the other sampling frames, even if there were other eligible children in a household.

Data collection
All data (interview and anthropometric) were collected by trained data collectors who had bachelor degrees or above. To maintain the data quality, demonstrations and pilot testing were conducted during the training period. Trained supervisors oversaw and monitored the data collectors during the field work. The supervisors also checked 5-10% of the anthropometric and interview data to ensure reliability.

Anthropometry data collection
The length or height of children was measured to the nearest 0.1 cm using the United Nations International Children's Emergency Fund (UNICEF) recommended wooden board with an upright wooden base and movable headpieces. Children ≥24 months were measured while standing upright while those less than 24 months in the recumbent position. Weight was measured to the nearest 0.1 kg using UNICEF's scale [22]. Immunization cards or home records of the date of birth, if available, were used to determine the age of the children. In the case of absence of these documents, mother's recall was taken using the local calendar and then converted to the Gregorian calendar. Adult weighing and height scales were used to measure maternal weight and height, respectively. Mothers were asked to remove shoes and heavy cloths before weight and height were measured. Weight and height were recorded to the nearest 0.1 kg and 1.0 cm, respectively.

Dietary data collection
Dietary data for children [7] and women [8] and household dietary and food insecurity data [9,23] were assessed using standard tools. Dietary intake was assessed for the preceding day (24 h). For children aged 6-23 months, seven food groups (grains, roots and tubers; legumes and nuts; dairy products (milk, yogurt, cheese); flesh foods (meat, fish, poultry and liver/organ meats); eggs; vitamin-A rich fruits and vegetables; and other fruits and vegetables) were included [7]. For children 24-59 months, an additional two food groups (oils and fats and independent categories of other fruits and vegetables) were included. The maximum dietary diversity for women of reproductive age (MDD-W) assessment includes 10 food groups (grains, white roots and tubers, and plantains; pulses (beans, peas and lentils); nuts and seeds; dairy; meat, poultry and fish; eggs; dark green leafy vegetables; other vitamin A-rich fruits and vegetables; other vegetables; other fruits). Consumption of food by any of the household member from any of 12 food groups in the last 24 h was also assessed and the household DDS was determined. The food groups were cereals; roots and tubers; vegetables; fruits; meat, poultry, offal; eggs; fish and seafood; pulses/legumes/nuts; milk and milk products; oils/ fats; sugar/honey and miscellaneous [23].

Other covariates
Data including socio-demographic (such as maternal age, maternal education, sex of the head of the household and paternal education), economic (household asset), environmental factors (such as water source and latrine type), health service utilization (such as place of delivery for index child), and household characteristics (such as the number of under-five children living in a household) were collected. Different socio-economic indicators were combined using principal component analysis to construct household wealth. The factor scores were divided into quantiles (poorest, poorer, middle, richer and richest) to indicate the relative socio-economic status of the  participants. The highest level of education achieved was categorized into no education, primary, and secondary and above. Water source was classified as piped, other improved and unimproved. The type of functional latrine used in the household was categorized into traditional pit latrine, improved latrine and no facility/bush/field.

Anthropometry and dietary data analyses
Height-for-age z score (HAZ), an indicator of linear growth, was compared with reference data from the World Health Organization (WHO) Multicentre Growth Reference Study Group, 2006 [24] using the ENA (Emergency Nutrition Assessment) SMART (Standardized Monitoring and Assessment of Relief and Transitions) 2011 software. Children whose HAZ is <−2 SD from the median of the WHO reference population were considered stunted (short for their age). In our analysis, both HAZ (continuous) and stunting (categorical; stunted = 1/ not stunted = 0) were used as outcome variables. Maternal body mass index (BMI) was calculated based on the measured weight (kg) and height (meters) (weight [kg]/ (height[meters]) 2 ). For children aged 6 to 23, and 24 to 59 months, the minimum acceptable DDS was defined as consuming food from four or more of the standardized set of seven (6 to 23 months) or nine (24 to 59 months) food groups on the preceding day of the survey [7,25]. For women (mothers/ caregivers of the index children), a threshold of at least 5 food groups of the 10 was considered acceptable [8]. The scores of household DDS were continuous, ranging from 0 to 12, based on whether any of the members of the household consumed any of the 12 food groups. Minimum household DDS was not determined because a dichotomous indicator has not been developed [8,23]. However, we assumed that consumption of food groups above the median number as adequate. Nine Household Food Insecurity Access Scale (HFIAS) generic questions were used with a dichotomized response (0 = no/1 = yes) to assess food insecurity [9]. Each of the questions were asked with a recall period of four weeks (30 days). If a respondent answers "yes" to any of the above nine questions, frequency-of-occurrence questions were asked to determine whether the condition happened rarely (1 = once or twice), sometimes (2 = three to ten times) or often (3 = more than ten times). The insecurity status was categorized into four groups (secured, mild, moderate, and severe) using the Food and Nutrition Technical Assistance (FANTA) algorithm [9].
Dietary patterns were identified by polychoric (tetrachoric) analysis-a family of factor analysis which uses a tetrachoric correlation matrix to construct latent variables from dichotomized (binary) observed data [26]. For each of the dietary patterns, factor scores were assigned for all study participants. Factor scores show the relative position of the study participants in each of the identified patterns, thus reflecting adherence to the patterns. Pattern-specific factor scores are calculated as the sum of the products of the factor loading coefficients and standardized daily consumption of food and nutrient groups related to the pattern. The factor scores were orthogonally (varimax) rotated to create less correlation among the patterns and to facilitate their interpretability. Participants were then assigned into tertiles (first [lowest adherence]; second; and third [highest adherence] tertiles) based on their factor scores. Eigenvalues (>1), scree plots, and interpretability of the factors were used to determine the number of dietary and nutrient patterns. Factor loadings (the correlation between each pattern and the food and nutrient groups) were calculated. Percentages of variances (the variations that were explained by the identified dietary and nutrient patterns) were also computed.

Statistical analyses
The chi-square (categorical variables), ANOVA (normally distributed continuous variables) and Kruskal-Wallis (continuous but not normally distributed) tests were used to compare differences of proportions, means and medians, respectively, between groups. Principal component analysis (PCA) was used to compute economic status (in quintiles) of households.
To assess the associations of household, maternal and child dietary diversity and patterns with HAZ and childhood stunting, β coefficients and the prevalence ratio (PR) with their corresponding 95% confidence intervals (CIs) were determined using multilevel linear and Poisson regression models, respectively [27]. Since the data were collected using a multi-stage cluster sampling technique, stunting could potentially be correlated in clusters (EAs). We, therefore, used a two-level model with individual factors as level 1 and geographical areas (EAs) at level 2 (random effects). A stepwise backward elimination of covariates in the models was conducted and potential factors were retained at p-value < 0.20. This method was used for both individual and community level factors. Dietary diversity and pattern scores were treated as categorical (model 1) and continuous (model 2) variables. Estimates of associations were adjusted for sociodemographic factors (child age, sex, maternal age and education, number of under-five children in a household), Fig. 2 Household, maternal and child dietary patterns and corresponding factor loadings and proportion of food groups. The colour gradation reflects how large and in which direction was the correlation between the food groups and the dietary patterns. Deep green colour refers a relatively higher correlation (a higher intake) of the food groups with the corresponding patterns. Deep red refers relatively a lower correlation (a lower intake) of the food groups with the corresponding dietary patterns maternal anthropometry (height and BMI), infant and young child feeding practices (exclusive breastfeeding) and household food security at level 1. At level 2, water source was included. Model fit was assessed using Akaike's (AIC) and Bayesian (BIC) information criteria. We tested interactions between DDSs, dietary patterns, other covariates with HAZ and stunting using multiplicative terms. We conducted sensitivity analysis: 1) by labelling missing values of covariates as "missing" and including in the models; 2) by including and excluding covariates (such as household wealth, paternal education, place of delivery and latrine type). Further, the association between joint classifications of tertiles of dietary patterns and HAZ was explored. Statistical analyses were performed using Stata version 14.1 (Stata Corporation, College Station, TX, USA). A 2-sided ttest value of P < 0.05 was considered statistically significant.

Participant characteristics
Almost half of the children (1805, 47.7%) were female. Only 3.1% of the households had a female head. In 2353 (90.8%) of households, there was only one child aged under five. The median maternal age was 29.0 years (IQR = 25.0, 35.0). More than half (2069; 54.6%) of the mothers were illiterate. The mean maternal BMI was 20.2 kg/m 2 (SD = 2.4). Almost two-thirds (2371; 62.6%) of mothers delivered the index child at home. The prevalence of stunting among children aged 0-59 months was 38.5% with a mean (SD) HAZ of 1.6 (1.8). A fifth (777; 20.5%) of the study participants had missing values of household income (Table 2). Figure 2 depicts household, maternal and child dietary patterns and corresponding factor loadings. For each, we identified three dietary patterns. Pattern 2 ("egg, meat, poultry and legume") and pattern 3 ("dairy, vegetable and fruit pattern") were similar for all groups. The "egg, meat, poultry and legume pattern" was characterized by a high intake of eggs, meats, legumes, cereals, oils, fats and sweets. The "dairy, vegetable and fruit based" pattern was characterized by a high intake of fruits tubers, roots, vegetables and milk and milk products. A "plantbased pattern" and a "dairy, vegetable and fruit pattern"    were identified for mothers. While these two patterns appeared to be similar, the "plant-based" pattern was however characterized by a high intake of grains, tubers and leafy vegetables. Individual food items used for food groupings and the proportion of food groups consumed by households, mothers and children are depicted in Additional file 1: Tables S1, S2 and S3.

Dietary patterns
DDSs, dietary patterns, household food security and stunting Table 3 and Additional file 1: Table S4 show the DDSs of households, mothers and children by whether or not stunting was present. Three-quarters (2859; 75.5%) of the mothers had a DDS less than five. Only 13.1% of children aged 6-23 months had a DDS greater than or equal to four. The proportion of children aged 24-59 months who had a DDS less than or equal to three was 61.7% (1411). The proportion of food secure households was 46.0% (1744). A marginally significant statistical difference in the proportion of stunting was found between those children with a household DDS < =5 and >5 (p = 0.045). The prevalence of stunting was significantly different by maternal DDS (p = 0.005). No significant difference in the prevalence of stunting was found between the DDSs of children aged 6-23 and 24-59 months. Except for the "fish, meat and miscellaneous" dietary pattern (p = 0.097), there were significant differences in the prevalence of stunting across the tertiles of the other household dietary patterns. The prevalence of stunting (41.6%) in the first tertile of household "dairy, vegetable and fruit" dietary pattern was higher compared to the third tertile (33.3%) (p < 0.001). Children in the third tertile of maternal and child "dairy, vegetable and fruit" dietary pattern were less likely to be stunted compared to the first and second tertiles (p = 0.002) ( Table 4).

Associations of DDSs and dietary patterns with HAZ and stunting
After adjusting for potential individual and community level factors, no significant associations between household, maternal and child DDSs with stunting (HAZ) was found. Children in the third tertile of the household "dairy, vegetable and fruit" pattern had a 0.16 (β = 0.16; 95% CI: 0.02, 0.30) increase in HAZ compared to those in the first tertile. Similarly, the prevalence of stunting among children in the third tertile of the pattern was lower (PR = 0.83; 95% CI: 0.72-0.95) compared to those in the first tertile. Those children in the second (β = −0.17; 95% CI: -0.31, −0.04) and third (β = −0.16; 95% CI: -0.30, −0.02) tertiles of maternal "egg, meat, poultry and legume" pattern had a significantly lower HAZ compared to those in the first tertile (Table 5 and Additional file 1: Table S5).
A 0.22 (β = 0.22; 95% CI: 0.06, 0.39) and 0.19 (β = 0.19; 95% CI: 0.04, 0.33) increase in HAZ was found for those in the third tertiles of the "dairy, vegetable and fruit" patterns for children 24-59 months and 6-59 months, respectively. The AIC and BIC were significantly lower for the household and maternal dietary patterns compared to the corresponding DDSs (Table 5). In the joint classification, children in the first tertile of maternal "egg, meat, poultry and legume" and the third tertile of child "dairy, vegetable and fruit" patterns had a 0.31 (β = 0.31; 95% CI: 0.05, 0.57) increase in HAZ compared to the respective third and first tertiles of the patterns. A 0.29 (β = 0.29; 95% CI: 0.07, 0.50) increase in HAZ was found for those children in the third tertiles of both the maternal and child "dairy, vegetable and fruit" dietary patterns, compared to the first tertiles of the patterns (Fig. 3).
The sensitivity analysis, which was undertaken by labelling missing values as "missing" and including them in the analysis, as well as incorporating additional covariates, did not change the findings materially. There were no significant interactions between the DDSs and dietary patterns, and the other covariates and HAZ (stunting). The interactions among dietary patterns for each of the levels (household, maternal and children) were also not significant (data not shown).

Discussion
For households, mothers and children, we computed the DDSs and identified three dietary patterns for each. At all levels, a dietary pattern characterized by high intake  Table 5 Adjusted β coefficients (95% confidence interval) for the associations of household, maternal and child dietary diversity scores and tertiles of dietary pattern scores with childhood height-for-age z score in Ethiopia, 2014 Adjusted β coefficient (95% confidence interval) for height-for-age z score The β coefficients were adjusted for both individual (maternal body mass index, age (if applicable), household food security, number of under-five children in a household, maternal education, maternal height, exclusive breast feeding) and community-level patterns (water source) P for trend was determined by including the tertiles of dietary patterns as continuous variables AIC Akaike's information criterion, BIC Bayesian information criterion *p < 0.05; **p < 0.01 of dairy, vegetables and fruits was positively associated with HAZ and inversely related to stunting. In addition, a maternal "plant-based" pattern (characterized by high intake of grains, white tubers, roots, plantains and dark leafy vegetables) was inversely and significantly associated with childhood stunting. Statistically non-significant associations were found between household, maternal and child DDSs and stunting. Dietary patterns show a better picture of eating habits compared to DDSs by reflecting mainly the behavioural aspect of food consumption that has a synergistic effect on health. Evidence shows that dietary patterns identified by factor analysis are associated with several noncommunicable diseases [28,29], highlighting the plausibility and validity of the approach. In low-and middle-income countries, a study also indicated that the application of factor analysis, using disaggregated food expenditure data to explore food consumption patterns, is an important approach to identify patterns and food groups that predict children's nutritional status [19]. The dietary patterns are defined based on the factor loadings of individual food items which contribute at different levels. Unlike DDSs, factor analysis is an a posteriori statistical analysis method that creates unrelated food patterns that could potentially be associated with an outcome [17]. In addition, using this approach, it is possible to assess the relative intake level of the individual food groups within a dietary pattern. Dietary patterns defined the overall characteristics of the dietary habits of the study groups. In this study, we found that the "dairy, vegetable and fruit" pattern was a common feature of household, maternal and childhood dietary habits.
The results indicate that a pattern characterised by a high intake of dairy, vegetables and fruits was positively associated with HAZ (inversely related to stunting). However, the proportion of households (32.7%), mothers (28.0%) and children (30.0%) consuming dairy products was low. In a recent study in Malawi, it was reported that frequent milk intake during pregnancy was positively associated with birth size [30]. In Ethiopia, a higher intake of cow's milk, in addition to cereals and/or legumes, was associated with a higher length-for-age zscores among children aged 5-11 months [31]. In our study area, cow's milk is the most commonly consumed type of dairy product. It is believed that the milk contains important nutrients, including protein, calcium and vitamin A, which are important for development and bone growth [32]. A systematic review of dairy consumption and physical growth has shown that a daily intake of 245 ml of milk is associated with 0.4 cm increase in height per annum compared to non-consumers [33]. It was also reported that low consumption of vegetables and fruits was associated with stunting and poor linear growth in children aged 6-23 months [34]. This implies that, in addition to available interventions to increase dietary diversity, targeting to increase accessibility and consumption of dairy, vegetables and fruits specifically, could have an important contribution to the reduction of stunting prevalence. a b Fig. 3 Multivariable adjusted β coefficient and 95% confidence interval of fractures in joint classified participants across nine strata formed with the tertiles of dietary patterns. The β coefficients were adjusted for both individual (maternal body mass index, age (if applicable), household food security, number of under-five children in a household, maternal education, maternal height, exclusive breast feeding) and community-level patterns (water source) Unlike other studies [14,15,35], we found a nonsignificant association of household, maternal and child DDSs with stunting. This difference could be explained by the differences in sample size, sampling and analysis methods. Particularly, with regard to the analysis methods, we used a multilevel Poisson regression model to determine the level of associations. This method allows for controlling for the geographical clustering effect of the samples, and accounting for potential correlations, under the premise that variations in childhood stunting could be due to both individual and community level factors [36]. Most studies that found a positive association between DDS and stunting did not consider this clustering effect [16,37]. The association may also be confounded by income [38,39]. However, in this study, even when we adjusted for income, the association still remained statistically non-significant.
A study by Daniels et al. suggested that the addition of portion size as part of data collection could improve the correlation of DDSs with nutrient adequacy [40]. Another study among Zambian infants suggested that although dietary diversity had a positive effect on linear growth, micronutrient adequacy among those who consumed fortified foods may be more accurately assessed using other food indicators [41]. In Eastern Kenya, a study reported that child DDS was not associated with childhood stunting [42]. In Ethiopia, while household food group expenditure index (identified using factor analysis) significantly predicted HAZ (β = 0.067; p = 0.03), dietary diversity was only marginally associated with HAZ (β = 0.037; p = 0.05) [19]. DDSs are important indicators of dietary quality in terms of micronutrient density and adequacy. However, DDSs only measure one dimension of dietary quality. Macronutrients (for instance, protein) also have an important role in growth and development in children [8]. Although DDS is an important approach to measure dietary quality, we recommend that the use of a posteriori dietary data analysis methods (such as factor analysis) as an alternative or complementary method, can give a further insight into the eating behaviours of a population group. These approaches are also important to understand the relative contribution of foods in a pattern that have a potential link with disease outcomes or nutritional status, eventually leading to identifying specific food items, which are most important in determining an outcome of interest (a disease or nutritional status).
Measures taken to ensure the quality of the data are a major strength of the study. Before, during and after the data collection, all possible quality control measures, including intensive training of data collectors, use of standard procedures and tools, intensive and supportive supervision and standardization of anthropometric measurements to minimize bias and associated errors were implemented. The use of qualitative dietary data without portion size and limited food items for the identification of dietary patterns in the factor analysis could be a limitation. Therefore, further validation studies are needed. In addition, due to the cross-sectional design, we cannot claim a cause-effect relationship between dietary patterns and stunting.

Conclusions
Identification of dietary patterns using a posteriori dietary analysis methods can be an alternative and feasible method of diet quality assessment in LICs as an alternative approach to DDSs. We found that, while DDSs are not significantly associated with HAZ (stunting), a dietary pattern characterized by a high intake of dairy, vegetables and fruits by households, mothers and children is positively associated with HAZ and inversely associated with stunting. These findings could be of importance in developing food-based interventions targeting households, mothers and children. In addition, the study suggests an alternative approach of analysing dietary data to determine dietary quality using an a posteriori method with the same data collected for DDSs. More research is warranted to confirm the findings.

Additional file
Additional file 1: Table S1. Components of food groups used for dietary pattern (household and maternal), 2014. Table S2. Components of food groups used for dietary pattern (children), 2014. Table S3. Proportion children, mothers and household by dietary diversity scores (DDS), 2014. Table S4. Dietary diversity and stunting in Ethiopia, 2014. Table S5. Prevalence ratios (95% confidence interval) for the associations of household, maternal and child dietary diversity scores and dietary patterns with childhood stunting in Ethiopia, 2014 (PDF 173 kb)