Discussion :

[Delphine Van Ruysten]

Bonjour,

Voilà, je dois coder un modèle de régression. D'habitude, j'y arrive sans trop de problèmes mais là j'ai un message d'erreur qui revient et je n'arrive pas à résoudre le problème.

Voici le modèle (et je dois interpréter les paramètres) :

http://imageshack.us/photo/my-images...515154134.png/

y: nombre de jours d'absence
z: sexe (0 ou 1)
x: age.

Et voici l'erreur que j'obtiens systématiquement :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
[4] ERREUR:  
  formule de modèle incorrecte dans ExtractVars
[5] ERREUR: the model GLM.1 is no longer available

Pourtant je mets bien l'élément "family=poisson(log)".

Comment puis-je résoudre ce problème?

Merci,

D.V.R.

[A. D.]

Bonjour,

Pour nous aider à vous aider (), il est toujours préférable de poster votre code (surtout qu'ici il semblerait que le modèle soit mal spécifié or difficile de dire ce qui coince sans savoir comment vous l'avez écrit...) et si possible, un extrait de données sur lequel on puisse reproduire le problème.

Merci !

Bonne continuation à vous !


Cordialement,

A.D.

[Delphine Van Ruysten]

Ah oui, excusez-moi

Donc voici les 5 premières lignes :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
 y sex age
 4   0  50
 0   0  21
 6   0  43
 0   0  30
 0   0  26

Et enfin le code que j'ai utilisé :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
GLM.1 <- glm(y ~ (sex)+(sex)*(age-39), family=poisson(log), data=Dataset) 
 
summary(GLM.1)

C'est assez trivial je sais, mais j'ai fait le même en mettant juste une normale à la place de la poisson et ça renvoie les estimateurs comme d'habitude. Or pour le coup, c'est une loi poisson qui est demandée.

Cela dit le modèle à suivre que j'ai mis en lien dans le 1er post semble exigeant et l'erreur vient peut-être du fait que j'ai mal retranscrit le modèle.

Merci de votre aide,

D.V.R.

[vchouraki]

Bonsoir,

Trop de parentheses peut-etre. Essayez comme ca :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
GLM.1 <- glm(y ~ sex + sex * I(age-39), family=poisson(log), data=Dataset) 
 
summary(GLM.1)

La fonction I() permet une operation arithmetique dans la formulation du modele.

HTH

Vincent

[Delphine Van Ruysten]

Merci beaucoup, je ne connaissais pas la fonction I()!

Le modèle fonctionne, je reçois les résultats.

Ce modèle correspond donc bien à la formule de mon 1ier post? Je pensais qu'il y aurait peut-être plus de difficultés.

Encore merci,

D.V.R.