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

Resumen Técnico: Dinámica SCAR del Consumo de Sustancias en Adolescentes

Planteamiento del Problema
El consumo de sustancias en adolescentes es un problema de salud pública complejo impulsado por la influencia de los pares, entornos sociales protectores y transiciones relacionadas con la escuela, como la desvinculación y el reingreso. Si bien existen modelos matemáticos para poblaciones adultas y dinámicas generales de la adicción, comparativamente pocos están adaptados a poblaciones adolescentes y escolares, donde los efectos de los pares son potentes y las transiciones educativas (abandono escolar, graduación, retorno) alteran directamente la población efectiva. Los modelos existentes a menudo no logran contabilizar simultáneamente el contagio conductual (iniciación, escalada) y los mecanismos de rotación específicos de la escuela. Este documento aborda la necesidad de un marco matemáticamente tratable que integre las transiciones conductuales impulsadas por pares con las dinámicas de desvinculación y reingreso escolar para comprender cómo estos mecanismos interactuantes generan resultados cualitativamente diferentes a nivel escolar.

Metodología
Los autores desarrollan un modelo de cuatro compartimentos Susceptible–Casual–Adicto–Resistente (SCAR) utilizando ecuaciones diferenciales ordinarias (EDO) adaptadas a un entorno de escuela secundaria. La población se divide en:

• $S(t)$: Estudiantes susceptibles (no usuarios, vulnerables conductualmente).
• $C(t)$: Usuarios casuales o experimentales.
• $A(t)$: Estudiantes con una implicación sostenida, problemática o de nivel de Trastorno por Consumo de Sustancias (TCS).
• $R(t)$: Estudiantes resistentes (no usuarios en entornos fuertemente protectores).

El modelo incorpora:

1. Rotación Demográfica: Entrada constante ($\mu$) y salida (graduación/partida) mediante tasas, con nuevos entrantes clasificados como susceptibles ($1-\rho$) o resistentes ($\rho$).
2. Influencia de Pares: Términos de transmisión de frecuencia dependiente ($\beta XY/N$) que representan la influencia que promueve el riesgo (iniciación desde $C$ o $A$ hacia $S$; escalada de $C$ a $A$) y la influencia protectora (desescalada de $A$ a $C$ o de $C$ a $S$; construcción de resistencia de $S$ a $R$).
3. Desvinculación y Reingreso: Una tasa $\delta$ a la cual los usuarios problemáticos ($A$) se desvinculan de la escuela. Una fracción $\phi$ de estos estudiantes reingresa a la población escolar (regresando a $S$), mientras que el resto ($1-\phi$) abandona el sistema modelado permanentemente.

El análisis emplea la teoría de sistemas dinámicos cualitativos, que incluye:

• Prueba de positividad y acotación de las soluciones.
• Derivación de un sistema escalado que describe la evolución de las proporciones ($x, y, z, w$) independiente del tamaño total de la población en ciertos regímenes.
• Análisis de equilibrio para estados libres de sustancias, libres de TCS e interiores (coexistencia).
• Derivación de umbrales de invasión ($R_y$ para el uso casual, $R_z^{(c)}$ para el uso problemático).
• Análisis de bifurcación de un solo parámetro para examinar cambios estructurales en la estabilidad del equilibrio y la multistabilidad.

Contribuciones Clave y Resultados

1. Regímenes Estructurales Definidos por el Reingreso ($\phi$):
   El modelo identifica una bifurcación crítica basada en el parámetro de retorno $\phi$.

• $\phi = 1$ (Población Conservada): La población escolar total es constante. Los equilibrios positivos en el sistema escalado corresponden a verdaderos equilibrios endémicos en las variables originales, permitiendo la coexistencia a largo plazo de todas las clases conductuales.
• $\phi < 1$ (Pérdida de Población): El uso problemático induce una pérdida neta de población ($dN/dt < 0$). En este régimen, los equilibrios positivos del sistema escalado representan "equilibrios de forma" (proporciones) pero no corresponden a verdaderos equilibrios endémicos en las variables originales, ya que la población total eventualmente disminuye.

2. Mecanismos de Umbral Distintos para la Iniciación y la Escalada:
   El análisis revela que el primer uso y la progresión hacia el uso problemático están gobernados por umbrales dinámicos separados:

• $R_y$ (Invasión del Uso Casual): Determina si los usuarios casuales pueden invadir un entorno libre de sustancias. Depende de la presión de iniciación ($\beta_{CS}$), el retorno protector ($\beta_{RC}$) y la población resistente base ($w_f$).
• $R_z^{(c)}$ (Invasión del Uso Problemático): Determina si los usuarios problemáticos pueden invadir un entorno donde el uso casual ya está establecido. Depende de la presión de escalada ($\beta_{AC}$), las tasas de eliminación ($\mu + \delta$) y la desescalada protectora ($\beta_{RA}$).
• Implicación: Prevenir la iniciación y prevenir la escalada son problemas de intervención distintos; suprimir uno no suprime automáticamente el otro.

3. Multiestabilidad y Bistabilidad:
   El modelo exhibe multiestabilidad, específicamente la bistabilidad entre un equilibrio libre de sustancias y un equilibrio interior estable de alto consumo.

• Subsistema Libre de TCS: Incluso sin uso problemático, el sistema puede converger a un estado libre de sustancias, a un estado estable de uso casual, o exhibir bistabilidad entre ambos, dependiendo del equilibrio entre la iniciación y la influencia protectora.
• Sistema Completo: Los ejemplos numéricos demuestran que el sistema completo puede soportar la bistabilidad entre un estado libre de sustancias y un estado interior estable donde todas las clases coexisten. El resultado final depende de las condiciones iniciales (atractor de cuenca) así como de los valores de los parámetros.

4. Análisis de Bifurcación de Palancas de Política:
   Los diagramos de bifurcación de un solo parámetro (para $\phi=1$) categorizan los parámetros en tres clases funcionales:

• Promotores de Riesgo: Aumentar $\beta_{AC}$ (escalada), $\beta_{AS}$ (reclutamiento por usuarios problemáticos) o $\beta_{CS}$ (reclutamiento por usuarios casuales) expande la región de equilibrios estables de alto uso. Notablemente, $\beta_{CS}$ actúa como un organizador estructural, ya que los usuarios casuales tanto reclutan susceptibles como suministran el grupo para futuros usos problemáticos.
• Protectores/Recuperación: Aumentar $\beta_{RA}$, $\beta_{RC}$, $\beta_{RS}$ o $\rho$ (entrada de resistentes) reduce o elimina la rama estable de alto uso, ampliando la cuenca de atracción para el estado libre de sustancias.
• Contextuales/Estructurales: La tasa de desvinculación $\delta$ y la rotación $\mu$ tienen efectos complejos. Si bien aumentar $\delta$ debilita matemáticamente la rama de alto uso, los autores advierten que esto puede reflejar el abandono escolar en lugar de la recuperación, destacando la capacidad del modelo para distinguir entre la eliminación vía intervención versus la eliminación vía desvinculación.

Significancia y Reivindicaciones
El artículo afirma que su contribución principal es una caracterización dinámica cualitativa del consumo de sustancias en adolescentes, más que una estimación empírica. La significancia radica en demostrar que los resultados a nivel escolar están gobernados no solo por umbrales de invasión, sino también por:

• Dinámicas de Escalada: El mecanismo distinto de pasar del uso casual al uso problemático.
• Multiestabilidad: La existencia de múltiples estados estables (bajo uso vs. alto uso) para conjuntos de parámetros idénticos, lo que implica que las condiciones iniciales y el contexto histórico (por ejemplo, grupos de uso temprano) determinan los resultados a largo plazo.
• Estructura de Reingreso: El papel crítico de si la desvinculación es reversible ($\phi$).

Los autores argumentan que una intervención eficaz requiere una estrategia estratificada que aborde la iniciación, la escalación, la influencia protectora y el reingreso exitoso. Enfatizan que las políticas deben tener como objetivo no solo reducir la prevalencia promedio, sino mantener el sistema alejado de los regímenes de parámetros y las regiones del espacio de estados que soportan resultados de alto uso sostenido. El marco se presenta como aplicable a diversas sustancias (cannabis, alcohol, nicotina) donde la influencia de los pares y las transiciones escolares son centrales, aunque la aplicación específica que motiva el modelo es el consumo de cannabis en adolescentes. El artículo reconoce las limitaciones, incluyendo la suposición de mezcla homogénea y la simplificación de las dinámicas de reingreso, sugiriendo que trabajos futuros incorporen estructuras de red y calibración de datos longitudinales.