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 min28,3 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ż 28,3 min. Zaokrąglij z do dwóch miejsc po przecinku.

μX = 12, σX = 1,5μY = 20, σY = 2X + Y > 28,3 min
P =
2.Rozkład normalnyTrudny
μ = 65x = ?90%

Rozkład wyników egzaminu jest normalny ze średnią 65 pkt i odchyleniem standardowym 10 pkt. Wyznacz wartość x, poniżej której leży 90% wszystkich wartości (90. percentyl).

μ = 65σ = 10P = 90%
x = pkt
3.Przedziały ufnościŚredni
n = 9, s = 8LU = ?x̄ = 51,5E99%

Próba 9 elementów ma średnią 51,5 i odchylenie standardowe s = 8. Populacja ma w przybliżeniu rozkład normalny. Wyznacz górną granicę przedziału ufności 99% dla średniej.

n = 9x̄ = 51,5s = 899%t(8) = 3,355
górna granica =
4.Przedziały ufnościTrudny
σ = 24, n = ?x̄ − Ex̄ + Ex̄E = 0,595%

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

σ = 24E = 0,595%
n =
5.Testy hipotez i przedziały ufnościŚredni
t −2,01102,011

Próba 49 studentów ma średni wynik 118,7 i odchylenie standardowe z próby 4. Wyznacz statystykę testową t dla H₀: μ = 120.

μ₀ = 120x̄ = 118,7s = 4n = 49
t =
6.Testy hipotez i przedziały ufnościTrudny
−1,961,96z = ?

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

x₁ = 105n₁ = 300x₂ = 48n₂ = 100
z =
7.Regresja liniowa i korelacjaŚredni
xy258368469578690790xy24660708090

Tabela podaje 6 par obserwacji x = liczba godzin nauki i y = wynik egzaminu. Wyznacz współczynnik korelacji Pearsona r.

n = 6
r =
8.Regresja liniowa i korelacjaTrudny
xy128222320419515614xy246152025x = 5,5

Tabela podaje 6 par obserwacji x = wiek samochodu w latach i y = cena w tys. €. Użyj prostej najmniejszych kwadratów ŷ = a + b·x. Wyznacz prognozę ŷ dla x = 5,5.

n = 6x = 5,5
ŷ =
9.Testy chi-kwadratŚredni
kaszelbrak kaszluSumapalący271340niepalący4357100Suma7070140

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 = 1402 × 2
χ² =
10.Testy chi-kwadratTrudny
zaliczoneniezaliczoneSumawiek 18–30332760wiek 31–50542680wiek 51+39140Suma12654180dfα = 0,05α = 0,0113,8416,63525,9919,21037,81511,34549,48813,277511,07015,086612,59216,812714,06718,475815,50720,090916,91921,6661018,30723,209

Tabela pokazuje liczebności obserwowane. Sprawdź na poziomie α = 0,05, czy obie zmienne są niezależne (bez poprawki na ciągłość). Oblicz statystykę testową χ², a potem porównaj ją z wartością krytyczną z tabeli.

N = 1803 × 2α = 0,05
χ² =

3Roztwory

1.Rozkład normalny

Odpowiedź: P = 0,9306

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 = (28,3 − 32)/2,5 = −1,48

Górny ogon

P = 1 − Φ(−1,48) = 1 − 0,06944 = 0,9306

Odpowiedź

P = 0,9306
2.Rozkład normalny

Odpowiedź: x = 77,82 pkt

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,9 → z = 1,282

Powrót do x

x = 65 + 1,282·10 = 77,82 pkt

Odpowiedź

x = 77,82 pkt
3.Przedziały ufności

Odpowiedź: górna granica = 60,45

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 = 8, t = 3,355

Margines błędu

E = 3,355·8/√9 = 8,947

Przedział

51,5 ± 8,947 → [42,55, 60,45]

Odpowiedź

górna granica = 60,45
4.Przedziały ufności

Odpowiedź: n = 8852

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·24/0,5)² = 8851

Zaokrąglij w górę

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

n = 8852

Odpowiedź

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

Odpowiedź: t = −2,275

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 = 4/√49 = 4/7 = 0,5714

Statystyka testowa

t = (118,7 − 120)/0,5714 = −2,275

Odpowiedź

t = −2,275
t0,025; df = 48 = 2,011

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

6.Testy hipotez i przedziały ufności

Odpowiedź: z = −2,317

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̂₁ = 105/300 = 0,35
p̂₂ = 48/100 = 0,48
p̂ = 153/400 = 0,3825

Błąd standardowy przy H₀

SE = √(0,3825·0,6175·(1/300 + 1/100)) = 0,05612

Statystyka testowa

z = (0,35 − 0,48)/0,05612 = −2,317

Odpowiedź

z = −2,317

|z| > 1,96, więc H₀ odrzuca się na poziomie 5 % (test dwustronny).

7.Regresja liniowa i korelacja

Odpowiedź: r = 0,9741

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̄ = 27/6 = 4,5
ȳ = 453/6 = 75,5

Sumy kwadratów i iloczynów

Sxx = Σ(x − x̄)² = 17,5
Sxy = Σ(x − x̄)(y − ȳ) = 117,5
Syy = Σ(y − ȳ)² = 831,5

Współczynnik korelacji

r = 117,5/√(17,5·831,5) = 0,9741

Odpowiedź

r = 0,9741
8.Regresja liniowa i korelacja

Odpowiedź: ŷ = 14,41

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̄ = 21/6 = 3,5
ȳ = 118/6 = 19,67

Sumy kwadratów i iloczynów

Sxx = Σ(x − x̄)² = 17,5
Sxy = Σ(x − x̄)(y − ȳ) = −46

Prosta regresji

b = −46/17,5 = −2,629
a = 19,67 − (−2,629)·3,5 = 28,87
ŷ = 28,87 − 2,629·x

Prognoza

ŷ = 28,87 + (−2,629)·5,5 = 14,41

Odpowiedź

ŷ = 14,41
9.Testy chi-kwadrat

Odpowiedź: χ² = 6,86

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

palący: 40·70/140 = 20, 40·70/140 = 20
niepalący: 100·70/140 = 50, 100·70/140 = 50

Statystyka chi-kwadrat

χ² = (27 − 20)²/20 + (13 − 20)²/20 + (43 − 50)²/50 + (57 − 50)²/50
χ² = 6,86

Odpowiedź

χ² = 6,86

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ź: χ² = 21,07

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

wiek 18–30: 42, 18
wiek 31–50: 56, 24
wiek 51+: 28, 12

Statystyka chi-kwadrat

χ² = (33 − 42)²/42 + (27 − 18)²/18 + (54 − 56)²/56 + (26 − 24)²/24 + (39 − 28)²/28 + (1 − 12)²/12
χ² = 21,07

Decyzja

df = 2, χ²crit = 5,991
21,07 > 5,991

Odpowiedź

χ² = 21,07

χ² jest większe od wartości krytycznej, więc H₀ zostaje odrzucona: zmienne wydają się powiązane.

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