dispersion in summary.glm ()

J'ai dirigé un glm.nb par

glm1<-glm.nb(x~factor(group))

avec groupe étant une catégorie et x une variable métrique. Lorsque j'essaie d'obtenir le résumé des résultats, j'obtiens des résultats légèrement différents, selon que j'utilise summary()ou summary.glm. summary(glm1)Donne moi

    ...
Coefficients:
                    Estimate Std. Error z value Pr(>|z|)  
    (Intercept)       0.1044     0.1519   0.687   0.4921  
    factor(gruppe)2   0.1580     0.2117   0.746   0.4555  
    factor(gruppe)3   0.3531     0.2085   1.693   0.0904 .
    ---
    Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

    (Dispersion parameter for Negative Binomial(0.7109) family taken to be 1)

alors que summary.glm (glm1) me donne

    ...
Coefficients:
                    Estimate Std. Error t value Pr(>|t|)  
    (Intercept)       0.1044     0.1481   0.705   0.4817  
    factor(gruppe)2   0.1580     0.2065   0.765   0.4447  
    factor(gruppe)3   0.3531     0.2033   1.737   0.0835 .
    ---
    Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

    (Dispersion parameter for Negative Binomial(0.7109) family taken to be 0.9509067)

Je comprends la signification du paramètre de dispersion, mais pas de la ligne

(Dispersion parameter for Negative Binomial(0.7109) family taken to be 0.9509067).

Dans le manuel, il est dit que ce serait la dispersion estimée, mais cela semble être une mauvaise estimation, car 0,95 n'est pas proche de 0,7109, ou la dispersion estimée est-elle quelque chose de différent du paramètre de dispersion estimé? Je suppose que je dois régler la dispersion dans summary.nb(x, dispersion=)quelque chose, mais je ne suis pas sûr, si je dois régler la dispersion sur 1 (ce qui donnera le même résultat que summary()si je devais insérer une estimation du paramètre de dispersion, dans ce cas à menais summary.nb(glm1, dispersion=0.7109)ou quelque chose d' autre? ou suis - je bien avec juste en utilisant le summary(glm1)?

summary.glm"negbin"summary.glmdispersionsummary.glm glm$\phi$$\phi$familyglm.nb"Negative Binomial(theta)"summary.glmsur le modèle monté par glm.nb, le code en

if (is.null(dispersion)) 
    dispersion <- if (object$family$family %in% c("poisson", 
        "binomial")) 
        1
    else if (df.r > 0) {
        est.disp <- TRUE
        if (any(object$weights == 0)) 
                warning("observations with zero weight not used for calculating dispersion")
            sum((object$weights * object$residuals^2)[object$weights > 
            0])/df.r
    }

"poisson""binomial"$\phi$summary.negbin

$\phi$dispersion

Deuxièmement, vous comprenez mal la sortie. Quand tu vois

Negative Binomial(0.7109)

$\hat{\theta}$$\phi$

$\phi$$\phi = 1$summary.negbin

summary(glm1, dispersion = 0.9509)

negbin$\phi$

$\theta$$\frac{1}{\theta}$$1$$\frac{1}{\theta}$$E$$Y$$E$$\mu E$$\mu$

attente

& variance

Comme le souligne @Momo, le paramètre de dispersion est une tout autre chose, que vous laisseriez varier pour effectuer une estimation de quasi-vraisemblance. Pour le modèle binomial négatif et le modèle (vrai) de Poisson, il est correctement fixé à une valeur de un.