在不估计整体平均斜率的情况下估计随机斜率是否合适?

机器算法验证 混合模式 lme4-nlme 斯坦
2022-03-14 12:29:27

我试图估计不同城市的个人(我的分组变量)对一些预测变量的反应是否存在差异。所以,在实践中,我有兴趣了解β来自每个城市。但是,我想使用随机斜率,以便这些估计值“缩小”到某个组均值。

我正在使用stan_glmerfrom包来拟合我的模型,它以与fromstanarm相同的方式接受输入glmerlme4

所以,问题#1。如果我不关心大坡度,可以这样指定模型:

stan_glmer(Y ~ (V1 + V2 + V3|city))

还是我需要先将 V1、V2、V3 包含为“固定效果”?例如

stan_glmer(Y ~ V1 + V2 + V3 + (V1 + V2 + V3|city))

其次,如果第二个模型是正确的选项(我有点怀疑),那么我如何解释输出?为了估计城市 1 中 V1 的斜率,我是否添加了整体βv1到城市 1βv1? 在这种情况下,您如何解释错误?

谢谢!

1个回答

将人口水平斜率固定为零的随机斜率拟合并不是不可能的——它在数学或统计上不是不合适的——但这是一个相当奇怪的模型,需要一些额外的理由。为什么你会期望跨城市的平均斜率恰好为零(这是忽略固定效应的模型所暗示的)?我看到拟合此类模型的唯一情况是

  • 作为(n)(诚然愚蠢的)空模型,用于对总体水平斜率的显着性进行似然比检验[在您使用贝叶斯方法时与您的情况无关]
  • 在根据实验设计效果为零的情况下,例如,当样本被随机分配到预处理条件下的测试和处理条件时(这将消除“之前”期间处理的固定效果,而不是固定的斜率,但想法是相似的)。

如果您有固定效应斜率和斜率的城市间变化,您确实需要将人口水平斜率添加到单个城市斜率偏差,并使用总和的后验分布进行推断 - 我不不确切知道这是如何在rstanarm相关的)中完成的……该tidybayes软件包可能会有所帮助。