Progettazione dello studio e popolazione di pazienti
Questo studio retrospettivo ha analizzato esami ecografici della tiroide effettuati tra il 13 giugno 2024 e il 13 gennaio 2025. Il protocollo dello studio è stato approvato dal Comitato Etico dell'Ospedale di Amicizia di Pechino, Università di Medicina Capitolare (numero di approvazione BFHHZS20240300) ed è stato condotto in conformità ai principi etici stabiliti nella Dichiarazione di Helsinki; il consenso informato è stato omesso a causa della natura retrospettiva dello studio. Il dataset comprendeva 68 noduli tiroidei provenienti da 63 pazienti (30 benigni e 38 maligni), e il singolo nodulo costituiva l'unità analitica. Le diagnosi di riferimento si basavano sulla citologia ottenuta mediante aspirato con ago sottile guidato da ecografia (FNAC). Il gruppo maligno comprendeva 38 noduli classificati come carcinoma tiroideo papillare alla FNAC.
Criteri di inclusione
I pazienti idonei dovevano soddisfare tutti i seguenti criteri: (1) aver eseguito un'ecografia convenzionale in scala di grigi e Doppler, seguita da CEUS con qualità dell'immagine soddisfacente che permettesse l'analisi successiva del flusso microvascolare; (2) aver ricevuto una diagnosi citopatologica di CPT confermata da biopsia con ago sottile, con o senza test positivo concomitante per la mutazione BRAFV600E; oppure aver ottenuto un risultato citologico di categoria III di Bethesda (atipia di significato indeterminato), ma con test positivo concomitante per la mutazione BRAFV600E; (3) aver ricevuto una diagnosi citopatologica di noduli benigni proliferativi, noduli adenomatosi o noduli follicolari benigni nella categoria II di Bethesda, senza test positivo concomitante per la mutazione BRAFV600E.
Criteri di esclusione
I pazienti sono stati esclusi dallo studio se era presente una delle seguenti condizioni: (1) immagini CEUS di scarsa qualità o segnale insufficiente dei microbolle che impediva una valutazione affidabile del flusso microvascolare; (2) assenza di risultati citologici o patologici; (3) diagnosi istopatologica di un sottotipo raro o speciale di carcinoma tiroideo (ad esempio varianti anaplastiche); oppure (4) reperti citologici suggestivi di neoplasia follicolare (Categoria Bethesda IV) o di qualsiasi lesione di origine follicolare indeterminata, indipendentemente dallo stato mutazionale; (5) tiroidite di Hashimoto
Acquisizione di ultrasuoni, CEUS e SRUS
Tutti gli esami sono stati eseguiti con il sistema ecografico di riferimento e una sonda lineare. I pazienti sono stati posizionati in decubito supino con il collo esteso. Il nodulo bersaglio è stato localizzato e misurato mediante ecografia in modo B; sono state registrate le microcalcificazioni e lo stesso piano di imaging centrato sulla lesione è stato utilizzato per valutare la vascolarizzazione intranodulare con Doppler colorato.
Dopo aver stabilito l'accesso endovenoso, il sistema è stato commutato su una modalità CEUS/URM con basso indice meccanico (URM indica imaging a Ultra-Resolution Microscopy). L'ecografia con mezzo di contrasto (CEUS) è stata eseguita mediante iniezione endovenosa in bolo di 1,2 mL di SonoVue, immediatamente seguita da un flush di 5 ml di soluzione fisiologica; il cronometro sullo schermo e l'acquisizione continua in cine sono stati avviati contestualmente alla somministrazione del bolo. La sonda è stata mantenuta in un piano fisso con pressione minima, e al paziente è stato chiesto di evitare di deglutire durante le fasi di afflusso e di eliminazione del mezzo di contrasto.
Per la CEUS quantitativa, è stata utilizzata una regione di interesse limitata alla lesione per ricavare i parametri tempo-intensità. Per la SRUS, il flusso di lavoro URM ha localizzato e tracciato i segnali delle microbolle dopo il controllo del movimento, generando misurazioni del rapporto vascolare, della complessità, della densità microvascolare, dell'indice di perfusione e della velocità del flusso. Le esportazioni delle immagini rappresentative mostravano le impostazioni di VSP 4, RES 2, CTR 3, SM 2, VEN 3 e CPT 10 s; le impostazioni corrispondenti non erano disponibili per i restanti esami.
Le 25 variabili di input del modello sono elencate nella Tabella 1 e raggruppate per modalità di acquisizione: età e sesso; microcalcificazioni in modo B; vascolarizzazione intranodulare in Doppler color; CEUS qualitativo; CEUS quantitativo; e 11 misurazioni microvascolari SRUS.
Selezione delle caratteristiche
Tutte e 25 le caratteristiche quantitative sono state mantenute nei modelli di apprendimento automatico. Sono state utilizzate solo caratteristiche numeriche; sono stati esclusi i nomi dei pazienti, i numeri di registrazione e i campi relativi alle dimensioni delle lesioni. Non è stata effettuata alcuna selezione delle caratteristiche basata sui dati, e gli stessi predittori predeterminati sono stati inseriti in ciascun classificatore.
Standardizzazione delle caratteristiche
Prima della cross-validation non è stata applicata alcuna trasformazione, imputazione o scala globale. La standardizzazione è stata applicata soltanto al classificatore SVM con funzione radiale mediante una pipeline StandardScaler. Lo scaler è stato addestrato sui noduli di training di ciascun fold e successivamente applicato ai noduli di validazione dello stesso fold. I classificatori basati su alberi hanno ricevuto le scale numeriche originali:

dove x è il valore originale della caratteristica, µ è la media del training-fold e σ è la deviazione standard del training-fold. La pre-elaborazione per fold ha impedito che le osservazioni di validazione contribuissero ai parametri di scala dell'SVM.
Partizionamento dei dati
La valutazione primaria delle prestazioni ha utilizzato una cross-validation stratificata a cinque ripetizioni (StratifiedGroupKFold) con shuffle = True e random_state = 42. Un identificatore del paziente ha definito 63 gruppi, e tutti i noduli dello stesso paziente sono stati assegnati allo stesso fold. Nessun paziente ha contribuito con noduli sia al sottoinsieme di addestramento che a quello di validazione di un fold.
Addestramento e validazione del modello
Protocollo di addestramento
Sono stati valutati cinque classificatori: foresta casuale (100 alberi; random_state = 42), SVM con funzione radiale di base (C = 1,0; gamma = scale; probability = True; pipeline StandardScaler; random_state = 42), albero decisionale (criterio di Gini; profondità illimitata; random_state = 42), XGBoost (100 stimatori; learning_rate = 0,3; max_depth = 6; subsample = 1,0; colsample_bytree = 1,0; random_state = 42) e boosting del gradiente (100 stimatori; learning_rate = 0,1; max_depth = 3; random_state = 42). Non è stata eseguita alcuna ricerca a griglia, ottimizzazione bayesiana, regolazione della soglia o selezione annidata del modello.
Validazione incrociata
All'interno di ciascuna delle cinque suddivisioni raggruppate per paziente, i modelli sono stati addestrati sui gruppi di pazienti rimanenti e valutati sui gruppi esclusi. Per ogni suddivisione sono stati calcolati accuratezza, sensibilità, specificità, precisione, punteggio F1 e area sotto la curva ROC (ROC-AUC), riassunti come media ± deviazione standard. Le predizioni OOF sono state raggruppate considerando tutti e 68 i noduli per generare una curva ROC convalidata tramite cross-validation e una matrice di confusione per ciascun modello.
Valutazione delle prestazioni
Metriche di valutazione
Le prestazioni del modello sono state valutate all'interno di ciascun fold di validazione raggruppato per paziente e a partire dalle predizioni OOF combinate, utilizzando le seguenti metriche:
Accuratezza: Proporzione complessiva di previsioni corrette

Sensibilità (Recupero): Proporzione di noduli maligni effettivi correttamente identificati

Specificità: proporzione di noduli benigni effettivi correttamente identificati

Precisione (Valore Predittivo Positivo): Proporzione di casi predetti come maligni che erano effettivamente maligni

F1-Score: media armonica tra precisione e richiamo

Area sotto la curva ROC (ROC-AUC): Misura della capacità del modello di discriminare tra noduli benigni e maligni a tutti i livelli di classificazione
dove TP = veri positivi (noduli maligni correttamente identificati), TN = veri negativi (noduli benigni correttamente identificati), FP = falsi positivi (noduli benigni classificati erroneamente come maligni) e FN = falsi negativi (noduli maligni classificati erroneamente come benigni).
Analisi della matrice di confusione
Le matrici di confusione OOF sono state generate raggruppando le previsioni effettuate per ciascun nodulo soltanto nel fold in cui il paziente corrispondente era stato escluso. Pertanto, ciascun nodulo ha ricevuto una previsione da un modello addestrato su noduli provenienti da altri pazienti.
Analisi dell'importanza delle caratteristiche
Per il modello Random Forest, i punteggi di importanza delle caratteristiche sono stati calcolati in base alla diminuzione media dell'impurità di Gini su tutti gli alberi decisionali. Sono state identificate e ordinate le 15 caratteristiche più importanti per determinare quali parametri microvascolari hanno contribuito in modo più significativo alle prestazioni di classificazione.
Analisi esplorativa con SHAP
È stata eseguita un'analisi SHAP OOF mirata utilizzando un algoritmo di spiegazione basato su permutazioni per la macchina a vettori di supporto (SVM). Per ciascun nodulo escluso, sono state utilizzate solo le osservazioni del set di addestramento corrispondente come distribuzione di riferimento (128 permutazioni antitetiche; seme casuale = 20260716). I valori SHAP medi assoluti hanno riassunto l'entità dei contributi, mentre i valori con segno hanno indicato la direzione. L'analisi aveva carattere esplorativo e non è stata utilizzata per inferire causalità, identificare biomarcatori indipendenti o definire soglie cliniche.
Analisi statistica
Analisi comparativa
Le prestazioni del modello sono state riassunte in modo descrittivo attraverso le cinque partizioni di validazione raggruppate per paziente e nelle predizioni OOF aggregate. Non è stato effettuato alcun test formale di ipotesi tra modelli né un ordinamento basato su P-value, poiché le partizioni sono correlate e la coorte è di piccole dimensioni.
Per i confronti di gruppo alla linea di base, la normalità delle variabili continue è stata valutata all'interno di ciascun gruppo di outcome utilizzando il test di Shapiro-Wilk. Di Welch t è stato utilizzato quando entrambi i gruppi erano compatibili con la normalità; altrimenti, un test di Mann-Whitney a due code U È stato utilizzato il test t. Le variabili categoriche sono state valutate mediante il test chi-quadrato di Pearson, con il test esatto di Fisher per tabelle 2 x 2 con dati scarsi. I valori di p erano bilaterali, esplorativi e non corretti (alfa = 0,05).
Riproducibilità
La riproducibilità è stata garantita da una definizione specificata del gruppo di pazienti, da un seme casuale fisso (42), da configurazioni fisse del classificatore e da un preelaborazione per fold. Gli identificatori sono stati utilizzati esclusivamente per il raggruppamento e non sono stati inseriti come predittori né esportati con gli output del modello.