Narysuj tabelę 2 × 2 i oznacz a, b, c i d, zanim zaczniesz liczyć.
Czułość i swoistość odczytuj w kolumnach, a wartości predykcyjne wzdłuż wierszy.
Dla innej chorobowości wyobraź sobie 10 000 osób i wypełnij tabelę za pomocą Se, Sp i chorobowości.
Dla ilorazu wiarygodności przelicz prawdopodobieństwo na szanse, pomnóż i przelicz z powrotem.
W badaniu klinicznym najpierw znajdź EER i CER; następnie RR = EER/CER, ARR = CER − EER i NNT = 1/ARR zaokrąglone w górę.
Statystyka opisowa
Średnia
x̄ = Σx/n
Odchylenie standardowe z próby
s = √(Σ(x − x̄)²/(n − 1))
Mediana
środkowa wartość uporządkowanych danych
Wynik z
z = (x − μ)/σ
Rozkład normalny
68% w zakresie ±1σ, 95% w zakresie ±2σ
Przedziały ufności
σ znane
x̄ ± z·σ/√n
σ nieznane
x̄ ± t·s/√n, df = n − 1
Wartość krytyczna z
90%: 1,645, 95%: 1,96, 99%: 2,576
Liczebność próby
n = (z·σ/E)², rounded up
Metoda
Wybierz między z (znane σ populacji) a t (tylko s z próby, z n − 1 stopniami swobody).
Wyznacz wartość krytyczną dla poziomu ufności.
Margines błędu: E = wartość krytyczna × odchylenie standardowe/√n.
Przedział: x̄ ± E.
Testy hipotez i przedziały ufności
Test z dla średniej
z = (x̄ − μ0)/(σ/√n)
Test t dla średniej
t = (x̄ − μ0)/(s/√n), df = n − 1
Dwustronny na poziomie 5%
odrzuć H0, jeśli |z| > 1,96
Przedział ufności
x̄ ± z·σ/√n
Metoda
Sformułuj H₀ i H₁ i ustal, czy test jest jednostronny, czy dwustronny.
Oblicz statystykę testową: (x̄ − μ₀)/(σ/√n) albo z s i t.
Porównaj z wartością krytyczną albo wyznacz wartość p.
Odrzuć H₀, jeśli statystyka przekracza wartość krytyczną (p poniżej α), i sformułuj wniosek słowami.
2Zadania do ćwiczenia
1.Statystyka w praktyce opartej na dowodachŚredni
Ćwiczenie (wymyślone wartości, nie stanowią wskazówek klinicznych). Badanie kliniczno-kontrolne porównuje 110 pacjentów z chorobą (przypadki) i 70 osób bez niej (kontrole). Tabela pokazuje, ilu w każdej grupie było narażonych na czynnik ryzyka. Oblicz iloraz szans OR dla narażenia. Zaokrąglij do dwóch miejsc po przecinku.
a = 43b = 20c = 67d = 50
OR =
2.Statystyka w praktyce opartej na dowodachTrudny
Ćwiczenie (wymyślone wartości, nie stanowią wskazówek klinicznych). Pacjent ma przedtestowe prawdopodobieństwo choroby 20%. Test o czułości 66% i swoistości 79% daje wynik: dodatnia. Oblicz potestowe prawdopodobieństwo choroby p′ w procentach. Zaokrąglij do jednego miejsca po przecinku.
przedtestowe 20%Se = 66%Sp = 79%dodatnia
p′ = %
3.Statystyka opisowaŚredni
Wyznacz odchylenie standardowe z próby s dla danych: 45, 28, 20, 54, 18, 12.
n = 6
s =
4.Statystyka opisowaTrudny
Wielkość ma rozkład normalny ze średnią 100 i odchyleniem standardowym 5. Jaka część P wartości leży powyżej 97,5?
μ = 100σ = 5x = 97,5
P =
5.Przedziały ufnościŚredni
Próba 12 elementów ma średnią 72,5 i odchylenie standardowe s = 2,5. Populacja ma w przybliżeniu rozkład normalny. Wyznacz margines błędu E przedziału ufności 95% dla średniej.
n = 12x̄ = 72,5s = 2,595%t(11) = 2,201
E =
6.Przedziały ufnościTrudny
Z góry nic nie wiadomo o frakcji wadliwych produktów. Ile produktów n trzeba skontrolować, aby przedział ufności 90 % dla frakcji miał margines błędu co najwyżej 0,05?
p = 0,5E = 0,0590%
n =
7.Testy hipotez i przedziały ufnościŚredni
Twierdzi się, że średnia populacji wynosi 80; odchylenie standardowe populacji jest znane i wynosi 15. Losowa próba 36 obserwacji ma średnią 85,9. Przetestuj H₀: μ = 80 przeciw H₁: μ > 80. Wyznacz wartość p p.
μ₀ = 80σ = 15n = 36x̄ = 85,9H₁: μ > μ₀
p =
8.Testy hipotez i przedziały ufnościTrudny
W badaniu klinicznym wyzdrowiało 40 z 80 pacjentów leczonych metodą A i 117 z 300 pacjentów leczonych metodą B. Korzystając z frakcji łącznej, wyznacz statystykę testową z dla H₀: p₁ = p₂.
x₁ = 40n₁ = 80x₂ = 117n₂ = 300
z =
3Roztwory
1.Statystyka w praktyce opartej na dowodach
Odpowiedź: OR = 1,6
Ilorazy wiarygodności i miary ryzyka: LR+ = czułość/(1 − swoistość), LR− = (1 − czułość)/swoistość. W badaniu klinicznym EER i CER to częstości zdarzeń z leczeniem i bez leczenia: RR = EER/CER, ARR = CER − EER, RRR = ARR/CER i NNT = 1/ARR, zaokrąglone w górę. W badaniu kliniczno-kontrolnym można użyć wyłącznie ilorazu szans OR = ad/(bc). Współczynnik zapadalności = nowe przypadki/osobo-czas.
Szanse narażenia
przypadki: 43/67 = 0,6418
kontrole: 20/50 = 0,4
Iloraz szans
OR = (43·50)/(20·67) = 1,604 ≈ 1,6
Odpowiedź
OR = 1,6
OR powyżej 1: narażenie jest częstsze wśród przypadków.
2.Statystyka w praktyce opartej na dowodach
Odpowiedź: p′ = 44 %
Wartości predykcyjne zależą od chorobowości: PPV = Se·p/(Se·p + (1 − Sp)(1 − p)) i NPV = Sp(1 − p)/(Sp(1 − p) + (1 − Se)p). Z ilorazami wiarygodności: szanse przedtestowe = p/(1 − p), szanse potestowe = szanse przedtestowe × LR, prawdopodobieństwo = szanse/(1 + szanse). Dla populacji N prawdziwie dodatnie = N·p·Se, a fałszywie dodatnie = N·(1 − p)·(1 − Sp). NNT = 1/(CER × RRR), zaokrąglone w górę.
Iloraz wiarygodności
LR+ = 0,66/(1 − 0,79) = 3,143
Szanse przedtestowe
0,2/(1 − 0,2) = 0,25
Szanse potestowe
0,25·3,143 = 0,7857
Prawdopodobieństwo potestowe
p′ = 0,7857/(1 + 0,7857) = 0,44 = 44%
Odpowiedź
p′ = 44%
3.Statystyka opisowa
Odpowiedź: s = 16,56
Statystyka opisowa: średnia to suma podzielona przez n; mediana to środkowa wartość uporządkowanych danych (średnia dwóch środkowych, gdy n jest parzyste). Odchylenie standardowe z próby to s = √(Σ(x − x̄)²/(n − 1)). Wynik z = (x − μ)/σ mówi, o ile odchyleń standardowych wartość leży od średniej.
Statystyka opisowa: średnia to suma podzielona przez n; mediana to środkowa wartość uporządkowanych danych (średnia dwóch środkowych, gdy n jest parzyste). Odchylenie standardowe z próby to s = √(Σ(x − x̄)²/(n − 1)). Wynik z = (x − μ)/σ mówi, o ile odchyleń standardowych wartość leży od średniej.
Standaryzuj
z = −0,5
Tablica rozkładu normalnego
Φ(−0,5) = 0,3085
P = 1 − 0,3085 = 0,6915
Odpowiedź
P = 0,6915
5.Przedziały ufności
Odpowiedź: E = 1,588
Przedział ufności dla średniej: x̄ ± margines błędu. Przy znanym σ populacji margines wynosi z·σ/√n; gdy znasz tylko odchylenie standardowe z próby s, użyj t z n − 1 stopniami swobody: t·s/√n. Wartości krytyczne: z = 1,645 (90%), 1,96 (95%), 2,576 (99%). Aby uzyskać margines E lub mniejszy, potrzebujesz n ≥ (z·σ/E)², zaokrąglone w górę.
Wartość krytyczna
df = n − 1 = 11, t = 2,201
Margines błędu
E = 2,201·2,5/√12 = 1,588
Odpowiedź
E = 1,588
6.Przedziały ufności
Odpowiedź: n = 271
Confidence interval for a proportion: p̂ ± z·√(p̂(1 − p̂)/n), where p̂ = x/n. To reach a margin E you need n ≥ z²·p(1 − p)/E², rounded up; without a prior estimate use p = 0,5, which gives the largest n.
Wyznacz n
n = 1,645²·0,5·0,5/0,05² = 270,6
Zaokrąglij w górę
n = 271
Odpowiedź
n = 271
7.Testy hipotez i przedziały ufności
Odpowiedź: p = 0,009137
Hypothesis test: the test statistic = (estimate − value under H₀)/(standard error under H₀). The p-value is the probability, assuming H₀ is true, of a statistic at least as extreme as the one observed: the area in one tail for a one-sided H₁, twice the tail area for a two-sided H₁. Reject H₀ when p ≤ α.
Statystyka testowa
SE = σ/√n = 15/√36 = 2,5
z = (85,9 − 80)/2,5 = 2,36
Pole ogona
Φ(2,36) = 0,9909
1 − Φ(2,36) = 0,0091
Wartość p
p = 0,009137
Odpowiedź
p = 0,009137
p ≤ 0,05, więc H₀ odrzuca się na poziomie 5 %.
8.Testy hipotez i przedziały ufności
Odpowiedź: z = 1,775
Hypothesis test: the test statistic = (estimate − value under H₀)/(standard error under H₀). The p-value is the probability, assuming H₀ is true, of a statistic at least as extreme as the one observed: the area in one tail for a one-sided H₁, twice the tail area for a two-sided H₁. Reject H₀ when p ≤ α.
Frakcje z prób
p̂₁ = 40/80 = 0,5
p̂₂ = 117/300 = 0,39
p̂ = 157/380 = 0,4132
Błąd standardowy przy H₀
SE = √(0,4132·0,5868·(1/80 + 1/300)) = 0,06196
Statystyka testowa
z = (0,5 − 0,39)/0,06196 = 1,775
Odpowiedź
z = 1,775
|z| < 1.96, więc nie ma podstaw do odrzucenia H₀ na poziomie 5 % (test dwustronny).