R 如何从拟isson模型中提取拟BIC（QBIC）？

我想通过使用准BIC（QBIC）作为信息标准拟合准泊松模型，使用GLMULTI进行特征选择。我使用以下代码：

glmulti(y ~ ., data=data,
               level=2, fitfunction=glm, family = quasipoisson, crit="useBIC", confsetsize=128, 
               method = "g", plotty = T)

其中useBIC来自AICcmodavg包。然后我收到错误消息：

Mean crit= NA
Error in plot.window(...) : need finite 'xlim' values
In addition: Warning messages:
1: In min(x) : no non-missing arguments to min; returning Inf
2: In max(x) : no non-missing arguments to max; returning -Inf
3: In min(x) : no non-missing arguments to min; returning Inf
4: In max(x) : no non-missing arguments to max; returning -Inf

对我来说，这表明该函数试图计算一个规则的BIC，但这不起作用，因为该模型没有可能性。我这么认为是因为当我在准泊松模型上调用

useBIC

时，我收到了NA

---

我最终想做的是使用GLMULTI和准泊松模型选择特征（变量）。我使用了GLMULTI拟合负二项模型和BIC，效果很好，但我想看看拟合准泊松模型是否会产生一组不同的特征。

在我看来，这更像是设置了一种动态绘图方法。是的，这是函数的一部分功能。但我不认为这个错误仅仅与绘图有关。它提醒我，有一次我错误地编码了一个自动轴刻度算法，产生了这种类型的错误。错误消息清楚地指向产生错误的

plot.window（）

函数。我不熟悉

glmulti

软件包，但该错误可能与一体式函数有关。我承认这可能是巧合，但你可能想排除这种可能性。我不知道

glmulti

包，但可能还有其他的可能性。也许这也是一个询问如何计算准bic的好方法。