使用 lmer - 混合模型仅为一个嵌套组指定随机斜率答案

【问题标题】：Specifying random slopes for only one nested group using lmer - mixed model使用 lmer - 混合模型仅为一个嵌套组指定随机斜率
【发布时间】：2019-11-15 02:06:45
【问题描述】：

我有多年和地点的作物产量数据。位置嵌套在一个省内。有5个省。我想指定一个模型，其中时间趋势（年份变量）因地点而异，但不因省而异。另外我想指定模型以显示位置嵌套在省内。这是在 lmer 中指定此类模型的正确方法吗？

lmer(y ~ x1 + x2 + x3 + (year|location) + (1|province/location), data = dat)

位置和省份编码为

province 1 -  location as 11, 12, 13, 14 etc
province 2 - location as 21, 22, 23, 24 etc 
province 3 - location as 31, 32, 33, 34 etc.

【问题讨论】：

您的state 变量是什么？这根本不在你的口头描述中......它应该是location吗？
抱歉，它应该是一个位置。我提供了有关我的位置和省份如何编码的更多信息

标签： r lme4 mixed-models

【解决方案1】：

如果您的位置是唯一编码，即在 A、B、C、D、E 省，您有 A1、A2、A3、...、B1、B2、B3 等州。 ., C1, C2, C3, ...，然后(year|location) 足以指定跨位置的截距和斜率的随机变化（根据对您问题的编辑，情况确实如此）。

如果不是（即您的位置标记为 1、2、3、...、1、2、3、... 1、2、3、...），那么您需要(year|province:location)，其中: 表示交互。

(year|province/location) 表示省内和省内位置之间的随机截距和斜率。
在任何一种情况下，(1|province/location) 看起来都是错误的/不必要的。
假设各省之间的截距和斜率恰好零变化似乎很奇怪。

【讨论】：