SCAR dynamics of adolescent substance use: peer influence, dropout, and bifurcation structure in a school-based model

技术摘要：青少年物质使用的 SCAR 动力学

问题陈述
青少年物质使用是一个复杂的公共卫生问题，受同伴影响、保护性社会环境以及与学校相关的转型（如脱离与回归）所驱动。虽然存在针对成年人群或一般成瘾动力学的数学模型，但专门针对青少年、基于学校环境的模型却相对较少，因为在这些环境中，同伴效应极强，且教育转型（辍学、毕业、回归）直接改变了有效人口规模。现有模型往往无法同时兼顾行为传染（起始、升级）与特定于学校的人口更替机制。本文旨在建立一个数学上易于处理的框架，将同伴驱动的行为转型与学校脱离及回归动力学相结合，以理解这些相互作用的机制如何产生性质迥异的学校层面结果。

方法论
作者开发了一个四组分易感–偶然–成瘾–抵抗（Susceptible–Casual–Addicted–Resistant, SCAR）模型，该模型利用适用于高中环境的常微分方程（ODEs）构建。人口被分为：

• $S(t)$： 易感学生（非使用者，具有行为脆弱性）。
• $C(t)$： 偶然或实验性使用者。
• $A(t)$： 具有持续性、问题性或达到物质使用障碍（SUD）水平参与的学生。
• $R(t)$： 抵抗型学生（处于强保护环境中的非使用者）。

该模型整合了以下内容：

1. 人口更替： 恒定的进入率 ($\mu$) 和退出率（毕业/离开），新进入者被分类为易感型 ($1-\rho$) 或抵抗型 ($\rho$)。
2. 同伴影响： 频率依赖型传播项 ($\beta XY/N$)，代表风险促进性影响（从 $C$ 或 $A$ 向 $S$ 的起始；从 $C$ 向 $A$ 的升级）和保护性影响（从 $A$ 向 $C$ 或 $C$ 向 $S$ 的降级；从 $S$ 向 $R$ 的抵抗力构建）。
3. 脱离与回归： 一个速率 $\delta$，描述了问题使用者 ($A$) 脱离学校的情况。一部分比例 $\phi$ 的学生会重新进入学校人口（回归至 $S$），而其余部分 ($1-\phi$) 则永久离开模拟系统。

分析采用了定性动力系统理论，包括：

• 解的正性和有界性证明。
• 推导出一个描述比例演化 ($x, y, z, w$) 的缩放系统，该系统在特定机制下独立于总人口规模。
• 平衡态分析，包括无物质使用、无 SUD 及内部（共存）状态。
• 推导入侵阈值（用于偶然使用的 $R_y$ 以及用于问题使用的 $R_z^{(c)}$）。
• 单参数分叉分析，用以研究平衡稳定性及多稳态的结构性变化。

核心贡献与结果

1. 由回归率 ($\phi$) 定义的结构机制：
   模型通过基于回归参数 $\phi$ 的关键分叉识别了不同的结构机制：

   • $\phi = 1$（守恒人口）： 总学校人口保持不变。缩放系统中的正平衡对应于原始变量中的真实流行平衡，允许所有行为类别长期共存。
   • $\phi < 1$（人口流失）： 问题性使用会导致净人口损失 ($dN/dt < 0$)。在此机制下，缩放系统的正平衡代表“形状平衡”（比例），但不对应于原始变量中的真实流行平衡，因为总人口最终会下降。

2. 起始与升级的截然不同的阈值机制：
   分析表明，首次使用和向问题性使用的进展受不同的动力学阈值控制：

   • $R_y$（偶然使用入侵）： 决定了偶然使用者能否入侵无物质使用的环境。它取决于起始压力 ($\beta_{CS}$)、保护性回归 ($\beta_{RC}$) 以及基准抵抗人口 ($w_f$)。
   • $R_z^{(c)}$（问题使用入侵）： 决定了在已存在偶然使用的情况下，问题使用者能否入侵。它取决于升级压力 ($\beta_{AC}$)、移除率 ($\mu + \delta$) 以及保护性降级 ($\beta_{RA}$)。
   • 启示： 防止起始与防止升级是两个不同的干预问题；抑制其中之一并不自动意味着可以抑制另一个。

3. 多稳态与双稳态：
   模型表现出多稳态特性，特别是存在于无物质使用平衡与稳定的内部高使用平衡之间的双稳态。

   • 无 SUD 子系统： 即使没有问题性使用，系统也可能收敛至无物质使用状态、稳定的偶然使用状态，或者根据起始与保护性影响之间的平衡表现出双稳态。
   • 全系统： 数值示例表明，全系统可以支持无物质使用状态与所有类别共存的稳定内部状态之间的双稳态。最终结果取决于初始条件（吸引域）以及参数值。

4. 政策杠杆的分叉分析：
   单参数分叉图（针对 $\phi=1$）将参数分为三类功能类别：

   • 风险促进型： 增加 $\beta_{AC}$（升级）、$\beta_{AS}$（问题使用者的招募）或 $\beta_{CS}$（偶然使用者的招募）会扩大稳定高使用平衡的区域。值得注意的是，$\beta_{CS}$ 起到了结构组织者的作用，因为偶然使用者既招募易感者，又为未来的问题使用者提供了来源池。
   • 保护/恢复型： 增加 $\beta_{RA}$、$\beta_{RC}$、$\beta_{RS}$ 或 $\rho$（抵抗型进入）会缩小或消除稳定的内部高使用分支，从而扩大无物质使用状态的吸引域。
   • 背景/结构型： 脱离率 $\delta$ 和更替率 $\mu$ 具有复杂的影响。虽然增加 $\delta$ 在数学上会削弱高使用分支，但作者提醒，这可能反映的是辍学而非康复，强调了该模型区分“通过干预移除”与“通过脱离移除”的能力。

意义与主张
本文声称其主要贡献在于对青少年物质使用的定性动力学特征描述，而非经验估计。其意义在于证明了学校层面的结果不仅受入侵阈值的支配，还受以下因素影响：

• 升级动力学： 从偶然使用转向问题性使用的独特机制。
• 多稳态： 存在多个稳定状态（低使用 vs 高使用）的可能性，这意味着初始条件和历史背景（例如早期使用集群）决定了长期结果。
• 回归结构： 脱离是否可逆 ($\phi$) 的关键作用。

作者认为，有效的干预需要一种分层策略，涵盖起始、升级、保护性影响以及成功的回归。他们强调，政策的目标不仅应是降低平均流行率，还应使系统远离支持持续高使用结果的参数机制和状态空间区域。该框架被认为适用于各种物质（如大麻、酒精、尼古丁），在这些领域中，同伴影响和学校转型是核心要素，尽管本文是以青少年大麻使用为具体的应用动机。论文承认了局限性，包括同质混合的假设以及对回归动力学的简化，并建议未来的工作应纳入网络结构和纵向数据校准。