Deze studie is goedgekeurd door de Ethische Commissie van het Jinhua Vrouwen- en Kinderziekenhuis (nr. 2024KY099). Geïnformeerde toestemming werd verkregen van de wettelijke voogden van de kinderen.
Patiëntgegevens
Deze retrospectieve studie omvatte kinderen met Mycoplasma pneumoniae pneumonie (MPP) die van november 2022 tot november 2023 werden opgenomen in het Kinderziekenhuis dat is aangesloten op het Jinhua Women and Children's Hospital. MPP werd vastgesteld op basis van ademhalingssymptomen, bevindingen van borstbeeldvorming en positieve Mycoplasma pneumoniae DNA- of RNA-testen. Algemene MPP (GMPP) verwijst naar kinderen die adequaat reageerden op standaard macrolidetherapie, terwijl RMPP werd gedefinieerd als aanhoudende koorts, verergerde symptomen of progressieve beeldvormingsafwijkingen na ten minste 7 dagen standaard macrolidtherapie10. De inclusiecriteria waren leeftijd <14 jaar, bevestigde Mycoplasma pneumoniae-infectie , luchtwegklachten en beeldvormingsbevindingen die overeenkwamen met longontsteking. De exclusiecriteria waren gemengde infectie binnen 10 dagen na het begin van de ziekte, leukemie, chronische longziekte, immunodeficiëntie, eerdere immunosuppressieve therapie, opname tijdens de herstelfase of onvolledige dossiers. De standaardtherapie verwees azitromycine 10 mg/kg eenmaal daags, oraal of intraveneus, met een maximale dosis van 500 mg/dag. Kinderen die niet aan de RMPP-criteria voldeden, werden geclassificeerd als algemene MPP (GMPP). In totaal werden 500 in aanmerking komende kinderen willekeurig verdeeld in een trainingscohort (n = 375) en een validatiecohort (n = 125) met R-softwareversie 4.1.2 (function: sample(), set.seed = 42). De validatiecohort was exclusief gereserveerd voor externe modelevaluatie en werd in geen enkele stap van modelconstructie gebruikt. Elektronische medische dossiergegevens werden geëxtraheerd met behulp van een gestandaardiseerde extractietemplate die a priori was ontwikkeld. Gegevensvelden omvatten demografie, klinische presentatie, laboratoriumwaarden, beeldvormende bevindingen en behandelgegevens. Twee getrainde onderzoekers haalden onafhankelijk de gegevens uit en de afwijkingen werden door consensus opgelost. De bevindingen van de borst-CT werden beoordeeld door radiologen die blind waren voor klinische groepering en laboratoriumresultaten. De betrouwbaarheid van CT-beeldinterpretatie tussen beoordelaars werd beoordeeld met behulp van Cohen's kappa-statistiek voor binaire beeldvormingsbevindingen (sputumproppen of pleurale effusie: ja/nee; longconsolidatie: ja/nee). Ten slotte werden 129 kinderen geclassificeerd als RMPP en 371 als algemene MPP.
Variabelen
Demografische, klinische, laboratorium- en beeldvormende variabelen werden uit elektronische medische dossiers gehaald met behulp van een vooraf gedefinieerd formulier. Klinische variabelen waren onder andere leeftijd, geslacht, koortsduur vóór opname, pieklichaamstemperatuur en hypoxemie. De koortsduur werd gedefinieerd als het interval van het begin van koorts tot opname. De piek van de lichaamstemperatuur was de hoogste geregistreerde temperatuur vóór of binnen 24 uur na opname. Hypoxemie werd gedefinieerd als perifere arteriële zuurstofsaturatie <92% op kamerlucht of de behoefte aan extra zuurstof. Laboratoriumvariabelen waren onder andere WBC, HB, PLT, CRP, ALB, ALT, CK-MB, LDH, D-dimeer, IL-6, IL-8, IL-10, IL-17, PCT en NE%. Nuchtere veneuze bloedmonsters werden binnen 24 uur na opname afgenomen, en de bemonsteringstijd ten opzichte van koortsstart en het starten van antibiotica werd gemeten wanneer beschikbaar. Mycoplasma pneumoniae DNA- of RNA-testen werden gebruikt voor etiologische bevestiging. De keeluitstrijkje realtime kwantitatieve PCR-test voor Mycoplasma pneumoniae DNA of RNA was positief. Kwantitatieve DNA/RNA-belasting werd niet opgenomen omdat gestandaardiseerde belastinggegevens niet beschikbaar waren voor alle patiënten. Een CT-scan van de borst die binnen 3 dagen voor of na opname werd uitgevoerd, werd beoordeeld op sputumproppen, pleuravocht en longconsolidatie. "Sputumproppen of pleurale effusie" werd als positief geregistreerd wanneer een van beide bevindingen aanwezig was. Longconsolidatie werd gedefinieerd als segmentale of lobaire parenchymale opaciteit op CT. Er werd een CT-scanner van 256 slices gebruikt. Kinderen werden in rugligging geplaatst en laag-dosis axiale borstscanning werd uitgevoerd met een detector van geschikte breedte onder ademhaling of terwijl de kinderen sliepen. De scandekking strekte zich uit van de longtop tot aan de longbasis om het gehele longparenchym te bedekken. De scanparameters werden als volgt ingesteld: buisspanning 100 kVp, slicedikte 5 mm, reconstructie slice dikte 1,25 mm, matrix 512 × 512, portaalrotatietijd 0,28 s, ruisindex 12. Afbeeldingen werden gereconstrueerd met behulp van ASIR-V en DLIR-algoritmen. Elke groep bevatte gereconstrueerde beelden met ASIR-V gewichten van 20%, 50% en 80%, evenals DLIR-L, DLIR-M en DLIR-H beelden. De verkregen beeldgegevens werden geïmporteerd in het werkstation. Alle beelden werden onafhankelijk dubbelblind beoordeeld door twee radiologen met meer dan vijf jaar klinische ervaring. Ze beoordeelden de aanwezigheid van abnormale longsignalen en somden de belangrijkste beeldvormende kenmerken samen. Voor gevallen met uiteenlopende interpretaties werd door wederzijdse discussie consensus bereikt.
Statistische analyse
Statistische analyses werden uitgevoerd met R-software versie 4.1.2. Categorische variabelen worden gepresenteerd als n (%) en vergeleken met behulp van de chi-kwadraattest of de exacte test van Fisher. Continue variabelen worden uitgedrukt als gemiddelde ± standaarddeviatie of mediaan (interkwartielbereik), volgens de verdeling. De normaliteit werd beoordeeld met behulp van de Shapiro-Wilk-test. Normaal verstreken variabelen werden vergeleken met de onafhankelijke steekproeven t-test, en niet-normaal verstreken variabelen met de Mann-Whitney U-test. Voor de modelbouw werd de basisvergelijkbaarheid tussen de trainings- en validatiecohorten bevestigd met behulp van χ2 - of Mann-Whitney U-tests voor alle kandidaatvoorspellers, waarbij geen significante verschillen tussen cohorten werden aangetoond (alle P > 0,05). In de trainingscohort werden variabelen met P < 0,05 in univariate logistische regressie ingevoerd in multivariabele logistische regressie. Multicollineariteit werd beoordeeld met behulp van variantie-inflatiefactoren (VIF's) die werden berekend met het "car"-pakket; alle geselecteerde voorspellers hadden VIF < 5, wat wijst op geen substantiële multicollineariteit. De lineariteitsaanname voor LDH werd onderzocht met behulp van logaritmische transformatie en beperkte kubieke splines (met het "rms"-pakket met 3 knopen geplaatst op het 10e, 50e en 90e percentiel); de lineaire vorm bleef behouden wanneer de modelfit niet verbeterde (likelihood ratio-test: P > 0,05 voor de niet-lineaire splineterm, ter ondersteuning van de lineaire specificatie). Een nomogram werd geconstrueerd uit het uiteindelijke multivariabele model met behulp van het "rms"-pakket (versie 6.3-0; functies: lrm voor logistische regressie, Predict voor waardevoorspelling en nomogram voor grafische weergave). Discriminatie werd geëvalueerd met behulp van ontvanger-bedrijfskarakteristiekcurves, oppervlakte onder de curve en C-index met 95% betrouwbaarheidsintervallen ("pROC"-pakket). De kalibratie werd beoordeeld met kalibratiecurves en 1.000 bootstrap-resamples. ("RMS"-pakket: Calibrate Function). De optimale voorspellingsgrens werd onafhankelijk binnen elke cohort bepaald met behulp van Youden's index via het pROC-pakket: 0,222 in de trainingscohort en 0,247 in de validatiecohort. De trainingsafgeleide cutoff (0,222) werd ook toegepast op de validatiecohort voor cross-cohort prestatievergelijking, in overeenstemming met standaard interne validatiepraktijken. Het 95% betrouwbaarheidsinterval voor de voorspelde kans werd bijvoorbeeld berekend met modelgebaseerde standaardfouten op de logitschaal. Beslissingscurveanalyse werd gebruikt om het netto voordeel over drempelkansen te evalueren ("RMDA"-pakket: decision_curve en plot_decision_curve functies). Een tweezijdige P < 0,05 werd als statistisch significant beschouwd. De volledige R-scripts, inclusief datapreprocessing, modelconstructie, validatie en nomogramgeneratie, zijn op verzoek beschikbaar bij de corresponderende auteur.
De ontbrekende gegevens waren minimaal (alle variabelen ontbraken <2%) en werden behandeld met mediane imputatie voor continue variabelen en modusimputatie voor categorische variabelen vóór de analyse. Alle categorische voorspellers werden gecodeerd als binaire (0/1) indicatorvariabelen. Continue voorspellers (koortsduur, piektemperatuur, LDH) werden in hun oorspronkelijke klinische eenheden behouden zonder categorisatie. Er werd geen variabeleselectie op basis van univariate screening uitgevoerd vóór multivariabele modellering; in plaats daarvan werden alle klinisch relevante kandidaten die in de literatuur waren geïdentificeerd overwogen, en degenen met P < 0,05 in de univariate analyse van de trainingscohort werden doorgeschakeld naar het multivariabele model.
Webcalculator-implementatie
Het webgebaseerde dynamische nomogram werd uitgerold met het "shiny" pakket (versie 1.7.4) in R. De gebruikersinterface is gebouwd met shiny::fluidPage(), shiny::sidebarLayout() en shiny::sliderInput() voor continue variabelen (koortsduur, piektemperatuur, LDH) en shiny::selectInput() voor binaire variabelen (sputumproppen of pleurale effusie, longconsolidatie, hypoxemie). De serverlogica riep de predict()-functie aan uit het aangepaste lrm-model om individuele risicoschattingen te berekenen, waarbij de resulterende kansen werden weergegeven via renderPlot() en renderText(). De applicatie werd gehost op ShinyApps.io (https://predictrmpp.shinyapps.io/RMPP/).