Kısa cevap "evet yapabilirsiniz" - ancak "büyük model" in Maksimum Olabilirlik Tahminlerini (MLE) her ikisine de takılan her iki modeldeki tüm ortak değişkenlerle karşılaştırmalısınız.
Bu, sorunuzu cevaplamak için olasılık teorisini elde etmenin "yarı resmi" bir yoludur
Örnekte, ve aynı tip değişkenlerdir (kesirler / yüzdeler) ve bu nedenle karşılaştırılabilirler. Aynı modeli her ikisine de uyduğunuzu varsayacağım. Yani iki modelimiz var: Y 2Y1Y2
l o g ( p 1 i
M1: Y1 ben~ B ı n ( n-1 ben, p1 ben)
l o g( p1 ben1 - p1 ben) = α1+ β1Xben
M2: Y2i~ B i n (n-2 i,p2 i)
l o g( p2 i1 - p2 i) = α2+ β2Xben
Böylece değerlendirmek istediğiniz hipoteziniz var:
'H0: β1> β2
Ve bazı verileriniz ve bazı önceki bilgileriniz (lojistik model kullanımı gibi) var. Yani olasılığı hesaplıyorsunuz:{ Y1 ben, Y2 i, Xben}ni = 1
P= Pr ( H0| { Y1 ben, Y2 i, Xben}ni = 1, Ben)
Şimdi , regresyon parametrelerinin hiçbirinin gerçek değerine bağlı değildir, bu nedenle marjinalleştirilerek kaldırılması gerekir.'H0
P= ∫∞- ∞∫∞- ∞∫∞- ∞∫∞- ∞Pr ( H0, α1, α2, β1, β2| { Y1 ben, Y2 i, Xben}ni = 1, Ben) dα1dα2dβ1dβ2
Hipotez, entegrasyon aralığını basitçe kısıtlar, bu yüzden:
P= ∫∞- ∞∫∞β2∫∞- ∞∫∞- ∞Pr ( α1, α2, β1, β2| { Y1 ben, Y2 i, Xben}ni = 1, Ben) dα1dα2dβ1dβ2
Olasılık veriler üzerinde koşullu olduğu için, her model için iki ayrı posterior içerecektir.
Pr ( α1, β1| { Y1 ben, Xben, Y2 i}ni = 1, Ben) Pr ( α2, β2| { Y2 i, Xben, Y1 ben}ni = 1, Ben)
Şimdi ve arasında doğrudan bir bağlantı olmadığı için , sadece bilinen üzerinden dolaylı bağlantılar, ikinci çıkacaktır. aynı ilk posteriordaki için de .Y1 benα2, β2XiY2i
Standart lojistik regresyon teorisinden ve tekdüze önceki olasılıklar varsayıldığında, parametreler için posterior, MLE'lere eşit ortalama ile V- ve ile belirtilen bilgi matrisine eşit varyans ile yaklaşık iki değişkenli normaldir - parametrelere bağlı olmayan, sadece MLE'ler. böylece bilinen varyans matrisi ile düz ileriye doğru normal integrallere sahipsiniz. herhangi bir katkı olmaksızın marjinalleşir (diğer herhangi bir "ortak değişken" gibi) ve olağan sonuçla kalırız (isterseniz türetmenin ayrıntılarını gönderebilirim, ancak oldukça "standart" şeyler):V1V2αj
P=Φ(β^2,MLE−β^1,MLEV1:β,β+V2:β,β−−−−−−−−−−−√)
Nerede sadece standart normal CDF olduğunu. Bu normal araç testinin olağan karşılaştırmasıdır. Ancak bu yaklaşımın her birinde aynı regresyon değişkenleri setinin kullanılmasını gerektirdiğini unutmayın. Birçok öngörücüye sahip çok değişkenli durumda, farklı regresyon değişkenlerine sahipseniz, integraller yukarıdaki teste etkili bir şekilde eşit olacaktır, ancak her iki modelden tüm ortak değişkenleri içeren "büyük model" deki iki betanın MLE'lerinden.Φ()