Jak działa test F-Snedecora?
Po co w ogóle to sprawdzać?
Test F-Snedecora sprawdza, czy dwie niezależne próby mają statystycznie istotnie
różną zmienność (wariancję), a nie tylko różną średnią. Typowe zastosowania
laboratoryjne: porównanie precyzji dwóch metod analitycznych, sprawdzenie, czy nowy
analityk ma podobną powtarzalność co doświadczony, porównanie dwóch urządzeń
pomiarowych, oraz walidacja metody. Test F bywa też krokiem pomocniczym przed
wyborem wariantu testu t dla dwóch prób niezależnych (klasyczny test t wymaga
zbliżonych wariancji obu grup, test Welcha - nie).
Wzór
F = s12 / s22,
df1 = n1 − 1, df2 = n2 − 1
gdzie s12 i s22
to wariancje próbkowe obu prób, a n1, n2
to ich liczności. W teście dwustronnym arkusz automatycznie umieszcza
WIĘKSZĄ z dwóch wariancji w liczniku (dzięki temu F zawsze ≥ 1, zgodnie z
powszechną konwencją i tablicami rozkładu F) - to, która próba (1 czy 2) faktycznie
trafiła do licznika, jest zawsze jawnie podane w wyniku. W teście
jednostronnym kolejność jest STAŁA (próba 1 zawsze w liczniku) -
przestawianie unieważniłoby kierunkowość takiego testu.
Uwaga o oznaczeniach: formalnie σ oznacza odchylenie
standardowe CAŁEJ populacji, a to, co faktycznie liczymy z Twoich pomiarów, jest tylko
jego OSZACOWANIEM na podstawie próby (poprawnie: s). Ten arkusz mimo to
konsekwentnie posługuje się symbolem σ - zarówno w hipotezach niżej, jak i w tabelach
wyniku PDF - bo to najbardziej naturalny symbol odchylenia standardowego dla laboranta
bez zaplecza statystycznego. Warto o tym rozróżnieniu pamiętać, jeśli porównujesz ten
raport z podręcznikiem lub inną publikacją, gdzie σ i s bywają rozróżniane ściślej.
Test dwustronny czy jednostronny?
Test dwustronny (H₀: σ₁² = σ₂², H₁: σ₁² ≠ σ₂²) sprawdza, czy
wariancje różnią się w dowolnym kierunku, i jest właściwym wyborem domyślnym, gdy
nie ma z góry założonego kierunku różnicy - typowe pytanie: "czy te dwie
metody/próbki mają podobną precyzję?".
Test jednostronny (H₀: σ₁² ≤ σ₂², H₁: σ₁² > σ₂²) sprawdza różnicę
tylko w jednym, wcześniej założonym kierunku - np. "czy metoda 1 ma GORSZĄ precyzję
(większą wariancję) niż referencyjna metoda 2?". Ma większą moc wykrywania różnicy w
tym kierunku, ale tylko wtedy, gdy kierunek został ustalony przed zebraniem
danych. Chcesz sprawdzić kierunek przeciwny (σ₁² < σ₂²)? Po prostu zamień, którą
próbę wpisujesz jako "1", a którą jako "2" - to dokładnie ten sam test, tylko z
odwróconymi etykietami, więc arkusz świadomie nie duplikuje tej opcji pod inną nazwą.
Dwa sposoby wprowadzania danych
Surowe obserwacje - każda wartość to JEDEN niezależny pomiar (bez
uśredniania powtórzeń analitycznych, w odróżnieniu od testu t-Studenta w tym
arkuszu!). Uśrednienie powtórzeń przed testem F zniekształciłoby wynik - test F
ocenia właśnie zmienność (precyzję), a uśrednienie tę zmienność częściowo usuwa. Z
tego samego powodu powtórzenia techniczne pojedynczej próbki (kilka odczytów tego
samego przygotowanego roztworu) nie powinny wchodzić do testu F jako osobne,
niezależne obserwacje - naruszałoby to założenie niezależności. Właściwym źródłem
danych jest seria NIEZALEŻNYCH oznaczeń (różne dni / różne przygotowania próbki /
różni analitycy - w zależności od tego, jaką zmienność akurat badasz).
Statystyki podsumowujące - podajesz wprost n oraz σ (lub wariancję) każdej
próby. Typowy przypadek, gdy ocena powtarzalności metody jest już gotowa w
dokumentacji walidacyjnej (seria niezależnych oznaczeń CRM podsumowana do n i σ).
Przełącznik σ/wariancja to wyłącznie kwestia interpretacji wpisywanej liczby - to
dokładnie ten sam test.
Decyzja
Obliczoną wartość F porównuje się z wartością krytyczną
Fkrit(α, df1, df2) rozkładu F-Snedecora. Jeżeli
F > Fkrit (równoważnie: p < α) - różnicę wariancji uznaje się
za statystycznie istotną. W przeciwnym razie nie ma podstaw, by odrzucić hipotezę
zerową o równości (lub nieprzewyższaniu) wariancji.
Wykres rozkładu F w raporcie
Raport PDF automatycznie zawiera wykres gęstości rozkładu F(df1, df2),
na którym zaznaczone są DWA obszary: pomarańczowy obszar krytyczny, odpowiadający
ZAŁOŻONEMU z góry poziomowi istotności (próg decyzji ustalony PRZED policzeniem testu),
oraz niebieski, kreskowany obszar odpowiadający wartości p WYNIKAJĄCEJ z
analizowanych danych. Pod wykresem program dodaje krótki opis (1–3 zdania),
napisany prostym językiem, tłumaczący wprost, co dokładnie widać na wykresie i jaki
płynie z tego wniosek – bez konieczności samodzielnego odczytywania wykresu czy
tablic rozkładu F.
Założenia testu
- Obie próby pochodzą z (w przybliżeniu) rozkładu normalnego.
- Pomiary w obrębie każdej próby, oraz między próbami, są niezależne od siebie.
- Wariancje są dodatnie i skończone (test nie ma sensu dla próby o zerowym rozrzucie).
- Zalecana liczebność próby to co najmniej n = 5 w każdej grupie – przy mniejszej
liczbie pomiarów oszacowanie zmienności jest mało stabilne, a trudniej też ocenić, czy
założenie normalności rozkładu w ogóle ma sens. Poniżej tego progu arkusz automatycznie
pokazuje w raporcie dodatkowe, żółte ostrzeżenie.
Uwaga: test F jest wyraźnie BARDZIEJ wrażliwy na odstępstwa od
normalności niż test t-Studenta - przy danych wyraźnie odbiegających od rozkładu
normalnego (zwłaszcza przy małych próbach) wynik testu F może być mylący. W takiej
sytuacji warto rozważyć alternatywne metody porównania zmienności (np. test Levene'a
lub test Browna-Forsythe'a, mniej wrażliwe na brak normalności), choć wykraczają one
poza ten arkusz.
Źródła
Fisher, R.A. (1924). On a Distribution Yielding the Error Functions of Several
Well Known Statistics. Proceedings of the International Congress of
Mathematics, Toronto, 2, 805–813.
Snedecor, G.W., Cochran, W.G. (1989). Statistical Methods (8th ed.). Iowa
State University Press.
PN-EN ISO 5725-2:2019 - Dokładność (poprawność i precyzja) metod pomiarowych i
wyników pomiarów - Część 2: Podstawowa metoda określania powtarzalności i
odtwarzalności znormalizowanej metody pomiarowej.