El protocolo del estudio fue aprobado por el Comité de Ética para la Investigación Biomédica que Involucra a Seres Humanos del Primer Hospital Popular de Jiashan, Hospital Afiliado a la Universidad de Jiaxing (número de aprobación: 2026 Research Approval No. 005; número de aceptación: LW2026004; fecha de aprobación: 13 de enero de 2026). Se dispensó el requisito de consentimiento informado individual porque este estudio retrospectivo utilizó datos clínicos con identificación eliminada. Antes del análisis, se eliminaron los identificadores directos, incluyendo nombre, número de identificación hospitalario, número de teléfono, dirección y número de identificación nacional. El manejo de los datos siguió los requisitos institucionales de protección de datos y la Declaración de Helsinki10.
Diseño del estudio y ensamblaje de la cohorte multicéntrica
Este estudio de cohorte retrospectivo multicéntrico utilizó datos bioquímicos y de coagulación al ingreso para construir un flujo de trabajo de puntuación reproducible para la estratificación temprana del riesgo de síndrome de disfunción orgánica múltiple (MODS) en pancreatitis aguda (AP). La AP se definió por al menos dos de los siguientes criterios: dolor abdominal típico, amilasa o lipasa sérica al menos tres veces por encima del límite superior normal, o hallazgos de imagen compatibles con AP. La ventana basal se definió como las primeras 24 horas tras el ingreso hospitalario. Para cada predictor de laboratorio, se utilizó el valor elegible más temprano dentro de esta ventana. No se utilizaron valores obtenidos tras el primer registro del inicio de MODS para la construcción de los predictores.
Se revisaron hospitalizaciones consecutivas por PA procedentes de los centros participantes desde el 1 de enero de 2023 hasta el 20 de diciembre de 2025, utilizando un registro unificado de cribado, una lista de elegibilidad, un diccionario de datos y un formulario de adjudicación de eventos finales. El registro de cribado registró el número total de ingresos revisados, ingresos duplicados, exclusiones por PA no confirmada, edad menor de 18 años, pancreatitis crónica, pancreatitis asociada a malignidad pancreática, hora de ingreso no fiable, datos de laboratorio a las 24 h no disponibles y registros de eventos finales incompletos, así como la cohorte analítica final. El registro de cribado bloqueado se utilizó para generar Figura 1. Los conteos de pacientes en la Figura 1 se verificaron frente al conjunto de datos analítico final antes del análisis.
Después de bloquear la cohorte analítica, los pacientes se asignaron a una cohorte de desarrollo y a una cohorte de validación reservada en una proporción de 7:3. La división se realizó en R versión 4.3.2 utilizando el paquete caret versión 6.0-94, con la semilla aleatoria fijada en 20260420. Se realizó un muestreo estratificado según el estado de MODS y el centro. Si la distribución por centro difería en más de 5 puntos porcentuales entre cohortes, se utilizaron estratos combinados de centro y MODS. La asignación de cohortes se guardó antes del ajuste del modelo y no se modificó posteriormente. Las reglas de preprocesamiento, configuraciones de imputación, ajuste del modelo, estimación de coeficientes, construcción del puntaje y selección del punto de corte se derivaron exclusivamente a partir de la cohorte de desarrollo. La cohorte de validación reservada se utilizó únicamente para la evaluación del modelo bloqueado.
Criterios de elegibilidad y procedimientos de inscripción
Los pacientes elegibles eran adultos de 18 años o más con diagnóstico confirmado de PA, tiempo de ingreso hospitalario rastreable, al menos un panel de laboratorio bioquímico o de coagulación dentro de las 24 horas posteriores al ingreso y registros hospitalarios completos para la adjudicación de los desenlaces. Se excluyeron a los pacientes si no se podía confirmar la PA, si presentaban pancreatitis por pancreatitis crónica o malignidad pancreática asociada, o si no se podía reconstruir de forma confiable el momento del primer ingreso hospitalario. Se excluyeron los casos de traslado cuando el momento del primer ingreso, el momento de la primera toma de laboratorio o el estado de disfunción orgánica temprana no eran rastreables entre instituciones. Para ingresos repetidos del mismo paciente durante el período del estudio, solo se conservó la primera hospitalización elegible. Para cada caso excluido se registró una causa principal de exclusión.
Resultados y adjudicación de los puntos finales
El resultado principal fue la SDOM durante la hospitalización, definida como disfunción de dos o más sistemas orgánicos durante la estancia hospitalaria. La disfunción orgánica se evaluó a partir de los registros respiratorios, cardiovasculares, renales, hepáticos, de coagulación y neurológicos en la historia clínica electrónica. Se aplicaron los mismos criterios preespecificados por sistema orgánico en todos los centros antes del desarrollo del modelo, y los criterios operativos utilizados para la adjudicación de los puntos finales se resumen en la Tabla Suplementaria 1.
Dos revisores capacitados evaluaron independientemente cada registro desidentificado utilizando notas de ingreso, notas de evolución, registros de unidad de cuidados intensivos, registros de procedimientos, informes de laboratorio, registros de signos vitales y resúmenes de alta. El formulario de adjudicación registró el código de estudio del paciente, el sistema orgánico afectado, el registro fuente, la fecha y hora de inicio, el valor de laboratorio o fisiológico de apoyo, la decisión del revisor y el estado final de adjudicación. Se aceptó como definitivo el acuerdo entre los revisores. Las discrepancias se resolvieron mediante un tercer revisor senior. Los resultados secundarios incluyeron ingreso a la unidad de cuidados intensivos, ventilación mecánica invasiva o no invasiva, terapia de reemplazo renal, mortalidad intrahospitalaria, necrosis infectada o procedimiento de drenaje y duración de la estancia hospitalaria. La duración de la estancia se calculó como la fecha de alta menos la fecha de ingreso más 1 día.
Predictores candidatos al momento del ingreso
Los predictores candidatos se limitaron a variables de rutina disponibles dentro de las 24 horas posteriores al ingreso. Para las pruebas repetidas dentro del período de 24 horas, se conservó el valor elegible más temprano. Las variables bioquímicas incluyeron recuento de glóbulos blancos, proteína C reactiva, procalcitonina cuando estuvo disponible, lactato, glucosa, calcio total, nitrógeno ureico en sangre, creatinina, alanina aminotransferasa, bilirrubina total, albúmina y deshidrogenasa láctica. Las variables de coagulación incluyeron recuento de plaquetas, tiempo de protrombina, relación normalizada internacional del tiempo de protrombina, tiempo de tromboplastina parcial activada, fibrinógeno y D-dímero.
Se recopilaron variables clínicas basales para la descripción del perfil de pacientes y el cálculo del puntaje de referencia. Estas variables incluyeron edad, sexo, índice de masa corporal, tiempo transcurrido desde el inicio de los síntomas hasta la admisión, etiología de la pancreatitis aguda, estado de tabaquismo, antecedentes de pancreatitis previa, comorbilidades, frecuencia cardíaca al ingreso, presión arterial media, estado de síndrome de respuesta inflamatoria sistémica, puntaje del Índice al Lado de la Cama para la Gravedad en la Pancreatitis Aguda, puntaje de Evaluación de la Fisiología Aguda y la Salud Crónica II y puntaje de Evaluación Secuencial de Insuficiencia Orgánica. Cuando tanto el tiempo de protrombina como el tiempo de protrombina con razón normalizada internacional estaban disponibles, se utilizó el tiempo de protrombina con razón normalizada internacional para la modelización, y se conservó el tiempo de protrombina para los informes descriptivos. El archivo de extracción contenía una fila por cada hospitalización y conservaba el código del centro, el código del paciente en el estudio, la hora de admisión, la hora de muestreo, el valor crudo, la unidad cruda, el valor armonizado, la unidad armonizada, la tabla de origen y el indicador de datos faltantes.
Armonización de laboratorios entre centros
Todas las variables de laboratorio se convirtieron a unidades objetivo preespecificadas antes del análisis: recuento de glóbulos blancos y recuento de plaquetas, 109/L; proteína C reactiva, mg/L; procalcitonina, ng/mL; lactato, glucosa, calcio total y nitrógeno ureico en sangre, mmol/L; creatinina y bilirrubina total, µmol/L; alanina aminotransferasa y deshidrogenasa láctica, U/L; albúmina y fibrinógeno, g/L; tiempo de protrombina y tiempo de tromboplastina parcial activada, segundos; relación internacional normalizada del tiempo de protrombina, relación adimensional; y D-dímero, mg/L de unidades equivalentes a fibrinógeno. Los valores de D-dímero reportados como µg/mL de unidades equivalentes a fibrinógeno se convirtieron a mg/L de unidades equivalentes a fibrinógeno mediante una conversión numérica 1:1. Los valores de D-dímero reportados como unidades de D-dímero se convirtieron únicamente cuando el laboratorio local proporcionó un factor de conversión validado. Los valores sin un factor de conversión validado se consideraron faltantes en el modelado.
Se revisaron las unidades de reporte, rangos de referencia, registros de ensayos y notas del analizador de cada centro antes de agruparlos. El archivo de armonización conservó la unidad original, la unidad objetivo, la regla de conversión y el revisor responsable de verificar la conversión. Los valores de laboratorio no se normalizaron según los rangos de referencia específicos de cada centro en el análisis principal, ya que la puntuación tenía la intención de utilizar valores clínicos absolutos. El código del centro se mantuvo durante la imputación y el análisis de sensibilidad. Se marcaron los valores negativos para pruebas que no pueden ser negativas, unidades incompatibles, marcas de tiempo fuera de la ventana de 24 horas y valores no respaldados por los registros originales. Solo se corrigió un valor cuando el registro fuente confirmó un error de transcripción o de conversión de unidades; de lo contrario, se estableció como faltante.
Preprocesamiento de datos y manejo de datos faltantes
La ausencia de datos se resumió por variable, centro y estado de MODS antes del ajuste del modelo. Los predictores candidatos se excluyeron del desarrollo del modelo si la ausencia general superaba el 30 % o si la ausencia superaba el 50 % en algún centro individual. Los pacientes se excluyeron del conjunto de datos de modelado si faltaba más del 40 % de los predictores candidatos o si no se podía determinar el resultado principal.
Los predictores continuos se ajustaron mediante winsorización en los percentiles 1.º y 99.º estimados en la cohorte de desarrollo. Los mismos umbrales se aplicaron a la cohorte de validación reservada. No se volvieron a estimar los umbrales en los datos de validación. Las variables con asimetría mayor a 2 tras la winsorización se transformaron utilizando el logaritmo natural; se usó log(x + 0,01) para variables que podían contener cero. Los predictores continuos se estandarizaron utilizando la media y la desviación estándar de la cohorte de desarrollo, y los mismos parámetros se aplicaron a los datos de validación.
Se realizó imputación múltiple utilizando ecuaciones encadenadas en R con el paquete mice versión 3.16.0. Se generaron veinte conjuntos de datos imputados con 20 iteraciones por conjunto de datos. Se utilizó el emparejamiento predictivo de medias para variables continuas, regresión logística para variables binarias y regresión logística politómica para variables categóricas nominales. El modelo de imputación incluyó predictores candidatos, código del centro, indicador de cohorte y estado de MODS para el desarrollo del modelo. Los resultados de validación no se utilizaron para modificar las reglas de preprocesamiento, la elegibilidad de variables ni los umbrales. Las estimaciones de regresión se combinaron utilizando las reglas de Rubin cuando se informaron. Las métricas de predicción se calcularon dentro de cada conjunto de datos imputados y se resumieron como estimaciones combinadas o medianas con rangos. Se prespecificaron dos análisis de sensibilidad: análisis de casos completos utilizando pacientes sin predictores del modelo final faltantes y análisis del valor basal utilizando el peor valor de laboratorio elegible dentro de las 24 horas en lugar del valor más temprano.
Análisis estadístico y desarrollo del modelo
Las características basales se resumieron por separado en las cohortes de desarrollo y de validación de retención. Las variables continuas se presentaron como media ± desviación estándar o mediana con rango intercuartílico, y las variables categóricas como n (%). Las variables bioquímicas y de coagulación al ingreso se compararon entre pacientes con y sin SDOM. Se utilizó la prueba t de Student o la prueba de suma de rangos de Wilcoxon para variables continuas, y la prueba de chi-cuadrado o la prueba exacta de Fisher para variables categóricas. Estos valores P fueron descriptivos y no se utilizaron por sí solos para la selección de predictores.
El desarrollo del modelo se limitó a la cohorte de desarrollo. Se utilizó la regresión logística penalizada como modelo principal para reducir el sobreajuste ante recuentos limitados de eventos y predictores de laboratorio correlacionados. El modelo se ajustó en R versión 4.3.2 utilizando el paquete glmnet versión 4.1-8 con familia binomial, predictores estandarizados y penalización mediante el operador de contracción absoluta mínima y selección (LASSO), especificada mediante alpha = 1. Para el ajuste de la penalización se empleó una validación cruzada de diez pliegues repetida 10 veces, con la semilla aleatoria fijada en 20260420. El parámetro de penalización se seleccionó utilizando la regla de un error estándar. Se guardaron la trayectoria de los coeficientes, los valores del error de validación cruzada, el parámetro de penalización seleccionado y los predictores retenidos. El flujo de trabajo del análisis se muestra en Figura 2. El software, los paquetes estadísticos, los sistemas de bases de datos y las plataformas de laboratorio utilizados en el protocolo se enumeran en la Tabla de Materiales.
No se utilizó la regresión logística sin penalización como modelo inferencial principal tras la selección con penalización. Si se informaron razones de momios para la interpretación clínica, se etiquetaron como estimaciones descriptivas posteriores a la selección. La discriminación se estimó mediante el área bajo la curva de característica operativa del receptor, con intervalos de confianza del 95 % calculados utilizando el método de DeLong en el paquete pROC versión 1.18.5. La calibración se evaluó mediante el intercepto de calibración, la pendiente de calibración, la puntuación de Brier y gráficos de calibración. Se empleó remuestreo bootstrap con 1 000 repeticiones en la cohorte de desarrollo para estimar el rendimiento corregido por optimismo. El análisis de curva de decisión no se mantuvo en el conjunto de análisis revisado; por lo tanto, las afirmaciones sobre beneficio clínico neto se eliminaron del manuscrito.
Derivación e interpretación del puntaje
El puntaje compuesto bioquímico-coagulación basado en puntos de admisión se derivó de los coeficientes del modelo penalizado bloqueado. Las variables continuas se utilizaron durante el ajuste del modelo. Las categorías con umbrales se introdujeron únicamente después de la finalización del modelo para facilitar el cálculo al pie del lecho. Los umbrales se seleccionaron según rangos de referencia clínicos establecidos, umbrales comúnmente utilizados para anormalidades de laboratorio o puntos de corte basados en coeficientes definidos en la cohorte de desarrollo antes de la validación. Los umbrales no se optimizaron en la cohorte de validación.
Se asignaron puntos enteros según la dirección y la magnitud relativa de los coeficientes fijos. Puntos más altos indicaban perfiles de laboratorio más adversos. La puntuación total se calculó sumando los puntos de los componentes retenidos. Las categorías de riesgo preespecificadas se utilizaron únicamente para la estratificación del riesgo; no se utilizaron para definir el SDMO. La puntuación tenía como objetivo identificar a los pacientes que podrían necesitar una observación más estrecha, evaluaciones de laboratorio repetidas o preparación para vías de soporte orgánico. No se utilizó como una regla diagnóstica independiente, un desencadenante de anticoagulación ni como sustituto de la reevaluación clínica.
Validación con conjunto de retención y evaluación comparativa
Se aplicaron al conjunto de validación con retención las reglas fijas de preprocesamiento, la estructura de imputación, los umbrales de winsorización, los parámetros de estandarización, los coeficientes del modelo, las reglas de puntuación y los puntos de corte preespecificados. No se realizó ninguna selección adicional de predictores, reajuste del modelo, recalibración de coeficientes ni optimización de umbrales. El desempeño de validación se resumió mediante indicadores de discriminación, calibración, puntuación de Brier, sensibilidad, especificidad, valor predictivo positivo y valor predictivo negativo.
Se calcularon herramientas de gravedad establecidas a partir de la misma ventana de ingreso cuando los datos estaban disponibles, incluyendo el Índice de Gravedad en Pancreatitis Aguda en la Cama, la puntuación de Evaluación de la Fisiología Aguda y la Salud Crónica II, la puntuación de Evaluación Secuencial de Insuficiencia Orgánica, la puntuación de Ranson, la Puntuación de Pancreatitis Aguda Inofensiva y los criterios del síndrome de respuesta inflamatoria sistémica. Las comparaciones utilizaron el mismo punto final, población de análisis y marco de validación. Se realizaron comparaciones pareadas del área bajo la curva utilizando las pruebas de DeLong.
Todos los gráficos de validación y los gráficos estratificados por puntuación se generaron a partir del conjunto de datos analíticos bloqueado. El tamaño de la muestra en cada panel de la figura se verificó en comparación con la cohorte analítica final, la cohorte de desarrollo, la cohorte de validación independiente o el estrato de puntuación preespecificado antes de la elaboración de las figuras. No se utilizaron métodos visuales diseñados para datos ómicos de alta dimensión, datos de conteo ecológicos o análisis de abundancia diferencial para datos de laboratorio clínico de rutina.