Zaproponowaną strukturę XAI zaimplementowano z wykorzystaniem zbioru danych Pima Indians Diabetes Dataset jako reprezentatywnego zbioru danych medycznych. Zbiór Pima Indians Diabetes Dataset posłużył jako jedyny zbiór danych do walidacji eksperymentalnej. Wszystkie raportowane wyniki predykcyjne, wyjaśnialności, statystyczne oraz dotyczące reprodukowalności zostały wygenerowane na podstawie tego zbioru. Zbiory TCGA-BRCA, TCGA-UCEC, MIMIC-III, ISIC 2019, HAM100, OhioT1DM i BraTS 2021 nie zostały poddane ocenie eksperymentalnej i zostały uwzględnione jedynie jako przykłady ilustrujące możliwość zastosowania ogólnego protokołu w różnych modalnościach danych medycznych. Kompletny zbiór danych został wykorzystany po usunięciu nieprawidłowych i zduplikowanych rekordów oraz przeprowadzeniu określonych operacji przetwarzania wstępnego przed opracowaniem modelu. Zbiór danych podzielono na niezależne podzbiory treningowe, walidacyjne i testowe, stosując zdefiniowaną strategię warstwowego podziału (stratified splitting). Zachowano niezależność próbek pomiędzy trzema podzbiorami, aby zminimalizować ryzyko wycieku danych.
Model predykcyjny skonfigurowano przy użyciu zdefiniowanej wcześniej architektury i hiperparametrów. Model wytrenowano za pomocą optymalizatora Adam z predefiniowanym współczynnikiem uczenia i rozmiarem partii przez maksymalnie 10 epok. Zastosowano mechanizm wczesnego zatrzymania (early stopping) w oparciu o stratę walidacyjną, a zachowano model odpowiadający najlepszym wynikom walidacji. W celu zwiększenia powtarzalności określono ziarna liczb pseudolosowych (random seeds) dla podziału zbioru danych, inicjalizacji modelu oraz powtórnych eksperymentów.
Wytrenowany model oceniono na niezależnym zbiorze testowym, wykorzystując wskaźniki dokładności, precyzji, czułości, swoistości, wynik F1, współczynnik korelacji Matthewsa (MCC), pole pod krzywą charakterystyki pracy odbiornika (AUC-ROC) oraz średnią precyzję (AP). Zaproponowany model porównano z wcześniej zdefiniowanymi modelami bazowymi, stosując identyczne procedury wstępnego przetwarzania, podziały danych oraz protokoły oceny. Na podstawie predykcji ze zbioru testowego wygenerowano krzywe charakterystyki pracy odbiornika (ROC), krzywe precyzji-czułości oraz macierz pomyłek.
W celu oceny stabilności wydajności przeprowadzono pięciokrotną walidację krzyżową na danych treningowych. Oszacowano 95% przedziały ufności dla głównych wskaźników wydajności oraz przeprowadzono zdefiniowane wcześniej porównania statystyczne między proponowanym modelem a modelami bazowymi.
Piokrotna walidacja krzyżowa pozwoliła uzyskać dokładność na poziomie 93,01 ± 0,48%, AUC-ROC wynoszące 0,954 ± 0,07, precyzję 92,42 ± 0,87%, czułość 93,02 ± 0,45% oraz wynik F1 na poziomie 92,72 ± 0,62%. 95% przedziały ufności bootstrap obliczone z 1 0 próbek wyniosły 0,897–0,973 dla dokładności, 0,890–0,970 dla precyzji, 0,80–0,963 dla czułości, 0,87–0,965 dla wyniku F1 oraz 0,931–0,984 dla AUC-ROC. Porównania statystyczne z modelami bazowymi CNN, Random Forest i SVM przeprowadzono przy użyciu testu McNemara dla dokładności, testu DeLonga dla AUC-ROC oraz parzystego testu bootstrap z 1 0 próbkami dla wyniku F1.
Pełna procedura szkolenia i ewaluacji została powtórzona w 10 niezależnych uruchomieniach z wykorzystaniem różnych ziarn losowości. Powtórzone uruchomienia dały wynik AUC-ROC na poziomie 0,957 ± 0,08, dokładność 93,24 ± 0,74% oraz wynik F1-score wynoszący 92,35 ± 0,92%, co wskazuje na stosunkowo niską zmienność między poszczególnymi wykonaniami. Metryki wydajności uzyskane w powtórzonych uruchomieniach zostały podsumowane za pomocą średniej i odchylenia standardowego. Zmienność między uruchomieniami została zbadana, aby określić, czy wydajność predykcyjna pozostawała stabilna przy wielokrotnym wykonaniu.
W niniejszym przykładzie praktycznym nie przeprowadzono walidacji zewnętrznej, ponieważ nie wykorzystano niezależnego zewnętrznego zbioru danych. W związku z tym wszystkie raportowane wyniki predykcyjne, statystyczne oraz wyniki powtarzalności uzyskano z Pima Indians Diabetes Dataset, korzystając z wcześniej zdefiniowanych partycji treningowych, walidacyjnych i niezależnych zbiorów testowych. Walidacja zewnętrzna została zachowana w protokole jako opcjonalna procedura dla zastosowań, w których dostępny jest niezależny zbiór danych pochodzący z innej instytucji, populacji, regionu geograficznego lub okresu czasu.
Wyjaśnienia modelu zostały wygenerowane przy użyciu technik wyjaśnialności stosowalnych do tabelarycznego zbioru danych Pima Indians Diabetes Dataset, w tym SHAP oraz LIME. Przeprowadzono ocenę globalnej istotności cech, a lokalne wyjaśnienia specyficzne dla pacjentów wygenerowano przy użyciu wydzielonych próbek testowych. Wyniki wyjaśnień zapisano wraz z odpowiadającymi im przewidywaniami modelu i wartościami cech, aby ułatwić powtarzalność oraz późniejszą interpretację.
Dla zachowania przejrzystości, niniejszy przykład praktyczny obejmował dziewięć głównych etapów procesu określonych w Tabeli 1. Walidacja zewnętrzna oraz ocena ekspercka kliniczna zostały zachowane jako opcjonalne moduły walidacyjne ogólnego protokołu. Walidacja zewnętrzna nie została przeprowadzona, ponieważ nie wykorzystano niezależnego zewnętrznego zbioru danych, a ocena ekspercka kliniczna nie została wykonana, ponieważ w niniejszym przykładzie praktycznym nie uwzględniono formalnego panelu oceny eksperckiej. W związku z tym moduły opcjonalne te nie są uważane za część dziewięciostopniowego procesu eksperymentalnego i nie są raportowane jako wyniki eksperymentalne.
Procedury wstępnego przetwarzania i partycjonowania zbioru danych pozwoliły na utworzenie ustandaryzowanego zbioru danych odpowiedniego do opracowania modelu. Usunięto duplikaty i niespójne rekordy, braki danych uzupełniono zgodnie z przyjętą strategią, cechy numeryczne ustandaryzowano, a następnie zbiór danych podzielono na niezależne podzbiory treningowe, walidacyjne i testowe. Charakterystyka otrzymanego zbioru danych oraz procedury wstępnego przetwarzania zostały zestawione w Tabeli 2. Ogólny schemat przygotowania i wstępnego przetwarzania zbioru danych z zakresu opieki zdrowotnej przedstawiono na Rysunku 2, natomiast schemat przygotowania eksperymentalnego, wstępnego przetwarzania, partycjonowania i oceny jakości zastosowany konkretnie do zbioru danych Pima Indians Diabetes Dataset pokazano na Rysunku 3. Końcowe środowisko obliczeniowe, architektura modelu oraz konfiguracja hiperparametrów wykorzystane do uzyskania raportowanych wyników zostały podsumowane w Tabeli 3.
Wykonanie sekcji 3 i 4 pozwoliło na uzyskanie zestandaryzowanego zbioru danych opieki zdrowotnej odpowiedniego do opracowania modelu. Procedury preprocessingu usunęły zdublowane i niespójne rekordy, uzupełniły brakujące wartości, zestandaryzowały cechy numeryczne, zakodowały zmienne kategoryczne tam, gdzie było to możliwe, oraz podzieliły zbiór danych na podzbiory treningowe, walidacyjne i testowe. Wynikowe dane zapewniły zestandaryzowane dane wejściowe wymagane do późniejszego opracowania modelu.
Po zakończeniu sekcji 5 i 6 skonfigurowany model wykazał stabilną zbieżność podczas treningu. Krzywe uczenia wykazały jednoczesny spadek strat treningowych i walidacyjnych oraz wzrost dokładności treningowej i walidacyjnej. Optymalizacja hiperparametrów poprawiła generalizację modelu, bez dowodów na istotne przeuczenie. W celu zapewnienia powtarzalności zoptymalizowany model został zarchiwizowany wraz z ostateczną architekturą, ustawieniami hiperparametrów, wersjami oprogramowania, ziarnami liczb pseudolosowych oraz wagami wytrenowanego modelu. Specyfikacje treningu modelu i wyniki optymalizacji podsumowano w Tabeli 4, a odpowiadające im krzywe uczenia dla zbioru treningowego i walidacyjnego przedstawiono na Rysunku 4.
W sekcji 7 oceniono zdolność predykcyjną modelu przy użyciu zestandaryzowanych metryk wydajności. Oceńne metryki klasyfikacji obejmowały dokładność (accuracy), precyzję (precision), czułość (recall), swoistość (specificity), wynik F1 (F1-score), MCC, AUC-ROC oraz AP. Zaproponowany model porównano z wcześniej zdefiniowanymi modelami bazowymi, stosując te same podziały zbioru danych, procedury wstępnego przetwarzania oraz protokół oceny. Porównanie wydajności predykcyjnej przedstawiono w Tabeli 5, natomiast krzywe ROC, krzywe precyzja-czułość, macierz pomyłek oraz porównawcze metryki wydajności zaprezentowano na Ryc. 5.
Zgodnie z Sekcją 8, w celu oceny interpretowalności modelu wygenerowano zarówno globalne, jak i lokalne wyjaśnienia jego predykcji. Wyjaśnienia modelu zostały opracowane przy użyciu metod SHAP i LIME dla tabelarycznego zbioru danych Pima Indians Diabetes Dataset, a w celu zilustrowania zmiany w predykcji przygotowano kontrfaktyczne wyjaśnienie specyficzne dla pacjenta. Metoda SHAP została wykorzystana do wygenerowania globalnej istotności cech oraz wizualizacji w postaci wykresów kaskadowych (waterfall) i zależności cech dla konkretnych pacjentów, natomiast LIME posłużył do generowania lokalnych wyjaśnień na podstawie próbek z wydzielonego zbioru testowego. Odpowiednie wyniki analizy wyjaśnialności przedstawiono na Rysunku 6.
Ostatni etap protokołu służył ocenie wiarygodności statystycznej i powtarzalności przepływu pracy. Cały proces powtórzono w 10 niezależnych uruchomieniach, stosując różne ziarna losowe (random seeds), przy zachowaniu tej samej strategii partycjonowania zbioru danych, parametrów przetwarzania wstępnego, architektury modelu, hiperparametrów, środowiska programistycznego oraz procedur ewaluacyjnych. Powtórzone uruchomienia dały AUC-ROC na poziomie 0,957 ± 0,08, dokładność (accuracy) 93,24 ± 0,74% oraz wskaźnik F1 na poziomie 92,35 ± 0,92%, co wskazuje na stosunkowo niską zmienność między kolejnymi wykonaniami.
Stabilność wyjaśnień nie została ilościowo oceniona w ramach niniejszej walidacji. Chociaż dla zbioru danych Pima Indians Diabetes Dataset wygenerowano wyjaśnienia SHAP i LIME, nie przeprowadzono oddzielnej analizy korelacji rang między seriami ani analizy pokrycia cech. W związku z tym nie raportuje się żadnych numerycznych wskaźników stabilności wyjaśnień ani zgodności atrybucji. Ilościowa ocena stabilności wyjaśnień została zachowana w protokole ogólnym jako opcjonalna procedura powtarzalności dla zastosowań, w których dostępne są powtarzalne wyniki wyjaśnień.
Porównania statystyczne oceniano przy użyciu zdefiniowanego progu istotności p < 0.05. W stosownych przypadkach stosowano dwustronne testy statystyczne, a odpowiadające im statystyki testowe, wartości p oraz 95% przedziały ufności raportowano w celu ilościowego określenia istotności statystycznej i niepewności zaobserwowanych różnic w wydajności. Wyniki eksperymentalnej walidacji statystycznej i powtarzalności, w tym wydajność walidacji krzyżowej, przedziały ufności, porównania statystyczne oraz zmienność między powtórzeniami, podsumowano w Tabeli 6. Odpowiadające im testy statystyczne i analizy powtarzalności zilustrowano na Rysunku 7.
Wyniki te dostarczyły dowodów eksperymentalnych na stabilność predykcyjną i powtarzalność w testowanych warunkach obliczeniowych oraz dla zastosowanego zbioru danych. Środowisko obliczeniowe, charakterystyka zbioru danych, procedury wstępnego przetwarzania, konfiguracja modelu, analizy statystyczne oraz ustawienia wyjaśnialności niezbędne do niezależnej replikacji zostały udokumentowane zgodnie z listą kontrolną powtarzalności przedstawioną w Tabeli 7. W zaprezentowanym w niniejszej pracy przykładzie walidacyjnym nie przeprowadzono oceny wygenerowanych wyników XAI przez ekspertów klinicznych.
Ogólnie zastosowanie protokołu pozwoliło na stworzenie zestandaryzowanego zbioru danych medycznych odpowiedniego do opracowania modelu AI oraz zoptymalizowanego modelu z wykorzystaniem zdefiniowanych ustawień hiperparametrów. Przeprowadzony proces umożliwił systematyczną ocenę wydajności predykcyjnej, generowanie wyjaśnień modelu przy użyciu odpowiednich technik XAI oraz statystyczną ocenę odporności i powtarzalności modelu. Wspólnie wyniki te wykazały możliwość implementacji i powtarzalność proponowanego procesu XAI w warunkach eksperymentalnych ocenianych w opracowanym przykładzie walidacyjnym.

Rysunek 1: Dziewięcioetapowy podstawowy schemat pracy w celu opracowania reprodukowalnych i wyjaśnialnych modeli AI w ochronie zdrowia. Schemat obejmuje (1) definicję zadania predykcyjnego w ochronie zdrowia, (2) konfigurację środowiska obliczeniowego, (3) pozyskanie i walidację zbioru danych, (4) wstępne przetwarzanie danych, (5) partycjonowanie zbioru danych, (6) trenowanie modelu predykcyjnego, (7) ocenę wydajności predykcyjnej, (8) analizę wyjaśnialności oraz (9) walidację statystyczną i ocenę reprodukowalności. Walidacja zewnętrzna oraz ocena ekspercka kliniczna są traktowane jako opcjonalne rozszerzenia protokołu i nie są włączone do dziewięcioetapowego podstawowego schematu pracy. Kliknij tutaj, aby zobaczyć powiększoną wersję tego rysunku.

Rysunek 2: Przebieg procesu przygotowania i wstępnego przetwarzania zbioru danych z zakresu opieki zdrowotnej. (A) Pozyskiwanie danych medycznych z dokumentacji klinicznej, obrazowania medycznego, sygnałów fizjologicznych oraz multimodalnych źródeł danych. (B) Integracja i wstępne przetwarzanie danych, w tym obsługa brakujących wartości, normalizacja, usuwanie szumów i augmentacja. (C) Podział zbioru danych na podzbiory treningowy, walidacyjny i testowy. (D) Ocena jakości przedstawiająca rozkład klas, normalizację cech oraz wyniki wstępnego przetwarzania przed opracowaniem modelu. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 3: Przebieg eksperymentalny przygotowania, wstępnego przetwarzania, partycjonowania i oceny jakości zbioru danych Pima Indians Diabetes Dataset. Schemat obejmuje pozyskanie zbioru danych, ocenę jakości danych, obsługę wartości nieprawidłowych i brakujących, standaryzację cech numerycznych, podział zbioru danych na podzbiory treningowy, walidacyjny i niezależny testowy oraz końcową weryfikację jakości przed opracowaniem modelu. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 4: Krzywe uczenia treningowego i walidacyjnego proponowanego modelu z wykorzystaniem zbioru danych Pima Indians Diabetes Dataset. (A) Strata treningowa i walidacyjna w całym okresie uczenia. (B) Dokładność treningowa i walidacyjna w całym okresie uczenia. (C) Powiększony widok straty w epokach 50–10. (D) Powiększony widok dokładności w epokach 50–10. Najlepszą wydajność walidacyjną uzyskano w 78. epoce, przy stracie walidacyjnej wynoszącej 0,142 i dokładności walidacyjnej na poziomie 94,6%. Wcześniejsze zatrzymanie (early stopping) zostało uruchomione w 8. epoce przy zastosowaniu cierpliwości (patience) wynoszącej 10 epok. Aby zobaczyć powiększoną wersję tej ryciny, kliknij tutaj.

Rycina 5: Eksperymentalna ocena wydajności predykcyjnej proponowanego frameworka XAI z wykorzystaniem niezależnego zbioru testowego (n = 154). (A) Krzywa charakterystyki operacyjnej odbiornika (ROC) przedstawiająca wydajność dyskryminacyjną proponowanego modelu XAI oraz modeli bazowych. (B) Krzywe precyzji i pełności (PR) przedstawiające wydajność precyzji oraz pełności proponowanego modelu XAI i modeli bazowych. (C) Macierz pomyłek proponowanego modelu XAI na niezależnym zbiorze testowym, wraz z odpowiadającą im dokładnością, czułością (pełnością) i swoistością. (D) Porównanie dokładności, precyzji, pełności, swoistości, miary F1-score, współczynnika korelacji Matthewsa (MCC), pola pod krzywą ROC (AUC) oraz średniej precyzji (AP) we wszystkich ocenianych modelach. Wszystkie wyniki ilościowe przedstawione na tej rycinie odpowiadają eksperymentowi walidacyjnemu na zbiorze danych Pima Indians Diabetes Dataset. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 6: Wyniki wyjaśnialności wygenerowane z predykcji niezależnego zbioru testowego dla proponowanego modelu wytrenowanego na zbiorze danych Pima Indians Diabetes Dataset. (A) Globalny wykres podsumowujący SHAP pokazujący rozkład wkładu cech w całym zbiorze testowym. (B) Wykres istotności cech SHAP oparty na średnich bezwzględnych wartościach SHAP. (C) Indywidualny dla pacjenta wykres kaskadowy (waterfall plot) SHAP pokazujący wkład poszczególnych cech w przewidywane prawdopodobieństwo cukrzycy. (D) Lokalna interpretacja LIME pokazująca wkład cech dla pacjenta testowego nr 125. (E) Wykres zależności SHAP pokazujący relację między wartościami glukozy a ich wkładem SHAP. (F) Indywidualna dla pacjenta interpretacja kontrfaktyczna pokazująca minimalne zmiany cech związane ze zmianą predykcji modelu. Skróty: SHAP = Shapley Additive exPlanations; LIME = Local Interpretable Model-agnostic Explanations. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 7: Eksperymentalna walidacja statystyczna i analiza powtarzalności zaproponowanego modelu z wykorzystaniem zbioru danych Pima Indians Diabetes Dataset. (A) Wyniki pięciokrotnej walidacji krzyżowej na zbiorze treningowym. (B) 95% przedziały ufności bootstrap dla wyników na niezależnym zbiorze testowym. (C) Porównania statystyczne z modelami bazowymi, zawierające zastosowany test statystyczny, statystykę testową, wartość p oraz poziom istotności. (D) Spójność wyników w 10 niezależnych uruchomieniach z wykorzystaniem różnych ziarn generatora liczb pseudolosowych. Do oceny sparowanej dokładności klasyfikacji wykorzystano test McNemara, do skorelowanego ROC-AUC test DeLonga, a do porównania miary F1-score zastosowano sparowany test bootstrap z 1 0 próbami ponownymi. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.
| Etap | Aktywność protokołu | Oczekiwany wynik | Status wdrożenia |
| 1 | Wybór i walidacja zbioru danych z zakresu opieki zdrowotnej | Zweryfikowany zbiór danych, cel predykcji, rozkład klas oraz informacje o jakości danych | Wykonano |
| 2 | Konfiguracja środowiska obliczeniowego | Zweryfikowany sprzęt, oprogramowanie, biblioteki, wersje i konfiguracja obliczeniowa | Wykonano |
| 3 | Preprocessing i kontrola jakości danych opieki zdrowotnej | Oczyszczony, przekształcony i zestandaryzowany zbiór danych | Wykonano |
| 4 | Podziel zbiór danych | Niezależne zbiory danych treningowych, walidacyjnych i testowych | Wykonano |
| 5 | Opracowanie i optymalizacja modelu predykcyjnego/XAI | Ostateczna architektura modelu i hiperparametry | Wykonano |
| 6 | Trenowanie i zapisywanie modelu | Wytrenowany model, punkt kontrolny, konfiguracja trenowania i logi | Wykonano |
| 7 | Ocena wydajności predykcyjnej i obliczeniowej | Metryki wydajności zbioru testowego i porównania wydajności | Wykonano |
| 8 | Generowanie i ocena wyników wyjaśnialności | SHAP/LIME i stosowne wyniki wyjaśnień | Wykonano |
| 9 | Przeprowadzenie walidacji statystycznej i oceny powtarzalności | Przedziały ufności, testy statystyczne, wyniki powtórzeń i miary powtarzalności | Wykonano |
Tabela 1: Podsumowanie dziewięcioetapowego schematu podstawowego protokołu zastosowanego w walidacji praktycznej z wykorzystaniem zbioru danych Pima Indians Diabetes Dataset. Tabela rozróżnia dziewięć etapów zaimplementowanych w przykładzie praktycznym dla zbioru Pima Indians Diabetes Dataset od walidacji zewnętrznej i oceny eksperckiej, które pozostają opcjonalnymi elementami ogólnego protokołu.
| Zbiór danych | Źródło danych | Zastosowanie kliniczne | Modalność danych | Podmioty/Rekordy | Próbki | Klasy | Rozkład klas | Trening/Walidacja/Test | Rola w niniejszym badaniu | Status walidacji | URL źródła |
| Pima Indians Diabetes Dataset | UCI Machine Learning Repository | Przewidywanie cukrzycy | Kliniczne dane tabelaryczne | 768 rekordów | 768 | 2 | 50 bez cukrzycy; 268 z cukrzycą | 60% / 20% / 20% | Główny zbiór danych do walidacji eksperymentalnej | Wykonano – pełny dziewięcioetapowy podstawowy przepływ pracy | https://archive.ics.uci.edu/dataset/34/pima+indians+diabetes |
| TCGA-BRCA | NCI Genomic Data Commons (GDC), projekt TCGA-BRCA | Klasyfikacja raka piersi | Kliniczne + histopatologiczne | 1 098 przypadków | Zależne od zbioru danych w zależności od typu danych | Zależne od punktu końcowego | Brak jednolitego rozkładu klas dla całego zbioru danych | Nie dotyczy – nie wykonano podziału eksperymentalnego | Wyłącznie przykład zastosowania protokołu | Nie użyto do walidacji eksperymentalnej | https://portal.gdc.cancer.gov/projects/TCGA-BRCA |
| TCGA-UCEC | NCI Genomic Data Commons (GDC), projekt TCGA-UCEC | Klasyfikacja raka endometrium | Kliniczne + histopatologiczne | Kohorta projektu; rozmiar zależy od wybranego typu danych i filtrów | Zależne od zbioru danych w zależności od typu danych | Zależne od punktu końcowego | Brak jednolitego rozkładu klas dla całego zbioru danych | Nie dotyczy – nie wykonano podziału eksperymentalnego | Wyłącznie przykład zastosowania protokołu | Nie użyto do walidacji eksperymentalnej | https://portal.gdc.cancer.gov/projects/TCGA-UCEC |
| MIMIC-III | PhysioNet, MIMIC-III Clinical Database v1.4 | Przewidywanie wyników klinicznych | Kliniczne szeregi czasowe | 46 520 pacjentów | 58 976 przyjęć na OIOM | Zależne od zadania | Brak jednolitego rozkładu klas dla całej bazy danych | Nie dotyczy – nie wykonano podziału eksperymentalnego | Wyłącznie przykład zastosowania protokołu | Nie użyto do walidacji eksperymentalnej | https://physionet.org/content/mimiciii/1.4/ |
| ISIC 2019 | International Skin Imaging Collaboration (ISIC) Challenge | Klasyfikacja zmian skórnych | Obrazy dermoskopowe | Nie dotyczy – zbiór obrazów | 25 31 obrazów treningowych | 8 kategorii treningowych | AKIEC 867; BCC 3 323; BKL 2 624; DF 239; MEL 4 52; NV 12 875; VASC 253; SCC 628 | Nie dotyczy – nie wykonano podziału eksperymentalnego | Wyłącznie przykład zastosowania protokołu | Nie użyto do walidacji eksperymentalnej | https://challenge.isic-archive.com/landing/2019/ |
| HAM10 | Harvard Dataverse / ISIC Archive | Klasyfikacja zmian skórnych | Obrazy dermoskopowe | 7 470 unikalnych zmian | 10 015 obrazów | 7 | AKIEC 327; BCC 514; BKL 1 09; DF 15; MEL 1 13; NV 6 705; VASC 142 | Nie dotyczy – nie wykonano podziału eksperymentalnego | Wyłącznie przykład zastosowania protokołu | Nie użyto do walidacji eksperymentalnej | https://doi.org/10.7910/DVN/DBW86T |
| OhioT1DM | Zbiór danych OhioT1DM, udostępniony publicznie do celów badawczych | Monitorowanie/przewidywanie cukrzycy | Kliniczne szeregi czasowe | 12 uczestników | 24 pliki XML w obrębie danych treningowych/rozwojowych i testowych | Przewidywanie ciągłego poziomu glukozy; nie jest to stałe zadanie klasyfikacji | Nie dotyczy | Pliki treningowe/rozwojowe i testowe zdefiniowane w zbiorze danych; nie wykonano tutaj podziału eksperymentalnego | Wyłącznie przykład zastosowania protokołu | Nie użyto do walidacji eksperymentalnej | https://pmc.ncbi.nlm.nih.gov/articles/PMC781904/ |
| BraTS 2021 | RSNA-ASNR-MICCAI BraTS 2021; CBICA/University of Pennsylvania | Segmentacja/klasyfikacja guzów mózgu | MRI | 2 040 przypadków w kohorcie konkursowej | 1 251 adnotowanych przypadków treningowych; cztery modalności MRI na przypadek | Zależne od zadania | Brak jednolitego rozkładu klas dla całego zbioru danych | Kohorty zdefiniowane w ramach konkursu; nie wykonano tutaj podziału eksperymentalnego | Wyłącznie przykład zastosowania protokołu | Nie użyto do walidacji eksperymentalnej | https://www.med.upenn.edu/cbica/brats2021/ |
Tabela 2: Charakterystyka zbiorów danych z zakresu opieki zdrowotnej i zastosowalność proponowanego protokołu. Tabela podsumowuje źródła danych, zastosowania kliniczne, modalności, charakterystykę próbek, rozkłady klas, strategie partycjonowania zbiorów danych, status walidacji oraz informacje o źródłach dla zbiorów danych z zakresu opieki zdrowotnej uwzględnionych w protokole. Zbiór danych Pima Indians Diabetes służy jako główny przykład praktyczny do walidacji protokołu.
| Element konfiguracyjny | Rzeczywiste środowisko walidacyjne | Status / Interpretacja |
| Zbiór danych | Zbiór danych o cukrzycy wśród Indian Pima | Rzeczywisty zbiór danych do walidacji eksperymentalnej |
| Algorytm / Model | Tabularny model predykcyjny do binarnej klasyfikacji cukrzycy | Użyj dokładnej nazwy architektury z zapisu eksperymentu, jeśli została ona udokumentowana osobno |
| Współczynnik uczenia | 0.001 | Układ eksperymentalny |
| Optymalizator | Adam | Ustawienie eksperymentalne |
| Wielkość partii | 32 | Układ eksperymentalny |
| Maksymalna liczba epok | 100 | Układ doświadczalny |
| Funkcja straty | Entropia krzyżowa | Ustawienia eksperymentalne |
| Funkcja aktywacji | ReLU | Układ eksperymentalny |
| Dropout | 0.5 | Układ eksperymentalny |
| Rozkład wag | 1e-4 | Układ eksperymentalny |
| Normalizacja wsadowa | Włączono | Układ doświadczalny |
| Augmentacja danych / Równoważenie klas | Włączono | Określ SMOTE tylko wtedy, gdy faktycznie został zastosowany w raportowanym przebiegu |
| Strategia walidacji | 5-krotna walidacja krzyżowa | Układ doświadczalny |
| Wczesne zatrzymanie | Cierpliwość (patience) = 10 epok | Układ doświadczalny |
| Ziarno losowe | 42 | Należy użyć dokładnie tej samej konfiguracji ziarna (seed), która została zapisana dla danego eksperymentu |
| Powtórne uruchomienia | 10 niezależnych uruchomień z wykorzystaniem różnych ziarn generatora liczb pseudolosowych | Analiza powtarzalności eksperymentalnej |
| Rozmiar obrazu wejściowego | Nie dotyczy | Zbiór danych Pima ma strukturę tabelaryczną |
| Wymiar cechy | 512 | Stosować wyłącznie w przypadku, gdy jest to ostateczna zaimplementowana reprezentacja cech |
| Wyjaśnialność | SHAP + LIME; wyjaśnienie kontrfaktyczne przedstawione na rysunku 6 | Rzeczywista przeprowadzona analiza XAI |
| Metryki ewaluacji | Dokładność, precyzja, pełność, swoistość, wynik F1, współczynnik MCC, AUC-ROC, AP | Zgłoszone metryki ewaluacji eksperymentalnej |
| Sprzęt | procesor GPU NVIDIA | Zastąp dokładnym modelem GPU, jeśli został on odnotowany |
| Oprogramowanie / Framework | Python 3.1; Scikit-learn; komponenty frameworka wykorzystane w implementacji | Należy używać dokładnych wersji pakietów, jeśli zostały one odnotowane |
Tabela 3: Środowisko obliczeniowe, architektura modelu i konfiguracja hiperparametrów wykorzystane do implementacji protokołu. Tabela określa platformę obliczeniową, środowisko programistyczne, architekturę modelu, konfigurację danych wejściowych, parametry optymalizacji, ustawienia regularyzacji oraz parametry odtwarzalności wykorzystane do implementacji proponowanego przepływu pracy XAI.
| Parametr | Reprezentatywna specyfikacja / wynik |
| Optymalizator | Adam |
| Szybkość uczenia | 0.01 |
| Wielkość partii | 32 |
| Maksymalna liczba epok | 10 |
| Cierpliwość wczesnego zatrzymania | 10 epochs |
| Najlepsza epoka | 78 |
| Epoka wczesnego zatrzymania | 8 |
| Najniższa strata walidacyjna | 0.142 |
| Najwyższa dokładność walidacyjna | 94.6% |
| Wynik uczenia | Strata treningowa i walidacyjna zmalała i zbiegła się |
| Wynik walidacji | Dokładność treningowa i walidacyjna wzrosła i ustabilizowała się |
| Wynik optymalizacji | Najlepsza wydajność modelu osiągnięta w epoce 78 |
| Kontrola przeuczenia | Wczesne zatrzymanie zapobiegło dalszej optymalizacji poza wybraną najlepszą epoką |
Tabela 4: Specyfikacje trenowania modelu i wyniki optymalizacji. Tabela zawiera podsumowanie optymalizatora, tempa uczenia, wielkości partii, maksymalnej liczby epok trenowania, wyników walidacji, zbieżności trenowania, rezultatu optymalizacji oraz strategii kontroli przeuczenia zastosowanych podczas opracowywania modelu.
| Model | Dokładność (%) | Precyzja (%) | Czułość (%) | Swoistość (%) | Wynik F1 (%) | MCC | AUC (ROC) | AP (PR) |
| Proponowany model XAI | 93.5 | 94.5 | 92.4 | 94.6 | 93.4 | 0.87 | 0.96 | 0.94 |
| Głębokie uczenie (CNN) | 89.1 | 89.8 | 8.1 | 90 | 8.9 | 0.78 | 0.92 | 0.9 |
| Random Forest | 84.3 | 85 | 83.2 | 85.7 | 84.1 | 0.67 | 0.86 | 0.81 |
| Maszyna wektorów nośnych (SVM) | 78.6 | 79.3 | 7.1 | 80 | 78.2 | 0.54 | 0.79 | 0.72 |
| Baseline (klasa większościowa) | 62.1 | 62.1 | 10 | 0 | 76.6 | 0 | 0.5 | 0.5 |
Tabela 5: Porównanie wydajności predykcyjnej ocenianych modeli. Tabela porównuje proponowany model XAI z ocenianymi modelami bazowymi z wykorzystaniem dokładności, precyzji, czułości, swoistości, wyniku F1, współczynnika korelacji Matthewsa, pola pod krzywą charakterystyki operacyjnej odbiornika oraz średniej precyzji.
| Analiza walidacyjna | Porównanie / Metryka | Test statystyczny | Statystyka testowa | Wartość p | Wynik eksperymentalny / 95% CI |
| Pięciokrotna walidacja krzyżowa | Dokładność | Opisowa (średnia ± SD) | — | — | 93,01 ± 0,48% |
| Pięciokrotna walidacja krzyżowa | AUC-ROC | Opisowa (średnia ± SD) | — | — | 0,954 ± 0,07 |
| Pięciokrotna walidacja krzyżowa | Precyzja | Opisowa (średnia ± SD) | — | — | 92,42 ± 0,87% |
| Pięciokrotna walidacja krzyżowa | Czułość | Opisowa (średnia ± SD) | — | — | 93,02 ± 0,45% |
| Pięciokrotna walidacja krzyżowa | Wskaźnik F1 | Opisowa (średnia ± SD) | — | — | 92,72 ± 0,62% |
| 95% przedział ufności bootstrap | Dokładność | Bootstrap (1 0 prób) | — | — | 0,935 [0,897, 0,973] |
| 95% przedział ufności bootstrap | Precyzja | Bootstrap (1 00 prób) | — | — | 0,930 [0,890, 0,970] |
| 95% przedział ufności bootstrap | Czułość | Bootstrap (1 0 prób) | — | — | 0,92 [0,80, 0,963] |
| 95% przedział ufności bootstrap | Wskaźnik F1 | Bootstrap (1 0 prób) | — | — | 0,926 [0,87, 0,965] |
| 95% przedział ufności bootstrap | AUC-ROC | Bootstrap (1 0 prób) | — | — | 0,958 [0,931, 0,984] |
| Proponowany model vs CNN | Dokładność (%) | Test McNemara | 9,32 | 0,023 | Istotne (α = 0,05) |
| Proponowany model vs CNN | AUC-ROC | Test DeLonga | 3,87 | 0,01 | Istotne (α = 0,05) |
| Proponowany model vs CNN | Wskaźnik F1 | Parowany bootstrap (1 0 prób) | 2,81 | 0,04 | Istotne (α = 0,05) |
| Proponowany model vs Random Forest | Dokładność (%) | Test McNemara | 14,27 | 0,002 | Istotne (α = 0,05) |
| Proponowany model vs Random Forest | AUC-ROC | Test DeLonga | 5,86 | <0,01 | Istotne (α = 0,05) |
| Proponowany model vs Random Forest | Wskaźnik F1 | Parowany bootstrap (1 00 prób) | 4,03 | 0,03 | Istotne (α = 0,05) |
| Proponowany model vs SVM | Dokładność (%) | Test McNemara | 16,08 | <0,01 | Istotne (α = 0,05) |
| Proponowany model vs SVM | AUC-ROC | Test DeLonga | 6,95 | <0,001 | Istotne (α = 0,05) |
| Proponowany model vs SVM | Wskaźnik F1 | Parowany bootstrap (1 0 prób) | 4,62 | <0,01 | Istotne (α = 0,05) |
| Powtarzalność serii pomiarowych (10 niezależnych uruchomień) | AUC-ROC | Opisowa (średnia ± SD) | — | — | 0,957 ± 0,08 |
| Powtarzalność serii pomiarowych (10 niezależnych uruchomień) | Dokładność (%) | Opisowa (średnia ± SD) | — | — | 93,24 ± 0,74% |
| Powtarzalność serii pomiarowych (10 niezależnych uruchomień) | Wskaźnik F1 (%) | Opisowa (średnia ± SD) | — | — | 92,35 ± 0,92% |
Tabela 6: Wyniki walidacji statystycznej i powtarzalności uzyskane z implementacji eksperymentalnej z wykorzystaniem zbioru danych Pima Indians Diabetes Dataset. Tabela podsumowuje wydajność pięciokrotnej walidacji krzyżowej, 95% przedziały ufności oparte na metodzie bootstrap, porównania statystyczne proponowanego modelu z modelami bazowymi oraz powtarzalność wielokrotnych uruchomień dla 10 niezależnych ziarn losowości. W ramach niniejszej walidacji nie przeprowadzono walidacji zewnętrznej ani oceny przez ekspertów klinicznych.
| Nr | Element listy kontrolnej | Wymagana dokumentacja | Zalecany zapis / weryfikacja |
| 1 | Środowisko oprogramowania | System operacyjny, język programowania i środowisko oprogramowania | Zapisz dokładne wersje systemu operacyjnego i języka programowania. |
| 2 | Wersje oprogramowania i bibliotek | Wersje głównych bibliotek i pakietów oprogramowania | Zapisz dokładne nazwy i wersje pakietów/bibliotek. |
| 3 | Specyfikacja sprzętowa | Specyfikacje CPU, GPU, RAM, pamięci masowej i akceleratorów | Zapisz użyty sprzęt obliczeniowy. |
| 4 | Identyfikacja zbioru danych | Nazwa zbioru danych, źródło, wersja i informacje o dostępie | Zapisz dokładne źródło/wersję zbioru danych oraz datę dostępu, jeśli dotyczy. |
| 5 | Charakterystyka zbioru danych | Wielkość próby i rozkład klas | Zapisz całkowitą liczbę próbek i rozkład klas dla każdego podziału. |
| 6 | Partycjonowanie zbioru danych | Strategia podziału na zbiór treningowy, walidacyjny i niezależny zbiór testowy | Udokumentuj proporcje/liczbę próbek w podziałach oraz stratyfikację. |
| 7 | Zapobieganie wyciekowi danych | Procedura stosowana w celu zapobieżenia wyciekowi informacji | Udokumentuj separację obiektów/próbek oraz szacowanie parametrów preprocessingu wyłącznie na zbiorze treningowym. |
| 8 | Parametry preprocessingu | Ustawienia czyszczenia, obsługi brakujących wartości, normalizacji, kodowania i transformacji | Zapisz wszystkie operacje preprocessingu i wartości parametrów. |
| 9 | Augmentacja danych | Metody augmentacji i ustawienia parametrów, jeśli dotyczy | Udokumentuj metody, prawdopodobieństwa i zakresy parametrów. |
| 10 | Architektura modelu | Końcowa architektura i konfiguracja modelu | Udokumentuj typ modelu, warstwy/komponenty, wymiary i funkcje aktywacji. |
| 1 | Hiperparametry | Hiperparametry treningu i optymalizacji | Zapisz szybkość uczenia, rozmiar partii (batch size), optymalizator, liczbę epok, wczesne zatrzymywanie (early stopping) i dostrojone parametry. |
| 12 | Ziarna losowości (Random seeds) | Ziarna losowości użyte do zapewnienia powtarzalności wykonania | Zapisz ziarna dla podziałów, inicjalizacji i powtórzonych uruchomień. |
| 13 | Konfiguracja treningu | Procedura treningowa i strategia wyboru modelu | Udokumentuj walidację, tworzenie punktów kontrolnych (checkpointing) i kryteria wyboru modelu. |
| 14 | Metryki ewaluacji | Predefiniowane predykcyjne i statystyczne metryki ewaluacji | Zapisz wszystkie raportowane miary wydajności. |
| 15 | Przedziały ufności | Przedziały niepewności dla miar wydajności | Udokumentuj procedurę szacowania przedziału ufności (CI) i próby bootstrapowe, jeśli dotyczy. |
| 16 | Testy statystyczne | Procedury statystyczne do porównania modeli | Udokumentuj test, statystykę, wartość p, próg istotności oraz założenia. |
| 17 | Analiza powtórzonych uruchomień | Niezależne wykonania z użyciem różnych ziaren losowości | Zapisz liczbę uruchomień oraz średnią ± SD kluczowych metryk. |
| 18 | Konfiguracja wyjaśnialności | Metody wyjaśnialności i ustawienia parametrów | Udokumentuj SHAP, LIME, Grad-CAM, mechanizmy uwagi (attention), analizy kontrfaktyczne lub inne odpowiednie ustawienia. |
| 19 | Ocena wyjaśnialności | Obiektywna ocena niezawodności i użyteczności wyjaśnień | Udokumentuj wierność (faithfulness), stabilność, niefidelność, lokalizację i ocenę ekspercką, jeśli dotyczy. |
| 20 | Artefakty modelu | Wagi wytrenowanego modelu i pliki konfiguracyjne | Zarchiwizuj końcowy wytrenowany model, architekturę/konfigurację oraz informacje o wyborze modelu. |
| 21 | Dostępność kodu | Kod wymagany do preprocessingu, treningu, ewaluacji i generowania wyjaśnień | Zapisz repozytorium lub informacje o kontrolowanym dostępie do kodu, jeśli dotyczy. |
| 2 | Dokumentacja wykonania | Polecenia, skrypty, pliki konfiguracyjne i instrukcje | Zapisz polecenia i konfigurację niezbędną do odtworzenia przepływu pracy (workflow). |
| 23 | Dokumentacja wyników | Predykcje, wyniki ewaluacji, wykresy i wyniki wyjaśnień | Zarchiwizuj wygenerowane wyniki i powiąż je z odpowiednim eksperymentem/uruchomieniem. |
| 24 | Weryfikacja powtarzalności | Weryfikacja spójności między niezależnymi wykonaniami | Porównaj wyniki powtórzonych uruchomień i zaraportuj predefiniowane miary zmienności. |
Tabela 7: Lista kontrolna odtwarzalności dla niezależnej replikacji proponowanego przepływu pracy XAI. Lista kontrolna określa wymagania obliczeniowe, dotyczące zbiorów danych, preprocessingu, rozwoju modelu, ewaluacji, walidacji statystycznej, wyjaśnialności oraz dokumentacji niezbędne do odtworzenia eksperymentalnego przepływu pracy.