Het studieprotocol is beoordeeld en goedgekeurd door de Ethische Commissie van het General Hospital of Northern Theater Command (Goedkeuringsnr. Ethics Y (2026) 75). Hoewel de klinische gegevens afkomstig waren van patiënten die tussen januari 2024 en december 2025 werden behandeld, werd de huidige studie uitgevoerd als een retrospectieve analyse van bestaande klinische dossiers. De in 2026 verkregen ethische goedkeuring omvatte de retrospectieve data-extractie, de-identificatie, analyse en publicatie van deze eerder verzamelde klinische gegevens. Er is geen prospectieve interventie of patiëntenwerving uitgevoerd vóór de ethische goedkeuring. De vereiste voor schriftelijke geïnformeerde toestemming is komen te vervallen vanwege het retrospectieve observationele onderzoeksontwerp. Alle procedures zijn uitgevoerd in overeenstemming met de Verklaring van Helsinki, en alle klinische gegevens zijn gede-identificeerd vóór de analyse. De huidige analyse is uitgevoerd als onderdeel van een goedgekeurd retrospectief neuro-interventioneel klinisch dataconductieonderzoek. De materialen en apparatuur die nodig zijn voor de hieronder beschreven procedures zijn samengevat in de Tabel met Materialen.
Patiëntenwerving en CIE-bepaling
Potentieel geschikte patiënten werden geïdentificeerd door het elektronisch medisch dossier-systeem te doorzoeken naar patiënten die tussen januari 2024 en december 2025 werden behandeld in de afdeling neurochirurgie van het General Hospital of Northern Theater Command. De initiële zoekopdracht werd beperkt op basis van opnameDatum, afdeling, diagnose van cerebrovasculaire ziekte en verslagen van neuro-interventionele procedures. De gedetailleerde workflow voor de screening van patiënten en de toewijzing aan de cohort is geïllustreerd in Figuur 2. Patiënten kwamen in aanmerking voor inclusie indien zij voldeden aan alle volgende criteria: (1) een leeftijd tussen 30 en 80 jaar; (2) diagnose van een cerebrovasculaire ziekte die neuro-interventionele evaluatie of behandeling vereist; (3) beschikbaarheid van volledige klinische dossiers en laboratoriumuitslagen; en (4) een postoperatieve CT-scan van het hoofd, uitgevoerd binnen 3 dagen na de operatie.
Patiënten werden uitgesloten als zij incomplete klinische gegevens of beeldvormingsgegevens hadden, geen endovasculaire interventie ondergingen, ernstige preoperatieve neurologische defecten vertoonden (modified Rankin Scale score ≥4), acute cerebrale infarcten vertoonden op preoperatieve diffusie-gewogen beeldvorming, chronische nierziekte in stadium G4 of hoger hadden volgens de Kidney Disease: Improving Global Outcomes (KDIGO) criteria, of ernstige hart- en vaatziekten hadden. De in- en uitsluitingscriteria werden opeenvolgend toegepast. Na de initiële elektronische zoekopdracht werd elk kandidaatrecord handmatig beoordeeld om de geschiktheid te bevestigen en, indien van toepassing, de reden voor uitsluiting te documenteren. Na screening werden geschikte patiënten in een verhouding van 8:2 willekeurig toegewezen aan trainings- en testcohorten, bestaande uit een trainingscohort (n = 128) en een testcohort (n = 3). Gestratificeerde steekproefname werd gebruikt om een consistente verdeling van CIE- en non-CIE-gevallen tussen de twee cohorten te behouden.
De diagnose van CIE was primair gebaseerd op de temporele associatie tussen het begin van de symptomen en de cerebrovasculaire interventieprocedures, in combinatie met uitsluitende beeldvormende beoordeling. Twee onafhankelijke specialisten in cerebrovasculaire interventies bekeken de postoperatieve klinische dossiers, het tijdstip van het begin van de symptomen, de neurologische manifestaties en de postoperatieve beeldvormingsresultaten voor alle vermoedelijke gevallen. Onenigheid werd via overleg opgelost totdat er consensus was bereikt. Patiënten werden als compatibel met CIE beschouwd indien er binnen 24 u na voltooiing van de neurointerventieprocedure sprake was van nieuw ontstane corticale blindheid, een veranderde bewustzijnsstaat (inclusief somnolentie, delirium of coma), insulten of focale neurologische uitval zoals hemiparese of afasie.
Na het optreden van de symptomen werden alternatieve oorzaken die soortgelijke neurologische manifestaties kunnen veroorzaken uitgesloten door postoperatieve CT- of MRI-beelden te beoordelen in combinatie met het klinische beloop. Postoperatieve CT van het hoofd werd uitgevoerd volgens het institutionele standaardprotocol voor cranial CT zonder contrast. De essentiële acquisitieparameters waren een buispanning van 80 kVp, een buisstroom van 260 mA, een plakdikte van 0,5 mm en standaard axiale reconstructie. De beelden werden door ervaren clinici beoordeeld op hersen- en botvensters. Bij de differentiaaldiagnose werden acute intracraniële bloedingen, subarachnoïdale bloedingen, intracerebrale bloedingen, nieuw ontstane cerebrale infarcten over een groot gebied, aanvallen-gerelateerde veranderingen, infecties en metabole stoornissen overwogen. CT-onderzoeken na het optreden van de symptomen toonden doorgaans focaal of diffuus cerebraal oedeem en hyperdense schaduwen in de cortex, subcorticale regio's of de subarachnoïdale ruimte, die hemorragische laesies nabootsten. MRI-bevindingen tonen gewoonlijk corticale zwelling, in het bijzonder van de temporopariëtale occipitale cortex. Bij sommige patiënten vertoonden CT- of MRI-onderzoeken echter geen duidelijke abnormaliteiten, ondanks compatibele klinische manifestaties.
Perioperatieve variabelenverwerving en constructie van CGR
Algemene demografische gegevens, waaronder leeftijd, geslacht en lichaamsgewicht, werden geëxtraheerd uit gestructureerde velden in het elektronisch medisch dossiersysteem. Variabelen betreffende de klinische geschiedenis, waaronder rookgeschiedenis, alcoholconsumptie, hypertensie, diabetes mellitus en coronaire hartziekten, werden geëxtraheerd uit opnameverslagen, dossiers van de medische voorgeschiedenis en ontslagdiagnoses, en werden gecontroleerd op consistentie. Laboratoriumgegevens werden geëxtraheerd uit het laboratoriuminformatiesysteem. Voor iedere patiënt werd het eerste veneuze bloedmonster gebruikt dat na opname en vóór de neurointerventionele procedure was afgenomen. De geëxtraheerde laboratoriumvariabelen omvatten serumcreatinine, de geschatte glomerulaire filtratiesnelheid (eGFR), totaalcholesterol en triglyceridengehalten. De eGFR werd berekend met de op creatinine gebaseerde vergelijking van de Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) en werd verkregen uit het laboratoriuminformatiesysteem van het ziekenhuis.
Omdat procedurele complexiteit, blootstelling aan contrastmiddelen en de locatie van de laesie de disruptie van de bloed-hersenbarrière en de retentie van contrastmiddelen tijdens neurointerventionele procedures kunnen beïnvloeden, werden perioperatieve procedurele variabelen ook in het voorspellende kader opgenomen. Procedure-gerelateerde variabelen omvatten laesielocatie, type procedure, procedureduur, type contrastmiddel en het totale contrastvolume. De laesielocatie en het type procedure werden bepaald op basis van operatieverslagen en angiografische gegevens. De procedureduur werd gedefinieerd als de tijd vanaf de arteriële punctie tot de voltooiing van de neurointerventionele procedure. Het totale contrastvolume werd gedefinieerd als het cumulatieve volume aan gejodeerd contrastmiddel dat vanaf het begin tot het einde van de procedure werd toegediend. De in deze studie gebruikte contrastmiddelen waren niet-ionische ICA, waaronder iodixanol (320 mg jodium/mL), toegediend volgens de institutionele klinische praktijk. Alle in deze studie opgenomen neurointerventionele procedures werden onder begeleiding van digitale subtractie-angiografie uitgevoerd door een ervaren neurointerventioneel team in het studiecentrum. Alle procedures werden uitgevoerd via de transfemorale arteriële benadering. Tijdens de interventie kregen patiënten intraoperatieve antistolling met heparine, waarbij de geactiveerde stollingstijd binnen het streefbereik van 250–30 s werd gehouden. De vitale functies werden gedurende de gehele procedure continu gemonitord.
Therapeutische interventies, waaronder aneurysma-embolisatie, angioplastiek en stentplaatsing, werden uitgevoerd op basis van de kenmerken van de laesie en het oordeel van de operator. Geschikte microkatheters, coils, stents of ballonnen werden geselecteerd op basis van de procedurele vereisten. Tijdens de angiografische beeldvorming en interventie werden contrastmiddelen toegediend via een hogedrukinjectie of een handmatige infusie met een constante snelheid. Het type contrastmiddel en het totale contrastvolume werden uit de procedureverslagen geëxtraheerd en, indien beschikbaar, gecontroleerd aan de hand van de anesthesie- of verpleegkundige dossiers. Verslagen met inconsistente of onvolledige informatie over het contrastvolume werden handmatig beoordeeld voordat ze werden opgenomen in de uiteindelijke analytische dataset.
Alle procedure-gerelateerde variabelen, inclusief de duur van de procedure, het type therapeutische interventie, het gebruik van apparatuur en de registraties van de toediening van contrastmiddelen, werden systematisch beoordeeld en geverifieerd op volledigheid en consistentie vóór de analyse. Om de onbalans tussen de contrastmiddellast en de renale klaringcapaciteit kwantitatief te evalueren, werd de CGR vastgesteld als een kernvoorspellende variabele en berekend volgens de volgende vergelijking:
(1)
Het totale contrastvolume werd geregistreerd in mL, en de eGFR werd geregistreerd in mL/min/1,73 m2. De eGFR-waarde die werd gebruikt voor de berekening van de CGR was de preoperatieve eGFR-waarde, verkregen uit de eerste veneuze bloedafname na opname en voorafgaand aan de neurointerventionele procedure. Voor elke patiënt werd de CGR berekend na verificatie van zowel het totale contrastvolume als de eGFR. De teller was het totale contrastvolume in mL en de noemer was de eGFR in mL/min/1,73 m2. Dezelfde berekeningsregel werd op alle patiënten toegepast vóór de ontwikkeling van het model.
Kenmerkselectie en ML-workflow
Om de analytische robuustheid te waarborgen, zijn alle klinische en laboratoriumvariabelen vóór de analyse onderworpen aan systematische kwaliteitscontroleprocedures. Voor elke kandidaatvariabele is het aandeel ontbrekende gegevens berekend voorafgaand aan de modelontwikkeling. In de definitieve analytische dataset werden geen ontbrekende waarden waargenomen bij de opgenomen predictoren of uitkomstlabels; daarom is geen enkele variabele uitgesloten vanwege ontbrekende gegevens en was multiple improvisatie niet vereist. Een gedetailleerd overzicht van de ontbrekende gegevens is opgenomen in Aanvullende Tabel 1.
Alle potentiële voorspellende variabelen werden aanvankelijk geëvalueerd met univariaat analyses om hun associaties met het optreden van CIE te beoordelen. Statistische methoden werden geselecteerd op basis van de distributiekenmerken van de gegevens. Variabelen met een normale verdeling werden uitgedrukt als gemiddelde ± standaarddeviatie en vergeleken met de t-toets voor onafhankelijke steekproeven. Variabelen met een niet-normale verdeling werden uitgedrukt als mediaan (interkwartielafstand, IQR) en geanalyseerd met de Mann–Whitney U-toets. Categorische variabelen werden gepresenteerd als frequenties en percentages en geanalyseerd met de chi-kwadraattoets of de exacte toets van Fisher. Statistische significantie werd gedefinieerd als p < 0.05.
De gehele dataset werd vóór de modelontwikkeling in een verhouding van 8:2 willekeurig verdeeld in trainings- en testcohorten. Om overfitting te minimaliseren en multicollineariteit te verminderen, werd least absolute shrinkage and selection operator (LASSO)-regressie in combinatie met 10-voudige kruisvalidatie toegepast voor de selectie van kenmerken. De selectie van kenmerken, de afstelling van hyperparameters en de modelontwikkeling werden uitsluitend uitgevoerd met het trainingscohort om datalekken te voorkomen. Het interne testcohort werd niet gebruikt tijdens de schatting van imputatieparameters, de selectie van kenmerken, de afstelling van hyperparameters of de modeltraining, en werd slechts één keer gebruikt voor de uiteindelijke prestatie-evaluatie. De optimalisatie van hyperparameters werd uitgevoerd met een grid search-strategie met 10-voudige kruisvalidatie binnen het trainingscohort. Specifiek werd het trainingscohort willekeurig verdeeld in 10 wederzijds exclusieve subsets. Tijdens elke iteratie werden negen subsets gebruikt voor de modeltraining en de overgebleven subset voor validatie. Dit proces werd 10 keer herhaald om ervoor te zorgen dat elke subset precies één keer als validatiecohort diende.
Alle computationele analyses werden uitgevoerd in een Python-omgeving. De ontwikkeling en evaluatie van de ML-modellen werden uitgevoerd met scikit-learn. De analysescrips die werden gebruikt voor datapreprocessing, kenmerkselectie, modeltraining, prestatie-evaluatie en DCA werden verstrekt als aanvullende codebestanden. Voor de modelimplementatie werden de geselecteerde kenmerken uit de LASSO-regressie gebruikt als inputvariabelen voor alle kandidaat-classificaties. Continue variabelen werden gestandaardiseerd met z-score normalisatie binnen de trainingscohort, en dezelfde schalingsparameters werden toegepast op de interne testcohort. Categorische variabelen werden gecodeerd volgens het vooraf gedefinieerde coderingsschema; de variabeledefinities en het coderingsschema zijn opgenomen in Aanvullende Tabel 2. Er werden vijf classificaties geïmplementeerd: Gaussian Naive Bayes, SVM, KNN, LightGBM en MLP. Voor SVM werden het kerneltype, de regularisatieparameter C en de kernelcoëfficiënt gamma afgestemd. Voor KNN werden het aantal buren en de afstandswegingsstrategie afgestemd. Voor LightGBM werden het aantal estimators, de leerfrequentie (learning rate), de maximale boomdiepte en het aantal bladeren afgestemd. Voor MLP werden de structuur van de verborgen laag, de activatiefunctie, de regularisatieparameter en de leerfrequentie afgestemd. Gaussian Naive Bayes werd geïmplementeerd met de gespecificeerde variance smoothing-instelling. De optimale hyperparametercombinatie voor elk model werd geselecteerd op basis van de cross-validatieprestaties binnen de trainingscohort. De uiteindelijke geoptimaliseerde hyperparametercombinaties die voor elk model zijn geselecteerd, zijn samengevat in Aanvullende Tabel 3. Na selectie van de hyperparameters werd elk definitief model opnieuw gefit op de volledige trainingscohort en één keer geëvalueerd in de interne testcohort.
Vergelijkingen van modelprestaties
Op basis van de geselecteerde predictoren werden vijf ML-modellen ontwikkeld en vergeleken. De prestaties van de modellen werden geëvalueerd in zowel de trainings- als de testcohort met behulp van de oppervlakte onder de receiver operating characteristic-curve (AUC), het 95% betrouwbaarheidsinterval (95% BI), sensitiviteit, specificiteit en nauwkeurigheid. Er werden receiver operating characteristic-curves gegenereerd om de discriminatieve prestaties van de verschillende modellen te vergelijken. De optimale classificatiedrempel voor elk model werd uitsluitend binnen de trainingscohort bepaald met behulp van de Youden-index. De afgeleide drempelwaarde werd vervolgens vastgesteld en ongewijzigd toegepast op de interne testcohort voor het berekenen van de drempelafhankelijke prestatiegegevens.
Decision curve-analyse
DCA werd uitgevoerd door het netto voordeel van elk model te berekenen over een reeks drempelwaarschijnlijkheden en de modelgebaseerde strategie te vergelijken met twee standaardstrategieën: alle patiënten behandelen en geen enkele patiënt behandelen. DCA werd gekozen omdat het een gelijktijdige beoordeling van de discriminatieprestaties en de potentiële klinische toepasbaarheid in perioperatieve besluitvormingssituaties mogelijk maakt. Om de reproduceerbaarheid te vergemakkelijken, werd de volledige workflow in de volgende volgorde uitgevoerd: patiëntidentificatie uit het elektronisch medisch dossier, toepassing van in- en uitsluitingscriteria, CIE-beoordeling door twee onafhankelijke specialisten, extractie en kruiscontrole van demografische, klinische geschiedenis, laboratorium- en proceduregerelateerde variabelen, CGR-berekening, beoordeling van ontbrekende gegevens en splitsing van de trainings- en testcohort, LASSO-gebaseerde kenmerkselectie binnen de trainingscohort, hyperparametertuning via 10-voudige kruisvalidatie, finale modelherfitting in de volledige trainingscohort, evaluatie van de interne testcohort, ROC-gebaseerde prestatiebeoordeling, bepaling van de drempelwaarde via de Youden-index en DCA.