Soru

Zorluk: ZorHücre Başına Tek Gözlemli Tasarım (Etkileşimsiz Model)

Bir tarım araştırmacısı, 44 farklı gübre türü (A faktörü) ve 55 farklı sulama yönteminin (B faktörü) bitki verimi üzerindeki etkilerini incelemek amacıyla iki yönlü varyans analizi (ANOVA) planlamaktadır. Zaman ve maliyet kısıtları nedeniyle, her bir gübre-sulama kombinasyonu (hücre) için yalnızca bir adet gözlem yapılabilmiştir (n=1n=1). Araştırmacı, test edilebilecek bir hata terimi kalabilmesi için analizini toplanabilir (etkileşimsiz) model üzerinden kurarak sadece ana etkileri test etmeye karar verir.

Ancak, araştırmacının bilmediği bir gerçek olarak, yığın (popülasyon) düzeyinde gübre türü ile sulama yöntemi arasında istatistiksel olarak anlamlı ve oldukça güçlü bir etkileşim etkisi bulunmaktadır.

Buna göre, popülasyonda var olan bu etkileşim etkisinin modele dahil edilememesi (toplanabilir modelin kullanılması), araştırmacının analiz sonuçlarını ve istatistiksel kararlarını nasıl etkiler?

  1. Hata kareler ortalaması (MSEMS_E) gerçekte olduğundan daha büyük tahmin edileceği için, faktörlerin ana etkilerine yönelik yapılacak F testlerinin istatistiksel gücü azalır.Cevap
  2. B
    Etkileşim etkisi, gruplar içi hata varyansından bağımsız bir kaynağa dönüşerek test edileceğinden, toplanabilir model kullanılmasına rağmen ana etkilerin anlamlılığı bu durumdan etkilenmez.
  3. C
    Kareler toplamlarının ayrıştırılmasındaki toplanabilirlik varsayımı çökeceği için, ihmal edilen etkileşim nedeniyle ana etkilerin kareler toplamları matematiksel olarak negatif değerler alır.
  4. D
    Hücre başına yalnızca tek gözlem bulunması, modeli matematiksel olarak sabit etkili modelden tesadüfi etkili (Model II) modele dönüştüreceğinden beklenen kareler ortalaması kuralları değişir.
  5. E
    Ana etkilerin serbestlik dereceleri toplamı genel serbestlik derecesini tamamen tüketeceğinden, modelde hata kareler ortalamasını (MSEMS_E) hesaplayabilmek için serbestlik derecesi kalmaz.

Cevap

Toplanabilir (etkileşimsiz) modelde, modele dahil edilmeyen etkileşim etkisi hata teriminin içerisine karışır ve hata kareler ortalamasını (MSEMS_E) gerçekte olduğundan daha büyük tahmin etmemize neden olur. Şişen hata terimi, ana etkiler için hesaplanan F istatistiklerinin küçülmesine yol açarak istatistiksel gücü azaltır.
Doğru yanıt, hata teriminin şişmesi nedeniyle F testinin gücünün düşeceğini belirten seçenektir. Hücre başına tek gözlem olan tasarımlarda etkileşimi tahmin etmeye yetecek veri bulunmadığından, modelde etkileşimin olmadığını (toplanabilirlik) varsaymak zorunludur. Ancak doğada/popülasyonda gerçek bir etkileşim varsa, modele alınmayan bu güçlü varyans kaynağı sistemde 'açıklanamayan varyans' yani hata (residual) olarak sınıflandırılır. Bu durumda hesaplanan hata kareler ortalaması (MSEMS_E), saf rastgele hatadan çok daha büyük olur. F formülünde (F=MSA/MSEF=MS_A/MS_E) bölen değer (payda) suni olarak büyüdüğü için F test değeri küçülür. Test değerinin küçülmesi, anlamlı farklılıkları tespit etme şansımızı (istatistiksel güç) azaltır.

Adım Adım Çözüm

1
Hücre başına tek gözlemli (n=1n=1) ANOVA tasarımının matematiksel özelliklerini belirlemek.
Gözlem sayısı yetersiz olduğundan etkileşim etkisini ayrı bir faktör olarak tahmin edebilecek serbestlik derecesi kalmaz. Bu nedenle araştırmacı, etkileşimin sıfır olduğunu varsayan toplanabilir modeli kurar.
Etkileşim modeline girmek için fazladan gözleme ihtiyaç vardır, aksi halde hata varyansı hesaplanamaz.
2
Popülasyonda gerçek bir etkileşim varsa toplanabilir modelin hata teriminin (MSEMS_E) nasıl etkilendiğini analiz etmek.
Araştırmacının modele koymadığı gerçek etkileşim varyansı doğrudan deneysel hata varyansının içine dahil olur. Beklenen kareler ortalaması (E(MSE)E(MS_E)), yalnızca deneysel hatayı (σ2\sigma^2) değil, hata ve etkileşim bileşenini birlikte yansıtır ve olması gerekenden yüksek çıkar.
Varyans analizi kuralı gereği açıklanamayan tüm kaynaklar hata (residual) teriminde birikir.
3
Şişmiş bir hata kareler ortalamasının (MSEMS_E) F testlerine ve analizin sonucuna etkisini değerlendirmek.
Ana etkilerin F istatistikleri, faktör kareler ortalamasının hata kareler ortalamasına bölünmesiyle (F=MSAnaEtki/MSEF = MS_{AnaEtki} / MS_E) bulunur. Paydadaki MSEMS_E gereksiz yere büyük olduğunda, bölme işleminin sonucu (F değeri) matematiksel olarak daha küçük çıkar.
Payda büyüdükçe oran küçülür.
4
Küçülen F istatistiğinin tipik hipotez testi hataları üzerindeki yansımasını yorumlamak.
Daha küçük bir F değeri, H0H_0 (Yokluk) hipotezini reddetmeyi zorlaştırır. Bu da gerçekte var olan ana etkileri tespit etme yeteneğini (istatistiksel gücü) düşürür ve yanlışlıkla H0H_0'ı kabul etme (II. Tip hata) olasılığını artırır.
İstatistiksel güç, gerçek bir etkiyi doğru şekilde tespit edebilme ihtimalidir.

Anahtar Kavram

İki Yönlü ANOVA'da Toplanabilirlik Varsayımı ve Etkileşimin Hata Terimiyle Karışması (Confounding)
Bu soruyu puanla