多层线性模型案例(stata代码和数据)数据格式dta格式
代码文件do文件,代码附有详细解释和说明
代码示例如下:
use 多层线性模型.dta, clear
*reml要求采用约束最大似然法估计参数,也可直接使用默认的极大似然估计mle,本例为了保持与HLM软件结果一致.
*nolog略去模型拟合的迭代信息
mixed mathach || school:, reml nolog
*检验截距是否存在显著校际差异y=β+rij
*方式1: Wald检验(Hox,2010)
disp "z = " 8.614081/1.078813
disp "P > |z| = " (1-normal(8.614081/1.078813))
*方式2: LR CHi2检验
quietly mixed mathach || school:, reml
estimates store m0
*简单模型,school不存在随机效应
quietly mixed mathach, reml
lrtest m0
*显著,拒绝简单模型,使用复杂的m0模型
*计算组内相关系数
qui mixed mathach || school:, reml nolog
estat icc
*值等于8.614081/(8.614081+39.14832)