[R] Régression linéaire avec variables quantitatives et qualitatives

[petitelectron]

Bonjour à tous,

Je débute dans le language R et je souhaite faire une régression linéaire en me basant sur 3 variables l'effet d'une varaiable (lnDON, quantitative) sur une autre (lnqPCR, quantitative) en rajoutant un effet (LABO, qualitative).

data<-data[ ,c("qPCR","DON", "LABO")]
data$lnqPCR<-log(data$qPCR)
data$lnDON<-log(data$DON)
plot(lnDON~lnqPCR, data=data)
data <- data[data$lnqPCR!=0, c("lnDON", "lnqPCR", "LABO") ]
data$LABO <-as.factor(as.character(data$LABO))
champi.lm<-lm(data$lnqPCR~data$lnDON*data$LABO)
summary(champi.lm)

et j'obtiens pour les deux dernieres étapes, les messages suivants :

> champi.lm<-lm(data$lnEF1a_Fgraminearum~data$lnV70*data$LABO)
Messages d'avis :
1: In model.response(mf, "numeric") :
  l'utilisation de type="numeric" avec une réponse de type facteur sera ignorée
2: In Ops.factor(y, z$residuals) : - not meaningful for factors
> summary(champi.lm)

Call:
lm(formula = data$lnEF1a_Fgraminearum ~ data$lnV70 * data$LABO)

Residuals:
Erreur dans quantile.default(resid) : factors are not allowed

Je suis dans la mouise ...
quelqu'un sait où est mon souci ??

Merci
-----

[Dlzlogic]

Bonjour,
Je veux bien essayer de vous aider, mais je ne connais pas R.
Il faudrait que vous décriviez plus précisément votre problème de régression (indépendamment de tout langage).
Quand on parle de régression on considère généralement que l'on a N couples de valeurs x et y liés par une relation telle que on peut trouver les paramètres A et B qui vérifient au mieux la relation y = A + B.x
Or, vous parlez de 3 variables ?

[defender19]

Salut,

Tu as une mauvaise syntaxe :

Si tu veux deux variables explicatives, c'est :

modele=lm(Y~var1+var2) et non "*"

[petitelectron]

Salut ! J'ai mis "*" car c'est ma var2 qui a un effet sur var1 !! c'est pas bon ?

[A voir en vidéo sur Futura]

[Scaeva]

Bonjour,
La syntaxe est bonne, le "*" est bien autorisé. Par contre, on ne pas faire une régression avec lm() de cette manière avec une variable qualitative. Le problème est là. Qu'est ce que vous voulez faire ? Une régression par modalité de la variable qualitative ?

[petitelectron]

En fait mes soucis venaient de mon jeu de données, j'avais des logarithmes de 0 et une variable ayant des valeurs quantitatives et qualitatives ... Donc je n'aurai pas pu aller très loin ces problèmes sont résolus


Mais maintenant que j'ai mes résultats d'anova je ne sais pas très bien comment les interpréter :S


Je voulais démontrer que deux nuage de points (data$lnV70 en fonction de data$lnEF1a_Fgraminearum) avaient bien la même pente et qu'il n'y avait pas d'effet de la variable LABO. Vu le résultat de mon summary(champi.lm) mes p-value < à 0,05 donc je peux conclure qu'il n'y a pas de différence significative, qu'en gros mes valeurs ne sont pas séparées par cet effet LABO ??

En fait je ne sais pas ce que me renseigne ces résultats c'est un peux du charabia pour moi...

[Dlzlogic]

Bonjour,
Si vous me transmettez vos données, je pourrai peut-être vous éclairer.
Vous pouvez me joindre sur mon site http://www.dlzlogic.com