《Frontiers in Public Health》：Bayesian spatio-temporal modeling of childhood anemia in Tanzania: spatial distribution, determinants, and causal impact of malaria
编辑推荐：
引言： 儿童贫血在中低收入国家仍是重大公共卫生挑战，损害儿童生长、认知与生存。在坦桑尼亚，贫血负担持续居高不下，但其空间格局、时间趋势以及与传染病（尤其是疟疾）的关联尚不十分清楚。
方法： 研究人员利用四轮坦桑尼亚人口与健康调查（2007/08–2022）数据
引言： 儿童贫血在中低收入国家仍是重大公共卫生挑战，损害儿童生长、认知与生存。在坦桑尼亚，贫血负担持续居高不下，但其空间格局、时间趋势以及与传染病（尤其是疟疾）的关联尚不十分清楚。
方法： 研究人员利用四轮坦桑尼亚人口与健康调查（2007/08–2022）数据，采用贝叶斯时空地统计框架并结合调查加权因果推断，绘制了贫血患病率地图，并通过平均因果效应（ACE）、归因分数（AF）和需治疗人数（NNT）估计量化了疟疾的影响。
结果： 贫血患病率仍处于高位（总体为60.1%），在西部、南部和湖区地带存在持续性聚集；患病率从2007/08年至2015/16年下降，随后在2022年再次上升。疟疾被识别为贫血的最强预测因子，且其效应随儿童年龄增长而增强：在6–23月龄儿童中比值比（OR）为2.89（95%可信区间CrI：2.25–3.72），在42.59月龄儿童中升至4.07（3.53–4.71）。疟疾的平均因果效应（ACE）在2015/16年达到峰值0.266，同年归因分数（AF）为6.4%，需治疗人数（NNT）最低（3.93），表明该时期疟疾预防效率最高。最低温度高于约18°C与贫血风险增加相关，而在已考虑疟疾状态后，降水未显示出可信的关联；中海拔集群（约1000–1600米）贫血风险较低，与这些海拔疟疾传播减少一致。母亲教育程度较低、家庭贫困、家庭规模较大和母亲年龄较小也与较高的贫血风险相关。
结论： 这些发现表明，坦桑尼亚儿童贫血受根深蒂固的空间异质性以及疟疾之间强烈的、年龄依赖的、随时间变化的因果关系的驱动，凸显了将疟疾控制与更广泛的社会、经济和环境干预相结合的地理靶向多部门策略的必要性。
论文解读文章
一、研究背景与问题
儿童贫血是全球最持久、最紧迫的公共卫生挑战之一，影响5岁以下儿童的生长、认知发育和生存。贫血以血红蛋白或红细胞水平低下为特征，限制组织氧输送，损害生理和认知功能。尽管有全球健康倡议和营养干预，贫血仍影响全球近40%的5岁以下儿童和三分之一的育龄女性，负担最重集中在撒哈拉以南非洲和南亚。
儿童贫血的病因是多因素的，涵盖营养、感染和社会经济决定因素。缺铁仍是首要原因，但在许多热带地区，疟疾、HIV和寄生虫感染等传染病通过破坏红细胞和干扰铁代谢放大贫血负担。这些生物医学原因进一步被贫困、粮食不安全、卫生条件差和医疗可及性受限等背景因素加剧，形成感染、营养不良和社会经济脆弱性的持续循环。
在坦桑尼亚，儿童贫血持续构成重大公共卫生关切。坦桑尼亚人口与健康调查和疟疾指标调查数据显示，6–59月龄儿童贫血患病率持续偏高：2007/08年为68.8%，2011/12年为56.5%，2015/16年为57.4%，2022年为59.7%，2岁以下儿童患病率始终最高。坦桑尼亚已实施多项多部门战略，如国家多部门营养行动计划（NMNAP）和国家疟疾控制规划，但患病率下降有限，提示现有干预可能未能充分捕捉和应对该疾病复杂而局部化的驱动因素。贫血分布高度异质，疟疾传播强度高、膳食多样性差和医疗可及性受限的地区患病率持续较高。然而，现有证据多依赖全国平均值，掩盖了地方差异，也未能捕捉时间动态。传统统计方法常忽略空间依赖和时间变化关系，降低了识别高风险区域或厘清社会经济、环境和感染决定因素交互效应的能力。因此，亟需能够揭示贫血风险细尺度变化并量化疟疾与贫血因果关系的整合时空方法。
二、研究设计与技术方法
研究人员利用2007/08、2011/12、2015/16和2022年四轮坦桑尼亚人口与健康调查（TDHS）数据，共纳入28,151名6–59月龄儿童。贫血状态依据DHS标准生物标志物方案判定：现场使用电池驱动的HemoCue光度计测量毛细血管血红蛋白浓度，按集群海拔进行DHS/CDC标准校正，校正后血红蛋白低于11.0 g/dL（WHO对6–59月龄儿童的阈值）即判定为贫血。疟疾状态使用现场快速诊断检测（RDT）确定。气候数据取自WorldClim全球气候数据库。
研究采用贝叶斯时空地统计建模框架，通过积分嵌套拉普拉斯近似（INLA）实现，空间成分使用随机偏微分方程（SPDE）方法以高斯马尔可夫随机场（GMRF）近似。模型以logistic回归为框架，纳入固定效应协变量、连续协变量的平滑非线性项（使用一阶或二阶随机游走先验）、调查集群和分层独立随机效应，以及遵循一阶自回归（AR(1)）过程的时空随机效应，其空间创新项采用Matérn协方差函数的高斯随机场。连续气候协变量离散化为25个等频区间以保证计算可行性。由于模型规模庞大（701个网格节点×4个调查期，共2,804个潜在空间维度），采用经验贝叶斯积分策略。DHS调查权重通过按比例缩放各观测的似然贡献加以应用。
因果分析方面，由于在反事实预测中整合潜在空间场在计算上不可行，因果量从独立的分年份分层logistic回归模型中估计。研究估计了三个流行病学指标：归因分数（AF）、需治疗人数（NNT）和平均因果效应（ACE）。因果识别依赖条件可交换性、条件阳性、一致性和时间性四个假设，其中时间性无法从横断面数据验证，被明确承认为研究局限。ACE通过G计算（标准化）估计，点估计和预测按DHS样本权重加权，95%置信区间通过分层集群bootstrap获得。
三、研究结果
3.1 研究人群描述
四轮调查共纳入28,151名儿童，其中269名（0.96%）所在集群无GPS坐标，被排除于需要空间位置或海拔的模型，剩余27,882名儿童用于分年份因果模型。儿童平均年龄为32.13个月（SD=15.56）。贫血患病率在整个研究期间保持高位，60.1%的儿童被分类为贫血。男孩略超半数（50.3%），大多数儿童（80.5%）居住在农村。疟疾阳性率为10.8%，2007/08年最高（13.2%），2022年最低（6.9%）。
3.2 空间分布与时间趋势
区域层面莫兰指数（Moran's I）检验统计量为0.495（p<0.0001），表明区域层面存在正向空间聚集。贫血患病率在西部、南部和湖区地带存在持续性聚集簇。时间趋势上，患病率从2007/08年到2015/16年下降，随后在2022年再次上升。
3.3 疟疾与贫血的关联及因果效应
疟疾被识别为贫血的最强预测因子，其效应随儿童年龄增长而增强：6–23月龄儿童中比值比（OR）为2.89（95% CrI：2.25–3.72），42.59月龄儿童中升至4.07（3.53–4.71）。疟疾的平均因果效应（ACE）在2015/16年达到峰值0.266，同年归因分数（AF）为6.4%，需治疗人数（NNT）最低（3.93），表明该时期疟疾预防效率最高。
3.4 环境与社会经济因素
最低温度高于约18°C与贫血风险增加相关；在已考虑疟疾状态后，降水未显示出可信关联。中海拔集群（约1000–1600米）贫血风险较低，与这些海拔疟疾传播减少一致。母亲教育程度较低、家庭贫困、家庭规模较大和母亲年龄较小均与较高贫血风险相关。
四、讨论与结论
研究结果表明，坦桑尼亚儿童贫血由根深蒂固的空间异质性以及疟疾之间强烈的、年龄依赖的、随时间变化的因果关系所驱动。疟疾预防在2015/16年效率最高，提示该时期干预投入可能产生了较好回报。环境因素中最低温度和海拔通过影响疟疾传播间接作用于贫血风险，而社会经济因素（母亲教育、家庭财富、家庭规模、母亲年龄）独立于疟疾发挥影响。
研究结论指出，需要采取地理靶向的多部门策略，将疟疾控制与更广泛的社会、经济和环境干预相结合。这些发现为在资源有限的环境中优化干预资源配置提供了证据基础，强调了同时解决感染性和结构性决定因素的重要性。
