大神求助!
毕业课题的一个问题,就是现在我利用stata循环产生了比如10个的多水平数据(单水平数据也是一个道理),然后依次对每个模型进行混合效应分析,那个如何把这10个混合模型的系数,比如x前的系数,置信区间提取出来?因为还想基于这10个模型的系数及一些性质进行分析!真实模拟是可能生存1000个多水平数据,所以还是得靠代码才可以!真心感谢!
谢谢,附带我的代码:
clear
set seed 1234
local i=1
preserve
while `i' <= 10{
set obs 6
generate school = _n
generate u_i = rnormal(0,2)
expand 10
bysort school: generate classroom = _n
generate u_ij = rnormal(0,3)
expand 16+int((25-16+1)*runiform())
bysort school classroom: generate child = _n
generate e_ijk = rnormal(0,5)
generate x_ijk=rnormal(0,1)
generate y = 70 + u_i + u_ij +5*x_ijk+ e_ijk
mixed y x, stddev
est store model`i'
local i = `i'+1
restore,preserve
}