Model for Assessment of Child Mortality under Different Parity: A Bayesian Swatch

In this study a model is conceptualized to measure the child mortality under different parity of women such that a better strategy can be formulated to bring down mortality rates. In the estimation of probability of child mortality some socio demographic variables are taken in consideration. The estimates are obtained under Bayesian procedure. Two different models are formulated for it and model fitting is observed by graphical approach along with the chi square test. First model is beta-binomial and second is binomial regression model. Second model shows the better fit on the data. The estimate of probability of child mortality at higher parities namely, parity 3, parity 4 and parity 5 were obtained as 0.06, 0.09 and 0.13 respectively on the basis of the second model.


Introduction
Children are important assets of a nation and reduction in infant and child mortality is one of the key factors which are taken into consideration while planning for the growth of the nation. Child mortality is also a good indicator of level and quality of health care as well as socio-economic condition of community. Recognizing the important role that child health plays in the overall health of societies, the Millennium Development Goals include a goal explicitly aimed at reducing child mortality by two-thirds, between 1990 and 2015.
Child mortality, also known as under-5 mortality, refers to the death of infants and children under the age of five. The trend in child mortality is declining over time but rate of decline is not fast. Child Mortality rate is the highest in low-income countries. Still many deaths in the third world go unnoticed since many poor families cannot afford to register their babies in the government registry.
Under five year mortality rate for the world was 91 per 1000 live birth in 1990 which decreased to 43 in 2015 (http://www.who.int/). The under five mortality for India was 126 in 1990 as compared to 48 in 2015 (http://www.childmortality.org/).
No doubt there is significant reduction in child mortality in India but as per UNICEF report, India is still in the list of high child mortality countries.
The infant and child mortality studies have long been of interest to demographers and person concerned with public health problem. The most common problem in such studies is the error in data of deaths during infancy and childhood. In this condition, development of stochastic models is a good choice to minimize the effect of such errors.
Perhaps the first attempt to model the infant mortality was made by Keyfitz [1], after that many authors contributed towards modeling of infant and child mortality [2] [3] [4] [5] [6]. The estimation of model parameters in above mentioned studies were based on classical estimation techniques. The classical estimation technique serves a lot of properties and requires less computational skill but it is unable to find estimates under complex modeling situations. On the other hand, Bayesian estimation approach is capable of doing estimation in complicated situations and provides direct interpretation of results.
The availability of computers based computational technique; Bayesian method of estimation is getting popularity as well as acceptance in all modeling situations. Various authors have been used this technique in the estimation of infant and child mortality [7] [8] [9] [10] [11].
The key features of Bayesian modeling are the likelihood function, which reflects information about the parameters contained in the data, and the prior distribution, which quantifies what, is known about the parameters before observing data. The prior distribution and likelihood are further combined to get the posterior distribution of parameter, which represents total knowledge about the parameters after the data have been observed. Simple summaries of this distribution can be used to isolate quantities of interest and ultimately to draw substantive conclusions whereas the classical modelling approaches uses only the likelihood function and hence sampling distribution of estimator, which reflects information about the parameters contained in the data.
The present study is an attempt to model child death in such a way that we get estimate of child mortality according to the parity of mothers. The aim behind it was to provide a platform for monitoring agencies to re-plan their strategies for bringing faster decline in child mortality. In present study an attempt is also made to find out the relationship of different socio-economic variables with child mortality in total reproductive life span of women according to different parity. Two models have been taken into account to estimate the mortality rate

Method
Let us consider the population of N women who have given birth to n children during their reproductive period. Define a random variable Z ij (i = 1, 2,  , N and j = 1, 2,  , n) as given below: with possible values 0, 1, 2, •••, n and it shows the number of child death experienced by i th women. If we assume that deaths of children experienced by i th women are independent and with same unknown probability say p i then i Y assumed to follow binomial distribution. Here n will denote parity (total children born in life span) of a women and let pi is probability of experiencing child death by i th women of parity n.

Model-I
The model -I is defined as below according to above explanation of the study: Let here n denotes parity of women and we have considered n = 3, 4 and 5, i denotes i th women of parity n (i − 1, 2,  , N). The prior distribution for parameter i p is taken as beta distribution with probability density function: Combining it by Bayes theorem the posterior probability density function will have a beta distribution. We use non-informative beta (1, 1) prior which is flat and reflect no prior information regarding the parameter being estimated.

Model-II
In this model, while estimating the probability of child death to particular women the consideration of different socio-economic and demographic variables was incorporated in model. For each n (n = 3, 4, 5) pi is unknown parameter of interest. Since p i is affected by various socio-economic variables. Then we consider pi as a function of a few explanatory variables X m , (m = 1,  , 4) and propose to use generalized linear model(GLM). The general linear model (logit model) is popularly used to see the association between Bernoulli response variable and the factors affecting it. In modeling by general linear model, the dependent variable is transformed into continuous form by using link function.
Dichotomous dependent variable like Z ij with probability of death of j th (i = 1, 2,  , n) child born to i th (i = 1, 2,  , N) women as p i , a link function called logit link function is used to make dependent variable continuous. It is actually expected value of logarithm of odds of experiencing a child death to not experiencing a child death. Thus GLM defines the relationship as Thus expected value of ( ) and if we take So, where, 1 2 , , , p X X X  (generally called explanatory variables) are the variables associated with ( ) log it i p known as logistic regression. The unknown regression coefficients in (4) are 0 1 2 , , , , p β β β β  which are estimated from the data and for the Bayesian approach they are estimated from their joint posterior distribution. The proposed model is linear for the regression coefficients on the logit scale.
In present study four explanatory variables are taken for the estimation of probability of child death experienced by women. The four explanatory variables are age at first birth (X 1 ), education of mother (X 2 ), religion (X 3 ) and type of house (X 4 ). The variable X 1 quantitative and continuous in nature and measured in years, and X 2 , X 3 , X 4 are categorical variables. X 2 has four levels as illiterate, primary, middle and above middle. There are two levels of X 3 as Hindu and non-Hindu. Finally the fourth explanatory variable X 4 has three levels as kaccha, semi pucca and pucca. Thus the proposed model is: where, Further, for the Bayesian analysis, priors are required for all the parameters before data in hand.
In present study independent non-informative priors have been used considering that we have no strong a priori idea about regression parameters and hence following prior distribution is considered: , , , p y β β β β and thus the joint posterior distribution can be obtained as: After obtaining the posterior distribution of regression coefficients we obtain the Bayesian estimate of regression coefficients from this distribution under squared error loss function with use of Monte Carlo Markov Chain process.
Here Gibbs sampling procedure is used for this purpose.

Data
To check the suitability of proposed methodology in estimating the situation of child mortality, data was taken from District Level Household and Facility Survey (DLHS-3, 2007-2008) [12]. A multi-stage stratified systematic sampling design was adopted in DLHS-3. In each district of country, 50 Primary Sampling Units (PSUs) which were census villages for rural areas and wards for urban ar-   Table 2 describes the distribution of women according to different explanatory variables taken in consideration for the estimation association of child death. Table 2 showed that out of 6688 women of parity 3, 3429 were illiterate; however, among literate women 939, 1017 and 1303 were having education up to primary, middle and above middle respectively. Among the women of parity 3, 5766 were Hindu and 1592, 3117 and 1979 were found to having kaccha, semi pucca and pucca house respectively. At parity 4, 4064 were illiterate, 898 attain primary, 675 were having education up to middle and 609 women had higher education. Out of all women of parity 4, 5109 were Hindu while 1137 were non-Hindu, further 1720 women were living in kaccha house, 3089 women in semi pucca and 1437 women were living in pucca house. Among the women of parity 5, 3946, 669, 452 and 273 were education level as illiterate, primary, middle and above middle respectively. 4184 women were Hindu while 1156 women were non-Hindu. In context of residential accommodation, 1600, 2750 and 990 women found to reside in kaccha, semipucca and pucca houses respectively. The distribution in Table 2 showed a decreasing trend in level of education when the parity of women increases. Most of the women found to reside in semipucca houses and belong to Hindu community.

Analysis
The Bayesian analysis of both of the proposed models was performed on the basis of data given in Table 1 and Table 2 for each parity separately on Win-BUGS software. To obtain the Bayes estimates of parameters of the models described in Equation (6) a programme was written in the language of WinBUGS.
After running this programme in WinBUGS we summarize the estimates after discarding 20,000 initial updates. Thus 40,000 updates were run after the initial burn in and MC error for each parameter was less than 5% of its standard deviation. In addition to controlled MC error, the Kernel density of each parameter was also found to be well in normal shape. After being confirmed with all the diagnostic about convergence of the parameters, the estimates were obtained.

Results
The results obtained from proposed model-I was tested for its fit by the help of graphs of observed frequency of death and expected frequency of death for each parity. Figure 1 shows the fit of parity 3, parity 4 and parity 5 respectively. It is clear that model-I is unable to capture the uncertainty among deaths of children of a mother of each considered parity. Apart from graphical fitting, an attempt was also made to check the fitting by calculating chi-square values for considered parity. The results are shown in Table 3. Table 3 Table 4. Table 4 Table 5, Table 6 and Table 7 for considered parity 3, parity 4 and parity 5 respectively.      Table 5 shows that the regression coefficient of AFB is negative for parity 3 but this value is very small (−0.00783), and HPD interval nearly covers equal area on both sides of 0 value (−0.1306, 0.1148). For the variable mother education all the three lower categories' regression coefficient are positive when the highest education category of mother is taken as reference category. The odd of child mortality is highest in illiterate category (1.47) and the 90% HPD interval for this regression coefficient comes under between 0 and above. The odds ratio for non Hindu is 0.62 as compared to Hindu which shows that Hindu has higher risk of child death as compared to other religions which are prevailing in Uttar Pradesh. The variable type of house does not give any clear pattern in determining child mortality. Table 6 describes the regression analysis for parity 4. The continuous variable AFB shows the positive regression coefficient with 90% HPD interval in positive. The odd of child mortality is higher by 1.12 times in unit increase in AFB as compared to mean level of AFB. Education status of mother shows that with decreasing of mother educational category the risk of child mortality becomes higher. But this result is not so conclusive because the magnitude and variation of regression coefficients are not strong. Again in parity 4 Hindu has higher risk of child mortality nearly in same magnitude as present in parity 3. The regression coefficient for kaccha (0.099) and semi-pucca (0.054) houses are positive though their 90% HPD intervals do not support these finding for conclusive statement.
For parity 5, the outcome of regression analysis is mentioned in Table 7. The table shows there is no role of AFB in child mortality determination. The negative association was found between mother education and child mortality which finding also exist in parity 3 and parity 4 as well. The finding for Hindu females remains same and having higher chance of experiencing child mortality. It is also see by table that the female belonging to pucca house has lesser chance of child mortality as compared to kachha and semipucca house.

Discussions
The aim of present study was to find out valid estimate of relationship between child mortality and selected socio-economic variable through mathematical modelling. Generalized linear model is used in Bayesian setup for this purpose, where risk of death was dependent on selected socio-economic factors. Data were used from DLHS-3 for UP state. Factors like education of mother were found decreasing effect on the risk of child death which is not significant but various studies have supported a direct causal relationship between mother's education and child mortality [3] [13]. This is due to the division of data on the basis of parity or due to consideration of total reproductive life span of women.
The factor religion shows significant increasing risk of child death in Hindus as compared to other religions. This result is quite hesitating to explain and no one can give a strong argument that, why this result occurred but it gave implication that Hindus are at higher risk of child mortality in comparison of other religion so a separate plan for Hindu's family are essential for reducing child mortality and it is also a point to think that the majority of this state is still have higher risk of child death. Although this finding supported by Bhalotra [14]. Bhalotra explored this hypothesis by using information reported by mothers on ideal family composition. They recognize, however, that adding this variable in their analyses does not shed more light on the Muslim mortality advantage. Nonetheless, they conclude that this may be due to data limitations and argue that this is a dimension of the Muslim mortality paradox that is worth exploring further.
Geruso and Spears [15] also conclude "In India, Muslims face significantly lower child mortality rates than Hindus, despite Muslim parents being poorer and less educated on average. Because observable characteristics would predict a Muslim disadvantage relative to Hindus, previous studies documenting this robust and persistent pattern have called it a 'puzzle' of Muslim mortality". Basu et al. [16] much-explored demographic relationship between fertility and child survival, they explain that "one would expect infant and child mortality to be higher among Muslims than Hindus in India. This expectation is buttressed by the fact that Muslim women also display so many of the other features of high mortality situations-relatively low levels of education, greater conservatism in medical matters, higher levels of poverty. And yet, the National Family Health Survey finds that infant and child mortality rates in India are significantly lower among Muslims. These results are intriguing not only because they fly against the expectations based on religious differentials in bio-demographic factors like parity and birth-spacing. They are surprising because, at least in the binary tables, they do not correlate particularly strongly with many of the other proximate determinants of infant and child mortality; determinants like ante-natal care, delivery conditions, exposure to infections and treatment of illness". Maitra [17] and Kravdal [18] briefly note that Muslims have lower mortality than Hindus, but do not address why that is the case. This issue quite complicated to explain, a detailed study is needed on that topic. The various studies already mentioned the child survival benefit in context of child mortality and infant mortality analysis by religion but these puzzles not specified for studies where child mortality for women considers for total reproductive life span. This is again a point of research in child mortality.
The factor type of house has no significant effect for parity-3, 4 but for parity 5 it shows significant effect. But there are various studies on infant or child mortality namely Gyimah [19], Hill [20] and Godson [21] in which generalized linear models (GLM) are used to find the factors related to infant or child mortality. According to above mentioned studies, mother education, birth order, and mother's age, type of locality, socioeconomic status and health facilities are found as the determinants of child or infant mortality. But in above mentioned study's authors were used the data of child mortality for some limited reproductive period of women and total reproductive life span of women was not studied.