Coder un modèle

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:

---

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.

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.

Ah oui, excusez-moi :oops:

Donc voici les 5 premières lignes :

Code:

---

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:

---

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.

Bonsoir,

Trop de parentheses peut-etre. Essayez comme ca :

Code:

---

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

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.