问题标签 [lmer]

我正在使用留一法来评估排除了一个数据点的模型预测该数据点（循环所有数据点）的效果。下面的代码已经成功地在基本相同的数据上运行，但 DV 略有不同，所以我很难理解为什么我会得到我得到的错误。这是相关的代码块：

dataFull2具有与 相同的列pennTrim，只是具有更多行。该功能的非常标准的使用predict()。我收到此错误：

关于可能导致此错误的任何想法？我可以使用基本相同的代码和相同的数据帧交换（来自不同来源的测量），并且代码没有错误logDur.PENN。logDur.Manual

尝试运行混合效果模型时出现以下错误：

这是在使用 caret 包中的 findLinearCombos 来查看是否有任何可线性组合的列之后：

关于如何删除此警告并找到不可线性分离的列的任何想法？

我想问题是 Intercept 列是线性组合？

我在网上冲浪并没有找到满意的答案。

如何为测试数据集中的每个观察结果从 lmer 对象生成预测区间？

interval与 predict 没有任何争论。sim显然不起作用，所以我需要使用 LMER 附带的 mcmc 函数并从分位数中提取？

当我尝试使用 lmer 进行预测时出现以下错误

这就是我的公式的样子

我能够总结模型，但我不能将它传递给预测函数。

我希望能够根据预测矩阵的列自动生成一个公式，然后将其传递给 lmer。我该怎么做？

我在 lmer 中拟合了一个模型，其中包括一个三向交互项，其中两个变量是分类变量，一个是连续变量。我正在尝试恢复所有斜率和截距的均值参数化，而不是 lmer 中默认的效果参数化，但我坚持正确的编码。例如（不包括随机效应），使用 iris 数据集我做了一个额外的分类变量（土壤），并用物种、萼片宽度和土壤拟合模型：

给出的输出

最后三个斜率值（Soil = Y）仍在效果参数化中，我无法计算出正确的编码来获得平均值。我认为这可能吗？任何建议将不胜感激。

我有两个没有随机效应的回归模型：一个是使用 lm 的 OLS，另一个包括使用 nle 的系数乘法。我希望为两者添加个人级别的随机效果。我已经设法使用 lme4 包为 OLS 函数执行此操作，但无法为乘法模型找到执行此操作的方法。

以下代码生成了一个与我正在处理的具有相似结构的数据集：

...

a、b 和 c 代表三个 1:5 维度尺度上的分数。a2 到 c5 是虚拟变量，代表相同比例的 2:5 水平。每个人 (id) 有 10 个观察值。val 是我希望使用回归模型预测的分数的代理。（但实际数据中的值可能与此处的结构不对应。）

我有两个没有随机效应的回归模型。一种是使用 12 个虚拟变量作为 val 的预测变量的常规 OLS：

第二个假设水平之间的相对距离对于三个维度（a，b，c）是共享的，但是维度在尺度方面不同。剩下 6 个系数（cA、cB、cC、cL2、cL3、cL4）+ 截距。

由于每个人有 10 个观察值，我们不能期望它们完全独立。因此，我希望在变量 id 定义的个体级别上添加随机效应。我找到了一种使用 lme4 包的方法：

问题是是否可以使用类似于乘法的回归模型（可能使用 lme4 或 nlme 包）对 id 变量添加随机效应？公式应该看起来像

有什么建议么？

我正在分析一些重复测量的药物试验数据，但我不确定在使用多面 ggplots 时如何绘制 lmer 结果。我已经从主数据集中绘制了各个斜率的初始图，但我正在按性别分别进行 lmer 分析。

使用公开可用的数据，与我的四个治疗组相比，只有两个治疗组，这是下面的可复制示例。它使用reshape2、lme4和ggplot2包。

有关数据集的信息在这里。

如何在两个方面从 lmer 绘制正确的摘要线，这两个方面的不同是Treatment？

在我的现实生活示例中，我还将分析女性，因此每个方面将绘制两组线。

我正在使用包将统计测试ReporteRs的R结果和绘图输出到Word 2010.

我正在使用包lmer中lme4做一些分析，我想将结果输出summary(foo.lmer)到Word。我希望捕获lmerMod对象中的所有内容并将其放入 Word。例如，这summary(foo.lmer)是我得到的结果：

如何使用固定宽度和硬返回存储此信息，以便我可以使用 ReporteRs 输出此信息？我可以复制和粘贴，但我正在尝试自动化输出，而且我经常重写输出的 Word 文档。

固定效果输出是最重要的部分，但是将所有输出复制到 Word 中会非常有用。

更新和解决 感谢在回复中使用 capture.output 的建议，这段代码对我有用，仍然使用 ReporteRs：

如果您想将您的回复作为答案，我可以接受并投票。

我在尝试适合 R 的嵌套模型时遇到了很大的问题。

我有 2 个条件的响应时间实验，每个条件 46 人，每个条件 32 个测量值。我希望将措施嵌套在人和人嵌套在条件中，但我无法让它发挥作用。

我认为应该有意义的代码是：

但是，我得到的只是一个错误：

不幸的是，我什至不知道从哪里开始寻找问题所在。有任何想法吗？请，请，请？=) 干杯！

当我使用 lmer 时有一个奇怪的行为：当我将使用 lmer 的拟合保存到一个对象中时，比如说fit0，使用 lmer，我可以查看摘要（输出未显示）：

如果我使用 save.image() 保存对象，关闭会话并重新打开它，摘要会给我：

如果我再次运行模型，我会得到预期的摘要，但如果我关闭会话则会丢失它。
怎么了？我怎样才能避免这个错误？

感谢帮助。

环境及版本：
Windows 7
R version 3.1.2 (2014-10-31)
GNU Emacs 24.3.1 (i386-mingw-nt6.1.7601)/ESS

这是一个最小的例子：

重新启动并添加库 lme4 后：