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

Technische Samenvatting: SCAR-dynamiek van adolescent middelengebruik

Probleemstelling
Adolescent middelengebruik is een complex volksgezondheidsprobleem dat wordt gedreven door peer-invloed, beschermende sociale omgevingen en schoolgerelateerde transities zoals desengagement en herintreding. Hoewel er wiskundige modellen bestaan voor volwassen populaties en algemene verslavingsdynamiek, zijn er relatief weinig modellen die specifiek zijn afgestemd op adolescenten en schoolpopulaties, waar peer-effecten krachtig zijn en educatieve transities (uitval, afstuderen, terugkeer) de effectieve populatie direct veranderen. Bestaande modellen falen vaak in het gelijktijdig rekening houden met gedragsmatige besmetting (initiatie, escalatie) en schoolspecifieke omloopsmechanismen. Dit artikel adresseert de behoefte aan een wiskundig hanteerbaar kader dat peer-gestuurde gedragsmatige transities integreert met schoolse desengagement en herintredingsdynamiek om te begrijpen hoe deze interagerende mechanismen kwalitatief verschillende schoolniveau-uitkomsten genereren.

Methodologie
De auteurs ontwikkelen een vier-compartimenten Susceptibel–Casual–Addicted–Resistant (SCAR) model met behulp van gewone differentiaalvergelijkingen (ODE's) die zijn aangepast aan een middelbare schoolomgeving. De populatie wordt verdeeld in:

• $S(t)$: Susceptibele studenten (niet-gebruikers, gedragsmatig kwetsbaar).
• $C(t)$: Casual of experimentele gebruikers.
• $A(t)$: Studenten met aanhoudend, problematisch of op SUD-niveau (Substance Use Disorder) betrokkenheid.
• $R(t)$: Resistente studenten (niet-gebruikers in sterke beschermende omgevingen).

Het model incorporeert:

1. Demografische Omloop: Constante instroom ($\mu$) en uitstroom (afstuderen/vertrek) ratio's, waarbij nieuwe instroom wordt geclassificeerd als susceptibel ($1-\rho$) of resistent ($\rho$).
2. Peer-invloed: Frequentie-afhankelijke transmissietermen ($\beta XY/N$) die risicoverhogende invloed vertegenwoordigen (initiatie van $C$ of $A$ naar $S$; escalatie van $C$ naar $A$) en beschermende invloed (de-escalatie van $A$ naar $C$ of $C$ naar $S$; opbouw van resistentie van $S$ naar $R$).
3. Desengagement en Herintreding: Een ratio $\delta$ waarbij problematische gebruikers ($A$) zich loskoppelen van de school. Een fractie $\phi$ van deze studenten keert terug in de schoolpopulatie (terugkeer naar $S$), terwijl de rest ($1-\phi$) het gemodelleerde systeem permanent verlaat.

De analyse maakt gebruik van kwalitatieve dynamische systeemtheorie, waaronder:

• Bewijs van positiviteit en begrensdheid van oplossingen.
• Afleiding van een geschaald systeem dat de evolutie van proporties ($x, y, z, w$) beschrijft, onafhankelijk van de totale populatieomvang in bepaalde regimes.
• Evenwichtsanalyse voor middelvrije, SUD-vrije en interieur (coexistentie) toestanden.
• Afleiding van invasiedrempels ($R_y$ voor casual gebruik, $R_z^{(c)}$ voor problematisch gebruik).
• Eén-parameter bifurcatieanalyse om structurele veranderingen in evenwichtsstabiliteit en multistabiliteit te onderzoeken.

Belangrijkste Bijdragen en Resultaten

1. Structurele Regimes Gedefinieerd door Herintreding ($\phi$):
   Het model identificeert een kritieke bifurcatie gebaseerd op de terugkeerparameter $\phi$.

• $\phi = 1$ (Geconserveerde Populatie): De totale schoolpopulatie is constant. Positieve evenwichten in het geschaalde systeem komen overeen met echte endemische evenwichten in de originele variabelen, wat de langdurige coexistentie van alle gedragsklassen mogelijk maakt.
• $\phi < 1$ (Populatieverlies): Problematisch gebruik induceert netto populatieverlies ($dN/dt < 0$). In dit regime vertegenwoordigen positieve evenwichten van het geschaalde systeem "vorm-evenwichten" (proporties), maar komen ze niet overeen met echte endemische evenwichten in de originele variabelen, aangezien de totale populatie uiteindelijk zal dalen.

2. Distinctieve Drempelmechanismen voor Initiatie en Escalatie:
   De analyse laat zien dat eerste gebruik en progressie naar problematisch gebruik worden beheerst door afzonderlijke dynamische drempels:

• $R_y$ (Invasie van Casual Gebruik): Bepaalt of casual gebruikers een middelvrije omgeving kunnen binnendringen. Het hangt af van de initiatiedruk ($\beta_{CS}$), de beschermende terugkeer ($\beta_{RC}$) en de basislijn van de resistente populatie ($w_f$).
• $R_z^{(c)}$ (Invasie van Problematisch Gebruik): Bepaalt of problematische gebruikers een omgeving kunnen binnendringen waar casual gebruik reeds gevestigd is. Het hangt af van de escalatiedruk ($\beta_{AC}$), de verwijderingsratio's ($\mu + \delta$) en de beschermende de-escalatie ($\beta_{RA}$).
• Implicatie: Het voorkomen van initiatie en het voorkomen van escalatie zijn afzonderlijke interventieproblemen; het onderdrukken van het één onderdrukt het ander niet automatisch.

3. Multistabiliteit en Bistabiliteit:
   Het model vertoont multistabiliteit, specifiek bistabiliteit tussen een middelvrije evenwichtstoestand en een stabiele interieur toestand met hoog gebruik.

• SUD-Vrij Subsysteem: Zelfs zonder problematisch gebruik kan het systeem convergeren naar een middelvrije staat, een stabiele staat van casual gebruik, of vertonen zij bistabiliteit tussen deze twee, afhankelijk van de balans tussen initiatie en beschermende invloed.
• Volledig Systeem: Numerieke voorbeelden demonstreren dat het volledige systeem bistabiliteit kan ondersteunen tussen een middelvrije staat en een stabiele interieur toestand waarin alle klassen coexisteren. De uiteindelijke uitkomst hangt af van de begincondities (het bassin van attractie) evenals van de parameterwaarden.

4. Bifurcatieanalyse van Beleidsinstrumenten:
   Eén-parameter bifurcatiediagrammen (voor $\phi=1$) categoriseren parameters in drie functionele klassen:

• Risicoverhogend: Het verhogen van $\beta_{AC}$ (escalatie), $\beta_{AS}$ (rekrutering door problematische gebruikers) of $\beta_{CS}$ (rekrutering door casual gebruikers) breidt de regio van stabiele evenwichten met hoog gebruik uit. Opvallend genoeg fungeert $\beta_{CS}$ als een structurele organisator, aangezien casual gebruikers zowel susceptibelen rekruteren als de bron vormen voor de toekomstige pool van problematisch gebruik.
• Beschermend/Herstel: Het verhogen van $\beta_{RA}$, $\beta_{RC}$, $\beta_{RS}$ of $\rho$ (resistente instroom) verkleint of elimineert de stabiele interieur tak met hoog gebruik, waardoor het bassin van attractie voor de middelvrije staat wordt vergroot.
• Contextueel/Structureel: De desengagement-ratio $\delta$ en de omloop $\mu$ hebben complexe effecten. Hoewel het verhogen van $\delta$ wiskundig gezien de tak met hoog gebruik verzwakt, waarschuwen de auteurs dat dit eerder een weerspiegeling is van uitval (dropout) dan van herstel, wat de bekwaamheid van het model benadrukt om onderscheid te maken tussen verwijdering via interventie versus verwijdering via desengagement.

Significantie en Claims
Het artikel claimt dat de primaire bijdrage een kwalitatieve dynamische karakterisering van adolescent middelengebruik is, in plaats van een empirische schatting. De significantie ligt in het aantonen dat uitkomsten op schoolniveau niet alleen worden beheerst door invasiedrempels, maar ook door:

• Escalatie-dynamiek: Het afzonderlijke mechanisme van beweging van casual naar problematisch gebruik.
• Multistabiliteit: Het bestaan van meerdere stabiele toestanden (laag gebruik versus hoog gebruik) voor identieke parameter sets, wat impliceert dat begincondities en historische context (bijv. clusters van vroeg gebruik) de langetermijnuitkomsten bepalen.
• Herintredingsstructuur: De cruciale rol van de vraag of desengagement omkeerbaar is ($\phi$).

De auteurs stellen dat effectieve interventie een gelaagde strategie vereist die zowel initiatie als escalatie, beschermende invloed en succesvolle herintreding aanpakt. Ze benadrukken dat beleid er niet alleen op moet gericht zijn het gemiddelde prevalentieniveau te verlagen, maar ook om het systeem weg te houden van de parameterregimes en de toestandsruimten die aanhoudend hoog gebruik ondersteunen. Het kader wordt gepresenteerd als toepasbaar op diverse middelen (cannabis, alcohol, nicotine) waar peer-invloed en schooltransities centraal staan, hoewel de specifieke motivatie voor het model de adolescenten-cannabisconsumptie is. Het artikel erkent beperkingen, zoals de aanname van homogene menging en de vereenvoudiging van herintredingsdynamiek, en suggereert dat toekomstig werk netwerkstructuren en longitudinale data-calibratie moet bevatten.