Discussion :

[habasque]

Bonjour,

Je souhaiterais représentant sur une carte des valeurs d'un champ Z ayant des coordonnées de type X,Y . Pour le champ Z, je cherche une méthode permettant de proposer une échelle de couleurs (l'équivalent de scatterplot et colorbar sous Matlab).
Pour le moment j'utilise le code suivant :

Code :

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

1
2
3
4
5
plot(dateset$Longitude,dataset$Latitude, xlab = 'Longitude', ylab = 'Latitude',xlim=c(-82,-70), ylim =c(-20,-5))
      par(new=TRUE)
      colors=cut(dataset$Depth, breaks=c(0, 10, 20, 30), labels=c("red", "green", "blue"))
      points(dataset$Longitude,dataset$Latitude,col=colors,pch=16)
      legend("topright", legend=c(0, 10, 20, 30),lwd=2,col=colors)

Quelqu'un aurait une idée, je préférerais une palette de couleur qui se créé automatiquement suivant les valeurs de Z ?
Merci d'avance,
Cdt,

[Sengar]

Salut,

Je connais mal lattice mais avec le package ggplot2 c'est très simple.
Exemple:

Code :

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

1
2
3
4
5
6
x<-1:10
y<-1:10
z<-c(1:3,1:3,1:4)
 
 
qplot(x,y,col=z,geom="point")

Et si tu veux changer le gradient de couleurs:

Code :

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

1
2
qplot(x,y,col=z,geom="point")+
scale_colour_gradient(low="red", high="blue")

Ce qui signifie que les couleurs vont passer graduellement du rouge (pour les petites valeurs de y ) au bleu (pour les plus grosses).

Il y a aussi d'autres types de gradients qui passent par plusieurs couleurs prédéfinies (pour ca voir l'aide de scale_colour_gradient) si tu es du genre perfectionniste.




Edit: Autre idée peut être plus pertinente que d'utiliser des couleurs: faire varier la taille des points proportionnellement à la valeur de z:

Code :

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

qplot(x,y,size=z,geom="point")

[habasque]

Merci c'est tout à fait ce que je recherchais ! A priori c'est un package de base !!

Voici le code :

Code :

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

1
2
3
4
5
x11(width=12,height=15,pointsize=12) 
p <-qplot(Longitude, Latitude, data=dataset, col=Depth, facets=Ano~Mes,xlim=c(-82,-70),ylim =c(-20,-5),xlab = 'Longitude', ylab = 'Latitude',main="Profundidad media de jurel por ano y mes - dia")+
  scale_colour_gradient(low="red", high="blue")
p + theme_bw()
lines(filename_shoreline)

Maintenant je souhaiterais superposer une ligne à chaque sous-plot mais en appelant 'lines' j'ai l'erreur suivante :

plot.new n'a pas encore été appelé

Faut il intégrer d'une manière ou d'une autre l'appel à lines dans la fonction qplot elle-même ?

Merci d'avance,

[Sengar]

Le mieux quand tu as plusieurs layers à superposer est de travailler avec ggplot() au lieu de qplot() qui est la version allégée: q pour quick.

Il faut juste comprendre que chaque couche d'un layer travaille à partir d'un data frame, et utilise les variables de ce data frame dès que tu fais un aes().

Pour ton problème, il faut aussi duppliquer ta courbe à appliquer à chacun des facet (si j'ai bien compris c'est la même pour chaque facet).
Tu la duppliques donc autant de fois que tu as de facet.

Je te donne le code ci dessous pour t'inspirer.

Ggplot2 demande un coup d'entrée pour tout comprendre, mais le jeu en vaut la chandelle (quand tu maitrises, tu fais des super graphes en deux lignes de code ^^).

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
x<-1:10
y<-1:10
z<-c(1:3,1:3,1:4)
t<-c(rep(1,5),rep(2,5))
data<-data.frame(x,y,z,t) #le data frame des données
 
 
#Ta courbe que je dupplique pour chaque facet: ici le facet depend de la variable t
x2<-1:10
y2<-sqrt(x2)
data2<-data.frame(x2,y2)
data2<-cbind.data.frame(data2,rep(t,each=nrow(data2))) # le data frame de la courbe
 
ggplot(data=data,aes(x=x,y=y))+geom_point(aes(col=z))+facet_grid(t~.)+
  geom_line(data=data2,aes(x=x2,y=y2))

[habasque]

Merci tout fonctionne !

Code :

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

1
2
3
 
p <- ggplot(data=dataset_dia,aes(x=Lon_S,y=Lat_S))+geom_point(aes(col=Depth_mean),xlab = 'Longitude', ylab = 'Latitude')+facet_grid(Ano~Mes)
p + geom_line(data=filename_shoreline,aes(x=Lon_S,y=Lat_S))