关于空模型：组内相关系数以及模型的选择

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

各位前辈，在下想咨询一个关于多水平模型的问题，我的目的是想建立一个社区-个体-调查时间的三水平空模型，看看建立三水平模型分析是否合适。在stata上写语句   . xtmelogit hamin2||COMMID:||IDind:
得到的结果是：

1. Refining starting values:
   
2. 
   
3. Iteration 0:   log likelihood = -2930.1485  
   
4. Iteration 1:   log likelihood = -2894.2116  
   
5. Iteration 2:   log likelihood = -2868.7771  
   
6. 
   
7. Performing gradient-based optimization:
   
8. 
   
9. Iteration 0:   log likelihood = -2868.7771  
   
10. Iteration 1:   log likelihood = -2868.6202  
    
11. Iteration 2:   log likelihood = -2868.6199  
    
12. Iteration 3:   log likelihood = -2868.6199  
    
13. 
    
14. Mixed-effects logistic regression               Number of obs      =      5430
    
15. 
    
16. --------------------------------------------------------------------------
    
17.                 |   No. of       Observations per Group       Integration
    
18. Group Variable |   Groups    Minimum    Average    Maximum      Points
    
19. ----------------+---------------------------------------------------------
    
20.          COMMID |      308          1       17.6         66           7
    
21.           IDind |     3998          1        1.4          3           7
    
22. --------------------------------------------------------------------------
    
23. 
    
24.                                                 Wald chi2(0)       =         .
    
25. Log likelihood = -2868.6199                     Prob > chi2        =         .
    
26. 
    
27. ------------------------------------------------------------------------------
    
28.       hamin2 |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
    
29. -------------+----------------------------------------------------------------
    
30.        _cons |  -2.272335   .1503214   -15.12   0.000    -2.566959    -1.97771
    
31. ------------------------------------------------------------------------------
    
32. 
    
33. ------------------------------------------------------------------------------
    
34.   Random-effects Parameters  |   Estimate   Std. Err.     [95% Conf. Interval]
    
35. -----------------------------+------------------------------------------------
    
36. COMMID: Identity             |
    
37.                    sd(_cons) |    1.31052   .1211253       1.09338    1.570783
    
38. -----------------------------+------------------------------------------------
    
39. IDind: Identity              |
    
40.                    sd(_cons) |   2.002017   .1565074      1.717612    2.333513
    
41. ------------------------------------------------------------------------------
    
42. LR test vs. logistic regression:     chi2(2) =   448.62   Prob > chi2 = 0.0000
    
43. 
    
44. Note: LR test is conservative and provided only for reference.
    
45. 
    
46. .

复制代码

我在书上看这种分三层的模型在随机效应结果处会出现三行sd，然后用每一层的σ除以三个相加的sd就能得到ICC，但为什么我的结果里面只有二三水平的sd？是默认第一水平sd为1？不知怎么算ICC，如果不算icc就无法知道三水平是否合适呀。
另外我也想请教关于两水平和三水平模型的比较，做lrtest这样可以吗？显然三水平模型的likelihood增加了，这样怎么说明三水平模型优于二水平呢？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：相关系数 Optimization observations Conservative observation 模型

你好，你这个问题解决了吗？

506232839 发表于 2015-9-12 22:08
你好，你这个问题解决了吗？

请问一下，您知道这个怎么做了吗，求解答啊

因为第一层logit的方差是固定的，3.2左右

ICC实际上是高水平单位方差占总方差的比重，来呈现处数据在高水平单位的聚集性（非独立性），三水平模型中残差是水平1单位的，个体是水平2单位，社区是水平3 单位，那么就是分别用水平2单位和水平3单位方差除以总方差，求得相应层的ICC,不需要关注水平1 残差的比重

506232839 发表于 2016-3-23 08:12
因为第一层logit的方差是固定的，3.2左右

请问对于广义线性多层模型，第一层的方差都是固定的吗？为π^2/3吗