Cette étude a été examinée et approuvée par le comité d’éthique de l’hôpital de cancérologie de Taizhou. Le consentement éclairé écrit de tous les participants a été obtenu avant la collecte des données, conformément aux normes éthiques institutionnelles et nationales. Les réactifs et les logiciels utilisés sont répertoriés dans la table des matériaux.
1. Vue d’ensemble du flux de travail
Le flux de travail de bout en bout pour le traitement des données cliniques structurées avec XGBoost est illustré à la figure 1 et comprend l’ingestion de fonctionnalités, la construction d’arbres, les mises à jour résiduelles, la régularisation et l’évaluation. Le flux de travail de prétraitement et de modélisation de séquence pour les données de surveillance à distance avec LSTM est illustré à la figure 2, englobant le rééchantillonnage en 10 minutes, la gestion des écarts, le fenêtrage sans chevauchement sur 7 jours, l’architecture réseau et l’inférence.
2. Recrutement des patients et données démographiques
Entre mars 2023 et janvier 2024, les patients ont été recrutés consécutivement dans les cliniques externes et les services d’hospitalisation de l’hôpital de cancérologie de Taizhou. Les critères d’inclusion étaient un diagnostic confirmé de BPCO selon les critères GOLD (VEMS 1/CVF < 70 %), une condition stable au moment de l’inscription et un âge compris entre 40 et 80 ans. Les critères d’exclusion comprenaient des comorbidités graves telles que le cancer du poumon, une maladie cardiovasculaire non contrôlée ou une déficience cognitive empêchant le consentement éclairé. Au total, 214 patients ont été recrutés (âge moyen de 63,7 ± 8,9 ans ; 68 % d’hommes ; 54 % de fumeurs actuels ou anciens).
3. Collecte des données
Les données cliniques structurées comprenaient la fonction pulmonaire (VEMS1, CVF), la saturation en oxygène du sang (SpO2), la durée du séjour à l’hôpital et la durée de la maladie. La fonction pulmonaire et la SpO2 ont été enregistrées en pourcentage, en jours de séjour à l’hôpital et en durée de la maladie en années. Les données de surveillance à distance ont été collectées à l’aide d’oxymètres portables certifiés et de trackers d’activité validés par rapport aux normes d’excellence des hôpitaux : les oxymètres de pouls ont été vérifiés avec les moniteurs Masimo Rad-97 (erreur absolue moyenne de 1,8 % dans la plage de 90 à 100 % de SpO2 ), les capteurs de fréquence cardiaque ont été validés par rapport à l’électrocardiographie (erreur moyenne de 2,3 bpm) et les compteurs de pas ont été calibrés sur un protocole de tapis roulant. Les données brutes ont été exportées sous forme de fichiers CSV horodatés à un intervalle d’échantillonnage de 10 minutes.
4. Justification de la taille de l’échantillon
Les seuils de 300 échantillons d’apprentissage pour XGBoost et d’environ 1000 séquences de 7 jours non chevauchantes pour LSTM ont été étayés par des expériences et des analyses de puissance/simulation. Une analyse de puissance a posteriori à l’aide de G*Power (version 3.1) a indiqué que 300 échantillons fournissaient une puissance de >80 % pour détecter une taille d’effet moyenne (f2 de Cohen = 0,15) à α = 0,05 en régression logistique. Les simulations ont montré qu’environ 1000 séquences étaient nécessaires pour une convergence stable du LSTM sur des séries temporelles physiologiques multivariées. Les preuves empiriques sont résumées dans le tableau 1 et visualisées dans les figures 3 et 4.
5. Prétraitement des données
Les entrées manquantes, identifiées comme des blancs, des valeurs sentinelles ou NaN, ont été imputées à l’aide de la moyenne de l’ensemble d’apprentissage pour chaque caractéristique afin d’éviter les fuites cibles :
(1)
où mj est le nombre de valeurs observées pour la caractéristique j. Il en a été de même
pour les kits de validation et de test. Pour éliminer le biais lié à l’échelle, les caractéristiques ont été normalisées avec la normalisation du score z à l’aide de paramètres d’entraînement uniquement :
(2)
où μj et σj sont la moyenne et l’écart-type de la caractéristique j estimés à partir des données d’entraînement. Pour les signaux de séries temporelles (SpO2, fréquence cardiaque, pas), les flux ont été alignés sur une cadence de 10 minutes ; Les écarts courts jusqu’à 30 minutes ont été comblés vers l’avant, et les écarts plus longs ont été lissés avec une moyenne mobile à trois points. Une fenêtre glissante de 7 jours avec 1008 pas de temps a ensuite été appliquée pour former des séquences qui ne se chevauchent pas afin d’éviter les fuites (voir figure 2).
6. Formation du modèle
Pour les données cliniques structurées, XGBoost a été ajusté via la recherche par grille sur le taux d’apprentissage (0,01-0,3), la profondeur maximale (3-10) et le nombre d’estimateurs (100-500) dans le cadre d’un objectif logistique binaire (flux de travail dans la figure 1). Pour les données de surveillance à distance, le LSTM comprenait deux couches LSTM cachées de 128 unités chacune, l’initialisation Xavier, un taux d’abandon de 0,5 et une couche de sortie sigmoïde ; l’optimisation a utilisé Adam avec un taux d’apprentissage de 0,001 implémenté dans TensorFlow (pipeline dans la figure 2). Le surajustement a été atténué par l’abandon et l’arrêt précoce (patience = 10), et le déséquilibre de classe a été corrigé avec SMOTE.
7. Stratégie d’évaluation
Les données structurées ont été évaluées par validation croisée stratifiée en 10 parties. Les données de séries chronologiques ont été évaluées à l’aide d’une validation préalable afin de préserver l’ordre temporel. Les mesures comprenaient l’exactitude, la précision, le rappel, la mesure F, l’aire sous la courbe ROC (AUC) et l’erreur quadratique moyenne (MSE). Les différences entre les modèles ont été évaluées à l’aide du test de Wilcoxon (bilatéral). Les chiffres comprennent des bandes de confiance à 95 % ou des barres d’erreur pour quantifier l’incertitude. Ces choix abordent directement les risques de fuite de séries chronologiques et le besoin de signification statistique demandé par les examinateurs.
8. Pertinence clinique et essais pilotes
Le modèle hybride a prédit le déclin de la SpO2 6 à 12 heures avant la manifestation clinique chez 78 % des patients suivis et a réduit la MSE de 42,5 % par rapport aux lignes de base. Dans le cadre d’un projet pilote de quatre semaines mené auprès de 20 patients atteints de BPCO, les scores de risque hebdomadaires du modèle concordaient avec les évaluations des médecins dans 85 % des évaluations, ce qui soutient le potentiel d’intégration clinique.
9. Logiciel et environnement
Les analyses ont été effectuées en Python 3.10.6 avec scikit-learn 1.2.2, XGBoost 1.7.5, TensorFlow 2.13 et PyTorch 1.13 sur Ubuntu 20.04 LTS avec un GPU NVIDIA RTX 3080, CUDA 11.6 et cuDNN 8.2. La version complète et les fournisseurs sont répertoriés dans la table des matériaux non numérotée.