Çalışma tasarımı ve hasta popülasyonu
Bu retrospektif çalışma, 13 Haziran 2024 ile 13 Ocak 2025 tarihleri arasında elde edilen tiroid ultrasonografisi incelemelerini analiz etmiştir. Çalışma protokolü, Capital Medical University, Beijing Friendship Hospital Etik Kurulu tarafından onaylanmış (onay No. BFHHZS20240300) ve Helsinki Deklarasyonu'nda belirtilen etik ilkeler doğrultusunda yürütülmüştür; retrospektif tasarım nedeniyle bilgilendirilmiş onamdan feragat edilmiştir. Veri seti, 63 hastadan alınan 68 tiroid nodülünden (30 benign ve 38 malign) oluşmuş olup, analiz birimi olarak nodüller esas alınmıştır. Referans teşhisler, ultrason eşliğinde ince iğne aspirasyon sitolojisine (İİAS) dayandırılmıştır. Malign grup, İİAS sonucunda papiller tiroid karsinomu olarak raporlanan 38 nodülden oluşmaktadır.
Dahil edilme kriterleri
Uygun hastaların aşağıdaki kriterlerin tamamını karşılaması gerekmiştir: (1) konvansiyonel gri tonlamalı ve Doppler ultrasonografiye tabi tutulmuş, ardından mikrodamar akım görüntüleme analizine izin veren yeterli görüntü kalitesine sahip CEUS uygulanmış; (2) eş zamanlı pozitif BRAFV600E mutasyon testi ile birlikte veya bu test olmaksızın ince iğne aspirasyon biyopsisi ile doğrulanmış PTC sitopatolojik tanısına sahip olanlar veya Bethesda Kategori III (belirsiz öneme sahip atipi) sitoloji sonucuna sahip olup aynı zamanda pozitif BRAFV600E mutasyon testi saptananlar; (3) eş zamanlı pozitif BRAFV600E mutasyon testi olmaksızın Bethesda Kategori II'de benign, proliferatif nodüller, adenomatöz nodüller veya benign foliküler nodül sitopatolojik tanısına sahip olanlar.
Hariç tutma kriterleri
Aşağıdaki koşullardan herhangi birinin mevcut olması durumunda hastalar çalışmadan hariç tutulmuştur: (1) güvenilir mikrovasküler akış değerlendirmesini engelleyen düşük kaliteli CEUS görüntüleri veya yetersiz mikrobaloncuk sinyali; (2) sitolojik veya patolojik sonuçların olmaması; (3) nadir veya özel bir tiroid karsinomu alt tipi (örneğin anaplastik varyantlar) histopatolojik tanısı; (4) mutasyon durumuna bakılmaksızın, foliküler neoplazm (Bethesda Kategori IV) düşündüren sitolojik bulgular veya belirsiz foliküler kökenli herhangi bir lezyon; (5) Hashimoto tiroiditi
Ultrasonda, CEUS ve SRUS görüntüleme
Tüm incelemeler, referans verilen ultrason sistemi ve lineer dizilimli bir prob ile gerçekleştirilmiştir. Hastalar, boyun ekstate edilmiş şekilde sırtüstü pozisyonda konumlandırılmıştır. Hedef nodül B-mod ultrason ile lokalize edilmiş ve ölçülmüştür; mikrokalsifikasyonlar kaydedilmiş ve intranodüler vaskülariteyi Color Doppler ile değerlendirmek için aynı lezyon merkezli görüntüleme düzlemi kullanılmıştır.
İntravenöz erişim sağlandıktan sonra, sistem düşük mekanik indeksli bir CEUS/URM moduna (URM, Ultra-Çözünürlüklü Mikroskopi görüntüleme anlamına gelir) geçirildi. CEUS, 1,2 mL SonoVue intravenöz bolus enjeksiyonunun ardından hemen 5 mL salin yıkaması yapılarak gerçekleştirildi; bolus uygulamasıyla birlikte ekrandaki zamanlayıcı ve sürekli cine görüntüleme başlatıldı. Prob, minimum basınçla sabit bir düzlemde tutuldu ve hastadan wash-in ve wash-out süreçleri sırasında yutkunmaktan kaçınması istendi.
Kantitatif CEUS için, zaman-yoğunluk parametrelerini türetmek amacıyla lezyonla sınırlı bir ilgi alanı kullanılmıştır. SRUS için, URM iş akışı hareket kontrolünden sonra mikrobaloncuk sinyallerini lokalize edip takip etmiş ve damar oranı, karmaşıklık, mikrovasküler yoğunluk, perfüzyon indeksi ve akış hızı ölçümlerini oluşturmuştur. Temsili görüntü çıktıları VSP 4, RES 2, CTR 3, SM 2, VEN 3 ve CPT 10 s ayarlarını göstermekteydi; geri kalan incelemeler için karşılık gelen ayarlar mevcut değildi.
25 model girdi değişkeni Tablo 1'de listelenmiş ve elde etme modalitesine göre şu şekilde gruplandırılmıştır: yaş ve cinsiyet; B-mod mikrokalsifikasyon; Renkli Doppler intranodüler vaskülarite; kalitatif CEUS; kantitatif CEUS ve 11 SRUS mikrovasküler ölçümü.
Öznitelik seçimi
25 kantitatif özelliğin tamamı makine öğrenimi modellerinde tutulmuştur. Yalnızca sayısal özellikler kullanılmış; hasta isimleri, kayıt numaraları ve lezyon boyutu alanları hariç tutulmuştur. Veriye dayalı herhangi bir özellik seçimi yapılmamış ve önceden belirlenmiş aynı öngörücüler her bir sınıflandırıcıya girilmiştir.
Öznitelik standardizasyonu
Çapraz doğrulamadan önce herhangi bir dönüşüm, imputasyon veya küresel ölçekleme uygulanmamıştır. Standardizasyon, yalnızca StandardScaler hattı aracılığıyla radyal tabanlı fonksiyon SVM'ye uygulanmıştır. Ölçekleyici, her katlamanın eğitim modülleri üzerinde eğitilmiş ve ardından o katlamanın doğrulama modüllerine uygulanmıştır. Ağaç tabanlı sınıflandırıcılar orijinal sayısal ölçekleri kullanmıştır:

burada x orijinal öznitelik değerini, µ eğitim katlaması ortalamasını ve σ eğitim katlaması standart sapmasını temsil eder. Katlama bazlı ön işleme, doğrulama gözlemlerinin SVM ölçekleme parametrelerine katkıda bulunmasını önlemiştir.
Veri bölümlendirme
Birincil performans değerlendirmesinde, shuffle = True ve random_state = 42 değerleri ile beş katmanlı StratifiedGroupKFold çapraz doğrulaması kullanılmıştır. Bir hasta tanımlayıcısı 63 grubu tanımlamış ve aynı hastadan alınan tüm nodüller aynı katmana atanmıştır. Hiçbir hasta, bir katmanın hem eğitim hem de doğrulama alt kümelerine nodül katkısında bulunmamıştır.
Model eğitimi ve doğrulaması
Eğitim protokolü
Beş sınıflandırıcı değerlendirildi: rastgele orman (100 ağaç; random_state = 42), radyal tabanlı fonksiyon SVM (C = 1.0; gamma = scale; probability = True; StandardScaler hattı; random_state = 42), karar ağacı (Gini kriteri; kısıtlanmamış derinlik; random_state = 42), XGBoost (100 tahminci; learning_rate = 0.3; max_depth = 6; subsample = 1.0; colsample_bytree = 1.0; random_state = 42) ve gradyan artırma (100 tahminci; learning_rate = 0.1; max_depth = 3; random_state = 42). Izgara araması, Bayesian optimizasyonu, eşik ayarlama veya iç içe model seçimi yapılmamıştır.
Çapraz doğrulama
Hasta bazlı olarak gruplandırılmış beş katlı çapraz doğrulamanın her birinde, modeller geri kalan hasta grupları üzerinde eğitilmiş ve ayrılan gruplar üzerinde değerlendirilmiştir. Her kat için doğruluk, duyarlılık, özgüllük, kesinlik, F1-skoru ve ROC-AUC hesaplanmış ve ortalama ± SD olarak özetlenmiştir. Her model için bir adet çapraz doğrulanmış ROC eğrisi ve bir adet hata matrisi oluşturmak amacıyla, OOF tahminleri toplam 68 nodül üzerinden birleştirilmiştir.
Performans değerlendirmesi
Değerlendirme metrikleri
Model performansı, her bir hasta gruplu doğrulama katlaması içinde ve havuzlanmış OOF tahminlerinden şu metrikler kullanılarak değerlendirilmiştir:
Doğruluk: Doğru tahminlerin genel oranı

Duyarlılık (Geri Çağırma): Gerçek malign nodüllerin doğru şekilde tanımlanma oranı

Özgüllük: Doğru şekilde tanımlanan gerçek benign nodüllerin oranı

Hassasiyet (Pozitif Öngörü Değeri): Öngörülen malign vakaların gerçekten malign olanlarının oranı

F1-Skoru: Kesinlik ve duyarlılığın harmonik ortalaması

Alıcı İşletim Karakteristik Eğrisi Altındaki Alan (ROC-AUC): Modelin, tüm sınıflandırma eşikleri boyunca benign ve malign nodülleri birbirinden ayırt edebilme yeteneğinin ölçümü
burada TP = gerçek pozitifler (doğru tanımlanmış malign nodüller), TN = gerçek negatifler (doğru tanımlanmış benign nodüller), FP = yanlış pozitifler (yanlışlıkla malign olarak sınıflandırılan benign nodüller) ve FN = yanlış negatifler (yanlışlıkla benign olarak sınıflandırılan malign nodüller) anlamına gelir.
Karışıklık matrisi analizi
OOF konfüzyon matrisleri, her nodül için yapılan tahminlerin yalnızca o hastanın dışarıda tutulduğu katmanda (fold) birleştirilmesiyle oluşturulmuştur. Böylece, her nodül için tahmin, diğer hastalardan alınan nodüllerle eğitilmiş bir model tarafından yapılmıştır.
Özellik önem analizi
Random Forest modeli için özellik önem skorları, tüm karar ağaçlarındaki ortalama Gini safsızlık azalmasına dayanarak hesaplanmıştır. Sınıflandırma performansına en önemli katkıyı sağlayan mikrovasküler parametreleri belirlemek amacıyla en önemli 15 özellik tanımlanmış ve sıralanmıştır.
Keşifsel SHAP analizi
SVM için permütasyon tabanlı bir açıklayıcı kullanılarak odaklanmış bir OOF SHAP analizi gerçekleştirilmiştir. Ayrılan her bir nodül için, arka plan dağılımı olarak yalnızca ilgili eğitim katlaması gözlemleri kullanılmıştır (128 antitetik permütasyon; rastgele tohum = 20260716). Ortalama mutlak SHAP değerleri katkıların büyüklüğünü özetlemiş ve işaretli değerler yönü belirtmiştir. Analiz keşifsel niteliktedir ve nedensellik çıkarımı yapmak, bağımsız biyobelirteçler tanımlamak veya klinik eşik değerler belirlemek amacıyla kullanılmamıştır.
İstatistiksel analiz
Karşılaştırmalı analiz
Model performansı, hasta gruplu beş doğrulama katlaması genelinde tanımlayıcı olarak ve havuzlanmış OOF tahminlerinde özetlenmiştir. Katlamaların birbiriyle ilişkili olması ve kohortun küçük olması nedeniyle, modeller arasında formal bir hipotez testi veya P-değeri tabanlı bir sıralama yapılmamıştır.
Başlangıç grubu karşılaştırmaları için, sürekli değişkenlerin normalliği her bir sonuç grubu içinde Shapiro-Wilk testi kullanılarak değerlendirildi. Her iki grup da normallik ile uyumlu olduğunda Welch t testi kullanıldı; aksi takdirde, iki yönlü Mann-Whitney U testi uygulandı. Kategorik değişkenler Pearson ki-kare testi ile, seyrek 2 x 2 tablolar için ise Fisher'ın kesin testi ile değerlendirildi. P değerleri iki yönlü, keşifsel ve düzeltilmemişti (alfa = 0,05).
Tekrarlanabilirlik
Tekrarlanabilirlik; belirlenmiş bir hasta grubu tanımı, sabit rastgele tohum (42), sabit sınıflandırıcı konfigürasyonları ve katman bazlı ön işleme ile desteklenmiştir. Tanımlayıcılar yalnızca gruplandırma için kullanılmış olup, öngördürücü olarak girilmemiş veya model çıktılarıyla birlikte dışa aktarılmamıştır.