2025-11-11T11:34:09.241880

LUME-DBN: Full Bayesian Learning of DBNs from Incomplete data in Intensive Care

Pirola, Stella, Grzegorczyk
Dynamic Bayesian networks (DBNs) are increasingly used in healthcare due to their ability to model complex temporal relationships in patient data while maintaining interpretability, an essential feature for clinical decision-making. However, existing approaches to handling missing data in longitudinal clinical datasets are largely derived from static Bayesian networks literature, failing to properly account for the temporal nature of the data. This gap limits the ability to quantify uncertainty over time, which is particularly critical in settings such as intensive care, where understanding the temporal dynamics is fundamental for model trustworthiness and applicability across diverse patient groups. Despite the potential of DBNs, a full Bayesian framework that integrates missing data handling remains underdeveloped. In this work, we propose a novel Gibbs sampling-based method for learning DBNs from incomplete data. Our method treats each missing value as an unknown parameter following a Gaussian distribution. At each iteration, the unobserved values are sampled from their full conditional distributions, allowing for principled imputation and uncertainty estimation. We evaluate our method on both simulated datasets and real-world intensive care data from critically ill patients. Compared to standard model-agnostic techniques such as MICE, our Bayesian approach demonstrates superior reconstruction accuracy and convergence properties. These results highlight the clinical relevance of incorporating full Bayesian inference in temporal models, providing more reliable imputations and offering deeper insight into model behavior. Our approach supports safer and more informed clinical decision-making, particularly in settings where missing data are frequent and potentially impactful.
academic

LUME-DBN: Vollständiges Bayessches Lernen von DBNs aus unvollständigen Daten in der Intensivmedizin

Grundinformationen

  • Papier-ID: 2511.04333
  • Titel: LUME-DBN: Full Bayesian Learning of DBNs from Incomplete data in Intensive Care
  • Autoren: Federico Pirola (Universität Milano-Bicocca), Fabio Stella (Universität Milano-Bicocca), Marco Grzegorczyk (Universität Groningen)
  • Klassifizierung: cs.LG (Maschinelles Lernen), cs.AI (Künstliche Intelligenz)
  • Veröffentlichungsdatum: 6. November 2025 (arXiv-Preprint)
  • Papierlink: https://arxiv.org/abs/2511.04333

Zusammenfassung

Dynamische Bayessche Netze (DBNs) finden zunehmend Anwendung im Gesundheitswesen, da sie komplexe zeitliche Beziehungen in Patientendaten modellieren können, während gleichzeitig Interpretierbarkeit gewährleistet wird – ein wesentliches Merkmal für klinische Entscheidungsfindung. Allerdings stammen bestehende Methoden zur Behandlung fehlender Werte in longitudinalen klinischen Datensätzen hauptsächlich aus der Literatur zu statischen Bayesschen Netzen und berücksichtigen nicht angemessen die zeitliche Natur der Daten. Diese Lücke beschränkt die Fähigkeit zur Quantifizierung zeitlicher Unsicherheit, was besonders in der Intensivmedizin kritisch ist, wo das Verständnis zeitlicher Dynamiken für die Modellzuverlässigkeit und Anwendbarkeit über verschiedene Patientengruppen hinweg entscheidend ist. Dieses Papier schlägt eine neue Gibbs-Sampling-basierte Methode zum Lernen von DBNs aus unvollständigen Daten vor, wobei jeder fehlende Wert als unbekannter Parameter behandelt wird, der einer Gaußverteilung folgt, und durch Sampling aus vollständigen bedingten Verteilungen eine prinzipiengestützte Imputation und Unsicherheitsschätzung ermöglicht wird.

Forschungshintergrund und Motivation

Kernproblem

Das Kernproblem dieser Forschung besteht darin, wie man dynamische Bayessche Netze effektiv lernen kann, wenn große Mengen fehlender Daten vorhanden sind, besonders in Anwendungen in der Intensivmedizin.

Bedeutung des Problems

  1. Klinische Dringlichkeit: In der Intensivmedizin ist eine zeitnahe und genaue Bewertung der Patientenerkrankungsentwicklung entscheidend für die Anleitung von Interventionsmaßnahmen
  2. Herausforderungen der Datenqualität: ICU-Daten sind häufig durch fehlende Werte, unregelmäßige Stichprobennahme und Messfehler beeinträchtigt
  3. Unsicherheitsquantifizierung: Traditionelle Methoden können die durch fehlende Werte eingeführte Unsicherheit nicht angemessen berücksichtigen, was zu verzerrten Parameterschätzungen führen kann

Einschränkungen bestehender Methoden

  1. Zeitliche Blindheit statischer Methoden: Bestehende Methoden zur Behandlung fehlender Daten stammen hauptsächlich aus statischen Bayesschen Netzen und berücksichtigen nicht die zeitliche Natur
  2. Unzulänglichkeiten frequentistischer Methoden: Traditionelle Imputations- oder frequentistische Methoden können die durch fehlende Werte eingeführte Unsicherheit möglicherweise nicht angemessen berücksichtigen
  3. Lokale Optimierungsprobleme: Algorithmen wie die strukturelle Erwartungsmaximierung (SEM) können in lokalen Optima steckenbleiben

Forschungsmotivation

Entwicklung eines vollständig Bayesschen Rahmens, der gleichzeitig Unsicherheit in Netzwerkstruktur, Parametern und fehlenden Werten behandeln kann, um zuverlässigere Unterstützung für klinische Entscheidungen zu bieten.

Kernbeiträge

  1. Theoretischer Beitrag: Herleitung geschlossener Lösungen für vollständig bedingte Verteilungen (FCDs) fehlender Werte in DBNs mit Nachweis ihrer Handhabbarkeit
  2. Methodische Innovation: Vorschlag des LUME-DBN-Algorithmus, der Gibbs-Sampling für die Imputation fehlender Daten mit MCMC-Strukturlernen kombiniert
  3. Experimentelle Validierung: Validierung der Methodeneffektivität an simulierten und echten ICU-Daten mit überlegener Rekonstruktionsgenauigkeit im Vergleich zu Methoden wie MICE
  4. Klinische Anwendung: Demonstration der Methode am PhysioNet-2012-Datensatz mit Aufdeckung aussagekräftiger zeitlicher Beziehungen in verschiedenen ICU-Typen

Methodische Details

Aufgabendefinition

Eingabe: Multivariate Zeitreihendaten mit fehlenden Werten D∈RN×k×(T+1)D \in \mathbb{R}^{N \times k \times (T+1)}, wobei NN die Anzahl der Stichproben, kk die Anzahl der Variablen und T+1T+1 die Anzahl der Zeitpunkte ist

Ausgabe: Posteriori-Verteilungsstichproben von DBN-Struktur, Parametern und fehlenden Werten

Einschränkungen: Annahme von Markov-Eigenschaft erster Ordnung und Abwesenheit momentaner Effekte

Modellarchitektur

DBN-Grundrahmen

Das DBN wird als kk unabhängige Bayessche lineare Regressionsmodelle (BLR) modelliert:

xit=β0(i)+∑j:(Xjt−1∈π(i))βj(i)xjt−1+ϵitx_i^t = \beta_0^{(i)} + \sum_{j:(X_j^{t-1} \in \pi(i))} \beta_j^{(i)} x_j^{t-1} + \epsilon_i^t

wobei π(i)\pi(i) die Menge der Elternknoten der Variablen XiX_i darstellt und ϵit∼N(0,σ(i)2)\epsilon_i^t \sim N(0, \sigma^2_{(i)}).

Priori-Verteilungseinstellung

  • Regressionskoeffizienten: β(i)∼N(μ(i),σ(i)2δ(i)2I)\beta^{(i)} \sim N(\mu^{(i)}, \sigma^2_{(i)}\delta^2_{(i)}I)
  • Rauschparameter: σ(i)2∼Inv-Gamma(a,b)\sigma^2_{(i)} \sim \text{Inv-Gamma}(a, b)
  • Unsicherheitsparameter: δ(i)2∼Inv-Gamma(αδ,βδ)\delta^2_{(i)} \sim \text{Inv-Gamma}(\alpha_\delta, \beta_\delta)
  • Größe der Elternknotenmenge: ∣π(i)∣∼Poisson(λ)|\pi(i)| \sim \text{Poisson}(\lambda)

Vollständig bedingte Verteilung fehlender Werte

Für den fehlenden Wert xit[MIS]x_i^t[MIS] der Variablen XiX_i zum Zeitpunkt tt ist die FCD:

P(xit[MIS]∣⋅)=N(μ∗,σ∗2)P(x_i^t[MIS] | \cdot) = N(\mu_*, \sigma^2_*)

wobei: σ∗2=(1σ(i)2+∑j:(Xit∈π(j))(βi(j))2σ(j)2)−1\sigma^2_* = \left(\frac{1}{\sigma^2_{(i)}} + \sum_{j:(X_i^t \in \pi(j))} \frac{(\beta_i^{(j)})^2}{\sigma^2_{(j)}}\right)^{-1}

μ∗=σ∗2⋅(μitσ(i)2+∑j:(Xit∈π(j))βi(j)(xjt+1−μ{−i}(j)(t+1))σ(j)2)\mu_* = \sigma^2_* \cdot \left(\frac{\mu_i^t}{\sigma^2_{(i)}} + \sum_{j:(X_i^t \in \pi(j))} \frac{\beta_i^{(j)}(x_j^{t+1} - \mu_{{\{-i\}}}^{(j)(t+1)})}{\sigma^2_{(j)}}\right)

Technische Innovationspunkte

  1. Einheitliche Imputationsstrategie: Entwurf von Gibbs-Schritten zur gemeinsamen Aktualisierung fehlender Werte über alle Regressionsmodelle hinweg
  2. Geschlossene Lösungsherleitung: Nachweis der Handhabbarkeit der FCD fehlender Werte, was effiziente MCMC-Inferenz ermöglicht
  3. Zeitliche Invarianz: Die FCD-Struktur ist zeitlich invariant bezüglich DBN-Parametern und erhöht die Recheneffizienz
  4. Flucht aus lokalen Optima: MCMC-Sampling ermöglicht das Entkommen aus lokalen Minima für genauere Netzwerkrekonstruktion

Experimentelles Setup

Datensätze

Simulierte Daten

  • Struktur: 10 unabhängige 10-Knoten-DBN-Strukturen mit maximal 5 Elternknoten pro Knoten
  • Zeitliche Länge: T∈{50,100,200}T \in \{50, 100, 200\}
  • Fehlquote: {10%,20%,30%,40%}\{10\%, 20\%, 30\%, 40\%\}
  • Parametereinstellung: Regressionskoeffizienten aus Uniform[0.2,0.8]\text{Uniform}[0.2, 0.8] gezogen, Rauschvarianz σ2=1\sigma^2 = 1

Echte Daten

  • Datenquelle: PhysioNet 2012 Challenge-Datensatz
  • Patientenzahl: 20.000+ erwachsene ICU-Patienten
  • Zeitfenster: Erste 48 Stunden des ICU-Aufenthalts
  • Variablenzahl: 11 klinische Variablen (Vitalzeichen, Blutparameter, physiologische Merkmale)
  • ICU-Gruppierung: MICU (34 Fälle), SICU (104 Fälle), CCU (114 Fälle), CSRU (62 Fälle)

Bewertungsmetriken

  • Strukturrekonstruktion: Fläche unter der Präzisions-Recall-Kurve (AUC-PR)
  • Konvergenzdiagnose: Potenzieller Skalierungsreduktionsfaktor (PSRF < 1,1)
  • Statistische Signifikanz: Gepaarter t-Test

Vergleichsmethoden

  • MICE: Multiple Imputation by Chained Equations
  • Temporal MICE: MICE-Variante mit verzögerten Prädiktorvariablen
  • Vollständige Daten: Als Leistungsobergrenze als Referenz

Implementierungsdetails

  • Sampling-Runden: 20.000 Runden mit 5.000 Burn-in-Runden
  • Aktualisierungshäufigkeit fehlender Werte: Aktualisierung alle 10 Runden (EM=10E_M = 10)
  • Kettenausdünnung: 1 von 5 Stichproben beibehalten zur Reduktion von Autokorrelation
  • Priori-Parameter: λ=1\lambda = 1, σ(i)2=δ(i)2=1\sigma^2_{(i)} = \delta^2_{(i)} = 1

Experimentelle Ergebnisse

Hauptergebnisse

Ergebnisse simulierter Daten

LUME-DBN übertrifft Baseline-Methoden signifikant in allen experimentellen Einstellungen:

  1. MICE-Leistung: Völliger Ausfall bei Fehlquoten über 20%, was seine Ineffektivität bei zeitlichen Daten widerspiegelt
  2. Temporal MICE: Bessere Leistung als MICE, aber immer noch signifikant schlechter als LUME-DBN
  3. LUME-DBN-Vorteil: Besonders hervorragende Leistung bei hohen Fehlquoten, mit geringerem Leistungsverlust im Vergleich zu vollständigen Daten bei großen Stichproben

Konvergenzanalyse

  • Strukturkonvergenz: Konvergenz innerhalb von 1,5k Runden bei allen Fehlquoten
  • Konvergenz fehlender Werte: 5k Runden erforderlich bei 40% Fehlquote
  • Konvergensstabilität: Konvergenzzeit nimmt mit Fehlquote zu, konvergiert aber letztendlich immer

ICU-Datenanalyse

Entdeckte klinische Beziehungen

  1. Selbstregulierungsschleifen: Starke interne Verbindungen zwischen Druckparametern (MAP, Sys, Dias) und atemwegsabhängigen Variablen (FiO2, PaCO2, PaO2, pH)
  2. Neurologische Wechselwirkungen: Vermindertes Bewusstsein führt zu erhöhter Herzfrequenz (CCU: GCS → HR)
  3. Hämodynamische Effekte: Blutdruck beeinflusst stark das Bewusstseinsniveau (medizinische Patienten: Dias, MAP → GCS)
  4. Körpertemperatur-Regulationsdynamik: Temperaturveränderungen während chirurgischer Genesung beeinflussen Urinausscheidung (Temp → Urine)
  5. Kardiopulmonale Rückkopplung: Niedriger Sauerstoffgehalt löst kompensatorische Herzfrequenzerhöhung aus (FiO2 → HR)

Normalisierte Auswirkungen

  • Lokale Normalisierung: Entdeckung weiterer ICU-spezifischer Beziehungen
  • Globale Normalisierung: Netzwerk zeigt mehr Gemeinsamkeiten, aber einige Beziehungen mangelt es an klinischen Belegen

Verwandte Arbeiten

Fehlende Daten in statischen Bayesschen Netzen

  • SEM-Algorithmus: Hard-EM-Variante effektiv bei begrenzten Daten, neigt aber zum Steckenbleiben in lokalen Optima
  • MCMC-Methoden: Neuere Sampling-Methoden können lokale Minima verlassen und genauere Rekonstruktionen erreichen

Dynamische Bayessche Netze

  • Bestehende Methoden: Verwenden hauptsächlich modellunabhängige Methoden wie MICE zur Behandlung fehlender Daten
  • Beitrag dieses Papiers: Erste Erweiterung von Sampling-Methoden auf die Behandlung fehlender Daten in DBNs

Klinische Anwendungen

  • Organversagenvorhersage: DBNs zur Vorhersage von Organversagensverläufen
  • Vorhersage physiologischer Veränderungen: Vorhersage physiologischer Veränderungen und Todesrisiko
  • Entscheidungsunterstützung: Bereitstellung interpretierbarer Entscheidungsunterstützung

Schlussfolgerungen und Diskussion

Hauptschlussfolgerungen

  1. Methodeneffektivität: LUME-DBN übertrifft bestehende Methoden sowohl bei der Strukturrekonstruktion als auch bei der Imputation fehlender Werte
  2. Klinische Relevanz: Entdeckte zeitliche Beziehungen haben klinische Bedeutung und unterstützen sicherere klinische Entscheidungen
  3. Unsicherheitsquantifizierung: Der vollständig Bayessche Rahmen bietet explizite Unsicherheitskodierung für Modell, Parameter und fehlende Werte

Einschränkungen

  1. Rechenkomplexität: MCMC-Sampling hat hohe Rechenkosten, erfordert Parallelisierungsoptimierung
  2. Annahme zufälliger Fehlwerte: Aktuelle Methode behandelt nur zufällige Fehlwerte; nicht-zufällige Fehlwertmuster in klinischen Daten erfordern weitere Forschung
  3. Stichprobengrößenbeschränkung: Bei kleinen Stichproben können einige Beziehungen möglicherweise nicht stabil genug sein
  4. Vorwissen: Bessere Integration klinischen Vorwissens zur Anleitung der Modellableitung erforderlich

Zukünftige Richtungen

  1. MNAR-Behandlung: Integration von Fehlwertgraph-Methoden zur Behandlung nicht-zufälliger Fehlwertmuster
  2. Nicht-homogene DBNs: Erweiterung auf global gekoppelte nicht-homogene DBNs zur Erfassung nicht-stationärer Beziehungen
  3. Gemischte Variablen: Behandlung gemischter kontinuierlicher und diskreter Variablentypen
  4. Echtzeit-Anwendungen: Entwicklung von Echtzeit-Klinischen-Entscheidungsunterstützungssystemen

Tiefgreifende Bewertung

Stärken

  1. Theoretische Strenge: Vollständige Herleitung geschlossener Lösungen für FCD fehlender Werte mit solider theoretischer Grundlage
  2. Methodische Innovation: Erste Anwendung vollständig Bayesscher Methoden auf DBN-Lernen aus fehlenden Daten
  3. Experimentelle Vollständigkeit: Validierung mit simulierten und echten Daten, Abdeckung verschiedener Fehlquoten und Stichprobengrößen
  4. Klinische Relevanz: Entdeckte Beziehungen haben klinische Bedeutung und validieren die praktische Anwendbarkeit der Methode
  5. Reproduzierbarkeit: Vollständige Algorithmusbeschreibung und Open-Source-Code bereitgestellt

Schwächen

  1. Recheneffizienz: Fehlende detaillierte Rechenzeit-Analyse und Optimierungsstrategien
  2. Frequentistische Vergleiche: Mangel an Vergleichen mit klassischen frequentistischen DBN-Lernmethoden
  3. Parameterempfindlichkeit: Unzureichende Sensitivitätsanalyse zur Hyperparameter-Auswahl
  4. Skalierbarkeit: Leistung bei größeren Netzwerken unbekannt

Auswirkungen

  1. Akademischer Beitrag: Bietet neuen theoretischen Rahmen für die Behandlung fehlender Daten in DBNs
  2. Praktischer Wert: Wichtige Anwendungsperspektiven in kritischen Bereichen wie dem Gesundheitswesen
  3. Methodische Universalität: Erweiterbar auf andere Bereiche, die Zeitreihen mit fehlenden Daten behandeln müssen

Anwendungsszenarien

  1. Gesundheitswesen: ICU-Überwachung, chronische Krankheitsmanagement, klinische Studienanalyse
  2. Finanzen: Zeitreihen-Risikomodellierung, Marktprognose
  3. Industrie: Gerätegesundheitsüberwachung, Qualitätskontrolle
  4. Umwelt: Klimamodellierung, Umweltüberwachung

Literaturverzeichnis

Das Papier zitiert 42 relevante Referenzen, die wichtige Arbeiten in mehreren Bereichen abdecken, einschließlich Bayessches Netzwerklernen, Behandlung fehlender Daten und medizinische Informatik, und bietet eine solide theoretische Grundlage für die Forschung.

Gesamtbewertung: Dies ist ein hochqualitatives Papier mit wichtigen methodischen Innovationen, das nicht nur theoretische Durchbrüche zeigt, sondern auch praktischen Wert in realen Anwendungen demonstriert. Obwohl es Raum für Verbesserungen bei Recheneffizienz und Methodenvergleichen gibt, sind seine Beiträge ausreichend, um die Entwicklung dieses Forschungsbereichs voranzutreiben.