- Research article
- Open Access
Socioeconomic and demographic factors determining the underweight prevalence among children under-five in Punjab
BMC Public Health volume 20, Article number: 1817 (2020)
Underweight prevalence continues to be major public health challenge worldwide, particularly in developing countries like Pakistan. This study is focused on socio-economic and demographic aspects of underweight prevalence among children under-five in Punjab.
In this study, several socioeconomic and demographic factors are considered using MICS-4 data-set. Only those variables which are usually described in the nutritional studies of children were picked. Covariates include: the age of children, sex of the children, age of mother, total number of children born to women, family wealth index quintile, source of drinking water, type of sanitation, place of residence, parents’ education and occupation. All Categorical variables are effect coded. The child’s age and the mother’s age are assumed to be nonlinear, geographical region is spatial effect, while other variables are parametric in nature.
Since, the response is binary, covariate comprises linear terms, nonlinear effects of continues covariates and geographic effects, so we have use Geo-additive models (based on Fully Bayesian approach) with binomial family under logit link. Statistical analysis is performed on Statistical package R using Bayes X and R2 Bayes X Libraries.
Underweight status of children was found to be positively associated with number of under-five children in household, total number of children ever born to women and age of mother when the child was born. Whereas, it negatively associated with place of residence, parent’s education and family wealth index quintile. On the regional effect, the Southern Punjab has higher prevalence of underweight compared to Central and Northern Punjab.
Similarity of our results with several other studies demonstrate that the Geo-additive models are an ideal substitute of other statistical models to analyze the underweight prevalence among children. Moreover, our findings suggest the Punjab Government, to introduce target-oriented programs such as poverty reduction and enhancement of education and health facilities for poor population and disadvantaged regions, especially Southern Punjab.
Underweight prevalence is a major cause of child mortality worldwide. Underweight children have lesser immunity against infections and have higher chance to die from common diseases, and those who survive are exposed to recurring illnesses and slower growth. Such children are more likely to have lower IQ, which not only affect their educational performance but also reduce their working abilities [1, 2].
According to MICS-2014 , the prevalence of underweight among under-five children in Punjab is 33.7%, which is higher than the overall the underweight prevalence in the country (Fig. 1). It is a point of serious concern, because Punjab is comparatively most developed province of Pakistan regarding health and educational facilities. Additionally, it is also the most densely populated province of country, so a minor change in the percentage can be of greater significance. Therefore, a question arises what are the factors behind such a high rate of underweight prevalence among children in the province.
To the best of our knowledge, not even a single research study has been conducted in Punjab so far, to explore and analyze the impact of factors influencing such a high proportion of underweight.
Alasfoor et al. (2007), Mengistu et al. (2013) and Rayhan et al. (2006), use bivariate and multivariate analysis to identify the causes of under-five malnutrition [5,6,7]. The key factors for under-five malnourishment were the education level of parents, income of the household and total number of under-five children in household. Sapkota & Gurung (2009) and Wolde et al. (2014) used logistic regression techniques to calculate the underweight prevalence among children under-five [8,9,10]. They discovered that the socioeconomic status of family, maternal education and gender, ethnicity and age of children were closely linked with underweight prevalence. Some other authors [11,12,13,14,15,16] found that age and gender of children, child’s malaria status, exclusive breastfeeding, child’s vaccination status, maternal education, parent’s access to media, ARI in children, poverty and the type of toilet used in household were strongly correlated with underweight.
Majority of previous studies on child’s underweight have usually focused on different socio-economic, demographic or health related factors but mostly have ignored spatial and nonlinear effects of covariates. Considering these aspects, Kandala et al. (2008), Lasisi et al. (2014) and Lasisi et al. (2014) use Bayesian geo-additive models to observe the pattern of underweight among the children under five year of age [14, 17, 18]. Their studies showed that underweight among children depends on residence, age of the child, maternal education, wealth status of the household and geographical zone.
We aim to study the impact of several socio-economic and demographic factors on underweight status of under-five children in Punjab, considering the spatial and nonlinear effects of covariates. This study is based on data-set of the MICS-2014 conducted by Punjab Bureau of Statistics with the collaboration of UNICEF.
Two-stage cluster sampling technique was used for selection of sample. At first stage 2050 clusters were randomly selected (called primary sampling units) from all 36 districts of Punjab through proportional allocation. On second stage systematic sample of 20 households (Called secondary sampling units) was randomly drawn from each of selected cluster . The information about one or more explanatory variable(s) used in this study was missing for 6789 out of 31,083 children interviewed, so these were discarded. Finally, the total sample size for the research is 10,774 women of reproductive age (15–49) in Gilgit Baltistan.
Underweight status of children (UWa), i.e. ‟weight-for-age” is taken as a dependent variable. In this regard, the weight and the height of each child in sample is calculated separately. These measures are then expressed as Z-scores from the median of the reference population.
Where, WAPi is the weight for age percentile of ith child, while μWAP and σWAP denote the median and standard deviation of WAP of reference population respectively. A child ‟i” is declared as an underweight if Zi < − 2, else Normal weight. So, our response variable (UWa) is binary with 0 & 1 refers to Normal weight and underweight children respectively.
We consider several socio-economic and demographic covariates as predictor, which are generally described in the nutritional studies of children. Covariates include; the age of children, sex of the children, age of mother, total number of children ever born to women, family wealth index quintile, source of drinking water, type of sanitation, place of residence (locality), parents’ education and occupation. All categorical variables are effect coded; Children’s age and Mother’s age (when child born) are assumed nonlinear; ‟Region” is geographic effect, while other covariates are parametric in nature (Additional file 1).
Since, our response is binary; covariates consist of usual linear terms, nonlinear effects of continuous covariates and geographic effects, such sort of models are called Geo–Additive Regression Models, which are special case of Structure Additive Regression Models. We use Geo–Additive Regression under logit link, because we are interested to find the relative risk of underweight at different levels of given covariates. For inference, we have used a Fully Bayesian Approach because we want to describe them probabilistically. In this approach, all unknown regression coefficients and the smoothing functions are considered as random variables and must be assigned with appropriate prior distributions (Additional file 2). This methodology is based on Markov priors and uses Markov Chain Monte Carlo techniques to draw samples from posterior, and for model checking it normally used the Deviance Information Criterion. For detail, see [19,20,21,22].
Statistical analysis is performed on R Language using BayesX and R2BayesX Libraries (Additional file 2).
Before executing the final model, a bivariate analysis was performed, in order to decide which variables were to be included in the final model. The association between underweight status of children and different socio-economic and demographic variables (based on Pearson Chi-square tests of independence at 5% level of significance) is given in Table 1. According to Table 1, All variables, except the gender of children, found significantly associated with the underweight status of children. It is, therefore, the gender of children is discarded from final model (Additional file 2).
According to summary of results (Table 2), going from rural to urban areas decreases the log odds of underweight by 0.13.
Underweight status of Children under-five in Punjab is negatively associated with education level of father and mother of children. The log odds of underweight decreased by 0.09, 0.17, 0.28 and 0.41 for children whose fathers attain primary, middle, secondary and higher level of education respectively w.r.t. those having no education at all. The log odds of underweight decreased by 0.08, 0.19, 0.25 and 0.59 for children whose mothers got primary, middle, secondary and higher education respectively, compared to those having no education (Table 2).
The risk of underweight is significantly affected by father’s occupation. Here, unemployed category is taken as a benchmark. The log odds of underweight are increased by 0.08 and 0.16 for the children whose fathers were officials and laborer respectively; in contrast, they decreased by 0.06 for child whose fathers were farmers and remained almost unchanged for children whose fathers owned any business entity (small or large) (Table 2).
Underweight status of Children under-five in Punjab is negatively correlated with household wealth index quintile. Going from the lowest wealth index quartile to second, middle, fourth and highest wealth index quantiles, decreases the log odds of underweight by 0.18, 0.35, 0.45 and 0.83 respectively (Table 2).
Total number of children ever born to a woman is positively associated with the risk of underweight. Addition of every single child into the total number of children ever born to a woman increases the log odds of underweight by 0.04 (Table 2).
The household’s access to improved sanitation facility decreases the log odds of underweight by 0.09 compared to unimproved sanitation facilities (Table 2).
Underweight status of children is inversely related to the age of their mothers. The risk of underweight found highest in the children whose mother’s age was less than twenty years at time of their birth, and it decline down gradually as the age of mother increase. The spread of underweight is almost same for mother’s age group of 20–42 years (when children born), relatively higher beyond this age group, and highest if mother age is grater then 45 years (Fig. 2).
Underweight status of children is directly but non-linearly correlated with the age of children (in months) as shown in Fig. 3 and Table 3. It shows that underweight status of children is normal from birth up to about 12 months of age, then it gradually inclines and got severe after 18 months of age. This pattern continues until the age of 20 months and turn out to be almost static up to the age of about 42 months, afterward falls sharply. The prevalence of underweight is revealed to be more pronounced in age group of 18–32 months (Fig. 3).
Exposure of child’s mother to mass media reduces the log odds of underweight by 0.02 in comparison with no access to media at all (Table 2).
Refer to Table 2, the prevalence of underweight is more noticeable in the Southern Punjab. Going from Southern Punjab to Central and Northern Punjab decreases the log odds of underweight by 0.09 and 0.39.
The underweight status of children is slightly but not significantly affected by mother’s occupation, source of drinking water and total number of under-five children in household (Table 2).
This study is focused on various socioeconomic and demographic aspects of underweight prevalence among under-five children in Punjab province of Pakistan and may facilitate the concerned administration and policy makers in formulation of apt polices to address it.
However, it exclude others factors related to health of children and mothers such as birth interval, health status of mother during pregnancy, place of delivery, weight of a child at birth, exclusive breast feeding, vaccination status of children, child morbidity status (fever, measles, diarrhea and ARI), type of food he/she received, mother’s BMI, maternal and eternal care and use of extra food during pregnancy etc. In addition, there are three commonly used measures of child nutritional status: Stunting (height-for-age), Wasting (weight-for-height) and Underweight (weight-for-age). For the sake of conciseness only the last mentioned has been used, in this study.
Since, quality, availability and accessibility of hygienic food, health care and other facilities (like water and sanitation, exposure to media etc.) are considerably lower in rural Punjab, therefore, the children from rural areas of Punjab found more like to be underweight than their urban counterparts. This finding is consistent with findings obtained by other authors [12, 14, 17].
Our finding discloses a strong positive association between underweight and education level of father and mother of children. It is because, the educated parents usually have more knowledge about diet and health of children. In addition, education may change their traditional beliefs about diseases and taking care of their children. It also found consistent with several other similar studies [5, 6, 10, 12, 14,15,16,17,18]. Furthermore, the risk of underweight is 18% lower in children whose mother possesses higher education than fathers with same educational qualification (Fig. 4). It might because, the mothers usually spend more time with children compared with their fathers.
The risk of underweight found more severe in children belongs to families with lowest wealth index quintile. Because, it is impossible for parents with a very small income to maintain the increasing need of food as well as health expenditures of children, which results undernourishment. On the other hand, wealthier parents are more likely to afford hygienic food, healthier living conditions and better health-care that improves the child nutrition and overall health. Many other studies on child nutrition illustrate the similar results [5, 10, 12,13,14,15,16].
Our study has found positive relationship between the number of children ever born to a woman and underweight prevalence. Same findings are revealed in a study carried out by Islam et al. (2013) . Generally, the number of children ever born to a woman (fertility) have an inverse effect on child nutrition and health, which results the child to be underweight. It is because, families with more children, experience more economic strain for food consumption; hence, the children from these families are more likely to suffer from malnutrition and underweight. In addition, up to 40% of people in Punjab fall in 1St and 2nd wealth index quintile, i.e. low-income quintiles (Table 4). Whereas, this percentage goes up to 45% for families with one or more children in household under-five year of age. Therefore, it is impossible, for such people to maintain their livelihood and fulfill their basic needs. Hence, addition of every newborn baby in family may further reduce the distribution of income on health and food of children, consequently the risk of underweight increase.
The underweight status of children found negatively associated with improvement in sanitation facilities. It is because, the improved sanitation reduces the risk of diseases; consequently, the risk of underweight is decreased. The source of drinking water revealed similar but comparatively lesser impact on underweight status of children. It might because, the 94% inhabitants of Punjab use improved source of drinking water while comparatively lesser proportion of population (66%) use improved sanitation facility . In addition, the ‟improved source” indicator is based on the categorization of water supplies by type of facility and does not consider the direct measurements of water quality, i.e. the water from improved source is not necessarily be free from contaminants .
The underweight status of children is negatively associated with mother’s access to media (Table 2). It could be because, the exposure to media may lead to educate the mothers about health and food of children as well as more effective allocation of income to devote on their health and nutrition. These results are consistent with findings obtained by other authors [12, 17].
Refer to Fig. 4, the underweight status of children is inversely related to age of mother when child born. The surprising result for the underweight status of children whose mother age is greater than forty years is probably due to smaller frequency of observations with very high dispersion among them. Mothers of only 2.8% of children had the age 40 years or more when they born as well as the variation among the risk of underweight relatively higher for this age group.
The reasons for nonlinear association of underweight status of children with their age are that: usually children born with approximately normal anthropometric status; but afterwards, the health status of the most children became worse due to different socio-economic and biological factors and until it improved and become stable at a low level. The highest risk of underweight during the age group of 18–32 months is might because, during the aforesaid age group, kids are usually weaned from breast as well as their mothers do lose their ability to produce enough milk to fulfill the nutritional requirement of growing children. These results are consistent with findings obtained by other authors [10, 11, 14].
The causes for such a high rate of underweight prevalence in Southern Punjab Southern Punjab also called Seraiki belt (consisting of 11 districts) could be:
Firstly, Southern Punjab is more deprived region of province, 41% of people residing in this region have lowest wealth index quintile (Table 5).
Secondly, 72% of people living in Southern Punjab belongs to rural areas with poor health, living and educational facilities which is highest in entire province (Table 6).
Thirdly, the statistics of parent’s education of under-five children is worse in Southern Punjab. In this region, 64% of mothers and 41% of fathers of under five children are illiterate (Table 7). This percentage is highest in entire province (Table 7).
Fourthly, 45% of the inhabitant of Southern Punjab have no access to improved sanitation which is also highest in entire province (Table 7).
Fifthly, about one-half of resident of this region has no access to mass media at all. It is not only highest compared with rest of Punjab but horrible to see in this era (Table 8).
In addition, Southern Punjab is considered as most deprived and neglected part of Punjab. Since 1970, South Punjab has not been receiving its due share from provincial resources. Development schemes for said region were diverted to Central and Northern parts of province. According to recent data, poverty and unemployment are relatively higher in this region [1, 3, 24,25,26,27]. Furthermore, unlike rest of Punjab province, People in Southern part of province, which are under the sway of powerful landed gentry and tribal sardars (Tribal chiefs), suffer from lack of opportunity to improve their livelihood.
Therefore, to-a-days, for fair distribution of resources and administrative reasons, the demand for creation of new autonomous province in south Punjab, i.e. Seraiki Province (Seraiki is a native language of this region), is under national debate and focus.
The results of this study reveal that the family’s wealth index quintile, mother education level, father education level, geographical region, locality and total number of children ever born to a woman (fertility) are significantly associated with underweight status of children under-five in Punjab. The results from this study are found to be consistent with other studies [6, 8, 11, 13].
In order to achieve the sustainable development goal Punjab needs to scale up target-oriented programs such as poverty reduction income-generating interventions and improvement of education and health facilities for poor population and deprived parts of province especially Southern Punjab.
Recommendations for researchers:
There are three commonly used measures of child nutritional status: Stunting, Wasting and Underweight. Only the last-mentioned measure has been used in this study. Hence, Stunting and wasting may also be used for said purpose. In addition, this study is specific to Punjab province; same may be don for other provinces/parts of Pakistan.
Acute Respiratory Illness
Body Mass Index
Millennium Development Goals
Multiple Indicator Cluster Survey
National Nutrition Survey
United Nation Children’s Emergency Fund
World Health Organization
MICS. Multipe Indicator cluster survey. Lahore: Punjab Bureau Of Statistics; 2014.
Unicef. The state of the World's children. New York: Oxford University Press; 1998.
PDHS. (2012–13). Pakistan demographic and health survey.
The World Bank. (2015). http://data.worldbank.org/indicator/SH.STA.MALN.ZS. Retrieved November 24, 2015, from The World Bank: http://www.worldbank.org/.
Alasfoor, D., Traissac, P., Gartner, A., & Delpeuch, F. (2007 , September). Determinants of persistent underweight among children, aged 6–35 months, after huge economic development and improvements in health Services in Oman. J Health Popul Nutr, 25(3), :359–369.
Khan R, M. I., & Hayat, M. S. Factors causing malnutrition among under-five children in Bangladesh. Pak J Nutr. 2006;5(6):558–62.
Mengistu, K., Alemu, K., & Destaw, B. (2013). Prevalence of malnutrition and associated factors among children aged 6-59 months at Hidabu Abote District, north Shewa, Oromia regional state. Gondar: Institute of Public Health, College of Medicine and Health Sciences, University of Gondar, Gondar, Ethiopia doi:https://doi.org/10.4172/2161-0509.T1-001.
Sapkota V, Gurung C. Prevalence and predictors of underweight, stunting and Waisting in under-five children. JNHR. 2009;7(15):120–6.
Wolde, T., Belache, T., & Birhanu, T. (2014). Prevalence of Undernutrition and determinant factors among preschool children in Hawassa, Southern Ethiopia Food Science and Quality Management, 29.
Anderson A, Bignell W, Winful S, Soyiri I, Steiner-Asiedu M. Risk factors for malnutrition among children 5-years and younger in the Akuapim-North District in theEastern region of Ghana. Curr Res J. 2010;2(3):183–8.
Islam MM, Alam M, Tariquzaman M, Kabir MA, Pervin R. Predictors of the number of under-five malnourished children in Bangladesh: application of the generalized Poisson regression model. BMC Public Health. 2013;13(11) Retrieved from http://www.biomedcentral.com/1471-2458/13/11.
Kisiangani I, Mbakaya C, Mekokha A, Magu D. Prevelence of Malunaration among Preschool Children (6–59 months) in Westren Province , Kenya. J Public Health Epidimol. 2014;6(11):398–406. https://doi.org/10.5897/JPHE2014.0660.
Lasisi KE, Adebayo S, Abdulhamid B. Analysis of Nigerian children underweight nutritional status using Geoadditive cox models with Gaussian and binomial links. Afr Dev Resourc Res Institut J. 2014;10(2) Retrieved from http://www.journals.adrri.org/.
Nzala S, Siziya S, Babaniyi O, Songolo P, Muula A, Rudatsikira E. Demographic, cultural and environment factors associated with frequency and severity of malnutrition among Zambian children less than five years of age. J Public Health Epidemiol. 2011;3(8):362–70 Retrieved from http://www.academicjournals.org.
Zewdu S. Magnitude and factors associated with malnutrition of children under-five years of age in rural Kebeles of Haramaya, Ethiopia. Harar Bull Health Sci. 2012:221–32.
Kandala, N.-B., Fahrmeir, L., Klasen, S., & Priebe, J. (2008). Geo-additive models of childhood Undernutrition in three sub-Saharan African countries. University of Warwick. Retrieved from https://doi.org/10.1002/psp.524.
Lasisi KE, Nwaosu S, Abdulazeez KA. Geoadditive cox models with Gaussian and binomial links for the analysis of wasting status of Nigerian children. J Nat Sci Res. 2014:4.
Belitz, C. (2007). Model selection in generalised structured additive regression models. MÄunchen.
Fahrmeir K, Lang. Penalized structured additive regression for space-time ata: a Bayesian perspective. Stat Sin. 2004;14:731–61.
Kneib, & Fahrmeir. (2005). Structured Additive Regression for Categorical Space-time Data: A mixed model approach. Biometrics, to appear.
Kneib, T. (2005). Mixed model based inference in structured additive regression. Munchen.
Unicef. (2015). Water and sanitation. Retrieved from Unicef: http://data.unicef.org/water-sanitation/water.html.
MICS. Multiple Indicator cluster survey. Lahore: Punjab Bureau Of Statistics; 2007.
MICS. (2011). Multiple Indicator cluster Srvey Punjab. Lahore: Punjab Bureau of Statistics.
NNS. (2011). National Nutrition Survey. Aga Khan University, Pakistan, Planing Comossion ,Planing and Development Division ,Government of Pakistan.
Unicef. Retrieved from Unicef: http://mics.unicef.org/surveys.
Filmer D, Pritchett LH. Estimating wealth effects without expenditure data – or tears: an application to educational enrolments in states of India. Demography. 2001;38(1):115–32.
Rutstein, S. (2008). The DHS wealth index: approaches for rural and urban areas. Demographic and Health Research Division, Macro International Inc..
Rutstein, S., & Johonson, K. (2004). The DHS wealth index. DHS comparative reports no. 6. Calverton, Maryland: ORC macro, MEASURE DHS.
Cheema A, Khalid L, Patnam M. The geography of poverty: evidence from the Punjab. Lahore J Econ. 2008:163–88.
Gazdar, H. (1999). Review of Pakistan Poverty Data, Monograph Series # 9. Sustainable Development Policy Institute Islamabad.
Jamal H. The Punjab indices of multiple Depriviation 2003–04 and 2007–08. Lahore: Planing & Development Deportment, Government of Punjab; 2011.
Malik, S. J. ( 2005). Agricultural growth and rural poverty: review of the evidence. Asian Development Bank.
Brezger K, Lang. BayesX: analyzing Bayesian structured additive regression models. J Stat Softw. 2005;14(11).
Sahlin K. Estimating convergence of Markov chain Monte Carlo simulations. Sweden: Mathematical Statistics, Stockholm University; 2011.
Brooks, S. P., & Roberts, G. O. (1998). Convergence assessment techniques for Markov chain Monte Carlo. Kluwer academic publishers, Dordrecht NL-3300 AA Netherlands.
Cowles MK, Carlin BP. Markov chain Monte Carlo convergence diagnostics: a comparative review. J Am Stat Assoc. 1996;91:883–904.
Plummer, M., Best, N., Cowles, K., & Vines, K. (2010). Coda:convergence diagnosis and output analysis for MCMC.
Authors are thankful to Mr. Muhammad Amin and Syed Wasim Abbas for their support in managing data files and development of map. Authors are also thankful to the reviewers and editor for their efforts and expertise to refine this paper.
Availability of data and material
MICS 2014 dataset is publicly available online , public access to the data base is open and no ethical approval is required to access and use this data.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
: Figure B1 ACF for step = 1.Auto correlation function plot at step = 1 created through R language
: Figure B2 ACF for step = 10. Auto correlation function plot at step = 10 created through R language
About this article
Cite this article
Farooq, R., Khan, H., Khan, M.A. et al. Socioeconomic and demographic factors determining the underweight prevalence among children under-five in Punjab. BMC Public Health 20, 1817 (2020). https://doi.org/10.1186/s12889-020-09675-5
- Bayes X
- Fully Bayesian approach
- Geo-additive models
- Markov chain Monte Carlo
- Wealth index quintile