韩国禁用高剧毒农药对自杀死亡率的长期有效性：一项为期12年的准实验研究

目的（Aims）：农药自毒（pesticide self-poisoning）占全球自杀死亡的14%–20%，主要集中在高剧毒农药（Highly Hazardous Pesticides, HHPs）仍可获取的农业地区。尽管韩国拥有先进的医疗体系，这仍构成重大挑战。对此，韩国自2011年11月起分阶段禁用含百草枯（paraquat）在内的HHPs。本研究旨在评估该禁令对农药自杀死亡率的长期影响。方法（Methods）：研究人员采用中断时间序列（Interrupted Time Series, ITS）研究设计，使用2004–2023年全国死亡数据，以带自回归积分移动平均误差（Autoregressive Integrated Moving Average errors, ARIMA errors）的模型考虑自相关性。敏感性分析包括替代干预时点、时期定义及结局指标。主要结局为农药相关自杀的月死亡例数。结果（Results）：研究期间记录268,869例自杀死亡，其中34,962例（13.0%）涉及农药中毒。农药自杀月均计数由禁令前244例降至禁令后81例。ITS分析显示禁令后无即时的水平变化（level change: −0.204; 95% CI: −0.989 to 0.581; P = 0.61），但前36个月下降加速显著（初始斜率改变 initial slope change: −0.122; 95% CI: −0.163 to −0.082; P < 0.001）。反事实（counterfactual）预测估算避免死亡10,846例（95% CI: −567 to 29,715），相当于模型估算减少48%。上述效应在所有敏感性分析中均稳健。结构性断点（structural breakpoint）分析识别到燃气中毒自杀在2008年8月上升，早于禁令三年余，提示无方法替代（method substitution）证据。农药与燃气中毒自杀的人群人口学特征存在显著差异。结论（Conclusions）：韩国HHPs禁令与随后12年间农药自杀死亡率的持续下降相关，且未见方法替代证据。研究结果支持通过农药法规实施手段限制（means restriction through pesticide regulation）是有效的自杀预防策略。

论文解读：韩国禁用高剧毒农药对自杀死亡率长期影响的中断时间序列准实验研究

该研究发表于《Epidemiology and Psychiatric Sciences》。全球每年自杀死亡超70万例，农药自毒（pesticide self-poisoning）在低和中低收入国家尤为突出，占全球自杀死亡约14%–20%；高剧毒农药（Highly Hazardous Pesticides, HHPs）致死率高且无特效解毒剂。手段限制（Means Restriction）——限制获取高致死性自杀手段——是循证自杀预防策略，已有国际证据表明国家级HHPs禁用可显著降低农药相关自杀且不出现完全的方法替代（method substitution）。韩国农业发达且农药易得，曾属OECD中自杀率最高国家之一，历史上农药摄入是常见自杀方式，尤以农村为多。2011年11月韩国启动百草枯（paraquat，致死率>70%、无解毒剂、占韩国农药自杀死亡多数）销售限制，2012年10月全面禁用，同期亦禁用数种有机磷及硫丹等WHO I–II级高毒杀虫剂。既往评估韩国该政策的研究随访仅至2013年或采用Joinpoint回归无法充分排除同期混杂，且对长期效应与方法替代的刻画不足。因此，研究人员开展此项长达12年随访（2004–2023）的准实验研究，以中断时间序列（Interrupted Time Series, ITS）结合差分法（Difference-in-Differences, DiD）探索性分析，检验HHPs禁令对农药自杀死亡率的影响并探查方法替代可能。研究结论为：韩国HHPs禁令与农药自杀死亡率持续12年下降相关，表现为禁令后前36个月下降斜率显著变陡而非即时水平骤降（符合禁售逐步消耗库存机制），反事实预测约避免一万余例死亡，断点分析显示燃气中毒上升早在2008年且与禁令无关，人口学特征差异大，无方法替代证据，支持农药法规手段限制是综合自杀预防的有效组成部分。

主要关键技术方法

研究人员使用韩国统计厅（Statistics Korea）微观数据综合服务（Microdata Integrated Service, MDIS）2004年1月至2023年12月全境≥15岁居民自杀死亡记录（ICD-10：X60–X84），主要结局为农药自杀（X68）月计数，次为上吊（X70）及总自杀，非农药自杀=总自杀-农药自杀。干预时点定为2011年12月（2011年11月23日宣布停续注册后首完整月）。主要分析为带ARIMA误差的ITS回归（Box–Cox优化平方根变换λ≈0.358稳定方差），建模基期趋势、干预后水平改变（level change）及初期（前36个月）与后期斜率改变（slope change），ARIMA阶数经AICc选定含季节成分（周期=12）。估算避免死亡数（deaths averted）为观测值与将干预系数置零的反事实预测之差。敏感性分析含变更干预时点（2012年5月、2012年12月）、变更初期时段长度（6/12/24/48/60个月）、两阶段干预模型（2011年12月销售限制为第一阶段、2012年11月全面禁令为第二阶段）、仅用干预前数据向前投影。方法替代探查采用Bai–Perron结构性断点分析（gas poisoning, X67）及农药与燃气中毒死者人口学χ²比较；探索性DiD以非农药或上吊自杀为对照。亚组按性别、年龄（15–29/30–49/50–64/≥65岁）、地区（首都圈vs非首都圈）、职业（农林业渔业AFF vs其他）分层。所有分析用R 4.4.2中forecast包完成。

研究结果（Results）

描述性统计（Descriptive statistics）：240个月内记录268,869例自杀死亡；干预前（2004年1月–2011年11月）106,386例（均值1,120/月），含农药自杀23,183例（均值244/月，占21.8%）；干预后（2011年12月–2023年12月）162,483例（均值1,121/月），含农药自杀11,779例（均值81/月，占7.2%），绝对降幅67%。禁令前三大方式为上吊（48.9%）、农药（21.8%）、高坠（14.1%）；禁令后变为上吊（50.9%）、高坠（16.4%）、燃气中毒（14.5%）、农药降至第四（7.2%）。

主要分析：中断时间序列（Primary analysis: ITS）：ITS显示干预后无显著即时水平改变（−0.204; 95% CI: −0.989 to 0.581; P = 0.61），但初期（前36个月）斜率改变显著（−0.122; 95% CI: −0.163 to −0.082; P < 0.001），第37个月起后期斜率改变不显著（0.012; 95% CI: −0.012 to 0.036; P = 0.33）。块Bootstrap置信区间较宽（−0.115; 95% CI: −0.214 to 0.006）。反事实比较估算避免农药自杀死亡10,846例（95% Bootstrap CI: −567 to 29,715），相当于较反事实减少47.9%（模型校正既有下降趋势后的保守估计）。年龄-性别标准化率分析结论一致。

亚组分析（Subgroup analyses）：男女均显著下降（男−0.100，女−0.071）；≥65岁降幅最大（−0.107），其次50–64岁（−0.060）、30–49岁（−0.046），15–29岁无显著变化（0.006; P = 0.55）；非首都圈降幅大于首都圈（分别为−0.114、−0.048）；农林业渔业从业者（AFF）、其他职业及职业未指定者均显著下降。

敏感性分析（Sensitivity analyses）：以年龄-性别标准化率替代计数的初期斜率改变仍显著（−0.005; 95% CI: −0.007 to −0.003）。延迟干预时点使效应量递减，以2012年12月（+12个月）为干预点时效应不再显著。不同初期时段长度（6–60个月）均检出显著下降。对数变换方向一致（月降低2.3%，36个月累计约57%）。仅用干预前拟合向前投影估算避免死亡11,473例（95% CI: 793–26,700），与主要分析吻合。两阶段干预模型中，第一阶段（2011年12月销售限制）斜率改变显著（−0.223; P = 0.011），第二阶段（2012年11月全面禁令）无额外显著水平或斜率改变，提示初期限制驱动主要效应，全面禁令起维持作用。

补充分析（Supplementary analyses）：燃气中毒（X67）Bai–Perron断点分析识别首个位点在2008年8月（月均从33升至141），早于禁令逾三年，与知名艺人木炭燃烧自杀事件（2008年10月）相近，表明燃气中毒上升独立于农药禁令。农药自杀死者54.8%为≥65岁，燃气中毒仅6.8%（Cramer's V = 0.55）；男性占比农药68.2% vs燃气81.4%（Cramer's V = 0.15）；地区分布相似。总自杀ITS见显著初期与后期斜率改变但无水平跳变。探索性DiD比较农药与上吊/非农药自杀虽违反平行趋势假定（农药禁令前已更快下降，反映农业人口缩减等），但差异显著，提供辅助证据——因事前趋势偏向压低估计，真实政策效应可能更大。同期其他自杀预防措施（2011年3月《自杀预防法》）无农药特异性，DiD框架可部分调整影响所有方式的因子，不太会解释农药特异性下降。

讨论与结论翻译（Conclusions）

韩国HHPs禁令与随后12年间农药相关自杀死亡率的持续下降相关联，表现为禁令实施后前三年下降加速而无即时水平跳变（符合停售后库存渐耗机制），反事实预测大量死亡得以避免；结构性断点分析显示燃气中毒上升早发生于禁令三年前且两类方法人群特征迥异，不支持方法替代；亚组中老年、农村及农业从业者降幅更明显。局限性含生态研究不能做个体因果推断、DiD平行趋势违反（但方向提示保守估计）、同期综合自杀预防措施不影响农药特异性发现。研究最长随访韩国百草枯禁令，结合ITS与DiD及多项敏感性分析结论稳健。综上所述：The ban on HHPs in South Korea was associated with a sustained decline in pesticide suicide mortality over 12 years, with no evidence of method substitution. These findings support means restriction through pesticide regulation as an effective suicide prevention strategy.（韩国禁用高剧毒农药与12年间农药自杀死亡率的持续下降相关，且无方法替代证据，支持通过农药法规实施手段限制是有效的自杀预防策略。）