lme4::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 或留空）的地方是 Fitting Linear Mixed-Effects Models using lme4 的第 7 页。这里的 table 与 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

lme4::lmer() 用 A -1？

lme4::lmer() With A -1?

r

statistics

lme4

random-effects