《Frontiers in Public Health》：National income inequality and suicide mortality across 170 countries, 2000–2021: within-country evidence and sensitivity to exposure uncertainty, with implications for health communication
编辑推荐：
引言：收入不平等是远端人群暴露，生态面板数据不能识别个体层面自杀机制。研究人员考察了国家可支配收入分配不平等的国内变化是否与自杀死亡率变化相关，并评估估计关联如何依赖于暴露不确定性、结局测量、样本构成、人口加权以及时间识别假设。
方法：研究人员将世界卫生组织(
引言：收入不平等是远端人群暴露，生态面板数据不能识别个体层面自杀机制。研究人员考察了国家可支配收入分配不平等的国内变化是否与自杀死亡率变化相关，并评估估计关联如何依赖于暴露不确定性、结局测量、样本构成、人口加权以及时间识别假设。
方法：研究人员将世界卫生组织(WHO)年龄标准化自杀率与标准化世界收入不平等数据库(SWIID)可支配收入基尼(Gini)估计值关联，覆盖170个国家2000–2021年共3,257个国家?年。主要设定为完全调整、未加权的双向固定效应(TWFE)模型，采用国家聚类标准误。研究人员使用SWIID 9.92版全部100条已完成序列，按Rubin规则与Barnard–Rubin自由度进行合并。敏感性分析考察人口加权、时间设定、替代结局、报告与测量质量子集及国家影响。研究人员还开展观测面板模拟，并在每基尼点±0.010对数率单位内做等效检验。
结果：后验均值代入估计为b = ?0.0157（95% CI，?0.0286至?0.0027；p = 0.018），而SWIID多重插补估计为b = ?0.00831（95% CI，?0.01688至0.00026；p = 0.057；缺失信息比例=0.173）。该敏感性区间跨过后验最小感兴趣效应|b| = 0.010。人口加权与对数人口加权估计分别为?0.00861（p = 0.193）与?0.00864（p = 0.047）。国家趋势、区域?年、一阶差分与动态估计更小且不精确。全部100条已完成数据系数与170个国家逐一剔除后验均值系数均为负，而WHO不确定区间较窄的子集更负。
结论：国家收入不平等与自杀死亡率存在小规模负向国内生态关联，其幅度不精确且依赖暴露处理、结局构建、样本构成与时间识别假设。研究不检验相对剥夺、地位焦虑、社会比较或社会凝聚等个体机制；其主要贡献是揭示国家层面关联的含义随估计量、暴露不确定性处理、结局测量、样本构成及时间识别假设而变化。
《Frontiers in Public Health》刊发的该项研究以2000–2021年170个国家为对象，关注国家层面可支配收入不平等（以基尼系数/Gini衡量）与自杀死亡率的生态关联。既有文献中，收入不平等与精神健康、自杀死亡率的关联方向并不一致；更严重的方法是学问题在于，汇总式跨国回归大多被国家间持久差异主导，而国家内年度变化所提供的信息很小。自杀作为罕见死亡事件，还受酒精环境、致死手段可及性、宗教与法律规范、死亡登记完整度等影响，因此生态研究很容易把测量误差、制度差异和缓慢共同趋势误读为机制性关联。研究人员开展本研究的核心目的不是证明“不平等导致自杀增多或减少”，而是把估计量、暴露不确定性、结局测量、样本构成和时间识别假设这些方法学因素系统化，说明同一组全球数据在不同设定下会给出怎样不同的国家层面关联。
为开展研究，研究人员使用WHO全球健康观察站MH_12年龄标准化男女合计自杀死亡率，匹配SWIID 9.92版100条已完成可支配收入Gini序列；协变量取自世界银行世界发展指标（World Development Indicators），包括按购买力平价计人均GDP（取对数）、失业率、城镇化率与互联网使用率。主要技术方法包括：双向固定效应（two-way fixed effects，TWFE）模型以消除国家不变与年度共同冲击；Mundlak/组内组间模型分离between-country与within-country估计量；把SWIID全部100条已完成序列通过Rubin规则与Barnard–Rubin自由度做多重插补（multiple imputation，MI）合并；人口加权、对数人口加权、国家特定线性趋势、WHO区域?年固定效应、一阶差分与动态固定效应等时间识别敏感性；以WHO死亡率数据库重叠子集、GBD 2021自伤结局、区间宽度与SWIID标准误分层做测量质量与报告子集分析；并做观测面板模拟与±0.010对数率单位的最小感兴趣效应（smallest effect size of interest，SESOI）等效检验。
3.1 面板覆盖与识别变异：研究人员描述3,257个国家?年面板，发现对数自杀率约95.4%方差来自国家均值间差异，Gini国内标准差仅1.44点，说明汇总估计主要比较长期自杀水平不同的国家，而非年度国家政策变化。
3.2 汇总、国家间与国家内估计：研究人员得到调整后汇总OLS接近零且不显著，混合模型给出正向但不显著的国家间系数，以及负向国家内系数；未加权TWFE完全调整系数为b = ?0.0157（p = 0.018），SWIID MI为b = ?0.0083（p = 0.057）。结论是国家内变动中，某国比自身通常水平更不平等的年份，自杀对数率略低，但幅度很小。
3.3 暴露不确定性与时间识别敏感性：研究人员比较后验均值代入与MI，MI使系数缩小约47%，但100条序列系数全为负。国家趋势、区域?年、一阶差分、动态模型系数更小且不显著，分数缺失信息（fraction of missing information，FMI）升高。模拟显示Rubin合并在某些非共轭情景下会衰减真实斜率，因此研究人员把?0.0157至?0.0083作为假设依赖的敏感性区间，而非正式边界。
3.4 面板构成、替代不平等指标与影响诊断：人口加权与对数人口加权未逆转负向点估计；报道重叠子集、较窄WHO区间、较小SWIID标准误、较高国内Gini变异子集更负。研究人员指出这既可能是测量误差减小，也可能是高收入与欧洲国家及更完整登记体系的选择效应。170个国家逐一剔除的后验均值系数均为负，区域剔除则区间变宽。
3.4节之后的结局稳健性显示，GBD 2021自伤结局同样为负向；在1,739个死亡率数据库重叠国家?年中，WHO建模结局与登记结局都更负。研究人员强调WHO与GBD共享登记输入，因此符号一致不是独立复制。
3.6 异质性与探索性背景分析：女性与男性斜率都负，男性点估计更负但性别交互不显著；收入组分斜率无统计显著差异；Hofstede个人主义不调节国家内斜率。结论是没有稳定文化调节证据。
讨论部分指出，研究提供的是小规模负向国内生态信号，既不是“无关联”也不是“不平等具有保护性”。汇总与国家间估计回答的是不同问题；高质量子集更负可能来自衰减减少，也可能来自发展、区域与统计能力的选择。研究人员强调，Gini不是可操作治疗，不能把国家层面系数翻译成个体心理获益或支持扩大不平等；公共卫生意义首先在于方法学：政策简报与比较研究应说明所用变异来自国家间还是国家内，并报告暴露不确定性、结局构建、样本构成与时间假设下的敏感性。未来研究需要把宏观分配变化联系到个体感知不平等、地位焦虑、流动预期、信任与求助行为，并用税收、转移支付、最低工资等政策准实验识别总效应与经不平等中介的效应。
结论部分译为：在2000–2021年170个国家中，数据支持可支配收入不平等与自杀死亡率之间存在小规模负向国家内生态关联，但其大小不精确且依赖假设。主要后验均值与SWIID MI估计介于b = ?0.0157至?0.0083之间；人口加权保留负向点估计，而趋势调整与短期模型给出更小且不精确的结果。更高质量与报告子集更负，但强烈选择于收入、区域与统计能力。研究不检验相对剥夺、地位焦虑、社会比较或社会凝聚等个体机制；这些只是可能的背景性解释，需个体层或多层数据验证。主要贡献在于表明：国家层面关联的含义会随估计量、暴露不确定性处理、结局测量、样本构成及时间识别假设而改变。
