Datenquelle
Diese Studie war eine retrospektive Analyse basierend auf der Medical Information Mart for Intensive Care IV (MIMIC-IV, Version 3.1) Datenbank, einer öffentlich zugänglichen Intensivpflegedatenbank, die auf PhysioNet (https://physionet.org/content/mimiciv/3.1/) gehostet wird. Die Datenbank wurde in Zusammenarbeit zwischen dem Massachusetts Institute of Technology Laboratory for Computational Physiology, dem Beth Israel Deaconess Medical Center und Philips entwickelt und enthält umfassende de-identifizierte klinische Daten, darunter demografische Informationen, Vitalzeichen, Labormessungen und Medikamentenakten.
Alle Patientendaten in MIMIC-IV sind gemäß dem Health Insurance Portability and Accountability Act (HIPAA) vollständig de-identifiziert, weshalb diese Studie keine Genehmigung des institutionellen Überprüfungsausschusses (IRB) benötigte. Die Nutzung der Datenbank wird durch eine Datennutzungsvereinbarung (DUA) geregelt, und der Zugang wird nur autorisierten Nutzern gewährt, die die erforderliche Schulung und Zertifizierung abgeschlossen haben. Die Autoren absolvierten die erforderliche Schulung, erhielten Zugangsdaten (PhysioNet-akkreditierter Nutzer) und führten diese Studie in Übereinstimmung mit allen relevanten Datennutzungsvorschriften durch.
Untersuchung von Population und Datenextraktion
Erwachsene Patienten (≥18 Jahre), die sich einem spinalen neurochirurgischen Eingriff unterzogen haben, wurden anhand der MIMIC-IV-Datenbank mittels verfahrensbezogener ICD-Codes identifiziert. Spinale neurochirurgische Eingriffe wurden definiert als ICD-9-Codes, beginnend mit 03 (Operationen am Rückenmark und am Rückenmarkskanal) und ICD-10-Codes, beginnend mit 00B, 00H, 00J, 00N, 00P, 00Q, 00R, 00S, 00T oder 00U (rückenmarksbezogene Verfahren). Bei Patienten mit mehreren Krankenhausaufenthalten wurde nur die erste Aufnahme mit einem qualifizierenden Verfahren beibehalten, um wiederholte Messungen derselben Person zu vermeiden.
Die Ergebnisvariable, Osteoporose, wurde anhand von ICD-9-Codes definiert, die mit 733 beginnen, und ICD-10-Codes, die mit M80, M81 oder M82 beginnen, die während der Indexaufnahme registriert wurden. Da sowohl Prädiktoren als auch Ergebnis aus demselben Krankenhausaufenthalt abgeleitet wurden, wurde die analytische Aufgabe in dieser Studie als Klassifikation innerhalb des Krankenhauses und nicht als prospektive Risikovorhersage dargestellt.
Insgesamt erfüllten 10.803 Patienten die Einschlusskriterien. Nach der stratifizierten Aufteilung wurde der Datensatz in einen Trainingssatz (n = 7.561), einen Validierungssatz (n = 1.621) und einen Testsatz (n = 1.621) unterteilt. Die Klassenverteilung war über alle Partitionen hinweg stark unausgewogen. Im Trainingsset waren 499 Patienten (6,60 %) osteoporose-positiv und 7.062 (93,40 %) osteoporose-negativ. Im Validierungssatz waren 107 Patienten (6,60 %) positiv und 1.514 (93,40 %) negativ. Der Testsatz zeigte die gleiche Verteilung, mit 107 positiven (6,60 %) und 1.514 negativen (93,40 %) Fällen.
Klinische Daten wurden mit relationaler Datenbankverwaltungssoftware von MIMIC-IV extrahiert. Extrahierte Variablen umfassten demografische Merkmale (z. B. Alter und Geschlecht), Labormessungen, Vitalzeichen und Medikamentenakten. Für jeden Patienten wurden Daten aus der Indexaufnahme für die anschließende Merkmalkonstruktion verwendet.
Um eine valide Modellbewertung sicherzustellen, wurde der Datensatz zunächst in Trainings-, Validierungs- und Testsätze mittels stratifizierter Stichprobe (70 %, 15 % und 15 %) aufgeteilt. Alle nachfolgenden Vorverarbeitungsschritte, die Verzerrungen einführen konnten, einschließlich Merkmalauswahl und Resampling, wurden ausschließlich innerhalb des Trainingssets durchgeführt. Das Ungleichgewicht der Klassen im Trainingssatz wurde durch eine Kombination aus Downsampling der Mehrheitsklasse und der Synthetic Minority Over-sampling Technique (SMOTE) behoben, während die Validierungs- und Testsätze die ursprüngliche Klassenverteilung beibehalten haben, um die reale Prävalenz widerzuspiegeln.
Feature-Engineering
Eine mehrstufige Feature-Engineering-Pipeline wurde implementiert, um die Modellleistung zu verbessern und gleichzeitig die Interpretierbarkeit zu erhalten. Für jeden Patienten wurden Medikamente-, Labor- und Vitalzeichenaufzeichnungen aus der Indexaufnahme zu Zusammenfassungsmerkmalen aggregiert (z. B. Durchschnittswerte für wiederholte Messungen und durchschnittliche Exposition für Medikamente).
Um die Funktionsknappheit zu verringern, wurden Kandidaten-Merkmalsgruppen mit mehr als 30 % fehlenden Werten ausgeschlossen. Um Datenlecks zu vermeiden, wurde dieser Filterschritt ausschließlich mit dem Trainingssatz durchgeführt, und derselbe Funktionssatz wurde dann auf die Validierungs- und Testsätze angewendet.
Anschließend wurde auf Feature-Gruppenebene ein Univariat-Merkmals-Screening mit dem Mann– Whitney-U-Test durchgeführt, um die Verteilungen zwischen osteoporose-positiven und osteoporose-negativen Patienten zu vergleichen. Feature-Gruppen wurden nach statistischer Signifikans bewertet, und die Top-k Gruppen (k = 20) wurden für die nachgelagerte Modellierung beibehalten. Dieses Auswahlverfahren wurde ebenfalls ausschließlich im Trainingsset durchgeführt, um Informationslecks zu verhindern.
Angesichts des relativ hohen Mangels zwischen den Kandidaten-Merkmalsgruppen bestand die Hauptaufgabe dieses Top-K-Auswahlschritts darin, Merkmale mit schwachen statistischen Signalen herauszufiltern, anstatt die Merkmalsdimensionalität strikt zu optimieren. In der Praxis war die Modellleistung nicht besonders empfindlich gegenüber dem exakten Wert von k innerhalb eines angemessenen Bereichs, was darauf hindeutet, dass dieser Schritt hauptsächlich als robustheitsorientiertes Screeningverfahren diente, um wenig informative Merkmale auszuschließen und klinisch relevante Signale zu erhalten.
Grundlegende demografische Merkmale wurden in allen Modellen beibehalten. Fehlende Werte in diesen demografischen Variablen wurden anhand des aus dem Trainingssatz berechneten Mittelwerts berechnet. Die mittlere Imputation wurde wegen ihrer Einfachheit und Stabilität in diesem großen Datensatz gewählt, obwohl sie die Varianz verringern und Verzerrungen einführen kann; Diese Einschränkung wird anerkannt.
Für ausgewählte Medikamente-, Labor- und Vitalzeichenmerkmale wurde das Fehlen einer aufgezeichneten Messung oder Exposition als Null kodiert, um eine konsistente numerische Eingabe für das Modelltraining zu ermöglichen. Wir erkennen an, dass dieser Ansatz wahre Nullwerte mit Fehlverhalten oder Nicht-Exposition vermischen kann, insbesondere bei Variablen, bei denen Null physiologisch nicht sinnvoll ist. Diese Kodierungsstrategie wurde jedoch in allen Stichproben konsistent angewandt und von der Z-Score-Standardisierung gefolgt, die skalenbedingte Verzerrungen mindert. Die potenzielle Auswirkung dieser Annahme wird als Einschränkung diskutiert.
Alle numerischen Merkmale wurden mittels Z-Score-Normalisierung standardisiert, basierend auf dem aus dem Trainingssatz geschätzten Mittelwert und Standardabweichung:
z = (x — μ)/σ
Dabei steht x für den rohen Merkmalswert, und μ und σ den Mittelwert und die Standardabweichung aus dem Trainingssatz bezeichnen. Kategoriale Variablen wurden numerisch codiert, um die Eingabeanforderungen der Deep-Learning-Modelle zu erfüllen.
Diese Feature-Engineering-Strategie ermöglichte eine systematische Reduzierung der Dimensionalität bei gleichzeitiger Erhaltung klinisch relevanter Bereiche, darunter Demografie, Nierenfunktion, hämatologische Indzes, Entzündungsmarker und metabolische Parameter. Detaillierte Ausgangsmerkmale der Studienpopulation über die Trainings-, Validierungs- und Testsets hinweg sind in Tabelle 1 dargestellt.
| Merkmale | Gesamt | Trainingsset | Validierungsmenge | Testsatz | P-Wert |
| (n=5238) | (n=1996) | (n=1621) | (n=1621) |
| Geschlecht | | <0,001 |
| Männlich | 2477 (47.29%) | 834 (41.78%) | 797 (49.17%) | 846 (52.19%) | |
| Weiblich | 2761 (52.71%) | 1162 (58.22%) | 824 (50.83%) | 775 (47.81%) | |
| Alter | 60.69 (48.29, 70.61) | 63.45 (52.00, 72.41) | 58.00 (46.00, 69.00) | 60.00 (46.00, 70.00) | <0,001 |
| Höhe | 165.23 (160.79, 171.21) | 164.25(161.41, 172.23) | 166.32 (161.85, 171.22) | 165.22 (161.31, 169.89) | <0,001 |
| Gewicht | 72.01 (56.77, 82.46) | 75.45 (63.22, 81.45) | 70.55 (60.22, 79.33) | 73.62 (65.22, 81.88) | <0,001 |
| Senna | 1.00 (-0.73, 8.60) | 1.00 (-0.30, 8.60) | 1.00 (-1.00, 8.60) | 1.00 (-1.00, 8.60) | <0,001 |
| Docusatnatrium | 100.00 (100.00, 100.00) | 100.00 (100.00, 100.00) | 100.00 (100.00, 100.00) | 100.00 (100.00, 100.00) | 0.001 |
| Heparin | 5000.00 (1655.28, 5000.00) | 5000.00 (-1.00, 5000.00) | 5000.00 (1600.00, 5000.00) | 5000.00 (3750.00, 5000.00) | 0.292 |
| Natriumchlorid 0,9 % Spülung | 3.00 (1.76, 3.00) | 3.00 (3.00, 3.00) | 3.00 (-1.00, 3.00) | 3.00 (3.00, 3.00) | <0,001 |
| Rote Blutkörperchen | 3.96 (3.51, 4.35) | 3.90 (3.49, 4.28) | 4.02 (3.53, 4.40) | 3.97 (3.50, 4.38) | <0,001 |
| Hämoglobin | 11.81 (10.51, 13.02) | 11.66 (10.47, 12.82) | 11.91 (10.58, 13.14) | 11.90 (10.49, 13.15) | <0,001 |
| RDW | 13.98 (13.22, 15.15) | 14.12 (13.35, 15.24) | 13.90 (13.10, 15.06) | 13.90 (13.18, 15.13) | <0,001 |
| Hämatokrit | 35.67 (31.95, 39.00) | 35.32 (31.63, 38.50) | 36.00 (32.28, 39.37) | 35.77 (32.01, 39.25) | <0,001 |
| Spezifisches Gewicht | 1.01 (1.01, 1.02) | 1.01 (1.01, 1.02) | 1.01 (1.01, 1.02) | 1.01 (1.01, 1.02) | <0,001 |
| Bicarbonat | 25.72 (24.04, 27.15) | 25.86 (24.17, 27.27) | 25.65 (23.92, 27.03) | 25.62 (24.00, 27.12) | <0,001 |
| Neutrophile | 68.67 (60.52, 76.23) | 69.41 (61.90, 76.52) | 68.10 (59.60, 75.88) | 68.34 (59.75, 76.23) | <0,001 |
| Kreatinin | 0.82 (0.68, 1.02) | 0.81 (0.68, 1.00) | 0.83 (0.69, 1.02) | 0.83 (0.68, 1.03) | <0,001 |
| Phosphat | 3.29 (2.92, 3.64) | 3.27 (2.93, 3.57) | 3.30 (2.93, 3.69) | 3.30 (2.90, 3.65) | 0.002 |
| Harnstoffstickstoff | 15.43 (12.08, 20.18) | 15.74 (12.50, 20.20) | 15.00 (11.64, 19.89) | 15.48 (12.00, 20.43) | <0,001 |
| Lymphozyten | 19.12 (12.22, 26.15) | 18.81 (12.49, 25.30) | 19.36 (12.40, 27.07) | 19.24 (11.72, 26.27) | 0.001 |
| Alanin-Aminotransferase (ALT) | 20.49 (10.84, 32.82) | 20.47 (11.67, 32.06) | 20.00 (10.00, 32.93) | 21.00 (10.67, 33.67) | 0.3283 |
| MCV | 90.90 (87.50, 94.32) | 91.03 (87.64, 94.51) | 90.64 (87.43, 94.00) | 91.00 (87.40, 94.40) | <0,001 |
| MCHC | 33.14 (32.37, 33.92) | 33.11 (32.37, 33.88) | 33.16 (32.38, 33.97) | 33.16 (32.35, 33.93) | <0,001 |
| Basophile | 0.42 (0.27, 0.61) | 0.41 (0.27, 0.59) | 0.45 (0.29, 0.63) | 0.41 (0.26, 0.60) | 0.037 |
Tabelle 1. Basischarakteristika der Studienpopulation, geschichtet nach Training, Validierung und Testsätzen.
Modellarchitektur
Um komplexe Beziehungen zwischen heterogenen tabellaren klinischen Merkmalen zu modellieren, verwendeten wir eine auf TabTransformer basierende Architektur, die mit dem PyTorch-Framework6 implementiert wurde. Das Modell wurde entwickelt, um sowohl kategorische als auch numerische Variablen zu verarbeiten und kontextuelle Interaktionen zwischen Merkmalen durch einen Selbst-Aufmerksamkeitsmechanismus zu erfassen.
Eingabedarstellung
Kategoriale Variablen wurden zunächst in dichte Einbettungen umgewandelt. Jedes kategoriale Merkmal wurde auf einen Einbettungsvektor der Dimension d(embedding_dim = 256) abgebildet. Numerische Merkmale wurden mittels Z-Score-Normalisierung standardisiert und dann durch eine lineare Transformationsschicht in denselben Einbettungsraum projiziert, was eine gemeinsame Verarbeitung mit kategorialen Merkmalen ermöglichte.
Die eingebetteten Merkmalsvektoren wurden zu einer Folge von Tokens verkettet:
X = [x1,x 2,... xn]
wobei jedes xi ∈ Rd ein eingebettetes Merkmal darstellt.
Transformatorencoder
Die verketteten Feature-Einbettungen wurden durch einen Stapel von Transformer-Encoder-Schichten geleitet (num_layers = 3). Jede Schicht bestand aus Mehrkopf-Selbstaufmerksamkeit, gefolgt von einem positionsweise Feedforward-Netzwerk.
Der Multi-Head-Selbstaufmerksamkeitsmechanismus ist definiert als:
—> Achtung(Q,K,V) = softmax
V
Dabei stehen Q, K und V für Abfrage-, Schlüssel- und Wertmatrizen. Mehrere Aufmerksamkeitsköpfe (num_heads = 8) wurden verwendet, um verschiedene Merkmalsinteraktionen zu erfassen.
Jeder Transformatorblock umfasste Multi-Head Self-Attention, Residualverbindung und Schichtnormalisierung, Feed-Forward-Netzwerk, Dropout (Rate = 0,243).
Merkmalaggregation und Klassifikation
Die Ausgabe der letzten Transformer-Schicht wurde abgeflacht und durch ein mehrschichtiges Perzeptron (MLP) zur Klassifikation geleitet. Das MLP bestand aus einer oder mehreren vollständig zusammenhängenden Schichten mit nichtlinearen Aktivierungsfunktionen.
Die finale Ausgangsschicht nutzte eine sigmoidale Aktivierungsfunktion, um die vorhergesagte Osteoporosewahrscheinlichkeit zu erzeugen:
—>
= σ(z)
wobei z die Logit-Ausgabe ist.
Verlustfunktion und Trainingsstrategie
Das Modell wurde mit binärem Kreuzentropieverlust trainiert:
—> L = −(1/N) Σ [ yi log(ŷi) + (1 − yi) log(1 − ŷi) ]
Wobei yi das wahre Label und
i die vorhergesagte Wahrscheinlichkeit bezeichnet.
Das Modell wurde mit dem Adam-Optimierer (Lernrate = 1,9e-4, Batchgröße = 64) für 100 Epochen optimiert. Frühzeitiges Stoppen wurde basierend auf der Validierungsleistung angewendet (Geduld = 15), um Überanpassungen zu verhindern. Zur Unterstützung der Reproduzierbarkeit wurden alle experimentellen Verfahren mit einem konstanten zufälligen Seed (Seed = 42) durchgeführt.
Modellinterpretierbarkeit
Um die Interpretierbarkeit zu verbessern, wurden SHapley Additive Explanations (SHAP) auf das trainierte Modell angewendet. SHAP-Werte wurden auf dem Testsatz berechnet, um den Beitrag jedes Merkmals zum vorhergesagten Ergebnis zu quantifizieren, was sowohl eine globale als auch individuelle Interpretation des Modellverhaltens ermöglichte.
Wie in Abbildung 1 dargestellt, fasst der Workflow den Prozess von der MIMIC-IV-Datenvorverarbeitung bis hin zur Modellkonstruktion und -klassifikation zusammen. Die TabTransformer-Architektur erfasst Interaktionen zwischen heterogenen tabellaren Merkmalen durch Selbstaufmerksamkeit, während SHAP-basierte Visualisierung interpretierbare Einblicke in Feature-Beiträge liefert. Dieser Rahmen ermöglicht eine intern validierte Analyse osteoporosebedingter Muster in einer heterogenen spinalen neurochirurgischen Kohorte.

Abbildung 1. Studiere Workflow und Modellarchitektur. Diese Abbildung veranschaulicht das gesamte Studiendesign, einschließlich Kohortenauswahl, Datenvorverarbeitung, Feature Engineering, Datensatzaufteilung und Modellentwicklung. Es fasst außerdem die für die Klassifikation verwendete TabTransformer-Architektur zusammen, einschließlich Merkmal-Einbettung, Transformer-Encoder-Schichten und der endgültigen Klassifikationsausgabe. Nur die erste qualifizierende Aufnahme pro Patient war enthalten. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.