Modeling Maternal Risk Factors Affecting Low Birth Weight Among Infants in Kenya
Nancy Wairimu Gathimba1, *, Anthony Wanjoya1, George Kipruto Kiplagat1, Levi Mbugua2, Koima Kibiwott3
1Department of Statistics and Acturial Sciences, Jomo Kenyatta University of Agriculture and Technology, Nairobi, Kenya
2Department of Statistics and Acturial Sciences, Technical University, Nairobi, Kenya
3Department of Statistics and Acturial Sciences, Kabarak University, Nakuru, Kenya
To cite this article:
Nancy Wairimu Gathimba, Anthony Wanjoya, George Kipruto Kiplagat, Levi Mbugua, Koima Kibiwott. Modeling Maternal Risk Factors Affecting Low Birth Weight Among Infants in Kenya. American Journal of Theoretical and Applied Statistics. Vol. 6, No. 1, 2017, pp. 22-31. doi: 10.11648/j.ajtas.20170601.13
Received: December 8, 2016; Accepted: December 20, 2016; Published: January 18, 2017
Abstract: The strong association of birth weight with infant mortality is the main focus of birth weight research, with the assumption that birth weight is a major determinant of infant survival. Studies on factors of low birth weight in Kenya have neglected the flexible approach of using smooth functions for some covariates in models. Such flexible approach reveals detailed relationship of covariates with the response. The study sought to investigate risk factors of low birth weight in Kenya by assuming a flexible approach for continuous covariates and geographical random effect. The study used semi parametric models to flexibly model the effects of selected covariates and spatial effects. Our spatial analysis is based on a flexible geo-additive model using the provinces as the geographic unit of analysis, which allows to separate smooth structured spatial effects from random effect. A Gaussian model for birth weight in grams and a binary logistic model for the binary outcome were fitted. Continuous covariates was modeled by the penalized (p) splines and spatial effects was smoothed by the two dimensional p-spline. The specific objectives of the study was to investigate factors of low birth weight in Kenya by taking into account the hierarchical nature of child birth weight factors using a Bayesian hierarchical model. The study used secondary data from Kenya 2014 demographic and health survey (KDHS) data. The study found that child birth order, mother weight and height are significant predictors of birth weight. Secondary education for mother, birth order categories 2-3 and 4-5, wealth index of richer family and mother height were significant predictors of child size at birth. The area associated with low birth weight was North Eastern and areas with increased risk to less than average size at birth were Central and Eastern. The study found support for the flexible modeling of some covariates that clearly have nonlinear influences. Nevertheless there was no strong support for inclusion of geographical spatial analysis. The spatial patterns and the maps generated could be used for targeting development efforts at a glance. These findings have important implications for targeting policy as well as the search for left-out variables that might account for these residual spatial patterns.
Keywords: Low Birth Weight, Geo-additive Models, Spatial Effect, Penalised Splines, Semi Parametric Bayesian Analysis
Background of the study
The birth weight of an infant is highly significant in two important aspects. First, it is strongly conditioned by the health and nutritional status of the mother and second, it is one of the most important determinants of the chances of the newborn to survive and experience healthy growth and development. According to World Health Organization around 3.3 million babies die within their first month of life every year and the proportion of neonatal deaths has increased in all regions of the world and is currently estimated at 41% , . 90% of these deaths occur in low and middle income countries and a third of all neonatal deaths occur in sub Saharan Africa. Low birth weight continues to be a significant public health problem globally and is associated with a range of both short- and long-term consequences. Overall, it is estimated that 15% to 20% of all births worldwide are Low Birth Weight, representing more than 20 million births a year. In 2012, the World Health Assembly Resolution 65.6 endorsed a Comprehensive implementation plan on maternal, infant and young child nutrition, which specified six global nutrition targets for 2025,Which include: 40% reduction in the number of children under-five who are stunted, 50% reduction of anaemia in women of reproductive age, 30% reduction in low birth weight, No increase in childhood overweight, Increasing the rate of exclusive breastfeeding in the first six months up to at least 50% and maintaining childhood wastage to less than 5%; . This resolution covers the third target: a 30% reduction of low birth weight. The purpose of this resolution is to increase attention to, investment in, and action for a set of cost-effective interventions and policies that can help Member States and their partners in reducing rates of low birth weight. The goal is to achieve a 30% reduction of the number of infants born with a weight lower than 2500 g by the year 2025. This would translate into a 3.9% relative reduction per year between 2012 and 2025 and a reduction from approximately 20 million to about 14million infants with low weight at birth. The 48 hours immediately following birth is the most crucial period for newborn survival. This is when the mother and child should receive follow-up care to prevent and treat illness. Prior to birth, the mother can increase her child's chance of survival and good health by attending antenatal care consultations, being immunized against tetanus, and avoiding smoking and use of alcohol. The prevalence of Low Birth Weight in Kenya is estimated to be 11 percent by World Health Organization and United Nations Children's Fund. The Kenya Demographic Health Survey of 2009 estimates Low Birth Weight to be 6 percent. Low Birth Weight is a major cause of morbidity and mortality in Kenya.  The current trends of infant and under five mortality rates in Kenya are declining. Child mortality has declined by 36 percent from 115 per 1000 in 2003 to 74 deaths per 1000 in 2008; , While infant mortality has declined by 32 percent from 77 deaths per 1000 in the 2003 to 52 deaths per 1000 in the 2008. Neonatal mortality changed from 33 deaths per 1000 reported in the 2003 to 31 deaths per 1000 reported in the 2008. The recorded decline of the infant and child mortality is an indicator of progress in achieving the fourth millennium development goal. Addressing challenges associated with newborn deaths in Kenya has the greatest potential of contributing to this progress. Low birth weight defined as less than 2500g is one of the biggest risk factors associated with neonatal deaths. Low birth weight babies may face serious health problems as newborns and are at increased risk of long-term disabilities. Researchers trail population trends in birth weight, assuming these to have implications for infant mortality trends; . In the United States of America, interventions to increase birth weight are recommended as a strategy to improve infant mortality, recommendations on Progress towards achieving the 1990 objectives for pregnancy and infant health. . There is an absence of data to explain the way in which low birth weight, small for gestational age, and preterm risks interact with neonatal mortality in high mortality burden settings. In the United States, where the neonatal mortality rate is relatively low at five per 1,000 live births, , mortality outcomes are reported to vary across groups of weight for gestational age, , with newborns born small for gestational age at 34–36 week gestation estimated to have a neonatal mortality risk as much as 44 times higher than the risk experienced by newborns with appropriate weight for gestational age, , . There are multiple causes of low birth weight, including early induction of labour or caesarean birth, multiple pregnancies, infections and chronic conditions such as diabetes and high blood pressure but the main causes are preterm birth and restricted growth in the womb resulting in babies who are referred to as "small for their dates"; . Premature birth is the primary cause of Low Birth Weight and it describes babies born before 37 completed weeks of gestation. The earlier a baby is born, the more likely it will weigh less. This is because the baby has less time in the mother’s uterus to grow and gain weight .
Pregnancy at extreme of reproductive age, both young teenagers and older women (after 40 years), are at high risk for giving birth to low birth weight babies, . His results of multivariate analysis showed a U-shaped relationship between maternal age and Low Birth Weight among whites, with the youngest (younger than 15) and the oldest (aged 40 and older) mothers being at high risk than 25-29 years old; older teenagers were not at any significantly higher risk. The consequences of low birth weight include fetal and neonatal mortality and morbidity, poor cognitive development and an increased risk of chronic diseases later in life; . The contribution of low birth weight to neonatal morbidity and mortality in Kenya cannot be ignored. Neonatal intensive care is not readily available because of its initial and running costs and where available the bed capacity is extremely low. Large rural district hospitals including those with pediatricians are poorly equipped to provide essential services to sick newborns and hence the need to implement simple, cost-effective and sustainable interventions to care for the special needs of newborns, . The need to focus on how to increase access to cost effective interventions that include control of the quality of infants born so as to decrease the burden and adopting simple strategies for the management of the high risk new born have been proposed, . Providing local solutions to public health problems have been found to be more acceptable and more likely to be implemented.
First univariate logistic regression was performed to select potential covariates for the multiple regression models. Covariates that were significant were considered as candidate variables for multiple regression models. Cross tabulations between categorical covariates and categorized birth weight (\ > 2.5kg or \ > = 2.5kg) and child size at birth was done to have percentage distributions of low birth weight and size at birth per covariate categories. The histogram of birth weight in kilograms was plotted to see the plausibility of Gaussian model.
The following multiple variable hierarchical model was fitted in R
Where g is the link function linking the mean of the response to the predictor
and is the area level random effect representing unmeasured contextual factors. In case of child birth weight as a response, the link function will be the identity link resulting in the Gaussian regression model. For the child size at birth as the response, the link function was the log of odds of less than average size at birth resulting in the binary logistic model. To take a more flexible approach, the continuous covariates and the area level random effects was modeled by the nonlinear smooth functions. This revealed their subtle influences which could not be shown if modeled parametrically. To reflect this flexible approach, (1) was changed to a geo-additive model
wherefor j=1, 2, 3,....., p are smooth functions expressing non linear relationship betweeen the response variable and the continuous covariate and is the area of the child random effect. The vector of coefficients determine the parametric relationship between the response and the categorical covariates. The smooth functions
was specified as Bayesian splines. This assumes approximating by the polynomial spline of degree l defined at equally spaced knots, 
Which are within the domain of the covariate xj. The Bayesian spline can be written as a linear combination of d =s+l basis functions, Bm, that is,
Now Bayesian estimation of the penalized spline (3) is equivalent in estimating model parameters
where first or second order random walk priors for the regression coefficients are assigned. A first order random walk prior for equidistant knots is given by:
where m = 2, 3,... d and a second order random walk prior for equidistant knots is given by: where m = 3, 4,..., d and ~ are random errors. The spatial effect was modeled by the tensor product of two dimensional p-spline defined as
Where () refers to the coordinates of the location of the data point, latitude and longitude, or location centroids based on the map. Note that
represents the effect of correlated unmeasured or unobserved location effects. The prior for
is based on spatial smoothness priors common in spatial statistics
 The most commonly used prior specification based on the four nearest neighbours was defined as: ~
For i, j = 2,..., k-1 with appropriate changes for corners and edges. Since model estimation will be by empirical Bayesian method, all variance parameters will be treated as unknown constants that will be estimated by restricted maximum likelihood (REML) method and hence their priors will not be given. The fixed effects was assigned diffuse priors. An advantage of the empirical Bayesian inference over full Bayesian inference is that questions about the convergence of Markov Chain Monte Carlo (MCMC) samples or sensitivity on hyper parameters do not arise .
Further more, a comparison of full Bayesian and empirical Bayesian approach in a simulation study, has shown empirical Bayesian approach yielding better point estimates, especially for Bernoulli distributed responses .
3. Results and Discussion
The data used are from the 2014 Kenya Demographic and Health Surveys. The DHS collect information on a nationally representative sample of women in child-bearing age (15-49). The questionnaire collects socioeconomic indicators for the respondent and her partner as well as the household she resides in, and then gathers a large amount of information on fertility patterns, health and care practices, health knowledge, and assesses the anthropometric status of all children of these women who were born within the past five years.
3.2. Descriptive Summaries
The percentage of low birth weight infants is higher among young mothers (aged 20 years or less) and in older mothers (aged 35 to 49 years) than in mothers aged 20 to 34 years (Table 1). The inverse relationship is also observed among wealth quintile. As education and household wealth increases, the percentage of low birth infants decreases. For instance, the percentage of low birth weight decreases from 9.4 percent among mothers with no education to 6.5 percent among mothers with more than secondary education. Likewise percentage of births in which infants weigh less than 2.5kg decreases from 8.5 percent among mothers in the lowest wealth quintile to 7 percent among mothers in the highest quintile. Among the regions (Table. 2) Nyanza region has the smallest proportion of low birth weight infants and the coastal region has the highest (3.5 and 12.7 percent respectively). Similar patterns in education and wealth quintile are seen for births categorized as very small and smaller than average as was seen for births less than 2.5kg.
To check the suitability of a Gaussian model, a histogram of birth weight was plotted. Fig 1 gives the histogram of birth weight in grams. The histogram show that birth weight is symmetrically distributed. This leads to a simpler model and analysis since the Gaussian assumption is more tenable.
|variable||Birth weight less than 2.5kg||Childsize very small||smallerthan average||Average or larger|
|Mother’s age at birth|
|Does not smoke||7.6||3.3||11.6||83.6|
|Birth weight||size of the child|
|less than 2.5kg||very small||smaller than average||Average or larger|
3.3. Empirical Bayesian
3.3.1. Fixed Effect
Fixed variables associated with child birth weight under the Gaussian model are birth order, mother weight, and mother height (Table 3). The birth order effects are positive which means children with higher birth order are associated with higher child birth weight than children of lower birth order. Positive effect of mother weight means that as mother weight increase child birth weight also increase. The positive effect for mother height means that the taller the mother the higher the child birth weight as well. Factors associated with child size at birth under binary logistic model are birth order 2–3 and 4–5, wealth index of richer family, mother education of secondary category, and mother height (Table 3). The effects of birth order 2–3 and 4–5 are negative which means children of birth order 2–3 and 4–5 are associated with lower risk of being small at birth compared to those with first birth order. Mother secondary education has a negative effect on child size at birth which means children of secondary.
|Variable||Logistic Coefficient(95%CI)||Gaussian Coefficient(95%CI)|
|Intercept||-0.481(-1.002 0.039)||2.851(2.645 3.058)|
|2 to 3||-0.351*(-0.657 -0.42)||-0.105*(0.04 0.196)|
|4 to 5||-0.530*(-0.951 -0.128)||0.138*(0.018 -0.259)|
|6+||-0.407(-0.901 0.090)||0.166*(0.017 -0.313)|
|45 to 70kg||-0.262(-0.547 0.023)||0.162*(0.064 0.263)|
|>70kg||-0.261(-0.774 0.251)||0.214*(0.065 0.363)|
|Primary||-0.201(-0.450 0.048)||-0.004(-0.091 0.082)|
|Secondary||-0.397(-0.789 0.017)||-0.20(-0.132 0.090)|
|Higher||-0.702(-2.194 0.793)||-0.160(-0.466 0.146)|
|Poor||-0.152(-0.413 0.112)||0.059(0.028 0.146)|
|Rich||-0.226(0.493 0.0044)||0.027(0.058 0.112)|
|Richer||-0.340*(-0.632 0.048)||0.047(-0.040 0.134)|
|Richest||-0.301(-0.643 0.040)||0.056(-0.039 0.152)|
|>150||-0.500*(-0.741 0.260)||0.126*(0.044 0.208)|
|Non -linear effects|
-means reference category
*means significant at 5% significance level
Negative effect of mother height on child birth size means that children of mothers whose height is equal to or greater than 150 cm are less likely to be smaller than average size at birth than children of mothers whose height is less than 150 cm.
3.3.2. Nonlinear Effects
Starting with the nonlinear effects to child birth weight (Fig 2 left), children of young mothers (aged 15 to 23 years) and older mothers (aged 35 to 49 years) are more likely to have low birth weight than children of mothers aged 23 to 35 years. These results are in line with the . Furthermore as number of antenatal visits for pregnancy increase, child birth weight also increases. With regard to nonlinear effects to child size at birth (Fig 2 right), children of mothers aged 15 to 25 years and children of mothers aged 35 to 49 years are prone to have small size at birth than children of mothers aged 25 to 35 years. Children whose mothers have less prenatal visits are prone to be small at birth. The study has found that, women at the reproductive ages of 25 years or less and 35 years and over are more prone to deliver low birth weight or small sized babies. (Fig 2 top left and Fig 2 top right). Mothers less than 25 years are actually prone to have physical and emotional maturity issues which may contribute to their elevated incidence of small size births or low birth weight infants. Their ignorance of how to take care of themselves during pregnancy works against child birth weight or size at birth. Accordingly, among mothers who are 35 years or older, there is a greater tendency to develop prenatal complications and a higher probability of inadequate nutrition, thus increasing their likelihood of delivering less than average size or low birth weight babies. The study has further shown that mothers whose prenatal visits are less than four are prone to have low birth weight or less than average size babies (Fig. 2 bottom left and Fig. 2 bottom right). Increased prenatal visits ensure mothers receive adequate diet literacy which helps improve child birth weight.
3.3.3. Spatial Effects
Spatial statistics extends traditional statistics to support the analysis of geographic data. It provides techniques to describe the distribution of data in the geographic space (descriptive spatial statistics),analyze the spatial patterns of the data (spatial pattern analysis), identify and measure spatial relationships (spatial regression), and create a surface from sampled data.
(a) Descriptive Spatial Analysis
Nyanza province has a higher percent of weight normal birth (greater than 2.5Kg), followed by Western and Rift valley (Fig. 3). Nairobi, Central and coast regions are associated by decreasing birth weight. With regard to child size at birth, Nyanza has the lowest percentage of very small size, followed by Nairobi, Eastern and North Eastern has the largest percent. (Fig. 4) and (Fig. 5).
(b) Spatial Regression Analysis
Figure 7 presents total residual spatial effects to child size at birth. North Eastern, Eastern, Coast, Western and Nyanza are associated with increased risk of child being small at birth(positive effects) while Rift valley and central are associated with decreased(negative effect) risk of child being small at birth. Rift valley and Western are associated with increased birth weight (Fig. 6) while coast and central provinces have a decreasing birth weight. The observed residual spatial heterogeneity may be due to unobserved factors not captured by the covariates in the models, and it is a matter of conjecture to identify them. Regional nutritional disparities may also explain the spatial heterogeneity of low birth weight.
The cause of regional nutritional differences can be natural disasters like floods, and varying climatic conditions. This study employed the use of geo-additive logistic model to study the relationship between birth weight and its risk factors. The geo-additive model allowed the mapping of residual spatial effects to child size and birth weight while accounting for non-linear covariate effects under the assumption of additiviness. This study found that child birth order, mother weight, mother height, mother education and family wealth are significant predictors of birth weight. The study is in consistent with the findings of the previous study done by, . The positive effect of mother's education on birth weight in this study is in consistent with other studies, ,  which found mother education as an important factor influencing birth weight. The finding of a positive effect of mother weight, and height on birth weight is also in line with, , , where they concluded that, such a relationship is because mother weight and height reflect food taken which has a direct influence on child birth weight. The mechanisms associated with small size at birth among the less educated, may include poor diet as a result of low dietary literacy . Furthermore limited education may also result in limited access to prenatal care, especially in settings where clients are expected to pay for service. Positive effect of family wealth on child birth weight may be due to the fact that wealth is associated with income level which determines kind of diet.
The study has also documented that women at the reproductive ages of 25 years or less and 35 years and over are more prone to deliver low birth weight or small sized babies (Fig 2 top left and Fig 2 top right). Mothers less than 25 years are actually prone to have physical andemotional maturity issues which may contribute to their elevated incidence of small size births or low birth weight infants. Their ignorance of how to take care of themselves during pregnancy works against child birth weight or size at birth. Accordingly, among mothers who are 35 years or older, there is a greater tendency to develop prenatal complications and a higher probability of inadequate nutrition, thus increasing their likelihood of delivering less than average size or low birth weight babies. The study has further shown that mothers whose prenatal visits are less than four are prone to have lowbirth weight or less than average size babies (Fig 2 bottom left and Fig 2 bottom right). Increased prenatal visits ensure mothers receive adequate diet literacy which helps improve child birth weight.
The study was not without weaknesses. Due to the cross-sectional nature of the data collection exercise, no temporal linkages can be made between birth weight or size at birth and any of the explanatory variables. Moreover, because the analysis was based on an existing data set, the study was limited to the use of variables found in the data set.
4. Conclusions and Recommendations
In this study we pooled the data from the 2014 demographic and Health surveys of Kenya to model the socio-demographic and spatial determinants of Low birth weight. The study found strong support for our approach of flexibly modeling some covariates that clearly have non-linear influences and for including a spatial analysis. The spatial effect of the low birth weight outcomes appeared weak across the various regions. The spatial analysis shows distinct spatial patterns that point to the influence of omitted variables with a strong spatial structure or possibly epidemiological processes that account for this spatial structure.
The maps generated could be used for targeting development efforts at a glance, or for exploring relationships between welfare indicators and others variables. For example, birth weight less than 2.5Kg map could be overlaid with maps of other types of data, say on size of birth very small. The visual nature of the maps may highlight unexpected relationships that would be overlooked in a standard regression analysis
We recommend that future related studies involving various adverse birth weight outcomes be traced at district level. This will enable minimal mismatch between risk factors and the low birth weight.
I am most grateful to Almighty God who through His infinite mercy and love guided me throughout the duration of writing this Project report. This work has been accomplished with the assistance of many people who deserve special mention for their support. My deep appreciation goes to my supervisor Dr. Anthony Wanjoya and Dr. Levi Mbugua who constantly told me it could be done and prompted me to continue when I had lost hope. I am particularly inspired by their dedication in seeing this work through. To them I say many thanks and God bless.