scenarios
library(tidyverse)
library(glue)
library(here)
library(raster)
library(knitr)
#library(bbnj)
devtools::load_all(here("../bbnj")) # devtools::install_local(force=T)
select = dplyr::select
dir_scenarios <- here("../bbnj/inst/app/www/scenarios")
dir_diffs <- here("../bbnj/inst/app/www/scenarios/diffs")
scenarios <- list(
s1 = "s01a.bio.now.mol50km",
s2a = "s02a.bio.alltime.mol50km",
s2b = "s02b.bio.future.mol50km",
s3 = "s03a.biofish.now.mol50km",
s4 = "s04a.biofish.alltime.mol50km")
diff_tbl <- function(x,...){
x %>%
mutate(
pct = sprintf("%.01f%%", pct*100)) %>%
kable(
format.args = list(big.mark=','),
align = c("r","l","r","r","r"), ...)
}
tibble(
s = names(scenarios),
scenario = unlist(scenarios)) %>%
kable()
s1 |
s01a.bio.now.mol50km |
s2a |
s02a.bio.alltime.mol50km |
s2b |
s02b.bio.future.mol50km |
s3 |
s03a.biofish.now.mol50km |
s4 |
s04a.biofish.alltime.mol50km |
For discussion: s1 - s2b
Biodiversity now (only) - Biodiversity future (only)
d <- scenarios_diff(scenarios[c("s1", "s2b")], dir_scenarios, dir_diffs)
diff_tbl(d$tbl)
-1 |
loss |
6,535 |
16,407,957 |
24.8% |
0 |
same |
14,003 |
35,158,473 |
53.2% |
1 |
gain |
5,785 |
14,524,871 |
22.0% |