Bu çalışma, Jinhua Kadın ve Çocuk Hastanesi Etik Komitesi tarafından onaylandı (No. 2024KY099). Bilgilendirilmiş onay, çocukların yasal vasilerinden alındı.
Hasta verileri
Bu retrospektif çalışma, Kasım 2022'den Kasım 2023'e kadar Jinhua Kadın ve Çocuk Hastanesi'ne bağlı Çocuk Hastanesi'nde hastaneye kaldırılan Mycoplasma pneumoniae pnömonisi (MPP) hastası çocukları içeriyordu. MPP, solunum semptomları, göğüs görüntüleme bulguları ve pozitif Mycoplasma pnömoniae DNA veya RNA testleri temelinde teşhis edildi. Genel MPP (GMPP), standart makrolid tedavisine yeterince yanıt veren çocukları ifade ederken, RMPP ise en az 7 gün standart makrolid tedavisinden sonra kalıcı ateş, şiddetli semptomlar veya ilerleyici görüntüleme anormallikleri olaraktanımlanmıştır 10. Dahil etme kriterleri <14 yaş, Mycoplasma pneumoniae enfeksiyonu, solunum semptomları ve pnömoni ile uyumlu görüntüleme bulguları idi. Dışlama kriterleri ise hastalık başlangıcından sonraki 10 gün içinde karma enfeksiyon, lösemi, kronik akciğer hastalığı, bağışıklık eksikliği, önceki immünsüpresif tedavi, iyileşme aşamasında yatış veya eksik kayıtlardı. Standart tedavi, günde bir kez ağızdan veya damar yoluyla 10 mg/kg azitromisin önerir, maksimum doz ise 500 mg/gündür. RMPP kriterlerini karşılamayan çocuklar genel MPP (GMPP) olarak sınıflandırıldı. Toplam 500 uygun çocuk, R yazılımı 4.1.2 sürümü (fonksiyon: sample(), set.seed = 42) kullanılarak rastgele bir eğitim grubu (n = 375) ve doğrulama kohortu (n = 125) olarak bölündü. Doğrulama grubu yalnızca dış model değerlendirmesi için ayrılmıştı ve model oluşturmanın hiçbir aşamasında kullanılmadı. Elektronik tıbbi kayıt verileri, önceden geliştirilmiş standart bir çıkarma şablonu kullanılarak çıkarıldı. Veri alanları arasında demografik bilgiler, klinik sunum, laboratuvar değerleri, görüntüleme bulguları ve tedavi kayıtları yer alıyordu. İki eğitimli araştırmacı verileri bağımsız olarak çıkardı ve tutarsızlıklar uzlaşmayla çözüldü. Göğüs BT bulguları, klinik gruplama ve laboratuvar sonuçlarına kör radyologlar tarafından değerlendirildi. BT görüntüleme yorumu için değerlendirme aracı arası güvenilirliği, ikili görüntüleme bulguları için Cohen'in kappa istatistiği kullanılarak değerlendirildi (balgam tıkacı veya plevral effuzyon: evet/hayır; akciğer konsolidasyonu: evet/hayır). Son olarak, 129 çocuk RMPP, 371 ise genel MPP olarak sınıflandırıldı.
Değişkenler
Demografik, klinik, laboratuvar ve görüntüleme değişkenleri, elektronik tıbbi kayıtlardan önceden tanımlanmış bir form kullanılarak çıkarıldı. Klinik değişkenler arasında yaş, cinsiyet, yatıştan önceki ateş süresi, zirve vücut sıcaklığı ve hipoksemi yer alıyordu. Ateş süresi, ateşin başlangıcından yatışa kadar olan aralık olarak tanımlanmıştır. En yüksek vücut sıcaklığı, kabul öncesi veya sonrasında 24 saat içinde kaydedilen en yüksek sıcaklıktı. Hipoksemi, periferik arteriyel oksijen doygunluğu <%92 oda havasında veya ek oksijen ihtiyacı olarak tanımlanmıştır. Laboratuvar değişkenleri arasında WBC, HB, PLT, CRP, ALB, ALT, CK-MB, LDH, D-dimer, IL-6, IL-8, IL-10, IL-17, PCT ve NE% vardı. Açlık venöz kan örnekleri kabul edildikten sonraki 24 saat içinde alındı ve ateş başlangıcı ile antibiyotik başlangıcına göre numune alma süresi kaydedildi. Mycoplasma pneumoniae Etiyolojik doğrulama için DNA veya RNA testi kullanıldı. Mycoplasma pneumoniae DNA veya RNA için boğaz sürüntüsü gerçek zamanlı nicel PCR testi pozitifti. Standartlaştırılmış yük verileri tüm hastalar için mevcut olmadığı için nicel DNA/RNA yükü dahil edilmedi. Kabul öncesi veya sonrası 3 gün içinde yapılan göğüs BT'si, balgam tıkacı, plevral efüzyon ve akciğer konsolidasyonu açısından incelendi. "Balgam tıkaları veya plevral efüzyon" herhangi bir bulgu bulunduğunda pozitif olarak kaydedildi. Akciğer konsolidasyonu, BT'de segmental veya lobar parenkimal opaklık olarak tanımlandı. 256 dilimli bir BT tarayıcı kullanıldı. Çocuklar sırtüstü pozisyona yerleştirildi ve düşük dozlu eksenel göğüs taraması, nefes tutma sırasında veya uyurken uygun genişlikte bir dedektörle yapıldı. Tarama kapsamı, akciğer ucundan akciğer tabanına kadar uzanarak tüm akciğer parenximasını kapsaydı. Tarama parametreleri şu şekilde ayarlandı: tüp voltajı 100 kVp, dilim kalınlığı 5 mm, yeniden yapılandırma dilimi kalınlığı 1,25 mm, matris 512 × 512, portal dönüş süresi 0,28 s, gürültü indeksi 12. Görüntüler ASIR-V ve DLIR algoritmaları kullanılarak yeniden oluşturuldu. Her grup, ASIR-V ağırlıkları %20, %50 ve %80 olan yeniden oluşturulmuş görüntülerin yanı sıra DLIR-L, DLIR-M ve DLIR-H görüntülerini içeriyordu. Elde edilen görüntü verileri iş istasyonuna aktarıldı. Tüm görüntüler, beş yıldan fazla klinik deneyime sahip iki radyolog tarafından çift kör bir şekilde bağımsız olarak incelendi. Anormal akciğer belirtilerinin varlığını değerlendirdiler ve ana görüntüleme özelliklerini özetlediler. Farklı yorumlara sahip vakalar için, karşılıklı tartışma yoluyla bir uzlaşmaya varıldı.
İstatistiksel analiz
İstatistiksel analizler R yazılımı 4.1.2 sürümü kullanılarak gerçekleştirilmiştir. Kategorik değişkenler n (%) olarak sunulur ve ki-kare testi veya Fisher'ın tam testi kullanılarak karşılaştırılır. Sürekli değişkenler, dağılıma göre ortalama ± standart sapma veya medyan (çeyreklerarası aralık) olarak ifade edilir. Normallik Shapiro-Wilk testi kullanılarak değerlendirildi. Normal dağılımlı değişkenler bağımsız örneklem t testiyle, normal dağılım olmayan değişkenler ise Mann-Whitney U testi kullanılarak karşılaştırıldı. Model oluşturulmadan önce, eğitim ve doğrulama kohortları arasındaki temel karşılaştırılabilirlik, tüm aday tahmin ediciler için χ2 veya Mann-Whitney U testleri kullanılarak doğrulanmış ve kohortlar arasında anlamlı farklar olmadığı gösterilmiştir (tüm P > 0.05). Eğitim grubunda, tek değişkenli lojistik regresyonda P < 0.05 olan değişkenler çok değişkenli lojistik regresyona girildi. Çok kollineerlik, "araba" paketiyle hesaplanan varyans enflasyon faktörleri (VIF) kullanılarak değerlendirildi; tüm seçilmiş tahmin edicilerin VIF < 5'e sahipti ve bu da belirgin bir multikollineerlik olmadığını gösteriyordu. LDH için doğrusal varsayımı, logaritmik dönüşüm ve kısıtlı kübik splinlar kullanılarak incelendi (10., 50. ve 90. pertiklerde 3 knotlu "rms" paketi kullanılarak); Model uyumu iyileşmediğinde doğrusal form korundu (olasılık oranı testi: P > doğrusal olmayan spline terimi için 0.05, doğrusal spesifikasyonu destekledi). Nihai çok değişkenli modelden "rms" paketi kullanılarak bir nomogram oluşturuldu (sürüm 6.3-0; fonksiyonlar: lojistik regresyon için lrm, değer tahmini için Predict ve grafiksel temsil için nomogram). Ayırtım, alıcı çalışma karakteristik eğrileri, eğrin altındaki alan ve C indeksi ile %95 güven aralıkları ("pROC" paketi) kullanılarak değerlendirildi. Kalibrasyon, 1.000 bootstrap yeniden örneklemesi ile kalibrasyon eğrileri kullanılarak değerlendirildi. ("RMS" paketi: Kalibrasyon fonksiyonu). Optimal tahmin kesme sınırı, her kohort içinde Youden indeksi kullanılarak pROC paketi üzerinden bağımsız olarak belirlendi: eğitim kohortunda 0.222, doğrulama kohortunda 0.247. Eğitimden türetilen sınır (0.222) ayrıca, standart iç doğrulama uygulamasıyla uyumlu olarak, çapraz kohort performans karşılaştırması için doğrulama kohortuna uygulandı. Örneğin, tahmin edilen olasılık için %95 güven aralığı, logit ölçeğinde model tabanlı standart hatalar kullanılarak hesaplandı. Karar eğrisi analizi, eşik olasılıklar arasında net faydayı değerlendirmek için kullanıldı ("rmda" paketi: decision_curve ve plot_decision_curve fonksiyonları). İki taraflı P < 0.05 istatistiksel olarak anlamlı kabul edildi. Veri ön işleme, model oluşturma, doğrulama ve nomogram üretimi dahil olmak üzere tam R betikleri, ilgili yazardan talep üzerine temin edilebilir.
Eksik veriler minimumdu (tüm değişkenler <%2 eksikti) ve analizden önce sürekli değişkenler için medyan atama ve kategorik değişkenler için mod ödümlemesi kullanılarak ele alındı. Tüm kategorik tahmin ediciler ikili (0/1) gösterge değişkenleri olarak kodlanmıştır. Sürekli tahmin ediciler (ateş süresi, tepe sıcaklık, LDH) orijinal klinik birimlerinde kategorize edilmeden korundu. Çok değişkenli modellemeden önce tek değişkenli taramaya dayalı bir değişken seçimi yapılmamıştır; bunun yerine, literatürde tanımlanan tüm klinik açıdan ilgili adaylar değerlendirilmiş ve eğitim kohortunun tek değişkenli analizinde P < 0.05 olanlar çok değişkenli modele yönlendirilmiştir.
Web hesaplayıcı uygulaması
Web tabanlı dinamik nomogram, R dilinde "parlak" paket (sürüm 1.7.4) kullanılarak dağıtıldı. Kullanıcı arayüzü, sürekli değişkenler (ateş süresi, tepe sıcaklığı, LDH) için shiny::fluidPage(), shiny::sidebarLayout() ve shiny::sliderInput() ile ikili değişkenler (balgam tıkacıları veya plevral efüzyon, akciğer konsolidasyonu, hipoksemi) için shiny::selectInput() ile oluşturulmuştur. Sunucu mantığı, uyumlu lrm modelinden predict() fonksiyonunu çağırarak bireysel risk tahminlerini hesapladı; ortaya çıkan olasılıklar renderPlot() ve renderText() ile render edildi. Başvuru ShinyApps.io (https://predictrmpp.shinyapps.io/RMPP/) kanalında yayınlandı.