Badanie to zostało ocenione i zatwierdzone przez Komisję Etyczną Szpitala Onkologicznego w Taizhou. Przed zebraniem danych uzyskano pisemną świadomą zgodę od wszystkich uczestników, zgodnie z instytucjonalnymi i krajowymi standardami etycznymi. Stosowane odczynniki i oprogramowanie są wymienione w Tabeli materiałów.
1. Przegląd przepływu pracy
Kompleksowy przepływ pracy przetwarzania ustrukturyzowanych danych klinicznych za pomocą XGBoost jest pokazany na rysunku 1 i obejmuje pozyskiwanie funkcji, konstrukcję drzewa, aktualizacje szczątkowe, regularyzację i ocenę. Przepływ pracy wstępnego przetwarzania i modelowania sekwencyjnego dla danych zdalnego monitorowania za pomocą LSTM pokazano na rysunku 2, obejmującym 10-minutowe ponowne próbkowanie, obsługę luk, 7-dniowe nienakładające się okna, architekturę sieci i wnioskowanie.
2. Rekrutacja pacjentów i dane demograficzne
Pacjenci byli rekrutowani kolejno z przychodni i oddziałów stacjonarnych w Taizhou Cancer Hospital w okresie od marca 2023 r. do stycznia 2024 r. Kryteriami włączenia były: potwierdzone rozpoznanie POChP zgodnie z kryteriami GOLD (FEV1/FVC < 70%), stabilny stan w momencie włączenia do badania oraz wiek od 40 do 80 lat. Kryteria wykluczenia obejmowały ciężkie choroby współistniejące, takie jak rak płuc, niekontrolowane choroby układu krążenia lub upośledzenie funkcji poznawczych uniemożliwiające świadomą zgodę. Do badania włączono łącznie 214 pacjentów (średnia wieku 63,7 ± 8,9 lat; 68% mężczyzn; 54% obecnych lub byłych palaczy).
3. Gromadzenie danych
Ustrukturyzowane dane kliniczne obejmowały czynność płuc (FEV1, FVC), wysycenie krwi tlenem (SpO2), długość pobytu w szpitalu i czas trwania choroby. Czynność płuc i SpO2 były rejestrowane w procentach, pobyt w szpitalu w dniach, a czas trwania choroby w latach. Dane ze zdalnego monitorowania zostały zebrane za pomocą certyfikowanych pulsoksymetrów noszonych na ciele i monitorów aktywności zweryfikowanych zgodnie ze szpitalnymi złotymi standardami: pulsoksymetry zostały sprawdzone krzyżowo z monitorami Masimo Rad-97 (średni błąd bezwzględny 1,8% w zakresie 90-100% SpO2 ), czujniki tętna zostały zweryfikowane pod kątem elektrokardiografii (średni błąd 2,3 bpm), a liczniki kroków zostały skalibrowane na protokole bieżni. Nieprzetworzone dane zostały wyeksportowane jako pliki CSV ze znacznikiem czasu w 10-minutowym interwale próbkowania.
4. Uzasadnienie wielkości próby
Progi 300 próbek treningowych dla XGBoost i około 1000 nienakładających się 7-dniowych sekwencji dla LSTM zostały poparte eksperymentami i analizami mocy/symulacji. Analiza mocy post hoc przy użyciu G*Power (wersja 3.1) wykazała, że 300 próbek dostarczyło >80% mocy do wykrycia średniej wielkości efektu (f2 Cohena = 0,15) przy α = 0,05 w regresji logistycznej. Symulacje wykazały, że około 1000 sekwencji było niezbędnych do stabilnej zbieżności LSTM w wielowymiarowych fizjologicznych szeregach czasowych. Dowody empiryczne podsumowano w tabeli 1 i zwizualizowano na rysunkach 3 i 4.
5. Wstępne przetwarzanie danych
Brakujące wpisy, zidentyfikowane jako puste miejsca, wartości wskaźnikowe lub NaN, zostały przypisane przy użyciu średniej zestawu treningowego dla każdej funkcji, aby uniknąć wycieku celu:
(1)
gdzie mj jest liczbą obserwowanych wartości dla cechy j. To samo
dotyczyło zestawów walidacyjnych i testowych. Aby usunąć odchylenie związane ze skalą, funkcje zostały ustandaryzowane z normalizacją wskaźnika z przy użyciu parametrów tylko do trenowania:
(2)
gdzie μj i σj są średnią i odchyleniem standardowym cechy j oszacowanymi na podstawie danych treningowych. W przypadku sygnałów szeregów czasowych (SpO2, tętno, kroki) strumienie zostały wyrównane do 10-minutowej kadencji; Krótkie luki do 30 minut były wypełniane do przodu, a dłuższe luki były wygładzane trzypunktową średnią kroczącą. Następnie zastosowano 7-dniowe okno przesuwne z 1008 krokami czasowymi, aby utworzyć nienakładające się sekwencje, aby zapobiec wyciekom (patrz rysunek 2).
6. Trening modelowy
W przypadku ustrukturyzowanych danych klinicznych XGBoost został dostrojony za pomocą wyszukiwania w siatce pod kątem szybkości uczenia się (0,01-0,3), maksymalnej głębokości (3-10) i liczby estymatorów (100-500) w ramach binarnego celu logistycznego (przepływ pracy na rysunku 1). W przypadku danych zdalnego monitorowania LSTM składał się z dwóch ukrytych warstw LSTM po 128 jednostek każda, inicjalizacji Xavier, współczynnika zaniku 0,5 i warstwy wyjściowej sigmoidalnej; optymalizacja używana przez Adama ze współczynnikiem uczenia się 0,001 zaimplementowanym w TensorFlow (potok na rysunku 2). Nadmierne dopasowanie zostało złagodzone poprzez rezygnację i wczesne zatrzymanie (cierpliwość = 10), a nierównowaga klas została naprawiona za pomocą SMOTE.
7. Strategia oceny
Dane ustrukturyzowane zostały ocenione za pomocą warstwowej 10-krotnej walidacji krzyżowej. Dane szeregów czasowych zostały ocenione z walidacją przejścia do przodu w celu zachowania porządku czasowego. Metryki obejmowały dokładność, precyzję, kompletność, miarę F, pole pod krzywą ROC (AUC) i błąd średniokwadratowy (MSE). Różnice między modelami oceniano za pomocą testu rang ze znakiem Wilcoxona (dwustronnego). Dane liczbowe obejmują 95% przedziały ufności lub słupki błędów w celu ilościowego określenia niepewności. Wybory te bezpośrednio odnoszą się do ryzyka wycieku szeregów czasowych i potrzeby istotności statystycznej wymaganej przez recenzentów.
8. Znaczenie kliniczne i badania pilotażowe
Model hybrydowy przewidywał spadek SpO2 na 6-12 godzin przed wystąpieniem objawów klinicznych u 78% monitorowanych pacjentów i zmniejszył MSE o 42,5% w stosunku do wartości wyjściowych. W czterotygodniowym programie pilotażowym z udziałem 20 pacjentów z POChP, tygodniowe modelowe oceny ryzyka zgadzały się z ocenami lekarzy w 85% przeglądów, co potwierdza potencjał integracji klinicznej.
9. Oprogramowanie i środowisko
Analizy przeprowadzono w Pythonie 3.10.6 z scikit-learn 1.2.2, XGBoost 1.7.5, TensorFlow 2.13 i PyTorch 1.13 na Ubuntu 20.04 LTS z procesorem graficznym NVIDIA RTX 3080, CUDA 11.6 i cuDNN 8.2. Pełne wersje i dostawcy są wymienieni w nienumerowanym Tabeli materiałów.