混合模型的 Lme4 语法

Question

我有 $m=4$ 组老鼠（即 group1、group2、group3、group4）。每组有不同数量的小鼠。我在 $l=4$ 不同的状态（即 state1、state2、state3、state4）对每组的每只老鼠测量参数 $(y)$。我想建立一个混合效应模型来分析组、状态和组*状态的影响，允许每只小鼠和每组内的可变性。$group_{m}$ 中的小鼠标有 id (1 ,2,3...,$group_{m}$) 的老鼠数量

$$y_{mln}=\mu +group_{m} +state_{l} +(group*state){ml}+b{ml}+\varepsilon_{mln} $$ $b_{ml}$ 对 $group_{m}$

中的第 n 个鼠标的随机效应

我的数据框有以下变量

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）并估计每个组和每只鼠标的效应的随机截距。