Analyzing factors associated with time to age at first marriage among women in Ethiopia: log logistic-gamma shared frailty model
BMC Women's Health volume 22, Article number: 191 (2022)
The main objective of this study is to fit Log logistic-Gamma shared frailty model for the determinant of time to age at first marriage among women in Ethiopia.
The data set in this study were obtained from Demography and Health survey conducted in Ethiopia in 2016. In this study, we used Log logistic-Gamma shared frailty model to account for the loss of independence that arises from the clustering of women in region of Ethiopia. A total of 12,066 women aged 15–49 in Ethiopia were included in this study.
Of all 12,066 women aged 15–49, 9466 (78.45%) were married and the median & mean age at first marriage for women living in Ethiopia were 17.2 years and 17.5 years respectively, while the minimum and maximum age at first marriage observed were 8 years and 49 years respectively.
The most significant contributing factors to delaying time to age at first marriage of women aged 15–49 in Ethiopia were increased education level of women, increased education level of the head, increased income, residing in urban and being followers of religion other than orthodox, catholic, protestant & Muslim. The heterogeneity of age at first marriage for women aged 15–49 among regions in Ethiopia was observed. The government of Ethiopia and the concerned bodies should revise the women's health policy and practice to reduce early marriage and give attention to women; illiterate, live in rural areas, and have illiterate and poor heads.
Time to age at first marriage is length of time from birth until the age by which the individuals have been married for the first time [1, 2]. For individuals, it is the length of time up to a new chapter of life [1, 3]. The couples may lead a new chapter of life by strengthening their economy, and they may have children over time thru birth or adoption [2, 3]. In Ethiopia, women tend to have tremendously less time to age at first marriage than men. The median age at first marriage is 16.5 years among women age 25–49 and 23.1 years among men age 25–49. 58% of women and only 9% of men age 25–49 marry before their 18th birthday [4, 5]. Our main motivational factors to conduct a depth study on time to age at first marriage is that the effects that they could have on reproductive health and socio-economic pertaining issues in Ethiopia. This is due to the fact that, time to age at first marriage, especially among girls is often associated with adolescent motherhood, school dropouts, maternal morbidity and mortality, and forfeited future life opportunities for the affected individual[1, 6]. In Ethiopia, a finding from study in Injibara town reveals that, 44.8% of women in Injibara town had first married below 18 years. In this area, the mean age at first marriage was 17 years, while the minimum and maximum age at first marriage was 9 and 23 years respectively . The result from another study in Ethiopia particularly, in Dabat district shows that 86.6% of women in Dabat district had first married below 18 years. These studies concluded that women in these areas had first married at earlier time .
Most study in Ethiopia conducted using logistic regression and Cox proportional hazard models to analysis the determinant of time to age at first marriage [4, 5]. But, those models have their own amiss and cannot be applied for multilevel survival data. For instance, logistic regression does not provide methods for survival data and a Cox proportional hazard model require independency assumption of samples [9, 10]. To overcome this situation, we used Log logistic-Gamma shared frailty model. This model permits the analysts to account for the loss of independence that arises from the clustering of subjects in higher level units. Similarly, it allows researchers to make valid inferences when examining the effect of both subject characteristics and cluster characteristics on the risk of the occurrence of the outcome. It is a model that incorporates cluster-specific random effects to account for within cluster homogeneity in outcomes . For instance, in this study women with in the same region share the same random effect with respect to marriage. As a result, here we used a Log logistic-Gamma shared frailty model for our two level data.
The general objective of this study is to analysis factors associated with time to age at first marriage among women in Ethiopia using advanced methodology. Distinguishing significant associated factors with time to age at first marriage among women in Ethiopia and estimating the clustering effect in modeling time to age at first marriage for the data set are the specific objectives of this study. Conducting this study has a great role in designing women's health policy and identifying priority health programs in Ethiopia. The study also set a solution for women’s health-related problems in Ethiopia.
Study design and setting
This study used a population-based cross-sectional design and data from 2016 Ethiopian Demographic and Health Survey (EDHS), which was gathered between January 18 and June 27, 2016. The purpose of the survey was to give estimates for demographic and health predictors of interest in Ethiopia's two administrative cities and nine geographical regions. The survey included 12,066 women between the ages of 15 and 49.
Variable in the study
Time to age at first marriage (dependent variable) is a measured value in years from birth to age at first marriage. On a sample of all women aged 15–49 in Ethiopia, we retrospectively observe the timing to first marriage since birth. Hence we have to consider two things. First, all cases with no observed events are right censored. Therefore the women who had not yet experienced the event of interest resulting in right censoring of the data. Second, in order to make censoring valid, we have to assume that all women marry before the age of 50.
Wealth index, media exposure, head/parents' occupations, women's education level, head/parents' education level, type of residence, religion, and Respondents' work status are the independent variables investigated in this study. In the frailty model we have considered regional state of the women as a clustering effect.
Method of data analysis
Multilevel survival analysis
Multilevel survival analysis is the statistical technique that can apply for clustered (grouped) survival times. Researchers often encounter grouped or multilevel data like individuals are nested within families, and families are nested within neighborhoods. In our study also encountered such kinds of data. For instance women aged 15–49 nested with in region i.e. we have a two level data. Analyzing such data requires special treatment because most multivariate models assume that observations are independent, and grouped data clearly violate this assumption. Statisticians and biomedical researchers identified adverse consequences of applying the Cox regression to grouped survival times [11, 12]. They noted that when the independent assumption of the Cox model is violated, the tests of statistical significance are biased and in ways that cannot be predicted beforehand . Since Log logistic-Gamma shared frailty model is one of the statistical models for multilevel survival analysis , we used it in this study.
Modeling frailty and shared frailty model
A model for multilevel survival data may consider the heterogeneity among individual failure times and the heterogeneity between groups of failure times. If we add a random effect from the proportional hazards model the resulting model is called the frailty model. The inclusion of this random effect in this model helps to measure the heterogeneity among individuals. The term “frailty” is used to express heterogeneity. In addition, to explicitly measure the heterogeneity among individuals, this model allows measuring truly the effects of covariates on the hazard function of individuals [14, 15]. In another way, if we add a group-specific random effect from the proportional hazards model the resulting model is called the shared frailty model. This model is used to know the heterogeneity among individual failure times and the heterogeneity between groups of failure times. This group-specific random effect can be considered as an unobserved covariate which is in common for that specific group. The term “shared frailty” comes from the concept that individuals within the same group share the same heterogeneity [9, 14, 15].
If G and ni represents the number of groups and individuals in ith group respectively, the hazard function of the jth individual in ith group is.
where wi is ith group unobservable covariates, ψ is regression coefficient for ith group unobservable covariates, Zij is the jth individual observable covariate vector in ith group. exp(ψwi) is the frailty of the ith group and if we replace it by ui, the hazard function can be simplified as:
The random variables u1, . . . ,uG assumed to have a common probability density function q.
This random effect model uses u as a random effect to represent group variation and the hazard function h0(t) exp(βZij) to represent individual variation. When a group's frailty value is high, the group experiences the event of interest at earlier times [9, 14, 15].
The gamma frailty model
In gamma frailty model, the frailty assumed to have a gamma distribution. It fits failure time very well. Mathematical tractability and analytical simplicity are the two main benefits of this model. Those benefits are the result of mathematical derivative-able of the Laplace transformation of the model [15, 16].
A frailty's density function in one parameter gamma distribution (θ) is.
If θ > 0 and u > 0, individuals in group experience the event of interest at earlier times (more frail). If u < 0, individuals in group experience the event of interest at latter times (less frail).The Kendall's Tau  measure the relationship between any two events times within group for the Gamma distribution and is given by:-
The Laplace transform of the model in (2.3) is expressed as; LP(s) = (1+θs)−1/θ.
Log logistic-gamma shared frailty model
Log-logistic distribution for Baseline hazard
The log-logistic distribution is a hybrid of the Gompertz and Gamma distributions, with the mean and variance values equal to one [20, 21].Various functional forms for f(t), S(t), h(t) & H(t) in log-logistic distribution are presented in (Table 1).
From a total of 12,066 women aged 15–49 in Ethiopia, 9466(78.45%) were married. The median & mean age at first marriage among women aged 15–49 were 17.2 years and 17.5 years respectively. In addition, the minimum and maximum age at first marriage observed were 8 years and 49 years respectively. The highest and lowest percentage of women aged 15–49 who were first married were in the Amhara region (11.14%) and in the Somali region (4.81%) respectively. Both median and mean age at first marriage among women in Addis Ababa was 19 years which was the highest relative to other regions. In contrast, the median and mean age at first marriage among women in the Amhara region was 14 years and 14.5 years respectively which were the lowest relative to other regions. Women aged 15–49 who were first married and reside in the rural part of Ethiopia were higher in percentage (57.97%) than those residing in the urban area (20.48%). Both median and mean age at first marriage among women in the urban part of Ethiopia were 18 years and 18.5 years respectively which were higher than women in the rural area (median = mean = 16 years).
Regards to the educational level of women, the percentage of illiterate women aged 15–49 who were first married was highest (50.39%), while the lowest for women at higher education level (2.54%). The median and mean age at first marriage among women aged 15–49 are 23 years and 22 years respectively for those women who achieved higher education, which is highest, while lowest for illiterate women (mean = 15, median = 17). The finding also shows that women aged 15–49 who were first married and followers of the Muslim religion in Ethiopia were highest in percentage (31.87%), while it is lowest for catholic religion followers (0.84%).
Both the mean and median age at first marriage among women aged 15–49 in Ethiopia was 18 years which is similar and highest for women who follow a religion other than orthodox, Muslim, catholic, and protestant, while it is lowest for the orthodox follower (mean = 16, median = 17).
The percentage of women aged 15–49 who were first married was highest for the women that having richest head (21.95%) and followed by the women that having poorest head (20.45%), while the lowest for those having a middle income head (11.78%). Median and mean age at first marriage were highest and similar for 15–49 aged women who having richest head in Ethiopia, while the lowest for those having poorest head (mean = 16.5, median = 17). The result also shows that the percentage of women aged 15–49 who were first married was highest for those women having an illiterate head (39.87%) and agriculturalist head(50.91%), while lowest for those having a better-educated head (4.91%) and professional head (4.46%). However, median age and mean age at first marriage are similar as well as highest (20 years) for those women having head who achieved higher education, while similar as well as lowest for those of having no educated head (16 years).
The percentage of women aged 15–49 who were first married was higher for media-accessible women (48.22%) than those who haven’t (30.23%). Median age and mean age at first marriage for women aged 15–49 are higher among those having better access to media and wealthier. The median age at first marriage among women aged 15–49 in Ethiopia was 18 years for women having professional and laborers head which were the highest, while the lowest for women having an agricultural head (16 years). Similarly, the mean age at first marriage was highest for women having a professional head (18.5%) and lowest for women having an agricultural head (16%).
Finally, the percentage of women aged 15–49 who were first married was higher for those who haven't a job (work) (50.53%) than those having work (27.92%). The median and mean age at first marriage is similar, which is 16 years irrespective of work status (Table 2).
Analysis based on log logistic-gamma shared frailty model
Analysis based on this model showed that place of residence of women, religion of women, education level of women, some categories at wealth index of household and head education were significant at 5% level of significance (P value < 0.05). In contrast work status of women, head occupation and access to media were not significant at 5% level of significance (P value > 0.05).
An odds ratio of greater than 1 indicates delaying the time to age at first marriage. Therefore women who attend primary school (ϕ = 1.034), secondary school (ϕ = 1.183) and higher education (ϕ = 1.325) had delayed age at first marriage by a factor of 1.034, 1.183 & 1.325 respectively than the reference category (No Education). Concerning religion of women, women who are follower of catholic; protestant, Muslim and other religion (religion except orthodox, catholic, protestant & Muslim) had delayed age at first marriage by a factor of 1.047, 1.05389, 1.05387 and 1.067 respectively when compared to the reference category (orthodox). Head education is also another significant covariate with odds ratio greater than 1 for all categories. As a result those women having head that attend primary school, secondary school and higher education had delayed age at first marriage by a factor of 1.016, 1.029 and 1.044 respectively when compared with the reference category (No Education). Categories of Significant covariates having odds ratio less than 1 imply that women characterized by those categories of the same covariate had less delayed age at first marriage relative to the reference category of the same covariates. For instance, women residing in rural area of Ethiopia have less expected survival time (duration of age at first marriage) than those residing in urban area (ref) of Ethiopia (ϕ = 0.967). Similarly, richer (ϕ = 0.970) and richest women (ϕ = 0.949) had less delayed age at first marriage than reference category (poorest women). Here, we see that poorer and middle categories are not statistically significant at 5% level of significance (P value > 0.05).
The value of the shape parameter in the Log logistic –Gamma shared frailty model was (γ = 0.131). This indicates that non monotonic hazard rates, specifically initially increasing and then decreasing rates. The variability (heterogeneity) of age at first marriage in women population of clusters (Region) estimated by this Log logistic –Gamma shared frailty model was θ = 0.088 and the dependence within region were about τ = 4.2% (Table 3).
At 5% level of significance, the results from Log logistic-Gamma shared frailty model shows that, residence of women, educational level of women, religion of women, wealth index of household and head educational level are the most significance factors in predicting time to age at first marriage, whereas work status of women, access to mass media, and head occupation are not. A study by Haloi  revealed that wealth index of household and head education is significant factors in affecting time to age at first marriage. This result coincides with our findings. Moreover, the result from our study shows that, women who had better educated parent have delayed age at first marriage than women who had least educated parent. With respect to wealth index, women who had richest parent have less delayed age at first marriage than women who had poorest parent. This might be explained from the fact that, the increment of education level and income of household head leads to the increment of the perception and knowhow of household head for the consequence of early marriage. In addition, in Ethiopian context, household head have a direct or indirect involvement in deciding the age of first marriage for their daughters. This implies wealthier and educator household head have a great role in delaying age at first marriage.
The finding from a study by Adebowale  revealed that residence of women, educational level of women and religion of women had a significant effect on age at first marriage. This result also coincides with our findings. In our finding, better educated women had delayed age at first marriage than least educated women. The possible reason might be the crucial role of education in creating and developing knowledge and attitude towards the impact of early marriage. Concerning to religion of women, orthodox follower women had less delayed age at first marriage than other religion follower’s women. This finding is in line with the fact that, preachers in orthodox religion of Ethiopia didn’t focused on religion education about marriage. So that, their followers might be in lack of awareness about the impact of early marriage and then have less delayed age at first marriage. With regards to residence of women, women residing in rural area of Ethiopia have less expected time to age at first marriage than those residing in urban area of Ethiopia. In Ethiopian context, parents in rural area take early marriage as a cultural set up to avoid sexual risk (unintended pregnancy, rape and abduction) of their children. In addition, women in rural area are less educated and have no decision power for their age at first marriage. According to Von  findings work status of women, educational level of women and residence of women had a significant effect on age at first marriage but not head occupation of women. This finding contradicts our finding with respect to work status. Access to mass media was found to have a significant effect on age at first marriage by Zahangir . But, in our study it was not. The findings from another study in Ethiopia reveals that, place of residence of women, head occupation, education level of women, work status of women, and head education have a significant effects on time to age at first marriage, while religion of women, access to media of women, and wealth index of a household were not have a significant effect . This seems a similar study with differed statistical analysis from the present study. But, the parameter estimate in the present study is more precise which leads to more accurate statistical inference than its counterpart. This is because the present study assumes a parametric distribution for the baseline hazard function of the model, while this previous study doesn’t assume. In another way, the findings from the present study are more useful for making effective decision relative to this previous study.
Strengths and limitations of this study
The study provides information on marriage in Ethiopian women by analyzing the impact of different covariates on time to age at first marriage using advanced methodology.
The study put a solution for women’s health-related problems in Ethiopia and gives input for designing women health-related policy and practice in Ethiopia.
The study verifies whether there is heterogeneity in time to age at first marriage for women in Ethiopia among regions in Ethiopia.
All low and middle-income countries may not have the same cultural norms, traditional beliefs, and moral rules for the early marriage of girls. In addition, they may not have a similar socio-economic condition. As a result, our study findings may not be conclude-able for those countries.
There may be recall bias when respondents respond to the questionnaire about time-to-age at first marriage.
In this study, we were considered various characteristics of women to analyze their effects on time-to-age at first marriage among women aged 15–49 in Ethiopia by using the Log logistic-Gamma shared frailty model. However, the most significant contributing factors to delaying time to age at first marriage of women aged 15–49 in Ethiopia were increased education level of women, increased education level of the head, increased income, residing in urban and being followers of religion other than orthodox, catholic, protestant & Muslim. In this study, the work status of women, access to mass media, and head occupation were not having a significant effect on time to age at first marriage. The model reveals the presence of a significant cluster (frailty) effect. This implies the heterogeneity of age at first marriage for women aged 15–49 among regions in Ethiopia. The government of Ethiopia and the concerned bodies should revise the women's health policy and practice to reduce early marriage and give attention to women; illiterate, live in rural areas, and have illiterate and poor heads.
Availability of data and materials
After a brief registration and explanation of the data's intended use from the online form of Ethiopian CSA, the data freely available at http://dhsprogram.com/data/dataset/Ethiopia_Standard-DHS_2016.cfm?flag=0
Central statistical agency
Demographic and health survey
Ethiopian demographic and health survey
Csa I. Central Statistical Agency (CSA)[Ethiopia] and ICF. Ethiopia Demographic and Health Survey, Addis Ababa. Ethiopia and Calverton. 2016.
Budu E, Ahinkorah BO, Seidu A-A, Hagan JE, Agbemavi W, Frimpong JB, et al. Child marriage and sexual autonomy among women in Sub-Saharan Africa: evidence from 31 demographic and health surveys. Int J Environ Res Public Health. 2021;18(7):3754.
Liang M, Simelane S, Chalasani S, Snow R. New estimations of child marriage: Evidence from 98 low-and middle-income countries. PLoS ONE. 2021;16(10): e0258378.
Gurmu E, Etana D. Age at first marriage and first birth interval in Ethiopia: analysis of the roles of social and demographic factors. African Popul Stud. 2014;28(3):1332–44.
Lindstrom DP, Kiros G-E, Hogan DP. Transition into first intercourse, marriage, and childbearing among Ethiopian women. Genus. 2009;65(2):45.
Ondenge K, Gvetadze R, Otieno G, Gust DA, McLellan-Lemal E. Factors associated with age of first sex among women screened for an observational contraceptive vaginal ring study in Kisumu, Kenya, 2014. Afr J Reprod Health. 2021;25(1):101–13.
Bezie M, Addisu D. Determinant of early marriage among married women in Injibara town, north west Ethiopia: community - based cross-sectional study. BMC Womens Health. 2019;19(1):1–6.
Alazbih N, Kaya AH, Mengistu M, Gelaye K. Age at first marriage and fertility decline in DABAT health and demographic surveillance system site, Northwest Ethiopia: decomposition analysis. Int J Womens Health. 2021;13:1197–206.
Hougaard P. Analysis of multivariate survival data. Berlin: Springer Science & Business Media; 2012.
Crowther MJ, Look MP, Riley RD. Multilevel mixed effects parametric survival models using adaptive Gauss-Hermite quadrature with application to recurrent events and individual participant data meta-analysis. Stat Med. 2014;33(22):3844–58.
Andersen PK, Gill RD. Cox’s regression model for counting processes: a large sample study. Ann Stat. 1982;1100–20.
Prentice RL, Williams BJ, Peterson AV. On the regression analysis of multivariate failure time data. Biometrika. 1981;68(2):373–9.
Wei L-J, Lin DY, Weissfeld L. Regression analysis of multivariate incomplete failure time data by modeling marginal distributions. J Am Stat Assoc. 1989;84(408):1065–73.
Clayton DG. A model for association in bivariate life tables and its application in epidemiological studies of familial tendency in chronic disease incidence. Biometrika. 1978;65(1):141–51.
Ebrahimi N, Zhang X, Berg A, Shinnar SBT-H of S. Frailty Model and its Application to Seizure Data. In: Advances in survival analysis. Elsevier; 2003. p. 495–517.
Wienke A. Frailty models in survival analysis. Cambridge: CRC Press; 2010.
Ibrahim JG, Chen M-H, Sinha D. Bayesian survival analysis. Berlin: Springer Science & Business Media; 2001.
Yu B. Estimation of shared gamma frailty models by a modified EM algorithm. Comput Stat Data Anal. 2006;50(2):463–74.
Bennett S. Analysis of survival data by the proportional odds model. Stat Med. 1983;2(2):273–7.
Bain LJ. Analysis for the linear failure-rate life-testing distribution. Technometrics. 1974;16(4):551–9.
Ragab A, Green J. On order statistics from thelog-logistic distribution and their properties. Commun Stat Methods. 1984;13(21):2713–24.
Haloi A, Limbu DK. Socio-economic factors influence the age at first marriage of Muslim women of a remote population from North-East India. Sociology. 2013;75:79.
Adebowale SA, Fagbamigbe FA, Okareh TO, Lawal GO. Survival analysis of timing of first marriage among women of reproductive age in Nigeria: regional differences. Afr J Reprod Health. 2012;16(4):95–107.
Von Elm B, Hirschman C. Age at first marriage in Peninsular Malaysia. J Marriage Fam. 1979;877–91.
Zahangir MS, Kamal MM. Several attributes linked with child marriage of females’ in Bangladesh. Int J Stat Syst. 2011;6(1):107–17.
Gobena MG, Berhan YB. Modeling the determinant of time to age at first marriage among women in Ethiopia using cox models with mixed effects. Reprod Health. 2022;19(1):1–6.
We have deepest gratitude to Central Statistical Agency of Ethiopia for making available the data used in this study.
Ethics approval and consent to participate
The data set used in this study was anonymized before we used it for data analysis. From the Ethiopian Central Statistical Agency, approval was sought and volition was given for this study. All methods in this study were carried out in accordance with relevant guidelines and regulations.
Consent for publication
Conflict of interests
The authors declare that as there is no conflict of interest.
About this article
Cite this article
Gobena, M.G., Alemu, Y.M. Analyzing factors associated with time to age at first marriage among women in Ethiopia: log logistic-gamma shared frailty model. BMC Women's Health 22, 191 (2022). https://doi.org/10.1186/s12905-022-01775-1