问题描述
假设我想用 lme4 计算随机模型效应。 这是一个示例模型以及通过明确指定从一个变量中获取随机模型效应的代码:
library(lme4)
model<-lmer(mpg ~ 1+ (1|wt) + (1|hp),data = mtcars,control=lmerControl(optCtrl=list(maxfun=20000)))
ranef(model)$wt
但是,我想在一个更灵活的函数中构建它,为此我需要粘贴我想要来自这样的向量的随机效果的效果:
randomeffectlist<-NULL
effectlist<-c("wt","hp")
for (y in length (effectlist)){
randomeffectlist[y]<-ranef(model)$effectlist[y]
}
在 $sign 之后,我想一个接一个地传递 factorlist 的元素,以便为每个效果获取 randomeffectlist 中的随机效果。在当前形式中,它不起作用。我想这是在 $ 符号之后传递一些东西的问题。如何更正此代码?
解决方法
计算要提取的元素时不能使用 $
:
ranef(model)$effectlist[y]
将逐字查找名为 "effectlist[y]"
的列表元素。相反,您应该使用双方括号运算符:
ranef(model)[[effectlist[y]]]
。
来自 [
的帮助页面:
"[[
和 $
都选择列表中的单个元素。主要区别在于 $
不允许计算索引,而 [[
允许。{{ 1}} 等价于 x$name
。另外,x[["name",exact = FALSE]]
的部分匹配行为可以使用精确参数来控制。"