Examkin

Soru

Zorluk: Çok zorBonferroni ve Sidak Düzeltmeleri

Bir kamu araştırma enstitüsünde görevli başistatistikçi, kuruma ait geniş çaplı bir veritabanı üzerinde 5050 farklı sosyoekonomik göstergeyi bağımsız örneklemler üzerinden karşılaştıran eş zamanlı hipotez testleri yürütecektir. Deney genelinde Grup Bazında Tip I Hata Oranının (FWER) α=0,05\alpha = 0,05 seviyesini aşmaması resmi bir gerekliliktir. Uzman ekip, kullanılacak çoklu karşılaştırma (post-hoc) prosedürü için Bonferroni ve Šidák düzeltmeleri üzerine aşağıdaki değerlendirmeleri yapmıştır:

I. Šidák düzeltmesi kullanıldığında her bir test için hesaplanan nominal anlamlılık düzeyi (αS\alpha_S), Bonferroni düzeltmesi ile hesaplanan değerden (αB\alpha_B) daima daha büyük olacaktır. Bu matematiksel eşitsizlik, Šidák yöntemini Bonferroni'ye göre testin gücünü artırma ve Tip II hata riskini azaltma açısından teorik olarak daha avantajlı kılar.
II. Bonferroni yöntemi yalnızca testlerin birbirinden tam bağımsız olduğu varsayımı altında FWER kontrolünü garanti edebilirken; Šidák yöntemi değişkenler arasındaki her türlü korelasyon (bağımlılık) yapısında Tip I hata oranının belirlenen seviyede kalmasını kesin olarak sağlar.
III. Eğer herhangi bir çoklu karşılaştırma düzeltmesi yapılmaksızın her bir gösterge için ayrı ayrı nominal α=0,05\alpha = 0,05 düzeyinde bağımsız testler gerçekleştirilirse, araştırma genelinde en az bir göstergede yanlışlıkla anlamlı fark bulunması olasılığı %90\%90'ın üzerine çıkacaktır.

Buna göre, uzman ekibin yaptığı değerlendirmelerden hangileri istatistiksel olarak kesinlikle doğrudur?

  1. A
    Yalnız I
  2. B
    Yalnız III
  3. C
    I ve II
  4. I ve IIICevap
  5. E
    I, II ve III

Cevap

Uzman ekibin değerlendirmelerinden I ve III numaralı ifadeler kesinlikle doğrudur.
I ve III numaralı ifadeleri içeren seçenek doğru cevaptır. Šidák düzeltmesi matematiksel yapısı gereği daima Bonferroni'den daha büyük bir anlamlılık düzeyi üretir, bu da testin gücünü nispeten artırır (I. ifade doğru). Düzeltme yapılmadan α=0,05\alpha=0,05 ile 5050 bağımsız test yapıldığında deney bazında hata olasılığı 1−(0,95)50≈%92,31 - (0,95)^{50} \approx \%92,3 olarak hesaplanır, ki bu %90\%90'ın üzerindedir (III. ifade doğru). Ancak Bonferroni bağımsızlık varsayımına ihtiyaç duymayan, Šidák ise kesinlikle bağımsızlık varsayan bir yöntemdir; dolayısıyla II. ifade teorik olarak tamamen hatalıdır.

Adım Adım Çözüm

1
I. ifadenin matematiksel analizini yapın.
Bonferroni için αB=αm\alpha_B = \frac{\alpha}{m}, Šidák için αS=1−(1−α)1/m\alpha_S = 1 - (1-\alpha)^{1/m} formülleri geçerlidir. m=50m=50 ve α=0,05\alpha=0,05 için αB=0,001\alpha_B = 0,001 iken, αS≈0,001025\alpha_S \approx 0,001025 bulunur. αS>αB\alpha_S > \alpha_B olduğundan, Šidák yöntemi null hipotezini reddetmeyi marjinal olarak daha kolay hale getirir, bu da gücü artırır ve Tip II hatayı azaltır. I. ifade doğrudur.
Yöntemlerin test güçlerini karşılaştırmak için hesapladıkları yeni eşik değerlerinin büyüklük sıralamasını bilmek gerekir.
2
II. ifadedeki teorik varsayımları inceleyin.
Bonferroni yöntemi Boole eşitsizliğine (P(∪Ai)≤∑P(Ai)P(\cup A_i) \leq \sum P(A_i)) dayandığı için testler arasındaki korelasyon veya bağımlılık yapısından bağımsız olarak daima FWER ≤α\leq \alpha garantisi verir. Šidák yöntemi ise P(HataYok)=(1−α)mP(Hata Yok) = (1-\alpha)^m eşitliğini kullandığı için testlerin tam bağımsız olması varsayımına dayanır. II. ifade durumu tam tersi yansıttığı için yanlıştır.
Post-hoc testlerinin veri yapısındaki bağımlılıklara karşı ne kadar dirençli olduklarını belirleyen temel istatistiksel varsayımları doğrulamak gerekir.
3
III. ifadedeki birikimli hata olasılığını (Düzeltilmemiş FWER) hesaplayın.
Hiçbir düzeltme yapılmadığında en az bir Tip I hata yapma olasılığı 1−(1−α)m1 - (1 - \alpha)^m formülü ile bulunur. 1−(1−0,05)50=1−(0,95)501 - (1 - 0,05)^{50} = 1 - (0,95)^{50} hesabı yaklaşık olarak 1−0,0769=0,92311 - 0,0769 = 0,9231 (\%92,31) değerini verir. Dolayısıyla hata olasılığı %90\%90'ın üzerine çıkacaktır. III. ifade doğrudur.
Çoklu karşılaştırmalarda alfa enflasyonunun ulaştığı boyutu matematiksel olarak kanıtlamak gerekir.

Anahtar Kavram

FWER Kontrolü ve Yöntem Varsayımları
Bu soruyu puanla

Konular

Bonferroni ve Sidak Düzeltmeleri

Bu konudan daha fazla soru