作者o1o3o1o31030 (主席)
看板Statistics
標題[問題] 關於羅吉斯-2LL
時間Thu Jul 21 16:41:16 2016
各位版大好
關於羅吉斯分析中的-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
→ celestialgod: null model fit出來beta0 = 0 07/22 00:22
→ celestialgod: 全部的y hat是0.5 07/22 00:22
→ celestialgod: LL = sum(nrow(data) * ln(0.5)) = 41.58883 07/22 00:23
→ celestialgod: 更正數字 -20.79442 07/22 00:23
→ celestialgod: 所以-2LL = 41.58884 07/22 00:24
→ celestialgod: 0.5的原因是 全部係數為0，所以exp(XB) = 1 07/22 00:24
→ celestialgod: P(Y=1|X) = P(Y=0|X) = 0.5 07/22 00:24
→ o1o3o1o31030: 了解! 但是這樣感覺怎麼有種怎樣都會顯著的感覺呢? 07/22 01:42
→ allen1985: 你所謂的“顯著”是? 07/22 02:34