Skip to content
Snippets Groups Projects
Commit c10d3991 authored by clblatte's avatar clblatte
Browse files

update rmd for checking opt data

parent 26630d56
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Main_cPouta.R
+ 7
− 7
View file @ c10d3991
...@@ -20,13 +20,13 @@ library(dplyr) ...@@ -20,13 +20,13 @@ library(dplyr)
path <- paste0(getwd()) path <- paste0(getwd())
### Running on cPouta or WIN local ### Running on cPouta or WIN local
on_cPouta <- FALSE on_cPouta <- TRUE
### Do the SQL queries for the DB (only needed first time!) ### Do the SQL queries for the DB (only needed first time!)
Process_SQL_Query <- TRUE Process_SQL_Query <- FALSE
### Export the defined columns from DB as CSV ### Export the defined columns from DB as CSV
Process_CSV_export_from_DB <- TRUE Process_CSV_export_from_DB <- FALSE
...@@ -34,7 +34,7 @@ Process_CSV_export_from_DB <- TRUE ...@@ -34,7 +34,7 @@ Process_CSV_export_from_DB <- TRUE
# Define name of genereated Optimization CSV file # Define name of genereated Optimization CSV file
# ------------------------------------------------ # ------------------------------------------------
outp_name <- "test_newCarbon" outp_name <- "rslt_RCP26_FIN_V11"
...@@ -53,8 +53,8 @@ if(on_cPouta == TRUE) { ...@@ -53,8 +53,8 @@ if(on_cPouta == TRUE) {
numCores <- detectCores() numCores <- detectCores()
numCores numCores
inputFolder <- "/media/volume/outp_rcp0/simulated" inputFolder <- "/media/volume/outp_rcp26/simulated"
outputFolder <- "/media/volume/outp_rcp0" outputFolder <- "/media/volume/outp_rcp26"
} else { } else {
# for Windows # for Windows
...@@ -219,7 +219,7 @@ columns <- paste0("id, ...@@ -219,7 +219,7 @@ columns <- paste0("id,
CARBON_STORAGE, CARBON_STORAGE,
CARBON_STORAGE_Update, CARBON_STORAGE_Update,
BM_Total, BM_total,
Carbon_soil, Carbon_soil,
CARBON_STORAGE_Ojanen, CARBON_STORAGE_Ojanen,
......
MultiForest/check_opt_data_rcp0_2.Rmd
+ 13
− 43
View file @ c10d3991
--- ---
title: "check_opt_data" title: "check_opt_data"
author: "CB" author: "CB"
date: "7/28/2020" date: "11/24/2020"
output: html_document output: html_document
--- ---
# RCP0 whole Finland - version 10 # RCP0 whole Finland - version 11
### (rslt_RCP0_FIN_V10.csv) ### (rslt_RCP0_FIN_V11.csv)
### (rslt_RCP0_FIN_V11_XY.csv)
```{r setup, include=FALSE} ```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE) knitr::opts_chunk$set(echo = TRUE)
...@@ -17,7 +18,7 @@ knitr::opts_chunk$set(echo = TRUE) ...@@ -17,7 +18,7 @@ knitr::opts_chunk$set(echo = TRUE)
on_cPouta <- TRUE on_cPouta <- TRUE
### define the of output file "outname_smooth.csv" ### define the of output file "outname_smooth.csv"
outname <- "rslt_RCP0_FIN_V10" outname <- "rslt_RCP0_FIN_V11"
### Set path to input files ### Set path to input files
...@@ -52,11 +53,11 @@ library(kableExtra) ...@@ -52,11 +53,11 @@ library(kableExtra)
```{r , echo=FALSE} ```{r , echo=FALSE}
poutaxy <- read.csv(paste0(datapath, "rslt_RCP0_FIN_V10_XY.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE) poutaxy <- read.csv(paste0(datapath, "rslt_RCP0_FIN_V11_XY.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE)
pouta0 <- read.csv(paste0(datapath, "rslt_RCP0_FIN_V10.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE) pouta0 <- read.csv(paste0(datapath, "rslt_RCP0_FIN_V11.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE)
#poutaxy <- read.csv(paste0(datapath, "rslt_test_XY.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE) #poutaxy <- read.csv(paste0(datapath, "test_newCarbon_XY.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE)
#pouta0 <- read.csv(paste0(datapath, "rslt_test.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE) #pouta0 <- read.csv(paste0(datapath, "test_newCarbon.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE)
pouta <- pouta0 %>% pouta <- pouta0 %>%
left_join(poutaxy[,c("id", "year", "regime", "standid", "CARBON_STORAGE")], by = c("id", "year", "regime")) left_join(poutaxy[,c("id", "year", "regime", "standid", "CARBON_STORAGE")], by = c("id", "year", "regime"))
...@@ -377,7 +378,7 @@ plot_carbon <- pouta %>% ...@@ -377,7 +378,7 @@ plot_carbon <- pouta %>%
boxplot(plot_carbon$CARBON_SINK, main="CARBON_SINK without the huge outlier") boxplot(plot_carbon$CARBON_SINK, main="CARBON_SINK without the huge outlier")
# before it was 1.250 million. NOW: instead periodical kg, we have yearly tonnes ( value / 5 / 1000) # before it was 1.250 million. NOW: instead periodical kg, we have yearly tonnes ( value / 5 / 1000)
carbon_sink_out <- pouta %>% filter(CARBON_SINK < -250 | CARBON_SINK > 250) %>% # source/sink of more than 200 t CO2/ha carbon_sink_out <- pouta %>% filter(CARBON_SINK < -250 | CARBON_SINK > 250) %>% # source/sink of more than 250 t CO2/ha/year
mutate(error = "CARBON_SINK") mutate(error = "CARBON_SINK")
# how many stands # how many stands
...@@ -611,7 +612,7 @@ out_all_spread[1:10,] %>% ...@@ -611,7 +612,7 @@ out_all_spread[1:10,] %>%
```{r , echo=TRUE} ```{r , echo=TRUE}
# 5 indicators are: ALL_MARKETED_MUSHROOMS, Harvested_V, i_Vm3, N_where_D_gt_40, V # 5 indicators are: ALL_MARKETED_MUSHROOMS, Harvested_V, i_Vm3, N_where_D_gt_40, V
unique(out_all[out_all$standid == 14505106,]$regime) # unique(out_all[out_all$standid == 14505106,]$regime)
# How does the input look like? # How does the input look like?
...@@ -630,7 +631,7 @@ unique(out_all[out_all$standid == 14505106,]$regime) ...@@ -630,7 +631,7 @@ unique(out_all[out_all$standid == 14505106,]$regime)
```{r , echo=TRUE} ```{r , echo=TRUE}
# 4 indicators are: CARBON_SINK, Harvested_V, i_Vm3, V # 4 indicators are: CARBON_SINK, Harvested_V, i_Vm3, V
unique(out_all[out_all$standid == 23271031,]$regime) # unique(out_all[out_all$standid == 23271031,]$regime)
# How does the input look like? # How does the input look like?
...@@ -649,7 +650,7 @@ unique(out_all[out_all$standid == 23271031,]$regime) ...@@ -649,7 +650,7 @@ unique(out_all[out_all$standid == 23271031,]$regime)
```{r , echo=TRUE} ```{r , echo=TRUE}
unique(out_all[out_all$error %in% c("i_Vm3", "ALL_MARKEDET_MUSHROOMS"),]$regime) # unique(out_all[out_all$error %in% c("i_Vm3", "ALL_MARKEDET_MUSHROOMS"),]$regime)
``` ```
...@@ -864,17 +865,6 @@ plotfkt("CAPERCAILLIE") ...@@ -864,17 +865,6 @@ plotfkt("CAPERCAILLIE")
## Indicator: HSI_CAPERCAILLIE
### according Väisänen, R. Metson (2008) (FIN only)
```{r , message=FALSE}
plotfkt("HSI_CAPERCAILLIE")
```
## Indicator: HAZEL_GROUSE ## Indicator: HAZEL_GROUSE
```{r , message=FALSE} ```{r , message=FALSE}
...@@ -885,26 +875,6 @@ plotfkt("HAZEL_GROUSE") ...@@ -885,26 +875,6 @@ plotfkt("HAZEL_GROUSE")
## Indicator: COMBINED_HSI
```{r , message=FALSE}
plotfkt("COMBINED_HSI")
```
## Indicator: weighted_mean_HSI27
```{r , message=FALSE}
plotfkt("weighted_mean_HSI27")
```
## Indicator: V_total_deadwood ## Indicator: V_total_deadwood
```{r , message=FALSE} ```{r , message=FALSE}
......
MultiForest/check_opt_data_rcp0_2_kyle.Rmd deleted 100644 → 0
+ 0
− 282
View file @ 26630d56
---
title: "check_opt_data"
author: "CB"
date: "7/28/2020"
output: html_document
---
# RCP0 whole Finland - version 7
### (rslt_RCP0_FIN_V7.csv)
```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE)
### Running on cPouta or WIN local
on_cPouta <- TRUE
### Set path to input files
if(on_cPouta == TRUE) {
# for pouta
datapath <- "/media/volume/outp_rcp0/"
path <- "/home/ubuntu/workspace/MultiForest2/"
# path <- paste0(getwd(),"/")
} else {
# for Windows
datapath <- "C:/MyTemp/r_SIMO_output/output/MF_FIN/RCP0_test/"
}
### load libraries
library(dplyr)
library(tidyr)
library(ggplot2)
library(data.table)
library(gridExtra)
library(egg)
library(kableExtra)
```
```{r , echo=FALSE}
poutaxy <- read.csv(paste0(datapath, "rslt_RCP0_FIN_V7_XY.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE)
pouta0 <- read.csv(paste0(datapath, "rslt_RCP0_FIN_V7.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE)
# poutaxy <- read.csv(paste0(datapath, "rslt_test_XY.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE)
# pouta0 <- read.csv(paste0(datapath, "rslt_test.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE)
pouta <- pouta0 %>%
dplyr::left_join(poutaxy[,c("id", "year", "regime", "standid", "CARBON_STORAGE")], by = c("id", "year", "regime"))
```
## How many stands?
```{r, echo=TRUE}
length(unique(pouta$id))
```
## Simulated regimes
```{r, echo=TRUE}
unique(pouta$regime)
```
# Check for outliers of certain indicators
Based on summary, the following are looked at due to high max-min values: </br>
- V </br>
- i_Vm3 </br>
- CARBON SINK </br>
## Indicator: V
```{r , echo=TRUE}
boxplot(pouta$V, main="Indicator: V")
V_out <- pouta %>% dplyr::filter(V > 1500) %>% # Standing volume above xx m3/ha
dplyr::mutate(error = "V")
# how many stands effected
length(unique(V_out$standid))
# ---------------
# which stands
# ---------------
unique(V_out$standid)
# which region
unique(V_out$region)
# which regime
unique(V_out$regime)
```
## Indicator: i_Vm3
```{r , echo=TRUE}
boxplot(pouta$i_Vm3, main="Indicator: i_Vm3")
i_Vm3_out <- pouta %>% filter(i_Vm3 > 50) %>% # annual volume increment above xx m3/ha
mutate(error = "i_Vm3")
# how many stands effected
length(unique(i_Vm3_out$id))
# ---------------
# which stands
# ---------------
unique(i_Vm3_out$standid)
# which region
unique(i_Vm3_out$region)
# which regime
unique(i_Vm3_out$regime)
```
### What are the high values
```{r , echo=TRUE}
i_Vm3_high <- pouta %>%
select(c("id", "standid", "year", "regime", "i_Vm3", "V","Harvested_V")) %>%
filter(id %in% unique(i_Vm3_out$id)) %>%
filter(i_Vm3 > 50 )
# first 20 rows
i_Vm3_high[1:20,] %>%
knitr::kable(caption = "FIRST 20 ROWS ONLY !!") %>%
kable_styling()
```
```{r , echo=TRUE}
## Indicator: CARBON SINK
```{r , echo=TRUE}
plot_carbon <- pouta %>%
filter(CARBON_SINK > -5000000)
boxplot(plot_carbon$CARBON_SINK, main="CARBON_SINK without the huge outlier")
carbon_sink_out <- pouta %>% filter(CARBON_SINK < -1000000 | CARBON_SINK > 1000000) %>% # source/sink of more than 1000 t CO2/ha
mutate(error = "CARBON_SINK")
# how many stands
length(unique(carbon_sink_out$id))
# ---------------
# which stands
# ---------------
unique(carbon_sink_out$standid)
# which regime
unique(carbon_sink_out$regime)
# which region
unique(carbon_sink_out$region)
```
### What are the high values
```{r , echo=TRUE}
carbononly <- pouta %>%
select(c("id", "standid", "year", "regime", "V", "Harvested_V" ,"CARBON_STORAGE", "CARBON_SINK")) %>%
filter(id %in% unique(carbon_sink_out$id))
carbononly %>%
filter(CARBON_SINK < -1000000 | CARBON_SINK > 1000000) %>%
knitr::kable(caption = "ALL ROWS !!") %>%
kable_styling()
```
### Plot some example stands
```{r , echo=TRUE}
carbonplot <- carbononly %>%
filter(regime %in% c("SA","initial_state") & id %in% 99037593 |
regime %in% c("BAU_30","initial_state") & id %in% 99034584 |
regime %in% c("CCF_2", "CCF_1", "CCF_3", "initial_state") & id %in% 99032541) %>%
mutate(deltaCarbonstorage = round(CARBON_SINK / 3.67, digits = 2))
# for scenarios BAU_30 and SA
p1 = ggplot(carbonplot[carbonplot$regime %in% c("BAU_30", "SA"),], aes(year, V)) +
geom_line() +
facet_grid(.~regime) +
theme(axis.text.x = element_text(angle = 90))
p2 = ggplot(carbonplot[carbonplot$regime %in% c("BAU_30", "SA"),], aes(year, CARBON_STORAGE)) +
geom_line() +
facet_grid(.~regime) +
theme(axis.text.x = element_text(angle = 90)) +
scale_y_continuous(labels = function(x) format(x, scientific = FALSE))
p3 = ggplot(carbonplot[carbonplot$regime %in% c("BAU_30", "SA"),], aes(year, deltaCarbonstorage)) +
geom_line() +
facet_grid(.~regime) +
theme(axis.text.x = element_text(angle = 90)) +
scale_y_continuous(labels = function(x) format(x, scientific = FALSE))
ggarrange(p1, p2, p3,
top = "BAU_30 99034584 standid: 9932854 - SA 99037593 standid: 10628286")
```
!!! HIGH Carbon sinks due to extremly high simulated Volumes !!!
```{r , echo=TRUE}
# for scenarios CC_2, _1 and _3
p4 = ggplot(carbonplot[carbonplot$regime %in% c("CCF_2", "CCF_1", "CCF_3"),], aes(year, V)) +
geom_line() +
facet_grid(.~regime) +
theme(axis.text.x = element_text(angle = 90))
p5 = ggplot(carbonplot[carbonplot$regime %in% c("CCF_2", "CCF_1", "CCF_3"),], aes(year, CARBON_STORAGE)) +
geom_line() +
facet_grid(.~regime) +
theme(axis.text.x = element_text(angle = 90)) +
scale_y_continuous(labels = function(x) format(x, scientific = FALSE))
p6 = ggplot(carbonplot[carbonplot$regime %in% c("CCF_2", "CCF_1", "CCF_3"),], aes(year, deltaCarbonstorage)) +
geom_line() +
facet_grid(.~regime) +
theme(axis.text.x = element_text(angle = 90)) +
scale_y_continuous(labels = function(x) format(x, scientific = FALSE))
ggarrange(p4, p5, p6,
top = "99032541 - standid 23271031")
```
MultiForest/check_opt_data_rcp26_2.Rmd
+ 10
− 41
View file @ c10d3991
--- ---
title: "check_opt_data" title: "check_opt_data"
author: "CB" author: "CB"
date: "7/28/2020" date: "11/24/2020"
output: html_document output: html_document
--- ---
# RCP26 whole Finland - version 10 # RCP26 whole Finland - version 11
### (rslt_RCP26_FIN_V10.csv) ### (rslt_RCP26_FIN_V11.csv)
### (rslt_RCP26_FIN_V10_XY.csv) ### (rslt_RCP26_FIN_V11_XY.csv)
```{r setup, include=FALSE} ```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE) knitr::opts_chunk$set(echo = TRUE)
...@@ -18,7 +18,7 @@ knitr::opts_chunk$set(echo = TRUE) ...@@ -18,7 +18,7 @@ knitr::opts_chunk$set(echo = TRUE)
on_cPouta <- TRUE on_cPouta <- TRUE
### define name of output file "outname_smooth.csv" ### define name of output file "outname_smooth.csv"
outname <- "rslt_RCP26_FIN_V10" outname <- "rslt_RCP26_FIN_V11"
### Set path to input files ### Set path to input files
...@@ -53,11 +53,11 @@ library(kableExtra) ...@@ -53,11 +53,11 @@ library(kableExtra)
```{r , echo=FALSE} ```{r , echo=FALSE}
poutaxy <- read.csv(paste0(datapath, "rslt_RCP26_FIN_V10_XY.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE) poutaxy <- read.csv(paste0(datapath, "rslt_RCP26_FIN_V11_XY.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE)
pouta0 <- read.csv(paste0(datapath, "rslt_RCP26_FIN_V10.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE) pouta0 <- read.csv(paste0(datapath, "rslt_RCP26_FIN_V11.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE)
#poutaxy <- read.csv(paste0(datapath, "rslt_test_XY.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE) #poutaxy <- read.csv(paste0(datapath, "test_newCarbon_XY.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE)
#pouta0 <- read.csv(paste0(datapath, "rslt_test.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE) #pouta0 <- read.csv(paste0(datapath, "test_newCarbon.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE)
pouta <- pouta0 %>% pouta <- pouta0 %>%
left_join(poutaxy[,c("id", "year", "regime", "standid", "CARBON_STORAGE")], by = c("id", "year", "regime")) left_join(poutaxy[,c("id", "year", "regime", "standid", "CARBON_STORAGE")], by = c("id", "year", "regime"))
...@@ -376,7 +376,7 @@ plot_carbon <- pouta %>% ...@@ -376,7 +376,7 @@ plot_carbon <- pouta %>%
boxplot(plot_carbon$CARBON_SINK, main="CARBON_SINK without the huge outlier") boxplot(plot_carbon$CARBON_SINK, main="CARBON_SINK without the huge outlier")
# before it was 1.250 million. NOW: instead periodical kg, we have yearly tonnes ( value / 5 / 1000) # before it was 1.250 million. NOW: instead periodical kg, we have yearly tonnes ( value / 5 / 1000)
carbon_sink_out <- pouta %>% filter(CARBON_SINK < -250 | CARBON_SINK > 250) %>% # source/sink of more than 1 t CO2/ha carbon_sink_out <- pouta %>% filter(CARBON_SINK < -250 | CARBON_SINK > 250) %>% # source/sink of more than 250 t CO2/ha/year
mutate(error = "CARBON_SINK") mutate(error = "CARBON_SINK")
# how many stands # how many stands
...@@ -733,17 +733,6 @@ plotfkt("CAPERCAILLIE") ...@@ -733,17 +733,6 @@ plotfkt("CAPERCAILLIE")
## Indicator: HSI_CAPERCAILLIE
### according Väisänen, R. Metson (2008) (FIN only)
```{r , message=FALSE}
plotfkt("HSI_CAPERCAILLIE")
```
## Indicator: HAZEL_GROUSE ## Indicator: HAZEL_GROUSE
```{r , message=FALSE} ```{r , message=FALSE}
...@@ -754,26 +743,6 @@ plotfkt("HAZEL_GROUSE") ...@@ -754,26 +743,6 @@ plotfkt("HAZEL_GROUSE")
## Indicator: COMBINED_HSI
```{r , message=FALSE}
plotfkt("COMBINED_HSI")
```
## Indicator: weighted_mean_HSI27
```{r , message=FALSE}
plotfkt("weighted_mean_HSI27")
```
## Indicator: V_total_deadwood ## Indicator: V_total_deadwood
```{r , message=FALSE} ```{r , message=FALSE}
......
MultiForest/check_opt_data_rcp45_2.Rmd
+ 10
− 41
View file @ c10d3991
--- ---
title: "check_opt_data" title: "check_opt_data"
author: "CB" author: "CB"
date: "7/28/2020" date: "11/24/2020"
output: html_document output: html_document
--- ---
# RCP45 whole Finland - version 10 # RCP45 whole Finland - version 11
### (rslt_RCP45_FIN_V10.csv) ### (rslt_RCP45_FIN_V11.csv)
### (rslt_RCP45_FIN_V10_XY.csv) ### (rslt_RCP45_FIN_V11_XY.csv)
```{r setup, include=FALSE} ```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE) knitr::opts_chunk$set(echo = TRUE)
...@@ -18,7 +18,7 @@ knitr::opts_chunk$set(echo = TRUE) ...@@ -18,7 +18,7 @@ knitr::opts_chunk$set(echo = TRUE)
on_cPouta <- TRUE on_cPouta <- TRUE
### define name of output file "outname_smooth.csv" ### define name of output file "outname_smooth.csv"
outname <- "rslt_RCP45_FIN_V10" outname <- "rslt_RCP45_FIN_V11"
### Set path to input files ### Set path to input files
...@@ -53,11 +53,11 @@ library(kableExtra) ...@@ -53,11 +53,11 @@ library(kableExtra)
```{r , echo=FALSE} ```{r , echo=FALSE}
poutaxy <- read.csv(paste0(datapath, "rslt_RCP45_FIN_V10_XY.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE) poutaxy <- read.csv(paste0(datapath, "rslt_RCP45_FIN_V11_XY.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE)
pouta0 <- read.csv(paste0(datapath, "rslt_RCP45_FIN_V10.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE) pouta0 <- read.csv(paste0(datapath, "rslt_RCP45_FIN_V11.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE)
#poutaxy <- read.csv(paste0(datapath, "rslt_test_XY.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE) #poutaxy <- read.csv(paste0(datapath, "test_newCarbon_XY.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE)
#pouta0 <- read.csv(paste0(datapath, "rslt_test.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE) #pouta0 <- read.csv(paste0(datapath, "test_newCarbon.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE)
pouta <- pouta0 %>% pouta <- pouta0 %>%
left_join(poutaxy[,c("id", "year", "regime", "standid", "CARBON_STORAGE")], by = c("id", "year", "regime")) left_join(poutaxy[,c("id", "year", "regime", "standid", "CARBON_STORAGE")], by = c("id", "year", "regime"))
...@@ -376,7 +376,7 @@ plot_carbon <- pouta %>% ...@@ -376,7 +376,7 @@ plot_carbon <- pouta %>%
boxplot(plot_carbon$CARBON_SINK, main="CARBON_SINK without the huge outlier") boxplot(plot_carbon$CARBON_SINK, main="CARBON_SINK without the huge outlier")
# before it was 1.250 million. NOW: instead periodical kg, we have yearly tonnes ( value / 5 / 1000) # before it was 1.250 million. NOW: instead periodical kg, we have yearly tonnes ( value / 5 / 1000)
carbon_sink_out <- pouta %>% filter(CARBON_SINK < -250 | CARBON_SINK > 250) %>% # source/sink of more than 1 t CO2/ha carbon_sink_out <- pouta %>% filter(CARBON_SINK < -250 | CARBON_SINK > 250) %>% # source/sink of more than 250 t CO2/ha year
mutate(error = "CARBON_SINK") mutate(error = "CARBON_SINK")
# how many stands # how many stands
...@@ -733,17 +733,6 @@ plotfkt("CAPERCAILLIE") ...@@ -733,17 +733,6 @@ plotfkt("CAPERCAILLIE")
## Indicator: HSI_CAPERCAILLIE
### according Väisänen, R. Metson (2008) (FIN only)
```{r , message=FALSE}
plotfkt("HSI_CAPERCAILLIE")
```
## Indicator: HAZEL_GROUSE ## Indicator: HAZEL_GROUSE
```{r , message=FALSE} ```{r , message=FALSE}
...@@ -754,26 +743,6 @@ plotfkt("HAZEL_GROUSE") ...@@ -754,26 +743,6 @@ plotfkt("HAZEL_GROUSE")
## Indicator: COMBINED_HSI
```{r , message=FALSE}
plotfkt("COMBINED_HSI")
```
## Indicator: weighted_mean_HSI27
```{r , message=FALSE}
plotfkt("weighted_mean_HSI27")
```
## Indicator: V_total_deadwood ## Indicator: V_total_deadwood
```{r , message=FALSE} ```{r , message=FALSE}
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment