Məzmuna keç
Educora
Universitet32 dəq81 / 82

Statistik nəticə çıxarma: qiymətləndirmə və fərziyyələrin yoxlanılması

Seçmədən baş məcmuya: qiymətləndirmələr, orta qiymət üçün etibarlılıq intervalları, p-qiyməti, əhəmiyyətlilik səviyyəsi və I/II növ səhvlərlə fərziyyənin yoxlanılması, t-kriteriyası, korrelyasiya və ən kiçik kvadratlar üsulu ilə xətti reqressiya.

Özünü yoxla
Bu dərsdə öyrənəcəksən
  • Seçmə üzrə orta qiymətin və dispersiyanın sürüşməmiş qiymətləndirmələrini hesablamaq
  • Orta qiymət üçün etibarlılıq intervalı qurmaq və onu düzgün şərh etmək
  • z- və ya t-kriteriyasını tətbiq etmək və p-qiymətini düzgün şərh etmək
  • Korrelyasiya əmsalını və ən kiçik kvadratlar üsulu ilə reqressiya xəttini tapmaq

Seçkidən əvvəl sorğu şirkəti 1000 nəfərlə danışır və «52% ± 3%» nəticəsini elan edir. Çay fabriki istehsal etdiyi milyonlarla qutunun hamısını deyil, 30-unu çəkir. Müəllim bilmək istəyir ki, yeni metod balları həqiqətən artırıb, yoxsa sinfin sadəcə bəxti gətirib. Statistik nəticə çıxarma belə suallara cavab verir: təsadüfi seçmə əsasında bütün baş məcmu haqqında etibarlı nəticəyə necə gəlmək olar və bu nəticəyə nə dərəcədə əmin ola bilərik.

Baş məcmu, seçmə və qiymətləndirmələr

Tərif
Baş məcmu və seçmə

Baş məcmu bizi maraqlandıran obyektlərin bütün çoxluğudur (bütün seçicilər, bütün qutular). Seçmə isə həqiqətən ölçdüyümüz hissədir; o, təsadüfi olmalıdır ki, hər obyektin seçilmə şansı eyni olsun. Baş məcmunu təsvir edən ədədlər (μ, σ, p) parametrlər, seçmədən hesablanan və onları qiymətləndirən ədədlər (x̄, s, p̂) isə statistikalar adlanır.

x̄ = (x₁ + x₂ + … + xₙ) / n, s² = ∑(xᵢ − x̄)² / (n − 1)x̄ = (x₁ + x₂ + … + xₙ) / n, s² = ∑(xᵢ − x̄)² / (n − 1)
burada:
  • x̄seçmə orta, μ-nü qiymətləndirir
  • s²seçmə dispersiya, σ²-nı qiymətləndirir
  • sseçmə orta kvadratik meyl
  • nseçmənin həcmi

Hər iki qiymətləndirmə sürüşməmişdir: bütün mümkün seçmələr üzrə orta hesabla M(x̄) = μ və M(s²) = σ².

Niyə n − 1? Meyllər eyni verilənlərə uyğunlaşdırılmış x̄-dən ölçülür, ona görə orta hesabla bir az kiçik alınır. Meyllərin cəmi həmişə 0-dır, yəni onlardan yalnız n − 1-i sərbəstdir (sərbəstlik dərəcələri). n əvəzinə n − 1-ə bölmək bu sürüşməni dəqiq aradan qaldırır.

Beş qutu çay

Qutuda 100 q çay olmalıdır. Beş qutunun kütləsi 98, 102, 101, 97 və 102 q-dır. x̄, s² və s-i tap.

Həllini göstər
x̄ = (98 + 102 + 101 + 97 + 102)/5 = 500/5 = 100 q.
Meyllər: −2, 2, 1, −3, 2 (cəmi 0 ✓); kvadratları: 4, 4, 1, 9, 4, cəmi 22.
s² = 22/(5 − 1) = 5,5 q², s = √5,5 ≈ 2,35 q.
5-ə bölsəydik, 4,4 q² alınar və səpələnmə azaldılmış olardı.

Orta qiymət üçün etibarlılıq intervalları

Tək bir x̄ ədədi μ-nü heç vaxt dəqiq tapmır. Etibarlılıq intervalı elə qurulan aralıqdır ki, seçmələrin 95%-ində μ-nü əhatə etsin. Mərkəzi limit teoreminə görə x̄ ≈ N(μ, σ²/n), ona görə 0,95 ehtimalla μ x̄-dən ən çoxu 1,96 standart xəta məsafəsindədir.

x̄ ± z* · σ/√nx̄ ± z* · σ/√n
burada:
  • z*kritik qiymət: 90% üçün 1,645, 95% üçün 1,96, 99% üçün 2,576
  • σ/√nσ/√nortanın standart xətası
  • σbaş məcmunun məlum orta kvadratik meyli

σ məlum olduqda (və ya n böyük olduqda — onda σ əvəzinə s götürmək olar) istifadə olunur.

x̄ ± t* · s/√n, df = n − 1x̄ ± t* · s/√n, df = n − 1
burada:
  • t*n − 1 sərbəstlik dərəcəli Styudent t-paylanmasının kritik qiyməti
  • sseçmə orta kvadratik meyl

σ məlum olmadıqda, xüsusən kiçik seçmələrdə istifadə olunur: t-paylanmasının «quyruqları» normalınkından qalındır, ona görə interval daha genişdir.

Etibarlılıq səviyyəsiz*t* (df = 4)t* (df = 9)t* (df = 29)
90%1,6452,1321,8331,699
95%1,9602,7762,2622,045
99%2,5764,6043,2502,756
İkitərəfli intervallar üçün kritik qiymətlər. n böyüdükcə t* z*-a yaxınlaşır.
Gündəlik dərs vaxtı

36 tələbədən ibarət təsadüfi seçmədə gündə orta hesabla x̄ = 2,4 saat dərs oxunur; əvvəlki sorğulardan σ = 0,9 saat məlumdur. μ orta qiyməti üçün 95%-lik etibarlılıq intervalı qur.

Həllini göstər
Standart xəta: σ/√n = 0,9/√36 = 0,9/6 = 0,15 saat.
Xəta həddi: 1,96 · 0,15 ≈ 0,29 saat.
İnterval: 2,4 ± 0,29, yəni [2,11; 2,69] saat.
Şərh: burada istifadə olunan üsul seçmələrin 95%-ində həqiqi ortanı tutur.
Kiçik seçmə: t-intervalı

Çay qutularının orta kütləsi üçün 95%-lik etibarlılıq intervalı qur (n = 5, x̄ = 100 q, s ≈ 2,35 q).

Həllini göstər
σ məlum deyil və n kiçikdir, ona görə df = 4 ilə t götürürük: t* = 2,776.
Standart xəta: s/√n = 2,345/√5 ≈ 1,049 q.
Xəta həddi: 2,776 · 1,049 ≈ 2,91 q.
İnterval: [97,09; 102,91] q. O, 100 q-ı əhatə edir, yəni verilənlər etiketə zidd deyil, amma n kiçik olduğu üçün interval genişdir.

Fərziyyələrin yoxlanılması

  1. 1
    Fərziyyələri yaz

    Sıfır fərziyyə H₀ «xüsusi heç nə yoxdur» iddiasıdır (μ = 500 q, metodun təsiri yoxdur). Alternativ fərziyyə H₁ şübhələndiyimiz haldır (μ < 500 q, metod kömək edir).

  2. 2
    Əhəmiyyətlilik səviyyəsini seç

    α — doğru H₀-ı səhvən rədd etmək riskidir; adətən α = 0,05 (bəzən 0,01).

  3. 3
    Kriteriyanın statistikasını hesabla

    O, verilənlərin H₀-dan standart xəta vahidlərində nə qədər uzaq olduğunu göstərir, məsələn, z və ya t.

  4. 4
    p-qiymətini tap

    H₀ doğru olduğu halda müşahidə olunan qədər və ya ondan da «ekstremal» nəticənin alınma ehtimalı.

  5. 5
    Qərar ver

    p ≤ α olarsa, H₀ rədd edilir (nəticə statistik əhəmiyyətlidir). p > α olarsa, H₀ rədd edilmir: verilənlər kifayət qədər güclü sübut deyil — amma bu, H₀-ın doğru olduğunu sübut etmir.

z = (x̄ − μ₀) / (σ/√n), t = (x̄ − μ₀) / (s/√n)z = (x̄ − μ₀) / (σ/√n), t = (x̄ − μ₀) / (s/√n)
burada:
  • μ₀H₀-ın iddia etdiyi orta qiymət
  • zσ məlum olduqda kriteriyanın statistikası (normal paylanma)
  • tσ məlum olmadıqda: n − 1 sərbəstlik dərəcəli t-paylanması
Çörəklər yüngüldürmü?

Çörəkxana iddia edir ki, çörəklərin orta kütləsi 500 q-dır, σ = 12 q. Müfəttiş 36 çörək çəkir və x̄ = 495 q alır. α = 0,05 səviyyəsində çörəklərin iddia olunandan yüngül olduğunu yoxla.

Həllini göstər
H₀: μ = 500 q; H₁: μ < 500 q (birtərəfli kriteriya).
z = (495 − 500)/(12/√36) = −5/2 = −2,5.
p-qiyməti = P(Z ≤ −2,5) = 1 − Φ(2,5) = 1 − 0,9938 ≈ 0,006.
p ≈ 0,006 < 0,05, deməli, H₀ rədd edilir: verilənlər çörəklərin orta hesabla 500 q-dan yüngül olduğuna güclü sübutdur.
Qeyd: bir çörək üçün 5 q kiçik fərqdir, lakin 36 çörəkdə standart xəta cəmi 2 q-dır.
QərarH₀ doğrudurH₀ yanlışdır
H₀ rədd edilirI növ səhv (ehtimalı α)Düzgün qərar (kriteriyanın gücü 1 − β)
H₀ rədd edilmirDüzgün qərarII növ səhv (ehtimalı β)
Kiçik α daha az «yalançı həyəcan», amma daha çox buraxılmış effekt deməkdir; seçmənin həcmini artırmaq hər iki səhvi azaldır.

t-kriteriyası. σ məlum olmadıqda onu s ilə əvəz edir və t-ni Styudent paylanması ilə (df = n − 1) müqayisə edirik. İki seçmə üçün t-kriteriyası iki qrupun, məsələn, yeni və köhnə metodla oxuyan siniflərin ortalarını müqayisə edir: ortaların fərqi onun standart xətasına bölünür. t-paylanmasını Dublindəki Guinness pivə zavodunda işləyən Uilyam Qosset «Student» təxəllüsü ilə 1908-ci ildə dərc etdirib.

Korrelyasiya və xətti reqressiya

r = ∑(xᵢ − x̄)(yᵢ − ȳ) / √(∑(xᵢ − x̄)² · ∑(yᵢ − ȳ)²)r = ∑(xᵢ − x̄)(yᵢ − ȳ) / √(∑(xᵢ − x̄)² · ∑(yᵢ − ȳ)²)
burada:
  • rPirsonun korrelyasiya əmsalı, −1 ≤ r ≤ 1
  • r ≈ ±1nöqtələr demək olar ki, bir düz xətt üzərindədir
  • r ≈ 0xətti asılılıq yoxdur

ŷ = b₀ + b₁x reqressiya xətti ən kiçik kvadratlar üsulu ilə seçilir: o, şaquli xətaların kvadratları cəmini, S(b₀, b₁) = ∑(yᵢ − b₀ − b₁xᵢ)² ifadəsini minimuma endirir. ∂S/∂b₀ və ∂S/∂b₁ xüsusi törəmələrini sıfıra bərabər etsək, iki xətti tənlik (normal tənliklər) alınır; onların həlli belədir:

b₁ = ∑(xᵢ − x̄)(yᵢ − ȳ) / ∑(xᵢ − x̄)², b₀ = ȳ − b₁ · x̄b₁ = ∑(xᵢ − x̄)(yᵢ − ȳ) / ∑(xᵢ − x̄)², b₀ = ȳ − b₁ · x̄
burada:
  • b₁bucaq əmsalı: x 1 vahid artanda y-in orta dəyişməsi
  • b₀sərbəst hədd (x = 0 olduqda ŷ)
  • ŷproqnozlaşdırılan qiymət

Xətt həmişə (x̄; ȳ) nöqtəsindən keçir, r² isə y-in dəyişkənliyinin xətt ilə izah olunan payıdır.

Dərs saatları və ballar

Beş tələbə: dərs saatları x = 1, 2, 3, 4, 5, test balları y = 52, 58, 65, 70, 80. Reqressiya xəttini və r-i tap, 6 saat üçün balı proqnozlaşdır.

Həllini göstər
x̄ = 3, ȳ = 325/5 = 65.
x − x̄: −2, −1, 0, 1, 2; y − ȳ: −13, −7, 0, 5, 15.
∑(x − x̄)(y − ȳ) = 26 + 7 + 0 + 5 + 30 = 68; ∑(x − x̄)² = 10; ∑(y − ȳ)² = 169 + 49 + 0 + 25 + 225 = 468.
b₁ = 68/10 = 6,8 bal/saat, b₀ = 65 − 6,8 · 3 = 44,6, deməli, ŷ = 44,6 + 6,8x.
r = 68/√(10 · 468) ≈ 68/68,41 ≈ 0,994, r² ≈ 0,99.
Proqnoz: ŷ(6) = 44,6 + 40,8 = 85,4 bal.
İnteraktiv
Simulyasiya yüklənir…
Ən kiçik kvadratlar xətti ŷ = 44,6 + 6,8x (a = b₀, b = b₁). O, (1; 52), (2; 58), (3; 65), (4; 70), (5; 80) nöqtələrinin yaxınlığından keçir; istənilən başqa a və ya b xətaların kvadratları cəmini 5,6-dan böyük edir.
Python
import numpy as np
from scipy import stats

tea = np.array([98, 102, 101, 97, 102])
se = tea.std(ddof=1) / np.sqrt(len(tea))
low, high = stats.t.interval(0.95, df=4, loc=tea.mean(), scale=se)
print(round(low, 2), round(high, 2))

t, p = stats.ttest_1samp(tea, 103)
print(round(t, 3), round(p, 3))

hours = [1, 2, 3, 4, 5]
score = [52, 58, 65, 70, 80]
b1, b0 = np.polyfit(hours, score, 1)
print(round(b1, 2), round(b0, 2), round(np.corrcoef(hours, score)[0, 1], 3))
▸ Gözlənilən nəticə
97.09 102.91
-2.86 0.046
6.8 44.6 0.994
scipy nümunələri yoxlayır: çay qutuları üçün t-intervalı, H₀: μ = 103 q fərziyyəsinin t-kriteriyası (t = −2,86, p = 0,046: α = 0,05-də rədd edilir, α = 0,01-də yox) və r ilə birlikdə reqressiya xətti.

Əsas fikirlər

  • x̄ və s² = ∑(xᵢ − x̄)²/(n − 1) μ və σ²-nın sürüşməmiş qiymətləndirmələridir; seçmə təsadüfi olmalıdır.
  • 95%-lik etibarlılıq intervalı: x̄ ± 1,96σ/√n, σ məlum deyilsə, x̄ ± t*·s/√n.
  • p ≤ α olarsa, H₀ rədd edilir; p-qiyməti H₀-ın doğru olması ehtimalı deyil.
  • I növ səhv: doğru H₀ rədd edilir (α); II növ səhv: yanlış H₀ rədd edilmir (β).
  • Reqressiya: b₁ = ∑(x − x̄)(y − ȳ)/∑(x − x̄)², b₀ = ȳ − b₁x̄; xətt (x̄; ȳ)-dən keçir.
  • r yalnız xətti əlaqəni ölçür, korrelyasiya isə səbəbiyyət demək deyil.

Özünü yoxla

10 sual. Hər düzgün cavab XP qazandırır.

1 / 10
n = 100, x̄ = 50, σ = 10. μ üçün 95%-lik etibarlılıq intervalı hansıdır?