Les posteurs les plus actifs de la semaine
Aucun utilisateur |
Sujets les plus vus
regression polynomial R
3 participants
Page 1 sur 1
regression polynomial R
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?
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
Re: regression polynomial R
Je connais pas la fonction polym.
Pour faire une régression polynomiale je fais quelque chose comme lm(y~x+I(x^2))
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
Re: regression polynomial R
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?
lm(y~x+I(x²)+z+I(z²)+x*z) c'est correcte?
leaticia- Nombre de messages : 16
Date d'inscription : 17/05/2011
Re: regression polynomial R
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.
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
Re: regression polynomial R
ç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))
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- Nombre de messages : 1156
Date d'inscription : 04/09/2009
Re: regression polynomial R
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 ??
merci encore pour votre aide
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 ??
merci encore pour votre aide
leaticia- Nombre de messages : 16
Date d'inscription : 17/05/2011
Re: regression polynomial R
Concrètement et en français tu cherches a établir quel modèle ?
droopy- Nombre de messages : 1156
Date d'inscription : 04/09/2009
Re: regression polynomial R
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
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
Re: regression polynomial R
Le terme X*Z pardon !!
leaticia- Nombre de messages : 16
Date d'inscription : 17/05/2011
Re: regression polynomial R
Tu as des variables explicatives binaires dont tu cherches a introduire le carré dans un modèle c'est bien ça ?
droopy- Nombre de messages : 1156
Date d'inscription : 04/09/2009
Re: regression polynomial R
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
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
Re: regression polynomial R
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.
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- Nombre de messages : 1156
Date d'inscription : 04/09/2009
Sujets similaires
» régression multiple et régression logistique
» Régression Logistique vs Régression Linéaire
» Regression / regression médiane
» Regression 3D
» Regression log, AUC=0.54
» Régression Logistique vs Régression Linéaire
» Regression / regression médiane
» Regression 3D
» Regression log, AUC=0.54
Page 1 sur 1
Permission de ce forum:
Vous ne pouvez pas répondre aux sujets dans ce forum
|
|