Constante de la loglik dans AIC()

[Pierre Bady]

re ...

Pour en avoir discuté avec David Anderson, je prend ça comme un résultat acquis et appuyé.

ok pour ce point.


sinon, pour revenir dans le topic initial, les familles de distribution classique ont une valeur 'aic' qui pourrait te renseigner sur leur mode de calcul ... dans la perspective d'un recalcul manuel comme le suggère Renaud.

Code : Tout sélectionner

> ?binomial
> binomial()$aic
function (y, n, mu, wt, dev)
{
    m <- if (any(n > 1))
        n
    else wt
    -2 * sum(ifelse(m > 0, (wt/m), 0) * dbinom(round(m * y),
        round(m), mu, log = TRUE))
}
> poisson()$aic
function (y, n, mu, wt, dev)
-2 * sum(dpois(y, mu, log = TRUE) * wt)
>

hth


pierre


PS: on est retourné sur le droit chemin ... enfin, je pense ... :p

[Logez Maxime]

Pour poursuivre le post de Pierre,

Les fonctions qui utilisent ces familles comme glm utilisent les fonctions aic de chaque famille mais pas seulement. Dans glm.fit tu trouves les lignes de commande suivante :

Code : Tout sélectionner

    aic <- family$aic
   aic.model <- aic(y, n, mu, weights, dev) + 2 * rank

Maxime

[Nicolas Péru]

ok...ça fait effectivement pas mal de chose à fouiller pour s'assurer du calcul de l'AIC.

En tout cas merci pour toutes ces précisions/échanges même si on a fait quelques détours HS...que le cirad nous pardonne !! :P
Mais c'est la faute de Pierre tout ça ! :D

Nicolas