1. Home
  2. Volume 25 Number 3, October 2025
  3. SAHAMONY

Simultaneous inference for empirical best predictors in generalized linear mixed models: A poverty study in West Java

NUR FITRIYANI SAHAMONY, KUSMAN SADIK, ANANG KURNIA

Abstract


Accurate poverty mapping at the district and municipal levels remains challenging due to small sample sizes in household surveys, which often result in unstable direct estimates. To address this issue, this study employs microdata from the 2023 National Socioeconomic Survey (SUSENAS) to estimate household-level poverty proportions across 27 districts and municipalities in West Java Province using a binomial Generalized Linear Mixed Model (GLMM) combined with the Empirical Best Predictor (EBP) and Simultaneous Confidence Intervals (SCI). The GLMM framework captures household characteristics and random area effects to account for spatial heterogeneity. Three SCI approaches—Bonferroni correction, Bootstrap-t, and the Simes procedure—were implemented to evaluate EBP uncertainty while controlling the family-wise error rate. Results reveal substantial disparities, with Tasikmalaya (21.7%), Bandung Barat (15.5%), and Cianjur (12.8%) consistently above the provincial average of (6.8%), while urban areas such as Cimahi, Bekasi, and Depok report poverty rates below 2%. All methods achieved full empirical coverage (ECP = 100%), although interval widths differed: Bonferroni produced the widest intervals (AIW = 44.99), Bootstrap-t yielded the narrowest and most efficient (AIW = 29.16), and Simes provided intermediate but highly consistent results (AIW = 33.24). These findings underscore the methodological importance of integrating GLMM, EBP, and SCI for small area estimation while offering practical insights for evidence-based policy development and poverty reduction strategies in Indonesia.


Keywords


Empirical best predictor; generalized linear mixed models; poverty estimation; small area estimation; simultaneous confidence intervals

References


. World Bank. “Poverty and Inequality Platform (PIP),” World Bank. Available: https://www.worldbank.org (accessed Jun. 30, 2025).

. Badan Pusat Statistik (BPS). 2022. Persentase Penduduk Miskin Maret 2022. (Jakarta: BPS-Statistics Indonesia).

. BPS Provinsi Jawa Barat. 2022. Profil Kemiskinan Provinsi Jawa Barat. (Bandung: BPS Provinsi Jawa Barat).

. BPS Provinsi Jawa Barat. “Kemiskinan Maret 2023 Provinsi Jawa Barat.” BPS Provinsi Jawa Barat. Available: https://jabar.bps.go.id (accessed Jun. 30, 2025).

. Rao J. N. K.; Molina I. 2015. Small area estimation. 2nd ed. (Hoboken, NJ: Wiley).

. Molina, I.; Rao, J. N. K. 2010. Small area estimation of poverty indicators. Can. J. Stat. 38(3) 369–385. DOI:https://www.jstor.org/stable/27896031.

. Reluga, K.; Lombardía, M.-J.; Sperlich, S. 2023. Simultaneous inference for empirical best predictors with a poverty study in small areas. J. Am. Stat. Assoc. 118(541) 583–595. DOI: https://doi.org/10.1080/01621459.2021.1942014

. Simes R. J. 1986. An improved Bonferroni procedure for multiple tests of significance. Biometrika. 73(3) 751–754.DOI: https://doi.org/10.1093/biomet/73.3.751.

. Zhang L.; Hall P.; Maiti T. 2022. Bootstrap t procedures for simultaneous confidence intervals for complex indicators in small area estimation. Stat. Probab. Lett. 188(1) 109446.

. Wang Y.; et al. 2022. Constructing statistical intervals for small area estimates under GLMM. BMC Med. Res. Methodol. 22 178.

. Frink N.; Schmid T. 2024. Small area estimation with generalized random forests: Estimating poverty rates in Mexico. arXiv preprint arXiv:2406.03861. Available: https://arxiv.org/abs/2406.03861 (accessed Jun. 24, 2025).

. Bugallo M.; et al. 2025. Small area estimation of labour force indicators under unit-level models. J. R. Stat. Soc. A (Stat. Soc.).188(1) 241–262.

. Butar F. B.; Lahiri P. 2003. On measures of uncertainty of empirical Bayes small-area estimators. J. Stat. Plan. Inference.112(1–2) 63–76.

. Hall P.; Maiti T. 2006. On parametric bootstrap methods for small area prediction. J R Stat Soc Series B Stat Methodol. 68(2) 221–238.

. Reluga; Katarzyna. 2020. Simultaneous and post-selection inference for mixed parameters. Doctoral Thesis. DOI: 10.13097/archive-ouverte/unige:138615.

. Romano J. P.; Wolf M. 2005. Exact and approximate stepdown methods for multiple hypothesis testing. J. Am. Stat. Assoc.100(469) 94–108.

. Hobza T.; Novák M.; Pustejovská P. 2018. The performance of the Bonferroni and step-down procedures under spatial dependence. Revstat Stat. J.16(1) 53–76.

. Chatterjee S.; Lahiri P.; Li H. 2008. Parametric bootstrap approximation to the distribution of EBLUP and related prediction intervals in linear mixed models. Ann. Stat. 36(3) 1221–1245.

. Benjamini Y.; Yekutieli D. 2001. The control of the false discovery rate in multiple testing under dependency. Ann. Stat. 29(4) 1165–1188.

. Pratiwi S.; Setianingrum L. 2022. Pola spasial dan tren kemiskinan di Indonesia tahun 2015–2020. J. Perenc. Pengemb. Kebijak. 2(3) 220–229.

. Suharyo W.; Akhmadi; Hastuti; Filaili R.; Budiati S.; Munawar W. 2005. Developing a poverty map for Indonesia: A tool for better targeting in poverty reduction and social protection programs – Book 4: Field verification. (Jakarta: The SMERU Research Institute).

. Bates D.; Mächler M.; Bolker B.; Walker S. 2015. Fitting linear mixed-effects models using lme4. J. Stat. Softw.67(1) 1–48.

. Muhammad F.; Muchtar M.; Sihombing P. R. 2024. The nexus of gender inequality and poverty rate in Indonesia. Perspektif: Perspektif: J. Ekon. Manaj. Univ. Bina Sarana Inform.22(1).

. Hasanujzaman M.; Omar M. A. 2022. Household and non-household factors influencing multidimensional energy poverty in Bangladesh: Demographics, urbanization and regional differentiation via a multilevel modeling approach. Energy Res. Soc. Sci.92 102083. DOI: https://doi.org/10.1016/j.erss.2022.102803.

. Pangeran M.; Arham M. A.; Dai S. I. 2022. The effect of education and health performance on poverty in Indonesia. Eur. J. Res. Dev. Sustain. 3(6) 119–124. https://www.scholarzest.com.

. Soseco T. 2021. Household size, education, and household wealth in Indonesia: Evidence from quantile regression. J. Ekon. Indones. 10(3) 281–297.DOI: https://doi.org/10.52813/jei.v10i3.72.

. Wahyuni W.; Gazali M.; Hidayaturrohman U. 2022. Pengelompokan dan pemetaan karakteristik kemiskinan di Provinsi Nusa Tenggara Barat menggunakan self organizing map (SOM) dan biplot. Syntax Literate. 7(11) 15587–15605.

. Bryant G.; Spies-Butcher B.; Stebbing A. 2024. Comparing asset-based welfare capitalism: Wealth inequality, housing finance and household risk. Hous. Stud. 39(2).DOI: https://doi.org/10.1080/02673037.2022.2056150.

. Fauziyah E.; Awang S. A.; Suryanto P.; Achmade B. 2025. Inequality and poverty of privately owned forests farmers in rural areas of Indonesia. For. Sci. Technol. 21(1) 1–14. DOI: https//doi.org/10.1080/21580103.2024.2409219.

. Fikri A. T.; Afifudin; Sekarsari R. W. 2023. Evaluasi implementasi program bantuan sosial dalam menekan angka pengemis di Kota Malang (studi kasus di Dinas Sosial Kota Malang). Respon Publik. 17(9) 11–21. https://jim.unisma.ac.id/index.php/rpp/article/view/21914.

. Tzavidis N.; Zhang L.-C.; Luna-Hernández Á.; Schmid T.; Rojas-Perilla N. 2018. From start to finish: A framework for the production of small area official statistics. J. R. Stat. Soc. A (Stat. Soc.) 181(4) 927–979.DOI: https://doi.org/10.1111/rssa.12364.

. Hall P.; Maiti T.; Parnell K. 2011. Confidence regions for empirical best linear unbiased predictors in small area estimation. The Annals of Statistics. 39(6) 2862–2882.

. Efron B.; Tibshirani R. J. 1994. Ann. Stat. (New York: Chapman & Hall/CRC).


Full Text: PDF

DOI: 10.24815/jn.v25i3.47540

Refbacks

  • There are currently no refbacks.