Badanie próbne
Próba obejmowała 881 uczestników z Hiszpanii (N = 325), Meksyku (N = 169), Gwatemali (N = 227) i Chile (N = 160), z których wszyscy byli rodzimymi użytkownikami języka hiszpańskiego. Próbę podzielono na dwie grupy: 451 osób w grupie z dysleksją (RD) oraz 430 osób w grupie czytelników o typowych osiągnięciach (NAR). Dzieci ze specjalnymi potrzebami edukacyjnymi – wymagające wsparcia i specyficznej uwagi edukacyjnej ze względu na dysfunkcje sensoryczne, problemy neurologiczne lub inne schorzenia – zostały wykluczone, ponieważ czynniki te są zazwyczaj stosowane jako kryteria wykluczające w przypadku trudności w nauce lub ciężkich zaburzeń behawioralnych, albo tymczasowo, albo na cały okres edukacji. Uczestnicy grupy RD zostali wybrani na podstawie kryteriów takich jak: IQ równe lub większe niż 80, centyl czasu czytania słów powyżej 75. centyla, centyl czasu czytania pseudosłów powyżej 75. centyla lub centyl poprawności czytania pseudosłów poniżej 25. centyla. Podobnie grupę NAR wybrano na podstawie porównywalnych kryteriów, z dodatkiem centyla zadania rozumienia tekstu PROLEC36 powyżej 25. centyla. Wszystkie centyle określono zgodnie z poziomem klasy uczestników. Klasyfikacja uczestników do tych grup opierała się na zadaniach czytania słów i czytania pseudosłów, z których każde składało się z oddzielnych bloków. W obu blokach uczestnicy otrzymali instrukcję, aby czytać na głos prezentowane bodźce werbalne tak szybko, jak to możliwe. Blok czytania słów obejmował 32 bodźce, natomiast blok czytania pseudosłów zawierał 48 bodźców. Znajomość słów kontrolowano zgodnie z kryteriami opisanymi przez Guzmana i Jimeneza37. Rejestrowano liczbę poprawnych odpowiedzi, błędów oraz czasy latencji, przy czym czasy latencji mierzono od pojawienia się każdego elementu do momentu rozpoczęcia odpowiedzi wokalnej ucznia. Blok czytania słów wykazał wysoką spójność wewnętrzną (α = .89), natomiast blok czytania pseudosłów wykazał jeszcze większą spójność wewnętrzną (α = .91). Rozkład uczniów z RD i NAR w poszczególnych klasach przedstawiał się następująco: w drugiej klasie było 30 uczniów NAR i 27 RD. W trzeciej klasie odnotowano 100 uczniów NAR oraz 88 RD. Czwarta klasa obejmowała 100 uczniów NAR i 116 RD. W piątej klasie było 100 uczniów NAR i 116 RD, natomiast w szóstej klasie 100 uczniów NAR i 104 RD. Aby ocenić różnice wiekowe w całej próbie oraz w obrębie każdego poziomu klasy, przeprowadzono jednoczynnikowe testy ANOVA. Wyniki nie wykazały istotnych różnic, a wszystkie wartości p przekroczały 0,05. Podobnie zbadano rozkład zmiennej płci za pomocą testu chi-kwadrat, gdzie wszystkie wartości p również przekroczyły 0,05.
Statystyka opisowa
Tabela 1 przedstawia statystyki opisowe (średnią i odchylenie standardowe) dla wieku oraz miar oceny. Dane zostały podzielone według płci, z rozróżnieniem między czytelnikami osiągającymi wyniki w normie a osobami z dysleksją w różnych krajach. Wyniki wykazały pozytywne i istotne statystycznie korelacje pomiędzy dużą liczbą miar oceny (patrz Tabela 2). Tabela przedstawia korelacje między miarami, obejmujące zakres od słabych do silnych. Korelacje poniżej 0,3 uznano za słabe, natomiast te powyżej 0,5 sklasyfikowano jako silne. Konkretne wartości korelacji podano wraz z ich odpowiednią klasyfikacją siły. Na przykład słabe korelacje (r < 0,3) obejmowały korelacje między blendowaniem a segmentacją (r = 0,354) oraz między usuwaniem a rozumieniem homofonów (r = 0,270). Umiarkowane korelacje (0,3 ≤ r < 0,5) zaobserwowano między strukturą gramatyczną a liczbą (r = 0,463) oraz słowami funkcyjnymi a strukturą gramatyczną (r = 0,512). Silne korelacje (r ≥ 0,5) są widoczne dla liczby i rodzaju (r = 0,642) oraz RAN obrazkowego i RAN kolorów (r = 0,442).
Tabela 1: Statystyki opisowe pomiarów. Tabela ta przedstawia średnią, odchylenie standardowe oraz inne statystyki opisowe dla poszczególnych pomiarów zawartych w narzędziu cyfrowym. AGE: wiek; NAR: czytelnicy osiągający standardowe wyniki (Normally Achieving Readers); RD: dysleksja/zaburzenia czytania; SEG: segmentacja; BLN: blendowanie; ISO: izolacja fonemów; DEL: usuwanie fonemów; HOM: rozumienie homofonów; R&S: rdzenie i przyrostki; GEN: rodzaj; NUM: liczba; FW: słowa funkcyjne; GRS: struktura gramatyczna; EXT: tekst ekspozycyjny; NAT: tekst narracyjny; VOC: dźwięczność artykulacji; MOA: sposób artykulacji; POA= miejsce artykulacji; DIG: RAN cyfr; LET: RAN liter; COL: RAN kolorów; PIC: RAN obrazków Prosimy kliknąć tutaj, aby pobrać tę tabelę.
Tabela 2: Współczynniki korelacji pomiarów.Tabela ta przedstawia współczynniki korelacji między różnymi pomiarami narzędzia cyfrowego, wskazując siłę i kierunek zależności między tymi pomiarami. SEG: Segmentacja; BLN: Blendowanie; ISO: Izolacja fonemów; DEL: Usuwanie fonemów; HOM: Rozumienie homofonów; R&S: Rdzenie i przyrostki; GEN: Rodzaj; NUM: Liczba; GRS: Struktura gramatyczna; FW: Słowa funkcyjne; EXT: Tekst ekspozycyjny; NAT: Tekst narracyjny; VOC: Głosowość artykulacji; MOA: Sposób artykulacji; POA= Miejsce artykulacji; DIG: RAN cyfr; LET: RAN liter; COL: RAN kolorów; PIC: RAN obrazków. Kliknij tutaj, aby pobrać tę tabelę.
Weryfikująca analiza czynnikowa
Przeprowadzono weryfikującą analizę czynnikową (CFA), aby ocenić proponowaną strukturę czynnikową baterii multimedialnej. Model obejmuje jeden czynnik drugiego rzędu oraz sześć zmiennych latentnych, z których każda reprezentuje odrębne moduły baterii multimedialnej. Konstrukty te obejmują świadomość fonologiczną (w tym zadania z segmentacji fonemicznej, izolacji, blendowania i usuwania), moduł morfo-ortograficzny (obejmujący zadania z rozumienia rdzeni, afiksów i homofonów), składnię (składającą się z zadań dotyczących rodzaju, liczby, struktury gramatycznej i słów funkcyjnych), percepcję mowy (obejmującą zadania dotyczące dźwięczności, sposobu artykulacji i miejsca artykulacji), rozumienie tekstu (obejmujące zadania z rozumienia tekstu ekspozycyjnego i narracyjnego) oraz szybkie automatyczne nazywanie (zadania z literami, kolorami, cyframi i obrazkami w ramach testu RAN). Aby zapewnić spójność wszystkich zadań w baterii multimedialnej (w zakresie od 0 do 1), zastosowano estymację z wykorzystaniem skalowania proporcji maksimum (POMS). Wyniki POMS obliczono dla zadań pomiarowych czasu (rdzenie i afiksy)38.
Analizy statystyczne i prezentacje graficzne wykonano przy użyciu oprogramowania R wersja 4.3.139, wykorzystując pakiety lavaan40, semTools41 oraz ggplot242. Dopasowanie modelu oceniono za pomocą różnych wskaźników dobroci dopasowania. Choć test chi-kwadrat dla dobroci dopasowania był istotny, χ2(df) = 632.01, p < .001, co sugeruje rozbieżność między modelem hipotetycznym a obserwowanymi danymi, należy zaznaczyć, że test chi-kwadrat jest wrażliwy na wielkość próby43. W związku z tym rozważono dodatkowe wskaźniki dopasowania.
Wskaźnik dopasowania porównawczego (CFI) wyniósł .961, co przekracza powszechnie przyjęty próg .95 i wskazuje na dobre dopasowanie do danych. Średniokwadratowy błąd aproksymacji (RMSEA) wyniósł .038, co wskazuje na ścisłe dopasowanie modelu do danych. Standaryzowana pierwiastkowa średnia kwadratowa reszt (SRMR) wyniosła .034, co znajduje się poniżej zalecanej wartości .05, sugerując idealne dopasowanie. Ładunki czynnikowe dla pozycji wynosiły od .36 do .81, co wskazuje na istotne ładunki dla ich odpowiednich czynników. Całkowita średnia wariancja wyekstrahowana (AVE) przekroczyła .50, co sugeruje odpowiednią trafność zbieżną, a całkowita niezawodność kompozytowa (CR) była powyżej .80, co wskazuje na dobrą niezawodność spójności wewnętrznej.
Podsumowując, wyniki CFA potwierdziły proponowaną strukturę czynnikową, wykazując dobry dopasowanie modelu, trafność zbieżną oraz rzetelność. Przedstawienie graficzne modelu znajduje się na Rysunku 1.

Rysunek 1: Analiza czynnikowa potwierdzająca. Na rysunku przedstawiono wyniki analizy czynnikowej potwierdzającej przeprowadzonej dla narzędzia cyfrowego, z zaznaczeniem ładunków czynnikowych oraz relacji między różnymi miarami przetwarzania poznawczego. Skróty: SIC= Sicole-R; PA= świadomość fonologiczna; SEG= segmentacja; BLN= blendowanie; ISO= izolacja; DEL= usuwanie; MO= przetwarzanie morfologiczne; HOM= rozumienie homofonów; R&S= rdzenie i afiksy; SYN= przetwarzanie syntaktyczne; GEN= rodzaj; NUM= liczba; GRS= struktura gramatyczna; FW= słowa funkcyjne; RC= rozumienie tekstu; EXT= tekst ekspozycyjny; NAT= tekst narracyjny; SP= percepcja mowy; VOC= dźwięczność artykulacji; MOA= sposób artykulacji; POA= miejsce artykulacji; RAN= szybkie nazywanie zautomatyzowane; NUM= RAN liczb; LET= RAN liter; COL= RAN kolorów; PIC= RAN obrazków. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.
Analizy niezmienniczości pomiarowej wykorzystano do sprawdzenia, czy struktura czynnikowa baterii multimedialnej była stabilna w zależności od płci. Testowanie niezmienniczości pomiarowej polega na serii porównań modeli, które definiują coraz bardziej rygorystyczne ograniczenia równości44. Przeprowadzono analizy statystyczne, stosując te same kryteria dopasowania modelu, co w poprzedniej analizie CFA. Stopniowo nakładano ograniczenia na parametry reprezentujące struktury: konfiguracyjną, metryczną (ładunki), skalarną (przechwyty) oraz ścisłą (residua)45. Słabe dopasowanie w którymkolwiek z tych modeli sugeruje, że aspekt poddany ograniczeniu nie funkcjonuje spójnie w różnych grupach. Stopień niezmienniczości określano łącznie, gdy χ2D wynosiło >0,05, a ΔCFI było <0,01. Podsumowanie wyników wskaźników modeli oraz różnic między nimi przedstawiono w Tabeli 3.
Tabela 3: Wskaźniki dopasowania modelu niezmienniczego.Tabela ta przedstawia wskaźniki dopasowania służące do testowania niezmienniczości modelu ze względu na płeć, oceniając, w jakim stopniu narzędzie cyfrowe zachowuje spójność swojej struktury i właściwości pomiarowych w grupach mężczyzn i kobiet. Wskaźniki te zapewniają niezawodność i trafność narzędzia, potwierdzając jego zdolność do spójnego pomiaru procesów poznawczych w różnych grupach płciowych. Kliknij tutaj, aby pobrać tę tabelę.
Model konfiguracyjny, w którym ograniczono jedynie konfigurację względną zmiennych w modelu, aby była taka sama w obu grupach, wykazał adekwatne dopasowanie do danych: χ2(292) = 745.970, p<.001, CFI = .964, SRMR = 0.034, RMSEA = .036, 90% CI (.033, .038). Model niezmienniczości metrycznej ograniczał konfigurację zmiennych oraz wszystkie ładunki czynnikowe, aby były stałe w poszczególnych grupach. Wskaźniki dopasowania były porównywalne z tymi z modelu konfiguracyjnego: χ2(310) = 768.56, p<.001, CFI = .963, SRMR = .037, RMSEA = .036, 90% CI (.033, .038). Niezmienniczość ładunków czynnikowych została potwierdzona przez nieistotne testy różnic oceniające podobieństwo modeli: χ2D (18) = 26.30, p =.09; ΔCFI = .001. W modelu niezmienniczości skalarnej konfiguracja, ładunki czynnikowe oraz średnie/punkty przecięcia wskaźników zostały ograniczone tak, aby były takie same dla każdej grupy. Wskaźniki dopasowania były mniej niż idealne: χ2(322) = 787,50, p<.001, CFI = .963, SRMS = .038, RMSEA = .035 (90% CI = .032, .038). Testy różnic oceniające podobieństwo modeli sugerowały występowanie niezmienniczości factorialnej: χ2D(12)=13.00, p = .369; ΔCFI = .001. Wreszcie, w modelu niezmienniczości ścisłej konfiguracja, ładunki czynnikowe, średnie/punkty przecięcia wskaźników oraz reszty zostały ograniczone tak, aby były takie same dla każdej grupy. Wskaźniki dopasowania były mniej niż idealne: χ2(341) =798.20, p<.001, CFI = .961, SRMS = .039, RMSEA = .0035 (90% CI = .032, .038). Niezmienniczość ścisła została potwierdzona przez nieistotne testy różnic oceniające podobieństwo modeli: χ2D(11) = 25.81, p<.001; ΔCFI =.002
Dokładność diagnostyczna
Aby ocenić dokładność i zdolność dyskryminacyjną baterii multimedialnej, przeprowadzono analizy krzywej charakterystyki operacyjnej odbiornika (ROC). Analiza ROC pomaga w wyznaczeniu optymalnego progu (wartości odcięcia) dla testu oceniającego w skali ciągłej, zapewniając równowagę między czułością (zdolnością do poprawnej identyfikacji wyników prawdziwie dodatnich) a swoistością (zdolnością do poprawnej identyfikacji wyników prawdziwie ujemnych). Dodatkowo analizę ROC wykorzystano do oceny zdolności testu do rozróżniania grup RD i NAR. W celu przeprowadzenia analizy obliczono wyniki Z dla każdego stopnia w każdym zadaniu baterii multimedialnej. Wynik całkowity stanowił suma wyników zeta. Analizy statystyczne i prezentacje graficzne wykonano przy użyciu programu R w wersji 4.3.1.41Pakiet pROC46 oraz ggplot242 wykorzystano pakiety. Pod względem dokładności diagnostycznej bateria multimedialna wykazała pole pod krzywą (AUC) wynoszące 94,39% [95% CI: 93,31%–96,24% (metoda DeLonga)] oraz czułość na poziomie 91,0 (Rycina 2Krzywe ROC w podziale na stopień zaawansowania wykazały następujące wskaźniki: stopień 2°, AUC = 96,195% [CI: 91,54%–100% (DeLong)], Se = 0,96, Sp = 0,90; stopień 3°, AUC = 95,3 [95% CI: 92,43%–98,18% (DeLong)], Se = 75,70, Sp = 72,34; 4° stopień, AUC=93,4 [95% CI: 90,2%–96,66% (DeLong)], Se=0,92, Sp=0,84; 5° stopień, AUC=95,9 [// 95% CI: 93,11%–98,75% (DeLong)], Se=0,90, Sp=0,95; 6° stopień, AUC=94,4 [95% CI: 91,11%–97,69% (DeLong)], Se=0,92, Sp=0,91.

Rycina 2: Analiza krzywej ROC. Na rysunku przedstawiono analizę krzywej ROC, która obrazuje dokładność diagnostyczną narzędzia cyfrowego poprzez zestawienie odsetka wyników prawdziwie dodatnich z odsetkiem wyników fałszywie dodatnich dla różnych wartości progowych. Jednostką na osi x krzywej ROC jest swoistość, a na osi y — czułość. Skróty: ROC = charakterystyka operacyjna odbiornika; AUC = pole pod krzywą. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.