作者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://webptt.com/cn.aspx?n=bbs/Statistics/M.1469090479.A.85B.html
※ 编辑: o1o3o1o31030 (111.80.51.185), 07/21/2016 18:05:30
1F：→ celestialgod: null model fit出来beta0 = 0 07/22 00:22
2F：→ celestialgod: 全部的y hat是0.5 07/22 00:22
3F：→ celestialgod: LL = sum(nrow(data) * ln(0.5)) = 41.58883 07/22 00:23
4F：→ celestialgod: 更正数字 -20.79442 07/22 00:23
5F：→ celestialgod: 所以-2LL = 41.58884 07/22 00:24
6F：→ celestialgod: 0.5的原因是 全部系数为0，所以exp(XB) = 1 07/22 00:24
7F：→ celestialgod: P(Y=1|X) = P(Y=0|X) = 0.5 07/22 00:24
8F：→ o1o3o1o31030: 了解! 但是这样感觉怎麽有种怎样都会显着的感觉呢? 07/22 01:42
9F：→ allen1985: 你所谓的“显着”是? 07/22 02:34