# INBO CODING CLUB
30 September 2021
Welcome!
## Share your code snippet
If you want to share your code snippet, copy paste your snippet within a section of three backticks (```):
As an **example**:
```
library(tidyverse)
```
(*you can copy paste this example and add your code further down*)
## Yellow sticky notes
Sticky notes are back! But still, some of us are remotely connected. Add you name here below and add a star for each challenge you could solve.
## Participants
Name | Challenges
--- | ---
Damiano Oldoni |
Emma Cartuyvels|**
Lynn Pallemaerts |**
Dirk Maes |
Tom De Dobbelaer |**
Jasmijn Hillaert |**
Cécile Herr |***
Adriaan Seynaeve |
Matthieu Chastel |
Marijke Thoonen |
Geert Spanoghe
Suzanna Lettens |**
Jolien Goossens |**
Maud Raman |
## Challenge 1
Cécile:
```
mypalette <- c("#391B15","#454961","#3F825A","#DE9250")
# Or from: https://sashamaps.net/docs/resources/20-colors/
mypalette <- c("black", "#4363d8", "#f58231", "#a9a9a9")
#' 2. How to use it for user case 1?
usercase1 <-
ggplot(area_biotopes_pararge,
aes(x = year,
y = meanArea,
color = biotope)) +
geom_point() +
geom_smooth() +
scale_color_manual(values = mypalette) +
labs(title = "Distribution of Pararge aegeria per biotope and region",
y = "Area (%)",
color = "BIOTOPE"
) +
facet_wrap(~ region)
usercase1
# 3. How to use it for user case 3?
usercase3 <-
ggplot(area_biotopes_2000,
aes(x = region, y = biotope)) +
geom_tile(aes(fill = meanArea)) +
scale_fill_gradientn(colours = mypalette) +
labs(title = "Average covered area (%) in 2000",
fill = "Area (%)") +
facet_wrap(~ species)
usercase3
# 4. How to invert the color scale of `usercase3`?
usercase3_inv <-
ggplot(area_biotopes_2000,
aes(x = region, y = biotope)) +
geom_tile(aes(fill = meanArea)) +
scale_fill_gradientn(colours = rev(mypalette)) +
labs(title = "Average covered area (%) in 2000",
fill = "Area (%)") +
facet_wrap(~ species)
usercase3_inv
```
Lynn
```
pal <- c("#554018","#7E5728","#AC6D3D","#DD8358")
usercase1 +
scale_color_manual(values=pal)
usercase3 +
scale_fill_continuous(type="viridis")
usercase3 +
scale_fill_continuous(type="viridis", direction=-1)
```
Emma:
```
#' 2. How to use it for user case 1?
usercase1 +
scale_color_manual(values = c("#56A16A","#7A9346","#998235","#AD6E3D"))
# 3. How to use it for user case 3?
usercase3 +
scale_fill_gradientn(colors = c("#56A16A","#7A9346","#998235","#AD6E3D"))
# 4. How to invert the color scale of `usercase3`?
usercase3 +
scale_fill_gradientn(colors = c("#AD6E3D","#998235","#7A9346","#56A16A"))
usercase3 +
scale_fill_gradientn(colors = rev(c("#56A16A","#7A9346","#998235","#AD6E3D")))
```
The package "ggsci" has a lot of predefined palettes, including colourblind-friendly ones
## Challenge 2
Cécile
```
#' 1. Display `usercase1` and `usercase2` next to each other
usercase1 + usercase2
#' 2. Display `usercase2` below `usercase1`
usercase1 / usercase2
#' 3. Display `usercase1` and `usercase2`in first row and `usercase3` and
#' `usercase4` in the second row creating a 2x2 plot grid
# this is the defaut behaviour for 4 plots, so this works:
usercase1 + usercase2 +usercase3 + usercase4
#' 4. `usercase1` and `usercase2` share same legend. How to avoid duplicating it in 1 and 2?
(usercase1 + usercase2 + plot_layout(guides = 'collect')) / (usercase3 + usercase4 + plot_layout(guides = 'keep'))
#' 5. What if you try to avoid duplicating the legend of `usercase3` and
#' `usercase4`?
# I will just take the maximum value and
# use it as max limit in both plots:
maxlim <- max(
area_biotopes_2000 %>% summarize(maxarea = max(meanArea)),
area_biotopes_2018 %>% summarize(maxarea = max(meanArea)))
usercase3 <- usercase3 +
scale_fill_continuous(limits = c(0, maxlim))
usercase4 <- usercase4 +
scale_fill_continuous(limits = c(0, maxlim))
(usercase1 + usercase2 ) / (usercase3 + usercase4) + plot_layout(guides = 'collect')
```
Matthieu
#' 1. Display `usercase1` and `usercase1` next to each other
usercase1 +
usercase2
#par(mfrow=c(1,2)) doesn't work?!
Lynn
```
#' 1. Display `usercase1` and `usercase2` next to each other
usercase1 + usercase2
#' 2. Display `usercase2` below `usercase1`
usercase1 / usercase2
#' 3. Display `usercase1` and `usercase2`in first row and `usercase3` and
#' `usercase4` in the second row creating a 2x2 plot grid
(usercase1 + usercase2) / (usercase3 + usercase4)
#' 4. `usercase1` and `usercase2` share same legend. How to avoid duplicating
#' it in 1 and 2?
(usercase1 + usercase2) +
plot_layout(guides="collect")
#' 5. What if you try to avoid duplicating the legend of `usercase3` and
#' `usercase4`? Can you understand why? How to solve it? Hint: _scale_ your
#' plots, check ggplot cheatsheet!
usercase3 + usercase4 +
plot_layout(guides="collect")
(usercase3 + scale_fill_continuous(limits=c(0,70))) +
(usercase4 + scale_fill_continuous(limits=c(0,70))) +
plot_layout(guides="collect")
```
Emma:
```
#' 1. Display `usercase1` and `usercase2` next to each other
usercase1 + usercase2
#' 2. Display `usercase2` below `usercase1`
usercase1/usercase2
#' 3. Display `usercase1` and `usercase2`in first row and `usercase3` and
#' `usercase4` in the second row creating a 2x2 plot grid
(usercase1|usercase2)/(usercase3|usercase4)
#' 4. `usercase1` and `usercase2` share same legend. How to avoid duplicating
#' it in 1 and 2?
usercase1/usercase2 +
plot_layout(guides = 'collect')
#' 5. What if you try to avoid duplicating the legend of `usercase3` and
#' `usercase4`? Can you understand why? How to solve it? Hint: _scale_ your
#' plots, check ggplot cheatsheet!
usercase3 + scale_fill_continuous(limits = c(0,75)) +
usercase4 + scale_fill_continuous(limits = c(0,75)) +
plot_layout(guides = 'collect')
```
## Challenge 3
Matthieu
#' 2. Remove the title and the legend of `usercase3`
usercase3 <-
ggplot(area_biotopes_2000,
aes(x = region, y = biotope)) +
geom_tile(aes(fill = meanArea)) +
#'' labs(title = "Average covered area (%) in 2000",
#'' fill = "Area (%)") +
theme(legend.position = "none")+
scale_fill_continuous(type="viridis", direction=-1)+
facet_wrap(~ species)
usercase3
## Bonus challenge
Sign in with Wallet
Connect another wallet