Predictors for the duration of breastfeeding among ethiopia women of childbearing age with babies; application of accelerate failure time and parametric shared frailty models | BMC Nutrition | Full Text

2022-09-24 04:28:04 By : Mr. Devlin Zhang

BMC Nutrition volume  8, Article number: 106 (2022 ) Cite this article

Duration of breastfeeding is the length of the time that infants who were initially breastfed continue to receive breast milk until weaning. The duration of breastfeeding is important for a child's health, growth, and development. However, the duration of breastfeeding decreases from time to time and further leads children to be exposed to malnutrition (stunting, wasting, and weight loss). Children who did not get enough breastfeeding are also exposed to different diseases. Previous studies used a simple survival model and didn’t see the shared frailty model on the variable of interest. Therefore, the current study aimed to investigate the factors affecting the duration of breastfeeding among Ethiopian women of reproductive age with babies.

A cross-sectional study design was conducted on 15,400 women of childbearing age with babies in nine regional states and two city administrations. The data source for the analysis was the 2016 EDHS data. The Cox-proportional hazard model, AFT, and parametric shared frailty models were conducted for the current investigation. Weibull-gamma shared frailty model was in favor of others for current data analysis.

Among the covariates, women living in urban area (Φ = 0.96; 95% CI; (0.94,0.97); p-value = 0.001), non-educated women(Φ = 1.03; 95% CI; (1.00,1.06); p-value = 0.039), primary educated women (Φ = 1.13; 95% CI; (1.11,1.15); p-value < 0.001), age of a child (Φ = 0.99; 95% CI; (0.76.0.99); p-value < 0.001) and non-smoker mothers (Φ = 1.60; 95% CI; (1.57, 1.63); p-value < 0.001),birth interval between 2–3 years(Φ = 1.02; 95% CI;(1.09, 1.25, p-value = 0.027), birth interval, > 3 years(Φ = 1.28; 95% CI; (1.06, 1.43); p-value < 0.01 significantly affected the duration of breastfeeding. The median survival time of breastfeeding of women of reproductive age with babies considered under study was 23.4 months. Clustering had a significant effect on the variable of interest.

Residence area, level of education, age of the child, smoking status of women, and birth interval of successive birth significantly affected the duration of breastfeeding in the current investigation. Hence, the health staff should conduct health-related education for young women, educated women, urban women, smoker women, and women with a shorter interval of birth to increase the women's attitude and awareness towards the use of long-duration of breastfeeding.

Breastfeeding is the act of breast milk transference from the mother to the infant. Breastfeeding is a unique source of nutrition that plays an important role in the growth, development, and survival of infants [1]. It satisfies a baby’s emotional needs and increases bonding between mother and infant, and it is considered to be physiologically, psychologically, and immunologically important [2]. Breastfeeding reduces the incidence of meningitis, malaria, asthma, respiratory diseases (such as pneumonia), ear infection, diarrhea, and urinary tract infection [3, 4].Breastfeeding also helps to reduce age-related disease (severe acute malnutrition (kwashiorkor) [5].Breastfeeding of an infant is not only crucial for optimal growth and development but also important for determinants of future physical and mental well-being because of the rapid growth and development of organs and tissues during the first year of life [6].

However, the duration of breastfeeding is declining in almost all parts of the world despite its nutritional and immunological benefits [7,8,9,10].An estimated number of 820,000 under-five children died because of the short duration of breastfeeding [11, 12].

In Ethiopia, the duration of breastfeeding decreases from 74 to 64% for infants of age group (0–1 month) and from 74 to 36% for infants of age group (4–5 months) [13]. The duration of breastfeeding decreases from 91 to 76% among children of the age group (18–23 months) [14]. Hence, the duration of breastfeeding is decreasing significantly and becomes below the WHO’s and UNICEF’s recommendations in the study area [15].

An insignificant number of researches are conducted on the issue concerned here, and most of the research conducted previously considered a simple survival model without considering the shared frailty term [16]. This gap became a motivation for the current investigation. To overcome, this gap, in the current study, the duration of breastfeeding of women were clustered by the region by introducing the shared frailty term in the usual survival model. The shared frailty term was also added to account for investigating the correlation that existed among the cluster of unobservable random effects. Hence, the current study used a parametric shared frailty model and accelerated failure time (AFT) models like Weibull, exponential, log-logistic, and log-normal distributions to compare and get the best model which fits the duration of breastfeeding of women data appropriately.

For the effective policy formulation and motivation of women for a longer duration of breastfeeding, it is crucial to study the various socio-economic and demographic factors affecting the duration of breastfeeding. Therefore, the objective of the current investigation was to identify socio-demographic, economic, and health-related factors affecting the duration of breastfeeding among Ethiopian women of reproductive age with babies.

The current study was conducted across the country, Ethiopia, which consists of nine regions and two city administrations namely; Tigray, Afar, Amhara, Oromia Somali, Benishangul-Gumuz, the Old Southern Nations Nationalities and Peoples (SNNP), Gambela and Harari regional states and two city administrations, namely Addis Ababa and Dire Dawa. The study population under the current investigation consists of all women of childbearing age.

The source of data for this study was reports of the Ethiopian Demographic and Health Survey (EDHS, 2016) which is obtained from the Central Statistical Agency (CSA) collected from January 18, 2016, to 27 June 2016. It is the fourth survey conducted in Ethiopia as part of the worldwide DHS project. The 2016 EDHS is designed to provide estimates for the health and demographic variables of interest.

Administratively, regions in Ethiopia are divided into zones, and zones, into administrative units, woredas. Each woreda is further subdivided into the lowest administrative unit, Kebele. Furthermore, each kebele was subdivided into census enumeration areas (EAs) or clusters. The 2016 EDHS sample was selected using a stratified, two-stage cluster sampling design (EDHS, 2016). Each region was stratified into urban and rural areas, yielding 21 sampling strata. Clusters are the sampling units for the first stage. The sample included 645 clusters (202 in urban areas and 443 in rural areas). Then, a fixed number of 28 households per cluster were selected with an equal chance of selection criteria from the newly created household listing. A total of 18,008 households were selected for the sample, of which 17,067 were interviewed. Among the interviewed households, 15,683 eligible women of reproductive age with babies were identified for individual interviews. Hence, this study considered the total number of 15,400 women of reproductive age with babies.

All women of reproductive age (15–49 years) with babies were eligible for the current investigation.

The duration of breastfeeding was considered as the dependent variable and defined as the period between the date of the beginning of breastfeeding and the time to stop breastfeeding before the end of the study.

The predictor variables considered under current investigation are indicated in Table1, in the results section.

In this study, SAS software (Version 9.4) was used for statistical data analysis. Semi-parametric methods such as Cox proportional hazards model were used for multivariable analysis. A comparison of the survival curve for each covariate was also conducted using the Kaplan–Meier method. The accelerated failure time model was also used to measure the direct effect of the explanatory variables on the survival time instead of a hazard. A shared frailty model was used to model the degree of correlation within groups assuming that all individuals in a given cluster share the same frailty; i.e., the subjects within a group are correlated because they share common frailty [17,18,19].

The existence of confounding was checked and controlled by comparing the estimated measure of the association before and after adjusting for confounding. In other words, the measure association was computed both before and after adjusting for the potential confounding factors. The difference between the two measures of association was 4% which is less than 10% and this indicates that there was little and insignificant confounding. The formula used for computing confounding is given by.

Magnitude of confounding = \(\frac{{OR}_{crude}-{OR}_{adjusted}}{{OR}_{crode}}\) , where OR stands for odds ratio.

Effect modification/ interaction effect occurs when the effect of a factor is different for different groups. The existence of effect modification was tested by considering the crude estimate of the association (odds ratio). The result indicates that there was very close to a weighted average of group-specific estimates of the association. In the current investigation, both stratum-specific estimates of the odds ratio went to one side of the crude odds ratio. Hence, the result of this study indicates that the crude odds ratio was less than the estimates of the odds ratio for the stratum-specific estimates which is an evidence for nonexistence of effect modification.

Multicollinearity in regression analysis occurs when two or more explanatory variables are highly correlated to each other, such that they do not provide unique or independent information in the regression model. Multicollinearity was tested in the current investigation using a metric known as Variance Inflection Factor (VIF). The VIF measures the correlation between predictor variables. The result in this regard revealed that there was no multicollinearity issue.

The baseline characteristics of covariates are indicated in Table 1. In Table 1, it is indicated that about 15,400 breastfed women of reproductive age with babies were included in this study. Among these, about 27.6% of the participants were censored, where the actual duration of breastfeeding is unknown and the remaining were events with known breastfeeding duration, 65.7% lived in rural areas, about 42.6% of the women were from poor households, 18.2% of them were from middle-income households and the rest were from high-income households. Table 1 also indicates that about 63.6% of the participants were women with female children and 65.4% of the participants had no job (unemployed). Out of 10,072 unemployed women, 27.7% of them were still breastfeeding and 72.3% of women stopped breastfeeding (event). Of a total of 5328 employed women, 27.5% of the women were still breastfeeding and 72.5% of women stopped breastfeeding.

Among the participants, about 92.1% were non-smokers, 17.9% of the women attained secondary and above education, 49.7% of the women attained primary education, and the rest were uneducated women. The religion category in Table 1 indicates that about 41.1% of the participants were Orthodox, 39.4% were Muslim, and 17.9% of them were Protestant. Regarding their regions, about 11.8% of them were from Addis Ababa city, 12.0% from Oromia, 11.7 from the Old SNNP, 11.0% from Amhara, about 7.1% of them were from Afar, 7.2% from Benishangul- Gumuz, 7.3% from Dire Dawa, 6.5% from Gambela, 5.8% from Harari, and about 10.8% of them were from Tigray region.

The non-parametric analysis was conducted using plots of the Kaplan Meir curves to the survival and hazard experience for the duration of breastfeeding as shown in Fig. 1 and Fig. 2 respectively. The result revealed that the survival plot decreased at an increasing rate at the beginning and decreased further from time to time. This implies that most women breastfed highly soon after birth. On the other hand, the hazard plot increased at an increasing rate at the beginning and increased as time increased.

Survival plot for the duration of breastfeeding

Hazard plot of duration of breastfeeding

The log-rank test was used at a 5% level of significance to validate the differences in the survival time of each factor. The difference between the probabilities of an event occurring at any time point was the null hypothesis that has been tested.

After making a comparison of the survivorship experience among groups of covariates, the next important step was model development. At the initial step, in the model-building process was conducted to identify sets of explanatory variables that had the potential for being included in the linear components of a multivariable proportional hazards model. In doing this, Cox proportional hazards regression model, including the univariable Cox proportional hazards models were fitted. From the univariate analysis; place of residence, wealth index of the family, sex of the child, child age, smoking status of women, birth interval, place of residence, wealth index, and women's education level were statistically significant for the variable of interest.

In the current investigation, the two basic assumptions of the Cox regression model, log-linearity and proportional hazards were tested as shown in Table 2. The log-linearity test revealed that the relationship between log hazard or log cumulative hazard and a covariate was linear. The proportional hazard test in this investigation indicates that the ratio of the hazard function for two individuals with different regression covariates does not vary with time.

The global fit test in Table 2 also shows that the Wald chi-square test statistic was significant which indicates that the proportional hazards assumption is violated. In other words, the plots of the duration of breastfeeding covariates were not parallel to each other.Thus, there is a violation of the proportional hazards assumption. This indicates that the residuals were not random, there is a systematic pattern and the smoothed plot looks not a straight line and there is some departure from the horizontal line. Thus, there is a violation of the proportional hazards assumption.

Since the proportional hazard assumptions were not satisfied, the accelerated failure time (AFT), including the univariable and multivariable analysis, the model should be conducted for the current data analysis. The univariable analyses were fitted for every covariate, considering AFT models for baseline characteristics of participants. For the duration of breastfeeding data, AFT models of Weibull, exponential, log-logistic, and log-normal distribution were considered. To select the best model for the current analysis, AFT models namely, Weibull, exponential, log-logistic and log-normal distributions were compared using AIC and BIC considering the model with the smallest AIC and BIC is the one that fits the data well, as shown in Table 3. Table 3 indicates that the Weibull distribution had the smallest AIC and BIC. Hence, it has been selected for univariate and multivariate data analysis in the current investigation.

In all univariable analyses of Weibull AFT models, age of women, smoking status, place of residence, wealth index of the household, sex of the child, women's education level, age of children, and birth interval of women were significantly associated with the duration of breastfeeding at 5% level of significance. The summary of the univariable analysis is given in Table 4. Variables that were significant at 5% of the level of significance in the univariable model were included in the multivariable data analysis. The multivariable data analysis under this study is indicated in Table 5. In Table 5, the main effect of the covariates namely place of residence, age of women, birth interval, women’s level of education, smoking status, and age of the child were considered as potential predictors for the variable of interest. The multivariable data analysis with the Weibull AFT models and corresponding AIC and BIC values were conducted as indicated in Table 5.

The result of the Weibull AFT model in Table 5 indicates that level of education, age of women, birth interval, place of residence, sex of the child, smoking status, and wealth status were statistically significant variables for the variable of interest.

To test the effect of regions on the variable of interest, the multivariable survival analysis was conducted including the Gamma shared frailty term, This was performed using the covariates; place of residence, education level of women, religion, employment status, age of the child, sex of a child, wealth index, women’s level of education, smoking status, and birth interval. In this study, the AIC and BIC criteria were considered to compare various candidates of parametric shared frailty models considering, the model with the smallest AIC and BIC as the best model.

The parametric baseline distributions namely Gamma frailty distribution, log-normal, and Inverse Gaussian distributions were fitted and compared by considering the regions of the women as frailty terms. The effect of random component (frailty) was significant for Weibull gamma shared frailty because of its smallest AIC and BIC values [20]. The final Weibull gamma shared frailty model is indicated in Table 6.

Table 6 indicates that the residential area had a significant effect on the duration of breastfeeding. Hence, the expected duration of breastfeeding for urban women was decreased by 4% as compared to rural women, keeping the other covariates constant(Φ = 0.96; 95% CI; (0.94,0.97); p-value = 0.001). This indicates that women in rural had a prolonged duration of breastfeeding than urban women.

The level of education had a significant effect on the variation in the duration of breastfeeding. Hence, comparing non-educated women with secondary and above, the expected duration of breastfeeding for non-educated women was increased by 3% as compared to secondary and above women keeping the other covariates constant(Φ = 1.03; 95% CI; (1.00,1.06); p-value = 0.039). Similarly, the expected duration of breastfeeding for primary educated women was increased by 13% as compared to secondary and above, keeping the other covariates constant(Φ = 1.13; 95% CI; (1.11,1.15); p-value < 0.001). This result indicates that the more educated women had a short duration of breastfeeding as compared to non-educated or less educated women.

The age of a child also played a significant role in the variation in the duration of breastfeeding. Hence, as the age of a child increased by one month, the expected duration of breastfeeding decreased by 1% keeping the other covariates constant (Φ = 0.99; 95% CI; (0.76, 0.99); p-value < 0.001). Hence, the increase in the age of a child leads to a decrease in the duration of breastfeeding.

The smoking status of mothers/women also played a significant role in the variation in the duration of breastfeeding. Comparing smoker women with non-smoker one, the expected duration of breastfeeding for non-smoker women was increased by 60% as compared to smoker women, keeping the other covariates constant. (Φ = 1.60; 95% CI; (1.57, 1.63); p-value < 0.001).

The birth interval between consecutive births had also a significant effect on the variation in the duration of breastfeeding. Hence, comparing a woman whose birth interval, is between 2–3 years with that of less than 2 years, the expected duration of breastfeeding for a woman whose birth interval, is 2–3 years was increased by 2% as compared to a woman whose birth interval, < 2 years, keeping the other conditions constant (Φ = 1.02; 95% CI; (1.09, 1.25); p-value < 0.027). Similarly, the expected duration of breastfeeding for a woman whose birth interval > 4 years was increased by 28% as compared to a woman whose birth interval < 2 years (Φ = 1.28; 95% CI; (1.06, 1.43); p-value < 0.01).

The age of women was also another significant variable for the variation of the duration of breastfeeding. Hence, the expected duration of breastfeeding for a woman whose age, 40–44 years was increased by 4% as compared to a woman of age, 15–19 years, given that the other covariates were constant (Φ = 1.041; 95% CI; (1.01, 1.22); p-value = 0.019) and the expected duration of breastfeeding for a woman of age,35–39 years was increased by 3% as compared to a woman whose age, 15–19 years (Φ = 1.030; 95% CI; (1.01, 1.55); p-value = 0.005). Hence, the more the age of the women leads to the longer duration of breastfeeding for the current investigation.

Place of residence has a significant role in the variation of the time-to-breastfeeding in Ethiopian women of reproductive age with babies. Women who lived in urban areas had a shorter duration of breastfeeding than women who lived in rural areas. Rural inhabitants have usually no access to getting a different alternative to feed their children [21]. The only option is the feeding of breast milk which may result in a longer duration for rural women as compared to urban. This result is supported by another previous study [22]. Most of the time rural women are housewives and they have a greater chance to get their child a day to give their breast milk.

A woman's age has a significant effect on the variable of interest. Women of older age have a longer duration of breastfeeding. This finding is supported by another study [23]. One of the previous studies reported that aged women are highly associated with a longer duration of breastfeeding [24]. The potential reason for this might be that woman in the menopause period has less chance of getting an extra child and she needs to extend her breastfeeding. Most of the time, young women, who made sexual intercourse, may have an extra child and they may be forced to stop her breastfeeding for the current child to prepare for the newly coming child [25].

Women's level of education has a significant effect on the variation of the duration of breastfeeding of Ethiopian women of reproductive age with babies. Educated women have a shorter duration of breastfeeding than non-educated or less educated women. The potential reason for this might be the fact that more educated women are government employed and they are forced to continue their regular time work for a short period (four months) after birth. This result is supported by many other previous studies [26,27,28]. Previously conducted research states that illiteracy has a positive association with the duration of breastfeeding [29].

The wealth index had a significant association with the duration of breastfeeding. Women of rich wealth status families had a shorter duration of breastfeeding than poor wealthy status. The potential reason for this might be the fact that rich women used another alternative for feeding their children. This finding is supported by another study [30]. It is reported that low-income status was associated with a longer duration of breastfeeding.

The smoking status of women has a significant effect on the duration of breastfeeding. Smoker women have a shorter duration of breastfeeding than non-smokers [31]. The potential reason for this might be that smoker women have lower milk production and let-down problems; this may be related to poor lactation management rather than physiological causes. And also smokes immediately after breastfeeding to cut down on the amount of nicotine in her milk during nursing [32]. The other reason might be the fact that smoker women may not have a strong physiological mechanism. This finding is supported by another previously conducted study [33]. Another similar study reported that smoking has a negative effect on the duration of breastfeeding [34].

The birth interval had a significant association with the duration of breastfeeding. Hence, women with a birth interval of 2–3 years had a longer duration of breastfeeding as compared to those women with a birth interval < 2 years. Furthermore, women with more than three years of birth interval have a longer duration of breastfeeding as compared to those women with breath intervals of fewer than two years. Mothers with short period birth intervals should give their milk to the newly coming child rather than the previous. This finding was supported by [35], who reported that breastfeeding was highly associated with birth interval.

The result of Weibull AFT and Weibull-gamma frailty models revealed that place of residence, age of women, level of women’s education, age of the child, smoking status, birth interval, and wealth of household was found to be significant predictors of the duration of breastfeeding of women of reproductive age with babies in Ethiopia. Among these significant predictors, the age of women, and the birth interval of consecutive births had a positive association with the duration of breastfeeding. On the other hand, the level of education, wealth status of women, smoking status of women, and urbanization had a negative association with the duration of breastfeeding.

In the current study, the clustering effect was significant (p-value = 0.000) in Weibull-gamma shared frailty model. This indicates that there was heterogeneity between the regions on the duration of breastfeeding. Weibull-gamma frailty model was in favor of the Weibull AFT model and well fitted the duration of breastfeeding data. The estimated median survival duration of breastfeeding for Ethiopian women of reproductive age with babies was found to be 23 months. The median time of duration of breastfeeding for rural women was 26 months. The concept behind Weibull–the gamma shared frailty model indicates that there was a frailty (clustering) effect on the duration of breastfeeding. Hence, regions played a significant role for the presence of heterogeneity and necessitate the use of shared frailty models.

This research was not without limitation; one of the limitations of current investigation whether the mothers delivered by cesarean or spontaneous was not identified. Hence, the duration of breast feeding by mode of delivery was not considered. Mainly mothers who delivered by cesarean section may fed her baby for shorter duration than mother delivered spontaneously.

As a recommendation, the health staff and those concerned bodies should encourage and advise women of childbearing age to prolong the duration of breastfeeding. The breastfeeding-promotion programs in Ethiopia should give special attention to young mothers, mothers who have a short duration (educated mothers, urban mothers, mothers with high wealth index, and smoker mothers). The health-related education and awareness creation should be conducted by the health staff about the importance of elongating the duration of breastfeeding. Further investigation, with the inclusion of additional variables not included here such as mode of delivery of women, are recommended to have additional information on factors affecting the duration of breastfeeding in the study area. Further studies should be also conducted in each region of Ethiopia and identify other factors that are not included in this study.

Andreas NJ, Kampmann B, Le-Doare KM. Human breast milk: A review on its composition and bioactivity. Early Human Dev. 2015;91(11):629–35.

Nascimento MBRD, Issler H. Breastfeeding: making the difference in the development, health and nutrition of term and preterm newborns. Rev Hosp Clin Fac Med Sao Paulo. 2003;58(1):49–60.

Benguigui Y, Stein F. Integrated management of childhood illness: an emphasis on the management of infectious diseases. Semin Pediatr Infect Dis. 2021;15(13):128. Elsevier.

Scott JAG, et al. Pneumonia research to reduce childhood mortality in the developing world. J Clin Investig. 2008;118(4):1291–300.

CAS  PubMed  PubMed Central  Google Scholar 

Gedefaw M, Berhe R. Determinates of childhood pneumonia and diarrhea with special emphasis to exclusive breastfeeding in north Achefer district, northwest Ethiopia: a case control study. Open J Epidemiol. 2015;5(02):107.

Miniello VL, et al. Complementary Feeding and Iron Status:“The Unbearable Lightness of Being” Infants. Nutrients. 2021;13(12):4201.

CAS  PubMed  PubMed Central  Google Scholar 

Delgado C, Matijasevich A. Breastfeeding up to two years of age or beyond and its influence on child growth and development: a systematic review. Cad Saude Publica. 2013;29(2):243–56.

Huang K, Atlas R, Parvez F. The significance of breastfeeding to incarcerated pregnant women: an exploratory study. Birth. 2012;39(2):145–55.

PubMed  PubMed Central  Google Scholar 

Larsen JS, Kronborg H. When breastfeeding is unsuccessful–mothers’ experiences after giving up breastfeeding. Scand J Caring Sci. 2013;27(4):848–56.

Turin CG, Ochoa TJ. The role of maternal breast milk in preventing infantile diarrhea in the developing world. Curr Trop Medicine Rep. 2014;1(2):97–105.

Holla-Bhar R, et al. Investing in breastfeeding–the world breastfeeding costing initiative. Int Breastfeed J. 2015;10(1):1–12.

Gupta A, et al. The world breastfeeding trends initiative: implementation of the global strategy for infant and young child feeding in 84 countries. J Public Health Policy. 2019;40(1):35–65.

Abdulahi M, Fretheim A, Magnus JH. Effect of breastfeeding education and support intervention (BFESI) versus routine care on timely initiation and exclusive breastfeeding in Southwest Ethiopia: study protocol for a cluster randomized controlled trial. BMC Pediatr. 2018;18(1):1–14.

Dibley MJ, Senarath U, Agho KE. Infant and young child feeding indicators across nine East and Southeast Asian countries: an analysis of National Survey Data 2000–2005. Public Health Nutr. 2010;13(9):1296–303.

Salim YM, Stones W. Determinants of exclusive breastfeeding in infants of six months and below in Malawi: a cross sectional study. BMC Pregnancy Childbirth. 2020;20(1):1–8.

Banerjee S, Wall MM, Carlin BP. Frailty modeling for spatially correlated survival data, with application to infant mortality in Minnesota. Biostatistics. 2003;4(1):123–42.

Balan TA, Putter H. A tutorial on frailty models. Stat Methods Med Res. 2020;29(11):3424–54.

PubMed  PubMed Central  Google Scholar 

Gebeyehu A. Survival analysis of time-to-first birth after marriage among women in Ethiopia: application of parametric shared frailty model. Jimma: Department of Statistics College of Natural Sciences, Jimma University; 2015.

Govindarajulu US, et al. Frailty models: applications to biomedical and genetic studies. Stat Med. 2011;30(22):2754–64.

Fentaw KD, et al. Time to first antenatal care visit among pregnant women in Ethiopia: secondary analysis of EDHS 2016; application of AFT shared frailty models. Arch Public Health. 2021;79(1):1–14.

Silver A. Twilight" is not good for maidens: gender, sexuality, and the family in stephenie meyer’s" twilight series. Stud Novel. 2010;42(1/2):121–38.

Amin T, Hablas H, Al Qader AA. Determinants of initiation and exclusivity of breastfeeding in Al Hassa Saudi Arabia. Breastfeeding Med. 2011;6(2):59–68.

Fallahzadeh H, Farajpour Z, Emam Z. Duration and determinants of birth interval in Yazd, Iran: a population study. Iran J Reprod Med. 2013;11(5):379.

PubMed  PubMed Central  Google Scholar 

Fisk CM, et al. Breastfeeding and reported morbidity during infancy: findings from the Southampton Women’s Survey. Matern Child Nutr. 2011;7(1):61–70.

Kimani-Murage EW, et al. Factors affecting actualisation of the WHO breastfeeding recommendations in urban poor settings in K enya. Matern Child Nutr. 2015;11(3):314–32.

Donath, S.M., L.H. Amir, and A.S. Team. Relationship between prenatal infant feeding intention and initiation and duration of breastfeeding: a cohort study. Acta Paediatr. 2003;92(3):352–6.

Brown A, Davies R. Fathers’ experiences of supporting breastfeeding: challenges for breastfeeding promotion and education. Matern Child Nutr. 2014;10(4):510–26.

PubMed  PubMed Central  Google Scholar 

Meedya S, Fahy K, Kable A. Factors that positively influence breastfeeding duration to 6 months: a literature review. Women Birth. 2010;23(4):135–45.

Li L, et al. Factors associated with the initiation and duration of breastfeeding by Chinese mothers in Perth, Western Australia. J Hum Lact. 2004;20(2):188–95.

Alemayehu T, Haidar J, Habte D. Determinants of exclusive breastfeeding practices in Ethiopia. Ethiop J Health Dev. 2009;23(1):122–31.

Giglia R, Binns CW, Alfonso H. Maternal cigarette smoking and breastfeeding duration. Acta Paediatr. 2006;95(11):1370–4.

Napierala M, et al. Tobacco smoking and breastfeeding: effect on the lactation process, breast milk composition and infant development A critical review. Environ Res. 2016;151:321–38.

Amir LH. Maternal smoking and reduced duration of breastfeeding: a review of possible mechanisms. Early Human Dev. 2001;64(1):45–67.

Amir LH, Donath SM. Does maternal smoking have a negative physiological effect on breastfeeding? Epidemiol Evid Birth. 2002;29(2):112–23.

Yohannes S, et al. Duration and determinants of birth interval among women of child bearing age in Southern Ethiopia. BMC Pregnancy Childbirth. 2011;11(1):1–6.

All the health staff at each government hospital is gratefully acknowledged for the data they supplied for our health research.

Authors declared that there are no conflicts of interest in this work.

This research did not receive any specific grant from funding agencies in the public, commercial, or not-from-profit sectors.

Department of Statistics, Debre Berhan University, Debre Berhan, Ethiopia

Getahun Mulugeta & Dagne Tesfaye

Department of Statistics, Bahir Dar University, Bahir Dar, Ethiopia

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

All authors made substantial contributions to conception and design, acquisition of data, took part in drafting it critically for important intellectual contact, agreed to submit to this journal, gave final approval of the version to be published, and agreed to be accountable for all accepts of the work. The author(s) read and approved the final manuscript.

Correspondence to Awoke Seyoum Tegegne.

To secure the confidentiality and compliance with the Declaration of Helsinki within the manuscript of patient-related data, the name of patients was not given to investigators, rather id number and important variables related to the current investigation were given to researchers. The waiver was done by Debre Birhan University Ethical approval committee, Ethiopia with reference number: RCS/1412/2021. Hence, the Debre Birhan University Ethical Committee approved and waived this study. The data used in the current investigation were secondary and in the data collection procedures, only id number of patients and important variables related to the current investigation were given to researchers. The name of the ethics committee who waived the informed consent for the current investigation is ‘’Debre Birhan University Ethical approval committee’’, Ethiopia with reference number: RCS/1412/2021.

Authors declared that there is no conflict of financial interest between the author and institutions.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Mulugeta, G., Tesfaye, D. & Tegegne, A.S. Predictors for the duration of breastfeeding among ethiopia women of childbearing age with babies; application of accelerate failure time and parametric shared frailty models. BMC Nutr 8, 106 (2022). https://doi.org/10.1186/s40795-022-00601-z

DOI: https://doi.org/10.1186/s40795-022-00601-z

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By using this website, you agree to our Terms and Conditions, California Privacy Statement, Privacy statement and Cookies policy. Manage cookies/Do not sell my data we use in the preference centre.

© 2022 BioMed Central Ltd unless otherwise stated. Part of Springer Nature.