r - 估计没有随机效应的固定效应

Question

可能是基本问题，但我不知道该怎么做。如果一个人想估计一个多级模型，但只有固定效应，而不是估计随机效应，那该怎么做呢？例如，这是我的模型：

model4 <- lmer(m ~ x + (  | id) , data = data_stacked2, REML = TRUE)

谢谢。

Answer 1

如果我理解正确，您想使用：

lmer(m ~ x + (1|id), data = data_stacked2, REML = TRUE)

更具体地说，上述模型具有固定x效应 ( ) 和随机效应 ( id)。如果您不需要估计随机效应，请lm(m ~ x, data)改用。

但是，如果您尝试运行具有随机斜率和固定截距的模型，您可能需要尝试以下公式

lmer(m ~ 1 + (0 + x|id), data = data_stacked2, REML = TRUE)

请澄清你的问题。