Plot gl and bn, subset to 28/2 and max scaler 7

Grendel/scenarios/numbers_report.Rmd

@@ -22,6 +22,8 @@ library(colorblindr)
 library(glue)
 
 date_stamp <- "2019-01-05"
+spring_cutoff <- 423
+scaler_cutoff <- 7
 ```
 
 Document compiled: `r Sys.time()`
@@ -77,19 +79,54 @@ default_numbers %>%
 left_join(numbers_per_day, ., 
           by = c("grain_dist", "day", "species"),
           suffix = c("", "_default")) %>% 
-  mutate(standarized = numbers/numbers_default) -> numbers_per_day
+  mutate(standarized = numbers/numbers_default) %>% 
+  filter(day <= spring_cutoff,
+         value <= scaler_cutoff) -> numbers_per_day
 ```
 
-```{r daily-avg-plot, echo=FALSE}
+```{r daily-avg-pf, echo=FALSE}
+sp <- "Pinkfoot"
 numbers_per_day %>% 
   filter(param == "GOOSE_BN_STARTNO_SCALER",
-         species == "pinkfoot") %>% 
+         species == tolower(sp)) %>% 
   ggplot(aes(day, standarized)) +
   geom_line(aes(color = factor(value))) + 
   scale_color_brewer(palette = "Set3", name = "Barnacle scaler") +
   hrbrthemes::theme_ipsum_rc(axis_title_size = 12) +
   labs(title = "Daily average numbers",
-       subtitle = "Pinkfeet numbers as a function of barnacle numbers",
+       subtitle = glue("{sp} numbers as a function of barnacle numbers"),
+       caption = glue("Numbers are average of {n_seasons} consecutive seasons"),
+       y = "proportion of default numbers") +
+  facet_grid(~grain_dist)
+```
+
+```{r daily-avg-bn, echo=FALSE}
+sp <- "Barnacle"
+numbers_per_day %>% 
+  filter(param == "GOOSE_BN_STARTNO_SCALER",
+         species == tolower(sp)) %>% 
+  ggplot(aes(day, standarized)) +
+  geom_line(aes(color = factor(value))) + 
+  scale_color_brewer(palette = "Set3", name = "Barnacle scaler") +
+  hrbrthemes::theme_ipsum_rc(axis_title_size = 12) +
+  labs(title = "Daily average numbers",
+       subtitle = glue("{sp} numbers as a function of barnacle numbers"),
+       caption = glue("Numbers are average of {n_seasons} consecutive seasons"),
+       y = "proportion of default numbers") +
+  facet_grid(~grain_dist)
+```
+
+```{r daily-avg-gl, echo=FALSE}
+sp <- "Greylag"
+numbers_per_day %>% 
+  filter(param == "GOOSE_BN_STARTNO_SCALER",
+         species == tolower(sp)) %>% 
+  ggplot(aes(day, standarized)) +
+  geom_line(aes(color = factor(value))) + 
+  scale_color_brewer(palette = "Set3", name = "Barnacle scaler") +
+  hrbrthemes::theme_ipsum_rc(axis_title_size = 12) +
+  labs(title = "Daily average numbers",
+       subtitle = glue("{sp} numbers as a function of barnacle numbers"),
        caption = glue("Numbers are average of {n_seasons} consecutive seasons"),
        y = "proportion of default numbers") +
   facet_grid(~grain_dist)