Źródło danych i projekt badania
Przeprowadzono retrospektywne badanie kohortowe na podstawie danych pochodzących z wersji 2.2 bazy danych Medical Information Mart for Intensive Care IV. Ta baza danych obejmuje dwa wewnętrzne systemy: kompleksową szpitalną elektroniczną dokumentację medyczną (EHR) oraz specyficzny dla OIOM-u system informacji klinicznej, obejmujący dane z lat 2008–2024 18. Badanie to przeprowadzono zgodnie z wytycznymi instytucjonalnymi. Baza danych MIMIC-IV uzyskała zgodę Massachusetts Institute of Technology oraz Beth Israel Deaconess Medical Center. Badanie to zostało przeanalizowane i zatwierdzone przez Komitet Etyki Badań Ludzkich Szpitala nr 2 w Ningbo (numer zatwierdzenia: YJ-NBEY-KY-2022-107-01, dzień 7 września 2022). Ze względu na retrospektywny projekt i anonimizowany charakter informacji zdrowotnych pacjenta, ta sama komisja etyczna zrezygnowała z wymogu świadomej zgody na podstawie tej samej zgody (YJ-NBEY-KY-2022-107-01). Dostęp do bazy danych zapewnił jeden z autorów (JQ L), który uzyskał niezbędne uwierzytelnienie i ukończył egzamin Collaborative Institutional Training Initiative (numer uwierzytelniający 60691748).
Ze względu na retrospektywny charakter tego badania, tradycyjne obliczenia wielkości próby ani analizy mocy nie zostały przeprowadzone. Jednak ilość dostępnych rekordów w bazie została uznana za wystarczającą, by spełnić cele badania. Ocena ta opiera się na zbadaniu związku między albumin-korygowaną luką anionową (ACAG) a wynikami śmiertelności w dużej i zróżnicowanej populacji pacjentów. Próbki wygodne pobrane z bazy danych skutecznie reprezentują populację intensywnej terapii, co wspiera analizy statystyczne stosowane w tym badaniu.
Uczestnicy
Badanie objęło wszystkich pacjentów z sepsą z bazy danych MIMIC-IV v2.2. Sepsa została zdefiniowana według kryteriów Sepsy-3.0, implementowanych za pomocą zapytań języka zapytań strukturalnych (SQL) do identyfikacji pacjentów z podejrzeniem infekcji oraz sekwencyjnej oceny niewydolności narządów (SOFA) zmiany wyniku ≥2 w ciągu 24 godzin19. Wszyscy pacjenci z tej grupy byli przyjmowani na OIT jako pierwszy punkt wejścia do szpitala; pacjenci przyjmowani na oddziały ogólne, a później przenoszeni na oddział intensywnej terapii, byli wykluczani, aby zapewnić jednorodną populację z wczesnym przyjęciem na OIOM. Ekstrakcja danych została przeprowadzona za pomocą SQL za pomocą narzędzia do zarządzania bazami danych20, a wszystkie wykorzystane oprogramowanie, bazy danych i narzędzia analityczne zostały podane w alfabetycznie uporządkowanej tabeli materiałów. Kryteria włączenia określały pacjentów z sepsą powyżej 18 roku życia, którzy po raz pierwszy trafili na OIOM. Kryteria wykluczenia obejmowały: (1) pacjentów poniżej 18 roku życia; (2) pobyty na OIOM-ie krótsze niż 48 godzin; (3) wielokrotne hospitalizacje na OIOM z powodu sepsy; oraz (4) niewystarczające dane, takie jak brakujące zapisy dotyczące albuminy lub luki anionowej (Rysunek 1).
Procedury badawcze i definicje
Ekstrakcja danych z MIMIC-IV była przeprowadzana za pomocą zapytań SQL. Wyodrębnione dane obejmowały demografię pacjentów (wiek, wzrost, wagę, płeć, status ubezpieczeniowy, rasę i stan cywilny), historię medyczną (w tym nadciśnienie, cukrzycę typu 2, niewydolność serca, zawał mięśnia sercowego, nowotwory złośliwe, przewlekłe choroby nerek, marskość wątroby, zapalenie wątroby, zapalenie płuc, przewlekłą obturacyjną chorobę płuc oraz hiperlipidemię) oraz wstępne wyniki badań laboratoryjnych. Wyniki te obejmowały liczbę białych i czerwonych krwinek, sód, potas, wapń, chlorki, glukozę, lukę anionową, pH krwi, czas trombinowy, fibrynogen, czas częściowy tromboplastyny, międzynarodowy stosunek normalizacji, cholesterol całkowity, aminotransferazę asparagianową oraz aminotransferazę alaninową. Dodatkowo, informacje o specjalnych metodach leczenia (takich jak ciągła terapia nerkowa i wentylacja mechaniczna), wynikach klinicznych (SOFA, ostra fizjologia i ocena zdrowia przewlekłego III [APACHE III], uproszczony wynik fizjologii ostrej II [SAPS II], wskaźnik ostrego nasilenia choroby Oxford [Oasis], wskaźnik współistniejących Charlsona [Charlson] oraz wyniki skali śpiączki Glasgow [GCS]) oraz wyniki kliniczne (długość pobytu w szpitalu, zebrano śmiertelność w szpitalach, pobyt na OIOM-ie oraz śmiertelność na OIOM-ie). Podczas czyszczenia danych wykluczono predyktory z brakującymi >30%.
Luka anionowa skorygowana albuminą (ACAG) została obliczona na podstawie pierwszych dostępnych pomiarów szczeliny anionowej w surowicy i albuminy w surowicy uzyskanej w ciągu 24 godzin od przyjęcia na OIOM. Jeśli w ciągu pierwszych 24 godzin dostępnych było wiele pomiarów, używano wartości najbliższych terminowi przyjęcia. Wzór na ACAG to:
Albumina - skorygowana luka anionowa = Luka anionowa + [40 - albumina serumowa (g/dL)] × 0,2521 (1)
Wyniki i miary
Główne wyniki tego badania to 30-dniowa śmiertelność hospitalizowana oraz 30-dniowa śmiertelność na OIOM-ie. Data indeksu analizy przeżycia to czas przyjęcia na OIOM. Pacjenci byli obserwowani przez 30 dni od tej daty indeksowej, niezależnie od statusu wypisu ze szpitala. Zdarzenia śmiertelne po wypisie były rejestrowane za pomocą rejestrów indeksu zgonów zabezpieczeń społecznych (SSDI) zintegrowanych z bazą danych MIMIC-IV, co zapewniało dokładne określenie zgonów poza szpitalem.
Analiza statystyczna
Zmienne ciągłe wyrażano jako średnią (odchylenie standardowe) lub medianę (zakres interkwartylowy), natomiast zmienne kategoryczne jako częstość i procent. Dla danych zgodnych z rozkładem normalnym stosowano testy t lub analizę wariancji (ANOVA). Dane niezgodne z rozkładem normalnym zostały przeanalizowane za pomocą testu Manna-Whitneya lub testu Kruskal-Wallis. Analiza przeżycia Kaplana-Meiera została wykorzystana do oceny częstości występowania zdarzeń końcowych na różnych poziomach ACAG, przy czym różnice oceniano za pomocą testu log-rank (Rysunek 2).
Model ryzyk proporcjonalnych Coxa został użyty do obliczenia wskaźników zagrożeń (HR) oraz 95% przedziałów ufności (CI) dla związku między ACAG a śmiertelnością. ACAG analizowano zarówno jako zmienną ciągłą, jak i według kwartyli. Dostosowania modelu zostały dokonane w następujący sposób: model 1 obejmował analizę jednowymiarową; Model 2 skorygowany o wiek, płeć, wzrost, wagę i rasę; oraz model 3 dodatkowo skorygowany o status ubezpieczeniowy, stan cywilny, liczbę białych krwinek i czerwonych krwinek, RDW, albuminę, chlorki, ALT, AST, SAF, wynik APACHE III, wynik SAPS II, wynik Oasis, wynik Charlsona oraz szereg chorób współistniejących, w tym nadciśnienie, cukrzycę typu 2, niewydolność serca, zawał mięśnia sercowego, nowotwory złośliwe, przewlekłe choroby nerek, ostrą niewydolność nerek, marskość wątroby, zapalenie wątroby, zapalenie płuc, udar, hiperlipidemia, ostre uszkodzenie nerek oraz przewlekła obturacyjna choroba płuc. Kolilinearność wśród kowarianów w modelu 3 oceniono za pomocą współczynnika wariancji inflacji (VIF), przy czym VIF > 5 wskazywał na istotną kolinearyczność. Nie wykryto istotnej kolineariczności (wszystkie VIF < 3).
Zastosowano model regresji ograniczonych sześciennych (RCS) z czterema węzłami, aby zbadać nieliniową zależność między bazowym ACAG a śmiertelnością zarówno w szpitalach, jak i na OIOM-ie. Analiza RCS konkretnie sprawdzała punkt zwrotny, w którym zmienia się nachylenie krzywej ryzyka śmiertelności. Analizy podgrupowe przeprowadzono w celu zbadania potencjalnych różnic między warstwami na podstawie wieku (≤70 lat vs. >70 lat), płci, indeksu masy ciała (BMI; <21,4, 21,4–31,2, ≥31,2 kg/m2), obecności nadciśnienia, cukrzycy typu 2, zawału mięśnia sercowego, przewlekłej choroby nerek, udaru, ostrego uszkodzenia nerek, a także stosowania ciągłej terapii nerkowymiany (CRRT) i wentylacji mechanicznej. Analizy te oceniały spójność prognostycznej wartości ACAG dla wyników głównych. Modele Coxa były również wykorzystywane w analizach podgrup do korygowania wszystkich zmiennych bazowych. Wartości p interakcji zostały obliczone w celu oceny modyfikacji efektu. Przetwarzanie i analizę danych przeprowadzono w wersji R 4.3.0 (RRID: SCR_001905), z istotnością statystyczną ustawioną na dwustronne p < 0,05. Do obsługi brakujących wartości zastosowano wielokrotną imputację za pomocą równań łańcuchowych (MICE) przy użyciu pakietu R "myszy" przy założeniu braku losowego. Przed imputacją wykluczono zmienne, których brakujące wartości przekraczały 30%.