Questo studio è stato esaminato e approvato dal Comitato Etico del Taizhou Cancer Hospital. Il consenso informato scritto è stato ottenuto da tutti i partecipanti prima della raccolta dei dati, in conformità con gli standard etici istituzionali e nazionali. I reagenti e il software utilizzati sono elencati nella Tabella dei Materiali.
1. Panoramica del flusso di lavoro
Il flusso di lavoro end-to-end per l'elaborazione di dati clinici strutturati con XGBoost è illustrato nella Figura 1 e include l'acquisizione di funzionalità, la costruzione di alberi, gli aggiornamenti residui, la regolarizzazione e la valutazione. Il flusso di lavoro di pre-elaborazione e modellazione delle sequenze per i dati di monitoraggio remoto con LSTM è illustrato nella Figura 2, che comprende il ricampionamento di 10 minuti, la gestione dei gap, la finestra non sovrapposta di 7 giorni, l'architettura di rete e l'inferenza.
2. Reclutamento dei pazienti e dati demografici
I pazienti sono stati reclutati consecutivamente da cliniche ambulatoriali e reparti di degenza presso il Taizhou Cancer Hospital tra marzo 2023 e gennaio 2024. I criteri di inclusione erano una diagnosi confermata di BPCO secondo i criteri GOLD (FEV1/FVC < 70%), una condizione stabile all'arruolamento e un'età compresa tra 40 e 80 anni. I criteri di esclusione includevano gravi comorbidità come il cancro ai polmoni, malattie cardiovascolari incontrollate o deterioramento cognitivo che impediva il consenso informato. Sono stati arruolati un totale di 214 pazienti (età media 63,7 ± 8,9 anni; 68% maschi; 54% fumatori o ex fumatori).
3. Raccolta dei dati
I dati clinici strutturati includevano la funzione polmonare (FEV1, FVC), la saturazione di ossigeno nel sangue (SpO2), la durata della degenza ospedaliera e la durata della malattia. La funzione polmonare e la SpO2 sono state registrate in percentuale, la degenza ospedaliera in giorni e la durata della malattia in anni. I dati del monitoraggio remoto sono stati raccolti utilizzando ossimetri indossabili certificati e tracker di attività convalidati rispetto agli standard ospedalieri: i pulsossimetri sono stati sottoposti a controllo incrociato con i monitor Masimo Rad-97 (errore assoluto medio 1,8% nell'intervallo SpO2 al 90-100%), i sensori di frequenza cardiaca sono stati convalidati rispetto all'elettrocardiografia (errore medio 2,3 bpm) e i contapassi sono stati calibrati su un protocollo di tapis roulant. I dati grezzi sono stati esportati come file CSV con data e ora a un intervallo di campionamento di 10 minuti.
4. Giustificazione della dimensione del campione
Soglie di 300 campioni di addestramento per XGBoost e circa 1000 sequenze di 7 giorni non sovrapposte per LSTM sono state supportate da esperimenti e analisi di potenza/simulazione. Un'analisi post hoc della potenza utilizzando G*Power (versione 3.1) ha indicato che 300 campioni hanno fornito >80% di potenza per rilevare una dimensione dell'effetto media (f2 di Cohen = 0,15) a α = 0,05 nella regressione logistica. Le simulazioni hanno mostrato che erano necessarie circa 1000 sequenze per una convergenza stabile della LSTM su serie temporali fisiologiche multivariate. L'evidenza empirica è riassunta nella Tabella 1 e visualizzata nella Figura 3 e nella Figura 4.
5. Pre-elaborazione dei dati
Le voci mancanti, identificate come spazi vuoti, valori sentinella o NaN, sono state imputate utilizzando la media del set di addestramento per ciascuna funzione per evitare perdite target:
(1)
dove mj è il numero di valori osservati per la funzione j. Lo stesso
è stato applicato ai set di convalida e di test. Per rimuovere il bias correlato alla scala, le funzionalità sono state standardizzate con la normalizzazione z-score utilizzando parametri di sola formazione:
(2)
dove μj e σj sono la media e la deviazione standard della funzione j stimata dai dati di addestramento. Per i segnali delle serie temporali (SpO2, frequenza cardiaca, passi), i flussi sono stati allineati a una cadenza di 10 minuti; Gli intervalli brevi fino a 30 minuti sono stati riempiti in avanti e gli intervalli più lunghi sono stati attenuati con una media mobile a tre punti. È stata quindi applicata una finestra scorrevole di 7 giorni con 1008 intervalli temporali per formare sequenze non sovrapposte per evitare perdite (vedi Figura 2).
6. Formazione del modello
Per i dati clinici strutturati, XGBoost è stato ottimizzato tramite la ricerca a griglia rispetto al tasso di apprendimento (0,01-0,3), alla profondità massima (3-10) e al numero di stimatori (100-500) in base a un obiettivo logistico binario (flusso di lavoro nella Figura 1). Per i dati di monitoraggio remoto, l'LSTM comprendeva due livelli LSTM nascosti con 128 unità ciascuno, inizializzazione Xavier, un tasso di dropout di 0,5 e un livello di uscita sigmoide; l'ottimizzazione ha utilizzato Adam con una velocità di apprendimento di 0,001 implementata in TensorFlow (pipeline nella Figura 2). L'overfitting è stato mitigato con l'abbandono e l'arresto anticipato (pazienza = 10) e lo squilibrio di classe è stato risolto con SMOTE.
7. Strategia di valutazione
I dati strutturati sono stati valutati con una convalida incrociata stratificata di 10 volte. I dati delle serie temporali sono stati valutati con la convalida walk-forward per preservare l'ordine temporale. Le metriche includevano accuratezza, precisione, richiamo, misura F, area sotto la curva ROC (AUC) e errore quadratico medio (MSE). Le differenze tra i modelli sono state valutate utilizzando il test del rango firmato di Wilcoxon (a due lati). Le cifre includono bande di confidenza o barre di errore del 95% per quantificare l'incertezza. Queste scelte affrontano direttamente i rischi di perdita di serie temporali e la necessità di significatività statistica richiesta dai revisori.
8. Rilevanza clinica e test pilota
Il modello ibrido ha previsto il declino della SpO2 6-12 ore prima della manifestazione clinica nel 78% dei pazienti monitorati e ha ridotto l'MSE del 42,5% rispetto ai valori basali. In un progetto pilota di quattro settimane con 20 pazienti con BPCO, i punteggi di rischio settimanali del modello concordavano con le valutazioni dei medici nell'85% delle revisioni, supportando il potenziale per l'integrazione clinica.
9. Software e ambiente
Le analisi sono state eseguite in Python 3.10.6 con scikit-learn 1.2.2, XGBoost 1.7.5, TensorFlow 2.13 e PyTorch 1.13 su Ubuntu 20.04 LTS con GPU NVIDIA RTX 3080, CUDA 11.6 e cuDNN 8.2. Il controllo completo delle versioni e i fornitori sono elencati nella tabella dei materiali non numerata.