我试图估计不同城市的个人(我的分组变量)对一些预测变量的反应是否存在差异。所以,在实践中,我有兴趣了解来自每个城市。但是,我想使用随机斜率,以便这些估计值“缩小”到某个组均值。
我正在使用stan_glmer
from包来拟合我的模型,它以与fromstanarm
相同的方式接受输入。glmer
lme4
所以,问题#1。如果我不关心大坡度,可以这样指定模型:
stan_glmer(Y ~ (V1 + V2 + V3|city))
还是我需要先将 V1、V2、V3 包含为“固定效果”?例如
stan_glmer(Y ~ V1 + V2 + V3 + (V1 + V2 + V3|city))
其次,如果第二个模型是正确的选项(我有点怀疑),那么我如何解释输出?为了估计城市 1 中 V1 的斜率,我是否添加了整体到城市 1? 在这种情况下,您如何解释错误?
谢谢!