Konular / Bölüm 6

Örnekten popülasyona

Örnekten hesapladığın ortalama, oran ya da varyansın popülasyondaki gerçek değeri hangi aralıkta yakaladığını, ne kadar güvenle söyleyebileceğini bilmek için.

Örnekten çıkan sayı, gerçeği ne kadar yakalar?

16 merinos kuzusunu tarttın, ortalama doğum ağırlığı 3,2 kg çıktı. Sürüdeki bütün kuzuların gerçek ortalaması da 3,2 kg mı? Neredeyse kesinlikle hayır. Başka 16 kuzu seçseydin 3,0 ya da 3,5 çıkabilirdi. Örnekten hesaplanan her sayı, popülasyondaki gerçek değerin bir tahminidir ve biraz yanılır.

Tahmin bu yüzden iki şekilde verilir:

  • Nokta tahmini: Tek bir sayı. “Ortalama doğum ağırlığı 3,2 kg.” Örnekten hesaplanan ortalama, varyans ya da oran popülasyon parametresinin nokta tahminidir. Popülasyon değerine tam eşit çıkması beklenmez.
  • Aralık tahmini (güven sınırları): Nokta tahmininin etrafına, gerçek değeri belli bir olasılıkla içine alan bir aralık çizilir. “Ortalama doğum ağırlığı %95 güvenle 2,0 ile 4,4 kg arasındadır.”

Aralık tahmini, nokta tahmininin söylemediği bir şeyi söyler: ne kadar emin olabileceğini. Dar bir aralık, tahminin hassas olduğunu; geniş bir aralık ise örneğin popülasyon hakkında az şey söylediğini gösterir.

İyi bir nokta tahmini nasıl olmalı?

ÖzellikAnlamıÖrnek
SapmasızlıkTahminin beklenen değeri popülasyon parametresine eşittir; tekrar tekrar örnek alınsa ortalamada doğru değeri bulur.Varyansta nn yerine n−1n - 1‘e bölmemizin sebebi budur.
KararlılıkÖrnek büyüdükçe tahmin popülasyon değerine yaklaşır.1000 kuzunun ortalaması 16 kuzununkinden daha güvenilirdir.
EtkinlikTahminin varyansı küçüktür, yani örnekten örneğe az oynar.Normal dağılan veride aritmetik ortalama, medyandan daha etkindir.
YeterlilikTahmin, örnekteki bilginin hepsini kullanır.Ortalama bütün değerleri kullanır; mod kullanmaz.

Güven aralığının mantığı

Aralığın genişliği üç şeye bağlıdır:

nokta tahmini±tablo deg˘eri⏟gu¨ven du¨zeyi×standart hata⏟yayılım ve n\text{nokta tahmini} \pm \underbrace{\text{tablo değeri}}_{\text{güven düzeyi}} \times \underbrace{\text{standart hata}}_{\text{yayılım ve } n}
  • Güven düzeyi: %99 güven istersen aralık %95’e göre genişler. Daha emin olmanın bedeli, daha az kesin bir cevaptır.
  • Verinin yayılımı: Değerler ne kadar dağınıksa aralık o kadar geniştir.
  • Örnek büyüklüğü: Standart hata n\sqrt{n}‘e bölündüğü için örnek büyüdükçe aralık daralır. Aralığı yarıya indirmek için dört kat örnek gerekir.

Güven sınırları P(−zα/2<z<zα/2)=1−αP(-z_{\alpha/2} < z < z_{\alpha/2}) = 1 - \alpha şeklinde yazılır. %95 güven için α=0,05\alpha = 0{,}05‘tir; aralık iki yanlı olduğu için bu yanılma payı iki kuyruğa paylaştırılır ve tablodan α/2=0,025\alpha/2 = 0{,}025 için değer okunur.

Hangi dağılımı kullanacağım?

Neyin aralığı?DurumDağılım
Bir ortalamaPopülasyon varyansı (σ2\sigma^2) biliniyorZ
Bir ortalamaσ2\sigma^2 bilinmiyor, n>30n > 30Z (σ yerine ss konur)
Bir ortalamaσ2\sigma^2 bilinmiyor, n≤30n \le 30t (n−1n - 1 serbestlik derecesi)
İki ortalama farkıKüçük örneklert (n1+n2−2n_1 + n_2 - 2 serbestlik derecesi)
Bir oran, iki oran farkıSayım verisiZ
VaryansKi-kare (n−1n - 1 serbestlik derecesi)

Neden 30? Küçük örneklerde ss, σ\sigma‘yı oldukça yanlış tahmin edebilir. t dağılımı bu ek belirsizliği hesaba katar: kuyrukları Z’den kalındır, tablo değeri büyüktür ve aralık genişler. Serbestlik derecesi arttıkça t, Z’ye yaklaşır. 30 serbestlik derecesinde %95 için t = 2,042, Z = 1,96; fark %4 kadardır ve bundan sonra iyice küçülür. “n > 30 ise Z” kuralı bu yüzden pratik bir kısayoldur. Bilgisayarla hesap yaparken t’yi her zaman kullanmak biraz daha doğrudur.

Kendi verinle hesapla

Aşağıdaki hesaplayıcıda ders notundaki altı örneğin hepsi yüklü. Birini seç, sonra güven katsayısını ya da gözlem sayısını değiştirip aralığın nasıl daralıp genişlediğine bak.

Neyin aralığı?
Ders notundaki örnekler:

Merinos kuzuların doğum ağırlığının standart sapması 2,4 kg olarak biliniyor. Sürüden şansa bağlı seçilen 16 kuzunun doğum ağırlığı ortalaması 3,2 kg. %95 güvenle popülasyon ortalamasının sınırları nedir?

Kullanılacak dağılım: Z

2,0244,376nokta tahmini 3,2

Popülasyon ortalamasının %95 güvenle 2,024 ile 4,376 arasında olduğu söylenebilir.

Güven sınırları

2,024 <μ<< \mu < 4,376

  1. Popülasyonun standart sapması (σ) bilindiği için Z dağılımı kullanılır.

  2. Güven katsayısı 1 − α = 0,95, yani α = 0,05. Aralık iki yanlı olduğu için α ikiye bölünür: α/2 = 0,025.

    P(−zα/2<z<zα/2)=1−α=0,95P(-z_{\alpha/2} < z < z_{\alpha/2}) = 1 - \alpha = 0{,}95
  3. Z tablosunda α/2 = 0,025 değerine karşılık gelen z değeri:

    z0,025=1,96z_{0{,}025} = 1{,}96
  4. Ortalamanın standart hatası:

    σxˉ=σn=2,416=0,6\sigma_{\bar{x}} = \dfrac{\sigma}{\sqrt{n}} = \dfrac{2{,}4}{\sqrt{16}} = 0{,}6
  5. μu¨s=xˉ+zα/2 σxˉ=3,2+1,96⋅0,6=4,376\mu_{\text{üs}} = \bar{x} + z_{\alpha/2}\,\sigma_{\bar{x}} = 3{,}2 + 1{,}96 \cdot 0{,}6 = 4{,}376
  6. μalt=xˉ−zα/2 σxˉ=3,2−1,96⋅0,6=2,024\mu_{\text{alt}} = \bar{x} - z_{\alpha/2}\,\sigma_{\bar{x}} = 3{,}2 - 1{,}96 \cdot 0{,}6 = 2{,}024

Formüller

Z ile bir ortalama (σ\sigma biliniyor ya da n>30n > 30):

μu¨s=xˉ+zα/2 σxˉ,μalt=xˉ−zα/2 σxˉ,σxˉ=σn\mu_{\text{üs}} = \bar{x} + z_{\alpha/2}\,\sigma_{\bar{x}}, \qquad \mu_{\text{alt}} = \bar{x} - z_{\alpha/2}\,\sigma_{\bar{x}}, \qquad \sigma_{\bar{x}} = \frac{\sigma}{\sqrt{n}}

t ile bir ortalama (σ\sigma bilinmiyor, n≤30n \le 30):

μu¨s, alt=xˉ±tα/2, n−1 Sxˉ,Sxˉ=S2n\mu_{\text{üs, alt}} = \bar{x} \pm t_{\alpha/2,\,n-1}\,S_{\bar{x}}, \qquad S_{\bar{x}} = \sqrt{\frac{S^2}{n}}

Bir oran:

pu¨s, alt=p^±zα/2p^q^np_{\text{üs, alt}} = \hat{p} \pm z_{\alpha/2}\sqrt{\frac{\hat{p}\hat{q}}{n}}

İki oran farkı:

(p1−p2)u¨s, alt=(p^1−p^2)±zα/2p^1q^1n1+p^2q^2n2(p_1 - p_2)_{\text{üs, alt}} = (\hat{p}_1 - \hat{p}_2) \pm z_{\alpha/2}\sqrt{\frac{\hat{p}_1\hat{q}_1}{n_1} + \frac{\hat{p}_2\hat{q}_2}{n_2}}

İki ortalama farkı (ortak varyansla):

(μ1−μ2)u¨s, alt=(xˉ1−xˉ2)±Sxˉ1−xˉ2 tα/2, n1+n2−2,Sxˉ1−xˉ2=Sp2(1n1+1n2)(\mu_1 - \mu_2)_{\text{üs, alt}} = (\bar{x}_1 - \bar{x}_2) \pm S_{\bar{x}_1 - \bar{x}_2}\,t_{\alpha/2,\,n_1+n_2-2}, \qquad S_{\bar{x}_1 - \bar{x}_2} = \sqrt{S_p^2\left(\frac{1}{n_1} + \frac{1}{n_2}\right)} Sp2=(n1−1)S12+(n2−1)S22n1+n2−2S_p^2 = \frac{(n_1 - 1)S_1^2 + (n_2 - 1)S_2^2}{n_1 + n_2 - 2}

Ki-kare ile varyans:

σu¨s2=S2(n−1)χA2,σalt2=S2(n−1)χB2\sigma^2_{\text{üs}} = \frac{S^2(n-1)}{\chi^2_A}, \qquad \sigma^2_{\text{alt}} = \frac{S^2(n-1)}{\chi^2_B}

Burada χA2\chi^2_A ve χB2\chi^2_B ki-kare tablosundan okunur. %95 güven ve 29 serbestlik derecesi için χA2=χ0,975; 292=16,05\chi^2_A = \chi^2_{0{,}975;\,29} = 16{,}05 (küçük değer, üst sınırı verir), χB2=χ0,025; 292=45,72\chi^2_B = \chi^2_{0{,}025;\,29} = 45{,}72 (büyük değer, alt sınırı verir). Ki-kare simetrik olmadığı için aralık da örnek varyansına göre simetrik değildir.

Sonucu doğru okumak

Aralık, tek tek değerlerin değil, popülasyon ortalamasının aralığıdır. Kuzu örneğinde 2,024–4,376 kg aralığı “kuzuların doğum ağırlığı bu aralıkta değişir” demek değildir; tek tek kuzuların ağırlığı çok daha geniş bir aralıkta dağılır. Doğru okuma: “Sürünün ortalama doğum ağırlığının %95 güvenle 2,024 ile 4,376 kg arasında olduğunu söyleyebiliriz.”

%95 ne demek? Aynı yöntemle yüz farklı örnekten yüz ayrı aralık hesaplasan, bunların yaklaşık 95’i gerçek popülasyon ortalamasını içine alır. Senin elindeki aralık o 95’ten biri olabilir de olmayabilir de; yöntem uzun vadede %95 oranında doğru sonuç verir.

Farkın aralığı sıfırı içeriyorsa: İki ortalama ya da iki oran arasındaki farkın aralığı sıfırı içine alıyorsa (alt sınır negatif, üst sınır pozitif), iki grup arasında gerçekte hiç fark olmayabilir demektir. Bu, aynı güven düzeyinde yapılan iki yanlı testte farkın önemsiz çıkmasıyla aynı anlama gelir. Negatif alt sınır sıfıra yuvarlanmaz; “birinci grup ikinciden biraz düşük de olabilir” demektir. Ayçiçeği örneğinde aralık yaklaşık −8,1 ile 20,1 cm arasıdır: A çeşidi B’den 20 cm daha büyük tabla da verebilir, 8 cm daha küçük de. Bu iki çeşit arasında bu veriyle bir fark gösterilemez.

Oranlarda: Oran negatif olamayacağı için hesaplanan alt sınır sıfırın altına düşerse 0 olarak okunur. Bu genellikle olay sayısının çok az olduğunu gösterir; o zaman normal yaklaşım zaten güvenilir değildir.