问题标签 [lmer]

我运行了一个模型，其中反应时间作为我的 DV 和 PWI 条件作为固定因素之一。

我将 contr.sum 用于所有固定因子。我运行以下模型来寻找三个不同 IV 的反应时间分数的差异（每个 2 个级别）。

输出如下：

输出（偏差编码）显示 PWI 条件的正 t 值为 2.067，其中非同源被编码 (-1) 和同源被编码 (+1)。因此，我将此解释为同源词的命名时间比非同源词要长得多。但是，这些结果与原始均值不匹配。非同源的手段比同源的手段长得多（就反应时间而言）。我该如何解释这个结果？

我正在尝试估计和模拟各种模型。我可以安装许多模型，sapply但不知何故，我无法访问输出。

我已经做到了，我可以通过进一步应用来使用模型sapply，例如，模拟：

但是，现在我需要提取fixefand ranefof sims。通过打印models或者sims它们看起来很像lmer输出，但它们不是。在尝试访问固定效果时，我收到这样的错误消息，就像使用 lmer 输出一样，这似乎是合乎逻辑的：

关于如何访问输出（或将其转换为能够访问它）的任何想法？

谢谢

这是输出sims：

我有一个关于从 R 中装有lmer( lme4) 的模型中提取随机效应的（原始）模型矩阵的问题。更具体地说，我想获得一个数据框或矩阵，其中包含随机效应项中涉及的所有变量. 问题更加复杂，因为该矩阵的某些条目为零。

我通常通过访问稀疏模型矩阵 ( Zt) 来提取这些矩阵getME，然后通过其维度将其转换为常规矩阵（见下文）。然而，只要（原始）模型矩阵包含零，这就会导致问题，因为Zt只包含非零元素。

这是一个示例，一个简单的混合效应模型，其中x1正常且x2包含五个正好为零的值：

我使用 拟合两个模型lmer，一个与x1另一个x2作​​为预测器：

在这里，我访问Zt了拟合模型对象的插槽。下面的代码演示了Zt不包含x2. 结果，我非常简单地转换为常规矩阵会引发错误。

这是我认为我可以做的。lmer似乎也保存了原始矩阵，只要只有一个簇变量，它似乎就可以很好地工作。

但是，该mmList插槽在网上的记录很差，我几乎没有提到人们将它用于编程。到目前为止，访问Zt似乎是更常见的选择。

是否可以从Zt即使原始模型矩阵包含零如果没有，那么我应该期待mmList什么？

我是混合效果模型的新手，我需要你的帮助。我在ggplot中绘制了下图：

但是，我想表示一个混合效应模型而不是lmin geom_smooth，因此我可以将其包含SITE为随机效应。

模型如下：

我已经包括TRTYEAR（治疗年份），因为我也对效果的模式感兴趣，对于某些群体来说，它可能会随着时间的推移而增加或减少。

接下来是我迄今为止从模型中提取绘图变量的最佳尝试，但只提取了TRTYEAR= 5、10 和 15 的值。

欢迎提出用完全不同的方法来表示这种分析的建议。我认为这个问题更适合 stackoverflow，因为它是关于 R 中的技术细节，而不是背后的统计数据。谢谢

我目前正在研究具有大量变量和交互项的多级模型。到目前为止，我已经使用 LMERConvenienceFunctions-package 的 plotLMER3d.fnc 函数来获取其中一些交互的 3D 图。不幸的是，当我向回归添加更多变量时，我收到一条错误消息：

fit3_plot <- lmer(OTIF ~ FtO* Variability + FtO* COLT2 + Products *FtO + FtO *Age4 + Products *FtO + Orders2*FtO + Age2*Variability +COLT2*Variability + COLT2*Age2 + COLT2*Products + COLT2*OpC + (1 | BU) + (1 | BU:RBU_SBU)，数据 = Data4)

plotLMER3d.fnc（fit3_plot，pred = “可变性”，intr = “COLT2”，plot.type = “persp”，phi = 25，theta=90）

gzfile（文件，“wb”）中的错误：无法打开连接

另外：警告信息：

In gzfile(file, "wb") : cannot open compressed file 'C:\Users\XXX\AppData\Local\Temp\RtmpKEw4FX/lmer___OTIF_FtO_Variability__FtO_COLT2__Products_FtO__FtO_Age4__Products_FtO__Orders2_FtO__Age2_Variability__COLT2_Variability__COLT2_Age2__COLT2_Products__COLT2_OpC__WWW1_BUWWW__WWW1_BU_RBU_SBUWWW___Data4___Variability_COLT2.rda', probable reason 'No such file or directory'

但是，当我从回归中删除前 6 个交互时，我能够获得一个图。

有人知道如何解决这个问题吗？我会很感激！

最好的，

基督教

使用这个虚拟数据集：

这有效：

但这不会：

对于后者，我只得到

从plot命令。

当我查看int.nlme和int.lmerwith 时str，它们看起来确实不同，但我无法弄清楚问题所在。非常感谢任何输入。

我正在使用lmerfrom lme4package 来计算方差分量的置信区间。

当我拟合模型时，会出现警告消息：

我搜索了很多以了解为什么会发生错误并最终做出决定there is difference between error and warning in the R world。

我想计算模型参数的置信区间并运行显示错误的代码：

是否有另一种方法可以使用 lmer 获得随机效应的 CI？

编辑 ：

我正在调用（ 1.1-9）中的lFormula函数来获得随机效应的稀疏模型矩阵的转置（）。问题是未使用的级别被丢弃在输出矩阵中。有没有办法解决这个问题？lmerlme4Zt

我意识到这个问题可能看起来很尴尬，所以这里是 GitHub Gist 上的一个示例，它试图拟合一个多成员模型，其中一个因子的级别比观察值多。