问题标签 [lmer]

我正在使用 HLMdiag 包中的工具来查看由lmer. dotplot_diag我对绘制厨师距离的功能有疑问。

这是我正在做的一个例子。

我得到以下情节：

问题Subject：有没有办法用它们的值而不是它们的索引来标记有影响的点？即，它将显示类似 M13 或 F10 的内容（Subject变量中的值，见下文）。

预先感谢您的帮助。

我对 lmer 做了一些更改。它可以正常工作，但是我无法摆脱运行程序时弹出的警告消息。我添加了以下选项，这些选项允许程序在不停止但带有警告消息的情况下运行。我相信这是 check.nobs.vs.rankZ = "warningSmall" 部分。我怎么能摆脱这个，有什么建议吗？谢谢你。

来自 R 的警告信息：

对于模拟研究，我对比了不同 LMEM 的重复测量能力。我想指定一个模型，其中允许随机截距和斜率相关，而其中一个不允许。但是，当我比较这两个模型时，它们似乎完全一样。当我在两个模型上运行方差分析时，没有相关性的模型甚至有一个更多的自由度，而它应该是相反的。

数据：

我试图解释 GLMM 中的自相关。我的响应变量是布尔值，它表示一组蜂巢的生命周期中是否存在 en 事件。我试图用一组描述每个巢的状态的数值变量来预测这种事件的概率。因此，我在具有嵌套作为随机效应的广义模型中使用了二项分布（使用 glmer()）。然而，这些事件是自相关的，所以我的残差中有一个非常可怕的模式。如果我在没有随机效应的错误中使用高斯分布，我会使用带有 gls() 的相关结构来估计一些相关参数，但这在这种情况下不起作用。我一直在寻找将这种自相关包含在 GLMM 中的方法，但我似乎没有做对。我发现可以使用函数 ts() 和 diff() 转换数据集，以允许模型包含响应的先前值作为预测变量。这适用于线性模型 lm()，使残差更好。

但是， lmer() 和 glm() 都不接受这些函数的输出。问题似乎是转换使某些值为负数，而 glm() 不接受二项式模型的负值（这是有道理的）。我已经读到可以考虑自相关 glm()，这已经是一种改进，但我无法让它发挥作用。我还读到 glmmPQL() 可以包含相关结构。该模型运行，但它并没有改善我的残差模式。它们似乎仍然是自相关的。

我尝试了不同的相关结构，但它们似乎都不起作用。

最后，我尝试了 dyn 包，它应该允许回归函数处理时间序列。但是再一次，该函数不会使用转换产生的值运行。

总而言之，我需要运行具有时间相关性的 GLMM，但我找不到这样做的方法。我将非常感谢一些帮助。

干杯！！！

我目前正在使用 lmer 和 lmerTest 进行分析。每次我在随机结构中添加效果时，在运行 summary() 时都会出现以下错误：

如果我离开这样的结构：

没有错误运行摘要（TRT5ToVerb.lmer2），返回AIC，BIC，logLik偏差，随机效应的估计，固定效应的估计及其对应的p值等，等等。

因此，尽管对象 TRT5ToVerb.lmer3 存在，但显然在我运行 lmerTest 时发生了一些事情。两者的唯一区别是随机结构：(1+Condition|Participant) vs. (1|Participant)

我的数据的一些特征：

1. Condition 和 Group 都是分类变量：Condition 包含 3 个级别，Group 2
2. 因变量 (TRT5ToVerb) 是连续的：它对应于以 ms 为单位的读取时间
3. 这是一个重复测量实验，每个参与者有 48 次观察（参与者 = 28）

我读过这个威胁，但我看不到明确的解决方案。是不是我必须将我的数据框转换为长格式？如果是这样，那么我如何在 lmer 中使用它？我希望不是这样。

谢谢！

免责声明：我既不是 R 专家，也不是统计学专家，所以请耐心等待。

我已经使用lmerTest和 using运行了以下模型lme4：

使用我在输入命令lmerTest 时收到以下错误：summary()

我看到这对其他用户来说已经是一个问题，并且一个用户能够绕过运行的问题lsmeans()。当我尝试 lsmeans 时，我得到了错误：

查看协方差矩阵时，我没有看到任何 NA。请注意，如果我简单地反转组因子中的对比，我就可以运行这个模型。我很难理解为什么会这样。

当我使用 lme4 而不是 lmerTest 运行相同的模型时，我能够获得 summary() 的所有输出，但没有 p 值（如预期的那样）。pvals.fnc 在 lme4 中已停产，我还没有找到替代方案。另外，最好以与我成功使用 lmerTest 的其他模型相同的方式估计模型 2 的 p 值。

有谁知道我此时应该做什么？任何帮助将非常感激！

我正在使用lme4R 中的包并尝试适应随机斜率和随机截距模型。如果有人能在我运行随机斜率和随机截距模型时帮助我理解这个错误以及如何处理这个错误，我会非常有帮助：

非常感谢

我想使用try()ortryCatch()或这样的函数来检测我的模型中是否存在名为“fit1”的错误。如果模型很好，我想使用“fit1”，否则我想使用“fit2”

你知道怎么做吗？我没有添加任何数据，因为我的问题可能很容易解决，但如果需要，我可以上传它们。

谢谢！

我正在尝试使用 pbkrtest 包中的函数 PBmodcomp 对两个 lmer 模型进行模型比较。但是我收到以下错误。

我的数据在这里可用：https ://www.dropbox.com/s/oweyw767qtpbqot/Data.txt

我正在运行的代码是

我该如何解决这个问题？

谢谢

我曾经使用下面的代码来计算lmer模型的标准化系数。但是，随着新版本的 lme 返回对象的结构发生了变化。

如何调整功能stdCoef.lmer使其与新lme4版本兼容？