ALL Metrics
-
Views
-
Downloads
Get PDF
Get XML
Cite
Export
Track
Research Article

Determining the lx and estimating the force of mortality for children under five in Bangladesh

[version 1; peer review: 2 approved with reservations]
PUBLISHED 03 Dec 2019
Author details Author details
OPEN PEER REVIEW
REVIEWER STATUS

Abstract

Background: Under-five mortality is a health indicator in population health and the health sciences. Force of death is a more accurate measure of mortality than others, which indicates the situation of mortality when time tends to zero. The purpose of this research is to construct a simulation for lx (the number of people living at exact age x) for under-five mortality in Bangladesh as a whole, and separately for rural and urban areas, and then estimating the mortality rate in regard to these matched models.
Methods: Secondary data were collected from Bangladesh Demographic and Health Survey 2007. A polynomial model was selected to match the lx values. To prove the accuracy of the model, the rule of cross-validity was applied.
Results: It has been shown that fourth degree polynomial models can be adjusted to lx values ​​with explanation of more than 94% variation. It was noted that the mortality rate shows a rapidly decreasing pattern for people aged 0-20 months, a monotonically increasing pattern for those aged 20-53.5 months and then it begins to decrease. It is found that the mortality rate in rural areas is higher than in urban areas of all ages.
Conclusions: It has been shown that the mortality rate in rural areas is higher than for urban areas of all ages, except for those aged 53.5 months. The health situation should therefore be improved to reduce mortality in rural Bangladesh.

Keywords

Children under five years, Force of mortality, Urban and rural mortality, F-test, Bangladesh

Abbreviations

ANOVA, Analysis of variance; BDHS, Bangladesh Demographic and Health Survey; BBS, Bangladesh Bureau of Statistics; CVPP, cross validity prediction power; IMR, infant mortality rates; EA, enumeration area; MDG 4: Millennium Development Goal 4; NGO, non-governmental organization

Introduction

The mortality rate of infants is a reflection of the well-being and nutritional status of children, and also indicates the socio-economic progress in a given country or region. In truth, it reflects the well-being of mothers, the availability and standard of care for well-being, and how easy it is to receive support in the community and/or country. Child mortality rates and, in particular, infant mortality rates (IMRs) dropped sharply in developed countries, but over the past five decades it has been difficult to lower rates in developing countries such as Bangladesh. It has been found that 29,000 children under the age of five die every day in Bangladesh1.

Five years of previous mortality studies of children under five (2010–2014) is 46 deaths per 1000 live births. Bangladesh achieved its Millennium Development Goal 4 (MDG 4) target for deaths under the age of five, decreasing the rate of child mortality to 48 per 1000 births by 2015. The current infant mortality rate is 38 deaths per 1000 live births, and the mortality rate for children is 8 per 1000 children. In childhood, the risk of death in the first month of life (28 deaths per 1000 live births) is almost three times higher than in the next 11 months (10 deaths per 1000 live births). It has also been found that mortality in the neonatal stage constitutes around 61% of all deaths under the age of five2.

It was noted that in 2007, around 9.2 million live births resulted in death of the child before their fifth birthday3, but in the next twelve months this number fell to 8.8 million4,5, with 41% of these deaths occurring in newborns5. The IMRs are 44.5 and 76 per 1,000 live births universal for more developed and less developed countries, respectively, indicating that infant death is still high in some vulnerable groups and areas, despite declining trends6. The highest under-five death rates are still observed in sub-Saharan Africa, while South Asia has the second highest rates in the world5,7. Many countries have not reached MDG 45, but Bangladesh is currently one of the few countries in the world, particularly in South Asia and sub-Saharan Africa, to have reached this goal. Bangladesh decreased death rates from 151 deaths per 1,000 live births in 1990 to 50 deaths per 1,000 live births in 20158. Current studies have revealed that the mortality risk index is 45 per 1,000 live births6 and that under-five mortality is between 719 and 6510 per 1,000 live births. Although the recent downward trend in infant mortality is notable11, it remains high in Bangladesh due to the high incidence of malnutrition and childhood diseases. One child dies every three or four minutes, resulting in 14 child deaths every hour in Bangladesh12. Deaths usually occur during the neonatal period13.

For this reason, a number of studies have been carried out on the mortality of infants and children and under five; however, in Bangladesh, there have been few studies investigating the mortality rate of infants and children under the age of five. The polynomial model can be used to approximate a complex nonlinear relationship, which is related to the polynomial degree, which depends on the data pattern. In addition, the response variable is estimated at any time within the domain of the matched model in terms of the explanatory variable, i.e. the respondent's age. Furthermore, mortality is important because it is considered a mortal situation when time tends to zero, but age-specific mortality is the mortality rate within a specific age group. Therefore, the main purpose of this article is to match some models to lx values for children under five in Bangladesh as a whole, as well as in rural and urban areas separately.

Methods

Data source

To achieve the above-mentioned goals of this study, secondary data on the lx value for age (in months), i.e. the group data for children under five in Bangladesh as a whole and separately for rural and urban areas, are from Bangladesh Demographic and Health Survey (BDHS) 200714. The 2007 BDHS was representative at national level and included the entire population that did not live in institutional housing units in the country. The survey was used as a sampling frame for the list of enumeration areas (EAs) prepared for the Population and Housing Census 2007, provided by Bangladesh Statistics Office (BBS). The BDHS study was conducted using a two-stage stratified sample of households. A total of 600 areas were initially selected, with 207 clusters in urban areas and 393 in rural areas. A complete list of households was then compiled in all selected areas to provide a sampling framework for the second phase of household selection. In the second phase of the sampling, a systematic sample of an average of 30 households per area was selected to provide statistically reliable estimates of the most important demographic and health variables for the country as a whole, for urban and rural areas individually, and for each of the seven divisions (Dhaka, Rajshahi, Rangpur, Chittagong, Khulna, Barisal and Sylhet). The total sample size is 8,721, of which 5,840 come from rural areas and 2,881 from urban areas during the ten years preceding the survey, in which the sample consists of eligible women, i.e. women who have ever been married and who are aged between 10 and 49 with at least one child. Their data are used as the raw data for this study.

Smoothing of age data

To better adjust the model, the age (in months) associated with the lx values in Bangladesh and rural and urban areas was smoothed separately using the 4253H smoothing method twice15. The smoothing method was done using Minitab 12.1. It should be noted that the smoothing method is performed to remove age as well as reporting errors before model fitting for better performance as well as better prediction. There are various methods of smoothing, such as the free-hand curve method, the moving average method and 4253H, the double method. The free-hand curve method does not have a mathematical background and the moving average method loses some observations. For this reason, the 4253H method is used, a double method for better efficacy that does not lose any observations.

Model fitting

The lx values for under-five mortality showed a non-linear distribution. Therefore, an n degree polynomial model is identified in this case and the model is given as:

y=a0+i=1naixi+u(1)

in which x represents the age in months that is used as explanatory variable in the model, y represents the lx value, a0 is a constant, ai is the coefficient of the xi (i = 1, 2, 3, ..., n) and u is the error term of the model16. Here, a suitable n is chosen for which the sum of the squares of the error is the lowest. These models were built using SPSS 22.

Model validation procedure

Cross-validation is a technique used in model selection to better estimate the test error of a predictive model. The idea behind cross-validation is to create a better assessment of a model’s predictive performance. For this, to check the legitimacy of these models over the population, the cross validity prediction power (CVPP), ρcv2, is employed at this juncture. The mathematical formula for the CVPP is specified below:

ρcv2=1(n1)(n2)(n+1)n(nk1)(nk2)(1R2)(2)

where n is the quantity of classes, k is the number of explanatory variables in the fitted model and R is the correlation between the observed and the fitted values of the response variable. The shrinkage coefficient of the fitted model is the absolute value of (ρcv2R2); where ρcv2 is the CVPP and R2 is the proportion of variation of this model. Furthermore, the stability situation of R2 of this model is (1 – shrinkage). The estimated CVPP corresponding to R2 and the results for the model fittings are summarized in Table 1. It has been shown that the CVPP can also be used as a model validation technique1726.

Table 1. The results of CVPP and information on model fittings.

ModelsnKParametersp-valuesR2ρcv2ShrinkageF-test
(calculated
and tabulated)
Model for
Bangladesh
115a0
a1
a2
a3
a4
<0.0001
0.001102
0.007429
0.020427
0.038182
0.945530.8217310.12379826.08
and
9.15 with (4,6) df
at 1% level of significance.
Significant
Model for
rural areas
115a0
a1
a2
a3
a4
<0.0001
0.001261
0.008829
0.023974
0.044073
0.94754 0.8283090.1192327.093
and
9.15 with (4,6) df
at 1% level of significance.
Significant
Model for
urban areas
115a0
a1
a2
a3
a4
<0.0001
0.000799
0.005044
0.014332
0.027845
0.94241 0.8115370.1308824.548
and
9.15 with (4,6) df
at 1% level of significance.
Significant

CVPP, cross validity prediction power; df, degrees of freedom.

The F-test

To identify the overall significance level of the formulated model and the significance of R2 of the model, the F-test is used in this paper. The formula for the F-test is given as:

F=R2(p1)(1R2)(np)with(p1,np)degreesoffreedom(df).

where p is the number of parameters of the fitted model, n is the number of cases and R2 is the coefficient of determination of the fitted model27.

Instantaneous force of mortality (μx)

The instantaneous force of mortality at age x is defined as the ratio of the instantaneous rate of decrease in lx to the value of lx. It is denoted by μx and, in the limiting case, it is given using differential calculus by:

μx=Ltt0lx+tlxt.lx=1lxddx(lx)=ddx(lnlx). [28, 29]

In brief, it is also called the force of mortality. In fact, the force of mortality is a more accurate measure in mortality studies than age-specific death rates. It measures the mortality situation when time tends to zero. This is an instantaneous measurement, rather than an interval measurement.

Results

The polynomial model is used to fit the age associated with the lx values for children under five in Bangladesh in rural and urban areas, and these fitted models are presented (Figure 1, Figure 2, Figure 3) below:

  • i) For Bangladesh, y=950.1666-39.3148x+2.1053x2-0.0449x3+0.0003x4; t-value: (42.33616) (-5.84934) (3.96269) (-3.12604) (2.64699)

  • ii) For rural areas, y=966.3290-26.1108x+1.3809x2-0.0294x3+0.0002x4; t-value: (63.16859) (-5.69949) (3.81335) (-3.00112) (2.54014)

  • iii) For urban areas, y=984.0909-13.2389x+0.7251x2-0.0155x3+0.0001 x4; t-value: (138.4389) (-6.2189) (4.3088) (-3.4093) (2.8858)

696ac7af-df86-43dc-8482-70297d33ea09_figure1.gif

Figure 1. Smoothed and predicted lx values for children under five in Bangladesh.

X axis represents age group and Y axis represents lx.

696ac7af-df86-43dc-8482-70297d33ea09_figure2.gif

Figure 2. Smoothed and predicted lx values for children under five in rural areas in Bangladesh.

X axis represents age group and Y axis represents lx.

696ac7af-df86-43dc-8482-70297d33ea09_figure3.gif

Figure 3. Smoothed and predicted lx values for children under five in urban areas in Bangladesh.

X axis represents age group and Y axis represents lx.

The results of the fitted models and the CVPP corresponding to the R2 of these fitted models of lx values for children under five in Bangladesh and in the rural and urban areas are presented in Table 2. It seems that the constructed models (i) to (iii) are cross-validated, and their shrinkages are also shown in Table 2. These imply that the fitted models (i) to (iii) are more than 81% stable, and this is demonstrated in Table 2.

Table 2. The observed, smoothed and predicted survival function (lx values) for children under five in Bangladesh.

Age in
Months
BangladeshRuralUrban
ObservedSmoothedPredictedObservedSmoothedPredictedObservedSmoothedPredicted
0100099295010009959661000997984
1.5835871896888913930947958966
4787791824855860882932932942
7751752764834834842917918922
10733733726822822816911912910
14.5725723700815813798910909902
20.5714713700805804796909909904
31.5702701717794794805908907912
41.5686688699784784791902904908
53.5677678657775776762902902895
60669669679770770776899899903

In addition, it can be seen that the parameters of these models are significant, with more than 94% of the explained variance, which is also shown in Table 3. The stability position for R2 of these matched models is over 87%. The calculated values for the F-test, i.e. the analysis of variance (ANOVA) for these matched models with (5, 7) df, are 26.08, 27.09 and 24.548 and are displayed in Table 2, while the corresponding tabular values are 9.15 for models (i) to (iii) at a significance level of 1%. Consequently, on the basis of these test statistics, it is concluded that the F-test is highly significant, and therefore, that these models are highly significant. That is why the fitting of the models is good. Hence, forecasts based on these models are shown in Table 1.

Table 3. Instantaneous force of mortality for children under five in Bangladesh, in rural and in urban areas.

Age in
Months
BangladeshRuralUrban
00.0413770.0270210.013453
1.50.0371720.0238280.011562
40.0298020.0186180.008661
70.0209940.0128620.005672
100.0130580.0079810.003283
14.50.0041750.0027190.000851
20.50.0010500.0004320.000426
31.50.0039450.0020200.001310
41.50.0154480.0078610.005014
53.50.0241070.0108860.008379
600.0182650.0065990.008007

It should be mentioned here that other more typical models, such as the linear model (R2 = 0.53), exponential (R2 = 0.55), square (R2 = 0.73) and the cubic model (R2 = 0.85), were also used for the Bangladesh data set, but they did not show a better fit due to the low coefficient of determination.

The ages associated with force of mortality for children aged under five in Bangladesh, and in rural and urban areas, are estimated from the fitted model of the lx values that are presented in Table 3 and shown in Figure 4 to Figure 6.

696ac7af-df86-43dc-8482-70297d33ea09_figure4.gif

Figure 4. Instantaneous force of mortality for children under five in Bangladesh.

X axis represents age in months and Y axis represents force of mortality.

696ac7af-df86-43dc-8482-70297d33ea09_figure5.gif

Figure 5. Instantaneous force of mortality for children under five in rural areas in Bangladesh.

X axis represents age in months and Y axis represents force of mortality.

696ac7af-df86-43dc-8482-70297d33ea09_figure6.gif

Figure 6. Instantaneous force of mortality for children under five in urban areas in Bangladesh.

X axis represents age in months and Y axis represents force of mortality.

In Table 3, as well as in the numbers described above, it can be seen that the mortality of children under five in Bangladesh, as well as in rural and urban areas, shows a rapidly decreasing pattern at 0–20 months, followed by a monotonically increasing pattern for age from 20 to 53.5 months and then it begins to decrease. In addition, it has been demonstrated that the mortality rate in rural areas is higher than for urban areas of all ages excluding those aged 53.5 months. All three cases show a similar pattern.

Discussion

Explanation of the empirical age patterns of death using these models is one of the most important topics in population studies, especially in demographics. In this study, the fourth degree polynomial models are well-suited to lx values for children ages under five, with an explanation for the huge proportion of variation. The mortality of infants and children in Bangladesh is significantly reduced, but the future increase in life expectancy requires an additional reduction in mortality30. This study estimates how age is related to the mortality of children under the age of five in Bangladesh as a whole and in rural and urban areas separately. We expect these recent mortality findings to convince the government, non-governmental organizations and those responsible for policy planning of the need for a planned increase in socio-demographic development and health care facilities. For this reason, the survival (lx) function taken from the data is selected as the raw material in this study. The mortality force is the ratio between the instantaneous rate of decrease in lx and the value of lx, which is obtained from the functional dependence between lx and age x. Mortality rate should be treated as the probability of death at a given moment, taking into account survival until that time. It was reported that the lx values for the Bangladeshi male population are in line with the four-parameter cubic polynomial model30 and it has been observed that the lx values for the Bangladeshi population are consistent with the third degree polynomial model31. However, in this study, efforts were made to focus on which mathematical model best suits the relationship between age and lx values for minors in Bangladesh as a whole, and in rural and urban areas separately. For this purpose, the polynomial model is chosen. It should be mentioned that a number of studies were carried out using the polynomial model3235. One of the limitations of the survey is that, except for a few questions, almost all information collected in BDHS surveys is subject to reporting errors.

Conclusions

In this study, it was found that the relationship between age and lx values for children under five in Bangladesh results from a four-level polynomial model containing five parameters. The mortality rate shows a rapidly decreasing pattern at 0–20 months and a monotonically increasing pattern between 20 and 53.5 months, and then begins to decrease in all cases. This research result has fundamental strategic implications, especially for educational and awareness program requirements for a real decrease in child mortality and compliance with public health interventions. The researcher recommends additional research, taking into account these unobserved issues that are probably associated with the death of infants and children, to better understand the relationship between family and public factors and the death of a child in Bangladesh. Intervention programs and policies should focus on this research to improve the health consequences for children, in order to achieve future improvement of the situation in Bangladesh as well as in other areas such as sub-Saharan Africa.

Ethical approval

Ethical approval for this study was not applicable, since ethical approval for the collection of data was previously approved for BDHS.

Data availability

Underlying data

The data from BDHS 2007 are free to access (https://dhsprogram.com/data/dataset/Bangladesh_Standard-DHS_2007.cfm?flag=0); however, before you can download data, users must register as a DHS data user for reasons laid out on the DHS website (https://www.dhsprogram.com/data/Registration-Rationale.cfm). Dataset access is only granted for legitimate research purposes (https://dhsprogram.com/data/new-user-registration.cfm).

Comments on this article Comments (0)

Version 1
VERSION 1 PUBLISHED 03 Dec 2019
Comment
Author details Author details
Competing interests
Grant information
Copyright
Download
 
Export To
metrics
Views Downloads
F1000Research - -
PubMed Central
Data from PMC are received and updated monthly.
- -
Citations
CITE
how to cite this article
Islam MR, Rana MS, Hossain MM et al. Determining the lx and estimating the force of mortality for children under five in Bangladesh [version 1; peer review: 2 approved with reservations] F1000Research 2019, 8:2045 (https://doi.org/10.12688/f1000research.21202.1)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
track
receive updates on this article
Track an article to receive email alerts on any updates to this article.

Open Peer Review

Current Reviewer Status: ?
Key to Reviewer Statuses VIEW
ApprovedThe paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approvedFundamental flaws in the paper seriously undermine the findings and conclusions
Version 1
VERSION 1
PUBLISHED 03 Dec 2019
Views
9
Cite
Reviewer Report 03 Aug 2020
Mahmudur Rahman, Consultant, International Centre for Diarrhoeal Disease Research, Bangladesh (icddr, b), Dhaka, Bangladesh 
Approved with Reservations
VIEWS 9
Determining the lx and estimating the force of mortality for children under five in Bangladesh
 
Review Comments:
 
Introduction:
  • Introductory paragraphs focus on IMR whereas this analysis is on Child mortality.
... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Rahman M. Reviewer Report For: Determining the lx and estimating the force of mortality for children under five in Bangladesh [version 1; peer review: 2 approved with reservations]. F1000Research 2019, 8:2045 (https://doi.org/10.5256/f1000research.23341.r66568)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
Views
7
Cite
Reviewer Report 03 Jul 2020
Duah Dwomoh, Biostatistics, School of Public Health, University of Ghana College of Health Sciences, Accra, Ghana 
Approved with Reservations
VIEWS 7
Abstract:
  • The authors indicated in the results section that mortality rate was monotonically increasing for those aged 20-53.5 months but in the conclusion it is restricted to only those aged 53.5 months. Authors should revise that
... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Dwomoh D. Reviewer Report For: Determining the lx and estimating the force of mortality for children under five in Bangladesh [version 1; peer review: 2 approved with reservations]. F1000Research 2019, 8:2045 (https://doi.org/10.5256/f1000research.23341.r65215)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.

Comments on this article Comments (0)

Version 1
VERSION 1 PUBLISHED 03 Dec 2019
Comment
Alongside their report, reviewers assign a status to the article:
Approved - the paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations - A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approved - fundamental flaws in the paper seriously undermine the findings and conclusions
Sign In
If you've forgotten your password, please enter your email address below and we'll send you instructions on how to reset your password.

The email address should be the one you originally registered with F1000.

Email address not valid, please try again

You registered with F1000 via Google, so we cannot reset your password.

To sign in, please click here.

If you still need help with your Google account password, please click here.

You registered with F1000 via Facebook, so we cannot reset your password.

To sign in, please click here.

If you still need help with your Facebook account password, please click here.

Code not correct, please try again
Email us for further assistance.
Server error, please try again.