$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Patienten und Studiendesign
Der erste Datensatz dieser Studie stammt aus der National Health and Nutrition Examination Survey (NHANES) Datenbank für den Zeitraum 2017–2018. Dieser spezielle Zeitraum wurde ausgewählt, weil die Umfragedaten Leberultraschall-Transient-Elastographiemessungen mit FibroScan-Technologie® enthielten. NHANES verwendet ein geschichtetes, mehrstufiges Probabilitätsstichprobensystem und dient als landesweite, bevölkerungsbasierte Umfrage, die alle zwei Jahre durchgeführt wird. Sie sammelt systematisch national repräsentative gesundheitsbezogene Daten zur nicht-institutionalisierten US-Bevölkerung, um den Ernährungs- und Gesundheitszustand der allgemeinen Zivilbevölkerung in den Vereinigten Staaten zu bewerten1. Die NHANES ist eine landesweit repräsentative Querschnittsstudie, die vom National Center for Health Statistics (NCHS) durchgeführt wird. Das Umfrageprotokoll wurde vom NCHS Research Ethics Review Board genehmigt, wobei von allen Teilnehmern dokumentierte informierte Zustimmung eingeholt wurde. Die Studie wurde gemäß den Erklärungen von Helsinki und Istanbul durchgeführt und vom Ethikkomitee des Ersten Angeschlossenen Krankenhauses der Medizinischen Universität Wenzhou (2016–246, 1. Dezember 2016) genehmigt, wobei von jedem Teilnehmer eine schriftliche informierte Zustimmung eingeholt wurde.
Der erste Datensatz bestand aus 5494 Personen aus der NHANES-Umfrage 2017–2018, die sich einer FibroScan-Untersuchung unterzogen hatten. Nach dem untenstehenden Ausschluss wurden insgesamt 2677 Teilnehmer in die Analyse einbezogen, darunter 718 Personen mit NAFLD und 1959 mit Nicht-NAFLD. Diese Teilnehmer wurden dann zufällig in einen Trainingssatz (n = 1785) und einen Testsatz (n = 892) eingeteilt. Der zweite Datensatz umfasste 582 Personen aus der Abteilung für Infektionskrankheiten am First Affiliated Hospital der Wenzhou Medical University (2018–2020). Nach Anwendung derselben Ausschlusskriterien wurden insgesamt 200 Personen eingeschlossen, darunter 159 Personen mit NAFLD und 41 Personen mit Nicht-NAFLD. Diese Kohorte wurde als unabhängiges Validierungsset verwendet. Das Studiendesign wurde entwickelt, um das Modell mit Multi-Center-Daten zu konstruieren und zu validieren, wodurch die Zuverlässigkeit und Generalisierbarkeit der Ergebnisse verbessert werden. Die klinischen Ausgangsmerkmale der NAFLD- und Nicht-NAFLD-Gruppen wurden anhand des Tabellenpakets 1 (Tabelle 1) zusammengefasst. Zusätzlich sind der Prozess der Patientenauswahl und der gesamte Studienablauf in Abbildung 1 dargestellt.
Diagnostische Kriterien und Ausschlusskriterien für NAFLD
Die NAFLD-Diagnose basierte auf folgenden Kriterien:16 Jahre oder älter: 18 Jahre oder älter, Teilnahme an einer transienten Elastographie (FibroScan) mit Alkoholkonsum begrenzt auf ≤140 g/Woche bei Frauen und ≤ 210 g/Woche bei Männern in den vorangegangenen 12 Monaten, ein kontrollierter Dämpfungsparameter (CAP)-Wert von ≥ 302 dB/m, gemessen mit dem FibroScan 502 V2 Touch-System (Echosens, Paris, Frankreich) mit entweder einer mittleren (M) oder extra-großen (XL) Sonde oder einer durch Leberbiopsie bestätigten Diagnose der Abteilung für Infektionskrankheiten des ersten angeschlossenen Krankenhauses der Wenzhou Medical University23.
Die Ausschlusskriterien für NAFLD sind wie folgt dargelegt:hoher Alkoholkonsum (durchschnittliche tägliche Aufnahme > 20 g bei Frauen und > 30 g bei Männern laut NHANES-Alkoholkonsum-Umfrage5), Vorhandensein von Hepatitis B oder C, HIV-Infektion, Autoimmunhepatitis, primäre galliäre Cholangitis, Morbus Wilson, langfristige Anwendung nichtsteroidaler entzündungshemmender Medikamente, Kalziumkanalblocker, Tamoxifen, Amiodaron, Kortikosteroide, Isoniazid oder Methotrexat, Schwangerschaft oder Stillen. und eine Diagnose von Leberkrebs oder einem anderen gutartigen oder bösartigen Tumor.
Datenerhebung und Variablenauswahl
Die potenziellen Prädiktorvariablen, die in dieser Studie enthalten sind, sind unten aufgeführt:
Demografische Merkmale (d. h. Alter und Geschlecht); Body-Mass-Index (BMI); CAP-Werte der Teilnehmer in der NHANES-Datenbank; Allgemeine biochemische Tests [d. h. Albumin (ALB), Globulin (GLO), Gesamtprotein (TP), Laktatdehydrogenase (LDH), Blutharnstoff-Stickstoff (BUN), Harnsäure (UA), Gamma-Glutamyl-Transferase (GGT), Triglyzerid (TG), Serum-Glukose (Glu), Serumkreatinin (SCr), Gesamtbilirubin (TBIL), Natrium (Na⁺), Chlorid (Cl⁻), Kalium (K⁺), Calcium (Ca), Bicarbonat (HCO₃), Gesamtcholesterin (TC), Aspartat-Aminotransferase (AST) und Alanin-Aminotransferase (ALT)]; Standardhämatologische Parameter [d. h. Rote Blutkörperchenzählungen (RBC), Weiße Blutkörperchenzählungen (WBC), Neutrophilzählungen (NEUT), Eosinophilzahlen (EOS), Lymphozytenzählungen (LYM), Monozytenzählungen (MON), Verteilungsbreite der roten Blutkörperchen (RDW) und Thrombozytenzahlen (PLT)]; Vorgeschichte von Bluthochdruck und Diabetes mellitus (DM). Unter den in der Studie eingeschlossenen Probanden wurden die diagnostischen Kriterien für Diabetes und Bluthochdruck aus einer früheren maschinell-learningbasierten Vorhersagemodellstudie mit Fokus auf NAFLD24 gewonnen.
Fehlende Datenverarbeitung und Merkmalsauswahl
Die in dieser Studie verwendeten Kohortendaten enthielten fehlende Werte. Das Ausschließen aller unvollständigen Datensätze würde nicht nur die Stichprobengröße der Analyse verringern, sondern auch die Datenqualität beeinträchtigen und potenziell die Vorhersageergebnisse verzerren. Daher wurden alle Daten mit fehlenden Werten über 20 % ausgeschlossen. Für Datensätze mit fehlenden Werten ≤20 % wurden je nach Datentyp unterschiedliche Imputationsmethoden angewendet: "norm" für stetige Variablen, "logreg" für binäre Klassifikationsvariablen und "polyreg" für Mehrklassenvariablen. Diese Imputationen wurden mit dem "mäus"-Paket in R für Mehrfachimputation25 durchgeführt. In dieser Studie wurden alle stetigen Variablen in binäre Variablen aufgeteilt, wobei optimale Klassifikationsschwellenwerte mittels Analyse der Empfänger-Betriebscharakteristik (ROC) bestimmt wurden. Konkret wurde der Grenzwert, der dem maximalen Youden-Index auf der ROC-Kurve entspricht, als optimales Klassifikationskriterium ausgewählt, wodurch die Klassifikationsleistung optimiert und gleichzeitig ein angemessenes Gleichgewicht zwischen Sensitivität und Spezifität gewahrt wird. Anschließend wurde die teilweise kleinste Quadrat-Diskriminantenanalyse (PLS-DA) eingesetzt, um die Daten effektiv zu gruppieren.
Zur weiteren Auswahl der Funktionen wurden die Personen zufällig in Trainings- und Testgruppen im Verhältnis 7:3 mit dem "Carart"-Paket eingeteilt. Zunächst wurde für die Merkmalsauswahl (LASSO) die Regression mit minimaler absoluter Schrumpfung und Auswahloperator verwendet, wobei 14 Schlüsselvariablen für die spätere Analyse identifiziert wurden. Anschließend wurde ChatGPT-4 angewendet, um jeder Variablen einen Wichtigkeitswert zuzuweisen. Um die Merkmalsbedeutung systematisch zu bewerten und gleichzeitig potenzielle Verzerrungen und stochastische Variationen zu kontrollieren, wurde ein standardisiertes Bewertungsprotokoll implementiert. Der spezifische Prompt für das Modell lautete: "Basierend auf der etablierten klinischen Literatur zur nicht-alkoholischen Fettlebererkrankung (NAFLD) sollten Sie für jede der folgenden 14 Variablen, die durch die LASSO-Regression identifiziert wurden, einen Wichtigkeitswert von 1 (niedrigster) bis 10 (höchster) zuweisen." Indem die Bewertung strikt auf die durch LASSO-Regression vorausgewählten Variablen beschränkt wurde, wurde das Risiko, halluzinierte oder irrelevante Merkmale einzubeziehen, minimiert. Um ein mögliches Datenleck und unbeabsichtigten Einsatz von Outcome-Prävalenz strikt zu verhindern, wurde das Sprachmodell vollständig auf den empirischen Datensatz verblindet. Die Bewertung war darauf beschränkt, bereits vorhandenes medizinisches Wissen über die Variablennamen zu synthetisieren. Dieses Verfahren verbessert traditionelle Methoden, wie die LASSO-Stabilitätsauswahl oder SHAP-basiertes Pruning, indem sichergestellt wird, dass rein datengetriebene Merkmale vor der endgültigen Modellintegration eine robuste pathophysiologische Plausibilität aufweisen. Zusätzlich wurde dieses Verfahren zur Minderung der Single-Response-Varianz zehnmal unabhängig davon iteriert. Der Durchschnittswert für jede Variable wurde über diese Iterationen hinweg berechnet, um eine objektive Priorisierung sicherzustellen. Die Variablen wurden dann in absteigender Reihenfolge basierend auf ihren Durchschnittswerten gerankt. Schließlich wurde die Auswahl auf nur jene Variablen mit einem durchschnittlichen Wichtigkeitswert von über 5 eingegrenzt, was zu 8 Schlüsselvariablen führte: SCr, Harnurinteil, GGT, Glu, Bluthochdruck, Diabetes, TG und BMI.
Entwicklung von AutoML-basierten Vorhersagemodellen für NAFLD
In dieser Studie wurde eine umfassende Suite klassischer und fortschrittlicher maschineller Lernalgorithmen mit H2O-AutoML für die effektive Diagnose von NAFLD integriert. Durch die Nutzung der AutoML-Fähigkeiten der H2O.ai-Plattform wurden maschinelle Lernanalysen für binäre Klassifikationsaufgaben durchgeführt. Zu den verwendeten Algorithmen gehörten eXtreme Gradient Boosting (XGBoost), Gradient Boosting Machine (GBM), Generalized Linear Model (GLM), Extremely Randomized Trees (XRT), Deep Learning (DL) und Stacked Ensemble. Diese Algorithmen wurden systematisch ausgewertet, um das optimale Modell für die Krankheitsdiagnose zu identifizieren. Um eine strenge methodologische Reproduzierbarkeit zu gewährleisten, wurden die Ausführungsparameter von H2O AutoML explizit definiert. Die automatisierte Suche war auf eine maximale Laufzeit von 11.687 Sekunden und maximal 302 Modelle beschränkt, wobei ein fester zufälliger Seed von 13 verwendet wurde. Interne Vorverarbeitungs-Flags beinhalteten die automatische Imputation von Restwerten mit Mittelwert-/Modusalgorithmen sowie die Zielkodierung für kategoriale Variablen mit hoher Kardinalität. Die ausgewählte optimale Architektur (bezeichnet als GBM_grid_1_model77) war eine Gradient-Boosting-Maschine mit folgenden spezifischen Hyperparametern: insgesamt 28 Bäume, eine maximale Baumtiefe von 5 und eine Lernrate von 0,1.
Um die Modellrobustheit zu verbessern und Überanpassungsrisiken zu mindern, wurde ein AutoML-Framework implementiert, das systematische Hyperparameter-Tuning- und Validierungsprotokolle integriert. Der Prozess begann mit einer automatisierten Untersuchung von 200 verschiedenen Modellkonfigurationen durch Hyperparameteroptimierung, wobei fünffache Kreuzvalidierung eingesetzt wurde, bei der der Trainingsdatensatz in fünf gegenseitig ausschließende Teilmengen aufgeteilt wurde. Während des iterativen Trainings nutzte jede Konfiguration vier Teilmengen (80 %) für den Modellbau, während eine Teilmenge (20 %) für die Validierung reserviert wurde, wobei diese Validierungsrolle sequentiell über alle Falten hinweg rotierte. Um die Recheneffizienz mit der Leistungsoptimierung auszubalancieren, wurde dynamisches frühes Stoppen basierend auf der Fläche unter der ROC-Kurve (AUC)-Metrik implementiert. Dieser Mechanismus unterbrach das Training, wenn AUC-Verbesserungen drei aufeinanderfolgende Zyklen unter eine Schwelle von 0,001 fielen, was sowohl auf die individuelle Modelloptimierung als auch auf den gesamten AutoML-Suchprozess angewendet wurde. Die endgültige Modellauswahl priorisierte Konfigurationen, die einen maximalen durchschnittlichen AUC über sowohl Trainings- als auch Validierungssets zeigten, während gleichzeitig eine konsistente Leistung zwischen diesen Sätzen und minimale metrische Varianz über Cross-Validierungsiterationen erforderlich waren. Dieser integrierte Ansatz gewährleistete eine optimale prädiktive Genauigkeit bei gleichzeitiger Erhaltung einer starken Generalisierung durch rigorose Validierungsprotokolle und automatisierte Optimierungsbeschränkungen.
Modellleistungsbewertung und Interpretation der Vorhersageergebnisse
Die Modellleistung und die Interpretation der Vorhersageergebnisse wurden umfassend mittels ROC-Kurven, F1-Werten und SHapley-additiver Erklärungsanalyse (SHAP) bewertet. Die Leistung des Modells und die Interpretation seiner Vorhersageergebnisse wurden umfassend mit ROC-Kurven, F1-Werten und SHapley Additive Explanation (SHAP)-Analyse bewertet. Anfangs wurden Vorhersagen auf dem Testdatensatz mit dem trainierten Modell generiert, und die vorhergesagten Wahrscheinlichkeiten für die positive Klasse (d. h. Klasse 1) wurden extrahiert (pred_prob). Die ROC-Kurve wurde mit dem pROC-Paket konstruiert, und die AUC des Modells sowie ihr 95%-Konfidenzintervall wurden berechnet. Der optimale Schwellenwert auf der ROC-Kurve wurde mittels Youdens J-Statistik (J = Sensitivität + Spezifität − 1) zur Bestimmung binärer Klassifikationslabels bestimmt. Basierend auf diesem Schwellenwert, der aus der ROC-Kurve abgeleitet wurde, wurden die vorhergesagten Wahrscheinlichkeiten in binäre Vorhersagelabels (0 oder 1) umgewandelt, und anschließend wurde eine Verwirrungsmatrix erzeugt. Der F1-Wert im Testsatz wurde mit der Verwirrungsmatrixfunktion berechnet, und eine Visualisierung der Verwirrungsmatrix wurde erstellt und gespeichert. Zusätzlich wurde das Modell auf einem unabhängigen externen Validierungsdatensatz aus 200 Fällen (vom First Affiliated Hospital der Wenzhou Medical University) weiter validiert. SHAP-Werte wurden mit dem "shapviz"-Paket analysiert, um die Auswirkungen jeder Variablen auf die Vorhersageergebnisse des Modells zu erläutern und so Einblicke in die Interpretation einzelner NAFLD-Likelihood-Vorhersagen zu liefern.
Statistische Methoden
"Statistische Analysen und Softwareentwicklung wurden mit R Version 4.2.3 (R Foundation for Statistical Computing, Wien, Österreich) durchgeführt. Kontinuierliche Variablen wurden zunächst mit dem Shapiro-Wilk-Test oder der visuellen Inspektion von Q-Q-Diagrammen auf Normalität bewertet. Normal verteilte Daten wurden als Mittelwert ± Standardabweichung (SD) präsentiert, und Vergleiche zwischen zwei unabhängigen Gruppen wurden mit unabhängigen Stichproben durchgeführt. Für mehrere Gruppenvergleiche wurde einseitig ANOVA zusammen mit post-hoc Tukey's HSD-Tests eingesetzt, wenn es angemessen war. Nicht normal verteilte kontinuierliche Daten wurden als Median zusammengefasst [Interquartilbereich (IQR), P25–P75], und Gruppenvergleiche wurden mit dem Mann-Whitney-U-Test für zwei unabhängige Gruppen oder dem Kruskal-Wallis-Test für mehrere Gruppen durchgeführt, gefolgt von Dunns Post-hoc-Test, falls nötig. Kategoriale Variablen wurden als Häufigkeiten und Prozentsätze (%) angegeben, und Vergleiche der Anteile zwischen den Gruppen wurden mit dem Chi-Quadrat-Test (χ-2-Test ) oder dem Fisher-Exakt-Test analysiert, wenn die erwarteten Zellzahlen unter 5 lagen. Das Signifikanzniveau wurde auf α = 0,05 (zweiseitig) festgelegt, und ein p-Wert < 0,05 wurde als statistisch signifikant angesehen.