首页
学习
活动
专区
工具
TVP
发布
精选内容/技术社群/优惠产品,尽在小程序
立即前往

什么时候应该在lmer()中指定REML = FALSE?

在使用lmer()函数进行线性混合模型拟合时,可以通过指定REML参数来选择估计固定效应的方法。REML代表了“Restricted Maximum Likelihood”,是一种用于估计模型参数的方法。

当REML参数设置为FALSE时,表示使用最大似然法(Maximum Likelihood,ML)来估计固定效应。ML方法在模型拟合时不考虑随机效应的方差,因此可能会导致估计值的偏差。

在以下情况下,可以考虑将REML参数设置为FALSE:

  1. 当样本量较大时,ML方法在计算过程中可能更加耗时。而如果样本量较大,那么使用REML方法可能更加高效。
  2. 当模型存在奇异性(singularity)或过度拟合(overfitting)的情况时,ML方法可能更适用。

需要注意的是,使用ML方法估计固定效应可能会导致随机效应的方差被低估,因此对于随机效应的解释可能不准确。在这种情况下,可以使用REML方法来更准确地估计随机效应的方差。

总结起来,当样本量较大、模型存在奇异性或过度拟合时,可以考虑在lmer()中指定REML = FALSE来使用最大似然法估计固定效应。

页面内容是否对你有帮助?
有帮助
没帮助

相关·内容

R语言用线性混合效应(多水平层次嵌套)模型分析声调高低与礼貌态度的关系|附代码数据

REML与ML 让我们从一个统计模型开始,指定(i)固定效应和(ii)各种随机效应的正态分布的变异和协方差。...在比较固定效应不同的模型时,你应该使用ML,而且你必须包括lmer(, REML=FALSE)。...在REML(REstricted ML)估计中,我们的主要兴趣是估计随机效应,而不是固定效应。想象一下,我们通过将上面第(i)组中的固定效应参数设置为某些合理的值来限制我们的参数空间。...这种平均法可以得到REML参数估计值和REML偏差值。因为这个过程对固定效应参数的关注度很低,所以它不应该被用来比较固定效应结构不同的模型。你应该在比较随机效应不同的模型时使用这个方法。...要做到这一点,你应该养成在运行模型比较时包括lmer(, REML=TRUE),并使用anova(, refit=FALSE)的习惯。

64000

R语言用线性混合效应(多水平层次嵌套)模型分析声调高低与礼貌态度的关系|附代码数据

REML与ML 让我们从一个统计模型开始,指定(i)固定效应和(ii)各种随机效应的正态分布的变异和协方差。...在比较固定效应不同的模型时,你应该使用ML,而且你必须包括lmer(, REML=FALSE)。...在REML(REstricted ML)估计中,我们的主要兴趣是估计随机效应,而不是固定效应。想象一下,我们通过将上面第(i)组中的固定效应参数设置为某些合理的值来限制我们的参数空间。...这种平均法可以得到REML参数估计值和REML偏差值。因为这个过程对固定效应参数的关注度很低,所以它不应该被用来比较固定效应结构不同的模型。你应该在比较随机效应不同的模型时使用这个方法。...要做到这一点,你应该养成在运行模型比较时包括lmer(, REML=TRUE),并使用anova(, refit=FALSE)的习惯。

64400
  • R语言LME4混合效应模型研究教师的受欢迎程度|附代码数据

    R中编码的基础知识。安装R软件包  lme4,和  lmerTest。...现在我们可以为数据中的100个不同类别绘制不同的回归线我们清楚地看到,外向性和受欢迎程度之间的关系在所有层级中并不相同,但平均而言,存在明显的正向关系。...                  alpha = as.factor(high_and_low)),              method = lm,              se     = FALSE...最后,我们在data = 命令后指定要使用的数据集summary(interceptonlymodel) #得到参数估计.## 通过REML进行线性混合模型拟合。...要在LMER中完成此操作,只需将随机斜率的变量添加到输入的随机部分。  (1|class)变成 (1+sex+extrav |class)。

    1K10

    BLUP is a good thing!

    BLUP值在育种中称为育种值。混合线性模型在基因组选择中,是GBLUP,ssGBLUP。 这个模型,是农业数据分析汇总应用最广泛的模型,值的好好学习! 至于LMM应该怎么学习?...数据描述 ❝睡眠剥夺研究中受试者每天的平均反应时间。第0天,受试者有正常的睡眠时间。从那天晚上开始,他们每晚只能睡3个小时。观察结果代表了每天对每个受试者进行的一系列测试的平均反应时间。...lme4的基本语法: library(lme4) mod1a = lmer(Reaction ~ Days + (1 | Subject), data=dat) summary(mod1a) 结果...: > summary(mod1a) Linear mixed model fit by REML ['lmerMod'] Formula: Reaction ~ Days + (1 | Subject...固定因子中,Days为数值协变量,截距的值为251.4,Days的值为10.46 「asreml:」 代码: library(asreml) mod1b = asreml(Reaction ~ Days

    58920

    R语言LME4混合效应模型研究教师的受欢迎程度|附代码数据

    p=11724文中本教程对多层_回归_模型进行了基本介绍(点击文末“阅读原文”获取完整代码数据) 。 介绍本教程期望: 多层_回归_模型的基础知识 。 R中编码的基础知识。...alpha = as.factor(high\_and\_low)), method = lm, se = FALSE...如果我们查看LMER函数的不同输入,则: “受欢迎程度”,表示我们要预测的因变量。 一个“〜”,用于表示我们现在给出了其他感兴趣的变量。(与回归方程式的'='相比)。...最后,我们在data = 命令后指定要使用的数据集 summary(interceptonlymodel) #得到参数估计.## 通过REML进行线性混合模型拟合。...要在LMER中完成此操作,只需将随机斜率的变量添加到输入的随机部分。 (1|class)变成 (1+sex+extrav |class)。

    79230

    R语言LME4混合效应模型研究教师的受欢迎程度

    现在我们可以为数据中的100个不同类别绘制不同的回归线 ? 我们清楚地看到,外向性和受欢迎程度之间的关系在所有阶层中并不相同,但平均而言,存在明显的正向关系。...group = class, size = as.factor(high_and_low), alpha = as.factor(high_and_low)), method = lm, se = FALSE...如果我们查看LMER函数的不同输入,则: “受欢迎程度”,表示我们要预测的因变量。 一个“〜”,用于表示我们现在给出了其他感兴趣的变量。(与回归方程式的'='相比)。 公式中表示截距的“ 1”。...最后,我们在data = 命令后指定要使用的数据集 summary(interceptonlymodel) #得到参数估计. ## 通过REML进行线性混合模型拟合。...要在LMER中完成此操作,只需将随机斜率的变量添加到输入的随机部分。 (1|class)变成 (1+sex+extrav |class)。

    1K20

    多水平模型、分层线性模型HLM、混合效应模型研究教师的受欢迎程度

    现在我们可以为数据中的100个不同类别绘制不同的回归线 我们清楚地看到,外向性和受欢迎程度之间的关系在所有层级中并不相同,但平均而言,存在明显的正向关系。...alpha = as.factor(high\_and\_low)), method = lm, se = FALSE...如果我们查看LMER函数的不同输入,则: “受欢迎程度”,表示我们要预测的因变量。 一个“〜”,用于表示我们现在给出了其他感兴趣的变量。(与回归方程式的'='相比)。 公式中表示截距的“ 1”。...最后,我们在data = 命令后指定要使用的数据集 summary(interceptonlymodel) #得到参数估计. ## 通过REML进行线性混合模型拟合。...要在LMER中完成此操作,只需将随机斜率的变量添加到输入的随机部分。 (1|class)变成 (1+sex+extrav |class)。

    1.5K20

    R语言实现混合模型

    噪声是我们模型中没有考虑的随机因素。而固定效应是那些可预测因素,而且能完整的划分总体。例如模型中的性别变量,我们清楚只有两种性别,而且理解这种变量的变化对结果的影响。...在R语言中我们使用mgcv包中的lmer函数来完成这项工作。首先载入faraway包以便读取psid数据集,然后加载mgcv包,再将年份数据中心化以方便解释模型,最后用lmer函数进行建模。...dat <- data.frame(rbind(dat.tx.a, dat.tx.b)) names(dat) <- c('measure.1', 'measure.2', 'measure.3') 根据指定的均值和协方差生成多元正态数据...结果显示如下: Linear mixed model fit by REML ['lmerMod'] Formula: measure ~ time + (1 | subject.id) Data...: dat REML criterion at convergence: 9750.6 Scaled residuals: Min 1Q Median 3Q

    4.4K70

    R语言广义线性混合模型GLMMs在生态学中应用可视化2实例合集|附数据代码

    另外,关于嵌套和交叉随机效应的问题,lme4包中的lmer函数支持拟合这些复杂的随机效应结构。你可以通过在公式中指定适当的随机效应项来实现这一点。...以下是这段代码的详细解释: # 实践2测试随机效应 # 第一个模型 # 第二个模型没有随机项,使用gls是因为它也可以通过REML拟合模型 # 似然比检验,对于小样本量来说可能不够精确...两种模型都考虑了随机效应,但使用了不同的估计方法(lme使用最大似然法,lmer使用REML或ML)。 参数自助法似然比检验:与之前的自助法类似,但这次是针对固定效应部分的模型进行比较。...在此步骤中,我们将计算将其增加到 20 年的影响。 modl2 <‐ extend power(modl2) 沿参数指定要扩展的变量,n 指定要替换它的级别。...该分析表明,该研究必须运行 16 年才能有≥80% 的功效来检测指定大小的影响。

    17410

    R语言广义线性混合模型GLMMs在生态学中应用可视化2实例合集|附数据代码

    另外,关于嵌套和交叉随机效应的问题,lme4包中的lmer函数支持拟合这些复杂的随机效应结构。你可以通过在公式中指定适当的随机效应项来实现这一点。...以下是这段代码的详细解释: # 实践2测试随机效应 # 第一个模型 # 第二个模型没有随机项,使用gls是因为它也可以通过REML拟合模型 # 似然比检验,对于小样本量来说可能不够精确...两种模型都考虑了随机效应,但使用了不同的估计方法(lme使用最大似然法,lmer使用REML或ML)。 参数自助法似然比检验:与之前的自助法类似,但这次是针对固定效应部分的模型进行比较。...在此步骤中,我们将计算将其增加到 20 年的影响。 modl2 <‐ extend power(modl2) 沿参数指定要扩展的变量,n 指定要替换它的级别。...该分析表明,该研究必须运行 16 年才能有≥80% 的功效来检测指定大小的影响。

    1.4K10

    BOLT-LMM用户手册笔记

    更准确地说: M = bim 文件中满足所有条件的 SNP 的数量: 未在任何 --exclude 文件中列出 通过 QC 过滤器筛选缺失 在 --modelSnps 文件中列出(如果已指定) N =...5.6 用户指定的筛选 要从分析中删除的单个值可以在一个或多个 --remove 文件中指定,其中列出了 FID 和 IID(每行一个单独的)。...请注意,即使指定了 --modelSnps 文件,基因型数据中的所有 SNP 仍会进行关联测试;只有混合模型中的随机效应被限制为 --modelSnps。...7.3 初始方差参数猜测 要指定一组用于启动 REML 迭代的方差参数(与 BOLT-REML 使用的默认过程相比,如果您有良好的初始猜测,这可能会节省时间),请使用 --remlGuessStr="string...将这些值视为上三角形矩阵的条目(对角线上有方差比例,对角线上有相关性),您应该在每个方差分量名称后指定这些 D(D+1)/2 值,方法是从左到右、从上到下读取它们。

    2.7K41

    GCTA如何计算GWAS中SNP的遗传力

    GWAS遗传力是显著SNP的解释百分比,具体可以参考我写的系列博客:GWAS分析中SNP解释百分比PVE | 第四篇,MLM模型中如何手动计算PVE?...1.2 --reml-alg 指定迭代方法(非必须) --reml-alg 0 # AI算法,默认算法 --reml-alg 1 # EM算法 1.3 --reml-priors 迭代初始值 指定迭代初始值...,当数据量较大时,指定初始值,会加快迭代速度。...--reml-priors 0.45 0.55 表示Va为0.45, Ve为0.55 1.4 --grm(必须) 指定GRM矩阵 --grm # 接二进制文件GRM的前缀 --grm-bin # 同上...2.2 基因型数据 plink的二进制文件 2.3 协变量 这里,示例数据中,没有提供协变量信息。如果提供,可以按照第一列是FID,第二列是IID,其它是协变量的方法整理数据。

    1.4K10
    领券