Entre quasipoisson et negative binomial

Postez ici vos questions, réponses, commentaires ou suggestions - Les sujets seront ultérieurement répartis dans les archives par les modérateurs

Modérateur : Groupe des modérateurs

Serge Faye
Messages : 34
Enregistré le : 14 Mar 2007, 09:49

Entre quasipoisson et negative binomial

Messagepar Serge Faye » 14 Mai 2011, 10:05

Bonjour,
pour étudier les tendances dans ma base, j'ai transformé le fichier d'individus en fichier de comptage. voici à quoi ressemble mon jeu de données
nbpst annee lage sexe activite
1 2001 -30 ans Homme Administration
1 2001 -30 ans Femme Administration
7 2001 -30 ans Femme Administration
15 2001 -30 ans Homme Administration
1 2001 -30 ans Homme Administration
1 2001 30-40 ans Femme Administration
1 2001 30-40 ans Homme Administration
1 2002 30-40 ans Femme Administration
2 2002 30-40 ans Homme Administration
7 2003 30-40 ans Femme Administration
4 2003 30-40 ans Homme Administration

Le tableau résume le nombre de problèmes (nbpst) par année, age, sexe et activité.

Je modélise nbpst (qui suit une loi de poisson) avec glm famille = poisson. Les variables explicatives sont annee, lage, sexe et activité.

Le rapport de la variance/moyenne est égal 46 indiquant une surdispersion. Surdispersion qui est confirmée par le test de Dean.

Donc je passe du poisson simple au quasipoisson et au negative binomial
Pour le quasipoisson j'ai :

(Dispersion parameter for quasipoisson family taken to be 26.17042)

Null deviance: 86850 on 3678 degrees of freedom
Residual deviance: 68188 on 3654 degrees of freedom
AIC: NA

Number of Fisher Scoring iterations: 6

et pour le negative binomial:
(Dispersion parameter for Negative Binomial(0.7512) family taken to be 1)

Null deviance: 5082.1 on 3678 degrees of freedom
Residual deviance: 3991.2 on 3654 degrees of freedom
AIC: 24429

Number of Fisher Scoring iterations: 1


Theta: 0.7512
Std. Err.: 0.0169

2 x log-likelihood: -24376.9570

Mes questions sont les suivantes:
1/Peut -on dire qu'avec le quasipoisson on confirme la présence de surdispersion (Dispersion parameter for quasipoisson family taken to be 26.17042)?
2/ Peut-on dire que le negative binomial modélise mieux mes données par rapport au quasipoisson (Dispersion parameter for Negative Binomial(0.7512) family taken to be 1)?
3/Compte tenue de la nature de mon fichier où la variable à expliquer est le résultat d'un dénombrement entre différents strates (qui sont mes variables explicatives)imbriqués les uns à la suite des autres est ce que la régression de poisson classique est le mieu adaptée?
Par avance, je vous remercie du temps que vous accorderez à mon message.
Cordialement,

Stéphane Laurent
Messages : 1557
Enregistré le : 05 Déc 2006, 19:07

Messagepar Stéphane Laurent » 15 Mai 2011, 09:18

Bonjour,
Est-ce toi qui a imposé le "1" comme paramètre de dispersion de la binomiale négative ?
Pour moi il n'y a aucune raison de modéliser un dénombrement par une loi de Poisson en général.

Serge Faye
Messages : 34
Enregistré le : 14 Mar 2007, 09:49

Messagepar Serge Faye » 15 Mai 2011, 09:27

Bonjour,
Je ne me suis pas imposé la règle du "1".
J'ai utilisé la modélisation de poisson par ce que ma variable etait du comptage et son histogramme suit la forme d'une distribution de poisson, à savoir en forme de L.

Stéphane Laurent
Messages : 1557
Enregistré le : 05 Déc 2006, 19:07

Messagepar Stéphane Laurent » 15 Mai 2011, 11:24

Ma remarque sur Poisson était une "réponse" à ta question 3. Ce "1" est une estimation ? J'étais surpris que ça tombe sur un nombre entier.

Serge Faye
Messages : 34
Enregistré le : 14 Mar 2007, 09:49

Messagepar Serge Faye » 15 Mai 2011, 14:49

Je vois mieux maintenant,
Je me demandais si un autre type de poisson (par exemple du poisson sur données regroupées) n'était pas mieux adapté, car je n'ai pas des individus, mais des comptages obtenus par regroupement. mes variables de regroupement étant mes variables explicative.
Pour le "1", c'est juste le résultat du modèle. je n'ai eu à faire aucune intervention.


Retourner vers « Questions en cours »

Qui est en ligne

Utilisateurs parcourant ce forum : Google [Bot] et 1 invité

cron