Stata MIXED命令提取个体截距与斜率：加入week²项报错求助

错误原因 #

当你在随机效应部分指定|| id: week week2时，Stata默认会包含个体随机截距，再加上week和week2的随机斜率，总共会生成3个随机效应项。但你的predict命令只指定了2个变量（u1 u0），变量数量不匹配，因此触发报错。

修正后的代码 #

// 拟合包含week和week2随机斜率的混合效应模型
mixed weight week week2 || id: week week2

// 提取3个随机效应项：u0=截距随机部分，u1=week随机斜率，u2=week2随机斜率
predict u0 u1 u2, reffects

// 生成个体特定的截距、week斜率、week2斜率
generate intercept = _b[_cons] + u0
generate slope_week = _b[week] + u1
generate slope_week2 = _b[week2] + u2

补充说明 #

如果不需要week2的随机斜率，仅保留截距和week的随机斜率，可以修改模型为：

mixed weight week week2 || id: week

这种情况下，week2仅作为固定效应项存在，随机效应只有截距和week的斜率，此时你原来的predict和生成变量的代码就可以正常运行。

内容的提问来源于stack exchange，提问作者bvowe