问题标签 [lmer]

我有一个mer用调用 to 创建的对象lmer()。

我可以得到随机效应，ranef()但我也想为每个随机效应获得相应数量的观察 - 有没有简单的方法可以做到这一点？

附加信息：

我可能没有在上面说得很清楚。例如，如果我有一个简单的 2 级模型，其中患者聚集在医院内并随机截取医院，我想提取每家医院的随机效应以及每家医院ranef() 内的患者数量。目前，我使用

这给了我类似的东西：

我想得到类似的东西：

为此，我一直在使用

然后cbind将结果从ranef(fullmodel)[[1]]

然而，这似乎不复杂且容易出错。

我想对我的数据进行分段增长回归，基于 Raudenbush/Bryk(1992/2002)，分层线性模型，p178 - 179。

我将不得不应用一个多层次模型来做到这一点。

我是我需要使用的 lme4-package。

我的数据包含每个 ID 一段时间内的销售收入。我有两个断点：sales>10 和 sales>25。我想分析断点前后的销售额增长情况。

每个ID都有两个变量：var1：123455555555（当销售额超过10（断点1）时，变量保持不变）var2：000000006789（当销售额超过25（断点2）时，变量在增加，在25之前是0)

我的多级模型具有以下形式：

Y= b0i + b1i*a1ti + b2i*a2ti + eti

（对不起，我不允许发布图片和公式：（）

a1ti 和 a2ti 是我指定片段的编码变量（var1 和 var2）。

第 1 级：每个人的重复观察
第 2 级：受试者之间的水平

也许有人可以解释一下我必须在我的 lmer()-Code 中指定级别 1 和级别 2 的位置以及我将如何获得我的两个部分的两个斜率？我是否必须在括号之间添加规格为 O + ... 或 1 + ... 才能仅获得截距等？

我想我的代码看起来像这样：

或这个：

谢谢！

最好的祝福亚历克萨

大家好，

当“lme4”中的“lmer”函数运行时，它会生成一个带有很多插槽的 S4 对象。我对这些插槽之一感兴趣，即模型@X，以及如何产生这个“X”插槽输出。我想尝试为我正在使用的不同模型函数（glmmPQL）重现这个输出，它不会自动产生这个'X'输出（仅供参考'lmer'产生一个'mer'类的对象，插槽'X'是一个模型固定效应的矩阵）。下面的代码显示了我在说什么。

我想弄清楚这个'X'数据是如何产生的？我通过在没有'()'的终端中编写'lmer'的代码来查看它，但我在那里找不到任何东西。我也试过 showMethod('lmer') 但它说 function 'lmer': 。

只是想知道是否有办法获取“X”插槽具体在做什么（或 S4 对象中的任何插槽）的源代码？或者有谁知道如何重现这个？非常感谢您的帮助和时间。

So here's my problem. I have a data set in R that I need to run a mixed effects model on. Here's the code:

When I try to run the script, it returns the following error:

I have tracked the error down to the "follows" parameter because when I just use stress_limit and word_position, it runs fine. If it helps, the data in "follows" are only 3 strings: n or l, consonant, vowel. I've tried replacing the spaces with _ but with no success. Is there something about the internal workings of the lmer() function that is preventing the use of "follows" in this case? Any help would be great!

For more info: intdiff contains numeric values, stress_limit is strings (Stressed or Unstressed) and word position is also strings (Word Medial or Word Initial).

EDIT: Here is a data sample that reproduces the error:

I tried the lme() function as well, but it returned this error:

The code in my original post is the exact code I'm using, minus the library(lme4) call, so I'm not leaving any information out I can think of.

My R version is 2.15.2

我正在尝试创建由两个分类变量解释的 MAM 的图形表示。我意识到我需要制作 2 个箱线图来显示每个解释变量的效果，因为两者之间的交互作用微不足道。我如何绘制箱线图，因为我似乎无法提取手段等，而且我不确定它是否是因为我使用的是 lmer。

到目前为止，这是我的代码....

我认为为了绘制这两个图表，我必须将 m4 用于上午或下午，将 m5 用于天气。

到目前为止我有这个但是后来卡住了，似乎无法在任何有帮助的地方找到答案，因为当我使用“$”时它会产生我不理解的错误。

谢谢

我正在尝试使用 lmer 函数来研究 3 种不同条件（cond=0、1、2）和目标的存在（target=False 或 True）之间的反应时间 (RT) 是否存在交互影响患者（患者）。

我写了以下等式：

我的问题是这个函数的默认截距是 cond = 0 和 target = False，而我希望截距是 cond= 0 和 target=True （以便查看 cond0*target= 之间是否存在显着差异真和 cond1*target=True)。

我将衷心感谢您的帮助。

这是我的输出

根据我的数据，选择的截距是cond0:targetTRUE，输出中的其他级别是cond1:targetFALSE和cond2:targetFALSE。

我正在尝试使用 lme4 库运行方差分析模型。我没有任何真实数据，所以我只是从 MASS 中提取了 Cars93 数据。虽然我可以让 lmer(...) 工作，但我似乎无法让 anova 表打印。

这是发生的事情：

任何帮助或建议将不胜感激。

是否有具有以下功能的 R 包：

(1) 模拟交互变量的不同值，(2) 绘制图表，展示交互中项的不同值对 Y 的影响，以及 (3) 与拟合 lmer( ) lme4 包的功能？

我查看了 arm、ez、coefplot2 和 fanovaGraph 包，但找不到我要找的东西。

我使用以下代码从数据库中删除了某些行：

这是一个虚拟变量，并且该特定虚拟变量的值为 1 的行已成功删除。（我在前后使用“dim”函数检查了数据库的尺寸；一切看起来都很正常。）

然而，在我这次使用新数据集“df2”运行回归模型后，我看到自由度急剧下降！这远远超过了删除的行数！

我想知道这怎么会发生。然后，我意识到新数据集有许多只有 NA 的行。在随机变量有缺失值的每一行，R 都有一整行的 NA 值。

在意识到上面的代码不是删除行的最佳方法后，我尝试了以下方法：

它似乎奏效了，因为我不再有同样的问题。但是你会说上面的这个新代码可能有一些我现在还没有真正意识到的（类似或其他）问题吗？

我正在处理 R 中blme 包中的 blmer() 函数，我的 sim 数据的 sigma.epsilon 不等于 1。但是，当我使用 blmer() 进行估计时，我发现在所有情况下输出残差方差都是 1。这与我最初的情况相反，因此结果并不好。我想知道是否有任何方法可以在这里调整剩余设置？例如通过设置任何先验信息？提前致谢！

我的代码在这里：

我的结果（见残差方差为 1.0000）：