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

Resumo Técnico: Dinâmica SCAR do Uso de Substâncias na Adolescência

Definição do Problema
O uso de substâncias na adolescência é um problema complexo de saúde pública impulsionado pela influência de pares, ambientes sociais protetivos e transições escolares, como o desengajamento e o reingresso. Embora existam modelos matemáticos para populações adultas e dinâmicas gerais de dependência, comparativamente poucos são adaptados para populações adolescentes em ambiente escolar, onde os efeitos de pares são potentes e as transições educacionais (evasão, graduação, retorno) alteram diretamente a população efetiva. Os modelos existentes frequentemente falham em contabilizar simultaneamente o contágio comportamental (iniciação, escalada) e os mecanismos de rotatividade específicos da escola. Este artigo aborda a necessidade de um arcabouço matematicamente tratável que integre transições comportamentais impulsionadas por pares com a dinâmica de desengajamento e reingresso escolar para compreender como esses mecanismos interativos geram resultados qualitativamente diferentes no nível escolar.

Metodologia
Os autores desenvolvem um modelo Suscetível–Casual–Viciado–Resistente (SCAR) de quatro compartimentos, utilizando equações diferenciais ordinárias (EDOs) adaptadas a um ambiente de ensino médio. A população é dividida em:

• $S(t)$: Estudantes suscetíveis (não usuários, comportamentalmente vulneráveis).
• $C(t)$: Usuários casuais ou experimentais.
• $A(t)$: Estudantes com envolvimento sustentado, problemático ou de nível de Transtorno por Uso de Substância (TUS).
• $R(t)$: Estudantes resistentes (não usuários em ambientes protetivos fortes).

O modelo incorpora:

1. Rotatividade Demográfica: Taxas constantes de entrada ($\mu$) e saída (graduação/partida), com novos entrantes classificados como suscetíveis ($1-\rho$) ou resistentes ($\rho$).
2. Influência de Pares: Termos de transmissão dependentes de frequência ($\beta XY/N$) representando influência promotora de risco (iniciação de $C$ ou $A$ para $S$; escalada de $C$ para $A$) e influência protetiva (desescalada de $A$ para $C$ ou $C$ para $S$; construção de resistência de $S$ para $R$).
3. Desengajamento e Reingresso: Uma taxa $\delta$ na qual usuários problemáticos ($A$) se desengajam da escola. Uma fração $\phi$ desses estudantes reingressa na população escolar (retornando para $S$), enquanto o restante ($1-\phi$) deixa o sistema modelado permanentemente.

A análise utiliza a teoria de sistemas dinâmicos qualitativos, incluindo:

• Prova de positividade e limitação das soluções.
• Derivação de um sistema escalonado que descreve a evolução das proporções ($x, y, z, w$) independente do tamanho total da população em certos regimes.
• Análise de equilíbrio para estados livres de substâncias, livres de TUS e estados interiores (coexistência).
• Derivação de limiares de invasão ($R_y$ para uso casual, $R_z^{(c)}$ para uso problemático).
• Análise de bifurcação de um parâmetro para examinar mudanças estruturais na estabilidade do equilíbrio e na multistabilidade.

Principais Contribuições e Resultados

1. Regimes Estruturais Definidos pelo Reingresso ($\phi$):
   O modelo identifica uma bifurcação crítica baseada no parâmetro de retorno $\phi$.

• $\phi = 1$ (População Conservada): A população escolar total é constante. Equilíbrios positivos no sistema escalonado correspondem a verdadeiros equilíbrios endêmicos nas variáveis originais, permitindo a coexistência de longo prazo de todas as classes comportamentais.
• $\phi < 1$ (Perda de População): O uso problemático induz uma perda líquida de população ($dN/dt < 0$). Neste regime, equilíbrios positivos do sistema escalonado representam "equilíbrios de forma" (proporções), mas não correspondem a verdadeiros equilíbrios endêmicos nas variáveis originais, pois a população total eventualmente declina.

2. Mecanismos de Limiar Distintos para Iniciação e Escalada:
   A análise revela que o primeiro uso e a progressão para o uso problemático são governados por limiares dinâmicos distintos:

• $R_y$ (Invasão do Uso Casual): Determina se usuários casuais podem invadir um ambiente livre de substâncias. Depende da pressão de iniciação ($\beta_{CS}$), do retorno protetivo ($\beta_{RC}$) e da população resistente basal ($w_f$).
• $R_z^{(c)}$ (Invasão do Uso Problemático): Determina se usuários problemáticos podem invadir um ambiente onde o uso casual já está estabelecido. Depende da pressão de escalada ($\beta_{AC}$), das taxas de remoção ($\mu + \delta$) e da desescalada protetiva ($\beta_{RA}$).
• Implicação: Prevenir a iniciação e prevenir a escalada são problemas de intervenção distintos; suprimir um não suprime automaticamente o outro.

3. Multiestabilidade e Bistabilidade:
   O modelo exibe multistabilidade, especificamente a bistabilidade entre um equilíbrio livre de substâncias e um equilíbrio interior estável de alto uso.

• Subsistema Livre de TUS: Mesmo sem o uso problemático, o sistema pode convergir para um estado livre de substâncias, um estado estável de uso casual ou exibir bistabilidade entre eles, dependendo do equilíbrio entre iniciação e influência protetiva.
• Sistema Completo: Exemplos numéricos demonstram que o sistema completo pode suportar a bistabilidade entre um estado livre de substâncias e um estado interior estável onde todas as classes coexistem. O resultado final depende das condições iniciais (atrator/bacia de atração) bem como dos valores dos parâmetros.

4. Análise de Bifurcação de Alavancas de Políticas:
   Diagramas de bifurcação de um parâmetro (para $\phi=1$) categorizam os parâmetros em três classes funcionais:

• Promotores de Risco: O aumento de $\beta_{AC}$ (escalada), $\beta_{AS}$ (recrutamento por usuários problemáticos) ou $\beta_{CS}$ (recrutamento por usuários casuais) expande a região de equilíbrios estáveis de alto uso. Notavelmente, $\beta_{CS}$ atua como um organizador estrutural, pois os usuários casuais tanto recrutam suscetíveis quanto fornecem o reservatório para o futuro uso problemático.
• Protetores/Recuperação: O aumento de $\beta_{RA}$, $\beta_{RC}$, $\beta_{RS}$ ou $\rho$ (entrada de resistentes) reduz ou elimina o ramo estável de alto uso, ampliando a bacia de atração para o estado livre de substâncias.
• Contextuais/Estruturais: A taxa de desengajamento $\delta$ e a rotatividade $\mu$ têm efeitos complexos. Embora aumentar $\delta$ matematicamente enfraqueça o ramo de alto uso, os autores alertam que isso pode refletir evasão em vez de recuperação, destacando a capacidade do modelo de distinguir entre remoção via intervenção versus remoção via desengajamento.

Significância e Alegações
O artigo alega que sua principal contribuição é uma caracterização dinâmica qualitativa do uso de substâncias na adolescência, e não uma estimativa empírica. A significância reside em demonstrar que os resultados ao nível escolar são governados não apenas por limiares de invasão, mas também por:

• Dinâmicas de Escalada: O mecanismo distinto de transição do uso casual para o uso problemático.
• Multiestabilidade: A existência de múltiplos estados estáveis (baixo uso vs. alto uso) para conjuntos de parâmetros idênticos, implicando que as condições iniciais e o contexto histórico (ex: agrupamentos de uso precoce) determinam os resultados de longo prazo.
• Estrutura de Reingresso: O papel crítico de se o desengajamento é reversível ($\phi$).

Os autores argumentam que uma intervenção eficaz requer uma estratégia em camadas que aborde a iniciação, a escalada, a influência protetiva e o reingresso bem-sucedido. Eles enfatizam que as políticas devem visar não apenas reduzir a prevalência média, mas manter o sistema longe de regimes de parâmetros e regiões do espaço de estados que suportam resultados de alto uso sustentado. O arcabouço é apresentado como aplicável a várias substâncias (cannabis, álcool, nicotina) onde a influência de pares e as transições escolares são centrais, embora a aplicação específica motivadora do modelo seja o uso de cannabis na adolescência. O artigo reconhece limitações, incluindo a suposição de mistura homogênea e a simplificação das dinâmicas de reingresso, sugerindo que trabalhos futuros incorporem estruturas de rede e calibração de dados longitudinais.