Multivariate small area modelling of undernutrition prevalence among under-five children in Bangladesh

District-representative data are rarely collected in the surveys for identifying localised disparities in Bangladesh, and so district-level estimates of undernutrition indicators – stunting, wasting and underweight – have remained largely unexplored. This study aims to estimate district-level prevalence of these indicators by employing a multivariate Fay–Herriot (MFH) model which accounts for the underlying correlation among the undernutrition indicators. Direct estimates (DIR) of the target indicators and their variance–covariance matrices calculated from the 2019 Bangladesh Multiple Indicator Cluster Survey microdata have been used as input for developing univariate Fay–Herriot (UFH), bivariate Fay–Herriot (BFH) and MFH models. The comparison of the various model-based estimates and their relative standard errors with the corresponding direct estimates reveals that the MFH estimator provides unbiased estimates with more accuracy than the DIR, UFH and BFH estimators. The MFH model-based district level estimates of stunting, wasting and underweight range between 16 and 43%, 15 and 36%, and 6 and 13% respectively. District level bivariate maps of undernutrition indicators show that districts in north-eastern and south-eastern parts are highly exposed to either form of undernutrition, than the districts in south-western and central parts of the country. In terms of the number of undernourished children, millions of children affected by either form of undernutrition are living in densely populated districts like the capital district Dhaka, though undernutrition indicators (as a proportion) are comparatively lower. These findings can be used to target districts with a concurrence of multiple forms of undernutrition, and in the design of urgent intervention programs to reduce the inequality in child undernutrition at the localised district level.