Statystyka 2: Wnioskowanie statystyczne: zestaw na powtórkę

Wzory, najważniejsze kroki metod dla każdego tematu i 10 zadań do ćwiczenia, z pełnymi rozwiązaniami na końcu.

Nowe zadaniaWróć do przedmiotu

W oknie drukowania wybierz „Zapisz jako PDF”, aby zachować kopię na swoim urządzeniu.

1Wzory i metody

Rozkład normalny

Wynik z
z = (x − μ)/σ
Poniżej x
P(X < x) = Φ(z)
Powyżej x
P(X > x) = 1 − Φ(z)
Pomiędzy
Φ(z2) − Φ(z1)
Percentyl
x = μ + z·σ

Metoda

  1. Zamień wartość na wynik standaryzowany: z = (x − μ)/σ.
  2. Odczytaj Φ(z), pole na lewo od z.
  3. Dla pola na prawo użyj 1 − Φ(z); dla przedziału odejmij dwa pola lewostronne.
  4. Dla percentyla licz wstecz: wyznacz z dla pola, a potem x = μ + z·σ.

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

  1. Wybierz między z (znane σ populacji) a t (tylko s z próby, z n − 1 stopniami swobody).
  2. Wyznacz wartość krytyczną dla poziomu ufności.
  3. Margines błędu: E = wartość krytyczna × odchylenie standardowe/√n.
  4. 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

  1. Sformułuj H₀ i H₁ i ustal, czy test jest jednostronny, czy dwustronny.
  2. Oblicz statystykę testową: (x̄ − μ₀)/(σ/√n) albo z s i t.
  3. Porównaj z wartością krytyczną albo wyznacz wartość p.
  4. Odrzuć H₀, jeśli statystyka przekracza wartość krytyczną (p poniżej α), i sformułuj wniosek słowami.

Regresja liniowa i korelacja

Sumy kwadratów
Sxx = Σ(x − x̄)² = Σx² − (Σx)²/n, Sxy = Σxy − ΣxΣy/n
Współczynnik kierunkowy
b = Sxy/Sxx
Wyraz wolny
a = ȳ − b·x̄
Korelacja
r = Sxy/√(SxxSyy)
Reszta
e = y − ŷ

Testy chi-kwadrat

Statystyka testowa
χ² = Σ(O − E)²/E
Test zgodności
E = n·p, df = k − 1 − (estimated parameters)
Niezależność
E = row total × column total / grand total, df = (r − 1)(c − 1)
Wartości krytyczne (α = 0,05)
df 1: 3,841, 2: 5,991, 3: 7,815, 4: 9,488
Decyzja
odrzuć H0, jeśli χ² > χ²crit

2Zadania do ćwiczenia

1.Rozkład normalnyŚredni
32 min35,8 minX + Y?

X ma rozkład normalny ze średnią 12 min i odchyleniem standardowym 1,5 min; zmienna Y jest niezależna i ma rozkład normalny ze średnią 20 min i odchyleniem standardowym 2 min. Wyznacz prawdopodobieństwo P, że X + Y jest większe niż 35,8 min. Zaokrąglij z do dwóch miejsc po przecinku.

μX = 12, σX = 1,5μY = 20, σY = 2X + Y > 35,8 min
P =
2.Rozkład normalnyTrudny
μ = 500x = ?99%

Rozkład czasów pracy baterii jest normalny ze średnią 500 h i odchyleniem standardowym 40 h. Wyznacz wartość x, poniżej której leży 99% wszystkich wartości (99. percentyl).

μ = 500σ = 40P = 99%
x = h
3.Przedziały ufnościŚredni
n = 6, s = 15,5L = ?Ux̄ = 10,5E95%

Próba 6 elementów ma średnią 10,5 i odchylenie standardowe s = 15,5. Populacja ma w przybliżeniu rozkład normalny. Wyznacz dolną granicę przedziału ufności 95% dla średniej.

n = 6x̄ = 10,5s = 15,595%t(5) = 2,571
dolna granica =
4.Przedziały ufnościTrudny
σ = 28, n = ?x̄ − Ex̄ + Ex̄E = 1,595%

Odchylenie standardowe populacji wynosi 28. Jak liczna musi być próba n, aby przedział ufności 95% dla średniej miał margines błędu co najwyżej 1,5?

σ = 28E = 1,595%
n =
5.Testy hipotez i przedziały ufnościŚredni
t −2,0302,03

Próba 36 studentów ma średni wynik 92,25 i odchylenie standardowe z próby 15. Wyznacz statystykę testową t dla H₀: μ = 100.

μ₀ = 100x̄ = 92,25s = 15n = 36
t =
6.Testy hipotez i przedziały ufnościTrudny
−1,961,96z = ?

W badaniu klinicznym wyzdrowiało 60 z 200 pacjentów leczonych metodą A i 18 z 80 pacjentów leczonych metodą B. Korzystając z frakcji łącznej, wyznacz statystykę testową z dla H₀: p₁ = p₂.

x₁ = 60n₁ = 200x₂ = 18n₂ = 80
z =
7.Regresja liniowa i korelacjaŚredni
xy139341544751954xy540455055

Tabela podaje 5 par obserwacji x = wydatki na reklamę w tys. € i y = sprzedaż w tys. €. Podane są ich sumy. Wyznacz z tych sum współczynnik korelacji Pearsona r.

n = 5Σx = 25Σy = 229Σx² = 165Σy² = 10655Σxy = 1225
r =
8.Regresja liniowa i korelacjaTrudny
xy1413485467579621165xy510405060ŷ = 48

Tabela podaje 6 par obserwacji x = wydatki na reklamę w tys. € i y = sprzedaż w tys. €. Użyj prostej najmniejszych kwadratów ŷ = a + b·x. Dla jakiej wartości x prosta przewiduje ŷ = 48?

n = 6ŷ = 48
x =
9.Testy chi-kwadratŚredni
poprawabrak poprawySumanowy lek9921120placebo5545100Suma15466220

Tabela kontyngencji pokazuje liczebności obserwowane. Oblicz statystykę chi-kwadrat χ² testu niezależności zmiennych wierszowej i kolumnowej. Nie stosuj poprawki na ciągłość.

N = 2202 × 2
χ² =
10.Testy chi-kwadratTrudny
dfα = 0,05α = 0,0113,8416,63525,9919,21037,81511,34549,48813,277511,07015,086612,59216,812714,06718,475815,50720,090916,91921,6661018,30723,209

Test niezależności chi-kwadrat przeprowadza się na tabeli kontyngencji o 2 wierszach i 4 kolumnach (bez sum). Korzystając z tabeli, wyznacz wartość krytyczną χ²crit na poziomie istotności α = 0,01.

α = 0,01
χ²crit =

3Roztwory

1.Rozkład normalny

Odpowiedź: P = 0,06426

Sumy i różnice zmiennych normalnych: jeśli X i Y są niezależne i mają rozkład normalny, to X + Y i X − Y również mają rozkład normalny. Średnie dodają się lub odejmują: E(X ± Y) = μX ± μY. Wariancje zawsze się dodają, także dla różnicy: Var(X ± Y) = σX² + σY², więc σ = √(σX² + σY²). Odchylenia standardowe nigdy nie dodają się bezpośrednio. Suma n niezależnych elementów o tym samym rozkładzie ma średnią n·μ i odchylenie standardowe σ·√n. Potem standaryzuj jak zwykle: z = (x − średnia)/σ i odczytaj Φ(z).

Rozkład sumy

μ = 32 min, σ = √(1,5² + 2²) = 2,5 min

Standaryzuj

z = (35,8 − 32)/2,5 = 1,52

Górny ogon

P = 1 − Φ(1,52) = 1 − 0,9357 = 0,06426

Odpowiedź

P = 0,06426
2.Rozkład normalny

Odpowiedź: x = 593,1 h

Rozkład normalny: zamień wartość x na wynik z, z = (x − μ)/σ, który mówi, o ile odchyleń standardowych wartość leży od średniej. Tablica rozkładu normalnego (czyli Φ(z)) podaje pole na lewo od z, czyli P(X < x). Pole na prawo to 1 − Φ(z), a pole między dwiema wartościami to różnica ich dwóch pól.

z dla percentyla

Φ(z) = 0,99 → z = 2,326

Powrót do x

x = 500 + 2,326·40 = 593,05 h

Odpowiedź

x = 593,05 h
3.Przedziały ufności

Odpowiedź: dolna granica = −5,769

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 = 5, t = 2,571

Margines błędu

E = 2,571·15,5/√6 = 16,27

Przedział

10,5 ± 16,27 → [−5,769, 26,77]

Odpowiedź

dolna granica = −5,769
4.Przedziały ufności

Odpowiedź: n = 1339

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ę.

Wyznacz n

n = (1,96·28/1,5)² = 1339

Zaokrąglij w górę

Liczebność próby jest liczbą całkowitą, a zaokrąglenie w dół dałoby za duży margines.

n = 1339

Odpowiedź

n = 1339
5.Testy hipotez i przedziały ufności

Odpowiedź: t = −3,1

Test dla średniej: statystyka testowa mierzy, o ile błędów standardowych średnia z próby leży od wartości z H₀: z = (x̄ − μ₀)/(σ/√n), gdy σ jest znane, t = (x̄ − μ₀)/(s/√n), gdy szacuje się je z próby. Przedział ufności to x̄ ± (wartość krytyczna)·(błąd standardowy).

Błąd standardowy

SE = 15/√36 = 15/6 = 2,5

Statystyka testowa

t = (92,25 − 100)/2,5 = −3,1

Odpowiedź

t = −3,1
t0,025; df = 35 = 2,03

|t| > 2,03, więc H₀ odrzuca się na poziomie 5 % (test dwustronny).

6.Testy hipotez i przedziały ufności

Odpowiedź: z = 1,265

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̂₁ = 60/200 = 0,3
p̂₂ = 18/80 = 0,225
p̂ = 78/280 = 0,2786

Błąd standardowy przy H₀

SE = √(0,2786·0,7214·(1/200 + 1/80)) = 0,0593

Statystyka testowa

z = (0,3 − 0,225)/0,0593 = 1,265

Odpowiedź

z = 1,265

|z| < 1.96, więc nie ma podstaw do odrzucenia H₀ na poziomie 5 % (test dwustronny).

7.Regresja liniowa i korelacja

Odpowiedź: r = 0,9794

Regresja metodą najmniejszych kwadratów: przy średnich x̄ i ȳ oraz Sxx = Σ(x − x̄)², Syy = Σ(y − ȳ)² i Sxy = Σ(x − x̄)(y − ȳ) prosta ŷ = a + b·x ma współczynnik kierunkowy b = Sxy/Sxx i wyraz wolny a = ȳ − b·x̄, więc zawsze przechodzi przez punkt (x̄, ȳ). Współczynnik korelacji Pearsona r = Sxy/√(Sxx·Syy) leży między −1 a 1; r² to część zmienności y wyjaśniona przez prostą.

Sumy scentrowane

Sxy = 1225 − 25·229/5 = 80
Sxx = 165 − 25²/5 = 40
Syy = 10655 − 229²/5 = 166,8

Współczynnik korelacji

r = 80/√(40·166,8) = 0,9794

Odpowiedź

r = 0,9794
8.Regresja liniowa i korelacja

Odpowiedź: x = 3,909

Regresja metodą najmniejszych kwadratów: przy średnich x̄ i ȳ oraz Sxx = Σ(x − x̄)², Syy = Σ(y − ȳ)² i Sxy = Σ(x − x̄)(y − ȳ) prosta ŷ = a + b·x ma współczynnik kierunkowy b = Sxy/Sxx i wyraz wolny a = ȳ − b·x̄, więc zawsze przechodzi przez punkt (x̄, ȳ). Współczynnik korelacji Pearsona r = Sxy/√(Sxx·Syy) leży między −1 a 1; r² to część zmienności y wyjaśniona przez prostą.

Średnie

x̄ = 36/6 = 6
ȳ = 319/6 = 53,17

Sumy kwadratów i iloczynów

Sxx = Σ(x − x̄)² = 70
Sxy = Σ(x − x̄)(y − ȳ) = 173

Prosta regresji

b = 173/70 = 2,471
a = 53,17 − 2,471·6 = 38,34
ŷ = 38,34 + 2,471·x

Wyznacz x

x = (48 − 38,34)/2,471 = 3,909

Odpowiedź

x = 3,909
9.Testy chi-kwadrat

Odpowiedź: χ² = 19,64

Testy chi-kwadrat: porównaj liczebności obserwowane O z liczebnościami E oczekiwanymi przy H₀ za pomocą χ² = Σ(O − E)²/E. Test zgodności: E = n·p dla każdej kategorii, df = k − 1 (minus jeden za każdy parametr estymowany z danych). Niezależność w tabeli r × c: E = suma wiersza × suma kolumny / suma ogólna, df = (r − 1)(c − 1). Odrzuć H₀, gdy χ² przekracza wartość krytyczną dla wybranego α.

Liczebności oczekiwane

nowy lek: 120·154/220 = 84, 120·66/220 = 36
placebo: 100·154/220 = 70, 100·66/220 = 30

Statystyka chi-kwadrat

χ² = (99 − 84)²/84 + (21 − 36)²/36 + (55 − 70)²/70 + (45 − 30)²/30
χ² = 19,64

Odpowiedź

χ² = 19,64

df = 1. Przy α = 0,05 wartość krytyczna wynosi 3,841, więc decyzja dotycząca H₀ (niezależności): odrzucona.

10.Testy chi-kwadrat

Odpowiedź: χ²crit = 11,35

Testy chi-kwadrat: porównaj liczebności obserwowane O z liczebnościami E oczekiwanymi przy H₀ za pomocą χ² = Σ(O − E)²/E. Test zgodności: E = n·p dla każdej kategorii, df = k − 1 (minus jeden za każdy parametr estymowany z danych). Niezależność w tabeli r × c: E = suma wiersza × suma kolumny / suma ogólna, df = (r − 1)(c − 1). Odrzuć H₀, gdy χ² przekracza wartość krytyczną dla wybranego α.

Stopnie swobody

df = (2 − 1)·(4 − 1) = 3

Tabela

χ²crit(df = 3, α = 0,01) = 11,345

Odpowiedź

χ²crit = 11,345

Odrzuć H₀, jeśli obliczone χ² jest większe od tej wartości.

ExamProblems · examproblems.com/courses/statistics-2-inference