Синтаксис Lme4 для смешанной модели

Question

У меня $ m = 4 $ группа мышей (например, group1, group2, group3, group4). У каждой группы разное количество мышей. Я измеряю параметр $ (y) $ на каждой мыши каждой группы при $ l = 4 $ разных состояниях (т.е. состояние 1, состояние 2, состояние 3, состояние 4). Я хотел бы создать модель смешанного эффекта для анализа влияния состояния группы, состояния и группы *, учитывающего изменчивость внутри каждой мыши и внутри каждой группы. Мыши в пределах $ group_ {m} $ помечены знаком id (1 , 2,3 ..., число мышей $ group_ {m} $)

$$ y_ {млн.} = \ Mu + group_ {m} + state_ {l} + (состояние группы *) { мл} + B {мл} + \ varepsilon_ {млн} $$ с $ B_ {мл} $ случайный эффект для п-го мыши в пределах $ group_ {M} $

Мой кадр данных имеет следующие переменные

value (num) 
state (factor: 4 levels) 
group (factor: 4 levels) 
id (within group) (num) 

Правильный ли соответствующий синтаксис?

lmer(value~group+state+group*state+(1|group)) 

Answer 1

Вы хотите, чтобы это

mouseID <- interaction(group, ID) 
lmer(value ~ group * state + (1|mouseID)) 

mouseID должен быть уникальным для каждой мыши.

Поскольку группа является фактором, вы не можете иметь ее как в фиксированной, так и в случайной части. Это привело бы к неидентифицируемой модели.

Answer 2

Я думаю, что вы ищете,

lmer(value ~ group*state + (1|group) + (1|id)) 

Эта модель оценивает фиксированный эффект группы и государства, а также взаимодействие между ними (R автоматически расширяется group*state к group + state + group*state) и оценки случайный перехват для эффекта каждой группы и для каждой мыши.