Régression linéaire sur deux nuages de points distincts

Bonjour à tous,

je cherche à rajouter deux régressions linéaires sur un plot où figure deux nuages de points différents.

je pensais à la fonction lines() ou abline() mais je ne sais pas laquelle est la meilleure (je crois que lines permet de faire plusieurs droites alors que ce n'est pas le cas de abline)
je voulais me baser sur les coefficients qui m'interessent. Seulement lorsque j'essaye avec lines ou abline aucune droites ne sort

Code:

---

1 2 3 4 5

> champi.lm<-lm(data$lnV70~data$lnEF1a_Fgraminearum+data$LABO) > coef(champi.lm) (Intercept) data$lnEF1a_Fgraminearum data$LABOMYCSA 6.1084405 0.6968371 3.9022808 > plot(lnV70~lnEF1a_Fgraminearum, data)

---

si vous avez des suggestions...

Merci

Bonjour,
Comment tu utilises les fonctions abline() ou lines() ? mets ton code si possible.

J'ai essayé plusieurs choses avec lines

En essayant de reprendre des points de ma courbe

Code:

---

1 2 3 4

plot(lnV70~lnEF1a_Fgraminearum, main="Corrélation entre la présence en Fusarium graminearum et la teneur en DON", xlab = "Présence de Fusarium graminearum (lnEF1a-Fgraminearum)", ylab = "Teneur en DON (lnV70)", col=as.factor(data$LABO), data=data) lines(x=-6.333561222, y=3.610917913, col="red")

---

ou en essayant juste avec une pente (ici celle de data$LABOMYCSA) :

Code:

---

1 2 3 4

plot(lnV70~lnEF1a_Fgraminearum, main="Corrélation entre la présence en Fusarium graminearum et la teneur en DON", xlab = "Présence de Fusarium graminearum (lnEF1a-Fgraminearum)", ylab = "Teneur en DON (lnV70)", col=as.factor(data$LABO), data=data) lines(3.9022808, col="red")

---

je n'ai pas de message d'erreur mais rien ne s'affiche sur mon plot...

Bonjour,

La fonction "lines" permet, comme son nom l'indique, de tracer une ligne/droite sur un graphe. Pour cela, il faut donc avoir, au minimum, deux points à relier l'un à l'autre. Il faut donc passer en paramètres de la fonction "lines", les coordonnées de deux points, par exemple :

Code:

---

1 2 3

plot( ... )   lines( x=c(1.2,3.5), y=c(5.2,1.3) )

---

Bonne continuation ! :)


Cordialement,


A.D.

Merci !

mais je me rends compte que la fonction lines() va être compliqué pour ce que je vais faire
il vaut mieux que j'utilise abline où je mets directement le coefficient de ma régression, comme ici :

Code:

---

1 2 3 4

plot(lnV70~lnEF1a_Fgraminearum, main="Corrélation entre la présence en Fusarium graminearum et la teneur en DON", xlab = "Présence de Fusarium graminearum (lnEF1a-Fgraminearum)", ylab = "Teneur en DON (lnV70)", col=as.factor(data$LABO), data=data) abline(a=coef(champi.lm)[1], b=coef(champi.lm)[2])

---

j'obtiens une droite de régression pour un de mes nuages de points mais je n'arrive pas à en faire plusieurs

Bonjour,
Est ce que tu as essayé quelque chose comme ça:

Code:

---

1 2	plot(Y~X, main="Scatterplot", ylab= "Y.nom",xlab="X.nom") abline(lm(Y~X), col="blue")

---

Et partir de la, tu pourra ajouter autant de "abline " que tu souhaite.
A defaut d'un jeu de données et les erreurs obtenus, ça sera toujours plus difficile d'y voir plus clair.

Ça y est ! j'ai trouvé.

Code:

---

1 2 3 4

plot(lnV70~lnEF1a_Fgraminearum, data)   abline(a=coef(champi.lm)[1], b=coef(champi.lm)[2]) abline(a=coef(champi.lm)[1] + coef(champi.lm)[3] ,b=coef(champi.lm)[2],col="red")

---

Je reprends entre crochet les coefficients de ma régression champi.lm ([1] pour l'intercept, [2] pour data$lnEF1a_Fgraminearum et [3]data$LABOMYCSA)

En espérant que ça puisse aider quelqu'un