Это исследование было одобрено Этическим комитетом Женской и детской больницы Цзиньхуа (No 2024KY099). Информированное согласие было получено от законных опекунов детей.
Данные пациентов
Это ретроспективное исследование включало детей с пневмонией Mycoplasma pneumoniae (MPP), госпитализированных в Детской больнице, аффилированной с Женской и детской больницей Цзиньхуа с ноября 2022 по ноябрь 2023 года. MPP был диагностирован на основе респираторных симптомов, результатов грудной визуализации и положительного теста Mycoplasma pneumoniae ДНК или РНК. Общий МПП (ГМПП) относится к детям, которые адекватно реагировали на стандартную макролидную терапию, тогда как РМПП определялся как устойчивая лихорадка, усугублённые симптомы или прогрессирующие нарушения визуализации после как минимум 7 дней стандартной терапиимакролидами (10). Критерии включения включали возраст <14 лет, подтверждённую инфекцию Mycoplasma pneumoniae , респираторные симптомы и результаты визуализации, соответствующие пневмонии. Критериями исключения были смешанная инфекция в течение 10 дней после начала заболевания, лейкемия, хроническое заболевание лёгких, иммунодефицит, предыдущая иммуносупрессивная терапия, госпитализация на этапе восстановления или неполные записи. Стандартная терапия включала азитромицин 10 мг/кг один раз в день, перорально или внутривенно, с максимальной дозой 500 мг/день. Дети, не соответствующие критериям RMPP, классифицировались как общие MPP (GMPP). Всего 500 подходящих детей были случайным образом разделены на обучающую когорту (n = 375) и валидационную когорту (n = 125) с использованием программы R версии 4.1.2 (функция: sample(), set.seed = 42). Когорта валидации была зарезервирована исключительно для внешней оценки модели и не использовалась ни на одном этапе построения модели. Данные электронных медицинских карт извлекались с использованием стандартизированного шаблона экстракции, разработанного заранее. Области данных включали демографию, клинические показатели, лабораторные показатели, результаты визуализации и записи о лечении. Два обученных исследователя независимо извлекли данные, и расхождения были устранены консенсусом. Результаты КТ грудной клетки оценивались радиологами, не обращая внимания на клинические группировки и лабораторные результаты. Надёжность интерпретации КТ-изображений между оценщиками оценивалась с использованием статистики каппа Коэна для бинарных результатов визуализации (мокроточные пробки или плевральный выпот: да/нет; консолидация лёгких: да/нет). Наконец, 129 детей были классифицированы как RMPP и 371 — как общий MPP.
Переменные
Демографические, клинические, лабораторные и визуализационные переменные извлекались из электронных медицинских записей с использованием заранее определённой формы. Клинические переменные включали возраст, пол, продолжительность температуры до госпитализации, пиковую температуру тела и гипоксемию. Продолжительность лихорадки определялась как интервал от начала лихорадки до госпитализации. Максимальная температура тела была самой высокой зафиксированной до или в течение 24 часов после госпитализации. Гипоксемия определялась как насыщение кислородом периферических артерий <92% в комнатном воздухе или необходимость дополнительного кислорода. Лабораторные переменные включали WBC, HB, PLT, CRP, ALB, ALT, CK-MB, LDH, D-димер, IL-6, IL-8, IL-10, IL-17, PCT и NE%. Образцы венозной крови натощак были взяты в течение 24 часов после поступления, а время отбора проб относительно начала лихорадки и начала приёма антибиотиков фиксировалось по мере доступности. Mycoplasma pneumoniae Для этиологического подтверждения использовался ДНК или РНК-тестирование. Количественный ПЦР-анализ Mycoplasma pneumoniae (РНК) в реальном времени показал положительный результат. Количественная нагрузка ДНК/РНК не была включена, поскольку стандартизированные данные о нагрузке были недоступны для всех пациентов. КТ грудной клетки, проведённая в течение 3 дней до или после госпитализации, была проверена на наличие пробок мокроты, плеврального выпота и консолидации лёгких. «Пробки мокроты или плевральный выпот» были отмечены как положительные при наличии любого из находок. Консолидация легких определялась как сегментарная или лобарная паренхимальная непрозрачность на КТ. Использовался КТ-сканер размером 256 срезов. Детей помещали в положение лежа на спине, а низкодозовое аксиальное сканирование грудной клетки проводилось с помощью детектора подходящей ширины при задержании дыхания или во время сна. Покрытие сканирования простиралось от вершины лёгких до основания лёгких, охватывая всю паренхиму лёгких. Параметры сканирования были установлены следующим образом: напряжение трубки 100 кВп, толщина среза 5 мм, толщина реконструкционного среза 1,25 мм, матрица 512 × 512, время вращения портала 0,28 с, индекс шума 12. Изображения были воссозданы с использованием алгоритмов ASIR-V и DLIR. Каждая группа содержала реконструированные изображения с весами ASIR-V 20%, 50% и 80%, а также изображения DLIR-L, DLIR-M и DLIR-H. Полученные данные изображений импортировались в рабочую станцию. Все изображения были независимо проверены двумя радиологами с более чем пятилетним клиническим опытом. Они оценили наличие аномальных признаков лёгких и резюмировали основные характеристики визуализации. В случаях с разными интерпретациями достигался консенсус путем взаимного обсуждения.
Статистический анализ
Статистический анализ проводился с использованием программного обеспечения R версии 4.1.2. Категориальные переменные представлены как n (%) и сравниваются с помощью теста хи-квадрата или точного теста Фишера. Непрерывные переменные выражаются как среднее ± стандартное отклонение или медиана (межквартильный диапазон) в зависимости от распределения. Нормальность оценивалась с помощью теста Шапиро-Уилка. Нормально распределённые переменные сравнивались с помощью теста t-независимых выборок, а ненормально распределенные переменные — с помощью теста Mann-Whitney U. До построения модели базовая сопоставимость между когортами обучения и валидации подтверждалась с помощью χ2 или тестов Манна-Уитни U для всех кандидатов-предикторов, при этом не было выявлено значимых различий между когортами (все P > 0,05). В обучающей когорте переменные с P < 0,05 в одномерной логистической регрессии были внесены в многомерную логистическую регрессию. Мультиколлинеарность оценивалась с использованием коэффициентов инфляции дисперсии (VIF), рассчитанных с помощью пакета «автомобиль»; все отобранные предикторы имели VIF < 5, что указывает на отсутствие значительной мультиколлинеарности. Предположение линейности для LDH изучалось с помощью логарифмического преобразования и ограниченных кубических сплайнов (с использованием пакета «rms» с 3 узлами, расположенными на 10-м, 50-м и 90-м перцентилях); линейная форма сохранялась, когда модель не улучшалась (тест на отношение правдоподобия: P > 0,05 для нелинейного сплайн-члена, поддерживающий линейную спецификацию). Номограмма была построена на основе окончательной многомерной модели с использованием пакета «rms» (версия 6.3-0; функции: lrm для логистической регрессии, Predict для прогнозирования значений и номограмма для графического представления). Дискриминация оценивалась с использованием кривых характеристик приёмника, площади под кривой и C-индекса с 95% доверительных интервалов («пакет pROC»). Калибровка оценивалась с помощью калибровочных кривых с 1000 повторных выборок bootstrap. ("пакет RMS: функция калибровки"). Оптимальный порог прогноза определялся независимо внутри каждой когорты с использованием индекса Юдена через пакет pROC: 0,222 в обучающей когорте и 0,247 в валидационной когорте. Порог, полученный на основе обучения (0,222), также применялся к когорте валидации для сравнения межкогортной производительности, что соответствует стандартной практике внутренней валидации. Например, 95% доверительный интервал для предсказанной вероятности был рассчитан с использованием стандартных ошибок на основе моделей по шкале логита. Анализ кривой принятия решения использовался для оценки чистой выгоды по пороговым вероятностям (пакет rmda: функции decision_curve и plot_decision_curve). Двусторонний коэффициент P < 0,05 считался статистически значимым. Полные R-скрипты, включая предобработку данных, построение модели, валидацию и генерацию номограмм, доступны у соответствующего автора по запросу.
Отсутствующие данные были минимальными (все переменные <2% отсутствуют) и обрабатывались с помощью медианной импутации для непрерывных переменных и модовой импутации для категориальных переменных до анализа. Все категориальные предикторы кодировались как бинарные (0/1) индикаторные переменные. Непрерывные предикторы (продолжительность лихорадки, пиковая температура, LDH) сохранялись в исходных клинических единицах без категоризации. Выбор переменных на основе одномерного скрининга до многомерного моделирования не проводился; вместо этого рассматривались все клинически значимые кандидаты, выявленные в литературе, а те, кто с P < 0,05 в одномерном анализе учебной когорты, были переведены в многомерную модель.
Реализация веб-калькулятора
Веб-ориентированная динамическая номограмма была развернута с помощью пакета «shiny» (версия 1.7.4) в R. Пользовательский интерфейс был построен с использованием shiny::fluidPage(), shiny::sidebarLayout() и shiny::sliderInput() для непрерывных переменных (продолжительность лихорадки, пиковая температура, LDH) и shiny::selectInput() для бинарных переменных (пробки мокроты или плевральный выпот, консолидация лёгких, гипоксемия). Серверная логика вызывала функцию predict() из модели Fitted LRM для вычисления индивидуальных оценок риска, при этом полученные вероятности отображались с помощью renderPlot() и renderText(). Приложение было размещено на ShinyApps.io (https://predictrmpp.shinyapps.io/RMPP/).