Les posteurs les plus actifs de la semaine


regression polynomial R

Aller en bas

regression polynomial R Empty regression polynomial R

Message par leaticia le Mar 14 Juin 2011 - 6:42

Bonjour,
pour ma régression polynômial d'ordre deux j'ai le message d'erreur suivant que je ne comprend pas:

reg<-lm(X~polym(Y,Z,degree=2)) avec y variables binaire {0,1} et Z,X variables continues

"Erreur dans poly(dots[[1L]], degree, raw = raw) :
'degree' doit être inférieur au nombre de points uniques
Calls: lm ... model.frame.default -> eval -> eval -> polym -> cbind -> poly"

dans le Help de la fonction polym l'argument degree of the polynomial. Must be less than the number of unique points.
je reformule ma question : c'est quoi les points uniques finalement?



leaticia

Nombre de messages : 16
Date d'inscription : 17/05/2011

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

regression polynomial R Empty Re: regression polynomial R

Message par popotam le Mar 14 Juin 2011 - 7:29

Je connais pas la fonction polym.
Pour faire une régression polynomiale je fais quelque chose comme lm(y~x+I(x^2))

popotam

Nombre de messages : 371
Date d'inscription : 27/09/2006

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

regression polynomial R Empty Re: regression polynomial R

Message par leaticia le Mar 14 Juin 2011 - 9:04

Merci beacoup pour ta réponse, mais pour la régression polynômial multiple je dois introduire tous les degrés ?
lm(y~x+I(x²)+z+I(z²)+x*z) c'est correcte?

leaticia

Nombre de messages : 16
Date d'inscription : 17/05/2011

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

regression polynomial R Empty Re: regression polynomial R

Message par popotam le Mar 14 Juin 2011 - 9:41

Jamais essayé mais je mettrais I(x*z).

Par contre dans ton 1er post tu parles d'une variable y qui est 0-1, donc y²=y.

popotam

Nombre de messages : 371
Date d'inscription : 27/09/2006

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

regression polynomial R Empty Re: regression polynomial R

Message par droopy le Mar 14 Juin 2011 - 9:44

ça dépend ... Ici dans ton modèle tu as bien mis tout les termes, dont les termes quadratiques, mais tu as aussi mis une interaction entre x et z c'est comme si tu avais écris : y~x+z+x:z+x²+z².
Si tu supposes qu'un effet additif alors il te faut écrire :
lm(y~x+I(x^2)+z+I(z^2))
qui est équivalent à :
lm(y~poly(x,2)+poly(z,2)) # attention polygone orthogonaux
# pour avoir les mêmes coefficients qu'avec la fonction I(x^2):
lm(y~poly(x,2,raw=TRUE)+poly(z,2,raw=TRUE))
droopy
droopy

Nombre de messages : 1126
Date d'inscription : 04/09/2009

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

regression polynomial R Empty Re: regression polynomial R

Message par leaticia le Mar 14 Juin 2011 - 11:42

Alors là je suis totalement perdu
I- pour la variables binaire l'écriture lm(y~x+I(x^2)) à l'air de fonctionner, merci

II- une regression
reg<-lm(y~x+I(x^2)+z+I(z^2)+I(z*x))
reg2<-lm(y~poly(x,2,raw=TRUE)+poly(z,2,raw=TRUE)+poly(X*Z,1,raw=TRUE))
reg3<-lm(y~polym(x,z,degree=2))donne le même R² mais des coefficients totalement différents ?? No

merci encore pour votre aide

leaticia

Nombre de messages : 16
Date d'inscription : 17/05/2011

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

regression polynomial R Empty Re: regression polynomial R

Message par droopy le Mar 14 Juin 2011 - 13:12

Concrètement et en français tu cherches a établir quel modèle ?
droopy
droopy

Nombre de messages : 1126
Date d'inscription : 04/09/2009

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

regression polynomial R Empty Re: regression polynomial R

Message par leaticia le Mar 14 Juin 2011 - 13:20

Re,
Je cherche à faire une régression polynômial multiples, de la variable Y sur X etZ.
et je souhaite faire apparaître le coefficient du terme Y*X .
Merci

leaticia

Nombre de messages : 16
Date d'inscription : 17/05/2011

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

regression polynomial R Empty Re: regression polynomial R

Message par leaticia le Mar 14 Juin 2011 - 13:21

Le terme X*Z pardon !!

leaticia

Nombre de messages : 16
Date d'inscription : 17/05/2011

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

regression polynomial R Empty Re: regression polynomial R

Message par droopy le Mer 15 Juin 2011 - 7:21

Tu as des variables explicatives binaires dont tu cherches a introduire le carré dans un modèle c'est bien ça ?
droopy
droopy

Nombre de messages : 1126
Date d'inscription : 04/09/2009

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

regression polynomial R Empty Re: regression polynomial R

Message par leaticia le Mer 15 Juin 2011 - 7:57

j'avoue que je m'exprime tres mal!!!
je récapitule:
j'ai une variable "X continue" a expliqué par les variables "Y continue" et "Z binaire"
en utilisant ceci;

reg1<-lm(X~polym(Y,Z,degree=2))
voila ce que j'ai comme message d'erreur:

Erreur dans poly(dots[[i]], degree, raw = raw) :
'degree' doit être inférieur au nombre de points uniques
Calls: lm ... model.frame.default -> eval -> eval -> polym -> cbind -> poly

mais popotam "que je remercie au passage" m'a proposer une autre façon de faire la régression qui est la suivante:

reg2<-lm(X~Y+Z+I(Y^2)+I(Y*Z)) qui fonctionne apparement !!


mais si je veux faire une régression avec que des variables continues avec les methodes que tu : propose
lm(X~Y+Z+I(Y^2)+I(Y*Z))
et
lm(x~poly(y,2,raw=TRUE)+poly(z,2,raw=TRUE)) j'e n'ai pas les mêmes résultats

Merci beaucoup pr tn aide

leaticia

Nombre de messages : 16
Date d'inscription : 17/05/2011

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

regression polynomial R Empty Re: regression polynomial R

Message par droopy le Mer 15 Juin 2011 - 8:10

c'est normal que tu n'es pas les mêmes résultats entre les deux méthodes puisque ce ne sont pas les mêmes modèles
Je n'avais pas vu que z était binaire, du coup ça n'a pas de sens de mettre son terme quadratique dans le modèle d'ou le message d'erreur que tu obtiens avec polym
le modèle construit dans reg2 suppose que tu cherches l'interaction entre y et z mais seulement sur le 1er degre de cette variable est-ce que c'est ce que tu souhaites ?
Sinon pour faire reg2 tu peux aussi le taper de cette manière :
lm(x~poly(y,2,raw=T)+y*z) mais c'est un peut moins propre que la méthode de popotam. Le I(quelque chose) stipule que la variable qui est dans le crochet doit être considérée telle quelle, parce que sinon le terme ^2 à une autre signification dans la syntaxe des modèles.
droopy
droopy

Nombre de messages : 1126
Date d'inscription : 04/09/2009

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

regression polynomial R Empty Re: regression polynomial R

Message par Contenu sponsorisé


Contenu sponsorisé


Revenir en haut Aller en bas

Revenir en haut

- Sujets similaires

 
Permission de ce forum:
Vous ne pouvez pas répondre aux sujets dans ce forum