lme4 4：：lmer()与A -1？

Question

我正在看一个已经编码的lmer模型，我不太明白-1是/正在做什么。代码看起来像fit = lmer(resids ~ -1 + (1|loc/time))

我相信(1|loc/time)部分可以等价地写成(1|loc) + (1|loc:time)，这是loc的随机截距，而time的随机截距在loc中变化。

现在我不太明白的部分是：-1，我认为这与平均值有关。我发现在那个位置使用-1 (相对于1或保留空白)的唯一地方是用lme4拟合线性混合效应模型的第7页。这里的表显示了它与offset(o)的结合，它用于“指定随机拦截具有先验已知手段”。所以，我的直觉是，去掉offset(o)和使用offset(0) (数字0而不是字符o)是一样的，这意味着先验的方法都是0。

这是正确的吗？

Answer 1

是的，这将模型的固定效应分量设置为零。虽然这是合法的，但这是一个相当奇怪的模型；我希望编写了它的人知道他们在做什么。我只能想出两个理由:你要给一个模型配上一个空的固定效果组件：

• 由于某些原因，您希望对拦截的重要性进行似然比检验(这是不寻常的，在大多数情况下，拦截并不具有特定的统计意义)。
• 你有一个特殊的实验设计，它的平均值被先验地知道为零(例如，你的响应变量是两个被随机交换的元素之间的某种差异)。

lmer(Reaction ~ -1 + (Days|Subject), sleepstudy)
Linear mixed model fit by maximum likelihood  ['lmerMod']
Formula: Reaction ~ -1 + (Days | Subject)
   Data: sleepstudy
      AIC       BIC    logLik  deviance  df.resid 
1840.7814 1853.5532 -916.3907 1832.7814       176 
Random effects:
 Groups   Name        Std.Dev. Corr
 Subject  (Intercept) 252.53       
          Days         11.93   0.88
 Residual              25.59       
Number of obs: 180, groups:  Subject, 18
No fixed effect coefficients