《Discover Public Health》：Cooling related household expenditure proxies and heatstroke mortality in an ecological panel study of Japanese prefectures from 2013 to 2022
编辑推荐：
背景：日本中暑（heatstroke，ICD-10 T67）死亡集中于≥65岁成人。冷却相关家庭电力支出（cooling-related household electricity expenditure）是否修饰温度—死亡关系，在不同代理操作化下仍不确定。方
背景：日本中暑（heatstroke，ICD-10 T67）死亡集中于≥65岁成人。冷却相关家庭电力支出（cooling-related household electricity expenditure）是否修饰温度—死亡关系，在不同代理操作化下仍不确定。方法：研究人员纳入47个日本都道府县生态面板，2013—2022年（n=470）。研究人员测试夏季—肩季电力支出比（summer–shoulder electricity expenditure ratio，SSR’），并测试来自家庭收支调查（Family Income and Expenditure Survey）的三个正向代理（碳酸饮料、冰淇淋、矿泉水）和两个负向代理（燃气、煤油）支出，采用泊松回归（Poisson regression）与双向固定效应（two-way fixed effects）、都道府县聚类稳健标准误（prefecture-clustered robust standard errors），并以≥65岁人口为偏移项（offset）。五个替代热指标用于敏感性分析。结果：温度异常（temperature anomaly）具有较大主效应（IRR=1.325，p<0.001），在所有热指标下稳健。主要温度×SSR’交互作用不显著（IRR=1.031，p=0.252）；碳酸饮料达到名义显著性（IRR=1.062，p=0.019），六个代理中五个方向一致。2020年固定效应不显著（IRR=1.10，p=0.39），尽管该时期温度异常最高（+1.57 °C），排除2020年的估计与全样本一致，表明极端热而非COVID-19驱动2020年死亡峰。以热带夜计数（tropical-night counts，最低温度≥25 °C的天数）替代日间峰值温度，与SSR’-电力产生显著负交互（IRR=0.917，p=0.017），提示冷却相关支出特异减弱夜间冷却失败对中暑死亡的影响。结论：分别通过碳酸饮料支出和热带夜暴露捕捉的行为与夜间热维度，值得作为中暑预防目标开展个体水平研究。
该研究关注气候变化背景下热事件频率与严重度上升，以及老年人热相关疾病和死亡风险不成比例的问题。日本中暑（heatstroke，ICD-10 T67）死亡构成增长的公共卫生负担，2018年和2020年等高温年死亡超过1500例，且集中于≥65岁成人。既往贝叶斯时空分析显示日本都道府县中暑死亡存在明显空间自相关，老年人口比例是主要驱动因素，东北和四国农村地区风险高于仅由温度预期的水平；全国公共卫生基础设施相对均一，但T67死亡在都道府县间差异不能完全由温度差异解释。空调（air conditioning，AC）在个体水平研究中是主要保护行为，跨国纵向研究估计AC普及解释日本、美国、加拿大和西班牙热相关死亡风险衰减的17%—20%，美国历史死亡数据、住院研究和气候适应综述提供一致证据。然而，日本都道府县层面AC使用代理有限；家庭收支调查中的冷却相关家庭电力支出（cooling-related household electricity expenditure）是年度可得的都道府县代理，但电力代理会混杂非冷却夏季用电，并可能对操作化选择敏感。温度变异与冷却强度在塑造中暑死亡中的交互作用尚未在全国都道府县纵向框架中正式检验，也未评估替代代理操作化下的稳健性；理解该交互对老龄化社会优先配置热适应干预十分重要。日本47个都道府县10年观测为面板生态推断提供自然实验。研究提出三个问题：RQ1，SSR’是否修饰温度异常与T67死亡关联；RQ2，长期护理保险（long-term care insurance，LTCI）严重度分布即轻度至重度比（mild-to-severe ratio，MSR）是否独立修饰温度—死亡关系；RQ3，观察到的修饰是否对冷却代理操作化稳健。RQ3因同行评审反馈指出单一支出代理验证脆弱而事后加入，并通过理论多代理框架检验：三个正控制（碳酸饮料、冰淇淋、矿泉水，均为独立于冷却基础设施的行为热响应指标）和两个负控制（燃气、煤油，均为冬季取暖标志，预期呈反向关系）。研究人员开展回顾性生态面板研究，以都道府县—年为分析单元，覆盖2013—2022年47个都道府县共470个观测。结论是：主要电力代理在都道府县聚类稳健推断下未显著修饰温度—中暑死亡关联；探索性多代理分析提示碳酸饮料支出和热带夜维度可能分别捕捉行为与夜间热维度；2020年死亡峰由极端热而非COVID-19独立冲击解释。该研究为热适应干预提供了基于代理的生态证据，并强调需个体水平冷却使用和每日热暴露数据确认。
关键技术方法概括（≤250字）：研究人员使用厚生劳动省中暑死亡公开计数、家庭收支调查e-Stat月度支出、LTCI状态报告和SSDSE-B人口估计；分析单元为都道府县—年，2013—2022年47个都道府县470个观测，共同完整样本n=460（Shiga县大津市因老年分层人口缺失全期排除）。主要采用双向固定效应泊松回归（Poisson regression with two-way fixed effects）、都道府县聚类稳健标准误（prefecture-clustered robust standard errors），以≥65岁人口为偏移项（offset）；构建SSR’及五个辅助支出比进行多代理正负控制；敏感性分析包括负二项（negative binomial）、HC1、剔除LTCI、仅都道府县固定效应、排除2020年及S6—S10五种替代热指标；分析使用Python 3.12 statsmodels 0.14.2。
研究结果：3.1描述性统计显示，T67死亡呈极端过离散（方差/均值58.9），2018年1577例、2020年1522例为高峰；SSR’范围0.70—1.36，均值0.99，SD=0.11；SSR’-电力与都道府县纬度负相关。3.2主要分析（RQ1）显示，温度×SSR’-电力交互为正但不显著（IRR=1.031，95% CI 0.979—1.085，p=0.252）；温度主效应大且显著（IRR=1.325，95% CI 1.232—1.424，p<0.001）；敏感性S2—S5方向与幅度类似，排除设定依赖。3.3多代理佐证分析（RQ3）显示，碳酸饮料交互达到名义显著（IRR=1.062，95% CI 1.010—1.116，p=0.019）；燃气和煤油IRR分别为0.998和0.974；六个代理中五个方向与理论预期一致；原文将这一模式视为描述性和假设生成性。3.4次要分析显示，MSR交互不显著（IRR=1.021，p=0.371），联合模型中两个交互均不显著，加入老年人口协变量后估计实质不变。3.5年份固定效应与2020年观察显示，相对2013年，2018年（IRR=1.382，p<0.001）、2019年（IRR=1.364，p<0.001）和2022年（IRR=1.703，p<0.001）残差年效应显著；2020年虽温度异常最高（+1.57 °C，39/47都道府县>1 °C）且T67死亡第二高（1522例），年固定效应不显著（IRR=1.101，p=0.392）。排除2020年的S5主交互（IRR=1.034，p=0.206，n=414）与全样本（IRR=1.031，p=0.252，n=460）实质相同，表明2020年死亡峰主要由极端热暴露解释，而非COVID-19独立年冲击；2018、2019、2022年残差可能反映湿度、湿球黑球温度（wet-bulb globe temperature，WBGT）、热浪时机或年度公共卫生沟通等未测因素。3.6敏感性分析显示，主交互方向与幅度在HC1、剔除LTCI、排除2020等设定下保持；替代热指标中温度主效应均显著（S6 IRR=1.232，S7=1.356，S8=1.260，S9=1.087，p=0.032，S10=1.250，p<0.001）；交互项S6—S9不显著，S10热带夜（tropical nights，日最低温≥25 °C天数）×SSR’-电力为显著负交互（IRR=0.917，95% CI 0.853—0.984，p=0.017），为五个替代热指标中唯一显著交互，提示较高SSR’都道府县中热带夜频率与T67死亡斜率更浅，且与碳酸饮料交互方向相反。
讨论部分总结：主要发现是主电力代理交互方向为正但未达常规显著，不能解释为个体冷却保护；多代理框架为事后扩展，正控制接近预期、碳酸饮料名义显著、负控制低于1，六中五方向一致，但原文强调其为描述性和假设生成性，未进行正式多重比较。碳酸饮料与电力支出不同，反映与冷却基础设施正交的人群行为热响应；辅助代理显著而主电力代理不显著提示生态代理可能不同捕捉热响应行为、残余热脆弱或真正冷却相关修饰，生态数据无法区分。2020年死亡峰方面，原文讨论社会隔离、口罩影响体温调节、户外活动减少影响热适应、医疗系统过载等可能机制，但条件于温度异常、都道府县、SSR’和LTCI率后2020年固定效应不显著，排除2020年主估计又与全样本几乎相同，两条证据收敛于极端热而非COVID-19驱动2020年T67死亡峰。热带夜信号为敏感性分析结果，原文将其解释为描述性而非确证性，并提出一种假设生成性解释：都道府县冷却相关电力支出至少部分捕捉持续夜间空调使用，从而减弱夜间冷却失败对中暑死亡的影响；碳酸饮料支出可能反映日间行为热响应。文献比较中，主电力代理零结果与Sera等、Barreca等、Ostro等个体AC保护证据不同，但与Gronlund等所述生态AC代理常未能检测个体保护的混合记录一致；Zeng等、Fujibe、Honda等、Ng等、Oka等、Ohashi和Iwamoto、Ma等日本研究支持都道府县热死亡关系并非直接比例于环境温度，冷却强度可能是修饰因素之一。MSR未修饰可能因高严重度LTCI用户不成比例机构化、冷却提供集中管制，以及MSR时间不变操作化限制组内变异信息。局限包括SSR’仅覆盖二人及以上家庭、Shiga县大津市老年人口缺失、MSR采用FY2022时间不变值、生态谬误（ecological fallacy）、2020年结构断点、缺乏个体AC拥有和使用数据、T67公开计数未按年龄分层、主要热暴露为夏季峰值温度且未纳入WBGT、未建模滞后结构及城市热岛效应等。
研究结论翻译：在这项2013—2022年47个日本都道府县生态面板研究中，主要基于电力的代理在都道府县聚类稳健推断下未显著修饰温度—中暑死亡关联。探索性多代理分析提示某些夏季支出比，尤其碳酸饮料支出，可能捕捉热响应行为或残余热脆弱而非冷却保护本身；这些发现应解释为基于代理的生态效应修饰，并需个体水平冷却使用和每日热暴露数据确认。次要公共卫生相关发现涉及2020年死亡峰：两条收敛证据——2020年固定效应不显著（IRR=1.10，p=0.39）和排除2020年敏感性S5（IRR=1.034，p=0.206，n=414）与全样本主估计（IRR=1.031，p=0.252，n=460）几乎相同——表明2020年日本T67死亡峰主要由当年极端热（平均温度异常+1.57 °C，39/47都道府县超过+1 °C）解释，而非独立COVID-19可归因风险冲击。相反，2018、2019、2022年显著残差年固定效应（IRR=1.38、1.36、1.70；均p<0.001）表明温度不完全解释这些高温年逐年变化，指向湿度、WBGT、热浪时机或年度公共卫生沟通等未测因素作为未来研究目标。第三个重要发现来自替代热指标敏感性分析：五个替代热暴露操作化中，仅热带夜计数（日最低温≥25 °C天数）产生显著温度×SSR’-电力交互且为负（IRR=0.917，95% CI 0.853—0.984，p=0.017）。这提示都道府县冷却相关支出特异性减弱夜间冷却失败对T67死亡的影响，构成独立于日间峰值温度的热暴露维度。尽管该发现为描述性和假设生成性，其指向夜间热暴露作为公共卫生适应目标，需在老龄化日本政策规划中与日间热保护并列考虑。
