r - 为什么我在 R 中的 glmer 模型需要这么长时间才能运行？

Question

我以前一直在 Statistica 中使用简单的统计数据，但我的硕士研究需要 R。我正在尝试运行以下代码来测试是否存在任何重要的交互，并且它一直在运行。如果我通过取出月份来简化模型，那么它就会运行，但从生物学上讲，月份很重要是有道理的，所以我真的很希望它能运行，包括月份作为一个因素。运行模型后，R studio 中的 parking 标志会持续数小时，这可能是什么原因？就像我说的那样，我很新，而且我自己学习这个真的很难。我正在处理存在/不存在数据(作为 %)，我将 cbind 作为我的因变量。到目前为止，这是我的 coad 的样子:

library(car)
library(languageR)
library(AICcmodavg)
library(lme4)

Scat <- read.csv("Scat2.csv", header=T)

attach(Scat)

names(Scat)

y <- cbind(Present,Absent)

ScatData <- glmer(y ~ Estate * Species * Month * Content * (1|Site) + Min + Max,family=binomial)

summary(ScatData)

一旦我开始运行实际模型，我什至不做总结，因为 R 还没有完成对实际模型结果的计算。我运行模型大约 4 小时，当我点击 parking 标志时，我收到了这条消息:

Warning message: In (function (fn, par, lower = rep.int(-Inf, n), upper = rep.int(Inf, : failure to converge in 10000 evaluations

如果您能就此事提供意见，我将不胜感激。

Answer 1

您的模型规范存在一些问题。你的模型

y ~ Estate * Species * Month * Content * (1|Site) + Min + Max

要求所有主要影响和地产、物种、月份、内容和地点的相互作用，这是非常复杂的。

此外，您已将网站指定为随机效应，并询问其与固定效应的交互作用。我不确定这是否可能，但这看起来肯定是错误的。您应该决定是希望网站成为固定效果还是随机效果。

如果您发布一个最小的可复制示例，我可以提供更具体的建议。