In this study, we took into account nine predictor variables: age group (5 years), education level of the respondent, husband's education level, wealth index, religion, residence type, desired last child, region, and contraceptive method, which were found to have significant effects on the dependent variable. The dependent variable was the ever-born child of Ethiopian women aged 15–49. This variable had 15 categories (1–15) in the main dataset. The data recorded the number of children born to women between the ages of 15 and 49 who contributed to the population by giving birth, with the lowest possible outcome being one child. Implies, there are no zero count observations, and the data characteristics indicate that this is a count data that is truncated at zero.
The average number of children per mother was 4.04, with a standard deviation of 2.48, indicating that the zero-truncated Poison regression, ZTPR model was inefficient in analyzing this data and that it was forwarded to another model called zero-truncated negative binomial regression, ZTNBR, which was used to treat the dataset well.
Descriptive Statistics Of Response Variable
The summary statistics of the response variable (children ever born) is present in Table 1.The Table 1 shows the summary statistics for the response variable (children ever born). The response variable (children ever born) has a range of 1–15. There were 5753 observations in this study, with the mean and standard deviation of the response variable being 4.04 and 2.48, respectively. The smallest number of children was one, and the largest number of children was fifteen. Furthermore, the summary statistic in Table 1 indicates an overdispersion in the number of children ever born. As a result, this overdispersion property demonstrates the suitability of zero-truncated negative binomial regression models (Mohebbi et al., 2014).
Table 1
Summary statistics for response Variable( Number children ever born)
Variable
|
Obs
|
Mean
|
Std.
|
Min
|
Max
|
NuCEB
|
5753
|
4.04
|
2.48
|
1
|
15
|
Figure 1 shows that the number of children born is consistent with a zero-truncated negative binomial distribution. Thus, rather than using the zero-truncated Poisson regression model, the number of children ever born dataset is estimated by over-dispersed existing models such as the zero-truncated negative binomial regression model.
Table 2
Frequency distribution of response variable (number of children ever born,NuCEB)
Total children ever born, NuCEB *
|
|
age in 5-year groups
|
N
|
%
|
Mean
|
Std
|
Min
|
Max
|
age in 5-year groups
|
15–19
|
296
|
5.15
|
1.42
|
0.66
|
1
|
4
|
20–24
|
1143
|
19.87
|
2.09
|
1.05
|
1
|
6
|
25–29
|
1858
|
32.30
|
3.42
|
1.70
|
1
|
10
|
30–34
|
1239
|
21.54
|
4.95
|
2.06
|
1
|
11
|
35–39
|
765
|
13.30
|
5.91
|
2.33
|
1
|
14
|
40–44
|
339
|
5.89
|
7.27
|
2.29
|
1
|
15
|
45–49
|
113
|
1.96
|
8.41
|
2.84
|
2
|
14
|
Region where mothers live
|
tigray
|
454
|
7.89
|
3.59
|
2.23
|
1
|
11
|
afar
|
652
|
11.33
|
4.07
|
2.35
|
1
|
13
|
amhara
|
511
|
8.88
|
3.75
|
2.32
|
1
|
12
|
oromia
|
719
|
12.50
|
4.41
|
2.74
|
1
|
14
|
somali
|
637
|
11.07
|
5.23
|
2.62
|
1
|
15
|
benishangul
|
530
|
9.21
|
4.25
|
2.53
|
1
|
12
|
snnpr
|
660
|
11.47
|
4.18
|
2.29
|
1
|
12
|
gambela
|
450
|
7.82
|
3.76
|
2.30
|
1
|
11
|
harari
|
447
|
7.77
|
3.75
|
2.32
|
1
|
15
|
addis adaba
|
291
|
5.06
|
2.20
|
1.22
|
1
|
8
|
dire dawa
|
402
|
6.99
|
3.73
|
2.62
|
1
|
14
|
Type of place of residence
|
urban
|
1328
|
23.08
|
3.18
|
2.25
|
1
|
13
|
rural
|
4425
|
76.92
|
4.29
|
2.49
|
1
|
15
|
highest educational level of mothers
|
no education
|
3149
|
54.74
|
4.98
|
2.49
|
1
|
15
|
primary
|
1823
|
31.69
|
3.16
|
2.08
|
1
|
12
|
secondary
|
480
|
8.34
|
2.33
|
1.46
|
1
|
9
|
higher
|
301
|
5.23
|
2.17
|
1.12
|
1
|
8
|
religion
|
orthodox
|
1612
|
28.02
|
3.45
|
2.22
|
1
|
12
|
catholic
|
32
|
0.56
|
3.25
|
1.87
|
1
|
7
|
protestant
|
1053
|
18.30
|
3.92
|
2.46
|
1
|
14
|
muslim
|
2974
|
51.69
|
4.39
|
2.57
|
1
|
15
|
traditional
|
62
|
1.08
|
4.48
|
2.13
|
1
|
11
|
other
|
20
|
0.35
|
4.85
|
2.50
|
1
|
9
|
sex of household head
|
male
|
4598
|
79.92
|
4.09
|
2.47
|
1
|
15
|
female
|
1155
|
20.08
|
3.84
|
2.52
|
1
|
15
|
wealth index combined
|
poorest
|
1964
|
34.14
|
4.69
|
2.50
|
1
|
15
|
poorer
|
994
|
17.28
|
4.42
|
2.56
|
1
|
13
|
middle
|
805
|
13.99
|
4.08
|
2.42
|
1
|
12
|
richer
|
738
|
12.83
|
3.72
|
2.38
|
1
|
15
|
richest
|
1252
|
21.76
|
2.87
|
1.99
|
1
|
13
|
|
Total
|
5753
|
100.00
|
4.04
|
2.48
|
1
|
15
|
Table 2 displays the total number of children born by sociodemographic variables, including counts, percentages, means, standard deviations, minimums, and maximums. It is truly the case that the number of children increases as the age of the mother increases. Mothers in the first age group, 15–19, had the lowest mean number of children ever born (1.46), the highest number of children born (4), and a standard deviation is 0.69. Apart from having the highest mean number of children (7.76, which is higher than the aggregate mean), mothers in the last age group (44–49) have a standard deviation of 2.71, and the maximum number of children ever born is 14. In terms of regional distribution, Somalia (6.67) and Oromia (6.33) had the highest mean number of children born, with standard deviations of 2.78 and 2.85, respectively. Furthermore, Addis Abeba city administration mothers had the lowest mean number of children (2.96) and the lowest standard deviation (1.64).
In terms of place of residence, mothers in the rural area had a higher mean number of children ever born (4.29) than those in the urban area (3.18), with a standard deviation of 2.49 and 2.25, respectively. It is a fact that, as the education level of mothers increases, the number of children born from them decreases. Evidently, mothers with no education level had the largest mean number of children ever born (4.98), whereas mothers with higher education had the smallest mean number of children ever born (2.17), with a corresponding standard deviation of 2.49 and 1.12, respectively. Surprisingly, the maximum number of children ever born from educated mothers is 8; whereas, from no educated mothers, the maximum number of children ever born is 15. The mean number of children ever born (4.09) is high when the sex of the household head is being male as compared to being female (3.84), with a standard deviation of 2.47 and 2.52, respectively.
Besides, the mean number of children ever born from Muslim mothers and traditional believers is 3.39 and 4.48, with a standard deviation of 2.57 and 2.13, respectively. The average number of children born to Protestant and Catholic mothers is 3.92 and 3.25, with standard deviations of 2.46 and 1.87, respectively. Regarding the wealth index, the average number of children born decreases as the wealth index rises from the poorest to the richest level. The mean number of children born when the mother's wealth index is the lowest(poorest) is 4.69, and the mean number of children born when the mother's wealth index is the highest(richest) is 2.87, with standard deviations of 2.50 and 1.99, respectively.
Estimating The Number Of Children Ever Born Using Ztpr And Ztnbr Models
Table 3 reveald the output of the ZTPR and ZTRNB regression models, estimated using Newton-Rapson iteration methods, is shown in Table 4 for the number of children ever born per woman. As a result of these analyses, the log number of children is expected to increase or decrease by approximately the corresponding coefficient in the column of coefficients in the ZTPR model for every one unit change in a unit of the significant predictors (coef.). In this model, variables with p-values less than 0.05 were considered statistically significant. We use incidence rate ratios (to interpret categorical data, which is important for explaining the change in percentage ) of significant predictors.
The childbearing tendency increases for Ethiopian mothers from the lower age group of 20–24 to 45–49. The mothers in the age groups of 20–24 and 25–29 had approximately 1.13 and 2.3 times more children than mothers in the age group of 15–19, respectively. Similarly, mothers aged 30–34 had 3.1 times more children, 35–39 had 3.67 times more children, 40–44 had 4.1 times more children, and 44–49 had 4.87 times more children than the reference age group of mothers. In terms of the regional distribution of the number of children ever born, mothers in Addis Abeba had 13%. In terms of religious affiliation, Protestant and Muslim mothers had 8% and 7% more children, respectively, than Orthodox mothers. Regarding the wealth index of mothers, the richest mothers had a 6% lower number of children ever born than the poorest mothers. fewer children than the Tigray reference group.
It is obvious that as mothers' education levels rise, the number of children born decreases. The result of zero truncated Poisson regression model support revealed that mothers with primary, secondary, and higher education levels had 5%, 15%, and 17% fewer children than mothers without education, respectively. In terms of religious affiliation, Protestant and Muslim mothers had 8% and 7% more children, respectively, than Orthodox mothers. Regarding the wealth index of mothers, the richest mothers had a 4% lower number of children ever born than the poorest mothers.
The number of children ever born increases as the number of living children increases and the number of births within the last five years increases, whereas the number of children ever born decreases as mothers' age at first birth decreases.
Table 3
the result of ZTPR and ZTNBR models
|
ZTPR
|
ZTNBR
|
Children_ever_born
|
IRR
|
SE
|
[95% CI]
|
IRR
|
SE
|
[95% CI]
|
Mother_age_group (15–19)
|
REF.
|
|
|
|
|
|
|
20–24
|
2.127***
|
0.189
|
1.786
|
2.532
|
2.127***
|
0.189
|
1.786
|
2.532
|
25–29
|
3.314***
|
0.289
|
2.794
|
3.931
|
3.314***
|
0.289
|
2.794
|
3.931
|
30–34
|
4.104***
|
0.363
|
3.451
|
4.882
|
4.104***
|
0.363
|
3.451
|
4.882
|
35–39
|
4.674***
|
0.423
|
3.913
|
5.582
|
4.674***
|
0.423
|
3.913
|
5.582
|
40–44
|
5.106***
|
0.476
|
4.253
|
6.130
|
5.106***
|
0.476
|
4.253
|
6.130
|
45–49
|
5.877***
|
0.581
|
4.841
|
7.134
|
5.877***
|
0.581
|
4.841
|
7.134
|
Region (Tigray)
|
REF.
|
Afar
|
0.963
|
0.042
|
0.883
|
1.050
|
0.963
|
0.042
|
0.883
|
1.050
|
Amhara
|
0.955
|
0.035
|
0.889
|
1.026
|
0.955
|
0.035
|
0.889
|
1.026
|
Oromia
|
0.972
|
0.038
|
0.900
|
1.050
|
0.972
|
0.038
|
0.900
|
1.050
|
Somali
|
1.027
|
0.045
|
0.943
|
1.118
|
1.027
|
0.045
|
0.943
|
1.118
|
Benishangul
|
1.030
|
0.041
|
0.953
|
1.114
|
1.030
|
0.041
|
0.953
|
1.114
|
SNNPR
|
1.000
|
0.040
|
0.925
|
1.081
|
1.000
|
0.040
|
0.925
|
1.081
|
Gambela
|
0.999
|
0.045
|
0.915
|
1.090
|
0.999
|
0.045
|
0.915
|
1.090
|
Harari
|
0.998
|
0.045
|
0.913
|
1.091
|
0.998
|
0.045
|
0.913
|
1.091
|
Addis Adaba
|
0.868**
|
0.053
|
0.769
|
0.979
|
0.868**
|
0.053
|
0.769
|
0.979
|
Dire Dawa
|
1.021
|
0.048
|
0.931
|
1.118
|
1.021
|
0.048
|
0.931
|
1.118
|
Residence (Urban)
|
REF.
|
Rural
|
1.007
|
0.026
|
0.958
|
1.058
|
1.007
|
0.026
|
0.958
|
1.058
|
Morher_Edu_level (No education)
|
REF.
|
Primary
|
0.953***
|
0.018
|
0.919
|
0.989
|
0.953***
|
0.018
|
0.919
|
0.989
|
Secondary
|
0.850***
|
0.034
|
0.785
|
0.920
|
0.85***
|
0.034
|
0.785
|
0.920
|
Higher
|
0.827***
|
0.045
|
0.744
|
0.920
|
0.827***
|
0.045
|
0.744
|
0.920
|
Religion(Orthodox)
|
REF.
|
Catholic
|
1.031
|
0.035
|
0.964
|
1.102
|
1.031
|
0.035
|
0.964
|
1.102
|
Protestant
|
1.075***
|
0.012
|
1.052
|
1.100
|
1.075***
|
0.012
|
1.052
|
1.100
|
Muslim
|
1.071***
|
0.011
|
1.049
|
1.094
|
1.071***
|
0.011
|
1.049
|
1.094
|
Traditional
|
1.009
|
0.029
|
0.953
|
1.069
|
1.009
|
0.029
|
0.953
|
1.069
|
Other
|
1.184***
|
0.055
|
1.081
|
1.296
|
1.184***
|
0.055
|
1.081
|
1.296
|
Sex_HH(Male)
|
REF.
|
Female
|
1.014
|
0.020
|
0.977
|
1.053
|
1.014
|
0.020
|
0.977
|
1.053
|
Wealth(Poorest)
|
REF.
|
Poorer
|
1.011
|
0.022
|
0.969
|
1.054
|
1.011
|
0.022
|
0.969
|
1.054
|
Middle
|
1.002
|
0.024
|
0.956
|
1.050
|
1.002
|
0.024
|
0.956
|
1.050
|
Richer
|
0.973
|
0.025
|
0.925
|
1.023
|
0.973
|
0.025
|
0.925
|
1.023
|
Richest
|
0.935**
|
0.030
|
0.877
|
0.997
|
0.935**
|
0.030
|
0.877
|
0.997
|
Living_children
|
1.142***
|
0.006
|
1.130
|
1.153
|
1.142***
|
0.006
|
1.130
|
1.153
|
Birthin_5_year
|
1.176***
|
0.012
|
1.152
|
1.200
|
1.176***
|
0.012
|
1.152
|
1.200
|
Age_at_first_birth
|
0.962***
|
0.002
|
0.958
|
0.966
|
0.962***
|
0.002
|
0.958
|
0.966
|
_cons
|
0.972
|
0.081
|
0.825
|
1.146
|
0.972
|
0.081
|
0.825
|
1.146
|
lnalpha
|
|
|
|
|
-17.419
|
19.409
|
-55.460
|
20.622
|
*** Significant At 0.01, ** Significant At 0.05, *significant At 0.1
Table 4
Model comparison using AIC and BIC
Model
|
N
|
ll(null)
|
ll(model)
|
df
|
AIC
|
BIC
|
ZTPR
|
5,753
|
-13055.1
|
-8750.103
|
34
|
17568.21
|
17794.56
|
ZTNBR
|
5,753
|
-12596.2
|
-8750.103
|
35
|
17570.21
|
17803.22
|
Table 4 displays the values of the well-known model selection criteria AIC, BIC, and log (L). Although both models produced nearly identical results, their standard errors were slightly different(Karimuzzaman et al., 2020). The AICs, BICs, and Log-Likelihoods for both models were computed and compared to determine the best fitted model. For the zero-truncated Poisson and negative binomial models, the AICs were 17568.21 and 17570.21, respectively. For the zero-truncated Poisson and Negative Binomial models, the BICs were 17794.56 and 17803.22, respectively. Based on AIC and BIC criteria, the zero truncated Poisson regression (ZTPR) model is preferable to the candidate model, the zero truncated negative binomial regression (ZTNBR). Furthermore, the overdispersion parameter delta (\(\delta\)) is not statistically significant, which means overdispersion is not a problem in the dataset. Therefore, this shows that a zero truncated Poisson regression model is appropriate to estimate the numbers of children ever born.