Discussion :

[smilounet]

Bonjour,
Je souhaite colorer la surface de plusieurs classes de x sur une régression loess.
Actuellement, je parviens à colorer mes classes en fonction de mon plot, or, j’aimerais que la surface colorée puisse s’arrêter au niveau de la courbe loess et non jusqu’à chaque point de mon plot.

Data = enjeux
X = score
Y = nb_taxon

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
#Plot et regression Loess
 
plot(nb_taxon~score, data=enjeux, type="l")
loessMod25 <- loess(nb_taxon~score, data=enjeux, span=0.25) 
smoothed25 <- predict(loessMod25) 
 
lines(smoothed25, x=enjeux$score, col="black")
 
#Coloration des classes sur la courbe de distribution
 
i <- score >= 6 & score <= 10
polygon(c(6,score[i],10), c(0,nb_taxon[i],0), col="lightblue")
 
c <- score >= 10 & score <= 15
polygon(c(10,score[a],15), c(0,nb_taxon[a],0), col="red")
 
b <- score >= 15 & score <= 20
polygon(c(15,score[b],20), c(0,nb_taxon[b],0), col="yellow")

Avez-vous des idées permettant de colorer les classes uniquement sur la surface de la courbe loess ?

Merci par avance !

Cdlt

[marou1991]

bonjour,

il faut remplacer nb_taxon par smoothed25 ca doit marcher !

Code :

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

1
2
3
4
5
6
7
8
9
 
i <- score >= 6 & score <= 10
polygon(c(6,score[i],10), c(0,smoothed25[i],0), col="lightblue")
 
c <- score >= 10 & score <= 15
polygon(c(10,score[a],15), c(0,smoothed25[a],0), col="red")
 
b <- score >= 15 & score <= 20
polygon(c(15,score[b],20), c(0,smoothed25[b],0), col="yellow")

[smilounet]

Effectivement ça fonctionne! Merci beaucoup.. !

Connaissez vous éventuellement une fonction permettant de colorer les classes de ma courbe sous ggplot?

[marou1991]

bonjour,

pour Ggplot ca sera plutôt avec la fonction geom_polygon

Code :

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

1
2
3
4
5
6
7
8
9
 
i <- score >= 6 & score <= 10
ggplot() + geom_polygon(mapping=aes(x=c(6,score[i],10), y= c(0,smoothed25[i],0)),fill = "lightblue") 
 
c <- score >= 10 & score <= 15
ggplot() + geom_polygon(mapping=aes(x=c(10,score[a],15), y= c(0,smoothed25[a],0)),fill = "red") 
 
b <- score >= 15 & score <= 20
ggplot() + geom_polygon(mapping=aes(x=c(15,score[b],20), y= c(0,smoothed25[b],0)),fill = "yellow")