Parametre Tahmin Yöntemleri

228 soru

Soru 81Soru

Bir finansal risk analisti, belirli bir küçük işletme kredisi türü için temerrüde düşme olasılığını (θ\theta) modellemektedir. Analist, geçmiş verilere dayanarak θ\theta parametresinin önsel (prior) dağılımını α=3\alpha=3 ve β=7\beta=7 parametreli Beta dağılımı olarak belirlemiştir. Bu kredi türünden rastgele seçilen n=20n=20 kredilik bir örneklem incelendiğinde, x=4x=4 kredinin temerrüde düştüğü gözlemlenmiştir. Kredilerin temerrüde düşme durumu birbirinden bağımsızdır.

Buna göre, karesel hata kayıp fonksiyonu (squared error loss) altında θ\theta parametresinin Bayes tahmin edicisinin aldığı değer aşağıdakilerden hangisidir?

Cevabı ve açıklamayı göster

Cevap: 730\frac{7}{30}

Cevap

Karesel hata kayıp fonksiyonu altında Bayes tahmin edicisi 730\frac{7}{30} olarak hesaplanır.
Karesel hata kayıp fonksiyonu kullanıldığında Bayes tahmin edicisi, sonsal (posterior) dağılımın beklenen değeridir. Önsel dağılım Beta(3,7)Beta(3, 7) ve örneklem bilgisi Binomial(20,θ)Binomial(20, \theta) ile x=4x=4 verildiğinde, sonsal dağılım Beta(3+4,7+204)=Beta(7,23)Beta(3+4, 7+20-4) = Beta(7, 23) olur. Bu dağılımın beklenen değeri 77+23=730\frac{7}{7+23} = \frac{7}{30}'dur.

Adım Adım Çözüm

1
Önsel dağılım ve olabilirlik fonksiyonunu belirle.
Önsel dağılım: θBeta(3,7)\theta \sim Beta(3, 7). Olabilirlik fonksiyonu: XBinomial(20,θ)X \sim Binomial(20, \theta) ve gözlem x=4x=4.
Sonsal dağılımı bulmak için Bayes teoreminin temel bileşenlerine ihtiyaç vardır.
2
Sonsal (posterior) dağılımın parametrelerini hesapla.
Beta-Binom eşlenik (conjugate) yapısından dolayı sonsal dağılım Beta(α+x,β+nx)Beta(\alpha + x, \beta + n - x) olur. αyeni=3+4=7\alpha_{yeni} = 3 + 4 = 7 ve βyeni=7+204=23\beta_{yeni} = 7 + 20 - 4 = 23. Sonsal dağılım: θxBeta(7,23)\theta|x \sim Beta(7, 23).
Bayes tahmin edicisi sonsal dağılım üzerinden hesaplanır.
3
Verilen kayıp fonksiyonuna göre Bayes tahmin edicisini belirle.
Karesel hata kayıp fonksiyonu (squared error loss) için Bayes tahmin edicisi, sonsal dağılımın beklenen değeridir (ortalamasıdır).
Risk fonksiyonunu minimize eden tahmin edici, karesel hata için sonsal ortalamaya eşittir.
4
Sonsal dağılımın beklenen değerini hesapla.
E(θx)=αyeniαyeni+βyeni=77+23=730E(\theta|x) = \frac{\alpha_{yeni}}{\alpha_{yeni} + \beta_{yeni}} = \frac{7}{7 + 23} = \frac{7}{30}.
Beta dağılımının beklenen değer formülü uygulanarak nihai sonuca ulaşılır.

Anahtar Kavram

Karesel Hata Kayıp Fonksiyonu ile Bayes Tahmini ve Beta-Binom Eşlenikliği
Soru 82Soru

Bir büyükşehir belediyesi, hizmet verdiği ilçelerin nüfus yoğunluğu (XX, bin kişi/km 2^2) ile kişi başına düşen aylık ortalama su tüketim miktarı (YY, ton) arasındaki ilişkiyi modellemek istemektedir. Kurulan basit doğrusal regresyon modeli Yi=β0+β1Xi+ϵiY_i = \beta_0 + \beta_1 X_i + \epsilon_i şeklindedir; burada ϵi\epsilon_i hata terimidir. Modelin parametreleri, normal denklemler kullanılarak en küçük kareler yöntemi ile tahmin edilecektir. Toplanan veri setine ait özet istatistikler aşağıda verilmiştir:

i=1n(XiXˉ)2=40 \sum_{i=1}^n (X_i - \bar{X})^2 = 40
i=1n(XiXˉ)(YiYˉ)=20 \sum_{i=1}^n (X_i - \bar{X})(Y_i - \bar{Y}) = -20
Xˉ=5,Yˉ=10 \bar{X} = 5, \quad \bar{Y} = 10

Buna göre, elde edilen en küçük kareler tahmin denklemine dayanarak, nüfus yoğunluğu X=7X = 7 olan bir ilçenin tahmini ortalama su tüketim miktarı (Y^\hat{Y}) aşağıdakilerden hangisidir?

Cevabı ve açıklamayı göster

Cevap: 9.0

Cevap

9.0
En küçük kareler (EKK) yönteminde regresyon doğrusunun eğim parametresi β^1\hat{\beta}_1, çapraz çarpımlar toplamının XX'in sapma kareler toplamına bölünmesiyle hesaplanır ve 2040=0.5\frac{-20}{40} = -0.5 olarak bulunur. Tahmin doğrusu her zaman örneklem ortalamalarından geçeceği için kesişim terimi β^0=Yˉβ^1Xˉ\hat{\beta}_0 = \bar{Y} - \hat{\beta}_1 \bar{X} şeklinde yazılır. Değerler yerine konulduğunda β^0=10(0.5)(5)=12.5\hat{\beta}_0 = 10 - (-0.5)(5) = 12.5 olur. Elde edilen tahmin denklemi Y^=12.50.5X\hat{Y} = 12.5 - 0.5X olarak biçimlenir. Verilen X=7X = 7 değeri bu denklemde yerine yazıldığında, Y^=12.53.5=9.0\hat{Y} = 12.5 - 3.5 = 9.0 sonucuna ulaşılır.

Adım Adım Çözüm

1
En küçük kareler yöntemi ile regresyon doğrusunun eğimini (β^1\hat{\beta}_1) hesapla.
β^1=(XiXˉ)(YiYˉ)(XiXˉ)2=2040=0.5\hat{\beta}_1 = \frac{\sum (X_i - \bar{X})(Y_i - \bar{Y})}{\sum (X_i - \bar{X})^2} = \frac{-20}{40} = -0.5
Normal denklemlerden elde edilen eğim formülü, ortak değişimin (kovaryansın) bağımsız değişkenin varyansına (sapma kareler toplamına) oranıdır.
2
Regresyon doğrusunun kesişim (sabit) terimini (β^0\hat{\beta}_0) hesapla.
β^0=Yˉβ^1Xˉ=10(0.5)(5)=10+2.5=12.5\hat{\beta}_0 = \bar{Y} - \hat{\beta}_1 \bar{X} = 10 - (-0.5)(5) = 10 + 2.5 = 12.5
En küçük kareler regresyon doğrusu her zaman (Xˉ,Yˉ)(\bar{X}, \bar{Y}) örneklem ortalamaları noktasından geçmek zorundadır.
3
Bulunan parametreleri regresyon denklemine yerleştirerek X=7X = 7 için tahmini değeri (Y^\hat{Y}) bul.
Y^=12.50.5(7)=12.53.5=9.0\hat{Y} = 12.5 - 0.5(7) = 12.5 - 3.5 = 9.0
Popülasyon için oluşturulan modelin doğrusal tahmin denklemi Y^=β^0+β^1X\hat{Y} = \hat{\beta}_0 + \hat{\beta}_1 X formundadır.

Anahtar Kavram

En Küçük Kareler Yöntemi ile Tahmin
Soru 83Soru

Bir araştırmacı, il bazında kişi başına düşen sağlık harcamalarını (YY), ilgili ilin ortalama yaş seviyesi (XX) ile açıklamak amacıyla, sabit terim içermeyen (orijinden geçen) Yi=βXi+ϵiY_i = \beta X_i + \epsilon_i modelini kurmuştur.

Araştırmacı, modeldeki hata terimlerinin varyansının XiX_i ile orantılı olduğunu (Var(ϵi)=σ2XiVar(\epsilon_i) = \sigma^2 X_i) varsayarak, standart en küçük kareler yöntemi yerine ağırlıklı bir minimizasyon problemi çözmeye karar vermiştir. Bu doğrultuda, katsayı tahminini aşağıdaki ağırlıklandırılmış hata kareler toplamı fonksiyonunu minimize ederek yapacaktır:

S(β)=i=1n1Xi(YiβXi)2S(\beta) = \sum_{i=1}^{n} \frac{1}{X_i} (Y_i - \beta X_i)^2

Buna göre, araştırmacının bu minimizasyon problemi sonucunda elde edeceği β\beta parametresinin En Küçük Kareler tahmin edicisi (β^\hat{\beta}) aşağıdakilerden hangisidir?

Cevabı ve açıklamayı göster

Cevap: i=1nYii=1nXi\frac{\sum_{i=1}^{n} Y_i}{\sum_{i=1}^{n} X_i}

Cevap

Minimizasyon problemi çözüldüğünde elde edilen doğru tahmin edici, gözlemlerin toplamlarının birbirine oranını ifade eden i=1nYii=1nXi\frac{\sum_{i=1}^{n} Y_i}{\sum_{i=1}^{n} X_i} değeridir.
Doğru cevap, ağırlıklandırılmış hata kareler toplamı fonksiyonunun β\beta'ya göre birinci türevinin alınıp sıfıra eşitlenmesiyle türetilir. Türev işlemi sırasında zincir kuralı gereği parantez içinin türevi olan Xi-X_i katsayısı ile baştaki ağırlık katsayısı olan 1/Xi1/X_i birbirini sadeleştirir. Geriye kalan (YiβXi)=0\sum (Y_i - \beta X_i) = 0 denklemi çözüldüğünde tahmin edicinin YiXi\frac{\sum Y_i}{\sum X_i} oranına eşit olduğu açıkça görülmektedir.

Adım Adım Çözüm

1
Minimize edilecek olan ağırlıklandırılmış hata kareler toplamı fonksiyonunu belirleyin.
S(β)=i=1n1Xi(YiβXi)2S(\beta) = \sum_{i=1}^{n} \frac{1}{X_i} (Y_i - \beta X_i)^2
En Küçük Kareler yöntemi, hata terimlerinin karesel toplamını en aza indiren parametre değerini bulmayı amaçlar.
2
S(β)S(\beta) fonksiyonunun β\beta'ya göre birinci türevini alın.
dS(β)dβ=i=1n1Xi2(YiβXi)(Xi)\frac{dS(\beta)}{d\beta} = \sum_{i=1}^{n} \frac{1}{X_i} \cdot 2(Y_i - \beta X_i) \cdot (-X_i)
Bir fonksiyonun minimum noktasını bulmak için fonksiyona göre birinci türev alınır.
3
Türev ifadesini sıfıra eşitleyerek gerekli sadeleştirmeleri yapın.
2i=1n[XiXi(YiβXi)]=0    2i=1n(YiβXi)=0-2 \sum_{i=1}^{n} \left[ \frac{X_i}{X_i} (Y_i - \beta X_i) \right] = 0 \implies -2 \sum_{i=1}^{n} (Y_i - \beta X_i) = 0
Türevi sıfıra eşitlemek normal denklemleri elde etmemizi sağlar. Ağırlık teriminin paydasındaki XiX_i ile zincir kuralından gelen Xi-X_i birbirini götürür.
4
Elde edilen normal denklemi β\beta için çözerek tahmin ediciyi (β^\hat{\beta}) yalnız bırakın.
i=1nYiβ^i=1nXi=0    β^i=1nXi=i=1nYi    β^=i=1nYii=1nXi\sum_{i=1}^{n} Y_i - \hat{\beta} \sum_{i=1}^{n} X_i = 0 \implies \hat{\beta} \sum_{i=1}^{n} X_i = \sum_{i=1}^{n} Y_i \implies \hat{\beta} = \frac{\sum_{i=1}^{n} Y_i}{\sum_{i=1}^{n} X_i}
Denklem çözüldüğünde en küçük kareler tahmin edicisi parametre formu elde edilmiş olur.

Anahtar Kavram

Ağırlıklı En Küçük Kareler Yönteminde (WLS) Normal Denklemlerin Türetilmesi

Alternatif Yöntem

Matematiksel beklentiler ve momentler kullanılarak, WLS (Ağırlıklı En Küçük Kareler) tahmininde ağırlıklandırılmış değişkenler oluşturulabilir: Yi=Yi/XiY_i^* = Y_i / \sqrt{X_i} ve Xi=Xi/Xi=XiX_i^* = X_i / \sqrt{X_i} = \sqrt{X_i}. Bu yeni değişkenler ile orijinden geçen regresyon yapıldığında β^=XiYi(Xi)2=Xi(Yi/Xi)(Xi)2=YiXi\hat{\beta} = \frac{\sum X_i^* Y_i^*}{\sum (X_i^*)^2} = \frac{\sum \sqrt{X_i} (Y_i / \sqrt{X_i})}{\sum (\sqrt{X_i})^2} = \frac{\sum Y_i}{\sum X_i} sonucu yine elde edilir.
Tahmini Süre:2m 0s
Soru 84Soru

Bir çevre mühendisi, bir göldeki cıva konsantrasyonunu (mikrogram/litre) ölçmektedir. Cıva konsantrasyonu logaritmik olarak normal dağıldığından, rastgele seçilen nn farklı su örneğinden elde edilen ölçümler (X1,X2,,XnX_1, X_2, \dots, X_n) Log-Normal dağılım ile modellenmiştir. Dağılımın temelini oluşturan varyans değerinin σ2=1\sigma^2 = 1 olduğu bilinmektedir. Olasılık yoğunluk fonksiyonu, bilinmeyen θ\theta parametresi (<θ<-\infty < \theta < \infty) için aşağıdaki gibi verilmiştir:

f(x;θ)=1x2πexp(12(lnxθ)2),x>0f(x; \theta) = \frac{1}{x \sqrt{2\pi}} \exp\left(-\frac{1}{2} (\ln x - \theta)^2\right), \quad x > 0

Buna göre, Neyman Çarpanlara Ayırma Teoremi (Neyman Factorization Theorem) dikkate alındığında, bilinmeyen θ\theta parametresi için yeterli (sufficient) istatistik aşağıdakilerden hangisidir?

Cevabı ve açıklamayı göster

Cevap: i=1nlnXi\sum_{i=1}^n \ln X_i

Cevap

θ\theta parametresi için yeterli istatistik i=1nlnXi\sum_{i=1}^n \ln X_i değeridir.
Örneklem için olabilirlik (likelihood) fonksiyonu yazıldığında L(θ;x)=i=1n1xi2πexp(12(lnxiθ)2)L(\theta; \mathbf{x}) = \prod_{i=1}^n \frac{1}{x_i \sqrt{2\pi}} \exp\left(-\frac{1}{2} (\ln x_i - \theta)^2\right) elde edilir. Üstel ifade içindeki tam kare açılımı yapılıp terimler düzenlendiğinde fonksiyon; parametreyi ve istatistiği birlikte içeren g(T(x),θ)=exp(θi=1nlnxinθ22)g(T(\mathbf{x}), \theta) = \exp\left(\theta \sum_{i=1}^n \ln x_i - \frac{n\theta^2}{2}\right) ile sadece gözlemlere bağlı olan h(x)=(2π)n/2(i=1nxi)1exp(12i=1n(lnxi)2)h(\mathbf{x}) = (2\pi)^{-n/2} \left(\prod_{i=1}^n x_i\right)^{-1} \exp\left(-\frac{1}{2} \sum_{i=1}^n (\ln x_i)^2\right) fonksiyonlarının çarpımı şeklinde yazılır. Neyman Çarpanlara Ayırma Teoremi gereği, gg fonksiyonunun içindeki T(x)=i=1nlnXiT(\mathbf{x}) = \sum_{i=1}^n \ln X_i terimi θ\theta için yeterli (sufficient) istatistiktir.

Adım Adım Çözüm

1
nn çaplı örneklem için ortak olasılık yoğunluk fonksiyonunu (olabilirlik fonksiyonunu) yazın.
L(θ;x)=i=1n[1xi2πexp(12(lnxiθ)2)]L(\theta; \mathbf{x}) = \prod_{i=1}^n \left[ \frac{1}{x_i \sqrt{2\pi}} \exp\left(-\frac{1}{2} (\ln x_i - \theta)^2\right) \right]
Yeterli istatistiği bulmanın ilk adımı, tüm örneklemin ortak dağılımını temsil eden olabilirlik fonksiyonunu oluşturmaktır.
2
Çarpım sembolünü üstel ifadenin içine toplam olarak geçirin ve tam kare açılımını uygulayın.
L(θ;x)=(2π)n/2(i=1nxi)1exp(12i=1n((lnxi)22θlnxi+θ2))L(\theta; \mathbf{x}) = (2\pi)^{-n/2} \left(\prod_{i=1}^n x_i\right)^{-1} \exp\left(-\frac{1}{2} \sum_{i=1}^n \left((\ln x_i)^2 - 2\theta\ln x_i + \theta^2\right)\right)
Neyman teoremini uygulayabilmek için üstel ifadedeki parametre ile gözlem değerlerinin ilişkisini açıkça görmek gereklidir.
3
Fonksiyonu g(T(x),θ)g(T(\mathbf{x}), \theta) ve h(x)h(\mathbf{x}) formunda iki çarpana ayırın.
L(θ;x)=exp(θi=1nlnxinθ22)g(T(x),θ)(2π)n/2(i=1nxi)1exp(12i=1n(lnxi)2)h(x)L(\theta; \mathbf{x}) = \underbrace{\exp\left(\theta \sum_{i=1}^n \ln x_i - \frac{n\theta^2}{2}\right)}_{g(T(\mathbf{x}), \theta)} \cdot \underbrace{(2\pi)^{-n/2} \left(\prod_{i=1}^n x_i\right)^{-1} \exp\left(-\frac{1}{2} \sum_{i=1}^n (\ln x_i)^2\right)}_{h(\mathbf{x})}
Teoreme göre, hem θ\theta'yı hem de istatistiği içeren gg fonksiyonunun içindeki T(x)T(\mathbf{x}) terimi yeterli istatistiktir.

Anahtar Kavram

Neyman Çarpanlara Ayırma Teoremi (Neyman Factorization Theorem)
Soru 85Soru

Bir sigorta şirketinin aktüerya departmanı, poliçe sahiplerinin hasar taleplerini incelemektedir. Bir poliçe sahibinin r=4r=4 adet hasar talebinde bulunana kadar geçen ay sayısı XX rassal değişkeni ile modellenmektedir. XX'in olasılık fonksiyonu, pp (0<p<10 < p < 1) herhangi bir ayda hasar talebi gelme olasılığı olmak üzere,

P(X=x;p)=(x13)p4(1p)x4,x=4,5,6,P(X=x; p) = \binom{x-1}{3} p^4 (1-p)^{x-4}, \quad x = 4, 5, 6, \dots

şeklinde verilmiştir.

Rastgele seçilen nn adet poliçe sahibi için 4. hasar talebine kadar geçen ay sayıları X1,X2,,XnX_1, X_2, \dots, X_n olarak gözlemlenmiştir.

Neyman Çarpanlara Ayırma Teoremi dikkate alındığında, pp parametresi için yeterli (süfisan) istatistik aşağıdakilerden hangisidir?

Cevabı ve açıklamayı göster

Cevap: i=1nXi\sum_{i=1}^n X_i

Cevap

pp parametresi için yeterli istatistik, gözlemlerin toplamı olan i=1nXi\sum_{i=1}^n X_i'dir.
Neyman Çarpanlara Ayırma Teoremi'ne göre olabilirlik fonksiyonu L(p;x)=g(T(x)p)h(x)L(p; \mathbf{x}) = g(T(\mathbf{x})|p) \cdot h(\mathbf{x}) şeklinde iki parçaya ayrılır. Örneklem verileri için olabilirlik fonksiyonu düzenlendiğinde parametre pp'yi içeren kısım p4n(1p)Xi4np^{4n}(1-p)^{\sum X_i - 4n} olarak bulunur. Bu ifade örnekleme yalnızca i=1nXi\sum_{i=1}^n X_i üzerinden bağlı olduğundan, doğru cevap gözlemlerin toplamını veren istatistiktir.

Adım Adım Çözüm

1
Verilen olasılık fonksiyonunu kullanarak nn birimlik örneklem için olabilirlik (likelihood) fonksiyonunu yazmak.
L(p;x)=i=1n[(xi13)p4(1p)xi4]L(p; \mathbf{x}) = \prod_{i=1}^n \left[ \binom{x_i-1}{3} p^4 (1-p)^{x_i-4} \right]
Yeterli istatistiği bulmak için öncelikle örneklemin ortak olasılık dağılımını temsil eden olabilirlik fonksiyonuna ihtiyaç vardır.
2
Olabilirlik fonksiyonundaki çarpım işlemini düzenleyerek ifadeyi sadeleştirmek ve üslü ifadeleri toplamak.
L(p;x)=p4n(1p)i=1n(xi4)i=1n(xi13)=[p4n(1p)i=1nxi4n]×[i=1n(xi13)]L(p; \mathbf{x}) = p^{4n} (1-p)^{\sum_{i=1}^n (x_i-4)} \prod_{i=1}^n \binom{x_i-1}{3} = \left[ p^{4n} (1-p)^{\sum_{i=1}^n x_i - 4n} \right] \times \left[ \prod_{i=1}^n \binom{x_i-1}{3} \right]
Neyman Çarpanlara Ayırma Teoremi'ni uygulayabilmek için parametreyi içeren terimler ile içermeyen terimleri belirgin şekilde ayırmak gerekir.
3
Neyman Çarpanlara Ayırma Teoremi kuralını (L=g(T(x)p)h(x)L = g(T(\mathbf{x}) | p) \cdot h(\mathbf{x})) uygulamak.
Parametreye bağlı kısım: g(T(x)p)=p4n(1p)i=1nxi4ng(T(\mathbf{x}) | p) = p^{4n} (1-p)^{\sum_{i=1}^n x_i - 4n}
Parametreden bağımsız kısım: h(x)=i=1n(xi13)h(\mathbf{x}) = \prod_{i=1}^n \binom{x_i-1}{3}
Teoreme göre, pp'ye bağlı olan gg fonksiyonu, örnekleme yalnızca yeterli istatistik T(x)T(\mathbf{x}) üzerinden bağlı olmalıdır.
4
gg fonksiyonu üzerinden yeterli istatistiği teşhis etmek.
T(x)=i=1nXiT(\mathbf{x}) = \sum_{i=1}^n X_i
Parametreye bağlı olan gg fonksiyonu incelendiğinde, örneklem verilerinin (gözlemlerin) fonksiyon içinde sadece toplam (i=1nxi\sum_{i=1}^n x_i) formunda yer aldığı görülmektedir.

Anahtar Kavram

Neyman Çarpanlara Ayırma Teoremi ile Yeterli İstatistik Bulma
Tahmini Süre:1m 30s
Soru 86Soru

Bir üniversitede uygulanan çevrimiçi bir sınavda, öğrencilerin sınavı tamamlama süreleri (dakika cinsinden) XX rassal değişkeni ile modellenmektedir. Çoğu öğrencinin süreyi sonuna kadar kullandığı gözlemlenmiş olup, XX'in olasılık yoğunluk fonksiyonu θ>0\theta > 0 (sınavın maksimum süresi) olmak üzere aşağıdaki gibi tanımlanmıştır:

f(x;θ)={4x3θ4,0<x<θ0,dig˘er durumlardaf(x; \theta) = \begin{cases} \frac{4x^3}{\theta^4}, & 0 < x < \theta \\ 0, & \text{diğer durumlarda} \end{cases}

Sınavı tamamlayan öğrenciler arasından alınan nn birimlik rastgele bir örneklemin ortalaması Xˉ\bar{X}'tir. Sınav sisteminin analiz modülü, bilinmeyen maksimum θ\theta süresini tahmin etmek için T=XˉT = \bar{X} istatistiğini kullanmaktadır.

Buna göre, TT tahmin edicisinin θ\theta parametresi için Hata Kareler Ortalaması (MSE) aşağıdakilerden hangisidir?

Cevabı ve açıklamayı göster

Cevap: 3n+275nθ2\frac{3n+2}{75n}\theta^2

Cevap

MSE(T)=3n+275nθ2MSE(T) = \frac{3n+2}{75n}\theta^2 ifadesidir.
Doğru Hata Kareler Ortalaması (MSE) değeri hesaplanırken, öncelikle örneklem ortalamasının beklenen değeri ve varyansı bulunur. Verilen olasılık yoğunluk fonksiyonunun integrali üzerinden E[X]=45θE[X] = \frac{4}{5}\theta ve Var(X)=275θ2Var(X) = \frac{2}{75}\theta^2 elde edilir. Buradan T=XˉT=\bar{X}'in yanlı bir tahmin edici olduğu görülür çünkü beklenen değeri parametrenin kendisine değil, 45\frac{4}{5}'ine eşittir. Yan (Bias) miktarı E[Xˉ]θ=15θE[\bar{X}] - \theta = -\frac{1}{5}\theta olur. MSE formülü MSE=Var(Xˉ)+Bias2MSE = Var(\bar{X}) + Bias^2 olduğundan, 275nθ2+(15θ)2=275nθ2+125θ2\frac{2}{75n}\theta^2 + (-\frac{1}{5}\theta)^2 = \frac{2}{75n}\theta^2 + \frac{1}{25}\theta^2 işlemi yapılır. Paydalar 75n75n değerinde eşitlendiğinde sonuç 3n+275nθ2\frac{3n+2}{75n}\theta^2 olarak eksiksiz bulunur.

Adım Adım Çözüm

1
Kitle dağılımının beklenen değerini (E[X]E[X]) ve varyansını (Var(X)Var(X)) hesaplayın.
E[X]=0θx4x3θ4dx=[4x55θ4]0θ=45θE[X] = \int_0^\theta x \frac{4x^3}{\theta^4} dx = \left[ \frac{4x^5}{5\theta^4} \right]_0^\theta = \frac{4}{5}\theta
E[X2]=0θx24x3θ4dx=[4x66θ4]0θ=23θ2E[X^2] = \int_0^\theta x^2 \frac{4x^3}{\theta^4} dx = \left[ \frac{4x^6}{6\theta^4} \right]_0^\theta = \frac{2}{3}\theta^2
Var(X)=E[X2](E[X])2=23θ2(45θ)2=504875θ2=275θ2Var(X) = E[X^2] - (E[X])^2 = \frac{2}{3}\theta^2 - \left(\frac{4}{5}\theta\right)^2 = \frac{50-48}{75}\theta^2 = \frac{2}{75}\theta^2
Örneklem ortalamasının özelliklerini türetebilmek için öncelikle popülasyonun (kitlenin) beklenen değeri ve varyansı bulunmalıdır.
2
T=XˉT = \bar{X} tahmin edicisinin beklenen değerini ve varyansını belirleyin.
E[Xˉ]=E[X]=45θE[\bar{X}] = E[X] = \frac{4}{5}\theta
Var(Xˉ)=Var(X)n=275nθ2Var(\bar{X}) = \frac{Var(X)}{n} = \frac{2}{75n}\theta^2
Örneklem ortalamasının dağılımına ait özellikler bağımsız örneklemler üzerinden ifade edilir.
3
Tahmin edicinin yan (bias) miktarını türetin ve karesini alın.
Bias(Xˉ)=E[Xˉ]θ=45θθ=15θBias(\bar{X}) = E[\bar{X}] - \theta = \frac{4}{5}\theta - \theta = -\frac{1}{5}\theta
Bias(Xˉ)2=(15θ)2=125θ2Bias(\bar{X})^2 = \left(-\frac{1}{5}\theta\right)^2 = \frac{1}{25}\theta^2
Tahmin edici yansız olmadığı için (E[Xˉ]θE[\bar{X}] \neq \theta), Hata Kareler Ortalaması hesabında yan miktarının karesine ihtiyaç vardır.
4
Hata Kareler Ortalaması (MSE) formülünü uygulayın.
MSE(Xˉ)=Var(Xˉ)+Bias(Xˉ)2=275nθ2+125θ2MSE(\bar{X}) = Var(\bar{X}) + Bias(\bar{X})^2 = \frac{2}{75n}\theta^2 + \frac{1}{25}\theta^2
Paydaları 75n75n üzerinde eşitlediğimizde (125θ2=3n75nθ2\frac{1}{25}\theta^2 = \frac{3n}{75n}\theta^2):
MSE(Xˉ)=3n+275nθ2MSE(\bar{X}) = \frac{3n+2}{75n}\theta^2
Bir tahmin edicinin genel performansı MSE üzerinden varyans ve yan karesi toplamı ile ifade edilir.

Anahtar Kavram

Tahmin ediciler için yansızlık durumunun kontrol edilmesi ve yanlı tahmin edicilerde performansın MSE (Hata Kareler Ortalaması = Varyans + Yan^2) kullanılarak ölçülmesi.
Soru 87Soru

Bir madencilik işletmesinde çıkarılan cevherlerdeki belirli bir mineralin yoğunluğu XX rassal değişkeni ile modellenmektedir. XX'in olasılık yoğunluk fonksiyonu aşağıdaki gibi verilmiştir:

f(x;θ)={3x2θ3,0<x<θ0,dig˘er durumlardaf(x; \theta) = \begin{cases} \frac{3x^2}{\theta^3}, & 0 < x < \theta \\ 0, & \text{diğer durumlarda} \end{cases}

Burada θ>0\theta > 0 bilinmeyen bir parametredir. Bu dağılımdan çekilen nn birimlik rastgele bir örneklem X1,X2,,XnX_1, X_2, \dots, X_n ve örneklem ortalaması Xˉ\bar{X} olmak üzere, θ\theta parametresini tahmin etmek için T=2XˉT = 2\bar{X} tahmin edicisi önerilmiştir.

Buna göre, TT tahmin edicisinin hata kareler ortalaması (MSE) aşağıdakilerden hangisidir?

Cevabı ve açıklamayı göster

Cevap: (5n+3)θ220n\frac{(5n+3)\theta^2}{20n}

Cevap

(5n+3)θ220n\frac{(5n+3)\theta^2}{20n} ifadesidir. Hata Kareler Ortalaması (MSE), tahmin edicinin varyansı ile yan miktarının karesinin toplamına eşittir.
Doğru seçenek, hata kareler ortalamasının (MSE) varyans ve yan miktarının karesinin toplamı (MSE=Var(T)+[Bias(T)]2MSE = Var(T) + [Bias(T)]^2) kuralına uygun olarak eksiksiz hesaplandığı ifadedir. Var(T)=3θ220nVar(T) = \frac{3\theta^2}{20n} ve Bias(T)2=θ24Bias(T)^2 = \frac{\theta^2}{4} değerleri payda eşitlenerek toplandığında (5n+3)θ220n\frac{(5n+3)\theta^2}{20n} sonucuna ulaşılır.

Adım Adım Çözüm

1
Kitlenin beklenen değerini (E[X]E[X]) ve varyansını (Var(X)Var(X)) hesapla.
E[X]=0θx3x2θ3dx=3θ3[x44]0θ=3θ4E[X] = \int_{0}^{\theta} x \frac{3x^2}{\theta^3} dx = \frac{3}{\theta^3} [\frac{x^4}{4}]_{0}^{\theta} = \frac{3\theta}{4}. E[X2]=0θx23x2θ3dx=3θ25E[X^2] = \int_{0}^{\theta} x^2 \frac{3x^2}{\theta^3} dx = \frac{3\theta^2}{5}. Var(X)=3θ25(3θ4)2=3θ280Var(X) = \frac{3\theta^2}{5} - (\frac{3\theta}{4})^2 = \frac{3\theta^2}{80}.
Örneklem ortalamasına dayalı bir istatistiğin özelliklerini bulmak için kitle momentlerine ihtiyaç vardır.
2
T=2XˉT = 2\bar{X} tahmin edicisinin beklenen değerini ve yan (bias) miktarını hesapla.
E[T]=2E[Xˉ]=2(3θ4)=3θ2E[T] = 2 E[\bar{X}] = 2(\frac{3\theta}{4}) = \frac{3\theta}{2}. Yan miktarı: Bias(T)=E[T]θ=3θ2θ=θ2Bias(T) = E[T] - \theta = \frac{3\theta}{2} - \theta = \frac{\theta}{2}.
MSE formülünde kullanılacak olan tahmin edicinin yanlılık durumunu belirlemek gereklidir.
3
T=2XˉT = 2\bar{X} tahmin edicisinin varyansını hesapla.
Var(T)=Var(2Xˉ)=4Var(Xˉ)=4(Var(X)n)=4(3θ280n)=3θ220nVar(T) = Var(2\bar{X}) = 4 Var(\bar{X}) = 4 (\frac{Var(X)}{n}) = 4 (\frac{3\theta^2}{80n}) = \frac{3\theta^2}{20n}.
MSE formülünün ikinci bileşeni olan varyans değerini bulmak gerekir.
4
MSE(T)=Var(T)+Bias(T)2MSE(T) = Var(T) + Bias(T)^2 formülünü uygula.
MSE(T)=3θ220n+(θ2)2=3θ220n+θ24=3θ2+5nθ220n=(5n+3)θ220nMSE(T) = \frac{3\theta^2}{20n} + (\frac{\theta}{2})^2 = \frac{3\theta^2}{20n} + \frac{\theta^2}{4} = \frac{3\theta^2 + 5n\theta^2}{20n} = \frac{(5n+3)\theta^2}{20n}.
Yanlı bir tahmin edicinin performansını değerlendirmek için toplam hata (MSE) formülünün kullanılması kuraldır.

Anahtar Kavram

Yanlı Tahmin Edicilerde Hata Kareler Ortalaması (MSE)
Tahmini Süre:2m 30s
Soru 88Soru

Bir veri merkezinde kullanılan yüksek kapasiteli belirli bir Katı Hal Sürücüsü (SSD) modelinin ömrü (yıl cinsinden) XX rastgele değişkeni ile modellenmektedir. XX'in olasılık yoğunluk fonksiyonu aşağıda verilmiştir:

f(x;β)=β2xeβx,x>0f(x; \beta) = \beta^2 x e^{-\beta x}, \quad x > 0

Burada β>0\beta > 0 bilinmeyen bir parametredir. Veri merkezindeki disklerden rastgele seçilen bağımsız ve özdeş dağılımlı (i.i.d.) nn adet diskin ömürleri X1,X2,,XnX_1, X_2, \dots, X_n olarak gözlemlenmiş ve örneklem ortalaması Xˉ=1ni=1nXi\bar{X} = \frac{1}{n} \sum_{i=1}^n X_i olarak hesaplanmıştır.

Buna göre, disk ömürlerinin varyansı olan τ(β)=2β2\tau(\beta) = \frac{2}{\beta^2} parametresi için en çok olabilirlik tahmin edicisi (MLE) aşağıdakilerden hangisidir?

Cevabı ve açıklamayı göster

Cevap: Xˉ22\frac{\bar{X}^2}{2}

Cevap

Xˉ22\frac{\bar{X}^2}{2}
Verilen olasılık yoğunluk fonksiyonu için örneklemin olabilirlik fonksiyonu L(β)=β2neβXiXiL(\beta) = \beta^{2n} e^{-\beta \sum X_i} \prod X_i şeklindedir. Log-olabilirlik fonksiyonunun türevi alınıp sıfıra eşitlendiğinde β\beta parametresinin en çok olabilirlik tahmin edicisi β^=2Xˉ\hat{\beta} = \frac{2}{\bar{X}} olarak bulunur. En çok olabilirlik tahmin edicilerinin değişmezlik (invariance) özelliği gereği, τ(β)=2β2\tau(\beta) = \frac{2}{\beta^2} fonksiyonunun tahmin edicisi, β\beta yerine β^\hat{\beta} yazılarak hesaplanır. Bu işlem sonucunda 2(2/Xˉ)2=Xˉ22\frac{2}{(2/\bar{X})^2} = \frac{\bar{X}^2}{2} değerine ulaşılır.

Adım Adım Çözüm

1
Olabilirlik (likelihood) fonksiyonunu oluşturma
L(β)=i=1nβ2XieβXi=β2n(i=1nXi)eβi=1nXiL(\beta) = \prod_{i=1}^n \beta^2 X_i e^{-\beta X_i} = \beta^{2n} \left(\prod_{i=1}^n X_i\right) e^{-\beta \sum_{i=1}^n X_i}
Bağımsız ve özdeş dağılımlı örneklem için ortak olasılık yoğunluk fonksiyonunu bulmak amacıyla.
2
Log-olabilirlik fonksiyonunu hesaplama
lnL(β)=2nlnβ+i=1nlnXiβi=1nXi\ln L(\beta) = 2n \ln \beta + \sum_{i=1}^n \ln X_i - \beta \sum_{i=1}^n X_i
Çarpım durumundaki ifadeyi toplama dönüştürerek parametreye göre türev alma işlemini kolaylaştırmak amacıyla.
3
Parametreye göre türev alıp sıfıra eşitleme
lnL(β)β=2nβi=1nXi=0    β^=2ni=1nXi=2Xˉ\frac{\partial \ln L(\beta)}{\partial \beta} = \frac{2n}{\beta} - \sum_{i=1}^n X_i = 0 \implies \hat{\beta} = \frac{2n}{\sum_{i=1}^n X_i} = \frac{2}{\bar{X}}
Log-olabilirlik fonksiyonunu maksimum yapan β\beta değerini (nokta tahmin edicisini) bulmak amacıyla.
4
Değişmezlik (invariance) özelliğini uygulama
τ(β)^=2β^2=2(2/Xˉ)2=24/Xˉ2=Xˉ22\widehat{\tau(\beta)} = \frac{2}{\hat{\beta}^2} = \frac{2}{(2/\bar{X})^2} = \frac{2}{4/\bar{X}^2} = \frac{\bar{X}^2}{2}
Bir parametrenin fonksiyonunun en çok olabilirlik tahmin edicisi, parametrenin tahmin edicisinin o fonksiyonda yerine konmasıyla elde edildiğinden.

Anahtar Kavram

En Çok Olabilirlik Tahmin Edicisi ve Değişmezlik Özelliği
Tahmini Süre:2m 30s
Soru 89Soru

Bir bölgedeki kamu projelerinin tamamlanma süreleri (ay cinsinden), α\alpha (şekil) ve β\beta (ölçek) parametreli bir Gamma dağılımı ile modellenmektedir. Bir araştırmacı, nn adet projenin tamamlanma sürelerini rastgele seçerek birinci örneklem momentini m1=1ni=1nXim_1 = \frac{1}{n} \sum_{i=1}^n X_i ve ikinci örneklem momentini m2=1ni=1nXi2m_2 = \frac{1}{n} \sum_{i=1}^n X_i^2 olarak hesaplamıştır.

Gamma dağılımının beklenen değeri E(X)=αβE(X) = \alpha\beta ve varyansı Var(X)=αβ2Var(X) = \alpha\beta^2 olduğuna göre, α\alpha parametresi için momentler yöntemi tahmin edicisi (α^\hat{\alpha}) aşağıdakilerden hangisidir?

Cevabı ve açıklamayı göster

Cevap: m12m2m12\frac{m_1^2}{m_2 - m_1^2}

Cevap

m12m2m12\frac{m_1^2}{m_2 - m_1^2} ifadesi α\alpha parametresi için doğru momentler yöntemi tahmin edicisidir.
Momentler yönteminde, parametre sayısı kadar kitle momenti örneklem momentine eşitlenerek denklem sistemi çözülür. Gamma dağılımı için teorik momentler hesaplandığında; E(X)=αβ=m1E(X) = \alpha\beta = m_1 ve E(X2)=Var(X)+E(X)2=αβ2+α2β2=m2E(X^2) = Var(X) + E(X)^2 = \alpha\beta^2 + \alpha^2\beta^2 = m_2 denklemleri kurulur. Birinci denklemden çekilen β=m1/α\beta = m_1/\alpha ifadesi ikinci denklemde yerine konup işlemler sadeleştirildiğinde doğrudan α=m12m2m12\alpha = \frac{m_1^2}{m_2 - m_1^2} sonucuna ulaşılır.

Adım Adım Çözüm

1
Kitle momentlerini dağılımın parametreleri cinsinden ifade edin.
E(X)=αβE(X) = \alpha\beta ve E(X2)=Var(X)+[E(X)]2=αβ2+(αβ)2=αβ2+α2β2E(X^2) = Var(X) + [E(X)]^2 = \alpha\beta^2 + (\alpha\beta)^2 = \alpha\beta^2 + \alpha^2\beta^2
Momentler yönteminde öncelikle teorik kitle momentleri yazılmalıdır.
2
Birinci ve ikinci kitle momentlerini örneklem momentlerine (m1m_1 ve m2m_2) eşitleyin.
1) αβ=m1\alpha\beta = m_1
2) αβ2+α2β2=m2\alpha\beta^2 + \alpha^2\beta^2 = m_2
Parametreleri tahmin edebilmek için teorik momentler ile hesaplanan örneklem momentleri eşitlenir.
3
Oluşan denklem sisteminde β\beta'yı yalnız bırakıp ikinci denklemde yerine yazın.
Birinci denklemden β=m1α\beta = \frac{m_1}{\alpha} bulunur. İkinci denklemde yerine yazılırsa: α(m1α)2+α2(m1α)2=m2\alpha\left(\frac{m_1}{\alpha}\right)^2 + \alpha^2\left(\frac{m_1}{\alpha}\right)^2 = m_2
Sistemdeki iki bilinmeyenden birini yok ederek α\alpha'yı tek bilinmeyen olarak elde etmek için.
4
Elde edilen son denklemi α\alpha için çözün.
m12α+m12=m2    m12α=m2m12    α=m12m2m12\frac{m_1^2}{\alpha} + m_1^2 = m_2 \implies \frac{m_1^2}{\alpha} = m_2 - m_1^2 \implies \alpha = \frac{m_1^2}{m_2 - m_1^2}
İstenilen α\alpha parametresinin tahmin edicisini açık formda matematiksel olarak yazmak için.

Anahtar Kavram

Momentler Yöntemi ile Tahmin (Method of Moments Estimation)
Soru 90Soru

Ekonometrik ve istatistiksel modellemelerde parametre tahmini için sıklıkla En Küçük Kareler (EKK) yöntemi tercih edilmektedir. Yi=β0+β1Xi+ϵiY_i = \beta_0 + \beta_1 X_i + \epsilon_i biçimindeki basit doğrusal modelde, XiX_i açıklayıcı değişken değerlerinin sabit (stokastik olmayan) olduğu kabul edilmektedir.

Elde edilen β^0\hat{\beta}_0 ve β^1\hat{\beta}_1 EKK tahmin edicilerinin kitlenin gerçek parametreleri için yansız (unbiased) birer tahmin edici olabilmesi adına, gözlemlenemeyen hata terimleri (ϵi\epsilon_i) üzerine konulması gereken temel varsayım aşağıdakilerden hangisidir?

Cevabı ve açıklamayı göster

Cevap: Hata terimlerinin beklenen değerinin sıfır olması, E(ϵi)=0E(\epsilon_i) = 0

Cevap

Hata terimlerinin beklenen değerinin sıfır olması, E(ϵi)=0E(\epsilon_i) = 0
Bir tahmin edicinin yansız olabilmesi için beklenen değerinin tahmin edilen parametreye eşit olması gerekir. EKK tahmin edicileri hesaplanırken, doğrusal kombinasyonlar üzerinden E(ϵi)=0E(\epsilon_i) = 0 olduğu kabul edilirse E(β^1)=β1E(\hat{\beta}_1) = \beta_1 sonucuna ulaşılır. Sabit varyans, normal dağılım veya sıfır kovaryans özellikleri yansızlık için değil; etkinlik, normallik veya tutarlılık gibi diğer istatistiksel özellikler için gereklidir.

Adım Adım Çözüm

1
EKK yönteminde eğim parametresi tahmin edicisinin formülünü yazmak.
β^1=(XiXˉ)(YiYˉ)(XiXˉ)2\hat{\beta}_1 = \frac{\sum (X_i - \bar{X})(Y_i - \bar{Y})}{\sum (X_i - \bar{X})^2}
Tahmin edicinin matematiksel yapısını anlamak, beklenen değerini hesaplamanın ilk adımıdır.
2
Tahmin ediciyi gerçek parametre cinsinden ifade etmek.
β^1=β1+wiϵiburada wi=XiXˉ(XiXˉ)2\hat{\beta}_1 = \beta_1 + \sum w_i \epsilon_i \quad \text{burada } w_i = \frac{X_i - \bar{X}}{\sum (X_i - \bar{X})^2}
Yansızlık incelemesi için tahmin ediciyi sabitler (wiw_i, β1\beta_1) ve rastgele değişkenler (ϵi\epsilon_i) olarak ayırmak gerekir.
3
Eşitliğin her iki tarafının beklenen değerini (E) almak.
E(β^1)=E(β1)+wiE(ϵi)=β1+wiE(ϵi)E(\hat{\beta}_1) = E(\beta_1) + \sum w_i E(\epsilon_i) = \beta_1 + \sum w_i E(\epsilon_i)
Yansızlık koşulu olan E(β^1)=β1E(\hat{\beta}_1) = \beta_1 eşitliğinin sağlanması için gereken şartı bulmak.
4
Yansızlık için gereken koşulu belirlemek.
E(ϵi)=0E(\epsilon_i) = 0 varsayımı altında denklem E(β^1)=β1+0=β1E(\hat{\beta}_1) = \beta_1 + 0 = \beta_1 halini alır.
Tahmin edicinin beklenen değerinin gerçek parametreye eşit olmasını sağlayan temel varsayım budur.

Anahtar Kavram

EKK Tahmin Edicilerinin Yansızlık Özelliği ve Gauss-Markov Varsayımları
Soru 91Soru

Türkiye İstatistik Kurumu (TÜİK) tarafından yürütülen bir hanehalkı bütçe anketinde, belirli bir zorunlu tüketim kalemine ait harcama tutarları, ortalaması θ\theta olan bir olasılık dağılımı ile modellenmektedir. Uzmanlar, ellerindeki büyük örneklemi kullanarak θ\theta parametresinin en çok olabilirlik tahmin edicisini (MLE) θ^\hat{\theta} olarak hesaplamıştır.

Çalışmanın ilerleyen aşamalarında, ekonomik bir refah endeksinin hesaplanabilmesi için θ\theta'nın doğrusal olmayan g(θ)=θ2+ln(θ)g(\theta) = \theta^2 + \ln(\theta) şeklindeki bir fonksiyonunun tahmin edilmesi gerekmiştir. Uzmanlar, En Çok Olabilirlik yönteminin özelliklerine dayanarak bu endeksin tahmin edicisini g(θ)^=θ^2+ln(θ^)\widehat{g(\theta)} = \hat{\theta}^2 + \ln(\hat{\theta}) olarak önermiştir.

Verilen bu duruma ve en çok olabilirlik tahmin edicilerinin genel özelliklerine göre, önerilen g(θ)^\widehat{g(\theta)} tahmin edicisi ile ilgili aşağıdakilerden hangisi kesinlikle doğrudur?

Cevabı ve açıklamayı göster

Cevap: Değişmezlik (invariance) özelliği gereği g(θ)g(\theta)'nın en çok olabilirlik tahmin edicisidir ve tutarlı (consistent) bir tahmin edici olmasına rağmen sonlu örneklemlerde genellikle yanlıdır (biased).

Cevap

Doğru cevap, tahmin edicinin değişmezlik özelliği gereği en çok olabilirlik tahmin edicisi olduğunu, tutarlı kaldığını ancak sonlu örneklemlerde genellikle yanlı olduğunu belirten seçenektir.
En Çok Olabilirlik Tahmin Edicilerinin (MLE) en önemli pratik özelliklerinden biri 'değişmezlik' (invariance) özelliğidir. Bu özelliğe göre, bir parametrenin fonksiyonunun MLE'sini bulmak için parametrenin kendi MLE'sini fonksiyonda yerine koymak yeterlidir. Asimptotik özellikler incelendiğinde; Sürekli Fonksiyon Teoremi uyarınca tutarlılık özelliği korunur (örneklem çapı sonsuza giderken gerçek değere yakınsar). Ancak doğrusal olmayan g(θ)=θ2+ln(θ)g(\theta) = \theta^2 + \ln(\theta) gibi fonksiyonlarda, Jensen eşitsizliği nedeniyle sonlu örneklemlerde yansızlık (unbiasedness) ortadan kalkar ve tahmin edici yanlı hale gelir.

Adım Adım Çözüm

1
En çok olabilirlik tahmin edicilerinin değişmezlik (invariance) özelliğini hatırla.
Eğer θ^\hat{\theta}, θ\theta'nın MLE'si ise, herhangi bir g(θ)g(\theta) fonksiyonunun MLE'si g(θ^)g(\hat{\theta}) olur.
Bu özellik, karmaşık fonksiyonların tahmin edilmesinde olabilirlik fonksiyonunun yeniden çözülmesini engeller.
2
Tahmin edicinin asimptotik özelliklerinden olan 'tutarlılık (consistency)' kavramının dönüşümlerde nasıl etkilendiğini değerlendir.
Sürekli Fonksiyon Teoremi (Continuous Mapping Theorem) gereğince, θ^\hat{\theta} tutarlı ise sürekli bir gg fonksiyonu için g(θ^)g(\hat{\theta}) da tutarlıdır.
En çok olabilirlik tahmin edicileri genel düzenlilik koşulları altında tutarlıdır, dolayısıyla g(θ)^\widehat{g(\theta)} da tutarlı olacaktır.
3
Tahmin edicinin 'yansızlık (unbiasedness)' özelliğinin doğrusal olmayan (non-linear) dönüşümlerde nasıl etkilendiğini incele.
g(θ)=θ2+ln(θ)g(\theta) = \theta^2 + \ln(\theta) doğrusal olmayan bir fonksiyondur. Jensen Eşitsizliği gereği, genelde E[g(θ^)]g(E[θ^])E[g(\hat{\theta})] \neq g(E[\hat{\theta}]) olur.
Bu nedenle, θ^\hat{\theta} yansız olsa bile, doğrusal olmayan bir fonksiyonun MLE'si genellikle küçük ve sonlu örneklemlerde yanlı (biased) olur, ancak asimptotik olarak yansızlığa yaklaşır.

Anahtar Kavram

En Çok Olabilirlik Tahmin Edicilerinin Değişmezliği (Invariance) ve Asimptotik Özellikler
Soru 92Soru

Bir finansal risk analisti, bir kredi portföyündeki müşterilerin temerrüde düşme durumlarını bağımsız Bernoulli denemeleri olarak modellemiş ve başarı olasılığını pp (0<p<10 < p < 1) olarak tanımlamıştır. Analist, bu dağılımın varyansı olan θ=p(1p)\theta = p(1-p) parametresini tahmin etmek amacıyla nn hacimli (n2n \geq 2) bir rastgele örneklem üzerinden aşağıdaki üç farklı tahmin ediciyi kurgulamıştır:

I. T1=Xˉ(1Xˉ)T_1 = \bar{X}(1 - \bar{X})
II. T2=nn1Xˉ(1Xˉ)T_2 = \frac{n}{n-1} \bar{X}(1 - \bar{X})
III. T3=X1(1X2)T_3 = X_1(1 - X_2)

Tahmin edicilerin asimptotik davranışları ve tutarlılık (consistency) özellikleri incelendiğinde aşağıdakilerden hangisi doğrudur?

Cevabı ve açıklamayı göster

Cevap: Yalnızca T1T_1 ve T2T_2 tutarlı tahmin edicilerdir; T3T_3 yansız olmasına rağmen varyansı asimptotik olarak sıfıra yakınsamadığı için tutarlı değildir.

Cevap

Yalnızca T1T_1 ve T2T_2 tutarlı tahmin edicilerdir; T3T_3 yansız olmasına rağmen varyansı asimptotik olarak sıfıra yakınsamadığı için tutarlı değildir.
Bir tahmin edicinin tutarlı olması için nn \to \infty iken parametreye olasılıkta yakınsaması gerekir. Xˉ\bar{X} tahmin edicisi Büyük Sayılar Yasası (WLLN) gereği pp'ye olasılıkta yakınsar (XˉPp\bar{X} \xrightarrow{P} p). Sürekli Dönüşüm Teoremi'ne göre g(x)=x(1x)g(x)=x(1-x) sürekli bir fonksiyon olduğundan Xˉ(1Xˉ)Pp(1p)\bar{X}(1-\bar{X}) \xrightarrow{P} p(1-p) olur, yani T1T_1 tutarlıdır. T2T_2 ise T1T_1'in asimptotik olarak 11'e yakınsayan n/(n1)n/(n-1) sabiti ile çarpılmış halidir ve Slutsky Teoremi gereği tutarlıdır. T3T_3 tahmin edicisi E(X1(1X2))=p(1p)E(X_1(1-X_2)) = p(1-p) olduğu için tam yansızdır ancak tahmini yalnızca iki gözleme dayandırır. Örneklem çapı nn sonsuza gitse bile verinin geri kalanını kullanmadığı için varyansı azalmaz ve sıfıra gitmez. Bu sebeple olasılıkta θ\theta'ya yığılma gerçekleşmez ve T3T_3 tutarlı değildir.

Adım Adım Çözüm

1
T1T_1, T2T_2 ve T3T_3 tahmin edicilerinin beklenen değerlerini inceleyerek yansızlık durumlarını belirlemek.
E(T1)=n1np(1p)E(T_1) = \frac{n-1}{n}p(1-p) olup yanlıdır (ancak limitinde asimptotik yansızdır). E(T2)=p(1p)E(T_2) = p(1-p) olup tam yansızdır. E(T3)=E(X1)E(1X2)=p(1p)E(T_3) = E(X_1)E(1-X_2) = p(1-p) olup tam yansızdır.
Tutarlılık analizi için öncelikle tahmin edicilerin asimptotik olarak parametreye yakınsayıp yakınsamadığı (beklenen değer limiti) kontrol edilmelidir.
2
Büyük Sayılar Yasası (WLLN) ve Sürekli Dönüşüm Teoremi kullanarak T1T_1 ve T2T_2'nin tutarlılığını değerlendirmek.
Büyük Sayılar Yasası gereği XˉPp\bar{X} \xrightarrow{P} p olur. Sürekli Dönüşüm Teoremi uygulanırsa Xˉ(1Xˉ)Pp(1p)\bar{X}(1-\bar{X}) \xrightarrow{P} p(1-p) elde edilir. Bu nedenle T1T_1 tutarlıdır. Slutsky Teoremi ile nn11\frac{n}{n-1} \to 1 olduğundan T2T_2 de aynı değere olasılıkta yakınsar ve tutarlıdır.
Olasılıkta yakınsama özellikleri, tahmin edicilerin tutarlılığını ispatlamak için yeterlidir.
3
T3T_3'ün varyansının nn \to \infty durumundaki limitini incelemek.
T3T_3 sadece ilk iki bağımsız gözleme (X1X_1 ve X2X_2) dayanır. Örneklem çapı nn ne kadar büyürse büyüsün yapısı değişmez. Varyansı sabittir ve sıfıra yakınsamaz. Bu nedenle olasılıkta yakınsama gerçekleşmez ve T3T_3 tutarlı değildir.
Yansız bir tahmin edicinin tutarlı olabilmesi için MSE'sinin (ve dolayısıyla varyansının) asimptotik olarak sıfıra gitmesi şarttır.

Anahtar Kavram

Tahmin Edicilerin Özellikleri: Tutarlılık (Consistency)
Soru 93Soru

Bir yenilenebilir enerji uzmanı, belirli bir bölgedeki rüzgâr türbinlerinin arızasız çalışma sürelerini (ay cinsinden) incelemektedir. Birbirinden bağımsız nn adet türbinin arızasız çalışma süreleri X1,X2,,XnX_1, X_2, \dots, X_n rastgele değişkenleri ile modellenmiş olup, bu değişkenlerin (0,θ)(0, \theta) aralığında sürekli düzgün (uniform) dağılıma sahip olduğu bilinmektedir. Burada θ>0\theta > 0 bilinmeyen bir parametredir.

Uzman, Neyman Çarpanlara Ayırma Teoremi'ni (Neyman Factorization Theorem) kullanarak θ\theta parametresi için bir yeterli (sufficient) istatistik belirlemek istemektedir.

Buna göre, aşağıdakilerden hangisi θ\theta parametresi için bir yeterli istatistiktir?

Cevabı ve açıklamayı göster

Cevap: Örneklemdeki en büyük gözlem değeri olan X(n)X_{(n)}

Cevap

θ\theta parametresi için yeterli istatistik, örneklemdeki en büyük gözlem değeri olan X(n)X_{(n)} istatistiğidir.
Düzgün dağılım gibi parametrenin tanım kümesini sınırlandırdığı durumlarda, yeterli istatistik genellikle uç değerlere (sıra istatistiklerine) bağlıdır. Olabilirlik fonksiyonu L(θx)=1θnI(X(n)θ)I(X(1)0)L(\theta | \mathbf{x}) = \frac{1}{\theta^n} I(X_{(n)} \leq \theta) I(X_{(1)} \geq 0) şeklinde yazıldığında, θ\theta'yı içeren çarpan yalnızca en büyük sıra istatistiği olan X(n)X_{(n)}'i barındırdığından, Neyman Çarpanlara Ayırma Teoremi gereği X(n)X_{(n)} yeterli istatistiktir.

Adım Adım Çözüm

1
Dağılımın olasılık yoğunluk fonksiyonunu (pdf) yazmak.
f(x;θ)=1θI(0xθ)f(x; \theta) = \frac{1}{\theta} I(0 \leq x \leq \theta)
Olabilirlik fonksiyonunu oluşturmak için kitle dağılımının matematiksel formuna ihtiyaç vardır. II gösterge (indicator) fonksiyonudur.
2
Örneklem için ortak olasılık yoğunluk (olabilirlik) fonksiyonunu oluşturmak.
L(θx)=i=1n1θI(0xiθ)=1θni=1nI(0xiθ)L(\theta | \mathbf{x}) = \prod_{i=1}^{n} \frac{1}{\theta} I(0 \leq x_i \leq \theta) = \frac{1}{\theta^n} \prod_{i=1}^{n} I(0 \leq x_i \leq \theta)
Neyman Çarpanlara Ayırma Teoremi'ni uygulayabilmek için nn çaplı bağımsız örneklemin ortak dağılımı yazılmalıdır.
3
Gösterge fonksiyonlarının çarpımını sıra istatistikleri cinsinden ifade etmek.
Tüm xix_i'lerin 00 ile θ\theta arasında olması, en küçük değerin 00'dan büyük veya eşit olması ve en büyük değerin θ\theta'dan küçük veya eşit olması demektir: I(0x(1))I(x(n)θ)I(0 \leq x_{(1)}) \cdot I(x_{(n)} \leq \theta). Böylece L(θx)=1θnI(x(n)θ)I(x(1)0)L(\theta | \mathbf{x}) = \frac{1}{\theta^n} I(x_{(n)} \leq \theta) \cdot I(x_{(1)} \geq 0) elde edilir.
Parametreye (θ\theta) bağlı sınırları, örneklemin uygun istatistikleriyle (minimum ve maksimum) ayrıştırmak gerekir.
4
Neyman Çarpanlara Ayırma Teoremi'ni uygulayarak g(T(x),θ)g(T(x), \theta) ve h(x)h(x) çarpanlarını belirlemek.
g(T(x),θ)=1θnI(x(n)θ)g(T(\mathbf{x}), \theta) = \frac{1}{\theta^n} I(x_{(n)} \leq \theta) ve h(x)=I(x(1)0)h(\mathbf{x}) = I(x_{(1)} \geq 0). Parametreye bağlı çarpan gg fonksiyonu X(n)X_{(n)} istatistiğini içerdiğinden yeterli istatistik T(x)=X(n)T(\mathbf{x}) = X_{(n)} olur.
Teoreme göre, L(θx)=g(T(x),θ)h(x)L(\theta | \mathbf{x}) = g(T(\mathbf{x}), \theta) \cdot h(\mathbf{x}) şeklinde yazılabiliyorsa T(x)T(\mathbf{x}) yeterli istatistiktir.

Anahtar Kavram

Neyman Çarpanlara Ayırma Teoremi ve Sıra İstatistikleri
Soru 94Soru

Bir meteoroloji istasyonunda ölçülen günlük rüzgar hızları, (0,θ)(0, \theta) aralığında bir düzgün dağılıma (Uniform dağılım) sahip rastgele değişkenler olarak modellenmektedir.

Bu istasyondan birbirinden bağımsız olarak alınan n=3n=3 günlük rüzgar hızı ölçümü X1,X2,X3X_1, X_2, X_3 olarak kaydedilmiştir. Bilinmeyen θ\theta parametresini tahmin etmek için başlangıçta yansız bir tahmin edici olan T=2X1T = 2X_1 önerilmiştir.

Ancak, Rao-Blackwell teoremi kullanılarak tam ve yeterli (complete and sufficient) istatistik olan S=max(X1,X2,X3)S = \max(X_1, X_2, X_3) yardımıyla bu tahmin edicinin iyileştirilebileceği fark edilmiştir.

Buna göre, TT tahmin edicisinin SS yeterli istatistiğine göre koşullu beklenen değeri (E[TS]E[T|S]) alınarak elde edilen gelişmiş (Rao-Blackwell) tahmin edici aşağıdakilerden hangisidir?

Cevabı ve açıklamayı göster

Cevap: 43S\frac{4}{3} S

Cevap

43S\frac{4}{3} S
Rao-Blackwell teoremine göre, yansız bir tahmin edicinin (T=2X1T=2X_1) tam ve yeterli istatistiğe (SS) göre koşullu beklenen değeri (T=E[TS]T^* = E[T|S]), yansız ve minimum varyanslı bir tahmin edici (UMVUE) verir. Düzgün dağılımda S=max(X1,X2,X3)S = \max(X_1, X_2, X_3) tam ve yeterli istatistiktir ve E[S]=34θE[S] = \frac{3}{4}\theta'dır. SS'yi yansız bir tahmin edici yapmak için 43\frac{4}{3} katsayısı ile çarpmak gerekir. Lehmann-Scheffe teoremi gereği bu tek yansız fonksiyon, aynı zamanda Rao-Blackwell tahmin edicisinin (E[TS]E[T|S]) kendisine eşittir.

Adım Adım Çözüm

1
T=2X1T = 2X_1 tahmin edicisinin yansız olup olmadığını kontrol etme
E[2X1]=2E[X1]=2(θ2)=θE[2X_1] = 2E[X_1] = 2(\frac{\theta}{2}) = \theta elde edilir. Demek ki TT tahmin edicisi yansızdır.
Rao-Blackwell teoreminin yansızlık özelliğini koruduğunu teyit etmek için.
2
Maksimum sıra istatistiğinin (SS) beklenen değerini bulma
n=3n=3 hacimli (0,θ)(0, \theta) aralığındaki düzgün dağılım için E[S]=nn+1θ=34θE[S] = \frac{n}{n+1}\theta = \frac{3}{4}\theta olarak hesaplanır.
Yeni tahmin edicinin SS'nin bir fonksiyonu olacağını ve yansız olması gerektiğini Lehmann-Scheffe teoreminden bildiğimiz için.
3
SS istatistiğini yansız bir tahmin ediciye dönüştürme
E[cS]=θE[c \cdot S] = \theta denklemini çözecek cc katsayısı aranır. c34θ=θ    c=43c \cdot \frac{3}{4}\theta = \theta \implies c = \frac{4}{3} bulunur.
Tam ve yeterli bir istatistiğin yansız fonksiyonu tektir ve bu aynı zamanda Rao-Blackwell tahmin edicisinin ta kendisidir.
4
Gelişmiş Rao-Blackwell tahmin edicisini (TT^*) yazma
T=43ST^* = \frac{4}{3} S sonucuna ulaşılır.
Tüm koşullar sağlandığı için UMVUE formülü belirlenmiş olur.

Anahtar Kavram

Rao-Blackwell Teoremi ve Lehmann-Scheffe Uygulaması
Soru 95Soru

Bir kamu laboratuvarında görevli bir gıda mühendisi, yeni geliştirilen bir gıda koruyucusunun belirli bir standart sıcaklıktaki etkinlik süresini (saat cinsinden) incelemektedir. Etkinlik süresini temsil eden rastgele değişkenin, θ\theta bilinmeyen bir parametre olmak üzere, (θ,θ+5)(\theta, \theta+5) aralığında sürekli düzgün (uniform) dağılıma sahip olduğu varsayılmaktadır. Mühendis, birbirinden bağımsız nn adet ölçüm yaparak X1,X2,,XnX_1, X_2, \dots, X_n rastgele örneklemini elde etmiştir.

Buna göre, Neyman Çarpanlara Ayırma Teoremi dikkate alındığında, θ\theta parametresi için yeterli (sufficient) istatistik aşağıdakilerden hangisidir?

(Not: X(1)=min(X1,,Xn)X_{(1)} = \min(X_1, \dots, X_n) ve X(n)=max(X1,,Xn)X_{(n)} = \max(X_1, \dots, X_n) olarak tanımlanmıştır.)

Cevabı ve açıklamayı göster

Cevap: (X(1),X(n))(X_{(1)}, X_{(n)})

Cevap

θ\theta parametresi için yeterli istatistik (X(1),X(n))(X_{(1)}, X_{(n)}) vektörüdür.
Neyman Çarpanlara Ayırma Teoremi'ne göre, olabilirlik fonksiyonu L(θ;x)=g(T(x)θ)h(x)L(\theta; \mathbf{x}) = g(T(\mathbf{x})|\theta) \cdot h(\mathbf{x}) biçiminde yazıldığında, T(x)T(\mathbf{x}) yeterli istatistiktir. Dağılım (θ,θ+5)(\theta, \theta+5) aralığında olduğu için olabilirlik fonksiyonu ancak tüm veriler bu aralıkta olduğunda sıfırdan farklıdır. Bu koşul, örneklemin en küçük değerinin θ\theta'dan büyük veya eşit (X(1)θX_{(1)} \ge \theta) ve örneklemin en büyük değerinin θ+5\theta+5'ten küçük veya eşit (X(n)θ+5X_{(n)} \le \theta+5) olmasıyla sağlanır. Parametreye bağlı bu şart g(T(x)θ)g(T(\mathbf{x})|\theta) çarpanını oluşturur ve bu fonksiyon veriye yalnızca X(1)X_{(1)} ve X(n)X_{(n)} üzerinden bağlı olduğu için yeterli istatistik bu iki sıra istatistiğinin oluşturduğu vektördür.

Adım Adım Çözüm

1
Sürekli düzgün dağılım için olasılık yoğunluk fonksiyonunu (pdf) gösterge (indicator) fonksiyonu ile tanımlayın.
f(x;θ)=15I(θxθ+5)f(x; \theta) = \frac{1}{5} I(\theta \le x \le \theta+5)
Tanım aralığı (destek kümesi) bilinmeyen parametreye bağlı olduğunda, yeterli istatistiği bulabilmek için sınır koşullarını matematiksel olarak ifade etmek zorunludur.
2
nn çaplı örneklem için olabilirlik (likelihood) fonksiyonunu oluşturun.
L(θ;x)=i=1n15I(θxiθ+5)=(15)ni=1nI(θxiθ+5)L(\theta; \mathbf{x}) = \prod_{i=1}^n \frac{1}{5} I(\theta \le x_i \le \theta+5) = (\frac{1}{5})^n \prod_{i=1}^n I(\theta \le x_i \le \theta+5)
Örneklemlerin bağımsızlığından dolayı ortak dağılım fonksiyonu, marjinal olasılık fonksiyonlarının çarpımıdır.
3
Gösterge fonksiyonlarının çarpımını, sıra istatistiklerini kullanarak tek bir ifadeye dönüştürün.
Tüm xix_i'lerin θ\theta ile θ+5\theta+5 arasında olması, en küçük değerin θ\theta'dan büyük ve en büyük değerin θ+5\theta+5'ten küçük olması demektir: L(θ;x)=(15)nI(θX(1))I(X(n)θ+5)L(\theta; \mathbf{x}) = (\frac{1}{5})^n I(\theta \le X_{(1)}) I(X_{(n)} \le \theta+5)
Ortak gösterge bölgesini parametre ve veri ayrımı yapabilmek için basitleştirmek gereklidir.
4
Neyman Çarpanlara Ayırma Teoremi'ni (L(θ;x)=g(T(x)θ)h(x)L(\theta; \mathbf{x}) = g(T(\mathbf{x})|\theta) \cdot h(\mathbf{x})) uygulayın.
g(T(x)θ)=I(θX(1))I(X(n)θ+5)g(T(\mathbf{x})|\theta) = I(\theta \le X_{(1)}) I(X_{(n)} \le \theta+5) ve h(x)=(15)nh(\mathbf{x}) = (\frac{1}{5})^n olarak belirlenir. gg fonksiyonu veriye sadece (X(1),X(n))(X_{(1)}, X_{(n)}) üzerinden bağlıdır.
Teoreme göre, olabilirlik fonksiyonunun parametreye bağlı çarpanı içindeki veri fonksiyonu T(x)T(\mathbf{x}), yeterli istatistiktir.

Anahtar Kavram

Sınırları parametreye bağlı dağılımlarda (düzgün dağılım gibi) yeterli istatistik bulunurken gösterge (indicator) fonksiyonlarının çarpımından doğan sıra istatistiklerinin (minimum ve maksimum) kullanımı.
Soru 96Soru

Bir tarım uzmanı, yeni bir gübre türünün tohumların çimlenme olasılığı (pp) üzerindeki etkisini incelemektedir. Bu amaçla rastgele seçilen 55 adet tohum ekilmiş ve her birinin çimlenme durumu (Xi=1X_i=1: çimlenme gerçekleşti, Xi=0X_i=0: çimlenme gerçekleşmedi) kaydedilmiştir. Gözlemlerin pp parametreli Bernoulli dağılımına sahip bağımsız ve özdeş rastgele değişkenler olduğu varsayılmaktadır. Parametre için T=X1T = X_1 yansız tahmin edicisi tanımlanmıştır.

Rao-Blackwell teoremi uygulanarak S=i=15XiS = \sum_{i=1}^5 X_i yeterli istatistiği yardımıyla elde edilen yeni tahmin edicinin (T=E[TS]T^* = E[T|S]) varyansının, başlangıçtaki TT tahmin edicisinin varyansına oranı kaçtır?

Cevabı ve açıklamayı göster

Cevap: 1/51/5

Cevap

Rao-Blackwell teoremi sonucunda elde edilen yeni tahmin edicinin varyansının orijinal tahmin ediciye oranı 1/51/5 olur.
Doğru cevap 1/51/5 oranıdır. X1,,X5X_1, \dots, X_5 i.i.d. Bernoulli(pp) olduğunda, pp için T=X1T=X_1 yansızdır ve S=XiS=\sum X_i yeterli istatistiktir. Rao-Blackwell teoremi gereği T=E[X1S]=S/5=XˉT^* = E[X_1|S] = S/5 = \bar{X} elde edilir. Xˉ\bar{X} tahmin edicisinin varyansı p(1p)/5p(1-p)/5 iken orijinal tahmin edici X1X_1'in varyansı p(1p)p(1-p)'dir. Bu iki değerin oranı 1/51/5 sonucunu verir.

Adım Adım Çözüm

1
Başlangıç tahmin edicisinin (TT) varyansını belirle.
Var(T)=Var(X1)=p(1p)Var(T) = Var(X_1) = p(1-p)
X1X_1 bir Bernoulli rastgele değişkeni olduğu için varyansı p(1p)p(1-p) formülü ile hesaplanır.
2
Rao-Blackwell teoremini uygulayarak yeni tahmin ediciyi (TT^*) bul.
T=E[X1i=15Xi]=15i=15Xi=XˉT^* = E[X_1 | \sum_{i=1}^5 X_i] = \frac{1}{5} \sum_{i=1}^5 X_i = \bar{X}
Bernoulli dağılımında, bir gözlemin toplam (yeterli istatistik) üzerindeki koşullu beklenen değeri örneklem ortalamasına eşittir.
3
Yeni tahmin edicinin (TT^*) varyansını hesapla.
Var(T)=Var(Xˉ)=p(1p)5Var(T^*) = Var(\bar{X}) = \frac{p(1-p)}{5}
Bağımsız ve özdeş dağılan rastgele değişkenlerin ortalamasının varyansı, kitle varyansının örneklem çapına (n=5n=5) bölünmesiyle bulunur.
4
Varyanslar arası oranı hesapla.
Var(T)Var(T)=p(1p)/5p(1p)=15\frac{Var(T^*)}{Var(T)} = \frac{p(1-p)/5}{p(1-p)} = \frac{1}{5}
Yeni tahmin edicinin varyansı ile orijinal varyans oranlanarak iyileşme miktarı bulunur.

Anahtar Kavram

Rao-Blackwell Teoremi ve Varyans Azaltma
Tahmini Süre:1m 30s
Soru 97Soru

Bir vergi denetmeni, bir bölgedeki yüksek gelirli mükelleflerin beyan ettikleri yıllık kazançları (bin TL) incelemektedir. Kazançların (XX), bilinen bir x0x_0 (xx0x \geq x_0) alt sınırı için aşağıdaki Pareto olasılık yoğunluk fonksiyonuna uygun olduğu varsayılmaktadır:

f(x;θ)=θx0θxθ+1,xx0,θ>0f(x; \theta) = \frac{\theta x_0^\theta}{x^{\theta+1}}, \quad x \geq x_0, \theta > 0

Rastgele seçilen nn adet mükellefin kazançları X1,X2,,XnX_1, X_2, \dots, X_n olduğuna göre, Neyman Çarpanlara Ayırma Teoremi'ne göre θ\theta parametresi için yeterli istatistik aşağıdakilerden hangisidir?

Cevabı ve açıklamayı göster

Cevap: i=1nXi\prod_{i=1}^{n} X_i

Cevap

Gözlemlerin çarpımı olan i=1nXi\prod_{i=1}^{n} X_i ifadesi, θ\theta parametresi için yeterli istatistiktir.
Olabilirlik fonksiyonu incelendiğinde, θ\theta parametresini içeren çarpanın verilere sadece gözlemlerin çarpımı (Xi\prod X_i) aracılığıyla bağlı olduğu görülmektedir. Neyman Çarpanlara Ayırma Teoremi'ne göre bu durum, gözlemlerin çarpımının parametre hakkında tam bilgi sağladığını ve dolayısıyla yeterli istatistik olduğunu kanıtlar.

Adım Adım Çözüm

1
Olabilirlik fonksiyonunu (L(θ)L(\theta)) kurun.
L(θ;x1,,xn)=i=1nθx0θxiθ+1=θn(x0θ)n(i=1nxi)θ+1L(\theta; x_1, \dots, x_n) = \prod_{i=1}^n \frac{\theta x_0^\theta}{x_i^{\theta+1}} = \frac{\theta^n (x_0^\theta)^n}{(\prod_{i=1}^n x_i)^{\theta+1}}
Yeterli istatistiği bulmanın ilk adımı örneklemin ortak olasılık yoğunluk fonksiyonunu yazmaktır.
2
Fonksiyonu parametreye bağlı ve bağımsız kısımlara ayırın.
L(θ;x)=[θnx0nθ(i=1nxi)θ][(i=1nxi)1]L(\theta; \mathbf{x}) = \left[ \theta^n x_0^{n\theta} \left( \prod_{i=1}^n x_i \right)^{-\theta} \right] \cdot \left[ \left( \prod_{i=1}^n x_i \right)^{-1} \right]
Neyman Çarpanlara Ayırma Teoremi, olabilirlik fonksiyonunun g(T(x),θ)h(x)g(T(\mathbf{x}), \theta) \cdot h(\mathbf{x}) şeklinde yazılmasını gerektirir.
3
Yeterli istatistiği (T(x)T(\mathbf{x})) belirleyin.
g(T(x),θ)=θnx0nθ(xi)θ    T(x)=i=1nXig(T(\mathbf{x}), \theta) = \theta^n x_0^{n\theta} (\prod x_i)^{-\theta} \implies T(\mathbf{x}) = \prod_{i=1}^n X_i
Parametreye (θ\theta) bağlı olan kısım sadece gözlemlerin çarpımı (TT) üzerinden verilere bağlıdır.

Anahtar Kavram

Neyman Çarpanlara Ayırma Teoremi uyarınca, olabilirlik fonksiyonu parametreye bağlı bir istatistik ve parametreden bağımsız bir gözlem fonksiyonu olarak çarpılabiliyorsa, o istatistik yeterlidir.

Alternatif Yöntem

Log-olabilirlik fonksiyonu üzerinden gidilerek ln(Xi)\sum \ln(X_i) ifadesinin de yeterli olduğu görülebilir; ancak şıklarda bu ifadenin karşılığı olan çarpım formu (Xi\prod X_i) yer almaktadır.
Tahmini Süre:1m 30s
Soru 98Soru

pp başarı olasılığına sahip bir Bernoulli kitlesinden n=2n=2 birimlik bağımsız bir rastgele örneklem (X1,X2X_1, X_2) seçilmiştir. pp parametresini tahmin etmek için T=X1T = X_1 yansız tahmin edicisi kullanılmaktadır. S=X1+X2S = X_1 + X_2 istatistiğinin pp parametresi için yeterli olduğu bilinmektedir.

Rao-Blackwell teoremi uygulanarak elde edilen geliştirilmiş tahmin edici T=E[TS]T^* = E[T|S] için varyans değeri aşağıdakilerden hangisidir?

Cevabı ve açıklamayı göster

Cevap: p(1p)2\frac{p(1-p)}{2}

Cevap

Geliştirilmiş tahmin edicinin varyansı p(1p)2\frac{p(1-p)}{2} ifadesine eşittir.
Rao-Blackwell teoremi uyarınca T=E[X1X1+X2]T^* = E[X_1 | X_1 + X_2] işlemi yapıldığında elde edilen tahmin edici Xˉ=(X1+X2)/2\bar{X} = (X_1+X_2)/2 olur. Bağımsız Bernoulli değişkenleri için örneklem ortalamasının varyansı Var(Xˉ)=p(1p)nVar(\bar{X}) = \frac{p(1-p)}{n} formülü ile hesaplanır; burada n=2n=2 olduğu için sonuç p(1p)2\frac{p(1-p)}{2} olur.

Adım Adım Çözüm

1
Orijinal tahmin edicinin varyansını hesapla
Var(T)=Var(X1)=p(1p)Var(T) = Var(X_1) = p(1-p)
Bernoulli dağılımında varyans p(1p)p(1-p) formülü ile hesaplanır.
2
Rao-Blackwell tahmin edicisini (TT^*) belirle
T=E[X1X1+X2]=X1+X22T^* = E[X_1 | X_1 + X_2] = \frac{X_1 + X_2}{2}
Simetri özelliğinden dolayı, toplam verildiğinde her bir gözlemin beklenen değeri ortalamaya eşit olur.
3
TT^* tahmin edicisinin varyansını hesapla
Var(T)=Var(X1+X22)=14[Var(X1)+Var(X2)]=2p(1p)4=p(1p)2Var(T^*) = Var\left(\frac{X_1 + X_2}{2}\right) = \frac{1}{4} [Var(X_1) + Var(X_2)] = \frac{2p(1-p)}{4} = \frac{p(1-p)}{2}
Bağımsız değişkenlerin toplamının varyansı, varyansların toplamına eşittir ve sabit çarpan dışarı karesi alınarak çıkar.

Anahtar Kavram

Rao-Blackwell teoremi, yansız bir tahmin edicinin yeterli bir istatistik üzerine koşullandırılmasıyla varyansın azaltılabileceğini (veya aynı kalacağını) ifade eder.
Soru 99Soru

Bir elektronik fabrikasında üretilen mikroçiplerin belirli bir yükleme altındaki maksimum çalışma gerilimi XX (volt cinsinden) rassal değişkeni ile modellenmektedir. XX değişkeninin olasılık yoğunluk fonksiyonu, θ>0\theta > 0 parametresi için aşağıda verilmiştir:

f(x;θ)={2(θx)θ2,0<x<θ0,dig˘er durumlarf(x; \theta) = \begin{cases} \frac{2(\theta - x)}{\theta^2}, & 0 < x < \theta \\ 0, & \text{diğer durumlar} \end{cases}

Bu kitle içinden seçilen nn birimlik bir örneklemin ortalaması olan Xˉ\bar{X} istatistiği, θ\theta parametresini tahmin etmek amacıyla kullanılmaktadır. Buna göre, Xˉ\bar{X} tahmin edicisinin hata kareler ortalaması (MSE) aşağıdakilerden hangisidir?

Cevabı ve açıklamayı göster

Cevap: (8n+1)θ218n\frac{(8n+1)\theta^2}{18n}

Cevap

(8n+1)θ218n\frac{(8n+1)\theta^2}{18n}
Yapılan hesaplamalar sonucunda kitlenin varyansının θ2/18\theta^2/18, tahmin edicinin yan miktarının ise 2θ/3-2\theta/3 olduğu görülmüştür. MSE=Var+Bias2MSE = Var + Bias^2 formülüne göre bu değerler yerine konulduğunda ((8n+1)θ2)/(18n)((8n+1)\theta^2)/(18n) sonucuna ulaşılır.

Adım Adım Çözüm

1
Kitlenin beklenen değerini (E[X]E[X]) hesaplayınız.
E[X]=0θx2(θx)θ2dx=2θ2[θx22x33]0θ=2θ2(θ32θ33)=θ3E[X] = \int_{0}^{\theta} x \frac{2(\theta - x)}{\theta^2} dx = \frac{2}{\theta^2} \left[ \frac{\theta x^2}{2} - \frac{x^3}{3} \right]_0^{\theta} = \frac{2}{\theta^2} \left( \frac{\theta^3}{2} - \frac{\theta^3}{3} \right) = \frac{\theta}{3}
Yan (bias) miktarını belirlemek için öncelikle tahmin edicinin beklenen değeri (E[Xˉ]=E[X]E[\bar{X}] = E[X]) bulunmalıdır.
2
Kitlenin varyansını (Var(X)Var(X)) hesaplayınız.
E[X2]=0θx22(θx)θ2dx=2θ2[θx33x44]0θ=θ26E[X^2] = \int_{0}^{\theta} x^2 \frac{2(\theta - x)}{\theta^2} dx = \frac{2}{\theta^2} \left[ \frac{\theta x^3}{3} - \frac{x^4}{4} \right]_0^{\theta} = \frac{\theta^2}{6}
Var(X)=E[X2](E[X])2=θ26(θ3)2=3θ22θ218=θ218Var(X) = E[X^2] - (E[X])^2 = \frac{\theta^2}{6} - \left(\frac{\theta}{3}\right)^2 = \frac{3\theta^2 - 2\theta^2}{18} = \frac{\theta^2}{18}
MSE hesaplamasında kullanılacak olan varyans bileşenini bulmak gerekir.
3
Xˉ\bar{X} tahmin edicisinin yan (bias) miktarını ve varyansını bulunuz.
Bias(Xˉ)=E[Xˉ]θ=θ3θ=2θ3Bias(\bar{X}) = E[\bar{X}] - \theta = \frac{\theta}{3} - \theta = -\frac{2\theta}{3}
Var(Xˉ)=Var(X)n=θ218nVar(\bar{X}) = \frac{Var(X)}{n} = \frac{\theta^2}{18n}
Örneklem ortalamasının özellikleri kitle özelliklerinden türetilir.
4
Hata Kareler Ortalaması (MSE) formülünü uygulayınız.
MSE(Xˉ)=Var(Xˉ)+[Bias(Xˉ)]2=θ218n+(2θ3)2=θ218n+4θ29=θ2+8nθ218n=(8n+1)θ218nMSE(\bar{X}) = Var(\bar{X}) + [Bias(\bar{X})]^2 = \frac{\theta^2}{18n} + \left(-\frac{2\theta}{3}\right)^2 = \frac{\theta^2}{18n} + \frac{4\theta^2}{9} = \frac{\theta^2 + 8n\theta^2}{18n} = \frac{(8n+1)\theta^2}{18n}
MSE, tahmin edicinin hem değişkenliğini (varyans) hem de doğruluğunu (yan) ölçen birleşik bir kriterdir.

Anahtar Kavram

Hata Kareler Ortalaması (MSE) = Varyans + (Yan)^2

Daha Fazla Pratik

Farklı bir tahmin edici olarak T=3XˉT = 3\bar{X} istatistiğinin yansız olup olmadığını kontrol edip MSE değerini hesaplayarak karşılaştırma yapabilirsiniz.
Tahmini Süre:2m 0s
Soru 100Soru

Bir barajdaki günlük su buharlaşma miktarı (milimetre cinsinden) XX rassal değişkeni, θ>0\theta > 0 olmak üzere aşağıdaki olasılık yoğunluk fonksiyonu ile modellenmektedir:

f(x;θ)={4x3θ4,0<x<θ0,dig˘er durumlarf(x; \theta) = \begin{cases} \frac{4x^3}{\theta^4}, & 0 < x < \theta \\ 0, & \text{diğer durumlar} \end{cases}

Bu kitleden çekilen nn birimlik bir rassal örneklem için örneklem ortalaması Xˉ\bar{X} tahmin edicisinin θ\theta parametresi için Hata Kareler Ortalaması (MSE) aşağıdakilerden hangisidir?

Cevabı ve açıklamayı göster

Cevap: (3n+2)θ275n\frac{(3n+2)\theta^2}{75n}

Cevap

Örneklem ortalaması Xˉ\bar{X} tahmin edicisinin θ\theta parametresi için Hata Kareler Ortalaması (3n+2)θ275n\frac{(3n+2)\theta^2}{75n} olarak bulunur.
Tahmin edicinin Hata Kareler Ortalaması (MSE), varyansı (2θ2/75n2\theta^2/75n) ile yan miktarının karesinin ([θ/5]2=θ2/25[-\theta/5]^2 = \theta^2/25) toplamıdır. Bu iki terim toplandığında payda eşitleme sonucu (3n+2)θ2/75n(3n+2)\theta^2/75n ifadesine ulaşılır.

Adım Adım Çözüm

1
XX rassal değişkeninin beklenen değerini (E[X]E[X]) hesaplayın.
E[X]=0θx4x3θ4dx=4θ4[x55]0θ=4θ5E[X] = \int_{0}^{\theta} x \cdot \frac{4x^3}{\theta^4} dx = \frac{4}{\theta^4} \left[ \frac{x^5}{5} \right]_{0}^{\theta} = \frac{4\theta}{5}
Tahmin edicinin beklenen değerini ve yan (bias) miktarını bulmak için kitlenin beklenen değerine ihtiyaç vardır.
2
XX rassal değişkeninin varyansını (Var(X)Var(X)) hesaplayın.
E[X2]=0θx24x3θ4dx=4θ4[x66]0θ=2θ23E[X^2] = \int_{0}^{\theta} x^2 \cdot \frac{4x^3}{\theta^4} dx = \frac{4}{\theta^4} \left[ \frac{x^6}{6} \right]_{0}^{\theta} = \frac{2\theta^2}{3} olduğundan; Var(X)=2θ23(4θ5)2=2θ2316θ225=50θ248θ275=2θ275Var(X) = \frac{2\theta^2}{3} - \left( \frac{4\theta}{5} \right)^2 = \frac{2\theta^2}{3} - \frac{16\theta^2}{25} = \frac{50\theta^2 - 48\theta^2}{75} = \frac{2\theta^2}{75}
Tahmin edicinin varyansını hesaplamak için kitlenin varyansının bilinmesi gerekir.
3
Tahmin edicinin beklenen değerini ve yan (bias) miktarını hesaplayın.
E[Xˉ]=E[X]=4θ5E[\bar{X}] = E[X] = \frac{4\theta}{5}
Yan(Xˉ)=E[Xˉ]θ=4θ5θ=θ5Yan(\bar{X}) = E[\bar{X}] - \theta = \frac{4\theta}{5} - \theta = -\frac{\theta}{5}
MSE formülünde yer alan yan (bias) bileşenini belirlemek için gereklidir.
4
Tahmin edicinin varyansını hesaplayın.
Var(Xˉ)=Var(X)n=2θ275nVar(\bar{X}) = \frac{Var(X)}{n} = \frac{2\theta^2}{75n}
MSE formülündeki varyans bileşenini belirlemek için gereklidir.
5
MSE=Var+Yan2MSE = Var + Yan^2 formülünü kullanarak sonucu bulun.
MSE(Xˉ)=2θ275n+(θ5)2=2θ275n+θ225MSE(\bar{X}) = \frac{2\theta^2}{75n} + \left( -\frac{\theta}{5} \right)^2 = \frac{2\theta^2}{75n} + \frac{\theta^2}{25}
MSE(Xˉ)=2θ2+3nθ275n=(3n+2)θ275nMSE(\bar{X}) = \frac{2\theta^2 + 3n\theta^2}{75n} = \frac{(3n+2)\theta^2}{75n}
Hata Kareler Ortalaması, varyans ve yan miktarının karesinin toplamına eşittir.

Anahtar Kavram

Hata Kareler Ortalaması (MSE), bir tahmin edicinin hem değişkenliğini (varyans) hem de doğruluğunu (yansızlık/yan) aynı anda ölçen pedagojik bir ölçüttür.

Daha Fazla Pratik

Tahmin edicinin yansız olması için örneklem ortalamasının hangi katsayı ile çarpılması gerektiğini düşünerek UMVUE kavramını araştırabilirsiniz.
Tahmini Süre:2m 30s
ÖncekiSayfa 5 / 12Sonraki
Parametre Tahmin Yöntemleri Alıştırma Soruları — KPSS İstatistik — Sayfa 5 | Examkin