6. W jakich przypadkach stosuje się jednoczynnikową analizę wariancji? Jakie założenia należy
sprawdzić przed użyciem tej techniki statystycznej? Podaj przykład zastosowania jednoczynnikowej
analizy wariancji.
Jednoczynnikowa analiza wariancji jest testem statystycznym służącym do porównywania średnich
w wielu populacjach. Stosujemy ją gdy chcemy ocenić czy występujące różnice są istotne
statystycznie gdy mamy jeden czynnik eksperymentalny.
Założenia:
·ð Każda populacja musi mieć rozkÅ‚ad normalny
·ð Pobrane do analizy próby sÄ… niezależne
·ð Próby pobrane z każdej populacji muszÄ… być losowymi próbami prostymi
·ð Wariancje w populacjach sÄ… równe
W przypadku, gdy założenia analizy wariancji nie są spełnione należy posługiwać się testem Kruskala-
Wallisa.
Przykład z zajęć:
Zbadać czy rodzaj paszy podawany kurczakom wpływa na przyrost ich wagi.
Mamy trzy rodzaje paszy: A, B oraz C.
Plan doświadczenia: układ całkowicie losowy ( kompletnie zrandomizowany) zrównoważony, czyli
jednakowa liczba kurczaków przypisana do trzech pasz.
Czynnik doświadczalny (inaczej eksperymentalny): rodzaj paszy (A, B, C); każda pasza jest obiektem
doświadczalnym, natomiast jednostki eksperymentalne to kurczaki.
µA- Å›redni przyrost dla paszy A
µB- Å›redni przyrost dla paszy B
µC- Å›redni przyrost dla paszy C
H : wszystkie Å›rednie grupowe sÄ… identyczne; µA= µB= µC
0
vs. H : istnieje różnica w przyrostach wagi dla przynajmniej dwóch pasz
1
Przed zastosowaniem testu dla równości średnich w analizie wariancji konieczne jest sprawdzenie
hipotezy o równości wariancji zmiennej odpowiedzi w poszczególnych grupach. Stosujemy różne
testy do badania równości wariancji, jednak najlepszy jest test Levene'a.
Jego założenia:
´2 A- wariancja Å›r odpowiedzi dla paszy A
´2B- wariancja Å›r odpowiedzi dla paszy B
´2C- wariancja Å›r odpowiedzi dla paszy C
H : wszystkie wariancje sÄ… identyczne
0
vs. H : przynajmniej jedna wariancja różni się od pozostałych
1
Gdy p-wartość > ą=0,05 to przyjmujemy hipotezę zerową, że wszystkie wariancje są równe (są
homogeniczne), czyli można stosować analizę wariancji.
Co więcej, warto również zastosować inny test dla zbadania wartości wariancji np. test Bartlett'a.
Musimy wczeÅ›niej sprawdzić normalność skÅ‚adnika losowego µ za pomocÄ… testu Shapiro-Wilka.
Gdy p-wartość > ą=0,05 to składniki losowe modelu są normalne (możemy zastosować test
Bartlett'a).
Gdy po zastosowaniu testu Bartlett a (bartlett.test(przyrost~pasza) p-wartość > ą=0,05 oznacza to,
że test ten także potwierdził jednorodność wariancji.
Następnie po wprowadzeniu formuły: summary (model) otrzymujemy podsumowanie naszego
modelu, które możemy przedstawić w formie tabeli analizy wariancji:
y
ródło zmienności Liczba stopni Suma Średni Statystyka p-
swobody kwadratów kwadrat testowa F wartość
Czynniki
eksperymentalne
(rodzaj paszy)
Czynnik losowy
(naturalna zmienność
wewnÄ…trz grup)
Gdy p-wartość jest < Ä… = 0,05 to należy odrzucić hipotezÄ™ zerowÄ… (H ) o tym, że µA= µB= µC, czyli że
0
średnie przyrosty dla każdej z pasz są takie same. Oznacza to, że pasze różnicują wartości przyrostu.
Ponieważ odrzucamy H o równości średnich to można przeprowadzić test post-hoc. Stosujemy
0
metodę HSD Tukey'a dla układu zrównoważonego tzw. test Turkey o formule:
TukeyHSD(aov(przyrost~pasza)).
Gdy p adj, czyli wartości prawdopodobieństwa testowego są > ą = 0,05 określa to, że wszystkie pary
rodzajów paszy różnią się statystycznie istotnie w zakresie średnich przyrostów wagi, gdy tylko jedna
p adj > ą = 0,05 (np. dla pary C-B) to tylko dla pasz B i C można odrzucić hipotezę o równości średnich.
Możemy zilustrować przedziały ufności na wykresie korzystając z formuły: plot
(TukeyHSD(aov(przyrost~pasza))).