问题标签 [lme4]

我将 lme4 用于线性混合效果模型

有 8 个演员，我有 3 个时期。现在我想绘制（使用 ggplot2）8 个不同方面的不同线条（这似乎并不困难）。但是，我不知道如何为同一方面的每个演员绘制三条不同的线（因为三个不同的时期）。我使用了来自http://www.sagepub.com/long/chapters/rcode/82689_10rc.txt的以下代码：

如何为三个时期绘制三条不同的线？

这是数据：

我正在尝试使用包中的dredge功能MuMIn，如下所示。因为我想为不同的数据集运行它，所以我使用了一个函数，它接受协变量的向量，并且公式是从函数内的这个向量创建的。

当我以这种方式运行它时，我得到符号不是子集的错误。我认为这与vars在公式中使用有关，尽管 lmer 运行良好

术语错误（as.formula（公式（x）））：在为函数“术语”选择方法时评估参数“x”时出错：错误：“符号”类型的对象不可子集

但这很好用，原因我不完全理解，但与疏浚功能的范围有关。（我在这里可能完全不合时宜）

我在分析中所做的更类似于此，其中 lmer 和 dredge 在函数内运行

Dredgin 只有 lmer1 给出了这个错误：

仅疏浚lmer2，出现以下错误

有没有办法解决这个问题？这是疏浚功能的必要限制吗？

LME4 对象的 getSummary 函数似乎不再适用于 memisc 包，并且将 LME4 模型结果输出到乳胶不再有效。但是，通过 stargazer 导出对象的唯一选择是不允许长表。有没有人想出修复 mtable 格式的方法，或者将这些结果输出到 longtable 的替代方法？

编辑：

这是一个可重现的例子。

这曾经有效，但我相信 lme4 现在会生成一个新对象，该对象不再被 mtable/getsummary 命令识别。

我有来自参与者 ( part) 的受试者内生理数据，他们都在三轮 ( round) 中查看了刺激物（阅读报纸），每轮有五篇论文 ( )，并且在每一轮中，报纸上的访问次数 ( )paper不定visit. 我有两个固定因素（CONDhier和CONDabund）加上交互作用来预测生理状态（例如，EDA），这通常是自回归的。我试图通过随机效应来考虑生理上的个体差异（让我们暂时只考虑截距），并且可能会因另一种随机效应而在回合中感到疲劳。

因此，我想在 R 中运行的模型在 SPSS 中是：

现在，我明白了，虽然lme交叉项做得不好，lmer（处理交叉项没有问题）不能使用不同的协方差结构。我可以运行简单的 lme 模型，例如

但更复杂的模型超出了我的范围。我已经读过 lme 中的交叉术语可以通过随机定义来完成，例如

但这似乎阻止了我的 AR1 结构和 part*round 的第二个随机截距。而且我不太确定它是否与我的 SPSS 语法相同。

那么，有什么建议吗？尽管有很多关于 lme 和 lmer 的不同著作，但我找不到一个同时具有交叉术语和 AR1 的著作。

（此外， lme 上的语法似乎很模糊：从几个不同的来源我了解到 | 将左侧的内容分组到右侧的内容下， / 构成嵌套项， ~1 是随机截距， ~x 是随机斜率, 和 ~1+x 两者都是, 但似乎至少有 : 和 -1 定义我在任何地方都找不到。有没有教程可以解释所有不同的定义？)

我可以lmer在没有此错误的旧计算机上使用，但新版本（昨天下载）给了我此错误消息：

响应此功能：

我真的很希望能够使用较新的版本，因为我已经读到它比旧版本具有更少的错误收敛问题。任何人都可以提供有关如何使其在没有上述错误的情况下运行的任何帮助，我们将不胜感激。

我的数据可以在这里找到：https ://dl.dropboxusercontent.com/u/16881915/Rfile.txt

我有两个逻辑回归模型lme4。当我尝试anova在它们上使用时出现错误：

我已经看到这个问题报告了类似的错误消息，但根据那里的答案，错误似乎是由于将长表达式作为参数传递给anova，我没有这样做。

如何在我的模型上运行方差分析？

我有一个我想的非常简单的问题。在一组参与者的纵向实验中，让每个人在 7 个不同的时间给其他人打分，比如说，10 个变量（例如“这个人很讨人喜欢”、“这个人很呆板”等等）。如果我想为一个变量/响应获得某种感知器和目标方差，我会使用：

在这里，我们有一个数据帧“x”的因变量“Var1”，其中指定了第一个时间点（这也是 x 的变量）。

到目前为止一切顺利，这工作得很好。

现在正如我所说，我有多个响应和多个时间点。因此，我想使用 a）“for”循环或 b）lapply 来一次获取所有模型。

无论哪种方式，我都必须以某种方式“索引”因变量，无论是指定列位置（x[,10]10 是 Var1 的假定位置）还是变量本身（x$Var1）或（至少有点奇怪）粘贴或打印将变量名称代入公式 ( col.names(c[10])。

我想说的是，这两种方法都不起作用。我总是收到关于不同可变长度的错误。但是，正如我所写，我使用的是完全相同的列！

你们中有人有运行多个 lmer 的经验吗？

欢迎和赞赏所有想法！我希望我不是太不清楚，如果您需要任何进一步的信息，我很乐意提供，尽我所能。

干杯，艾尔

我有一个数据集，我正在使用 lme4 拟合混合模型回归。

接下来我想绘制结果拟合，因此使用来自http://glmm.wikidot.com/faq的示例代码，我为原始值生成了预测。

这很容易用 ggplot 绘制：

但如图所示，ggplot 绘制的回归线的截距与估计的回归截距相比是关闭的；至少 Nhc 线的情况是这样，它显然是负截距，而常见的估计截距是 0.54。

我的错误很容易（我认为）：使用 geom_smooth 导致的拟合与回归所说的不同，因为 geom_smooth 每次处理都单独获取数据，并且在没有拟合模型的上下文的情况下以面值计算。但我不知道如何正确绘制线条。

glmer我有一个使用in 函数拟合的 GLMM 对象，R并且想要执行 k 折交叉验证。对于简单的 GLM，我使用CVbinary了 pkg 中的函数，DAAG如下所示。

但是，当将 IndID 的随机项添加到模型中时，模型的 S4 类会导致错误（如下）glmer。

我一直在网上寻找，但无法找到CVbinary与 S4 对象类似的功能，但在手动编码之前想在这里仔细检查一下。

简而言之，（假设我正确解释了R错误）是否有对 S4 对象执行 k 折交叉验证的函数？

使用lmer我收到以下警告：

使用后产生此错误anova(model1, model2)。我试图使这个可重现，但如果我dput数据并再试一次，我的错误不会在 dput 数据上重现，尽管原始数据帧和新数据帧具有完全相同的str.

如果在干净的会话中再次尝试，并且错误重现，并且再次丢失dput

我知道我在这里没有给人们太多的工作，就像我说我很想重现这个问题一样。Cayone阐明了这个警告？