单季错误占用模型的健康包装无标记和

Question

您好我正在努力为单季制作适合度分析，使用未标记的包装并使用包装AICcmodavg进行评估。您可以将原始模型下载为RDS here

library(unmarked)
library(AICcmodavg)
BestMylu <- readRDS("best2.My.Lu2.rds")
obs.boot <- mb.gof.test(BestMylu, nsim = 5000)

我收到以下错误：

  

data.frame（det.hist，preds.psi）中的错误：     参数意味着行数不同：123,111   另外：警告信息：   1：由于缺少相应的协变量，一些观察结果已被丢弃。   2:12个网站由于缺少数据而被丢弃。   3：由于缺少相应的协变量，一些观察结果已被丢弃。   由于数据丢失，已丢弃4:12个站点。

我知道我得到的是因为原始data.frame中的NA数据，我可以回去删除那些行，但我不得不重做25个不同的模型，我宁愿克服这个错误。

有没有办法克服这个函数中的错误，还是可以使用另一个函数来获得合适的好处？

Answer 1

此错误更像是您的响应变量中缺少数据 - 因此，当执行实际响应和预测响应之间的比较时，行数不相等 - 检查列车上的响应变量行数以及预测输出中的行数，看看它们不相等的原因。