sensitivity.Rmd

---
title: "MCT sensitivity analyses"
author: "Beni Stocker"
date: "2/15/2021"
output: html_document
---

```{r setup}
library(tidyverse)
library(segmented)

source("R/calc_cwd_lue0_v2.R")
source("R/calc_cwd_lue0_byilon.R")
source("R/calc_cwd_et0_byilon.R")
# source("R/selgmented.R")
# source("R/pscore.test.R")
```

## Change point detection

Perform change point detection in SIF ~ CWD regression.

```{r}
## get data for one longitude slice, based on SIF
calc_cwd_lue0_byilon(5000, drop_data = FALSE, dirn = "~/mct/data/df_cwd_lue0_2")
load("~/mct/data/df_cwd_lue0_2/df_cwd_lue0_5000.RData")

## get data for one longitude slice, based on EF
calc_cwd_et0_byilon(5000, drop_data = FALSE, dirn = "~/mct/data/df_cwd_et0_2/")
load("~/mct/data/df_cwd_et0_2/df_cwd_et0_5000.RData")
```

Look at one cell.
```{r}
# df$lat %>% sort() 

# caused problem: 18.175

df_data <- df %>% 
  filter(lat == 37.175) %>% 
  select(data) %>% 
  unnest(data)

df_inst <- df %>% 
  filter(lat == 37.175) %>% 
  select(data_inst) %>% 
  unnest(data_inst)

# # out <- calc_cwd_lue0(df_data, df_inst, nam_lue = "SIF", do_plot = TRUE)
# out <- calc_cwd_lue0(df_data, df_inst, nam_lue = "et", do_plot = TRUE)
# out$gg
# 
# # out <- calc_cwd_lue0(df_data, df_inst, nam_lue = "nSIF", do_plot = TRUE)
# out <- calc_cwd_lue0(df_data, df_inst, nam_lue = "fet", do_plot = TRUE)
# out$gg
```

Steps above are now implemented as
```{r}
out <- calc_cwd_lue0(df = df_data, inst = df_inst, nam_lue = "fet", do_plot = TRUE, verbose = TRUE)
out$gg + ylim(0,0.6)
out$cwd_lue0
```