[問題] 關於羅吉斯-2LL

看板Statistics作者 (主席)時間9年前 (2016/07/21 16:41), 9年前編輯推噓0(009)
留言9則, 3人參與, 最新討論串1/1
各位版大好 關於羅吉斯分析中的-2LL 關於null model的-2LLo計算方法應該是如何呢? 用下方資料羅吉斯分析的結果得到 B0=-3.604 B1=0.359 B2=0.538 B3=-0.449 我僅將常數B0=-3.604,B1~B3令為零,帶入羅吉斯模型得到-2LLo=109.73 但是統計軟體跑出的值-2LLo=41.589 是不是我對null model的理解錯誤了呢? 請教各位版大了,感謝! PS 我用的公式 LL=sum(yi*ln(pi)+(1-yi)*ln(1-pi)), 其中yi為實際值(1 or 0), pi為羅吉斯預測機率 --------------------------------------------------------- 編號 Y X1 X2 X3 1 1 13.6 10.3 6.2 2 1 4.2 16.5 12 3 1 12.5 9.1 12.7 4 0 7.1 6.2 11.1 5 1 7.2 10.8 9.5 6 1 13.5 10.1 5.7 7 0 8.7 10 12.9 8 1 14.1 8.3 12.1 9 0 6.4 11.2 13.8 10 1 9.5 12.6 7.6 11 0 8.5 6.9 8.2 12 0 10.8 5.1 10.5 13 0 13.1 9.9 18.1 14 1 11 8 11.8 15 1 12.9 9.9 8.7 16 0 2.1 7 11.2 17 0 12.8 6.9 10.4 18 1 6.7 10.2 3.9 19 0 11.6 4.4 11 20 1 14 9.2 8.9 21 1 10.4 9.9 7.9 22 0 6.2 7.2 17 23 0 7.6 7.1 13.6 24 1 9.1 10.4 7.4 25 0 7.6 6.7 12.2 26 0 8.8 3.1 11.2 27 0 9 12.9 10.8 28 1 11.8 13.8 11.3 29 1 5.1 6.7 8.5 30 0 7.5 8.5 6 -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 111.80.51.185 ※ 文章網址: https://www.ptt.cc/bbs/Statistics/M.1469090479.A.85B.html ※ 編輯: o1o3o1o31030 (111.80.51.185), 07/21/2016 18:05:30

07/22 00:22, , 1F
null model fit出來beta0 = 0
07/22 00:22, 1F

07/22 00:22, , 2F
全部的y hat是0.5
07/22 00:22, 2F

07/22 00:23, , 3F
LL = sum(nrow(data) * ln(0.5)) = 41.58883
07/22 00:23, 3F

07/22 00:23, , 4F
更正數字 -20.79442
07/22 00:23, 4F

07/22 00:24, , 5F
所以-2LL = 41.58884
07/22 00:24, 5F

07/22 00:24, , 6F
0.5的原因是 全部係數為0,所以exp(XB) = 1
07/22 00:24, 6F

07/22 00:24, , 7F
P(Y=1|X) = P(Y=0|X) = 0.5
07/22 00:24, 7F

07/22 01:42, , 8F
了解! 但是這樣感覺怎麼有種怎樣都會顯著的感覺呢?
07/22 01:42, 8F

07/22 02:34, , 9F
你所謂的“顯著”是?
07/22 02:34, 9F
文章代碼(AID): #1Na8glXR (Statistics)