Este estudio fue aprobado por el Comité de Ética del Hospital de Mujeres y Niños de Jinhua (Nº 2024KY099). Se obtuvo el consentimiento informado de los tutores legales de los menores.
Datos de pacientes
Este estudio retrospectivo incluyó a niños con neumonía por Mycoplasma pneumoniae (MPP) hospitalizados en el Hospital Infantil Afiliado al Hospital de Mujeres y Niños de Jinhua desde noviembre de 2022 hasta noviembre de 2023. La MPP se diagnosticó basándose en síntomas respiratorios, hallazgos de imagen torácica y pruebas positivas de ADN o ARN por Mycoplasma pneumoniae . La MPP general (GMPP) se refiere a los niños que respondieron adecuadamente a la terapia estándar con macrólidos, mientras que la RMPP se definió como fiebre persistente, síntomas agravados o anomalías progresivas por imagen tras al menos 7 días de terapia estándar conmacrólidos 10. Los criterios de inclusión fueron edad <14 años, infección confirmada por Mycoplasma pneumoniae , síntomas respiratorios y hallazgos de imagen compatibles con neumonía. Los criterios de exclusión fueron infección mixta dentro de los 10 días posteriores al inicio de la enfermedad, leucemia, enfermedad pulmonar crónica, inmunodeficiencia, terapia inmunosupresora previa, ingreso durante la fase de recuperación o registros incompletos. La terapia estándar se refería a azitromicina 10 mg/kg una vez al día, por vía oral o intravenosa, con una dosis máxima de 500 mg/día. Los niños que no cumplían los criterios RMPP se clasificaban como MPP general (GMPP). Un total de 500 niños elegibles se dividieron aleatoriamente en una cohorte de entrenamiento (n = 375) y una cohorte de validación (n = 125) usando la versión 4.1.2 del software R (función: sample(), set.seed = 42). La cohorte de validación se reservó exclusivamente para la evaluación externa del modelo y no se utilizó en ninguna etapa de la construcción del modelo. Los datos de los historiales médicos electrónicos se extrajeron utilizando una plantilla de extracción estandarizada desarrollada a previo. Los campos de datos incluían demografía, presentación clínica, valores de laboratorio, hallazgos de imagen y registros de tratamiento. Dos investigadores formados extrajeron los datos de forma independiente y las discrepancias se resolvieron por consenso. Los hallazgos de la tomografía computarizada torácica fueron evaluados por radiólogos cegados ante la agrupación clínica y los resultados de laboratorio. La fiabilidad entre evaluadores para la interpretación de imágenes en TC se evaluó utilizando la estadística kappa de Cohen para los hallazgos binarios de imagen (tapones de esputo o derrame pleural: sí/no; consolidación pulmonar: sí/no). Finalmente, 129 niños fueron clasificados como RMPP y 371 como MPP general.
Variables
Se extrajeron variables demográficas, clínicas, de laboratorio e imágenes de los registros médicos electrónicos utilizando un formulario predefinido. Las variables clínicas incluyeron edad, sexo, duración de la fiebre antes del ingreso, temperatura corporal máxima e hipoxemia. La duración de la fiebre se definió como el intervalo desde el inicio de la fiebre hasta el ingreso. La temperatura máxima corporal fue la más alta registrada antes o dentro de las 24 horas posteriores a la admisión. La hipoxemia se definió como la saturación periférica de oxígeno arterial <92% en el aire ambiente o la necesidad de oxígeno suplementario. Las variables de laboratorio incluían WBC, HB, PLT, CRP, ALB, ALT, CK-MB, LDH, D-dímero, IL-6, IL-8, IL-10, IL-17, PCT y NE%. Se recogieron muestras de sangre venosa en ayunas en las 24 horas posteriores a la admisión, y el tiempo de muestreo relativo al inicio de la fiebre y el inicio del antibiótico se registró cuando estaba disponible. Mycoplasma pneumoniae Se utilizó pruebas de ADN o ARN para confirmación etiológica. El test cuantitativo de PCR en tiempo real por hisopo de garganta para el ADN o ARN de Mycoplasma pneumoniae fue positivo. No se incluyó la carga cuantitativa de ADN/ARN porque no estaban disponibles datos estandarizados para todos los pacientes. Se revisó el TAC de tórax realizado dentro de los 3 días antes o después del ingreso para detectar tapones de esputo, derrame pleural y consolidación pulmonar. "Tapones de esputo o derrame pleural" se registró como positivo cuando alguno de los dos hallazgos estaba presente. La consolidación pulmonar se definió como opacidad parenquimatosa segmentaria o lobar en la TC. Se utilizó un escáner CT de 256 cortes. Los niños se colocaron en posición supina y se realizó un escaneo axial de tórax de baja dosis con un detector de ancho adecuado durante la suspensión de la respiración o mientras los niños dormían. La cobertura de la ecografía se extendía desde el ápice pulmonar hasta la base pulmonar para cubrir todo el parénquima pulmonar. Los parámetros de barrido se establecieron de la siguiente manera: voltaje de la válvula 100 kVp, grosor de corte 5 mm, grosor de corte de reconstrucción 1,25 mm, matriz 512 × 512, tiempo de rotación de la pórtica 0,28 s, índice de ruido 12. Las imágenes se reconstruyeron utilizando algoritmos ASIR-V y DLIR. Cada grupo contenía imágenes reconstruidas con pesos ASIR-V del 20%, 50% y 80%, así como imágenes DLIR-L, DLIR-M y DLIR-H. Los datos de imagen adquiridos se importaban a la estación de trabajo. Todas las imágenes fueron revisadas de forma independiente y en doble ciego por dos radiólogos con más de cinco años de experiencia clínica. Evaluaron la presencia de signos pulmonares anormales y resumieron las principales características de imagen. Para los casos con interpretaciones discrepantes, se alcanzó un consenso mediante discusión mutua.
Análisis estadístico
Los análisis estadísticos se realizaron utilizando la versión 4.1.2 del software R. Las variables categóricas se presentan como n (%) y se comparan usando la prueba del qui-cuadrado o la prueba exacta de Fisher. Las variables continuas se expresan como media ± desviación estándar o mediana (rango intercuartílico), según la distribución. La normalidad se evaluó mediante la prueba de Shapiro-Wilk. Las variables normalmente distribuidas se compararon usando la prueba t de muestras independientes, y las variables distribuidas no normalmente usando la prueba U de Mann-Whitney. Antes de la construcción del modelo, la comparabilidad inicial entre las cohortes de entrenamiento y validación se confirmó usando pruebas χ2 o Mann-Whitney U para todos los predictores candidatos, sin mostrar diferencias significativas entre cohortes (todas las P > 0,05). En la cohorte de entrenamiento, se introdujeron variables con P < 0,05 en regresión logística univariante en la regresión logística multivariable. La multicolinealidad se evaluó utilizando factores de inflación de varianza (VIFs) calculados con el paquete "car"; todos los predictores seleccionados presentaron < de 5 de FIV, lo que indica que no hay multilinealidad sustancial. La suposición de linealidad para LDH se examinó usando transformación logarítmica y splines cúbicos restringidos (usando el paquete "rms" con 3 nudos situados en los percentiles 10, 50 y 90); la forma lineal se mantuvo cuando el ajuste del modelo no mejoró (prueba de razón de verosimilitud: P > 0,05 para el término spline no lineal, apoyando la especificación lineal). Se construyó un nomograma a partir del modelo multivariable final usando el paquete "rms" (versión 6.3-0; funciones: lrm para regresión logística, Predict para predicción de valores y nomograma para representación gráfica). La discriminación se evaluó utilizando curvas características de funcionamiento del receptor, área bajo la curva y índice C con intervalos de confianza del 95% (paquete "pROC"). La calibración se evaluó utilizando curvas de calibración con 1.000 remuestras bootstrap. (paquete "RMS": Calibrar función). El corte óptimo de predicción se determinó de forma independiente dentro de cada cohorte utilizando el índice de Youden mediante el paquete pROC: 0,222 en la cohorte de entrenamiento y 0,247 en la cohorte de validación. El corte derivado del entrenamiento (0,222) también se aplicó a la cohorte de validación para la comparación de rendimiento entre cohortes, coherente con la práctica estándar interna de validación. El intervalo de confianza del 95% para la probabilidad predicha, por ejemplo, se calculó utilizando errores estándar basados en modelos en la escala logit. Se utilizó el análisis de la curva de decisión para evaluar el beneficio neto a través de probabilidades umbral (paquete "rmda": funciones decision_curve y plot_decision_curve). Se consideró estadísticamente significativo un P bilateral < 0,05. Los scripts R completos, incluyendo preprocesamiento de datos, construcción de modelos, validación y generación de nomogramas, están disponibles del autor correspondiente bajo solicitud.
Los datos faltantes eran mínimos (todas las variables <2% faltantes) y se manejaron mediante imputación mediana para variables continuas y imputación de modo para variables categóricas antes del análisis. Todos los predictores categóricos se codificaron como variables indicadoras binarias (0/1). Los predictores continuos (duración de la fiebre, temperatura máxima, LDH) se mantuvieron en sus unidades clínicas originales sin categorización. No se realizó selección de variables basada en el cribado univariado antes del modelado multivariable; en su lugar, se consideraron todos los candidatos clínicamente relevantes identificados en la literatura, y aquellos con P < 0,05 en el análisis univariante de la cohorte de entrenamiento avanzaron al modelo multivariable.
Implementación de calculadoras web
El nomograma dinámico basado en web se desplegó usando el paquete "shiny" (versión 1.7.4) en R. La interfaz de usuario se construyó con shiny::fluidPage(), shiny::sidebarLayout() y shiny::sliderInput() para variables continuas (duración de la fiebre, temperatura pico, LDH) y shiny::selectInput() para variables binarias (tapones de esputo o derrame pleural, consolidación pulmonar, hipoxemia). La lógica del servidor invocaba la función predict() del modelo lrm ajustado para calcular estimaciones individuales de riesgo, con las probabilidades resultantes representadas mediante renderPlot() y renderText(). La aplicación se alojaba en ShinyApps.io (https://predictrmpp.shinyapps.io/RMPP/).