Statistika

Příklad 10: Rozsah výběru pro podíl

Vzorcovník zatím není na e-learningu. Vyučující slíbili, že ho brzy nahrají a že bude v učebně i při testu. Až bude venku, upravím ruční výpočty tak, aby odkazovaly na vzorce z něj.

Zadání

Strana 10 v PDF

Stanovte potřebný rozsah výběru pro šetření, jehož předmětem je zjištění, zda se obyvatelé Šluknova zamýšleli nad tím, že by se z města vystěhovali. Z předvýzkumu bylo zjištěno, že podíl těch, kteří o tom uvažovali již někdy, je 20 %. Odhady chceme provádět se spolehlivostí 1-α = 0,95 a připouštíme chybu ve výši 2 %.

Zadání příkladu, pokračování:

Jak by se požadovaný rozsah výběru změnil, kdybychom byli ochotni připustit chybu 4 %?

Zadání příkladu, pokračování:

Jak by se požadovaný rozsah výběru změnil, kdybychom chtěli odhadnout maximální podíl těch, kteří uvažovali o vystěhování se ze Šluknova?

Co známeHodnota
podíl z předvýzkumu p0,2
spolehlivost 1 − α0,95
přípustná chyba Δ0,02 (1. pokračování: 0,04)

Typ proměnné a co z toho plyne

Postup s tahákem

Mít po ruce: tahák se vzorci, str. 1. Vzorec i kvantily u jsou na taháku. Procedura Sample Size Determination ve Statgraphicsu dá jiná čísla (1 671 místo 1 537), do testu se píše výpočet podle vzorce.

1. Kvantily

Tahák: Hladiny na Z rozdělení, hladiny spolehlivosti

oboustranný: u0,975 = 1,960     maximální podíl (jednostranný): u0,95 = 1,645

2. Dosadit

Tahák: Odhad minimálního rozsahu výběru, podle proporce

Na taháku je překlep: ve vzorci stojí p(p − 1), správně je p(1 − p). S p(p − 1) by vyšlo záporné číslo.

Δ = 2 %:   n ≥ u²0,975 · p(1 − p) / Δ² = 1,96² · 0,2 · 0,8 / 0,02² = 1 536,64 → 1 537
Δ = 4 %:   n ≥ 1,96² · 0,16 / 0,04² = 384,16 → 385
maximální podíl, Δ = 2 %:   n ≥ 1,645² · 0,16 / 0,02² = 1 082,41 → 1 083

Δ dosadit jako desetinné číslo (2 % = 0,02), stejně jako p. Rozsah se zaokrouhluje vždy nahoru.

3. Kontrola v Sample Size Determination (nepovinné)

Tools > Sample Size Determination > One Sample, Binomial Proportion, 0,2. V Options Absolute Error = 0,02 (nebo 0,04), Confidence level 95 %. Pro maximální podíl Alternative Hypothesis = Less Than.

Options: Absolute Error 0,02, Not Equal
Options: Absolute Error 0,02, Not Equal
Statgraphics: n = 1671. Ručně vychází 1 537
Statgraphics: n = 1 671. Podle vzorce vychází 1 537

Katedra: výsledek se liší, protože „zde jsou používány kvantily rozdělení F“. Na testu platí výpočet podle vzorce.

Ruční řešení po krocích

a) Základní případ: Δ = 2 %, oboustranný odhad

Krok 1: vzorec (RUČNĚ).
n ≥ u21−α/2 · p(1 − p) / Δ2
Krok 2: kvantil (RUČNĚ z tabulek).
u0,975 = 1,96
Krok 3: dosazení (RUČNĚ).
n ≥ 1,962 · 0,2 · 0,8 / 0,022 = 3,8416 · 0,16 / 0,0004 = 1 536,64

n = 1 537 respondentů

b) 1. pokračování: chyba 4 %

Krok 4 (RUČNĚ).
n ≥ 1,962 · 0,16 / 0,042 = 0,614656 / 0,0016 = 384,16

n = 385 respondentů, zhruba čtvrtina původního rozsahu.

Chyba je ve jmenovateli na druhou: dvojnásobná chyba = čtvrtinový rozsah (1 536,64 / 4 = 384,16).

c) 2. pokračování: maximální podíl (jednostranný odhad), Δ = 2 %

Krok 5 (RUČNĚ). Maximální podíl = pravostranný interval, kvantil u1−α. Vzorec podle přednášky a řešeného příkladu 13_5:
n ≥ u21−α · p(1 − p) / Δ2
u0,95 = 1,645
n ≥ 1,6452 · 0,16 / 0,0004 = 2,706025 · 0,16 / 0,0004 = 1 082,41

n = 1 083 respondentů, o 454 méně než u oboustranného odhadu.

Jednostranný odhad potřebuje menší výběr, protože u0,95 = 1,645 < u0,975 = 1,96.
Zadání se ptá „v prvním případě“ jen u příkladu 9. Tady není řečeno, jestli maximální podíl počítat s chybou 2 %, nebo 4 %. Hlavní výsledek je pro původní 2 %. Kdyby se myslelo 4 %: 1,645² · 0,16 / 0,0016 = 270,6 → 271 respondentů.
Kvantil u0,95 se v materiálech katedry zaokrouhluje různě: v 13_5 (rozsah výběru) 1,645, v 13_1 a 13_3 (intervaly) 1,64. S 1,64 by vyšlo 1 075,84 → 1 076. Tady je použito 1,645 jako v 13_5, které řeší přesně tenhle typ úlohy. S přesným kvantilem 1,64485 vyjde 1 082,22 → 1 083.
oboustranný, Δ = 2 % oboustranný, Δ = 4 % maximální podíl, Δ = 2 % 1 5373851 083 potřebný počet respondentů
Minimální rozsah výběru ve třech variantách příkladu 10.

Interpretace

Abychom podíl obyvatel, kteří uvažovali o vystěhování, odhadli s chybou nejvýše 2 procentní body a spolehlivostí 95 %, musíme se zeptat alespoň 1 537 lidí. Když připustíme chybu 4 procentní body, stačí 385 lidí. Pokud nám jde jen o maximální podíl (s chybou 2 %), stačí 1 083 lidí.

Na co se můžou zeptat v testu

Co dělat, když z předvýzkumu podíl neznáme?

Dosadit p = 0,5, kdy je p(1 − p) = 0,25 největší a rozsah vyjde „na jistotu“: 1,96² · 0,25 / 0,0004 = 2 401 respondentů. To je odhad podle obecného vzorce, v řešených příkladech katedry tato varianta není.

Jak se rozsah změní při spolehlivosti 99 % (Δ = 2 %)?

u0,995 = 2,576: 2,576² · 0,16 / 0,0004 = 2 654,3 → 2 655 respondentů. Vyšší spolehlivost potřebuje větší výběr.

Kolik by stačilo při výběru bez vracení z 6 500 obyvatel?

nBV = 6 500 · 1 537 / (6 500 + 1 537) = 1 243,1 → 1 244 respondentů. Tady je úspora velká, protože výběr tvoří skoro čtvrtinu populace.

Proč u podílu nikdy nepoužíváme t?

Vzorec pro π stojí na aproximaci normálním rozdělením (proto podmínka n · p(1 − p) > 9). Rozptyl se počítá přímo z p, nic se zvlášť neodhaduje ze směrodatné odchylky, takže t se tu nepoužívá.