美国社区健康食品零售销售份额（集约边际）与多维健康结局的异质性关联——基于双/去偏机器学习的全国人口普查区证据

《Food Policy》：The intensive margin of food access and its association with health outcomes: evidence of heterogeneous effects across U.S. Communities

　　

研究人员将关注点从外延边际（extensive margin，即食品零售点数量）转向集约边际（intensive margin，以健康食品零售商销售份额衡量），考察本地食品环境与健康的关系。利用大型全国企业数据集并与人口普查区（census tract）层面健康结局合并，研究人员构建了人口普查区食品环境的基于销售额的健康度指数（sales-based healthfulness index），并估计其与居民肥胖、糖尿病、高血压、心理健康问题和抑郁的关系。研究人员应用带交叉拟合（cross-fitting）的双/去偏机器学习（double/debiased machine learning, DML）灵活调整高维社会经济和建成环境协变量，并将集约边际获取与总体零售商密度分离。研究发现，较高的健康食品销售份额与更好的健康结局正相关；健康食品销售份额每增加10个点，与肥胖患病率降低0.87个百分点、糖尿病患病率降低0.19个百分点、高血压患病率降低0.61个百分点相关，心理健康问题和抑郁也有类似改善。这些关联在基线食品环境间存在异质性：心血管代谢结局（cardiometabolic outcomes）在低获取人口普查区反应最强，而心理健康改善在中等和高获取社区更大。

研究背景与问题：美国饮食相关慢性病是紧迫的公共卫生挑战，政策关注本地零售食品环境及居民可 readily、可负担购买食品的营养质量，Make America Healthy Again（MAHA）倡议进一步提升了对食品环境的关注。更健康的食品环境虽被普遍视为理想目标，但其与健康的关联可能因地点和结局不同而不均，理解何处、对谁最重要有助于有效配置有限政策资源。既有研究多以“外延边际”（extensive margin）即食品零售点的数量或存在衡量食品环境，零售点计数给大型超市与小便利店相同权重，不能反映销售量与相对市场存在，即“集约边际”（intensive margin）；计数反映潜在获取（potential access），而非销售指标捕捉的实现市场活动。由此存在两个缺口：本地零售食品市场实际组成及销售指标可能带来的食品环境重分类受关注有限；多数研究估计单一平均关系，对不同社区异质性关注不足。研究人员因此分析食品环境集约边际与肥胖、糖尿病、高血压、心理健康问题、抑郁的关联，并以销售份额重构改良零售食品环境指数（modified Retail Food Environment Index, mRFEI），估计平均关联与异质性。其意义在于把食品获取从门店计数扩展至销售构成，为食品政策评估提供补充信息；但该研究为观察性分析，估计值应解释为稳健关联而非因果效应。

关键技术方法：数据来源包括National Establishment Time-Series（NETS）2019年企业记录以构建2019年销售基于mRFEI，2010年NETS记录构建基线mRFEI作为控制；CDC 500 Cities and Places 2019提供人口普查区健康结局，2010年American Community Survey（ACS）提供社会经济、人口与地理控制变量，共覆盖70,102个census tracts。零售商按North American Industry Classification System（NAICS）分类为健康或不健康类型；mRFEI=100×健康零售商销售额/健康与不健康零售商总销售额。研究人员采用DML与10折交叉拟合，使用Lasso、随机森林、梯度提升和神经网络选择最佳样本外预测，并以部分线性模型残差化估计；进一步采用连续局部DML、广义倾向得分（generalized propensity score, GPS）和核平滑估计剂量—反应关系及条件平均处理效应（conditional average treatment effect, CATE），并按低获取、食物沼泽和健康食品环境分样本分析。

研究结果：
5.1 Baseline DML estimates（基线DML估计）：全样本DML显示，mRFEI对五种健康结局的估计系数均为负且统计显著。肥胖系数为?0.0869，即mRFEI增加10个点与肥胖患病率降低约0.87个百分点相关；糖尿病系数为?0.0193，对应降低约0.19个百分点；高血压系数为?0.0613，对应降低约0.61个百分点；心理健康问题系数为?0.0219，抑郁系数为?0.0216，约对应0.22个百分点降低。结局模型第一阶段交叉拟合R²为0.542—0.681，表明观察到的tract特征共同解释了健康变异的较大份额；控制变量用于剔除混杂变异，不应解释为结构参数。结论是较高健康食品销售份额与更好的心血管代谢和心理健康结局正相关。

5.2 Heterogeneity across baseline food environments（基线食品环境异质性）：按mRFEI≤52.42定义为低获取（low access）、52.4254.09定义为健康食品环境，分样本DML并以两样本z检验比较系数。肥胖在低获取tract的系数为?0.0975，大于健康食品环境的?0.0361，差异显著；高血压在低获取为?0.0910，大于健康食品环境的?0.0368，差异显著；糖尿病在低获取为?0.0163，但三类环境间差异均不显著（p≥0.131）；心理健康问题在低获取tract关联较弱，在食物沼泽和健康食品环境更强；抑郁在低获取tract最强，在食物沼泽较小。结论是异质性具有结局特异性：心血管代谢结局在低获取环境反应最强，心理健康问题在中等和高获取环境更强，抑郁则表现相反排序。

5.3 Nonlinear dose–response and CATE profiles（非线性剂量—反应与CATE曲线）：针对家庭收入低于中位数的tract，连续局部DML估计CATE随预测mRFEI变化的曲线。低收入tract中，肥胖和糖尿病在低获取环境由很低mRFEI向中段上升时改善最明显；高血压在低获取环境于较高预测mRFEI仍有较大负向关联，健康食品环境曲线也呈下降趋势；心理健康问题在食物沼泽和健康食品环境中段显示明显改善；抑郁在低获取环境中低至中段改善最强，在食物沼泽和健康食品环境较高预测mRFEI处趋近于零且较不稳定。极端mRFEI区域不确定性增大，中央区间推断更可靠。结论是关系明显非线性，统一销售份额基准可能忽略社区差异，政策推断需按结局和基线食品环境区分。

讨论部分总结：该研究为观察性分析，因果解释依赖条件独立/无混杂假设，而该假设不能完全验证，因此估计值应解释为稳健条件关联而非因果效应；但DML在较少参数假设下对高维关系去偏，仍对食品政策和公共卫生具有信息价值。食品获取政策常关注社区是否有足够健康食品店，而销售份额提示还应关注零售活动的强度与构成：若健康零售商已存在但销售份额较低，应考察食品选择、价格、交通和居民购物模式等限制；新店投资仍需评估其提供内容、能否吸引持续购买以及如何改变区域销售构成。非线性结果表明统一销售份额基准可能忽略社区基线条件和健康脆弱性的差异。按结局看，肥胖和高血压在低获取tract关联最强，糖尿病方向相同但跨环境差异不显著，心理健康问题在食物沼泽和健康食品环境关联最大，抑郁则在低获取tract最强。政策评估除零售商计数外，应追踪销售在不同零售商类型间的分布，并单独衡量新鲜及其他营养食品是否更可得、更可负担且被居民购买；销售基于mRFEI不能替代个体产品层面消费测量。

研究结论部分翻译：美国饮食相关慢性病仍然常见，相关资源负担巨大；成人肥胖、高血压和糖尿病广泛存在，是医疗利用的重要贡献因素，因此预防相关环境即使发生增量变化也可能具有有意义的总体影响。同时，食品零售获取在空间上仍不均衡，许多社区尤其是农村和低收入地区面临不可忽视的距离与交通约束，附近零售选择的构成可能影响日常购买机会。本研究利用DML分析人口普查区层面食品获取的集约边际及其与不同社区健康结局的差异化关联，并进一步分析低收入census tracts中健康结局的CATE。研究的核心贡献是以销售为基础重新界定“食品获取”：简单门店计数将大型超市和小型便利店视为等价，尽管二者对居民实际可购买食品的贡献不同；集约边际mRFEI通过健康零售商在当地食品销售中的份额捕捉这一区别。基线DML估计与该机制一致：在灵活调整tract特征并将集约边际获取与总体零售商密度分离后，较高的健康食品销售份额与美国本土相连地区观察数据中较低的肥胖、糖尿病和高血压，以及心理健康问题和抑郁改善系统相关。研究解释受观察性研究的标准识别限制；DML改善高维混杂调整并通过交叉拟合减少过拟合偏倚，但不能保证消除与零售商构成和健康均相关的未观察tract层面因素所致偏倚。未来研究可将集约边际测量与准实验变异（如政策诱导的销售激励）结合，并利用面板结构分离tract内变化与时间不变差异。