大神求助！
毕业课题的一个问题，就是现在我利用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
}