Задание случайных уклонов только для одной вложенной группы с использованием смешанной модели lmer

Question

У меня есть данные о урожайности за несколько лет и в разных местах.Места вложены в провинцию.Есть 5 провинций.Я хочу указать модель, в которой временная тенденция (переменная года) зависит от местоположения, но не от провинции.Также я хочу указать модель для отображения местоположений в провинции.Это правильный способ указания такой модели в лмер.

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.

Answer 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) выглядит неправильно / не нужно.
• Кажется странным предположить, что существует ровно ноль различий в пересечениях и склонах между провинциями.