如何对lmer()中的随机部分使用update()？
在R语言中，lmer()函数是用于拟合线性混合效应模型的函数，其中包含固定效应和随机效应。update()函数可以用于更新模型对象，包括lmer()函数返回的模型对象。
对lmer()中的随机部分使用update()的方法是通过指定新的随机效应来更新模型。具体步骤如下：
- 首先，使用lmer()函数拟合一个初始的线性混合效应模型，并将其存储在一个对象中，例如model。
- 然后，使用update()函数来更新模型对象。在update()函数中，通过指定formula参数来指定新的随机效应。新的随机效应应该使用“~”符号来表示，并与原始模型中的随机效应进行组合。例如，如果原始模型中的随机效应为(1|group)，新的随机效应可以为(1|group) + (1|subject)。
- 最后，将更新后的模型对象存储在一个新的对象中，例如updated_model。
需要注意的是，update()函数只能用于更新模型的随机部分，不能用于更新模型的固定部分。
下面是一个示例代码：
# 拟合初始的线性混合效应模型
model <- lmer(y ~ x + (1|group), data = data)
# 更新模型的随机部分
updated_model <- update(model, . ~ . + (1|subject))在这个示例中，初始模型中的随机效应为(1|group)，更新后的模型中的随机效应为(1|group) + (1|subject)。
对于lmer()函数中的其他参数和用法，可以参考R语言的官方文档或相关教程。
腾讯云相关产品和产品介绍链接地址：
- 腾讯云计算服务：https://cloud.tencent.com/product
- 腾讯云数据库：https://cloud.tencent.com/product/cdb
- 腾讯云服务器：https://cloud.tencent.com/product/cvm
- 腾讯云人工智能：https://cloud.tencent.com/product/ai
- 腾讯云物联网：https://cloud.tencent.com/product/iotexplorer
- 腾讯云移动开发：https://cloud.tencent.com/product/mobdev
- 腾讯云存储：https://cloud.tencent.com/product/cos
- 腾讯云区块链：https://cloud.tencent.com/product/baas
- 腾讯云元宇宙：https://cloud.tencent.com/product/vr
相关·内容
我想使用update()函数来更新我的模型的随机部分，特别是添加一个随机效果。大多数示例(help("update")、help("update.formula")、)都集中在模型的固定部分。在下面的示例中，我如何使用update()从fm0转到fm1？library(lme4)
(fm0
浏览 0提问于2017-07-19得票数 7
回答已采纳
library(lme4)
假设我们也对从fit中提取原始公式(即固定效应和随机效应部分)感兴趣，我们想把它传递给第二个称为fit2的LMM。显然，将terms(fit)传递给lmer()的formula参数不起作用。> fit2 <- lmer(terms(fit), data = sleepstudy)
浏览 7提问于2022-01-11得票数 6
回答已采纳
我试图建立一个随机混合效应模型，以CCT作为结果，以时间*行为作为主要的预测交互作用项。我很难运行这个模型，因为这是我的代码，它一直在说“未使用的参数”错误消息。我使用变量ID作为随机效应变量，这是我的代码。谢谢!model.set$f <- with(model.set, ID)
m <- lmer(CCT ~ time*behavior, random = ~ 1|f, data = model.set)
浏览 2提问于2022-04-18得票数 0
1回答
有没有办法修改(覆盖)lmer模型中的随机效果？对于固定效果，有一个名为my_lmer@beta的插槽，我可以使用以下命令更改固定效果：对于随机效果，有没有类似的方法？lmer-object是否已经包含随机效果，或者是稍后由ranef()计算的。
浏览 5提问于2016-08-30得票数 0
回答已采纳
以M为固定效应，年份嵌套为随机效应，用线性混合模型考察M对Com的影响。如果我在包装'nlme‘中使用lme没有问题。fit <- lme(Com ~ M + random = ~ 1|plot/year, data)
# Error: number of
浏览 1提问于2019-12-12得票数 0
我从40只不同的老鼠身上收集了一次连续的测量结果。他们被分成八个围栏，在四种不同的光处理中的一种下(每种处理一式两份)。我的实验测试了两种老鼠，雄性和雌性，但这些比例并不相等。我主要感兴趣的是光处理之间的差异，而不是外壳之间的差异，尽管我最终也想评估这一点。我的数据最适合嵌套ANOVA模型还是线性混合效应模型？公式应该是什么样子的？我知道这不可能是双向方差分析，因为子组的大小不相等
浏览 18提问于2019-04-17得票数 0
回答已采纳
我想要符合以下模式：但这并不是趋同的。我尝试过使用Nelder方法和L B来修改控制。上面的简单模型产生了：checkConv中的警告(attr(opt，“派生”)，opt$par，ctrl = control$ch
浏览 3提问于2022-06-14得票数 0
我目前正在测试是否应该在我的lmer模型中包含某些随机效应。我用了anova函数。到目前为止，我的过程是用函数调用lmer()和REML=TRUE (默认选项)来拟合模型。然后，我在两个模型上调用anova()，其中一个模型包含要测试的随机效应，而另一个模型没有。然而，众所周知，anova()函数用ML来修改模型，但是在新版本的anova()中，您可以通过设置选项refit=FALSE来阻止anova()这样做。da
浏览 5提问于2014-04-06得票数 13
回答已采纳
我正在试图弄清楚如何从一个简单的lmer模型中组织随机效果的ranef.mer列表，该模型只有一个随机截获和一个变量(性别)。fit.b <- lmer(Math ~ 1 + Sex + (1+Sex|SchoolID), data=pisa_com, REML=FALSE)
我已经使用qqmath绘制了随机效果图，但是我需要能够通过集群编号(在本例中是学校)来标记每个随机效
浏览 1提问于2013-11-07得票数 1
1回答
我要模拟来自多层Logistic分布的数据，每次1000次，每次估计参数，并计算估计值的平均值。但是，在glmer()函数中，结果似乎不能像lm()函数那样提取，比如lm(y~x)$coef。如何从glmer()函数中提取结果？replicate(10,simu(30,5,-1,.3,.3,.3,.13,1,0))如何从
浏览 10提问于2015-07-09得票数 1
回答已采纳
1回答
我正在看一个已经编码的lmer模型，我不太明白-1是/正在做什么。代码看起来像fit = lmer(resids ~ -1 + (1|loc/time))
我相信(1|loc/time)部分可以等价地写成(1|loc) + (1|loc:time)，这是loc的随机截距，而time的随机截距在loc中变化。现在我不太明白的部分是：-1，我认为这与平均值有关。我发现在那个
浏览 2提问于2021-06-21得票数 1
回答已采纳
3回答
我有一个公式列表，我想在一个循环中使用一个函数。这些公式有的是随机效应模型，有的是直截了当的线性模型。我希望该函数检测模型是否包含随机效应，如果包含随机效应，则使用lmer()来拟合模型。否则，它应该使用lm()。对于如何检查此条件有何建议(除了将公式转换为字符串和检查括号外)？在这个阶段，他们有相同的课程，所以我不能只是检查。我也可以使用错误处理来捕捉当lmer()从模型
浏览 5提问于2020-12-08得票数 3
回答已采纳
1回答
我试图在R中创建一个包装函数，它将向使用lmer()的混合效应线性回归模型提供一系列参数(响应、预测器、最大随机效应、参数)。我遇到的问题是如何使包装器函数具有足够的灵活性，以适应不同数量的参数。我不知道如何做到这一点，没有硬编码参数的函数(x1，x2，x3等)。当我们考虑最大随机效应，包括随机斜率和拦截(slope|intercept)的<em
浏览 3提问于2017-03-09得票数 1
回答已采纳
1回答
我有一个70k+行的大型数据集和几个变量为data.Also的列，我有一个列，需要使用的因素超过5000个。任何帮助都将不胜感激。
浏览 4提问于2020-04-29得票数 3
回答已采纳
我试着用一种前后梯度的方法来衡量工业发展对生物的影响。我在R中使用线性混合模型方法，并且很难指定合适的模型，特别是随机效应。响应变量(例如，关键物种丰度的变化)将作为距离撞击边缘的距离的函数来衡量，使用在固定距离沿多条横断面(“梯度”)从扰动边缘辐射出来的地块来衡量。这里的效果是线性的，而不是曲线的，但是你明白了。周期和距离之间的相互作用(相当于
浏览 0提问于2013-09-07得票数 2
回答已采纳
3回答
我需要从standard error的输出中提取lmer的方差分量。library(lme4)以下是差异部分的估计数：
s2 <- VarCorr(model在第12章，随机因素的实验，由道格拉斯C.蒙哥马利(DouglasC.Montgomery)著的“”一书中，在本章末尾，示例12-2是由SAS完成<em
浏览 2提问于2015-07-29得票数 3
回答已采纳
1回答
我使用R中的lmer()函数(lme4包)来分析一项纵向研究，在其中我测量了120名受试者，6次。首先，我指定了这样一个模型：model1 = lmer(DV ~ 1 + X1*X2 + (1+X1|SubjectID), REML="false")
X1是一个时变变量因为这些主题嵌套在几个团队中，所以我被建议在团队级别(级别-3)中包含一个随机拦截。然而，我只发现如何同时包含<em
浏览 3提问于2016-11-12得票数 3
回答已采纳
2回答
0.16666667每个主题为errorType (O或P)提供两个数据点，并且每个主题处于条件G (N=30)或N (N=33)中。我对主效应和交互都很感兴趣。然而，当我使用lmer时，我得到了一组完全不同的结果：
> lmer(errors ~ Condition * errorType + (1 | subject),此外，lmer结果与glm结果完全相同，这让我相信有些地方
浏览 4提问于2014-01-07得票数 2
我不明白这两个模型之间的区别(我使用的是lmer4包) Model.1<-lmer(Y~ X1+ X2+ X3 data=Data, REML=FALSE)
+(1|Subject) +(X1|Subject)+ (X2|Subject),
d
浏览 29提问于2020-10-16得票数 0
1回答
在开始使用固定项之前，我正在尝试选择GLMM中最好的随机效果结构。为了做到这一点，我包括了所有的固定效应及其相互作用(超越最优模型)，然后我尝试使用不同的随机因素组合。我正在使用公式lmer()。用REML对模型进行估计。然后，我获得每个模型的AIC()，并对它们进行比较。
但我也想知道模型在没有随机效应的情况下的AIC。我读过了，然后我应该<e
浏览 15提问于2015-01-04得票数 1
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
云直播
对象存储活动推荐
腾讯云开发者
