准泊松模型的 Akaike 准则 (AIC)

Question

我在 R 中运行准泊松 glm 模型，我想计算 AIC 标准。我运行下面的代码，但我不确定它是否正确。另外，我发现 AIC 减小而自由度增加这一事实很奇怪。

fqaic <- function(model) {
  loglik <- sum(dpois(model$y, model$fitted.values, log = TRUE))
  phi <- summary(model)$dispersion
  qaic <- -2*loglik + 2*summary(model)$df[3]*phi
  return(qaic)
}

如有任何建议，我将不胜感激。

Answer 1

这里有一些注释，表明除其他外，bbmle

、

AICcmodavg

和

MuMIn

包中存在计算 qAIC 的现有方法。 qAIC 的正确公式是

-2 * logLik/dispersion + 2 * df

Answer 2

这是对最高评价答案（作者：Ben Bolker）的回答。由于我刚刚创建了帐户，因此无法直接发表评论。

我读过 Ben Bolker 的文档“Dealing with quasi-models in R”，但我很难理解它。

它说：“对模型进行两次拟合，一次使用常规似然模型（族=二项式、泊松等），一次使用准变体——从前者中提取对数似然，从后者中提取离散参数”

这是否意味着我可以使用上面公式中的这些参数（-2 * logLik/离散度 + 2 * df）来计算拟泊松模型的 qAIC？

提前非常感谢！