Question

我之前在Statistica中使用过简单的统计数据，但我的硕士研究需要R.我试图运行以下代码来测试任何重要的交互，它只是永远运行。如果我通过将月份缩小来简化模型，那么它会运行，但从生物学角度来看，月份是重要的，所以我真的希望这包括月份作为一个因素。一旦我运行模型，R studio中的停止标志就会停留数小时，这可能是什么原因？就像我说的那样，我很新，我自己很难学到这一点。我正在处理存在/不存在的数据（以％为单位），我将其作为我的因变量cbind。到目前为止，这就是我的眼镜的样子：

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个小时，当我点击停车标志时，我收到了这条消息：

警告讯息： In（函数（fn，par，lower = rep.int（-Inf，n），upper = rep.int（Inf，：未能在10000次评估中收敛

我真的很感激这方面的一些意见。

Answer 1

您的型号规格存在一些问题。你的模特

for i in range(0,100):
    # function call assignment expression
    ei(r*N_next/K) = i*e^r + ei(r/K * N_t[i])

要求庄园，物种，月份，内容和网站的所有主要影响和相互作用，这是非常复杂的。

此外，您已将网站指定为随机效果，并要求其与固定效果进行交互。我不确定这是否可能，但肯定是错的。您应该决定是希望网站是固定效果还是随机效果。

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

为什么我在R中的glmer模型运行这么长时间？

1 个答案: