卢旺达儿童发育迟缓的时空变异与风险因素：基于地理加性模型的十年追踪研究

《BMC Public Health》：Spatio-temporal variation of child stunting and associated risk factors in Rwanda

编辑推荐：为解决卢旺达5岁以下儿童发育迟缓率持续高于国家目标的问题，研究人员利用2010-2020年三轮人口健康调查数据，采用贝叶斯地理加性模型分析时空变异规律。研究发现儿童性别、出生顺序、母亲身高等13项因素持续显著，且西部省份空间效应增强。该研究为制定区域精准干预策略提供多部门协同治理依据，对实现SDGs营养目标具有重要实践意义。

　　在撒哈拉以南非洲地区，儿童发育迟缓（Stunting）始终是严峻的公共卫生挑战。尽管全球5岁以下儿童发育迟缓率从2000年的33.1%降至2020年的22.0%，但卢旺达2020年仍高达33%，远高于非洲平均水平（30.7%）和国家设定的19%目标。更令人困惑的是，在实施"每户一牛"（Girinka）等本土解决方案后，不同地区效果差异显著——如2015年Nyabihu地区患病率达59%，而Kicukiro仅17%。这种地理异质性提示传统干预策略可能忽视区域特异性风险因素的动态变化。为此，卢旺达大学应用统计系Francois Niragire领衔的多学科团队，在《BMC Public Health》发表研究，首次揭示风险因素的时空演变规律。

研究团队采用贝叶斯地理加性逻辑回归模型（Bayesian geo-additive logistic regression），分析2010、2014/15、2019/20三轮卢旺达人口健康调查（RDHS）中11,422名儿童的横断面数据。通过马尔可夫随机场（MRF）先验处理空间自相关，结合P样条（P-splines）拟合非线性效应，构建包含结构化空间效应（f_struc）和非结构化效应（f_unstruc）的混合模型。数据预处理阶段保留0-59月龄儿童的身高年龄Z评分（HAZ），以WHO标准定义HAZ<-2SD为发育迟缓。

主要发现如下：

1. ?时间维度变异：母亲年龄效应呈现"V型→水平S型→倒U型"的演变，2010年年轻/高龄母亲风险更高，2020年转为中年母亲风险突出。母亲教育程度的影响在2020年减弱（POR=0.92, 95%CI:0.739-1.172），而农村居住风险持续上升（2010年POR=1.28→2020年1.38）。

2. ?空间异质性：西部Rubavu（14）、Nyabihu（15）等地区空间效应增强，后验概率图显示2020年这些区域总空间效应显著（p<0.05）。相反，东部Gatsibo（26）等地的空间效应十年间衰减68%。

3. ?稳定风险因素：男性儿童风险始终高于女性（POR≈1.5）；海拔>2000m地区风险增加64-79%；母亲身高<145cm使风险提高3倍以上（POR≈0.25）。

4. ?新现保护因素：2015年后，4次以上产检（ANC）的保护作用增强（POR从0.75降至0.66），可能与"村厨房"（Igikoni cy'Umudugudu）营养教育项目推广相关。

讨论指出，海拔相关的空间效应可能反映高原地带农业产出受限与缺氧代谢障碍的双重压力。而西部风险集聚提示需关注"牛群分配"政策在山区实施中的适应性障碍。该研究创新性地将时空交互项纳入健康地理分析框架，为卢旺达2024年营养战略提供三方面启示：（1）对Nyabihu等热点地区实施海拔适应的作物改良计划；（2）针对2-3岁儿童窗口期强化辅食补充；（3）建立母亲身高-出生体重联合筛查机制。这种"地理精准营养"（Geographically Precision Nutrition）策略，有望突破当前营养干预的"平原偏好"困局。