Wykład 5: Testy parametryczne
Biometria i
Biostatystyka
Wykład 5: Testy parametryczne
Biometria i
Biostatystyka
Testy Istotności
Czemu służą?
Formułowanie hipotez
Statystyka testowa
P-wartości
Istotność/znamienność statystyczna
Testowanie wartości średniej w populacji
Dwustronne testy istotności i przedziały
ufności
P-wartości vs. ustalone α
Testy Istotności
2 najczęściej stosowane formy
wnioskowania statystycznego
Przedział ufności
Ocena parametru populacji
Test istotności
Ocena słuszności konkretnych stwierdzeń
o populacji dokonana na podstawie
dowodów dostarczonych przez dane.
Testy istotności jako
narzędzie wnioskowania
statystycznego
Test istotności
Formalna procedura porównywania
zebranych danych z hipotezą, której
prawdziwość chcemy ocenić.
Hipoteza: stwierdzenie o parametrach
populacji lub modelu
Wyniki testu są wyrażone w postaci
prawdopodobieństwa, które mierzy
jak hipoteza się zgadza z danymi
Testowanie hipotez
Hipoteza zerowa
Stwierdzenie będące przedmiotem
testu istotności
Test istotności jest tak skonstruowany,
by ocenić siłę dowodów przeciw niej
Stwierdzenie w postaci “brak wpływu”
lub “brak różnic”
Skrót: H
0
Testowanie hipotez, cd.
Hipoteza alternatywna
Stwierdzenie, które podejrzewamy o
prawdziwość jako alternatywę dla H
0
Skrót: H
a
Hipotezy zawsze odnoszą się do populacji lub
modelu, nie do poszczególnych wyników.
Hipotezy testów parametrycznych dotyczą
parametrów populacji.
Istnieją jednostronne lub dwustronne
hipotezy alternatywne.
Statystyki testowe
Test opiera się na statystyce szacującej parametr,
o którym mówi hipoteza. Zwykle jest to taka sama
estymata, którą wykorzystuje się przy
konstruowaniu przedziałów ufności dla badanego
parametru. Kiedy H
0
jest prawdziwa, oczekujemy,
że estymacja daje wartości bliskie tym
wynikającym z H
0
.
Wartości estymat dalekie od wartości parametru
opisanej przez H
0
dostarczają dowodu przeciw H
0
.
Hipoteza alternatywna determinuje, który kierunek
uważamy za sprzyjający H
a
.
Statystyki testowe, cd.
Statystyki testowe
Mierzą zgodność hipotezy zerowej i
danych.
Używane są do obliczania
prawdopodobieństwa, które
potrzebujemy do testu istotności.
Zmienne losowe ze znanym
rozkładem.
P-wartości
Test istotności znajduje
prawdopodobieństwo otrzymania
wyniku tak skrajnego lub bardziej,
niż w danej chwili obserwowany
.
Skrajny: daleki od tego, którego
byśmy się spodziewali gdyby hipoteza
zerowa była prawdziwa
P-wartości, cd.
Prawdopodobieństwo, obliczane przy
założeniu prawdziwości H
0
, że statystyka
testowa da wynik tak skrajny lub bardziej niż
w danej chwili obserwowany jest nazywane P-
wartością testu.
Im mniejsza P-wartość, tym mocniejszy
dowód przeciw H
0
dostarczony przez dane
Oblicza się ją używając funkcji rozkładu
gęstości prawdopodobieństwa statystyki
testowej (np. standardowego rozkładu
normalnego dla statystyki testowej z).
Istotność Statystyczna
Jeśli poziom P jest mniejsza lub równa
niż α, mówimy że dane są znamienne
statystycznie na poziomie istotności α.
Przykład:
poziom istotności 0.01 znaczy
że mamy dowód tak silny, że pojawiałby
się tylko w 1% obserwacji, gdyby
hipoteza zerowa była rzeczywiście
prawdziwa.
Co mówi P-wartość równa 0.03?
Krok 1. Korzystając z generatora liczb
losowych N(0,1) wygeneruj dwa
zbiory pomiarów po 16 pomiarów w
każdym
Krok 2. Zapamiętaj P-wartości
uzyskane z zastosowanego testu
statystycznego (np. testu t)
Liczba powtórzeń kroków 1-2:
N=25000
Przykład obliczeniowy
N = 1262
Test istotności - kroki
Można ocenić istotność dowodów
przeciw hipotezie zerowej dostarczonych
przez dane wykonując czynności:
1.
Sformułuj hipotezę zerową H
0
i hipotezę
alternatywną H
a
. Test jest przeznaczony do
oceny siły dowodów przeciw H
0
. H
a
jest
stwierdzeniem które zaakceptujemy jeśli
dowody pozwolą nam odrzucić H
0
.
2.
Oblicz wartość statystyki testowej. Ta
statystyka zwykle określa jak daleko dane są
od H
0
.
Test istotności - kroki
3.
Znajdź P-wartość dla obserwowanych danych. Jest
to prawdopodobieństwo (przyjmując, że hipoteza
zerowa jest prawdziwa), że statystyka będzie się
opowiadać przeciw hipotezie co najmniej tak silnie
jak z tymi danymi.
4.
Sformułuj wniosek. Wybierz poziom istotności α (jak
bardzo dowód przeciw H
0
uznajesz za decydujący).
Jeśli P-wartość jest mniejsza lub równa α - wniosek,
że hipoteza alternatywna jest prawdziwa; jeśli jest
większa od α - wniosek, że dane nie dostarczają
wystarczająco silnego dowodu do odrzucenia H
0
.
Twój wniosek jest podsumowaniem badań
wykonanych za pomocą testu istotności.
Testy dla wartości średniej w
populacji
Mamy prostą próbę losową o liczności n
wylosowanej z normalnej populacji o
nieznananej średniej μ ale ze znaną wariancją
σ
2
. Chcemy sprawdzić hipotezę że μ ma
określoną wartość, np. μ
0
.
Hipoteza zerowa: H
0
: μ = μ
0
.
Statystyka testowa:
jeśli hipoteza zerowa jest prawdziwa, to statystyka z
ma standardowy rozkład normalny.
n
x
z
/
0
Testy dla wartości średniej w
populacji
Alternatywna hipoteza: Jeśli jest
jednostronna to np. H
a
: μ > μ
0
P-wartość: prawdopodobieństwo, że
standardowa normalna zmienna
losowa Z przybierze wartość co
najmniej taką jak obserwowana z.
Możemy podobnie wnioskować dla
innych hipotez alternatywnych.
Dwustronne testy
istotności i przedziały
ufności
Poziom α dwustronnego testu
istotności odrzuca hipotezę H
0
: μ =
μ
0
dokładnie wtedy, gdy wartość μ
0
znajduje się poza przedziałem
ufności dla dla μ obliczonym dla
poziomu 1-α .
P-wartości vs. ustalone α
P-wartość to najmniejszy poziom α, na
którym dane są istotne.
Znając P-wartość możemy ocenić istotność
na każdym poziomie.
To daje więcej informacji niż sprawdzanie
odrzucony-lub-nie na ustalonym poziomie
istotności.
Wartość z*, taka że P{Z>z
*
} jest równa
zadanej liczbie a, 0<a<1, jest nazywana
wartością krytyczną standardowego
rozkładu normalnego.
Wykorzystywanie i
nadużywanie testów
Wybieranie poziomu istotności
Wybierz poziom α z góry jeśli musisz
dokonać decyzji.
Nie ma sensu, jeśli chcesz jedynie opisać siłę
swoich dowodów.
Jeśli stosujesz test istotności z ustalonym α
żeby podjąć decyzję, wybierz α pytając jak
silny dowód jest potrzebny żeby odrzucić H
0
.
To też zależy od tego jak wiarygodna jest hipoteza
zerowa.
Wybór poziomu istotności
Jeśli H
0
reprezentuje hipotezę, w
którą każdy wierzył przez lata, będzie
potrzebny mocny dowód (małe α)
żeby je obalić.
Siła dowodu potrzebnego do
odrzucenia H
0
zależy od
konsekwencji podjęcia takiej decyzji.
Kosztowne: silny dowód
Wybór poziomu istotności
Lepiej podawać P-wartość, która
pozwala każdemu z nas decydować
indywidualnie czy mamy
wystarczająco silne dowody.
Nie ma ostrej granicy między
„istotny” a „nieistotny” a jedynie
rosnąca siła dowodu przy malejącej
P-wartości.
Czego istotność
statystyczna nie oznacza?
Istotność statystyczna to nie to samo co
praktyczna istotność.
Przykład: hipoteza bez zależności jest odrzucana
Nie oznacza mocnego związku, ale że jest silny dowód
na istnienie jakiegoś związku
Kilka punktów odstających może się
przyczynić do zaobserwowania wysokiej
istotności wyniku testu, jeśli ślepo stosuje
się testy istotności.
Punkty odstające mogą również zniwelować
istotność.
Nie ignoruj braku
istotności
Jeśli badacz ma wyraźny powód
podejrzewać, że zależność istnieje a potem
nie może znaleźć istotnego dowodu nań, to
może to być ciekawą wiadomością ---
czasami ciekawszą niż to gdyby
potwierdzono związek na poziomie
istotności 5%.
Ukrywanie negatywnych wyników może
skazać innych badaczy na poszukiwanie
zależności, która nie istnieje.
Wnioskowanie statystyczne nie
jest słuszne dla wszystkich
danych
Formalne wnioskowanie
statystyczne nie koryguje błędnie
przeprowadzonego eksperymentu.
Losowość próbkowania gwarantuje,
że prawa probabilistyki mają
zastosowanie w naszych testach
istotności oraz wyznaczanych
przedziałach ufności.
Nie szukaj na siłę
znamienności
statystycznej!
Można także wnioskować bez
uwzględniania istotności
statystycznej.
Jeśli zdecydujesz jakiego wyniku
szukasz, zaprojektuj eksperyment
lub próbę, która Cię do niego
doprowadzi i zastosuj test istotności
żeby ocenić wagę dowodów.
Nie szukaj na siłę
znamienności
statystycznej!
Ponieważ udane próby szukania
zjawisk naukowych zwykle kończą
się znalezieniem istotności
statystycznej, uczynienie samej
istotności obiektem poszukiwań
jest zbyt nęcące.
Nie szukaj na siłę
znamienności
statystycznej!
Jeśli już masz hipotezę, zaprojektuj
badanie tak, żeby otrzymać wynik,
który uważasz za istniejący.
Jeśli wynik badania jest
statystycznie istotny, masz już
rzeczywisty dowód.
Moc i wnioskowanie
Badanie użyteczności przedziału
ufności
Poziom ufności: mówi nam jak
niezawodna jest ta metoda przy
wielokrotnych powtórzeniach
eksperymentu.
Margines błędu: mówi nam, jak czuła
jest ta metoda lub jak bardzo przedział
ogranicza szacowanie parametru
Moc i wnioskowanie
Badanie użyteczności testów istotności
przy ustalonym α.
Poziom istotności: mówi nam jak
wiarygodna jest ta metoda w użyciu
Moc testu
: mówi nam o zdolności testu do
wykrywania tego, że hipoteza zerowa jest
fałszywa
Mierzona jako prawdopodobieństwo że test
odrzuci hipotezę zerową kiedy alternatywna jest
prawdziwa.
Im wyższe prawdopodobieństwo, tym
bardziej czuły jest test.
Moc testu
Prawdopodobieństwo, że test
istotności przy ustalonym α odrzuci
H
0
, kiedy alternatywna wartość
parametru jest prawdziwa, jest
nazywane mocą testu do
wykrywania tej alternatywy.
Obliczanie mocy testu
Sformułuj H
0
, H
a
(konkretną
alternatywę, którą chcemy wykryć) i
poziom istotności α.
Znajdź wartości , które
spowodują że odrzucimy H
0
.
Oblicz prawdopodobieństwo
zaobserwowania tych wartości , dla
których alternatywa jest prawdziwa.
x
x
Zwiększanie mocy testu
Zwiększ α. 5%-owy test istotności
będzie miał większą szansę odrzucenia
alternatywy niż 1%-owy, ponieważ jest
wymagana mniejsza siła dowodu.
Weź taką alternatywę, która jest dalej
od μ
0
. Wartości μ, które są w H
a
ale
leżą blisko do hipotetycznej wartości μ
0
są cięższe do wykrycia (mniejsza moc)
niż wartości μ, które są daleko od μ
0
.
Zwiększanie mocy testu,
cd.
Zwiększ rozmiar próby. Więcej danych
dostarczy więcej informacji o więc jest
większa szansa odróżnienia wartości μ.
Zmniejsz σ. To daje taki sam efekt jak
zwiększanie rozmiaru próby: więcej
informacji o μ. Poprawienie procesu
pomiarów i ograniczenie uwagi na
subpopulacje to dwa najpopularniejsze
sposoby na zmniejszenie σ.
x
Dwa typy błędów
Podczas przeprowadzania testów istotności
musimy przyjąć jedną hipotezę a drugą
odrzucić.
Mamy nadzieję na trafność decyzji, ale nie
zawsze jest to dobry wybór.
2 typy niewłaściwych decyzji:
Jeśli odrzucimy H
0
(przyjmiemy H
a
) kiedy w
rzeczywistości H
0
jest prawdziwe, jest to błąd
pierwszego rodzaju (odniesiony do p-wartości).
Jeśli przyjmiemy H
0
(odrzucimy H
a
) kiedy w
rzeczywistości H
a
jest prawdziwe, jest to błąd
drugiego rodzaju (odniesiony do mocy).
Ogólnie przyjęte
postępowanie przy
testowaniu hipotez
Zdefiniuj H
0
i H
a
tak samo jak do testu
istotności.
Popatrz na problem jak na decyzję -
prawdopodobieństwa popełnienia błędów
I-szego i II-ego rodzaju są powiązane.
Błędy pierwszego rodzaju są poważniejsze.
Wybierz α (poziom istotności) i rozważ
testy tylko takie, gdzie
prawdopodobieństwo popełnienia błędu I-
szego rodzaju nie jest większe od α.
Ogólnie przyjęte
postępowanie przy
testowaniu hipotez
Spośród testów wybierz ten, który ma
jak najmniejsze prawdopodobieństwo
popełnienia błędu II-ego rodzaju (czyli
jak największą moc). Jeśli to
prawdopodobieństwo jest zbyt duże,
będziesz musisz wziąć większą liczbę
prób żeby zmniejszyć ryzyko błędu.
Jednostronny test t dla
jednej próby
Niech prosta próba losowa (PPL) o liczności n
jest losowana z populacji o nieznanej wart.
oczekiwanej μ. Żeby zweryfikować hipotezę
że H
0
: μ = μ
0
na podstawie PPL, oblicz
statystykę t
Zmienna losowa T ma rozkład t(n-1), P-
wartość dla testu H
0
przeciw H
a
: μ > μ
0
wynosi a przeciw H
a
: μ < μ
0
wynosi
n
/
s
x
t
0
t
T
P
)
t
T
(
P
Test t - przykład 1
Niech PPL o liczności n jest losowana z
populacji o nieznanej wartości
oczekiwanej μ.
[114, 123.3, 116.7, 129.0, 118, 124.6, 123.1, 117.4, 111,
121.7, 124.5, 130.5]
średnia próby = 121.15
odchylenie standardowe próby = 5.89
Żeby zweryfikować hipotezę że H
0
: μ =
μ
0
=120 na podstawie PPL liczności n,
oblicz statystykę t dla jednej próby
6764
.
0
12
/
89
.
5
120
15
.
121
n
/
s
x
t
0
Test t - przykład 1
W kategoriach zmiennej losowej T
z rozkładem t(n-1), P-wartość dla
testu H
0
przeciw H
a
: μ > μ
0
wynosi
t
T
P
stopień swobody
Wartość krytyczna dla
α=0.95
Jeśli t
1
=0.6 to cdf(t
1
)=0.71967 więc p
1
=1-
0.71967=0.28033
Jeśli t
2
=0.7 to cdf(t
2
)=0.75077 więc p
2
=1-
0.75077=0.24923
więc p
2
< p < p
1
Test t - przykład 1
Z rozkładu t Studenta
otrzymujemy
tcdf(0.6764,11) = 0.743621
p = 1 – 0.743621 =
0.256379
Test t - przykład 1
Test t - przykład 1
0.7436
p=0.2564
Test t - przykład 1
Końcowy wniosek:
Nie możemy odrzucić H
0
mówiącej, że wartość średnia w
populacji wynosi μ=120 i przyjąć
H
a
z p=0.2564
Kolejny przykład jednostronnego
testu t dla jednej próby
Żeby zweryfikować hipotezę H
0
: μ = 135
na podstawie PPL liczności n, obliczmy
statystykę t dla jednej próby:
W kategoriach zmiennej losowej T z
rozkładem t(n-1), P-wartość dla testu na
H
0
przeciw
H
a
: μ < μ
0
wynosi
n
/
s
135
x
t
t
T
P
Test t - przykład 2
[114, 123.3, 116.7, 129.0, 118, 124.6,
123.1, 117.4, 111, 121.7, 124.5, 130.5]
średnia próby = 121.15
odchylenie standardowe próby = 5.89
Żeby zweryfikować hipotezę że
H
0
: μ = μ
0
=135, obliczmy
statystykę t
1456
.
8
12
/
89
.
5
135
15
.
121
n
/
s
x
t
0
Z rozkładu t Studenta
otrzymujemy
P(T ≤ t)
=
tcdf(-8.1456,11) =2.7e-6
p =
0.0000027 < 0.000005
Test t - przykład 2
Test t - przykład 2
p=2.7e-6
t = -8.1456
Test t - przykład 2
Końcowy wniosek:
Odrzucamy H
0
mówiącą, że wartość
średnia w populacji wynosi μ=135 i
przyjmujemy H
a
: μ < μ
0
z p <
0.000005
Dwustronny test t dla jednej
próby
Żeby zweryfikować hipotezę że H
0
: μ =
115
na podstawie PPL liczności n,
obliczmy statystykę t
W kategoriach zmiennej losowej T z
rozkładem t(n-1), P-wartość dla testu H
0
przeciw
H
a
: μ ≠ μ
0
wynosi
n
/
s
115
x
t
|
t
|
|
T
|
P
Test t - przykład 3
[114, 123.3, 116.7, 129.0, 118, 124.6,
123.1, 117.4, 111, 121.7, 124.5, 130.5]
średnia próby = 121.15
odchylenie standardowe próby = 5.89
Żeby zweryfikować hipotezę że
H
0
: μ = μ
0
=115 obliczmy
statystykę t
6170
.
3
12
/
89
.
5
115
15
.
121
n
/
s
x
t
0
• Ponieważ wartość krytyczna
t
0.05
=2.2010 a nasza obserwowana
wartość jest wyższa, odrzucamy
hipotezę zerową na poziomie α=0.05.
• Z rozkładu t otrzymujemy
P(|T| > |t|)
= 2*(1-
tcdf(3.6170,11)) =
0.0040
Więc dokładnie p =
0.0040
Test t - przykład 3
Test t - przykład 3
p=0.0040
-t = -3.6170
t = 3.6170
Test t - przykład 3
Końcowy wniosek:
odrzucamy H
0
że wartość średnia w
populacji wynosi μ=115 i
przyjmujemy H
a
: μ ≠ μ
0
z p =
0.0040
Pary obserwacji - test t
W badaniu par obserwacji wyniki
są łączone w pary i porównywane
w jej obrębie.
Przykład: wyniki przed i po kursie
Kluczowe punkty odnośnie
par obserwacji
•
Analiza par obserwacji jest
konieczna kiedy mamy dwa
pomiary lub obserwacje każdego
obiektu i chcemy zbadać zmianę
jednej względem drugiej.
Zazwyczaj obserwacje w pewnym
sensie są pomiarami „przed” i
„po”.
Kluczowe punkty odnośnie
łączonych par
•
W każdej parze odejmuje się
pomiar „przed” od pomiaru „po”.
•
Analizuje się rozkład różnic
stosując przedziały ufności i testy
istotności dla jednej próby.
Odporność procedur t
Wnioskowanie statystyczne jest
nazywane odpornym, jeśli wymagane
metodyka obliczeń jest nieczuła na
naruszenie przyjętych założeń.
Procedury t są dość odporne na
odstępstwa od normalności rozkładu
populacji z wyjątkiem obserwowania
silnej skośności rozkłądu lub
występowania punktów odstających.
Odporność procedur t
Większe próby to większa dokładność
oszacowania P-wartości oraz wartości
krytycznych dystrybucji t, gdy rozkład
populacji nie jest normalny.
W przypadku małolicznych prób, zanim
przejdzie się do testów t należy
narysować wykres qq lub ramkowy w
celu sprawdzenia skośności i punktów
ostających.
Praktyczne wskazówki do
wnioskowania z jednej
średniej
Liczność próby: mniej niż 15: Zastosuj
procedury t jeśli dane mają w przybliżeniu
rozkład normalny. Jeśli ich rozkład jest daleki
od normalnego, albo zaobserwowaliśmy punkty
odstające, nie używaj t.
Liczność próby co najmniej 15: Można
zastosować procedury t, chyba że istnieją
punkty odstające lub rozkład jest bardzo skośny
Duże próby: Można stosować procedury t,
nawet dla skośnych rozkładów, ale liczność
musi być duża: ponad 40
Porównywanie dwóch
średnich
Problemy z dwiema próbami
Cel wnioskowania: porównanie
odpowiedzi w dwóch grupach.
Przyjmuje się że każda grupa to próba
z oddzielnej populacji.
Odpowiedzi w każdej z grup są
niezależne.
Pojęcia
Populacja
Zmienna
Średnia
Odchylenie standardowe
Liczność próby
Średnia próby
Odchylenie standardowe próby
Statystyka z dla dwóch
prób
Naturalny estymator różnicy μ
1
- μ
2
jest
różnicą dwóch średnich prób,
Żeby wnioskować z tej statystyki
musimy znać jej rozkład próbkowania.
Jeśli rozkłady obu populacji są normalne,
rozkład
jest także normalny.
2
1
x
x
2
1
x
x
2
2
2
1
2
1
2
1
,
n
n
N
Statystyka z dla dwóch
prób
Niech jest średnią prostej próby losowej o
liczności n
1
z populacji o rozkładzie N(μ
1,
σ
1
)
a to średnia PPL o liczności n
2
z populacji
o rozkładzie N(μ
2,
σ
2
). Wtedy taka statystyka
ma standardowy rozkład funkcji gęstości
prawdopodobieństwa N(0,1).
1
x
2
x
2
2
2
1
2
1
2
1
2
1
n
n
x
x
z
Procedura t dla dwóch
prób
Ten test zakłada że populacje, z
których wzięto próbki mają
równe wariancje
Procedura t dla dwóch
prób
Jeśli odchylenia standardowe
populacji nie są znane,
przybliżamy je na podstawie
odchylenia standardowego próby.
Stosuj właściwy wzór, zależny od
liczności prób: n
1
i n
2
.
Procedura t dla dwóch
prób
Kiedy liczności prób są różne i
przynajmniej jedno z n
1
i n
2
jest
małe (<30) stosuj:
Taka statystyka ma rozkład t Studenta
z k=n
1
+ n
2
– 2 stopniami swobody.
2
1
2
1
2
1
2
2
2
2
1
1
2
1
2
1
n
n
n
n
2
n
n
s
)
1
n
(
s
)
1
n
(
x
x
t
Procedura t dla dwóch
prób
Kiedy n
1
i n
2
są sobie równe
(niezależnie od liczności) zastosuj:
Taka statystyka ma rozkład t Studenta
z k=2(n – 1) stopniami swobody.
n
s
s
x
x
t
2
2
2
1
2
1
2
1
Procedura t dla dwóch
prób
Kiedy n
1
i n
2
są różne, ale duże (≥30)
stosuj:
Taka statystyka ma rozkład t Studenta
z k=n
1
+ n
2
- 2 stopniami swobody.
2
2
2
1
2
1
2
1
2
1
n
s
n
s
x
x
t
Test istotności t dla dwóch
prób
Niech prosta próba losowa o liczności n
1
jest losowana z populacji o rozkładzie
normalnym z nieznaną wartością
oczekiwaną μ
1
i jest niezależna od
drugiej próby o liczności n
2
, losowanej z
innej populacji o nieznanej wartości
oczekiwanej μ
2.
Zakładamy że odchylenia standardowe
obu populacji są równe.
Test istotności t dla dwóch
prób
Żeby zweryfikować hipotezę H
0
: μ
1
= μ
2
przeciw H
a
: μ
1
≠ μ
2
, obliczymy
odpowiednią statystykę t i użyjemy
P-wartości lub wartości
krytycznych dla rozkładu t
Studenta z odpowiednim stopniem
swobody.
Test istotności t dla
dwóch prób - przykład
Niech PPL o liczności n
1
jest losowana z
populacji o nieznanej wartości oczekiwanej μ
1
.
[114, 123.3, 116.7, 129.0, 118, 124.6, 123.1,
117.4, 111, 121.7, 124.5, 130.5]
średnia próby = 121.15
odchylenie standardowe próby = 5.89
Kolejna PPL o liczności n
2
jest losowana z
populacji o nieznanej wartości oczekiwanej μ
2
.
[120, 125.5, 126, 125.5, 128.5, 125, 128, 116,
122, 121, 117, 125]
średnia próby = 123.29
odchylenie standardowe próby = 4.08
Test istotności t dla
dwóch prób - przykład
0383
.
1
n
s
s
x
x
t
2
2
2
1
2
1
2
1
=0 z założenia
p = 0.3104
Wniosek: nie możemy odrzucić hipotezy, że wartości
średnie obu populacji są równe na poziomie istotności
α=0.05
Ponieważ n
1
=n
2
=12
Przedział ufności dla dwóch
prób
Niech PPL liczności n
1
jest losowana z
normalnej populacji o nieznanej wart.
oczekiwanej μ
1
i kolejna PPL o liczności n
2
jest
losowana z kolejnej normalnej populacji o
nieznanej wart. oczekiwanej μ
2
. Dla dużych n
1
i n
2
przedział ufności dla μ
1
- μ
2
jest dany przez
gdzie statystyka t ma n
1
+n
2
-2 stopni swobody.
2
2
2
1
2
1
2
1
2
2
2
1
2
1
2
1
*
,
*
n
s
n
s
t
x
x
n
s
n
s
t
x
x
Odporność procedur dla
dwóch prób
Procedury t dla dwóch prób są bardziej
odporne niż metody t dla jednej próby.
Jeśli liczności obu prób są równe a rozkłady
porównywanych populacji mają podobne
kształty, otrzymujemy dobrą dokładność.
Jeśli kształty rozkładów populacji są różne,
potrzeba prób o większej liczności.
Planując badanie oparte na dwóch
próbach, powinieneś zwykle wybierać
równe liczności tych prób.
Test Bartlett’a
homogeniczności
(jednorodności) wariancji
Testowanie hipotezy zerowej
mówiącej że dwie populacje
reprezentowane przez dwie próby
mają równe wariancje.
Hipoteza alternatywna mówi że
dwie wariancje są różne
Test jednorodności wariancji
Bartletta
Oblicz logarytm naturalny
wariancji każdej z prób.
Zsumuj liczby stopni swobody
Oblicz średnią ważoną wariancji i
jej logarytm naturalny
a
i
95
9
...
12
17
)
1
n
(
185167
.
2
112459
.
0
ln
,
112459
.
0
)
1
n
(
s
)
1
n
(
s
a
i
a
2
i
i
2
Test jednorodności wariancji
Bartletta
Oblicz ważoną sumę logarytmów
wariancji każdej z prób
Oblicz statystykę
227545
.
229
s
ln
)
1
n
(
a
2
i
i
636721
.
21
s
ln
)
1
n
(
s
ln
)
1
n
(
X
a
2
i
i
2
a
i
2
Test jednorodności wariancji
Bartletta
Jeśli hipoteza zerowa jest
prawdziwa, to skorygowana
statystyka X
2
ma w przybliżeniu
rozkład χ
2
o a-1 stopniach
swobody.
914
.
20
034566
.
1
63672
.
21
C
X
X
034566
.
1
)
1
n
(
1
1
n
1
)
1
a
(
3
1
1
C
2
2
kor
a
a
i
i
Test jednorodności wariancji
Bartletta
475
.
18
2
]
7
[
01
.
0
Wartość krytyczna dla a-1=7 stopni
swobody oraz alfa = 0.01 wynosi
Wartość skorygowanej statystyki wynosi
20.914. Wnioskujemy więc, że wariancje są
niejednorodne.
Przybliżony test t Welch’a
równości wartości średnich
dwóch populacji o różnych
wariancjach
Przybliżony test t Welch’a
Ten test oblicza przybliżoną t-
wartość, taką
’
, dla której krytyczna
wartość jest liczona jako średnia
ważona pojedynczych wartości
krytycznych t odpowiadających
stopniom swobody dwóch prób.
Przybliżony test t Welch’a
Dwupróbowa statystyka t
’
jest
równa:
2
2
2
1
2
1
2
1
2
1
n
s
n
s
x
x
't
Przybliżony test t Welch’a
Wartość krytyczna t
α’
dla błędu
pierwszego rodzaju jest liczona
jako:
2
2
2
1
2
1
2
2
2
]
[
1
2
1
]
[
n
s
n
s
n
s
t
n
s
t
't
2
1
Przybliżony test t Welch’a
Ta statystyka NIE MA rozkładu t ale
możemy przybliżyć rozkład t’ dla
dwóch prób stosując rozkład t(k) z
przybliżeniem liczby stopni
swobody k.
k równa się mniejszą z dwóch liczb:
n
1
–1 oraz n
2
–1.
Zadanie domowe