《Epidemiology and Psychiatric Sciences》：Community mental health workforce capacity and official early psychiatric readmission indicators in Italy: a longitudinal ecological panel study
编辑推荐：
目的：劳动力人数与早期再入院指标常并列出现在卫生系统仪表盘中，这助长了“更多人员应伴随更少返院”的假设。研究人员检验了前一年医疗-护理人员密度（medical-nursing headcount density）的年度变化是否与意大利官方供方侧指标——即公立精
目的：劳动力人数与早期再入院指标常并列出现在卫生系统仪表盘中，这助长了“更多人员应伴随更少返院”的假设。研究人员检验了前一年医疗-护理人员密度（medical-nursing headcount density）的年度变化是否与意大利官方供方侧指标——即公立精神科病房符合条件的一般精神病学出院后30天内非计划再入院（unplanned readmission within 30 days）相关，并次要检验7天内再入院。
方法：利用国家精神卫生、医院出院和人口数据，研究人员构建了一个2015—2024年意大利21个领土单元的回顾性面板（210个单元-年）。暴露变量为上一年年末每10万成人中医疗和护理人员数。研究人员使用泊松伪极大似然法（Poisson pseudo-maximum likelihood, PPML）对30天再入院计数建模，以符合条件的出院数作为偏移量（offset），并纳入领土单元和年份固定效应。主要推断采用限制性野聚类得分自举法（restricted wild cluster score bootstrap），含9,999次Rademacher抽样和21个聚类。次要和敏感性分析均预先设定，源数据空白未予插补，分析材料和脚本已公开存档。
结果：全部210个观测的结果计数和符合条件的出院分母均完整；暴露变量可得205个观测。每10万成人每增加10名医疗-护理人员，发病率比（incidence rate ratio, IRR）为1.006（95% CI 0.934–1.034；p = 0.8095），与相对变化从降低6.6%到增加3.4%的范围相容。在线性PPML敏感性分析中，IRR范围为0.995至1.042，所有置信区间均包含1。在31次逐一剔除单元或年份的重新估计中，IRR范围为0.974至1.044，每个置信区间均包含1。非线性检验* p* = 0.8572。7天IRR为1.032（95% CI 0.929–1.075；p = 0.6743），而加权线性模型估计为+0.215个百分点（95% CI ?1.159至+0.704；p = 0.7203）。
结论：数据未提供证据表明医疗-护理人员密度的逐年变化与官方30天再入院指标之间存在可检测的纵向生态关联。这并不意味着劳动力不重要：人员数并不能反映全职等效可用性（full-time-equivalent availability）、部署或转介过程，且供方侧设计无法支持患者层面或因果推断。领土政策应将充足资源与其部署、出院后组织以及照护提供所处的社会、经济和健康背景的整体系统视角相结合。
论文解读文章
一、研究背景与问题
精神科早期再入院是精神卫生服务指标中最具直观说服力、也最易被过度解读的指标之一。其在行政数据中的可观测性，使再入院和出院后随访成为国际基准比较的组成部分（Hermann等，2006）。然而，早期返院反映的是入院和出院过程、临床严重程度、既往服务使用、社会境况、随访可得性及地方照护系统等多重因素，既非纯粹的住院质量度量，也非对个体患者再入院概率的估计（Durbin等，2007）。六国欧洲登记研究进一步表明，数据库互操作性和指标定义限制着跨国可比性（Katschnig等，2019）。因此，监测需要明确的分母、时间窗口和分析层级，观察到的变异不应归因于照护的单一组成部分。
劳动力配置为此指标提供了一个看似合理但间接的路径。有效的可得性可支持出院规划、及时随访、用药连续性、信息交换和住院至社区过渡期的危机应对（Vigod等，2013；Donisi等，2016）。医疗和护理人员被选为核心临床劳动力群体，因其在国家精神卫生信息系统中纵向报告，并参与上述过渡功能。意大利提供了特别有信息量的场景：其以社区为导向的系统通过区域结构和精神卫生部门组织，国家系统报告服务资源和精神科住院活动。年度报告将意大利精神卫生信息系统（Sistema Informativo per la Salute Mentale, SISM）的劳动力人数与官方医院出院记录（Scheda di Dimissione Ospedaliera, SDO）早期再入院指标并列呈现，但二者在同一领土内的变化是否同步尚不清楚。本研究检验了滞后医疗-护理人员密度与官方供方侧30天指标之间的纵向生态关联，并将对应的7天指标预设为次要结局。
二、研究设计与主要结论
研究人员开展了一项2015—2024年意大利21个稳定领土单元的回顾性纵向生态面板研究，共210个地区/自治省-年观测。分析单元为地区/自治省-年。研究使用常规收集的公开可得的聚合行政和统计数据，未访问、链接或分析任何个体层面记录。在拟合任何推断模型之前，研究人员锁定了一份统计分析计划（Statistical Analysis Plan, SAP），规定了分析人群、结局和暴露定义、模型层级、源数据空白处理、主要估计量、推断程序以及次要和敏感性分析。锁定计划附有结局定义与可比性备忘录、源数据审计、冻结分析面板和SHA-256数据冻结清单。
主要结局为官方SDO报告的符合条件的一般精神科出院后30天内非计划再入院的年度地区/省层面计数。符合条件的普通精神科出院年数作为描述性指标的分母和回归模型中的偏移量。主要暴露为SISM在t ? 1年12月31日的医疗人员和护士人数，按t年1月1日每10万成人居民计算。确认性关联采用PPML估计，以30天再入院计数为结局，以符合条件的普通精神科出院数的自然对数为偏移量，并纳入地区/省和结局年份固定效应。推断以地区/省层面聚类，形成21个领土聚类。主要小聚类推断采用限制性野聚类得分自举法，在零假设下施加约束，9,999次Rademacher权重抽样，并通过双侧得分检验反转获得95%置信区间。
主要结果为：每10万成人每增加10名医疗-护理人员，30天再入院的IRR为1.006（得分反转95% CI 0.934–1.034；限制性野得分p = 0.8095）。伴随的聚类稳健t(20) 95% CI为0.966–1.048（p = 0.7476）。主要得分区间对应相对变化从?6.6%到+3.4%。因此，数据未提供证据表明滞后医疗-护理人员密度的变化与官方供方侧30天指标之间存在可检测的纵向生态关联。7天结局的IRR为1.032（95% CI 0.929–1.075；p = 0.6743），加权线性固定效应模型估计为+0.215个百分点（95% CI ?1.159至+0.704；p = 0.7203）。所有预设敏感性分析、31次逐一剔除单元或年份的重新估计以及非线性检验均支持同一解释。
三、关键研究方法与技术
本研究采用的主要技术方法包括：第一，构建冻结的纵向生态面板，覆盖2015—2024年意大利21个领土单元（19个地区加2个自治省），共210个单元-年观测，数据来源于意大利精神卫生信息系统年度报告、医院出院记录和意大利国家统计局（Istituto Nazionale di Statistica, ISTAT）居民人口数据；第二，采用泊松伪极大似然法（PPML）作为主要估计量，以符合条件的普通精神科出院数的对数为偏移量，纳入领土单元和年份双向固定效应；第三，使用限制性野聚类得分自举法处理21个小聚类推断问题，以9,999次Rademacher权重抽样和得分检验反转获得置信区间；第四，在建模前锁定统计分析计划（SAP），对源数据进行文档审计，构建SHA-256数据冻结清单，全部分析脚本化并公开存档以保证可审计性；第五，预设了排除特定年份、改变滞后结构、安慰剂时序、加入情境协变量、单元特异性线性趋势、加权线性模型、逐一剔除单元/年份的敏感性分析以及限制性三次样条非线性检验。
四、研究结果
描述性结果：冻结面板包含210个地区/自治省-年观测。结局计数和符合条件的出院分母完整；滞后医疗-护理人员密度可得205个观测。符合条件的普通精神科出院数从2015年的95,530例下降至2020年的71,827例，2024年回升至81,135例。官方30天早期再入院指标在2016年最高（17.69%），2020年降至13.54%，2024年为14.33%。7天指标范围为2017年的7.00%至2023年的8.36%，2024年为8.12%。年度中位滞后医疗-护理人员密度范围为每10万成人32.6（2023年）至36.2（2016年）。
主要确认性分析：PPML分析包含205个观测和21个聚类。每10万成人每增加10名医疗-护理人员，30天再入院的IRR为1.006（得分反转95% CI 0.934–1.034；p = 0.8095），对应相对变化从?6.6%至+3.4%。数据未提供证据表明存在可检测的纵向生态关联。
预设敏感性分析：排除2020年得IRR 1.007（0.930–1.036；p = 0.8180）；排除2020–2021年得IRR 1.002（0.929–1.030；p = 0.9184）；排除2015年得IRR 1.012（0.954–1.052；p = 0.5059）；两年滞后得IRR 0.995（0.937–1.039；p = 0.7250）；领先/安慰剂时序得IRR 1.042（0.937–1.082；p = 0.4682）；情境调整模型得IRR 1.016（0.955–1.055；p = 0.4024）；单元特异性线性趋势得IRR 0.995（0.955–1.035；p = 0.5388）。所有置信区间均包含零值1.00。
次要7天结局：每10万成人每增加10名医疗-护理人员，7天再入院的IRR为1.032（得分反转95% CI 0.929–1.075；p = 0.6743）。加权线性固定效应模型估计30天指标变化为+0.215个百分点（95% CI ?1.159至+0.704；p = 0.7203）。
影响性分析：21次逐一剔除地区/省重新估计中，IRR范围为0.974至1.015（p = 0.4867–0.9957）；10次逐一剔除结局年份重新估计中，IRR范围为0.978至1.044（p = 0.3273–0.9912）。每个得分反转置信区间均包含1。限制性三次样条非线性检验p = 0.8572。
职业特异性次级分析：每10万成人每增加10名人员，医生IRR为1.062（0.859–1.161；p = 0.5501；Holm调整p = 1.0000），护士IRR为0.994（0.901–1.044；p = 0.8318；Holm调整p = 1.0000），精神科医生IRR为1.263（0.901–1.582；p = 0.2510；Holm调整p = 0.7530）。191观测共同样本上，医疗-护理、医生、护士和精神科医生IRR分别为1.005、1.059、0.991和1.263。
描述性事后检查：原始暴露标准差为9.64，双向固定效应后残差标准差为4.33（SD比0.45；残差方差占比20.2%）。单元内中位范围为11.33（IQR 7.34–14.88）。条件最小可检测效应端点为IRR 0.944和1.059。观测和拟合总再入院数均为125,189；杠杆率最高为REG_12在2018年（0.501），该单元-年已包含在影响性重新估计中。
五、讨论与结论
讨论部分指出，早期精神科再入院受出院组织、随访、用药连续性、信息传递、社会支持、危机应对和住院-社区协调共同塑造。过渡性干预是复杂安排而非孤立的投入要素。全职等效人员配置（full-time equivalent staffing）将比人员数更好地考虑合同工作时间和兼职角色，但仍是结构性度量：若无空缺、跨场所和班次部署、技能组合、工作量、人员连续性和照护质量信息，全职等效本身对再入院的预测能力可能仍有限。精神科住院使用同样多维，反映结构性资源、流行病学变化、诊断病例组合、人口构成、自愿和强制路径、临床复杂性、社区照护连续性、社会决定因素和地方组织。
测量挑战并非意大利特有，但国际证据必须在暴露和结局定义的水平上进行比较。六国欧洲登记研究要求实质性协调后才能比较；韩国和日本索赔研究将住院人员配置与再入院关联，但其患者层面设计和服务情境不同于本研究的区域社区劳动力面板。意大利患者层面研究仍针对不同层级。国际相关要点不是可移植的意大利效应估计，而是共同的测量问题：劳动力供给、部署照护和医院使用不能假定在不同司法管辖区或分析层级间等效。
研究结论部分翻译如下：在2015—2024年领土面板中，数据未提供证据表明人员数与供方侧再入院之间存在可检测的纵向生态关联。这并不使劳动力变得无关紧要。它表明名义人员数和供方侧再入院不应替代部署的临床能力、照护连续性或领土绩效。一个更有信息量的政策框架需要充足资源和对人员时间、可及性、患者和病例组合、社会和经济条件、身体医疗保健及精神科过渡过程如何共同组织的多维监测。
