Rmd: maximum H_dom per watershed and its regime

quick_summary.Rmd

@@ -25,7 +25,7 @@ path <- paste0(getwd(),"/")
 #               "CC85"
 
 
-sim_variant = "without"
+sim_variant = "CC85"
 
 
 library(tidyr)
@@ -117,6 +117,37 @@ meanV_CCF <- rslt[rslt$regime %in% c("CCF_1", "CCF_2", "CCF_3", "CCF_4"), ] %>%
 
 ```
 
+### Development of average h_dom for certain regimes
+```{r}
+meanH_dom <- rslt[rslt$regime %in% c("BAU", "SA", "CCF_2", "BAUwGTR", "SA_DWextract") , ] %>% 
+   group_by(year, regime, gpkg) %>% 
+   summarise(meanH_dom = mean(H_dom, na.rm = TRUE )) %>% 
+   ggplot(aes(year, meanH_dom)) + 
+   geom_line( aes(color = regime)) +
+   theme(axis.text.x = element_text(angle = 90)) +
+   scale_x_continuous(breaks=c(2016, 2026, 2036, 2046, 2056, 2066, 2076, 2086, 2096, 2106)) +
+   facet_wrap(. ~gpkg)
+ plot(meanH_dom)
+```
+
+
+ 
+### What is the maximum H_dom per watershed, and which regime is causing it
+
+```{r, warning=FALSE, message=FALSE}
+require(data.table)
+dt <- data.table(rslt)
+max <- dt[ , max(H_dom, na.rm = TRUE ), by = gpkg]
+max <- max %>% rename(H_dom = V1 )
+
+maxH_dom <- rslt %>% 
+   semi_join(max, by = c("gpkg","H_dom")) %>% 
+   select(gpkg, id, year, regime, H_dom)
+
+kable(maxH_dom) %>% kable_styling()
+ 
+```
+
 
 ### Development of average V_total_deadwood (m3/ha) for certain regimes
 ```{r}