Bu olabilirlik oranlarını doğru hesapladım mı?

R için ez paketinin yazarıyım ve ANOVA'ların çıkışına otomatik olasılık oranlarını (LR'ler) otomatik olarak hesaplamak için bir güncelleme üzerinde çalışıyorum. Fikir, ANOVA'nın başardığı etkinin testine benzer her etki için bir LR sağlamaktır. Örneğin, ana etki için LR temel etkisi içeren bir modeli için bir boş modeli karşılaştırmasını temsil eder, bir etkileşim için LR temel etkileri de içeren bir modele göre, her iki bileşen temel etkileri içeren bir modelin karşılaştırmasını temsil eder ve etkileşimleri vb.

Şimdi, LR hesaplama benim anlayış gelen Glover & Dixon ( PDF temel hesaplamaları yanı sıra karmaşıklığı için düzeltmeler ve karşı eke kapsar), BORTOLUSSI & Dixon ( Ek PDF kapakları hesaplamaları tekrarlayan ölçümler değişkeni içeren). Anlayışımı test etmek için, dfs ve SS'leri örnek bir ANOVA'dan (sahte veriler kullanılarak 2 * 2 * 3 * 4 tasarımından üretilen) alan ve her efekt için LR'nin hesaplanması adımlarını içeren bu elektronik tabloyu geliştirdim .

Böyle bir hesaplamaya biraz daha güven duyan biri bir göz atabilir ve her şeyi doğru yaptığımdan emin olsaydı gerçekten çok memnun olurum. Soyut kodu tercih edenler için, işte ezANOVA () güncellemesini uygulayan R kodu (bkz. Özellikle satır 15-95).

SS değerlerinden LR'nin hesaplanmasıyla ilgili akıl yürütme oldukça adil olmasına rağmen, en küçük kareler yöntemi eşdeğerdir, ancak olasılık tahmini ile aynı değildir. (Fark, örneğin, en küçük kareler yaklaşımında (n-1) bölünen ve maksimum olabilirlikte n'ye bölünen se'nin hesaplanmasında gösterilebilir. ).

Bunun bazı sonuçları var: LR'yi, olasılık ile orantılı olduğundan hesaplayabilirsiniz , ancak bu, anova modelinizin kendisinin olasılığını vermez. Sadece oran hakkında bir şeyler söyler. AIC olasılık olarak klasik olarak tanımlandığından, AIC'yi istediğiniz gibi kullanıp kullanamayacağınızdan emin değilim.$\frac{1}{s}$

Ben elektronik tabloya baktım, ama içinde "düzeltilmemiş LR" değerleri (ben de tam olarak orada hesaplamak için ne çalıştığını tam olarak takip etmiyorum) bana pek olası görünmüyor.

Bir yan notta, LR testinin gücü, sadece istediğiniz modelleri kontrastlayabilmenizdir, hepsi için bunu yapmak zorunda kalmazsınız (bu da çok test hatasını azaltır). Bunu her dönem için yaparsanız, LR'niz tamamen bir F testine eşittir ve en az kareler söz konusu olduğunda, sayısal olarak bile aynı şeyi bildiğim kadarıyla.

Miliniz değişebilir, ancak iki farklı çerçevenin konseptlerini karıştırmaktan hiç emin olamadım (örn. Maksimum olasılığa karşı en küçük kareler). Şahsen, F istatistiklerini rapor eder ve LR'yi modelleri karşılaştırmaya izin veren bir işlevde uygularım (örneğin, tam olarak bunu yapan lme modelleri için anova işlevi).

Benim 2 sentim.

PS: Kodunuza baktım, ama tüm değişkenleri gerçekten anlayamadım. Kodunuzu yorumları kullanarak ek açıklama yaparsanız, bu hayatı yeniden kolaylaştırır. EXCEL sayfası da anlaşılması en kolay yöntem değildir. Ondan bir şey yapıp yapamayacağımı görmek için daha sonra tekrar kontrol edeceğim.