Discussion :

[davidhac]

Bonjour à tous,

Voici mon jeu de données ainsi que le code pour réaliser un diagramme à bâton :

Code :

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

1
2
3
4
5
6
7
8
9
10
 
df <- data.frame(bloc = c(1, 2, 3),
                 g  = c(20,22,18))
 
ggplot(df) +
  aes(x = bloc, y = g, fill = g, na.rm = T) +
  geom_bar(stat='identity') +
  labs(x = "Bloc", y = "G (m2/ha)", fill = "surface terrière", title = "Surface terrière par bloc")+
  geom_text(aes(label = g), position = position_stack(vjust = 0.5), size = 4, angle = 90) +
  scale_fill_gradient2(low = "yellow", high = "red", na.value = NA)

Voici le second jeu de données que je souhaiterai rajouter sur le diagramme :

Code :

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

1
2
 
cible <- data.frame(type_peuplement= c("chenaie","futaies_feuillues", "pineraies_melezins","sapinieres_pessieres"," douglasaies"), ymin = c(15,16,23,33,34), ymax = c(17,20,27,38,40))

Je souhaiterai afficher le dataframe "cible" en arrière plan sous forme de bandes horizontales de couleurs différentes et plus ou moins transparentes les 5 types de peuplements. La bande doit pouvoir faire toute la largeur du diagramme. Par exemple : pour la "chênaie", la bande irait en ordonnée de 15 à 17 (ymin, ymax) et ferait la même longueur que l'axe des abscisses. De plus j'aurai souhaité faire apparaître "type_peuplement" sur la légende.

Dans un monde parfait, l'ajout de ces bandes horizontales pourrait se faire sur tout type de dataframe ayant la même structure que"df" mais quel que soit la longueur de l'abscisse....



Je ne sais pas si ma demande est très clair et si elle est réalisable mais si vous avez d'éventuelles idées...

PS : j'ai essayé avec geom_rect(), mais je ne parvients pas à obtenir de résultat

Merci,
David

[olivier.decourt]

Bonjour David.
Pour le moment, voici une solution pas très jolie et très "à la main". Il n'y a pas de légende pour le type de peuplement, c'est à rajouter.
Ça peut sans doute s'automatiser, de manière plus ou moins élégante.

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
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
library(tidyverse)
 
df <- data.frame(bloc = c(1, 2, 3),
                 g  = c(20,22,18))
cible <- data.frame(type_peuplement= c("chenaie","futaies_feuillues", 
                                       "pineraies_melezins","sapinieres_pessieres",
                                       "douglasaies"), 
                    ymin = c(15,16,23,33,34), 
                    ymax = c(17,20,27,38,40))
 
cible2 <- pivot_wider(cible,
                      names_from=type_peuplement,
                      values_from = c(ymin, ymax))
 
df2 <- df %>% 
         left_join(cible2, by=character())
 
ggplot(df2)+
  geom_rect(aes(ymin=ymin_chenaie, ymax=ymax_chenaie),
            fill=viridis::magma(5)[[1]], alpha=0.5) +
  geom_rect(aes(ymin=ymin_futaies_feuillues, ymax=ymax_futaies_feuillues),
            fill=viridis::magma(5)[[2]], alpha=0.5) +
  geom_rect(aes(ymin=ymin_pineraies_melezins, ymax=ymax_pineraies_melezins),
            fill=viridis::magma(5)[[3]], alpha=0.5) +
  geom_rect(aes(ymin=ymin_sapinieres_pessieres, ymax=ymax_sapinieres_pessieres),
            fill=viridis::magma(5)[[4]], alpha=0.5) +
  geom_rect(aes(ymin=ymin_douglasaies, ymax=ymax_douglasaies),
            fill=viridis::magma(5)[[5]], alpha=0.5) +
  aes(x = bloc, y = g, fill = g, xmin=-Inf, xmax=Inf, label = g)  +
  geom_bar(stat='identity', orientation = "x") +
  labs(x = "Bloc", y = "G (m2/ha)", fill = "surface terrière", title = "Surface terrière par bloc")+
  geom_text(aes(y=g/2), size = 4, angle = 90) +
  scale_fill_gradient2(low = "yellow", high = "red", na.value = NA)

[statnmap]

Merci Olivier pour la piste.
Je propose une utilisation un peu plus automatisée.
Il faut ruser pour ne pas qu'il y ait deux utilisations de "fill" différente. Et pour cela, on peut utiliser un vecteur de couleurs pré-établi.
L'avantage c'est qu'on peut en mettre autant qu'on veut maintenant. Mais attention aux choses qui sont figées par l'utilisation de `cible$type_peuplement` par exemple, et à celles qui sont plus souples car dans l'aes().
Voici ma proposition:

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
19
20
21
22
23
24
25
26
library(tidyverse)
 
df <- data.frame(bloc = c(1, 2, 3),
                 g  = c(20,22,18))
 
cible <- tibble(
  type_peuplement= c("chenaie","futaies_feuillues", 
                     "pineraies_melezins","sapinieres_pessieres",
                     "douglasaies"), 
  ymin = c(15,16,23,33,34), 
  ymax = c(17,20,27,38,40),
  colour = viridis::magma(5)) %>% 
  # Arrange est necessaire pour la légende
  arrange(type_peuplement)
 
cible_colour_scale <- setNames(cible$colour, cible$type_peuplement)
 
ggplot(data = df) +  
  geom_rect(data = cible,
            aes(xmin=-Inf, xmax=Inf, ymin = ymin, ymax = ymax, colour = type_peuplement), fill = cible$colour) +
  scale_colour_manual(values = cible_colour_scale) +
  geom_col(orientation = "x", aes(x = bloc, y = g, fill = g)) +
  labs(x = "Bloc", y = "G (m2/ha)", fill = "surface terrière", title = "Surface terrière par bloc")+
  geom_text(aes(x = bloc, y = g/2, label = g), size = 4, angle = 90) +
  scale_fill_gradient2(low = "yellow", high = "red", na.value = NA) +
  guides(colour = guide_legend(override.aes = list(fill = cible$colour)))

Ce qui donne ça:

[davidhac]

Merci pour votre aide.

David