问题标签 [lme4]

我正在做 lmer 并希望 R 忽略NA单元格。我已经把这个命令：

但是，在摘要中lmer，它仍然给我与 excel 原始文件中的行数相同的行数（在 obs 数中），所以似乎NA公式中没有排除这些行。您对此有什么建议吗？

我在这里和其他网站上尝试了一些建议，但似乎没有一个有效

我的模型看起来像

cont1和cont2都是连续的解释变量，random是8个水平左右的随机因子

对于标准 glm (ie)

我曾使用此代码进行预测

等等

我已经能够使用 glmer 获得预测，但我无法获得 cont2 的每个级别的预测，例如在标准 glm 中。

我试图复制这里建议的代码glmer - 使用二项式数据进行预测（cbind 计数数据）

但是我收到错误消息

如果有人能建议我如何解决这个问题，那将不胜感激

干杯，迪伦

我正在尝试运行多元线性混合模型，并且需要使用远程工作站来减少计算时间。

当我在我的个人计算机（R 版本 2.15.1、lme4 版本 0.999999-0、64 位 Unix）上从 lme4 运行 glmer() 时，我的模型运行正常。

当我切换到删除工作站（R 版本 3.0.2、lme4 版本 1.0-6、64 位 Linux）时，我的模型无法运行，并且收到错误和警告消息：

有趣的是，如果我对任一因变量运行 lmer()，模型将运行而不会出现错误或警告。

这是示例代码：

这是我个人电脑的输出：

以下是使用较新版本 lme4 的远程工作站的输出：

我想知道模型是否形成不正确，或者一个版本的 lme4 是否存在并发症。非常感谢任何想法或建议。谢谢！！

我无法理解 glmer 与泊松模型和 difflsmeans 之间结果的一些差异。这两个函数都来自 lmerTest 包。基本上，glmer 告诉我这两个系数在 p < 0.05 时很重要，但是当我使用 difflsmeans 时，它会给出不同的结果。我正在使用泊松链接和偏移量来模拟计数数据，并包括治疗的固定效应和批次实验的随机效应。

在其他使用高斯链接的分析中，lmer 和 difflsmeans 都给出了相同的结果。

将 difflsmeans 用于具有泊松链接的混合模型是否有效？

还有其他方法可以检查系数的重要性吗？（我知道这已经被问过，但我的意思是在这个分析的背景下）

提前致谢

来自 GLMER 的结果

来自差异的结果意味着

谁能告诉我为什么斜率系数在从具有随机斜率的 lmer 模型中提取的斜率系数与从适合同一数据集的 lmList 模型中提取的斜率系数之间存在偏差？

谢谢...

我正在处理金融/经济数据，以防您想知道以下某些系数的大小...我的一般问题与 R 中线性随机效应模型输出的参数系数的模拟有关。我是尝试使用模型系数和 R 中同一模型的方差-协方差 (VCOV) 矩阵生成 beta 系数的随机样本。我的问题是：为什么我会收到以下关于预期值平方根的错误mvtnorm{} 包中的 rmvnorm() 函数？我该如何处理这个警告/问题？

警告消息：在 sqrt(ev$values) 中：产生了 NaN

当我谷歌以下搜索字符串：“rmvnorm，”警告消息：在 sqrt(ev$values) 中：NaNs 产生，”我看到： http ://www.nickfieller.staff.shef.ac.uk/sheff-only /mvatasksols6-9.pdf 在第 4 页上，此错误表示“负特征值。”虽然，我从概念上或实际上不知道负特征值是什么，也不知道为什么会在这种情况下产生它们。

第二个搜索结果：[ http://www.r-tutor.com/r-introduction/basic-data-types/complex 2表示出现这个错误是因为试图取-1的平方根，即“不是一个复杂的值”（你不能取 -1 的平方根）。

问题仍然存在，随机生成的 beta 发生了什么，如何纠正？

sessionInfo() R 版本 3.0.2 (2013-09-25) 平台：x86_64-apple-darwin10.8.0 (64-bit)

使用以下软件包/版本 mvtnorm_0.9-9994、lme4_1.1-5、Rcpp_0.10.3、Matrix_1.1-2-2、lattice_0.20-23

我在尝试嵌入lmer函数时遇到了问题。这是一个使用来自lexdec. 如果我lmer直接在数据框上运行，是没有问题的。例如，假设我想查看词汇决策任务中的阅读时间是否作为Trial的函数而不同。有两种类型的词刺激，“动物”（例如“狗”）和“植物”（例如“樱桃”）。我可以计算动物词的混合效果模型：

但是，如果我将 lmer 模型嵌入到函数中（比如不要为每个级别的类键入相同的命令），我会收到一条错误消息。你知道为什么吗？任何建议将不胜感激！

我有一个模型：

dummy个人水平随机效应在哪里解释过度分散，即factor(1:nrow(dat))

运行此程序时，我收到以下我不明白的错误。这是否意味着我过度拟合了我的模型？

checkNlevels(reTrms$flist, n = n, control) 中的错误：每个分组因子的水平数必须 < 观察数

但是，当我与家人一起运行此模型时，poisson我没有收到此错误，例如

我知道个体水平的随机效应在高斯 GLMM 中甚至可能没有意义，但我想知道泊松示例是否对我隐藏了某些东西，表明模型过度拟合？

因此，我查看了许多类似主题的帖子，但似乎都不是我所需要的，或者我根本不了解他们提供的解决方案......所以就这样......

我用 lme4 运行了一个混合效应模型来查看一些黑猩猩数据。我有两个因素（攻击率；交配率）会影响我的依赖（喂食时间）。

我想制作两个散点图，显示每个预测变量和结果变量之间的关系，但我想画一条线，它来自模型估计（而不是 (lm(y ~ x )) 类型，它只给出一个简单的回归线，而不是基于完整 LMM 的回归线）。

我有一种感觉，这只有 ggplot2 才有可能，但我还没有真正弄清楚如何做到这一点。花了一天的大部分时间浏览书籍和论坛，我希望这可能有一个相当直接的答案，如果有人知道他们在做什么。

感谢您提前提供任何提示！

亚历克斯

我正在使用glmultiglmer 的包装器（二项式），它summary是：

This is glmulti 1.0.7, Apr. 2013. Length Class Mode 0 NULL NULL

按照在这个线程上所做的事情，虽然这是 for lmer， glmulti 在使用带有 lme4 的遗传算法时无限期运行，我得到与上面相同的结果。难道是因为版本已经改变并且包装必须以不同的方式进行？以下是虚拟代码（从上面的链接中提取）：

lme4版本：1.1.5

glmulti版本：1.0.7

《R 版本 3.0.2 (2013-09-25)》

解决方案

这有效：

packageVersion('lme4') '1.1.5'</p>

packageVersion('glmulti') '1.0.7'</p>

版本：3.1.0

仅供参考：来自包维护者：

“fitfunc 必须是函数的名称，因此您的其他调用（包括 glmulti 调用中的函数定义）无法正常工作。”

“你将 glmulti 的第一个参数命名为‘公式’，它必须是未命名的或‘y’……对不起。但是 y 是一个公式（如果传递一个字符串，它只是因变量）。”