james9609 发表于 2013-7-27 22:21 
你第二层单位数量是多少,如果第二层单位个数比较少的话,可以会使得迭代计算失败,这与是否加入虚拟变量应 ...
感谢 @james9609 老师!我把情况再仔细解释一下。
样本学生来自同一所大学、不同专业的大一新生。hschool 变量是大一新生所毕业的
高中是否为重点中学。
第二层是专业(major),我将这所大学的所有专业按照学科属性分成了6个大类。
研究设计一个随机系数模型,也就是在第二层没有解释变量,只有一个常数项和随机项,也可以说是“随机系数模型”:
第一层: GPA = a+ bX1+ cX2+ dX3 + f sex + g hschool + e
第二层:没有解释变量,只在第一层解释变量中加入随机效应
a= α0 + γ0 (γ0为随机项)
f = α1 + γ1 (γ1为随机项)
g = α2 + γ2 (γ2为随机项)
c= α3 + γ3 (γ3为随机项)
将第二层代入第一层,这样第一层的常数项、X2、sex 和hscool 三个解释变量就存在随机效应了。
我写的命令: xtmixed GPA x1 x2 x3 i.sex i.hschool || major : R.sex R.hshool x2 , variance mle
报错提示: “invalid level specification”
我发现stata在第二层当中如果有虚拟变量的话,就不可以再加入其它任何变量,有且只能有这一个虚拟变量了。
现在的问题是,
第二层有多个虚拟变量、或者第二层同时有虚拟变量和连续变量 的时候,stata的命令该怎么写呢?