Discussion :

[Fnouch]

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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

[arm3366]

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

[Fnouch]

J'ai essayé plusieurs choses avec lines

En essayant de reprendre des points de ma courbe

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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...

[A. D.]

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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

Bonne continuation !


Cordialement,


A.D.

[Fnouch]

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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

[arm3366]

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

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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.

[Fnouch]

Ça y est ! j'ai trouvé.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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