环境暴露与疟疾负担的国家间不平等（2000–2023）：制度能力、经济发展和结构性脆弱性的全球面板证据

《Global Public Health》：From environmental exposure to disease burden: Institutional capacity, structural vulnerability, and global malaria inequality, 2000–2023

　　

在环境条件相似的国家之间，疟疾负担（malaria burden）差异显著。研究人员分析了2000–2023年间疟疾发病率（malaria incidence）、制度能力（institutional capacity）、经济发展（economic development）、结构性脆弱性（structural vulnerability）与选定环境指标（environmental indicators）之间的国家层面关联。生态面板（ecological panel）包含186个国家的4,382个国家-年（country-year）完整观测。年固定效应模型（year-fixed-effects models）将log(incidence + 1)与标准化预测变量（standardised predictors）相联系，并采用国家聚类不确定性（country-clustered uncertainty）。结构性脆弱性与发病率呈正相关（β = 0.931，p = .007），而人均GDP（GDP per capita）与发病率呈负相关（β = ?2.014，p < .001）。调整后，制度能力和选定环境指标无统计学显著性。纳入国家和年份固定效应后，结构性脆弱性的关联减弱至β = 0.406（p = .232），而GDP仍呈负相关。排除与疟疾相关的ND-GAIN指标后，年份效应模型中的正关联仍然保留；排除整个卫生部门（health sector）以及改变部门覆盖范围会影响研究结果。GDP表现出相对一致的关联；结构性脆弱性的关联取决于模型设定（specification）和复合指标构建（composite construction）。国家层面聚合（national aggregation）、来源溯源（source provenance）不完整和潜在混杂（confounding）限制了结果解释。这些生态学关联不能建立因果通路（causal pathways）。需要利用次国家层面（subnational）的环境、干预和卫生系统数据来研究所提出的机制。

研究背景与问题：疟疾负担在全球分布不均，环境条件相似国家之间差异显著；2000年以来全球发病率下降，但许多低中收入国家进展停滞。温度、降雨和湿度影响媒介与寄生虫，但传播还依赖土地利用、住房、卫生系统、贫困、监测和干预覆盖。历史上欧洲和北美在适宜传播地区消除疟疾，伴随社会经济发展、基础设施改善和有效公共卫生措施，提示环境适宜性转化为疾病负担取决于社会与制度条件。制度能力（institutional capacity）可能通过监测、项目协调、采购、供应链、病媒控制可及性、及时诊断、持续治疗、社区卫健康工作者、疫情应对、问责融资和环境风险适应影响疟疾控制；结构性脆弱性（structural vulnerability）可能限制资源和制度能力利用。既有研究关注贫困、经济发展与疟疾负担，Sempungu et al. (2023)分析撒哈拉以南非洲治理、发展援助和疟疾结局。本研究因此分析2000–2023年全球发病率面板，联合考虑结构与环境指标，区分国家间模式与国家内关联，并检验脆弱性构建敏感性。

研究开展与意义：研究人员构建生态面板（ecological panel），源面板含4,848个国家-年、202个国家；经校正和列表删除后，全模型纳入4,382个国家-年、186个国家。研究以年固定效应模型（year-fixed-effects models）估计log(incidence + 1)与标准化制度能力（WGI composite）、log人均GDP（GDP per capita）、ND-GAIN结构性脆弱性、温度异常和年降水的关联，采用国家聚类不确定性；并以国家和年份双向固定效应、滞后预测、IHS变换、剔除疟疾相关ND-GAIN指标和卫生部门等作敏感性分析；探索性交互模型检验环境关联是否随制度或结构背景变化；VIF诊断共线性；未使用多重插补。该研究意义在于区分跨国结构梯度与国家内变化证据，并显式评估设定、复合构建和部门覆盖敏感性。

关键技术方法（≤250字）：研究使用2000–2023年186国4,382个国家-年生态面板；样本来源包括GBD（Global Burden of Disease）疟疾发病率估计、世界银行WDI/WGI（World Development Indicators/Worldwide Governance Indicators）、ND-GAIN（Notre Dame Global Adaptation Initiative）脆弱性复合指标、FAOSTAT温度异常和世界银行CCKP降水。以年固定效应OLS（ordinary least squares）估计log(incidence + 1)与标准化预测变量的关联，采用国家聚类稳健推断；以双向固定效应、滞后预测、IHS（inverse hyperbolic sine）变换、复合指标剔除和交互模型作敏感性分析，并计算VIF（variance inflation factor）。

研究结果：
**Analytical sample（分析样本）**：源面板为4,848个国家-年、202个国家；ND-GAIN校正和列表删除后，全模型为4,382个国家-年、186个国家；排除发生在所有区域，东亚太平洋损失19.6%，低收入国家-年损失20.7%。该分析确定全模型样本及其非随机缺失特征。
**Descriptive statistics（描述统计）**：疟疾发病率均值5,331.4例/10万（SD 11,754），中位0.0（IQR 0.0–1,436.6），53.4%国家-年为零发病，最大60,059例/10万；负担集中于撒哈拉以南非洲，并在南亚、东南亚和拉美部分地区集中。描述统计显示发病率强烈正偏和地理集中。
**Bivariate associations（双变量关联）**：Spearman相关显示，发病率与人均GDP（ρ = ?.734）、制度能力（ρ = ?.643）负相关，与结构性脆弱性（ρ = .664）正相关；温度异常ρ = ?.226，降水ρ = .157；结构预测变量互相关：GDP与制度能力ρ = .799，GDP与脆弱性ρ = ?.787，制度能力与脆弱性ρ = ?.617。双变量分析表明结构指标关联强于选定环境指标，且结构预测变量重叠限制调整模型分离。
**Institutional-capacity quartile gradients（制度能力四分位梯度）**：按WGI四分位，Q1均值11,871例/10万（中位3,033.5，零发病17.5%），Q2均值7,961.1（中位71.94，零30.6%），Q3均值1,387.9（中位0，零68.1%），Q4均值43.47（中位0，零96.5%）；Q1/Q4均值比约273.1，但受极值和小Q4均值影响，Q4中位为零使中位比未定义。四分位比较显示描述性梯度，非因果阈值。
**Multivariable panel regression（多变量面板回归）**：气候-only模型中温度异常β = ?1.011（95% CI ?1.374至?0.648，p < .001），降水β = 0.389（p = .177），调整R² = .069；结构-only模型中log GDP β = ?2.094（p < .001），脆弱性β = 1.090（p = .002），制度能力β = ?0.075（p = .850），调整R² = .548；全模型中脆弱性β = 0.931（95% CI 0.260–1.602，p = .007），log GDP β = ?2.014（95% CI ?2.878至?1.151，p < .001），制度能力β = ?0.183（p = .651），温度异常β = ?0.191（p = .166），降水β = 0.358（p = .140），调整R² = .549；指数化拟合比约2.54和0.13每1 SD，非算术均值比或因果风险比。全模型显示GDP负关联相对一致，脆弱性正关联，制度能力和选定环境指标调整后不显著。
**Sensitivity and robustness analyses（敏感性和稳健性分析）**：国家和年份双向固定效应后，脆弱性从β = 0.931（p = .007）减至β = 0.406（95% CI ?0.262至1.075，p = .232），GDP仍负关联β = ?1.337（p < .001）；原始归档系列中脆弱性从1.358减至0.486（p = .185）；排除ID_HEAL_02保留年效应正关联β = 0.995（p = .005），双向固定效应β = 0.372（p = .257）；排除整个卫生部门在同4,382观测中削弱关联β = 0.313（p = .310；双向β = 0.364，p = .220）；在3,655观测、153国且五个非卫生部门可用时，无卫生部门年效应β = 0.746（p = .015），双向β = 0.386（p = .239）；滞后和IHS模型保留正年效应，完整WGI结果相似；双向调整R² = .969。敏感性分析表明GDP关联较一致，脆弱性关联依赖设定、复合构建和部门覆盖。
**Exploratory interaction analyses（探索性交互分析）**：年效应模型中温度异常×制度能力β = 0.613（p < .001），×log GDP β = 0.537（p < .001），×脆弱性β = ?0.166（p = .177）；方向和幅度随双向固定效应及环境指标变化。交互分析为初步探索，不足以确立缓冲、放大或调节机制。

讨论总结：研究覆盖186国2000–2023年，经济发展与疟疾发病率呈相对一致负关联；脆弱性关联依赖模型设定、复合构建和部门覆盖，主模型正关联在双向固定效应后减弱且不显著。制度能力的双变量相关和四分位梯度显示治理质量与疟疾负担密切，但调整相关结构预测变量后估计不精确，不能识别独立关联；监测、项目协调、采购、供应连续性、社区卫健康工作者整合、疫情应对和问责仍为合理通路但未直接测量。结构性脆弱性模式不稳定；排除疟疾相关指标后年效应正关联保留，不能假定直接重叠解释；卫生部门排除依赖部门覆盖；复合指标描述结构条件但不识别疟疾服务、住房、融资或供应链机制。GDP关联需考虑反向因果：贫困可能经住房、预防和就医增加脆弱性，疟疾可降低生产力、阻碍人力资本积累并约束发展；未测政治、地理和历史因素可能影响二者，滞后估计不能解决。温度异常和年降水的弱调整关联不意味气候不重要；国家指标仅捕捉选定维度，国家平均掩盖局部传播、干预覆盖和卫生系统绩效。交互结果对缓冲或放大支持有限，社会生态韧性是解释框架而非已测机制。四分位描述负担梯度而非阈值，Q1/Q4均值比受极值影响。与Sempungu et al. (2023)比较需考虑覆盖、结局标准化、暴露定义、协变量及国家间/国家内估计目标；Deane et al. (2025)的乌干达定性证据提示生计相关暴露、性别责任和预防约束等局部通路，但不能作为国家系数因果证据。历史制度遗留、债务与财政约束、捐助优先、商品依赖和卫生技术可及等外部决定因素未测量，属背景假设。国家平均未解决国内公平，需次国家环境、干预、卫生系统及纵向数据。政策上，与项目证据结合，支持评估制度与结构投资作为病媒控制、化学预防、诊断、有效治疗和疫苗的补充，但本研究不识别其效果或相对有效性。局限包括观察设计不能因果、年初固定效应不能控制未测稳定国家特征、双向固定效应不能解决时变混杂和反向因果、国家聚合掩盖局部、环境指标不全面、完整病例选择可能影响结果、WGI和GDP及ND-GAIN测量限制、共线性限制精度、原始复合发布不可恢复及GBD发布标签未验证等；优势是24年186国多指标、显式比较跨国梯度与国内估计及替代脆弱性构建、冻结输入和可执行代码支持计算复现。结论翻译：人均GDP与全国疟疾发病率呈相对一致负关联。脆弱性从年效应模型β = 0.931减至国家和年份固定效应β = 0.406（p = .232）。这一不精确的国内估计不能确立脆弱性关系不存在。仅指标排除和整个卫生部门排除针对不同构念，部门覆盖进一步影响解释。核心贡献是区分广泛结构梯度与跨设定和测量仍稳定的关联。制度能力四分位不识别阈值，探索性交互不确立调节，弱调整气候关联不能推广到整个疟疾生态。制度加强和可靠提供仍是生物医学干预的相关补充，尽管其效果未在此识别。需要次国家环境、干预和卫生系统数据，结合处理混杂和反向因果的设计来研究这些通路。