如何决定将 REML 设置为 True 还是 False?

机器算法验证 r lme4-nlme
2022-03-18 16:08:00

我找到了一个网页告诉 lmer:

如果您的随机效应是嵌套的,或者您只有一个随机效应,并且您的数据是平衡的(即每个因子组中的样本量相似),请将 REML 设置为 FALSE,因为您可以使用最大似然。如果您的随机效应被交叉,请不要设置 REML 参数,因为无论如何它默认为 TRUE。

我的lmer模型中有 2 个随机效应。一种是嵌套的:

(1|Random1A/Random1B) + (1|Random2)

我应该设置REMLTRUE(默认情况下)还是FALSE

1个回答

在我(并非完全不知情)看来,您从网页和您收到的评论中得到了一些有问题的建议。

  • 您可以随时使用 REML(或 ML)(无论随机效应结构如何 - 单个与多个、平衡与不平衡、交叉与嵌套)
  • 在简单的情况下(平衡/嵌套/等),可以证明 REML 提供方差分量的无偏估计(但不能提供例如标准偏差或对数标准偏差的无偏估计)
  • 如果使用 REML 而不是 ML 拟合,则无法比较固定效应不同的模型;这就是为什么评论者建议REML=FALSE您在尝试进行模型选择时使用
  • 但是,我不建议您首先进行模型选择,当然如果您要依赖条件置信区间和 p 值(即,在不考虑影响的情况下分析改装的“最小充分”模型)型号选择)

从我在Fox et al 2015的章节中

当您对随机效应方差的大小感兴趣时,通常最好使用 REML(如果可用),但当您通过假设检验或信息论标准(如 AIC)比较具有不同固定效应的模型时,则永远不要使用它。