Names of the restructered GPKGs

gpkg <- unique(rslt$gpkg)
print(gpkg)
##  [1] "MV_Hartola"   "MV_Kitee"     "MV_Korsnas"   "MV_Parikkala"
##  [5] "MV_Pori"      "MV_Pyhtaa"    "MV_Raasepori" "MV_Simo"     
##  [9] "MV_Vaala"     "MV_Voyri"

Simulated regimes

##  [1] "SA"         "CCF_1"      "CCF_2"      "CCF_3"      "CCF_4"     
##  [6] "BAUwT"      "BAUwT_GTR"  "BAUwT_m5"   "BAUwT_5"    "BAUwT_10"  
## [11] "BAUwT_15"   "BAUwT_30"   "BAUwoT_m20" "BAU_m5"     "BAU"       
## [16] "BAUwGTR"    "BAUwoT"     "BAU_5"      "BAU_10"     "BAUwoT_10" 
## [21] "BAU_15"     "BAU_30"

Number of simulated stands per GPKG and their size

stands <- rslt %>% 
  group_by(gpkg) %>% 
  summarise(simulated_stands = n_distinct(id),
            min_stand_size = min(AREA),
            max_stand_size = max(AREA),
            mean_size = mean(AREA))

kable(stands) %>% kable_styling()
gpkg simulated_stands min_stand_size max_stand_size mean_size
MV_Hartola 5 0.771 2.467 1.5679802
MV_Kitee 5 0.496 2.404 1.1994000
MV_Korsnas 5 0.289 4.861 2.0297640
MV_Parikkala 5 0.152 1.200 0.8710000
MV_Pori 5 0.112 0.637 0.3596000
MV_Pyhtaa 5 0.498 2.448 1.2284896
MV_Raasepori 5 0.068 1.784 0.2184878
MV_Simo 5 0.141 4.792 2.6578507
MV_Vaala 5 0.454 4.682 3.0449601
MV_Voyri 5 0.355 1.238 0.6504268

Development of average stand volume (m3/ha) for certain regimes

meanV <- rslt[rslt$regime %in% c("BAU", "SA", "CCF_1") , ] %>% 
   group_by(year, regime, gpkg) %>% 
   summarise(meanV = mean(V)) %>% 
   ggplot(aes(year, meanV)) + 
   geom_line( aes(color = regime)) +
   facet_wrap(. ~gpkg)
 plot(meanV)

Average cash flow (Euro/ha) for certain regimes

Cash flow = The sum of all revenues and costs for a specific forest stand

meanCash <- rslt[rslt$regime %in% c("BAU", "SA", "CCF_1") , ] %>% 
   group_by(year, regime, gpkg) %>% 
   mutate(cash_flow = ifelse(is.na(cash_flow), 0, cash_flow)) %>% 
   summarise(meanCash = mean(cash_flow), na.rm = TRUE) %>% 
   ggplot(aes(year, meanCash)) + 
   geom_line( aes(color = regime)) +
   facet_wrap(. ~gpkg)
 plot(meanCash)