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
Masy worków cementu są niezależne i mają rozkład normalny ze średnią 25 kg i odchyleniem standardowym 0,5 kg. Wyznacz prawdopodobieństwo P, że łączna masa 10 z nich przekroczy 248,4 kg. Zaokrąglij z do dwóch miejsc po przecinku.
μ = 25 kgσ = 0,5 kgn = 10> 248,4 kg
P =
2.Rozkład normalnyTrudny
Rozkład wyników testu IQ jest normalny ze średnią 100 pkt i odchyleniem standardowym 15 pkt. Przeciągnij wartość graniczną tak, aby zacieniowane dolne 95% leżało na lewo od niej. Znajdź wartość graniczną x.
μ = 100 pktσ = 15 pkt95%
x = pkt
3.Przedziały ufnościŚredni
Próba 25 elementów ma średnią 90,5 i odchylenie standardowe s = 13. Populacja ma w przybliżeniu rozkład normalny. Wyznacz dolną granicę przedziału ufności 99% dla średniej.
n = 25x̄ = 90,5s = 1399%t(24) = 2,797
dolna granica =
4.Przedziały ufnościTrudny
Średnice 11 toczonych wałków mają odchylenie standardowe z próby 7,3 µm. Zakładając, że średnice mają rozkład normalny, wyznacz σ: dolna granica przedziału ufności 95 % dla odchylenia standardowego populacji. Przy df = 10 2,5 % rozkładu χ² leży poniżej 3,247, a 2,5 % powyżej 20,483.
n = 11s = 7,3 µmχ² = 3,247, 20,483
σ = µm
5.Testy hipotez i przedziały ufnościŚredni
Twierdzi się, że średnia populacji wynosi 500; odchylenie standardowe populacji jest znane i wynosi 10. Losowa próba 100 obserwacji ma średnią 500,82. Przetestuj H₀: μ = 500 przeciw H₁: μ > 500. Wyznacz wartość p p.
μ₀ = 500σ = 10n = 100x̄ = 500,82H₁: μ > μ₀
p =
6.Testy hipotez i przedziały ufnościTrudny
W badaniu klinicznym wyzdrowiało 75 z 300 pacjentów leczonych metodą A i 27 z 150 pacjentów leczonych metodą B. Korzystając z frakcji łącznej, wyznacz statystykę testową z dla H₀: p₁ = p₂.
x₁ = 75n₁ = 300x₂ = 27n₂ = 150
z =
7.Regresja liniowa i korelacjaŚredni
Tabela podaje 5 par obserwacji x = temperatura w °C i y = liczba sprzedanych lodów. Wyznacz współczynnik determinacji r² prostej najmniejszych kwadratów.
n = 5
r² =
8.Regresja liniowa i korelacjaTrudny
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 = 2,5.
n = 6x = 2,5
ŷ =
9.Testy chi-kwadratŚredni
Tabela kontyngencji pokazuje liczebności obserwowane. W teście niezależności chi-kwadrat wyznacz liczebność oczekiwaną E dla komórki „pracownicy, samochód”.
suma wiersza = 74suma kolumny = 50N = 129
E =
10.Testy chi-kwadratTrudny
Test zgodności chi-kwadrat sprawdza, czy liczebności w 4 kategoriach mają rozkład Poissona, którego średnią estymuje się z tych samych danych. 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,8438
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 całkowitej
μT = 250 kg, σT = 0,5·√10 = 1,581 kg
Standaryzuj
z = (248,4 − 250)/1,581 = −1,01
Prawdopodobieństwo
P = 1 − Φ(−1,01) = 0,8438
Odpowiedź
P = 0,8438
2.Rozkład normalny
Odpowiedź: x = 124,6 pkt
Od pola z powrotem do z: Φ(z) to pole na lewo od z. Dla pola lewego ogona A znajdź w tablicy z, dla którego Φ(z) = A. Pole prawego ogona A oznacza Φ(z) = 1 − A. Pole środkowe A między −z a z zostawia (1 − A)/2 w każdym ogonie, więc Φ(z) = (1 + A)/2. Dla wartości x użyj x = μ + zσ.
Pole po lewej
Φ(z) = 0,95
z z tablicy
z ≈ 1,64
Z powrotem do pierwotnych jednostek
x = 100 + 1,64·15 = 124,6 pkt
Odpowiedź
x = 124,6 pkt
3.Przedziały ufności
Odpowiedź: dolna granica = 83,23
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 = 24, t = 2,797
Margines błędu
E = 2,797·13/√25 = 7,272
Przedział
90,5 ± 7,272 → [83,23, 97,77]
Odpowiedź
dolna granica = 83,23
4.Przedziały ufności
Odpowiedź: σ = 5,101 µm
Przedział dla wariancji: (n − 1)s²/σ² ma rozkład χ² z n − 1 stopniami swobody. Przedział 95 % dla σ² rozciąga się od (n − 1)s² podzielonego przez górną wartość χ² (2,5 % powyżej niej) do (n − 1)s² podzielonego przez dolną wartość χ² (2,5 % poniżej niej); przedział dla σ to pierwiastki z obu granic. Nie jest on symetryczny względem s.
Suma kwadratów
(n − 1)s² = 10·7,3² = 532,9
Granice dla σ²
532,9/20,483 = 26,02
532,9/3,247 = 164,1
Granica dla σ
σ = √26,02 = 5,101 µm
Odpowiedź
σ = 5,101 µm
Przedział 95 % dla σ: [5,101, 12,81] µm.
5.Testy hipotez i przedziały ufności
Odpowiedź: p = 0,2061
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 = 10/√100 = 1
z = (500,82 − 500)/1 = 0,82
Pole ogona
Φ(0,82) = 0,7939
1 − Φ(0,82) = 0,2061
Wartość p
p = 0,2061
Odpowiedź
p = 0,2061
p > 0,05, więc nie ma podstaw do odrzucenia H₀ na poziomie 5 %.
6.Testy hipotez i przedziały ufności
Odpowiedź: z = 1,672
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̂₁ = 75/300 = 0,25
p̂₂ = 27/150 = 0,18
p̂ = 102/450 = 0,2267
Błąd standardowy przy H₀
SE = √(0,2267·0,7733·(1/300 + 1/150)) = 0,04187
Statystyka testowa
z = (0,25 − 0,18)/0,04187 = 1,672
Odpowiedź
z = 1,672
|z| < 1.96, więc nie ma podstaw do odrzucenia H₀ na poziomie 5 % (test dwustronny).
7.Regresja liniowa i korelacja
Odpowiedź: r² = 0,9611
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̄ = 105/5 = 21
ȳ = 514/5 = 102,8
Sumy kwadratów i iloczynów
Sxx = Σ(x − x̄)² = 40
Sxy = Σ(x − x̄)(y − ȳ) = 180
Syy = Σ(y − ȳ)² = 842,8
Współczynnik determinacji
r² = 180²/(40·842,8) = 0,9611
Odpowiedź
r² = 0,9611
Około 96,11% zmienności y jest wyjaśnione przez prostą.
8.Regresja liniowa i korelacja
Odpowiedź: ŷ = 18,86
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
ȳ = 96/6 = 16
Sumy kwadratów i iloczynów
Sxx = Σ(x − x̄)² = 17,5
Sxy = Σ(x − x̄)(y − ȳ) = −50
Prosta regresji
b = −50/17,5 = −2,857
a = 16 − (−2,857)·3,5 = 26
ŷ = 26 − 2,857·x
Prognoza
ŷ = 26 + (−2,857)·2,5 = 18,86
Odpowiedź
ŷ = 18,86
9.Testy chi-kwadrat
Odpowiedź: E = 28,68
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ść oczekiwana
E = 74·50/129 = 28,68
Odpowiedź
E = 28,68
10.Testy chi-kwadrat
Odpowiedź: χ²crit = 9,21
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 = 4 − 1 − 1 = 2
Tabela
χ²crit(df = 2, α = 0,01) = 9,210
Odpowiedź
χ²crit = 9,210
Odrzuć H₀, jeśli obliczone χ² jest większe od tej wartości.