Last updated: 2021-08-02

Checks: 7 0

Knit directory: emlr_obs_v_XXX/

This reproducible R Markdown analysis was created with workflowr (version 1.6.2). The Checks tab describes the reproducibility checks that were applied when the results were created. The Past versions tab lists the development history.

R Markdown file: up-to-date

Great! Since the R Markdown file has been committed to the Git repository, you know the exact version of the code that produced these results.

Environment: empty

Great job! The global environment was empty. Objects defined in the global environment can affect the analysis in your R Markdown file in unknown ways. For reproduciblity it’s best to always run the code in an empty environment.

Seed: set.seed(20200707)

The command set.seed(20200707) was run prior to running the code in the R Markdown file. Setting a seed ensures that any results that rely on randomness, e.g. subsampling or permutations, are reproducible.

Session information: recorded

Great job! Recording the operating system, R version, and package versions is critical for reproducibility.

Cache: none

Nice! There were no cached chunks for this analysis, so you can be confident that you successfully produced the results during this run.

File paths: relative

Great job! Using relative paths to the files within your workflowr project makes it easier to run your code on other machines.

Repository version: defe32c

Great! You are using Git for version control. Tracking code development and connecting the code version to the results is critical for reproducibility.

The results in this page were generated with repository version defe32c. See the Past versions tab to see a history of the changes made to the R Markdown and HTML files.

Note that you need to be careful to ensure that all relevant files for the analysis have been committed to Git prior to generating the results (you can use wflow_publish or wflow_git_commit). workflowr only checks the R Markdown file, but you know if there are other scripts or data files that it depends on. Below is the status of the Git repository when the results were generated: 


Ignored files:
    Ignored:    .Rhistory
    Ignored:    .Rproj.user/

Unstaged changes:
    Modified:   code/Workflowr_project_managment.R
    Modified:   data/auxillary/params_local.rds

Note that any generated files, e.g. HTML, png, CSS, etc., are not included in this status report because it is ok for generated content to have uncommitted changes.

These are the previous versions of the repository in which changes were made to the R Markdown (analysis/eMLR_model_fitting.Rmd) and HTML (docs/eMLR_model_fitting.html) files. If you’ve configured a remote Git repository (see ?wflow_git_remote), click on the hyperlinks in the table below to view the files as they were in that past version.

File Version Author Date Message
html 127b801 jens-daniel-mueller 2021-07-24 Build site.
html 1f65ef1 jens-daniel-mueller 2021-07-23 Build site.
html 912d90e jens-daniel-mueller 2021-07-23 Build site.
html 2477316 jens-daniel-mueller 2021-07-23 rebuild: surface dcant mapping seperate
html c9ccc00 jens-daniel-mueller 2021-07-22 Build site.
html 426b2df jens-daniel-mueller 2021-07-21 Build site.
html 971ce87 jens-daniel-mueller 2021-07-13 Build site.
Rmd 89d9fcb jens-daniel-mueller 2021-07-13 complete revision
html c18a8b1 jens-daniel-mueller 2021-07-09 Build site.
Rmd 795a3b0 jens-daniel-mueller 2021-07-09 complete revision
html 910d64a jens-daniel-mueller 2021-07-02 Build site.
html 1cbf907 jens-daniel-mueller 2021-07-02 Build site.
html 955f5dc jens-daniel-mueller 2021-06-29 Build site.
html 9480f94 jens-daniel-mueller 2021-06-29 Build site.
Rmd af271b5 jens-daniel-mueller 2021-06-29 plot spatial residual ditribution
html 147cb07 jens-daniel-mueller 2021-06-29 Build site.
Rmd 80ae560 jens-daniel-mueller 2021-06-29 plot spatial residual ditribution
html d054c39 jens-daniel-mueller 2021-06-24 Build site.
html 4fcd066 jens-daniel-mueller 2021-06-24 Build site.
html 3b6387f jens-daniel-mueller 2021-06-18 Build site.
html 3052a6c jens-daniel-mueller 2021-06-15 Build site.
html 3cb64bd jens-daniel-mueller 2021-06-15 Build site.
html c6b3da6 jens-daniel-mueller 2021-06-14 Build site.
html 48c73fc jens-daniel-mueller 2021-06-14 Build site.
html 439ee80 jens-daniel-mueller 2021-06-11 Build site.
html 33ffcab jens-daniel-mueller 2021-06-10 Build site.
html 7e1f407 jens-daniel-mueller 2021-06-10 Build site.
html 2cbe18c jens-daniel-mueller 2021-06-10 added zonal mean section control plots
html 69c79d0 jens-daniel-mueller 2021-06-08 Build site.
html 1772903 jens-daniel-mueller 2021-06-07 Build site.
html 594ed9a jens-daniel-mueller 2021-06-04 Build site.
html db7df0e jens-daniel-mueller 2021-06-04 rebuild without overlapping eras
html 2edc791 jens-daniel-mueller 2021-06-04 Build site.
html 207339d jens-daniel-mueller 2021-06-03 Build site.
html 315710b jens-daniel-mueller 2021-06-03 include anomalous changes
html be90356 jens-daniel-mueller 2021-06-02 Build site.
html d37a85d jens-daniel-mueller 2021-05-31 Build site.
html 4b7a5ee jens-daniel-mueller 2021-05-28 Build site.
Rmd a61014f jens-daniel-mueller 2021-05-28 test run with cstar target variable
html 12b455a jens-daniel-mueller 2021-05-27 Build site.
Rmd ee13efb jens-daniel-mueller 2021-05-27 optional source of local params fully implemented
html 8c736a6 jens-daniel-mueller 2021-05-27 Build site.
Rmd dc8e4e1 jens-daniel-mueller 2021-05-27 optional source of local params implemented
html 25bd183 jens-daniel-mueller 2021-05-26 Build site.
Rmd e2b3434 jens-daniel-mueller 2021-05-26 generate seperate output
html b79cb2d jens-daniel-mueller 2021-05-20 Build site.
Rmd 8adf952 jens-daniel-mueller 2021-05-20 parallelisation tested, not succesful
html 62bd574 jens-daniel-mueller 2021-05-20 Build site.
html 7c56c39 jens-daniel-mueller 2021-05-19 Build site.
html 0de759e jens-daniel-mueller 2021-05-13 Build site.
html 52e7583 jens-daniel-mueller 2021-05-12 Build site.
html 969e631 jens-daniel-mueller 2021-05-12 Build site.
html d2a83bc jens-daniel-mueller 2021-04-16 Build site.
html c0a47df jens-daniel-mueller 2021-04-16 Build site.
html 50290e8 jens-daniel-mueller 2021-04-16 overlapping eras
html a00ec94 jens-daniel-mueller 2021-04-16 Build site.
html b6fe355 jens-daniel-mueller 2021-04-16 Build site.
html 81b7c6d jens-daniel-mueller 2021-04-16 Build site.
html ddec5b7 jens-daniel-mueller 2021-04-15 Build site.
html 29edae5 jens-daniel-mueller 2021-04-14 Build site.
html 099d566 jens-daniel-mueller 2021-04-14 Build site.
html bb44686 jens-daniel-mueller 2021-04-14 Build site.
html bf40480 jens-daniel-mueller 2021-04-13 Build site.
html 9f31fe3 jens-daniel-mueller 2021-04-13 Build site.
html c98be84 jens-daniel-mueller 2021-04-13 Build site.
html 338dd3c jens-daniel-mueller 2021-04-09 Build site.
html a79ca2c jens-daniel-mueller 2021-04-09 included model data
html 7831fd3 jens-daniel-mueller 2021-04-09 Build site.
html 20b26d4 jens-daniel-mueller 2021-04-08 Build site.
html c0895f8 jens-daniel-mueller 2021-04-07 Build site.
html 156d5b7 jens-daniel-mueller 2021-04-07 Build site.
html eb827c9 jens-daniel-mueller 2021-04-07 Build site.
html d3c3fde jens-daniel-mueller 2021-03-26 Build site.
Rmd 88fa519 jens-daniel-mueller 2021-03-26 included model data
html 594d503 jens-daniel-mueller 2021-03-26 Build site.
Rmd e4ef0f6 jens-daniel-mueller 2021-03-26 included model data
html 06bfaf7 jens-daniel-mueller 2021-03-26 Build site.
Rmd 6d2102a jens-daniel-mueller 2021-03-26 included model data
html 2eb881a jens-daniel-mueller 2021-03-26 Build site.
Rmd a9e976a jens-daniel-mueller 2021-03-26 included model data
html 03b6009 jens-daniel-mueller 2021-03-23 removed before copying template
html de6256b jens-daniel-mueller 2021-03-23 Build site.
html f155edd jens-daniel-mueller 2021-03-23 Build site.
html 380d215 jens-daniel-mueller 2021-03-21 Build site.
html 33b385b jens-daniel-mueller 2021-03-20 Build site.
Rmd 6c5618b jens-daniel-mueller 2021-03-20 rerun with vif threshold
html 108d6b2 jens-daniel-mueller 2021-03-20 Build site.
Rmd 2b99ba0 jens-daniel-mueller 2021-03-20 included vif removel option
html 330dcd0 jens-daniel-mueller 2021-03-20 Build site.
Rmd 4df8990 jens-daniel-mueller 2021-03-20 added vif calculation and plots
html 83a13de jens-daniel-mueller 2021-03-20 Build site.
html cf98c6d jens-daniel-mueller 2021-03-16 Build site.
html a1d52ff jens-daniel-mueller 2021-03-15 Build site.
html 0bade3b jens-daniel-mueller 2021-03-15 Build site.
html 27c1f4b jens-daniel-mueller 2021-03-14 Build site.
html af75ebf jens-daniel-mueller 2021-03-14 Build site.
html 5017709 jens-daniel-mueller 2021-03-11 Build site.
html 585b07f jens-daniel-mueller 2021-03-11 Build site.
html 6482ed7 jens-daniel-mueller 2021-03-11 Build site.
html 85a5ed2 jens-daniel-mueller 2021-03-10 Build site.
html 00688a1 jens-daniel-mueller 2021-03-05 Build site.
html 6c0bec6 jens-daniel-mueller 2021-03-05 Build site.
html 3c2ec33 jens-daniel-mueller 2021-03-05 Build site.
html af70b94 jens-daniel-mueller 2021-03-04 Build site.
Rmd c9cf1fd jens-daniel-mueller 2021-03-04 rebuild with NA in Cant replaced by 0
html 27ae473 jens-daniel-mueller 2021-02-24 Build site.
Rmd 7f77d91 jens-daniel-mueller 2021-02-24 removed log10 color scale
html fec3558 jens-daniel-mueller 2021-02-24 Build site.
Rmd 9ebedac jens-daniel-mueller 2021-02-24 latitudinal residual plots
html 4bc00ea jens-daniel-mueller 2021-02-24 Build site.
Rmd de11bfe jens-daniel-mueller 2021-02-24 clean up purrr approach and residual plots
html 42eca5d jens-daniel-mueller 2021-02-24 Build site.
Rmd 06a2f3b jens-daniel-mueller 2021-02-24 purrr residual plots by basin
html a1ba577 jens-daniel-mueller 2021-02-24 Build site.
Rmd 9ae7d87 jens-daniel-mueller 2021-02-24 loop residual plots by basin
html 071743d jens-daniel-mueller 2021-02-24 Build site.
Rmd c45672c jens-daniel-mueller 2021-02-24 added residual plots
html ac1a836 jens-daniel-mueller 2021-02-24 Build site.
Rmd 5f655e0 jens-daniel-mueller 2021-02-24 added plots back to after switching to map aproach
html b03fbd3 jens-daniel-mueller 2021-02-24 Build site.
Rmd c69736b jens-daniel-mueller 2021-02-24 added plots back to after switching to map aproach
html 86406d5 jens-daniel-mueller 2021-02-24 Build site.
Rmd 1b3c171 jens-daniel-mueller 2021-02-24 introduced purrr::map to model fitting, rebuild all
html 3d3b4cc jens-daniel-mueller 2021-02-23 Build site.
Rmd cbfc388 jens-daniel-mueller 2021-02-23 introduced purrr::map to model fitting
html 7b672f7 jens-daniel-mueller 2021-01-11 Build site.
html 33ba23c jens-daniel-mueller 2021-01-07 Build site.
Rmd 0ad30ba jens-daniel-mueller 2021-01-07 removed GLODAP gamma filter, target variable mapped by eras+era
html 318609d jens-daniel-mueller 2020-12-23 adapted more variable predictor selection
html 9d0b2d0 jens-daniel-mueller 2020-12-23 Build site.
html 0aa2b50 jens-daniel-mueller 2020-12-23 remove html before duplication
html 39113c3 jens-daniel-mueller 2020-12-23 Build site.
Rmd bef9220 jens-daniel-mueller 2020-12-23 rebuild before sensitivity test
html 2886da0 jens-daniel-mueller 2020-12-19 Build site.
html 02f0ee9 jens-daniel-mueller 2020-12-18 cleaned up for copying template
html 965dba3 jens-daniel-mueller 2020-12-18 Build site.
html 5d452fe jens-daniel-mueller 2020-12-18 Build site.
Rmd ca65bf5 jens-daniel-mueller 2020-12-18 rebuild after final cleaning
html 7bcb4eb jens-daniel-mueller 2020-12-18 Build site.
html d397028 jens-daniel-mueller 2020-12-18 Build site.
Rmd 7e1b1c0 jens-daniel-mueller 2020-12-18 rebuild without na predictors
html 7131186 jens-daniel-mueller 2020-12-17 Build site.
Rmd 737d904 jens-daniel-mueller 2020-12-17 rebuild without na predictors
html 22b07fb jens-daniel-mueller 2020-12-17 Build site.
html a84ff3c jens-daniel-mueller 2020-12-17 Build site.
Rmd 40369db jens-daniel-mueller 2020-12-17 model selection criterion added
html 5b48ef5 jens-daniel-mueller 2020-12-17 Build site.
Rmd e6ed2bc jens-daniel-mueller 2020-12-17 plotted model results
html f3a708f jens-daniel-mueller 2020-12-17 Build site.
Rmd 7c8ace9 jens-daniel-mueller 2020-12-17 new MLR fitting routine, rmse corrected
html e4ca289 jens-daniel-mueller 2020-12-16 Build site.
Rmd 3d5a3e2 jens-daniel-mueller 2020-12-16 new MLR fitting routine
html 158fe26 jens-daniel-mueller 2020-12-15 Build site.
html 7a9a4cb jens-daniel-mueller 2020-12-15 Build site.
html 61b263c jens-daniel-mueller 2020-12-15 Build site.
html 4d612dd jens-daniel-mueller 2020-12-15 Build site.
html e91cebd jens-daniel-mueller 2020-12-15 Build site.
Rmd d7992c4 jens-daniel-mueller 2020-12-15 eMLR target variable selection
html 953caf3 jens-daniel-mueller 2020-12-15 Build site.
html 42daf5c jens-daniel-mueller 2020-12-14 Build site.
Rmd 923aa7f jens-daniel-mueller 2020-12-14 rebuild with new path and auto folder creation
html 984697e jens-daniel-mueller 2020-12-12 Build site.
html 3ebff89 jens-daniel-mueller 2020-12-12 Build site.
html ba112d3 jens-daniel-mueller 2020-12-11 Build site.
Rmd 91b2b29 jens-daniel-mueller 2020-12-11 selectable basinmask, try 5
html b01a367 jens-daniel-mueller 2020-12-09 Build site.
Rmd 71c63b0 jens-daniel-mueller 2020-12-09 rerun with variable predictor assignment
html 24a632f jens-daniel-mueller 2020-12-07 Build site.
html 92dca91 jens-daniel-mueller 2020-12-07 Build site.
html 6a8004b jens-daniel-mueller 2020-12-07 Build site.
html 70bf1a5 jens-daniel-mueller 2020-12-07 Build site.
html 7555355 jens-daniel-mueller 2020-12-07 Build site.
html 143d6fa jens-daniel-mueller 2020-12-07 Build site.
html abc6818 jens-daniel-mueller 2020-12-03 Build site.
Rmd 992ba15 jens-daniel-mueller 2020-12-03 rebuild with variable inventory depth
html c8c2e7b jens-daniel-mueller 2020-12-03 Build site.
Rmd 83203db jens-daniel-mueller 2020-12-03 calculate cant with variable inventory depth
html 090e4d5 jens-daniel-mueller 2020-12-02 Build site.
html 7c25f7a jens-daniel-mueller 2020-12-02 Build site.
html ec8dc38 jens-daniel-mueller 2020-12-02 Build site.
html c987de1 jens-daniel-mueller 2020-12-02 Build site.
html f8358f8 jens-daniel-mueller 2020-12-02 Build site.
html b03ddb8 jens-daniel-mueller 2020-12-02 Build site.
Rmd 9183e8f jens-daniel-mueller 2020-12-02 revised assignment of era to eras
html 22d0127 jens-daniel-mueller 2020-12-01 Build site.
html 0ff728b jens-daniel-mueller 2020-12-01 Build site.
html 91435ae jens-daniel-mueller 2020-12-01 Build site.
Rmd 17d09be jens-daniel-mueller 2020-12-01 auto eras naming
html cf19652 jens-daniel-mueller 2020-11-30 Build site.
Rmd 0895ad6 jens-daniel-mueller 2020-11-30 rebuild with all plot output
Rmd 2842970 jens-daniel-mueller 2020-11-30 cleaned for eMLR part only
html 196be51 jens-daniel-mueller 2020-11-30 Build site.
Rmd 7a4b015 jens-daniel-mueller 2020-11-30 first rebuild on ETH server
Rmd bc61ce3 Jens Müller 2020-11-30 Initial commit

1 Version ID

The results displayed on this site correspond to the Version_ID: v_XXX

2 Required data

Required are:

  • cleaned and prepared GLODAPv2_2020 file
GLODAP <-
  read_csv(paste(path_version_data,
                 "GLODAPv2.2020_MLR_fitting_ready.csv",
                 sep = ""))

3 Predictor combinations

Find all possible combinations of following considered predictor variables:

  • sal, temp, aou, oxygen, nitrate, silicate
# the following code is a workaround to find all predictor combinations
# using the olsrr package and fit all models for one era, slab, and basin

i_basin <- unique(GLODAP$basin)[1]
i_era   <- unique(GLODAP$era)[1]

# subset one basin and era for fitting
GLODAP_basin_era <- GLODAP %>%
  filter(basin == i_basin, era == i_era,
         data_source == "obs")

i_gamma_slab <- unique(GLODAP_basin_era$gamma_slab)[3]
print(i_gamma_slab)

# subset one gamma slab
GLODAP_basin_era_slab <- GLODAP_basin_era %>%
  filter(gamma_slab == i_gamma_slab)

# fit the full linear model, i.e. all predictor combinations
lm_full <- lm(paste(
  params_local$MLR_target,
  paste(params_local$MLR_predictors, collapse = " + "),
  sep = " ~ "
),
data = GLODAP_basin_era_slab)

# fit linear models for all possible predictor combinations
# unfortunately, this functions does not provide model coefficients (yet)
lm_all <- ols_step_all_possible(lm_full)

# convert to tibble
lm_all <- as_tibble(lm_all)

# format model formula
lm_all <- lm_all %>% 
  select(n, predictors) %>% 
  mutate(model = str_replace_all(predictors, " ", " + "),
         model = paste(params_local$MLR_target, "~", model))

# remove helper objects
rm(i_gamma_slab,
   i_era,
   i_basin,
   GLODAP_basin_era,
   GLODAP_basin_era_slab,
   lm_full)

4 Apply predictor threshold

Select combinations with a total number of predictors in the range:

  • Minimum: 2
  • Maximum: 5
lm_all <- lm_all %>% 
  filter(n >= params_local$MLR_predictors_min,
         n <= params_local$MLR_predictors_max)

This results in a total number of MLR models of:

  • 56

5 Fit all models

Individual linear regression models were fitted for the chosen target variable:

  • cstar_tref

as a function of each predictor combination. Fitting was performed separately within each basin, era, and slab. Model diagnostics, such as the root mean squared error (RMSE), were calculated for each fitted model.

# prepare nested data frame
GLODAP %>%
  filter_all(any_vars(is.na(.)))

GLODAP_nested <- GLODAP %>% 
  # filter(gamma_slab %in% unique(GLODAP$gamma_slab)[c(5,6)]) %>%
  nest(data = !c(gamma_slab, era, basin, data_source))

# expand with model definitions
GLODAP_nested_lm <- expand_grid(
  GLODAP_nested,
  lm_all#[1:3,]
)

# fit models and extract tidy model output
GLODAP_nested_lm_fit <- GLODAP_nested_lm %>% 
  mutate(
    fit = map2(.x = data, .y = model, 
               ~ lm(as.formula(.y), data = .x)),
    tidied = map(fit, tidy),
    glanced = map(fit, glance),
    augmented = map(fit, augment),
    vif = map(fit, ols_vif_tol)
  )

6 Tidy models

# extract glanced model output (model diagnostics, such as AIC)
GLODAP_glanced <- GLODAP_nested_lm_fit %>%
  select(-c(data, fit, tidied, augmented, vif)) %>%
  unnest(glanced) %>% 
  rename(n_predictors = n)

# extract tidy model output (model coefficients)
GLODAP_tidy <- GLODAP_nested_lm_fit %>%
  select(-c(data, fit, glanced, augmented, vif)) %>%
  unnest(tidied)

# extract augmented model output (fitted values and residuals)
GLODAP_augmented <- GLODAP_nested_lm_fit %>% 
  select(-c(data, fit, tidied, glanced, vif)) %>% 
  unnest(augmented) 

# print(object.size(GLODAP_augmented), units = "MB")

# extract VIC from output
GLODAP_glanced_vif <- GLODAP_nested_lm_fit %>%
  select(-c(data, fit, tidied, augmented, glanced)) %>%
  unnest(vif)

# calculte max vif per model
GLODAP_glanced_vif_max <- GLODAP_glanced_vif %>% 
  group_by(data_source, era, basin, gamma_slab, model, n) %>% 
  summarise(vif_max = max(VIF)) %>% 
  ungroup()

# calculate RMSE from augmented output
GLODAP_glanced_rmse <- GLODAP_augmented %>%
  group_by(data_source, era, basin, gamma_slab, model) %>%
  summarise(rmse = sqrt(c(crossprod(.resid)) / length(.resid))) %>%
  ungroup()

# add RMSE and vif_max to glanced output
GLODAP_glanced <- full_join(GLODAP_glanced, GLODAP_glanced_rmse)
GLODAP_glanced <- full_join(GLODAP_glanced, GLODAP_glanced_vif_max)
rm(GLODAP_glanced_rmse)
rm(GLODAP_glanced_vif_max)


# extract input data
GLODAP_data <- GLODAP_nested_lm_fit %>%
  select(-c(fit, tidied, glanced, augmented, vif)) %>%
  unnest(data)

# append input data with augmented data
GLODAP_augmented <- bind_cols(
  GLODAP_data,
  GLODAP_augmented %>% select(.fitted, .resid)
)

rm(GLODAP, GLODAP_nested, GLODAP_nested_lm, GLODAP_nested_lm_fit, lm_all,
   GLODAP_data)

7 Prepare coeffcients

Coefficients are prepared for the mapping of Cant and the chosen target variable.

7.1 VIF threshold

To avoid multicollinearity among predictors, models were excluded with a VIF above:

  • 500

After removing models affected by multicollinearity, the targeted number of MLRs (10) was undercut in following fitting units:

GLODAP_glanced_clean <- GLODAP_glanced %>% 
  filter(vif_max <= params_local$vif_max)

GLODAP_glanced_clean %>% 
  group_by(data_source, era, basin, gamma_slab) %>% 
  count() %>% 
  filter(n < params_local$MLR_number)
# A tibble: 0 x 5
# Groups:   data_source, era, basin, gamma_slab [0]
# … with 5 variables: data_source <chr>, era <chr>, basin <chr>,
#   gamma_slab <chr>, n <int>

7.2 Predictor selection

Within each basin and slab, the following number of best linear regression models was selected:

  • 10

The criterion used to select the best models was:

  • rmse

The criterion was summed up for two adjacent eras, and the models with lowest summed values were selected.

# calculate RMSE sum for adjacent eras
lm_all <- GLODAP_glanced_clean  %>%
  select(data_source, basin, gamma_slab, model, era, AIC, BIC, rmse) %>% 
  arrange(era) %>% 
  group_by(data_source, basin, gamma_slab, model) %>% 
  mutate(rmse_sum = rmse + lag(rmse),
         aic_sum = AIC + lag(AIC),
         bic_sum = BIC + lag(BIC)
         ) %>% 
  ungroup() %>% 
  select(-c(era)) %>% 
  drop_na()

# subset models with lowest summed criterion
# chose which criterion is applied

if (params_local$MLR_criterion == "aic") {
  lm_best <- lm_all %>%
    group_by(data_source, basin, gamma_slab) %>%
    slice_min(order_by = aic_sum,
              with_ties = FALSE,
              n = params_local$MLR_number) %>%
    ungroup() %>%
    arrange(basin, gamma_slab, model)
} 

if (params_local$MLR_criterion == "bic") {
  lm_best <- lm_all %>%
    group_by(data_source, basin, gamma_slab) %>%
    slice_min(order_by = bic_sum,
              with_ties = FALSE,
              n = params_local$MLR_number) %>%
    ungroup() %>%
    arrange(basin, gamma_slab, model)
}

if (params_local$MLR_criterion == "rmse") {
  lm_best <- lm_all %>%
    group_by(data_source, basin, gamma_slab) %>%
    slice_min(order_by = rmse_sum,
              with_ties = FALSE,
              n = params_local$MLR_number) %>%
    ungroup() %>%
    arrange(basin, gamma_slab, model)
}


# print table
lm_best %>% 
  kable() %>%
  add_header_above() %>%
  kable_styling() %>%
  scroll_box(width = "100%", height = "400px")
data_source basin gamma_slab model AIC BIC rmse rmse_sum aic_sum bic_sum
obs Atlantic (-Inf,26] cstar_tref ~ sal + aou + oxygen + nitrate 2590.5074 2613.9559 8.0403150 15.942817 5733.041 5781.132
mod Atlantic (-Inf,26] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 2363.1882 2390.0727 7.3564632 11.148659 4204.439 4257.959
obs Atlantic (-Inf,26] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 2591.4166 2618.7732 8.0284077 15.896938 5732.082 5788.188
mod Atlantic (-Inf,26] cstar_tref ~ sal + aou + oxygen + silicate 2361.6788 2384.7226 7.3617099 11.158631 4201.756 4247.631
obs Atlantic (-Inf,26] cstar_tref ~ sal + aou + oxygen + silicate 2665.8737 2689.3222 8.9072731 17.646280 5898.762 5946.853
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate 2340.9786 2364.0225 7.1435141 10.906996 4175.183 4221.057
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate 2586.0867 2609.5352 7.9921665 15.868970 5725.696 5773.786
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 2326.5994 2353.4839 6.9754563 10.738841 4162.786 4216.307
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 2587.2239 2614.5805 7.9828030 15.833248 5725.823 5781.929
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + silicate 2324.7698 2347.8136 6.9771845 10.745933 4159.902 4205.777
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + silicate 2660.7146 2684.1631 8.8450542 17.550126 5890.109 5938.199
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + nitrate + silicate 2325.7240 2348.7679 6.9868682 10.767349 4162.920 4208.795
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + oxygen 2671.7509 2691.2913 9.0031172 17.967623 5925.517 5965.592
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + oxygen + nitrate 2347.8127 2370.8566 7.2148261 10.976829 4181.756 4227.630
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + oxygen + nitrate 2584.8491 2608.2976 7.9787385 15.850311 5723.861 5771.952
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 2327.6156 2354.5001 6.9857672 10.747533 4163.517 4217.037
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 2586.0519 2613.4085 7.9701016 15.818606 5724.429 5780.535
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + oxygen + silicate 2326.0362 2349.0800 6.9900392 10.758060 4161.040 4206.915
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + oxygen + silicate 2660.0722 2683.5207 8.8373376 17.537265 5888.936 5937.026
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + silicate 2337.7926 2356.9958 7.1312098 10.987632 4186.195 4224.424
mod Atlantic (26,26.5] cstar_tref ~ aou + oxygen + nitrate + silicate 4725.8177 4753.4738 5.7978577 11.462711 10464.064 10520.588
obs Atlantic (26,26.5] cstar_tref ~ sal + aou + nitrate 6365.8110 6389.8007 8.3966348 17.066480 13776.748 13825.443
mod Atlantic (26,26.5] cstar_tref ~ sal + aou + nitrate + silicate 4839.4335 4867.0896 6.2591790 12.449391 10738.738 10795.262
obs Atlantic (26,26.5] cstar_tref ~ sal + aou + nitrate + silicate 6367.2679 6396.0555 8.3940902 16.963085 13756.008 13814.443
mod Atlantic (26,26.5] cstar_tref ~ sal + aou + oxygen + nitrate 4744.5758 4772.2319 5.8716092 12.142941 10667.524 10724.047
obs Atlantic (26,26.5] cstar_tref ~ sal + aou + oxygen + nitrate 6148.0204 6176.8080 7.4274305 15.347642 13373.942 13432.376
mod Atlantic (26,26.5] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 4638.1170 4670.3825 5.4577865 10.986155 10334.075 10400.019
obs Atlantic (26,26.5] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 6147.6651 6181.2506 7.4176747 15.255594 13353.987 13422.161
obs Atlantic (26,26.5] cstar_tref ~ sal + aou + oxygen + silicate 6426.5113 6455.2989 8.6762366 17.085564 13776.354 13834.789
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate 4752.9219 4780.5780 5.9047242 11.945972 10607.991 10664.515
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate 6234.0452 6262.8329 7.7926806 15.762586 13472.901 13531.336
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + silicate 4637.6093 4669.8747 5.4559193 10.436647 10144.128 10210.072
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + silicate 6235.4451 6269.0307 7.7900713 15.713566 13464.223 13532.397
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + silicate 6458.5820 6487.3696 8.8329093 17.402059 13847.359 13905.794
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + oxygen + nitrate 4758.3178 4785.9739 5.9262335 12.010675 10626.325 10682.848
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + oxygen + nitrate 6293.9865 6322.7741 8.0577494 16.139921 13561.769 13620.204
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 4648.5864 4680.8518 5.4964263 10.488920 10159.390 10225.334
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 6294.2800 6327.8656 8.0500798 16.100134 13555.829 13624.003
mod Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate 4711.0534 4738.7095 5.7404608 11.220878 10389.191 10445.715
mod Atlantic (26,26.5] cstar_tref ~ temp + oxygen + nitrate + silicate 4707.1904 4734.8465 5.7255374 11.172046 10374.057 10430.580
mod Atlantic (26.5,26.75] cstar_tref ~ aou + oxygen + nitrate 4765.1225 4788.3674 5.2635411 11.646350 11739.102 11787.196
mod Atlantic (26.5,26.75] cstar_tref ~ aou + oxygen + nitrate + silicate 4722.0832 4749.9771 5.1122185 11.487005 11695.386 11753.099
obs Atlantic (26.5,26.75] cstar_tref ~ sal + aou + nitrate + silicate 7160.0518 7189.6230 8.0170983 16.480869 16279.278 16339.782
obs Atlantic (26.5,26.75] cstar_tref ~ sal + aou + oxygen + nitrate 6946.6770 6976.2482 7.2216503 14.799494 15782.634 15843.137
obs Atlantic (26.5,26.75] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 6944.4863 6978.9860 7.2068448 14.717161 15759.510 15830.097
obs Atlantic (26.5,26.75] cstar_tref ~ sal + aou + oxygen + silicate 7216.3053 7245.8765 8.2410253 16.213439 16182.305 16242.809
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate 4753.8862 4781.7801 5.2186112 11.554422 11714.139 11771.851
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate 6917.0813 6946.6525 7.1177383 14.609469 15723.757 15784.261
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + silicate 4694.4329 4726.9758 5.0149826 11.348350 11655.865 11723.196
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + silicate 6916.6964 6951.1961 7.1094301 14.531844 15701.557 15772.144
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate 7248.5485 7278.1197 8.3721841 16.487629 16260.105 16320.609
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + oxygen + nitrate 4763.8831 4791.7771 5.2525097 11.663554 11749.256 11806.968
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + oxygen + nitrate 6919.0273 6948.5985 7.1245247 14.626204 15729.103 15789.607
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 4684.8043 4717.3472 4.9838058 11.393759 11671.814 11739.146
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 6920.0467 6954.5465 7.1211042 14.559852 15710.539 15781.127
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + oxygen + silicate 7248.4822 7278.0534 8.3719123 16.510544 16267.349 16327.852
mod Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate 4753.3825 4776.6275 5.2236710 11.716233 11763.642 11811.736
mod Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate 4692.5214 4720.4153 5.0152701 11.503015 11703.202 11760.915
mod Atlantic (26.5,26.75] cstar_tref ~ temp + oxygen + nitrate 4762.1210 4785.3659 5.2533188 11.818717 11796.120 11844.214
mod Atlantic (26.5,26.75] cstar_tref ~ temp + oxygen + nitrate + silicate 4683.5780 4711.4719 4.9863037 11.548885 11718.664 11776.377
obs Atlantic (26.75,27] cstar_tref ~ aou + oxygen + nitrate 7837.1212 7862.0494 9.0374453 18.406155 18260.303 18311.538
mod Atlantic (26.75,27] cstar_tref ~ aou + oxygen + nitrate + silicate 6105.9957 6136.0738 3.7570489 8.758703 15028.452 15090.293
obs Atlantic (26.75,27] cstar_tref ~ aou + oxygen + nitrate + silicate 7787.4800 7817.3938 8.8241359 18.136363 18195.440 18256.921
mod Atlantic (26.75,27] cstar_tref ~ aou + oxygen + silicate 6935.5608 6960.6259 5.4623304 11.022558 16167.488 16219.021
mod Atlantic (26.75,27] cstar_tref ~ sal + aou + nitrate + silicate 6328.6878 6358.7659 4.1531010 9.532688 15465.449 15527.289
obs Atlantic (26.75,27] cstar_tref ~ sal + aou + oxygen + nitrate 7622.6157 7652.5295 8.1762633 16.437025 17689.352 17750.834
obs Atlantic (26.75,27] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 7597.1736 7632.0731 8.0731384 16.271960 17644.475 17716.203
mod Atlantic (26.75,27] cstar_tref ~ sal + oxygen + nitrate + silicate 6784.4557 6814.5338 5.0986273 10.788885 16086.392 16148.232
mod Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + nitrate 6449.1607 6479.2387 4.3844907 11.055763 16219.037 16280.877
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + nitrate 7584.7512 7614.6651 8.0343139 16.134338 17595.525 17657.007
mod Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + silicate 6022.4224 6057.5135 3.6151087 8.498083 14876.229 14948.376
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + silicate 7550.0191 7584.9186 7.8989649 15.898609 17527.279 17599.007
mod Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate 6975.3877 7005.4658 5.5561160 11.091739 16196.267 16258.108
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + oxygen + nitrate 7577.1235 7607.0374 8.0060182 16.090436 17582.405 17643.886
mod Atlantic (26.75,27] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 6035.1316 6070.2227 3.6358454 8.531073 14896.312 14968.459
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 7537.4777 7572.3772 7.8532769 15.828271 17505.948 17577.676
obs Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate 7846.5466 7871.4749 9.0769310 18.462845 18274.954 18326.188
mod Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate 6060.0417 6090.1198 3.6801460 8.683943 14983.759 15045.599
obs Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate 7801.3963 7831.3101 8.8811181 18.219592 18217.372 18278.853
mod Atlantic (26.75,27] cstar_tref ~ temp + oxygen + nitrate + silicate 6131.7604 6161.8385 3.8008664 8.829289 15069.920 15131.761
mod Atlantic (27,27.25] cstar_tref ~ aou + oxygen + nitrate + silicate 6070.9059 6101.4763 2.9833466 6.149017 14449.928 14512.865
obs Atlantic (27,27.25] cstar_tref ~ aou + oxygen + nitrate + silicate 10567.2634 10598.9144 9.3551879 18.374151 24581.960 24647.018
obs Atlantic (27,27.25] cstar_tref ~ sal + aou + nitrate 10544.4593 10570.8351 9.2880388 18.454182 24619.801 24674.016
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate 6147.3795 6177.9499 3.0794504 7.005205 15226.641 15289.578
obs Atlantic (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate 10541.7993 10573.4503 9.2730639 18.362530 24586.631 24651.689
obs Atlantic (27,27.25] cstar_tref ~ sal + aou + oxygen + nitrate 10221.5927 10253.2437 8.2998628 16.880651 24043.550 24108.608
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 5849.8722 5885.5376 2.7198515 5.884190 14229.524 14302.951
obs Atlantic (27,27.25] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 10205.7529 10242.6791 8.2487509 16.551118 23902.057 23977.959
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + oxygen + silicate 7153.7669 7184.3373 4.6738961 8.642014 16267.955 16330.892
mod Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + nitrate 6622.7139 6653.2843 3.7502473 8.526377 16339.993 16402.931
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + nitrate 10207.5449 10239.1959 8.2595888 16.742341 23985.032 24050.091
mod Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + silicate 5700.7108 5736.3762 2.5567474 5.552390 13902.092 13975.519
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + silicate 10187.7909 10224.7171 8.1976068 16.316762 23797.738 23873.639
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + oxygen + nitrate 10214.6597 10246.3107 8.2799620 16.733091 23978.609 24043.667
mod Atlantic (27,27.25] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 5764.0525 5799.7179 2.6247797 5.696967 14047.535 14120.962
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 10194.9238 10231.8500 8.2178786 16.277083 23776.190 23852.091
mod Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate 6629.0487 6654.5240 3.7632288 8.644390 16415.111 16467.559
mod Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate 5707.1343 5737.7046 2.5656920 5.659429 14011.362 14074.300
obs Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate 10568.8727 10600.5237 9.3604026 18.372574 24580.654 24645.712
mod Atlantic (27,27.25] cstar_tref ~ temp + oxygen + nitrate + silicate 5768.4879 5799.0583 2.6317922 5.853371 14204.476 14267.414
mod Atlantic (27.25,27.5] cstar_tref ~ aou + oxygen + nitrate + silicate 8786.1542 8818.5249 3.5816596 6.170343 18823.626 18889.938
obs Atlantic (27.25,27.5] cstar_tref ~ aou + oxygen + nitrate + silicate 11796.8317 11829.3408 8.3138742 16.059708 27794.456 27861.425
obs Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate 11809.0032 11836.0941 8.3493098 16.055758 27781.116 27836.923
mod Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate 8723.0850 8755.4557 3.5129500 6.338903 19131.512 19197.823
obs Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate 11802.3313 11834.8403 8.3276079 16.029908 27773.961 27840.930
obs Atlantic (27.25,27.5] cstar_tref ~ sal + aou + oxygen + nitrate 11652.7531 11685.2622 7.9620369 15.544581 27552.113 27619.081
obs Atlantic (27.25,27.5] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 11618.2698 11656.1971 7.8753332 15.312483 27430.336 27508.466
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate 9260.4142 9292.7849 4.1432625 7.317711 20160.741 20227.052
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate 11648.9404 11681.4495 7.9529313 15.528379 27543.981 27610.950
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + silicate 8428.4364 8466.2022 3.2070382 5.383373 17733.972 17811.335
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + silicate 11617.0346 11654.9619 7.8724141 15.292563 27418.546 27496.676
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + oxygen + nitrate 9217.6411 9250.0117 4.0891896 7.090860 19881.236 19947.547
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + oxygen + nitrate 11653.6371 11686.1462 7.9641495 15.542929 27550.706 27617.675
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 8353.2657 8391.0315 3.1338460 5.128509 17290.085 17367.449
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 11624.4605 11662.3878 7.8899786 15.315923 27429.572 27507.702
mod Atlantic (27.25,27.5] cstar_tref ~ temp + aou + nitrate 9274.4649 9301.4404 4.1637374 7.340729 20176.179 20231.438
mod Atlantic (27.25,27.5] cstar_tref ~ temp + aou + nitrate + silicate 8488.1961 8520.5667 3.2684498 5.446619 17795.296 17861.608
obs Atlantic (27.25,27.5] cstar_tref ~ temp + aou + nitrate + silicate 11799.2387 11831.7477 8.3198822 16.066167 27797.131 27864.100
mod Atlantic (27.25,27.5] cstar_tref ~ temp + oxygen + nitrate 9228.8051 9255.7806 4.1057556 7.108155 19891.427 19946.686
mod Atlantic (27.25,27.5] cstar_tref ~ temp + oxygen + nitrate + silicate 8409.2187 8441.5894 3.1901242 5.189165 17353.311 17419.623
mod Atlantic (27.5,27.75] cstar_tref ~ aou + oxygen + nitrate 12914.0300 12942.6368 4.2249068 7.223679 27447.937 27506.380
mod Atlantic (27.5,27.75] cstar_tref ~ aou + oxygen + nitrate + silicate 12370.9195 12405.2476 3.7441091 6.171028 25686.006 25756.138
obs Atlantic (27.5,27.75] cstar_tref ~ aou + oxygen + nitrate + silicate 10724.9212 10756.9895 7.7006700 15.212360 23674.634 23739.946
obs Atlantic (27.5,27.75] cstar_tref ~ sal + aou + nitrate 10678.1072 10704.8308 7.5900071 14.984513 23566.606 23621.033
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate 12456.4625 12490.7906 3.8157708 7.214896 27715.437 27785.569
obs Atlantic (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate 10670.6690 10702.7373 7.5669041 14.911018 23535.416 23600.728
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + oxygen + nitrate 12788.4476 12822.7757 4.1071155 7.096136 27305.560 27375.692
obs Atlantic (27.5,27.75] cstar_tref ~ sal + aou + oxygen + nitrate 10605.0457 10637.1140 7.4082026 14.590801 23386.044 23451.357
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 12093.1721 12133.2215 3.5190214 5.945744 25409.792 25491.613
obs Atlantic (27.5,27.75] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 10600.4002 10637.8133 7.3923183 14.567911 23379.724 23455.921
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + nitrate 12866.0520 12900.3801 4.1783670 7.299828 27633.324 27703.456
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + nitrate 10598.3998 10630.4681 7.3923173 14.554244 23368.544 23433.856
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + nitrate + silicate 12167.1255 12207.1750 3.5771748 6.133272 25783.440 25865.260
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + nitrate + silicate 10595.7136 10633.1266 7.3811364 14.536343 23364.323 23440.520
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + oxygen + nitrate 10597.1664 10629.2347 7.3893728 14.546127 23364.590 23429.902
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 12300.5706 12340.6200 3.6845519 6.411677 26290.587 26372.408
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 10595.6954 10633.1085 7.3810931 14.533290 23362.720 23438.917
mod Atlantic (27.5,27.75] cstar_tref ~ temp + aou + nitrate + silicate 12517.8494 12552.1775 3.8680401 6.437819 26162.965 26233.097
obs Atlantic (27.5,27.75] cstar_tref ~ temp + aou + nitrate + silicate 10722.7728 10754.8411 7.6953282 15.198221 23668.072 23733.384
mod Atlantic (27.5,27.75] cstar_tref ~ temp + oxygen + nitrate + silicate 12669.9960 12704.3241 4.0006961 6.756213 26717.772 26787.904
mod Atlantic (27.75,27.85] cstar_tref ~ aou + oxygen + nitrate 7346.6814 7372.8306 3.4531242 5.871906 15082.874 15136.150
mod Atlantic (27.75,27.85] cstar_tref ~ aou + oxygen + nitrate + silicate 6642.8461 6674.2251 2.6739103 4.531041 13494.278 13558.209
obs Atlantic (27.75,27.85] cstar_tref ~ aou + oxygen + nitrate + silicate 8778.3248 8809.4053 6.8164722 14.153364 20049.442 20112.974
mod Atlantic (27.75,27.85] cstar_tref ~ sal + aou + nitrate + silicate 7300.2815 7331.6606 3.3930973 5.936251 15206.746 15270.678
obs Atlantic (27.75,27.85] cstar_tref ~ sal + aou + nitrate + silicate 8772.6255 8803.7059 6.8016940 14.176498 20060.751 20124.283
mod Atlantic (27.75,27.85] cstar_tref ~ sal + aou + oxygen + nitrate 7281.8705 7313.2496 3.3705385 5.752463 14968.529 15032.461
mod Atlantic (27.75,27.85] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 6642.9575 6679.5664 2.6720813 4.523417 13485.901 13560.487
obs Atlantic (27.75,27.85] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 8774.1125 8810.3730 6.8003655 14.126853 20042.547 20116.668
obs Atlantic (27.75,27.85] cstar_tref ~ sal + oxygen + nitrate + silicate 8772.1355 8803.2159 6.8004251 14.163379 20054.955 20118.486
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate 7585.7969 7617.1760 3.7629032 6.422394 15642.374 15706.305
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + silicate 6909.4463 6946.0552 2.9429474 5.017795 14134.907 14209.494
obs Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + silicate 8774.6213 8810.8818 6.8016831 14.147662 20051.823 20125.944
obs Atlantic (27.75,27.85] cstar_tref ~ sal + temp + nitrate + silicate 8787.3474 8818.4278 6.8399329 14.185912 20062.550 20126.081
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 7021.8999 7058.5088 3.0653312 5.275877 14459.970 14534.556
obs Atlantic (27.75,27.85] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 8773.3864 8809.6468 6.7984853 14.144203 20050.472 20124.592
mod Atlantic (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate 6925.8460 6957.2250 2.9626322 5.057378 14181.339 14245.270
obs Atlantic (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate 8778.4723 8809.5527 6.8168550 14.167813 20055.911 20119.442
obs Atlantic (27.75,27.85] cstar_tref ~ temp + nitrate + silicate 8788.0804 8813.9807 6.8470552 14.198490 20063.733 20116.676
mod Atlantic (27.75,27.85] cstar_tref ~ temp + oxygen + nitrate + silicate 7036.7701 7068.1492 3.0841252 5.310710 14497.102 14561.033
obs Atlantic (27.75,27.85] cstar_tref ~ temp + oxygen + nitrate + silicate 8778.2736 8809.3540 6.8163392 14.167646 20055.868 20119.400
mod Atlantic (27.85,27.95] cstar_tref ~ aou + oxygen + nitrate 6291.1302 6316.3851 3.6785612 6.892110 14334.222 14386.217
mod Atlantic (27.85,27.95] cstar_tref ~ aou + oxygen + nitrate + silicate 6098.0726 6128.3786 3.3804472 6.011588 13522.094 13584.488
obs Atlantic (27.85,27.95] cstar_tref ~ aou + oxygen + nitrate + silicate 12970.1635 13003.4895 7.2068201 14.302669 30063.599 30131.936
obs Atlantic (27.85,27.95] cstar_tref ~ oxygen + nitrate + silicate 13037.0884 13064.8601 7.3381036 14.628990 30265.618 30322.566
mod Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate 6089.0574 6119.3633 3.3672686 6.360859 13913.927 13976.320
mod Atlantic (27.85,27.95] cstar_tref ~ sal + aou + oxygen + nitrate 6288.8098 6319.1157 3.6716816 6.822539 14272.710 14335.104
mod Atlantic (27.85,27.95] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 6016.2729 6051.6298 3.2599102 5.879047 13428.092 13500.884
obs Atlantic (27.85,27.95] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 12856.5520 12895.4323 6.9918638 13.990593 29882.309 29962.036
mod Atlantic (27.85,27.95] cstar_tref ~ sal + oxygen + nitrate + silicate 6164.9276 6195.2335 3.4797993 6.695287 14211.894 14274.287
obs Atlantic (27.85,27.95] cstar_tref ~ sal + oxygen + nitrate + silicate 13034.1905 13067.5165 7.3286959 14.612359 30259.709 30328.046
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + silicate 6048.5271 6083.8840 3.3057871 6.022938 13574.458 13647.250
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + silicate 12870.3451 12909.2254 7.0171686 14.060961 29928.552 30008.278
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + nitrate + silicate 12899.7133 12933.0393 7.0750583 14.182605 30001.477 30069.814
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 6071.1550 6106.5119 3.3383568 6.105939 13654.206 13726.998
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 12872.8223 12911.7027 7.0217231 14.073588 29936.820 30016.547
mod Atlantic (27.85,27.95] cstar_tref ~ temp + aou + nitrate + silicate 6167.0827 6197.3886 3.4830501 6.230688 13725.669 13788.063
obs Atlantic (27.85,27.95] cstar_tref ~ temp + aou + nitrate + silicate 13000.1296 13033.4556 7.2636063 14.424083 30139.406 30207.743
obs Atlantic (27.85,27.95] cstar_tref ~ temp + nitrate + silicate 13006.5689 13034.3406 7.2796796 14.459860 30157.739 30214.686
mod Atlantic (27.85,27.95] cstar_tref ~ temp + oxygen + nitrate + silicate 6193.0269 6223.3328 3.5224239 6.322233 13810.036 13872.429
obs Atlantic (27.85,27.95] cstar_tref ~ temp + oxygen + nitrate + silicate 13004.1862 13037.5122 7.2713280 14.439837 30149.131 30217.468
mod Atlantic (27.95,28.05] cstar_tref ~ aou + oxygen + nitrate 7570.9017 7596.4634 5.2704382 9.423926 17906.267 17959.350
mod Atlantic (27.95,28.05] cstar_tref ~ aou + oxygen + nitrate + silicate 7571.1608 7601.8348 5.2667006 9.156298 17670.112 17733.812
obs Atlantic (27.95,28.05] cstar_tref ~ aou + oxygen + nitrate + silicate 15666.2283 15700.5457 7.8200380 14.039514 35560.554 35631.033
mod Atlantic (27.95,28.05] cstar_tref ~ sal + aou + nitrate + silicate 7739.5111 7770.1851 5.6406906 9.475782 17787.206 17850.906
mod Atlantic (27.95,28.05] cstar_tref ~ sal + aou + oxygen + nitrate 7557.9264 7588.6004 5.2383739 9.284809 17800.452 17864.152
obs Atlantic (27.95,28.05] cstar_tref ~ sal + aou + oxygen + nitrate 15603.9709 15638.2884 7.7126878 13.889949 35456.589 35527.067
mod Atlantic (27.95,28.05] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 7559.9261 7595.7124 5.2383732 8.948138 17488.948 17563.265
obs Atlantic (27.95,28.05] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 15505.7029 15545.7399 7.5428851 13.546344 35185.547 35267.771
obs Atlantic (27.95,28.05] cstar_tref ~ sal + oxygen + nitrate 15627.8896 15656.4875 7.7571991 14.060958 35602.648 35661.380
obs Atlantic (27.95,28.05] cstar_tref ~ sal + oxygen + nitrate + silicate 15589.3571 15623.6746 7.6877035 13.955640 35531.216 35601.694
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate 7667.1251 7697.7990 5.4767362 9.540863 17925.497 17989.197
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate 15624.5625 15658.8799 7.7480296 13.980279 35531.454 35601.932
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate + silicate 7661.9493 7697.7356 5.4607451 9.179873 17600.127 17674.444
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate + silicate 15529.0039 15569.0409 7.5820084 13.642158 35266.405 35348.630
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate 15527.1955 15561.5130 7.5823311 13.644507 35264.644 35335.123
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + oxygen + nitrate 15618.3952 15652.7126 7.7374275 13.958222 35514.020 35584.498
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 7695.3765 7731.1628 5.5356374 9.263047 17641.632 17715.950
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 15529.1952 15569.2322 7.5823306 13.643981 35268.113 35350.338
mod Atlantic (27.95,28.05] cstar_tref ~ temp + aou + nitrate + silicate 7662.0829 7692.7569 5.4654948 9.396001 17799.034 17862.734
mod Atlantic (27.95,28.05] cstar_tref ~ temp + oxygen + nitrate + silicate 7694.7789 7725.4529 5.5388020 9.488910 17849.798 17913.498
mod Atlantic (28.05,28.1] cstar_tref ~ aou + oxygen + nitrate 3416.1968 3440.1361 1.6505810 3.112891 8517.242 8567.456
mod Atlantic (28.05,28.1] cstar_tref ~ aou + oxygen + nitrate + silicate 3395.1899 3423.9170 1.6293129 3.077630 8471.023 8531.280
obs Atlantic (28.05,28.1] cstar_tref ~ sal + aou + nitrate 13268.3993 13296.2541 7.3626079 13.215347 30341.802 30399.121
obs Atlantic (28.05,28.1] cstar_tref ~ sal + aou + nitrate + silicate 13260.9950 13294.4208 7.3447932 13.196713 30335.649 30404.432
mod Atlantic (28.05,28.1] cstar_tref ~ sal + aou + oxygen + nitrate 3409.4278 3438.1549 1.6424422 3.093519 8490.650 8550.906
mod Atlantic (28.05,28.1] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 3383.8318 3417.3467 1.6170903 3.042036 8415.627 8485.926
obs Atlantic (28.05,28.1] cstar_tref ~ sal + oxygen + nitrate 13252.2071 13280.0619 7.3319616 13.142358 30286.721 30344.040
obs Atlantic (28.05,28.1] cstar_tref ~ sal + oxygen + nitrate + silicate 13245.3168 13278.7426 7.3151898 13.125583 30281.828 30350.610
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate 3441.0056 3469.7327 1.6719400 3.144107 8563.063 8623.320
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate 13268.7239 13302.1496 7.3594309 13.191258 30324.956 30393.739
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + silicate 3410.0633 3443.5782 1.6411793 3.084458 8478.035 8548.335
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + silicate 13260.4026 13299.3993 7.3398899 13.171473 30318.411 30398.657
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + oxygen + nitrate 3457.5336 3486.2607 1.6875898 3.173421 8605.737 8665.993
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + oxygen + nitrate 13251.9286 13285.3544 7.3276596 13.136946 30287.419 30356.202
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 3423.6914 3457.2063 1.6538356 3.108640 8514.175 8584.474
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 13246.5571 13285.5538 7.3137583 13.123037 30284.041 30364.287
mod Atlantic (28.05,28.1] cstar_tref ~ temp + aou + nitrate + silicate 3427.4295 3456.1566 1.6591937 3.137211 8560.710 8620.966
obs Atlantic (28.05,28.1] cstar_tref ~ temp + oxygen + nitrate 13303.8677 13331.7225 7.4301855 13.314809 30406.370 30463.689
mod Atlantic (28.05,28.1] cstar_tref ~ temp + oxygen + nitrate + silicate 3442.2020 3470.9290 1.6730678 3.164889 8601.791 8662.048
obs Atlantic (28.05,28.1] cstar_tref ~ temp + oxygen + nitrate + silicate 13290.5873 13324.0131 7.4009962 13.278757 30388.840 30457.623
mod Atlantic (28.1,28.15] cstar_tref ~ aou + oxygen + nitrate 4519.3947 4545.1872 1.3989049 2.494642 9422.979 9475.723
mod Atlantic (28.1,28.15] cstar_tref ~ aou + oxygen + nitrate + silicate 4510.8461 4541.7972 1.3931749 2.463829 9341.402 9404.695
obs Atlantic (28.1,28.15] cstar_tref ~ sal + aou + nitrate + silicate 12794.5180 12827.2672 9.6495246 15.534421 28756.774 28824.465
mod Atlantic (28.1,28.15] cstar_tref ~ sal + aou + oxygen + nitrate 4487.0136 4517.9647 1.3803152 2.431360 9257.676 9320.968
mod Atlantic (28.1,28.15] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 4485.4399 4521.5495 1.3783971 2.404634 9180.713 9254.554
obs Atlantic (28.1,28.15] cstar_tref ~ sal + oxygen + nitrate + silicate 12802.5573 12835.3064 9.6719194 15.554291 28762.668 28830.359
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + nitrate 4507.8459 4538.7970 1.3915495 2.457752 9324.901 9388.193
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + nitrate 12736.1417 12768.8908 9.4884550 15.465035 28775.664 28843.355
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + nitrate + silicate 4505.4928 4541.6024 1.3891944 2.429035 9243.432 9317.273
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + nitrate + silicate 12732.4480 12770.6553 9.4728899 15.350820 28690.784 28769.757
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + oxygen + nitrate 4517.2243 4548.1754 1.3966368 2.474741 9370.248 9433.540
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + oxygen + nitrate 12719.8762 12752.6253 9.4440567 15.400839 28742.814 28810.505
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 4514.1867 4550.2963 1.3939018 2.443558 9282.567 9356.408
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 12715.5685 12753.7758 9.4268953 15.282099 28654.543 28733.516
obs Atlantic (28.1,28.15] cstar_tref ~ temp + aou + nitrate 12738.8407 12766.1316 9.5013202 15.479826 28777.973 28834.382
mod Atlantic (28.1,28.15] cstar_tref ~ temp + aou + nitrate + silicate 4540.8863 4571.8374 1.4095550 2.506384 9449.699 9512.991
obs Atlantic (28.1,28.15] cstar_tref ~ temp + aou + nitrate + silicate 12733.2826 12766.0317 9.4806357 15.363259 28693.608 28761.299
obs Atlantic (28.1,28.15] cstar_tref ~ temp + oxygen + nitrate 12721.6309 12748.9218 9.4542871 15.411507 28742.936 28799.345
mod Atlantic (28.1,28.15] cstar_tref ~ temp + oxygen + nitrate + silicate 4551.0857 4582.0368 1.4151601 2.524854 9497.680 9560.972
obs Atlantic (28.1,28.15] cstar_tref ~ temp + oxygen + nitrate + silicate 12715.4908 12748.2399 9.4321220 15.296436 28660.235 28727.926
mod Atlantic (28.15,28.2] cstar_tref ~ aou + oxygen + nitrate 5665.4747 5693.6633 0.9453751 1.792059 12584.326 12642.126
mod Atlantic (28.15,28.2] cstar_tref ~ aou + oxygen + nitrate + silicate 5392.2630 5426.0893 0.8847151 1.686929 12015.512 12084.872
obs Atlantic (28.15,28.2] cstar_tref ~ aou + oxygen + nitrate + silicate 6286.1409 6314.8680 8.3126732 14.546784 13114.489 13172.950
mod Atlantic (28.15,28.2] cstar_tref ~ sal + aou + oxygen + nitrate 5461.4626 5495.2889 0.8995910 1.719108 12202.423 12271.783
mod Atlantic (28.15,28.2] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 5383.5635 5423.0275 0.8824371 1.684582 12008.341 12089.261
obs Atlantic (28.15,28.2] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 6282.9596 6316.4745 8.2884299 14.478304 13098.368 13166.572
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate 5533.6271 5567.4534 0.9153709 1.745607 12346.269 12415.629
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate 6249.6213 6278.3484 8.1432977 14.370038 13075.488 13133.948
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate + silicate 5459.6927 5499.1567 0.8987742 1.712180 12161.365 12242.285
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate + silicate 6247.3094 6280.8243 8.1235285 14.267313 13047.037 13115.241
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + oxygen + nitrate 5592.7289 5626.5552 0.9285003 1.767266 12461.748 12531.108
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + oxygen + nitrate 6239.9109 6268.6380 8.0988453 14.267148 13045.995 13104.455
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 5514.0426 5553.5066 0.9106223 1.731244 12264.430 12345.350
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 6237.5019 6271.0168 8.0787420 14.172499 13020.076 13088.281
obs Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate 6247.6950 6271.6342 8.1436360 14.371871 13072.065 13120.782
mod Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate 5479.5462 5513.3725 0.9035196 1.717598 12183.779 12253.139
obs Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate 6245.4045 6274.1315 8.1239639 14.271545 13044.429 13102.889
obs Atlantic (28.15,28.2] cstar_tref ~ temp + oxygen + nitrate 6238.0426 6261.9818 8.0994466 14.269694 13042.788 13091.505
mod Atlantic (28.15,28.2] cstar_tref ~ temp + oxygen + nitrate + silicate 5535.7820 5569.6083 0.9158463 1.737247 12289.405 12358.765
obs Atlantic (28.15,28.2] cstar_tref ~ temp + oxygen + nitrate + silicate 6235.6262 6264.3533 8.0793079 14.176688 13017.448 13075.908
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate 3105.5197 3132.8964 0.5818610 1.374117 10255.650 10313.075
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate 4635.6685 4658.7889 5.2202982 10.428215 10925.085 10972.858
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate 3105.8018 3138.6539 0.5815778 1.373827 10257.881 10326.792
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate 4632.8355 4660.5799 5.2035721 10.399716 10919.621 10976.949
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + oxygen + nitrate 3027.4388 3060.2908 0.5688023 1.355144 10134.456 10203.367
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + oxygen + nitrate 4637.6619 4665.4063 5.2202751 10.424649 10927.686 10985.013
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 3010.6226 3048.9500 0.5657768 1.350530 10107.469 10187.864
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 4634.7830 4667.1515 5.2033908 10.394718 10921.672 10988.553
mod Atlantic (28.2, Inf] cstar_tref ~ sal + oxygen + nitrate 4322.8905 4350.2672 0.8216310 1.640854 11674.521 11731.946
obs Atlantic (28.2, Inf] cstar_tref ~ sal + oxygen + nitrate 4640.4293 4663.5496 5.2368267 10.486123 10946.038 10993.811
mod Atlantic (28.2, Inf] cstar_tref ~ sal + oxygen + nitrate + silicate 4254.0925 4286.9446 0.8053073 1.611956 11514.605 11583.515
obs Atlantic (28.2, Inf] cstar_tref ~ sal + oxygen + nitrate + silicate 4636.9606 4664.7050 5.2178447 10.462703 10942.839 11000.166
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate 3033.1749 3066.0269 0.5697279 1.356575 10144.065 10212.976
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate 4635.4090 4663.1534 5.2124716 10.407099 10921.598 10978.925
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + silicate 3016.6692 3054.9966 0.5667473 1.352146 10118.467 10198.862
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + silicate 4633.3522 4665.7206 5.1984495 10.383125 10917.617 10984.499
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + oxygen + nitrate 3064.0270 3096.8791 0.5747320 1.359751 10160.914 10229.824
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + oxygen + nitrate 4625.6372 4653.3816 5.1787597 10.358813 10906.078 10963.405
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 3040.8732 3079.2006 0.5706488 1.354222 10128.660 10209.055
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 4624.5449 4656.9134 5.1681369 10.343485 10905.126 10972.008
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + oxygen + nitrate 20526.9477 20562.5181 9.7517664 25.314444 48819.573 48891.926
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 16473.2780 16513.3553 9.1564282 18.069649 38025.142 38107.231
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 20445.5641 20487.0630 9.6063503 24.873640 48610.035 48694.446
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate 16443.7652 16478.1172 9.1009856 17.135159 37373.547 37443.909
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate 19885.9730 19921.5435 8.6881281 23.746070 47954.666 48027.019
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 16001.5818 16041.6591 8.2509552 15.814785 36573.094 36655.183
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 19821.4513 19862.9502 8.5846147 23.327223 47748.400 47832.811
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate 16473.2682 16507.6202 9.1604520 17.508926 37632.224 37702.586
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + nitrate 16711.5017 16740.1284 9.6593580 18.363211 38117.412 38176.047
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + nitrate + silicate 16658.4835 16692.8355 9.5427518 18.212454 38042.915 38113.278
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + oxygen + nitrate 16434.3617 16468.7137 9.0821131 17.057779 37320.495 37390.857
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + oxygen + nitrate 19712.0358 19747.6063 8.4200644 23.358570 47726.641 47798.994
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 15988.9007 16028.9780 8.2278901 15.697977 36485.936 36568.025
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 19658.1688 19699.6676 8.3357320 22.954478 47527.813 47612.224
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + oxygen + silicate 16452.5117 16486.8637 9.1185747 17.411664 37571.716 37642.078
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + silicate 16674.4463 16703.0729 9.5806666 18.269439 38070.000 38128.635
obs Indo-Pacific (-Inf,26] cstar_tref ~ temp + aou + nitrate 20747.6439 20777.2860 10.1510169 26.030918 49175.323 49235.617
obs Indo-Pacific (-Inf,26] cstar_tref ~ temp + aou + nitrate + silicate 20512.9362 20548.5067 9.7271783 25.007932 48681.394 48753.746
obs Indo-Pacific (-Inf,26] cstar_tref ~ temp + oxygen + nitrate 20363.2271 20392.8692 9.4717113 25.017134 48646.318 48706.612
obs Indo-Pacific (-Inf,26] cstar_tref ~ temp + oxygen + nitrate + silicate 20182.6989 20218.2693 9.1652726 24.164376 48224.800 48297.153
mod Indo-Pacific (26,26.5] cstar_tref ~ aou + oxygen + nitrate + silicate 16404.5727 16439.4712 6.5842704 13.535750 35744.605 35815.293
obs Indo-Pacific (26,26.5] cstar_tref ~ aou + oxygen + nitrate + silicate 18892.2217 18927.2541 9.9945045 28.138069 43374.420 43445.149
mod Indo-Pacific (26,26.5] cstar_tref ~ aou + oxygen + silicate 16510.8366 16539.9187 6.7295087 13.697817 35862.787 35921.693
obs Indo-Pacific (26,26.5] cstar_tref ~ oxygen + nitrate + silicate 19005.7717 19034.9654 10.2247204 28.372925 43487.419 43546.360
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + aou + nitrate + silicate 18898.9624 18933.9949 10.0077909 28.140583 43377.794 43448.524
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 16175.1543 16215.8692 6.2842430 12.772692 35120.420 35202.889
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 18876.0647 18916.9359 9.9588037 28.086963 43355.448 43437.966
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + aou + oxygen + silicate 16266.8998 16301.7983 6.4040980 12.894196 35211.628 35282.316
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + oxygen + nitrate + silicate 18894.0250 18929.0575 9.9980573 28.146091 43377.619 43448.348
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + silicate 16266.8131 16307.5280 6.4014053 13.047370 35350.144 35432.613
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + silicate 18900.3075 18941.1786 10.0064991 28.089374 43365.515 43448.032
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate 16383.7105 16418.6090 6.5566457 13.208021 35469.725 35540.413
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + silicate 16578.0411 16612.9396 6.8185232 13.663764 35829.426 35900.114
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + silicate 18905.7656 18940.7981 10.0212183 28.391214 43458.263 43528.992
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 16288.9393 16329.6542 6.4300138 13.094122 35387.963 35470.432
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 18892.2969 18933.1681 9.9907138 28.047900 43349.446 43431.964
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + oxygen + silicate 16398.4108 16433.3093 6.5760991 13.244089 35498.786 35569.474
mod Indo-Pacific (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate 16454.8474 16489.7459 6.6513210 13.668347 35848.900 35919.587
obs Indo-Pacific (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate 18970.8565 19005.8889 10.1506017 28.301657 43455.394 43526.123
obs Indo-Pacific (26,26.5] cstar_tref ~ temp + oxygen + nitrate + silicate 18963.3769 18998.4093 10.1356497 28.260626 43439.765 43510.494
obs Indo-Pacific (26.5,26.75] cstar_tref ~ aou + oxygen + nitrate 18744.8503 18774.0734 9.5037134 25.485112 41835.375 41894.202
mod Indo-Pacific (26.5,26.75] cstar_tref ~ aou + oxygen + nitrate + silicate 14364.0992 14398.2327 6.4675572 12.442869 29895.922 29964.804
obs Indo-Pacific (26.5,26.75] cstar_tref ~ aou + oxygen + nitrate + silicate 18538.6438 18573.7116 9.1238312 24.935758 41572.448 41643.041
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + aou + nitrate + silicate 18803.5550 18838.6228 9.6098870 25.462155 41851.394 41921.987
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + aou + oxygen + nitrate 14468.4149 14502.5484 6.6238731 12.446297 29874.786 29943.669
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 14303.8919 14343.7143 6.3761016 11.886079 29445.307 29525.670
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + aou + oxygen + silicate 14304.5582 14338.6917 6.3799949 11.911693 29463.016 29531.898
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + oxygen + nitrate + silicate 18749.5060 18784.5738 9.5086596 25.328530 41786.077 41856.670
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + silicate 14326.6874 14366.5098 6.4094638 12.070075 29598.644 29679.007
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + silicate 18662.0409 18702.9534 9.3434401 25.169411 41702.735 41785.094
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate 14328.0044 14362.1378 6.4143328 12.089287 29610.210 29679.092
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + silicate 14344.9860 14379.1195 6.4393186 12.142150 29650.908 29719.791
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + silicate 18672.1843 18707.2521 9.3656964 25.359467 41768.971 41839.564
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 14324.3700 14364.1924 6.4060642 12.057931 29588.844 29669.207
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 18654.5607 18695.4732 9.3297568 25.143980 41691.165 41773.523
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + oxygen + silicate 14325.0779 14359.2114 6.4100368 12.079795 29602.849 29671.732
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + silicate 14378.1319 14406.5765 6.4913398 12.194492 29682.327 29739.729
obs Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate 18671.8939 18706.9617 9.3651634 25.197899 41712.942 41783.535
obs Indo-Pacific (26.5,26.75] cstar_tref ~ temp + nitrate + silicate 18688.4276 18717.6508 9.3992322 25.405462 41787.504 41846.331
obs Indo-Pacific (26.5,26.75] cstar_tref ~ temp + oxygen + nitrate + silicate 18662.2388 18697.3066 9.3474644 25.166142 41698.394 41768.987
obs Indo-Pacific (26.75,27] cstar_tref ~ aou + nitrate + silicate 22355.5355 22385.9475 7.6342330 21.723824 51417.686 51479.006
obs Indo-Pacific (26.75,27] cstar_tref ~ oxygen + nitrate + silicate 22216.1155 22246.5275 7.4715847 21.542424 51268.747 51330.066
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + aou + nitrate + silicate 22094.0355 22130.5299 7.3297486 21.400948 51148.850 51222.433
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + aou + oxygen 17281.1004 17310.5177 6.2718712 12.101316 35847.710 35907.017
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + aou + oxygen + nitrate 17274.2714 17309.5721 6.2614437 12.087814 35839.804 35910.972
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 17275.4375 17316.6217 6.2604597 12.079937 35836.067 35919.097
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + aou + oxygen + silicate 17283.0890 17318.3896 6.2718576 12.096516 35846.907 35918.076
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + oxygen + nitrate + silicate 22079.4639 22115.9583 7.3132694 21.275226 51078.568 51152.151
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + nitrate 17377.2015 17412.5021 6.3840939 12.554431 36277.253 36348.422
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + silicate 17379.0165 17420.2007 6.3838714 12.551525 36278.531 36361.561
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + silicate 22077.4821 22120.0589 7.3087728 21.178155 51031.034 51116.881
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate 22102.8917 22139.3861 7.3397823 21.415978 51160.244 51233.827
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + oxygen 17342.8016 17372.2188 6.3452297 12.453888 36182.264 36241.571
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + oxygen + nitrate 17338.2981 17373.5988 6.3374573 12.444680 36178.389 36249.558
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 17339.0774 17380.2615 6.3359994 12.442834 36180.798 36263.828
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 22073.9998 22116.5766 7.3048426 21.167717 51024.197 51110.044
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + oxygen + silicate 17344.7339 17380.0346 6.3451488 12.453653 36186.048 36257.217
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate 22101.0968 22137.5912 7.3377477 21.419642 51161.342 51234.925
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate 22119.1400 22149.5520 7.3605002 21.656837 51285.416 51346.736
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + oxygen + nitrate + silicate 22098.6331 22135.1275 7.3349558 21.397491 51149.045 51222.628
mod Indo-Pacific (27,27.25] cstar_tref ~ aou + oxygen + nitrate + silicate 21279.6571 21315.9172 7.3669503 14.061546 42952.897 43025.698
mod Indo-Pacific (27,27.25] cstar_tref ~ oxygen + nitrate + silicate 21278.3693 21308.5860 7.3677931 14.062695 42949.908 43010.575
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate 21291.7568 21328.0168 7.3812813 14.091822 42980.517 43053.317
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate 24957.1168 24994.4827 6.7753920 18.174947 57106.779 57182.161
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + oxygen + nitrate + silicate 21277.4147 21313.6748 7.3642975 14.027868 42920.350 42993.150
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + oxygen + nitrate + silicate 24878.2190 24915.5849 6.7043585 18.009422 56958.446 57033.827
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + silicate 21224.0536 21266.3570 7.2991050 13.998881 42904.340 42989.274
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + silicate 24429.7115 24473.3049 6.3127906 17.503293 56426.973 56514.918
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + silicate 24501.5527 24538.9186 6.3753675 17.769498 56647.245 56722.626
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate 24492.6456 24530.0114 6.3677864 17.636929 56546.324 56621.706
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 21144.1819 21186.4853 7.2060648 13.827929 42748.156 42833.091
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 24422.4555 24466.0490 6.3066747 17.500288 56422.035 56509.980
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + oxygen + silicate 21159.0829 21195.3429 7.2256528 13.904905 42817.354 42890.154
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + oxygen + silicate 24490.3075 24527.6733 6.3657978 17.755865 56633.024 56708.405
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate 21306.2802 21342.5403 7.3985197 14.105570 42991.646 43064.447
obs Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate 24699.3102 24736.6761 6.5460297 18.134969 56986.422 57061.803
obs Indo-Pacific (27,27.25] cstar_tref ~ temp + nitrate + silicate 24720.8096 24751.9478 6.5666107 18.274040 57090.780 57153.598
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + oxygen + nitrate + silicate 21228.3904 21264.6504 7.3065377 13.936434 42838.273 42911.073
obs Indo-Pacific (27,27.25] cstar_tref ~ temp + oxygen + nitrate + silicate 24693.3499 24730.7158 6.5408199 18.122563 56975.280 57050.661
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + oxygen + silicate 21239.9825 21270.1992 7.3225062 14.010225 42904.518 42965.185
obs Indo-Pacific (27.25,27.5] cstar_tref ~ aou + nitrate + silicate 25879.1438 25910.6725 5.9129828 17.837357 59937.905 60001.344
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate 25386.9966 25424.8310 5.5627876 17.178047 59218.309 59294.435
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + oxygen + nitrate + silicate 25335.0951 25372.9295 5.5272312 17.141074 59165.341 59241.468
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate 23058.4643 23095.7267 5.5468250 10.444741 48690.465 48765.868
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + silicate 21182.9342 21226.4070 4.2976018 7.896199 44191.128 44279.098
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + silicate 25273.3657 25317.5059 5.4838826 17.088901 59098.972 59187.787
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate 21181.6424 21218.9047 4.2980154 7.922051 44247.837 44323.240
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate 25679.1419 25716.9763 5.7672387 17.483688 59586.231 59662.357
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate 25285.1811 25323.0155 5.4932508 17.165698 59159.399 59235.526
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + oxygen + nitrate 23041.7521 23079.0145 5.5342408 10.421252 48654.766 48730.170
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 21107.5282 21151.0010 4.2537842 7.815188 44027.225 44115.196
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 25276.2003 25320.3404 5.4858034 17.099152 59108.075 59196.890
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + oxygen + silicate 21105.7594 21143.0218 4.2539178 7.843039 44089.493 44164.896
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + oxygen + silicate 25668.2816 25706.1160 5.7595056 17.473459 59573.509 59649.636
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + silicate 22595.1556 22632.4179 5.2083287 9.232847 46554.181 46629.584
obs Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + silicate 25818.9192 25856.7535 5.8676995 17.671288 59790.740 59866.867
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + silicate 22593.8719 22624.9238 5.2088357 9.237990 46560.705 46623.541
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + oxygen + nitrate + silicate 22511.1658 22548.4282 5.1492150 9.123039 46362.214 46437.617
obs Indo-Pacific (27.25,27.5] cstar_tref ~ temp + oxygen + nitrate + silicate 25831.3753 25869.2097 5.8767365 17.678334 59801.722 59877.849
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + oxygen + silicate 22509.2734 22540.3254 5.1492902 9.129433 46371.857 46434.693
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate 25011.1111 25048.9484 5.3022630 16.639159 61392.659 61469.262
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + oxygen + nitrate 25387.5318 25425.3692 5.5545471 16.756772 61656.103 61732.707
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + oxygen + nitrate + silicate 25101.7307 25139.5680 5.3619305 16.775643 61547.121 61623.724
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou 21601.7629 21633.2481 3.5678356 6.807589 45765.961 45829.674
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + nitrate 21590.4965 21628.2788 3.5619416 6.801175 45755.199 45831.656
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + nitrate 25289.5388 25327.3762 5.4877372 16.711346 61576.140 61652.744
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + nitrate + silicate 20449.5857 20493.6650 3.0890675 5.979204 43554.643 43643.842
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + nitrate + silicate 24875.7702 24919.9138 5.2130957 16.421858 61151.858 61241.229
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate 20457.6868 20495.4691 3.0929586 6.012468 43654.882 43731.339
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate 21239.8295 21277.6118 3.4096291 6.391423 44633.554 44710.011
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + oxygen + nitrate 21589.6373 21627.4196 3.5615603 6.802277 45758.603 45835.059
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + oxygen + nitrate 25253.7773 25291.6146 5.4635562 16.710594 61560.092 61636.696
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 20425.5958 20469.6751 3.0798457 5.960864 43501.231 43590.430
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 24852.0080 24896.1516 5.1978212 16.428120 61146.243 61235.614
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + oxygen + silicate 20432.3631 20470.1453 3.0832126 5.990731 43591.242 43667.698
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate 21594.5345 21626.0197 3.5646229 6.713673 45494.360 45558.073
obs Indo-Pacific (27.5,27.75] cstar_tref ~ temp + aou + nitrate + silicate 25181.8030 25219.6404 5.4152118 16.675021 61498.847 61575.451
obs Indo-Pacific (27.5,27.75] cstar_tref ~ temp + oxygen + nitrate 25351.2956 25382.8267 5.5311135 16.806757 61679.626 61743.462
mod Indo-Pacific (27.5,27.75] cstar_tref ~ temp + oxygen + nitrate + silicate 21968.6673 22006.4496 3.7338346 6.806633 45642.285 45718.742
obs Indo-Pacific (27.5,27.75] cstar_tref ~ temp + oxygen + nitrate + silicate 25166.0210 25203.8583 5.4046685 16.677166 61493.713 61570.317
obs Indo-Pacific (27.75,27.85] cstar_tref ~ aou + oxygen + nitrate 19116.4209 19146.6951 5.0268311 11.837270 38843.751 38903.980
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + aou + nitrate 11006.2240 11035.2452 2.2802229 4.486519 21727.113 21785.098
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + aou + nitrate + silicate 10915.6709 10950.4964 2.2375747 4.370104 21473.764 21543.347
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + aou + oxygen + nitrate 19090.6118 19126.9408 5.0046838 11.720788 38737.536 38809.810
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + nitrate + silicate 11016.1709 11045.1922 2.2848545 4.571284 21909.948 21967.933
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + oxygen + nitrate + silicate 10922.5909 10957.4164 2.2407356 4.421157 21588.313 21657.895
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate 10996.4192 11031.2447 2.2747384 4.422709 21589.477 21659.059
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate 19029.7083 19066.0374 4.9565203 11.644027 38651.422 38723.697
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + silicate 10917.5223 10958.1521 2.2375068 4.340113 21409.136 21490.315
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + silicate 19027.0479 19069.4318 4.9528539 11.605898 38620.237 38704.557
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + nitrate + silicate 10992.4918 11027.3173 2.2729167 4.559176 21887.907 21957.489
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + oxygen + nitrate 10996.5903 11031.4158 2.2748178 4.429252 21604.207 21673.789
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + oxygen + nitrate 19002.4287 19038.7578 4.9350977 11.608787 38611.923 38684.197
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 10919.1993 10959.8290 2.2382724 4.349080 21429.680 21510.859
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 18999.4614 19041.8453 4.9312069 11.571511 38581.326 38665.646
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate 19040.6788 19070.9530 4.9667386 11.712892 38711.977 38772.205
mod Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate 11287.6024 11322.4279 2.4139537 4.564980 21887.549 21957.131
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate 19042.3731 19078.7022 4.9664975 11.703963 38708.058 38780.333
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + oxygen + nitrate 19011.8447 19042.1189 4.9440514 11.670787 38666.113 38726.342
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + oxygen + nitrate + silicate 19013.3721 19049.7011 4.9436804 11.662234 38662.450 38734.725
obs Indo-Pacific (27.85,27.95] cstar_tref ~ aou + oxygen + nitrate + silicate 23377.8721 23415.6214 4.5227292 11.257637 50292.522 50368.101
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate 11327.7923 11363.1065 2.0317348 3.971636 22101.537 22171.996
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + aou + oxygen + nitrate 23391.7124 23429.4617 4.5305801 11.196396 50222.961 50298.541
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + oxygen + nitrate + silicate 11413.6863 11449.0006 2.0648170 4.067805 22352.889 22423.348
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate 11186.8203 11222.1345 1.9785843 3.882403 21863.499 21933.958
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate 23312.5445 23350.2938 4.4858554 11.129747 50117.147 50192.727
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + silicate 11057.5000 11098.6999 1.9303243 3.758198 21525.714 21607.916
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + silicate 23283.0115 23327.0523 4.4681645 11.111995 50089.541 50177.716
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + oxygen + nitrate 11198.7441 11234.0583 1.9830255 3.893080 21892.329 21962.788
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + oxygen + nitrate 23281.2848 23319.0341 4.4683175 11.095121 50065.059 50140.638
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 11072.0082 11113.2082 1.9355977 3.773150 21567.527 21649.730
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 23257.5242 23301.5651 4.4539164 11.080713 50043.290 50131.466
mod Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate 11587.4342 11616.8627 2.1341946 4.101615 22432.006 22490.722
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate 23321.4952 23352.9530 4.4920154 11.192289 50192.446 50255.429
mod Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate + silicate 11521.9684 11557.2826 2.1072906 4.012707 22202.983 22273.442
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate + silicate 23311.6768 23349.4261 4.4853677 11.174066 50170.642 50246.222
mod Indo-Pacific (27.85,27.95] cstar_tref ~ temp + oxygen + nitrate 11606.4125 11635.8410 2.1418245 4.114155 22463.871 22522.587
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + oxygen + nitrate 23291.3335 23322.7913 4.4750692 11.156091 50139.011 50201.994
mod Indo-Pacific (27.85,27.95] cstar_tref ~ temp + oxygen + nitrate + silicate 11541.8812 11577.1954 2.1151959 4.027891 22242.606 22313.065
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + oxygen + nitrate + silicate 23285.0719 23322.8212 4.4704386 11.141038 50122.123 50197.702
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + oxygen + nitrate 30821.6235 30854.3796 4.7549248 11.371053 72564.397 72630.900
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + aou + oxygen + nitrate 30675.3269 30714.6341 4.6872552 11.285972 72386.862 72466.666
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou 11234.0221 11263.2609 2.1668040 4.242799 22971.034 23029.830
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate 11185.0717 11220.1583 2.1453486 4.166873 22779.017 22849.572
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate 30609.6693 30648.9766 4.6576033 11.262135 72332.314 72412.118
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate + silicate 11008.9821 11049.9165 2.0720095 3.985971 22306.609 22388.923
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate + silicate 30263.3119 30309.1703 4.5033891 11.107919 71987.953 72081.058
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate 11129.7498 11164.8364 2.1222927 4.150051 22740.494 22811.049
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate 11140.3367 11175.4232 2.1266856 4.142478 22718.790 22789.344
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + oxygen 11234.9048 11264.1436 2.1671776 4.243309 22972.275 23031.070
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + oxygen + nitrate 11187.9033 11222.9899 2.1465354 4.168872 22784.040 22854.594
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + oxygen + nitrate 30600.5821 30639.8893 4.6535141 11.268117 72342.448 72422.252
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 11013.2117 11054.1460 2.0737219 3.990400 22318.578 22400.891
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 30308.7375 30354.5960 4.5232054 11.136685 72050.461 72143.566
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + oxygen + silicate 11130.1720 11165.2586 2.1224677 4.150389 22741.354 22811.908
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + silicate 11143.4461 11172.6849 2.1288090 4.161039 22764.208 22823.003
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + nitrate 30811.0064 30843.7624 4.7500478 11.371129 72563.220 72629.723
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + nitrate + silicate 30727.9407 30767.2480 4.7111526 11.328905 72475.810 72555.614
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + oxygen + nitrate 30814.7428 30847.4988 4.7517636 11.381189 72582.843 72649.346
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + oxygen + nitrate + silicate 30758.3905 30797.6978 4.7250387 11.353796 72527.220 72607.023
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + aou + nitrate 21676.7700 21707.9785 4.1994617 10.275853 51083.718 51147.052
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + aou + nitrate + silicate 8305.9265 8339.3923 2.0205219 3.911490 17006.699 17074.099
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + aou + nitrate + silicate 21668.3030 21705.7532 4.1936760 10.259527 51061.418 51137.418
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + oxygen + nitrate 21711.6072 21742.8157 4.2187760 10.301915 51128.676 51192.010
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + oxygen + nitrate + silicate 8286.7084 8320.1742 2.0106101 3.914408 17016.056 17083.457
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + oxygen + nitrate + silicate 21712.9020 21750.3522 4.2183841 10.298729 51127.781 51203.782
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate 21552.7194 21590.1697 4.1303134 10.206595 50961.502 51037.503
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + silicate 7963.2444 8002.2878 1.8499466 3.655202 16469.985 16548.619
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + silicate 21510.3461 21554.0380 4.1062431 10.172023 50905.353 50994.021
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate 8150.6945 8184.1603 1.9418367 3.922809 17047.971 17115.372
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate 8068.5108 8101.9766 1.9014270 3.780990 16743.717 16811.118
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + oxygen + nitrate 21553.2763 21590.7266 4.1306164 10.210925 50968.102 51044.102
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 7964.1414 8003.1848 1.8503712 3.654487 16468.214 16546.848
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 21529.9176 21573.6095 4.1168423 10.192528 50939.807 51028.474
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + oxygen + silicate 8150.7023 8184.1681 1.9418406 3.923081 17048.550 17115.951
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + nitrate + silicate 8018.0417 8051.5075 1.8770293 3.777420 16739.819 16807.220
obs Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + nitrate + silicate 21697.0890 21734.5392 4.2096070 10.292922 51116.422 51192.422
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + nitrate + silicate 8118.5818 8146.4699 1.9269315 3.882691 16959.726 17015.894
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + oxygen + nitrate + silicate 8019.8446 8053.3104 1.8778954 3.777513 16739.904 16807.304
obs Indo-Pacific (28.05,28.1] cstar_tref ~ temp + oxygen + nitrate + silicate 21697.0611 21734.5113 4.2095915 10.308749 51140.115 51216.115
obs Indo-Pacific (28.1,28.2] cstar_tref ~ aou + nitrate + silicate 19525.2233 19556.0841 3.8063492 9.408290 48982.532 49045.654
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + nitrate + silicate 13918.4913 13955.3262 1.8415611 3.855428 30717.172 30791.713
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + nitrate + silicate 19527.2201 19564.2531 3.8063475 9.396658 48967.050 49042.797
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + nitrate + silicate 13921.9321 13952.6278 1.8430239 3.857533 30721.138 30783.255
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + oxygen + nitrate + silicate 13887.8993 13924.7342 1.8333574 3.843787 30673.048 30747.588
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + oxygen + nitrate + silicate 19541.3094 19578.3424 3.8139276 9.395435 48966.370 49042.117
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + oxygen + silicate 13936.3134 13967.0092 1.8468962 3.895873 30869.918 30932.035
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + nitrate 19474.7144 19511.7474 3.7782317 9.432237 49020.724 49096.471
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + nitrate + silicate 13536.1366 13579.1106 1.7411047 3.719575 30196.340 30283.304
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + nitrate + silicate 19475.6915 19518.8967 3.7776860 9.367648 48916.937 49005.308
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate 13805.5517 13842.3866 1.8114559 3.857299 30729.030 30803.571
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + nitrate + silicate 13670.8145 13707.6494 1.7761835 3.765724 30373.220 30447.761
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + oxygen + nitrate 19470.2653 19507.2983 3.7758589 9.401234 48968.697 49044.444
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 13536.5787 13579.5527 1.7412171 3.719292 30195.195 30282.159
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 19468.5687 19511.7738 3.7738885 9.351293 48888.737 48977.109
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + oxygen + silicate 13805.6470 13842.4819 1.8114811 3.857915 30731.411 30805.951
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate 13830.3502 13861.0459 1.8185545 3.908157 30919.489 30981.606
obs Indo-Pacific (28.1,28.2] cstar_tref ~ temp + aou + nitrate + silicate 19496.8343 19533.8672 3.7900511 9.383629 48942.140 49017.887
obs Indo-Pacific (28.1,28.2] cstar_tref ~ temp + oxygen + nitrate 19479.2530 19510.1139 3.7817217 9.417359 48992.766 49055.888
obs Indo-Pacific (28.1,28.2] cstar_tref ~ temp + oxygen + nitrate + silicate 19475.3653 19512.3983 3.7785790 9.367684 48913.175 48988.922
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + nitrate + silicate 19469.3787 19509.8818 1.1296116 2.682498 49772.808 49855.344
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + silicate 19513.6617 19547.4143 1.1337593 2.744006 50406.107 50474.886
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + oxygen + nitrate + silicate 19466.0708 19506.5739 1.1293157 2.690242 49853.636 49936.171
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + oxygen + silicate 19477.2173 19510.9700 1.1304920 2.739329 50355.385 50424.165
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + nitrate 6029.3811 6059.3005 3.9029656 9.465579 15233.221 15294.882
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + nitrate + silicate 19129.8757 19177.1294 1.0994725 2.651065 49421.720 49518.011
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + nitrate + silicate 6023.0495 6057.9554 3.8879675 9.449952 15228.557 15300.495
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + oxygen 6027.3138 6057.2332 3.8992387 9.445375 15222.455 15284.116
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + oxygen + nitrate 5994.0156 6028.9215 3.8361520 9.382250 15191.137 15263.075
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + oxygen + silicate 6026.9630 6061.8690 3.8950052 9.439017 15222.981 15294.919
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + silicate 19286.8253 19327.3285 1.1133991 2.721413 50158.655 50241.190
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + nitrate 6028.1603 6053.0931 3.9043710 9.468227 15230.654 15282.039
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + nitrate + silicate 19252.4074 19292.9106 1.1103687 2.745977 50401.486 50484.022
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + nitrate + silicate 6022.3018 6052.2212 3.8902182 9.453130 15226.298 15287.960
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + oxygen + nitrate 6025.7861 6055.7055 3.8964870 9.460200 15230.205 15291.866
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 19128.7504 19176.0041 1.0993745 2.651977 49431.191 49527.483
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 6020.6726 6055.5786 3.8836996 9.446521 15226.622 15298.560
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + oxygen + silicate 19286.8878 19327.3910 1.1134046 2.721664 50161.209 50243.745
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + silicate 19288.4906 19322.2432 1.1137223 2.750602 50448.221 50517.001
obs Indo-Pacific (28.2,28.3] cstar_tref ~ temp + aou + oxygen + nitrate + silicate 5987.7598 6022.6657 3.8250783 9.474188 15238.839 15310.777
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + nitrate + silicate 8611.9568 8648.3522 0.9338469 1.852705 19669.419 19743.784
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 8375.8178 8418.2791 0.8995695 1.814739 19401.979 19488.738
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + oxygen + nitrate + silicate 2077.5484 2105.1664 3.6039883 9.317084 5797.060 5855.292
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + oxygen + silicate 8496.6013 8532.9967 0.9170827 1.934470 20397.273 20471.637
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + nitrate + silicate 9095.3748 9125.7043 1.0078153 1.927358 20157.003 20218.973
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + oxygen + nitrate + silicate 8373.9247 8410.3201 0.8995846 1.817601 19423.796 19498.161
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + oxygen + nitrate + silicate 2077.1141 2100.7867 3.6113817 9.403065 5810.638 5860.550
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + nitrate 2126.2453 2149.9178 3.8512519 9.351248 5799.205 5849.118
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + nitrate + silicate 8362.3523 8404.8136 0.8976693 1.813010 19390.061 19476.819
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + nitrate + silicate 2079.0467 2106.6647 3.6110630 9.057992 5742.644 5800.875
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + silicate 8489.0363 8525.4316 0.9159938 1.933529 20390.904 20465.269
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + silicate 2128.6726 2152.3451 3.8635072 9.335626 5795.677 5845.589
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + nitrate + silicate 9094.9748 9131.3702 1.0074355 1.922818 20121.058 20195.423
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + nitrate + silicate 2077.0525 2100.7250 3.6110903 9.125931 5753.171 5803.084
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + oxygen + nitrate 2125.0397 2148.7122 3.8451793 9.345892 5798.152 5848.065
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 8375.4469 8417.9082 0.8995171 1.814853 19403.111 19489.869
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 2079.0423 2106.6603 3.6110424 9.051050 5741.149 5799.380
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + oxygen + silicate 8496.3378 8532.7332 0.9170447 1.934673 20398.965 20473.330
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + oxygen + silicate 2128.5982 2152.2708 3.8631311 9.329448 5794.359 5844.272
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + silicate 2126.6793 2146.4064 3.8635412 9.402873 5805.991 5847.585
obs Indo-Pacific (28.4, Inf] cstar_tref ~ sal + aou + nitrate + silicate 355.4346 368.8392 2.9145566 8.474782 1495.882 1528.445
obs Indo-Pacific (28.4, Inf] cstar_tref ~ sal + temp + aou 373.5685 384.7390 3.3723701 8.645043 1492.900 1520.035
obs Indo-Pacific (28.4, Inf] cstar_tref ~ sal + temp + aou + nitrate 353.4495 366.8541 2.8729324 8.125372 1473.397 1505.959
obs Indo-Pacific (28.4, Inf] cstar_tref ~ sal + temp + aou + nitrate + silicate 355.0500 370.6887 2.8646268 8.106124 1476.246 1514.236
obs Indo-Pacific (28.4, Inf] cstar_tref ~ sal + temp + aou + silicate 375.1966 388.6013 3.3632948 8.633319 1496.347 1528.909
obs Indo-Pacific (28.4, Inf] cstar_tref ~ sal + temp + oxygen + nitrate 354.0168 367.4214 2.8847665 8.223173 1479.808 1512.371
obs Indo-Pacific (28.4, Inf] cstar_tref ~ sal + temp + oxygen + nitrate + silicate 355.4579 371.0966 2.8731066 8.179582 1481.090 1519.079
obs Indo-Pacific (28.4, Inf] cstar_tref ~ temp + aou + nitrate 358.6077 369.7783 3.0258877 8.527043 1493.211 1520.346
obs Indo-Pacific (28.4, Inf] cstar_tref ~ temp + aou + nitrate + silicate 360.5409 373.9455 3.0244225 8.525124 1497.114 1529.676
obs Indo-Pacific (28.4, Inf] cstar_tref ~ temp + oxygen + nitrate + silicate 360.8928 374.2974 3.0321442 8.639047 1504.350 1536.913

7.3 Target variable coefficients

A data frame to map the target variable is prepared.

# extend best model selection with coefficients

lm_best <- left_join(
  lm_best,
  GLODAP_tidy %>% select(data_source, basin, gamma_slab, era, model, term, estimate)
)

7.4 dcant coeffcients

A data frame of coefficient offsets is prepared to facilitate the direct mapping of dcant.

# subtract coefficients of adjacent era  
lm_best_dcant <- lm_best %>%
  select(-c(AIC, BIC, rmse, aic_sum, bic_sum, rmse_sum)) %>% 
  arrange(era) %>%
  group_by(data_source, basin, gamma_slab, model, term) %>%
  mutate(delta_coeff = estimate - lag(estimate)) %>%
  ungroup() %>%
  arrange(basin, gamma_slab, model, term) %>%
  drop_na() %>%
  select(-c(era,estimate))

# pivot to wide format
lm_best_dcant <- lm_best_dcant %>%
  pivot_wider(values_from = delta_coeff,
              names_from = term,
              names_prefix = "delta_coeff_",
              values_fill = 0)
# create table of target variable coefficients in wide format
lm_best_target <- lm_best %>% 
  select(-c(AIC, BIC, rmse, aic_sum, bic_sum, rmse_sum)) %>% 
  pivot_wider(names_from = "term",
              names_prefix = "coeff_",
              values_from = "estimate",
              values_fill = 0
              )
GLODAP_augmented_best <- left_join(
  lm_best_target %>% select(data_source, basin, gamma_slab, era, model),
  GLODAP_augmented
)

7.5 Write files

lm_best_target %>%
  write_csv(paste(path_version_data,
                  "lm_best_target.csv",
                  sep = ""))

lm_best_dcant %>%
  write_csv(paste(path_version_data,
                  "lm_best_dcant.csv",
                  sep = ""))

GLODAP_augmented_best %>%
  write_csv(paste(path_version_data,
                  "lm_residuals_best_models.csv",
                  sep = ""))

8 Model diagnotics

8.1 Selection criterion vs predictors

The selection criterion (rmse) was plotted against the number of predictors (limited to 2 - 5).

8.1.1 All models

GLODAP_glanced %>%
  group_split(era) %>%
  # head(1) %>%
  map(
    ~ ggplot(data = .x,
             aes(
               as.factor(n_predictors),
               !!sym(params_local$MLR_criterion),
               col = data_source
             )) +
      geom_hline(yintercept = c(0, 10)) +
      geom_boxplot() +
      facet_grid(gamma_slab ~ basin) +
      scale_color_brewer(palette = "Set1") +
      ylim(c(0, NA)) +
      labs(x = "Number of predictors",
           title = paste("Era:", unique(.x$era)))
  )
[[1]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
9480f94 jens-daniel-mueller 2021-06-29
d054c39 jens-daniel-mueller 2021-06-24
4fcd066 jens-daniel-mueller 2021-06-24
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
d3c3fde jens-daniel-mueller 2021-03-26
594d503 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
ac1a836 jens-daniel-mueller 2021-02-24
b03fbd3 jens-daniel-mueller 2021-02-24
3d3b4cc jens-daniel-mueller 2021-02-23
7b672f7 jens-daniel-mueller 2021-01-11
33ba23c jens-daniel-mueller 2021-01-07
318609d jens-daniel-mueller 2020-12-23
9d0b2d0 jens-daniel-mueller 2020-12-23
0aa2b50 jens-daniel-mueller 2020-12-23
2886da0 jens-daniel-mueller 2020-12-19
02f0ee9 jens-daniel-mueller 2020-12-18
7bcb4eb jens-daniel-mueller 2020-12-18
7131186 jens-daniel-mueller 2020-12-17
5b48ef5 jens-daniel-mueller 2020-12-17
f3a708f jens-daniel-mueller 2020-12-17 

[[2]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
9480f94 jens-daniel-mueller 2021-06-29
d054c39 jens-daniel-mueller 2021-06-24
4fcd066 jens-daniel-mueller 2021-06-24
3b6387f jens-daniel-mueller 2021-06-18
c6b3da6 jens-daniel-mueller 2021-06-14
439ee80 jens-daniel-mueller 2021-06-11
33ffcab jens-daniel-mueller 2021-06-10
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
d3c3fde jens-daniel-mueller 2021-03-26
594d503 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26

8.1.2 Best models

left_join(lm_best_target %>% select(basin, gamma_slab, era, model),
          GLODAP_glanced) %>%
  group_split(era) %>%
  # head(1) %>%
  imap(
  ~ ggplot(data = .x,
         aes("",
             !!sym(params_local$MLR_criterion),
             col = data_source)) +
  geom_hline(yintercept = c(0, 10)) +
  geom_boxplot() +
  facet_grid(gamma_slab ~ basin) +
  scale_color_brewer(palette = "Set1") +
  ylim(c(0, NA)) +
  labs(x = "Number of predictors pooled",
           title = paste("Era:", unique(.x$era)))
  )
[[1]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
9480f94 jens-daniel-mueller 2021-06-29
d054c39 jens-daniel-mueller 2021-06-24
4fcd066 jens-daniel-mueller 2021-06-24
c6b3da6 jens-daniel-mueller 2021-06-14
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
d3c3fde jens-daniel-mueller 2021-03-26
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
ac1a836 jens-daniel-mueller 2021-02-24
b03fbd3 jens-daniel-mueller 2021-02-24
3d3b4cc jens-daniel-mueller 2021-02-23
7b672f7 jens-daniel-mueller 2021-01-11
33ba23c jens-daniel-mueller 2021-01-07
318609d jens-daniel-mueller 2020-12-23
9d0b2d0 jens-daniel-mueller 2020-12-23
0aa2b50 jens-daniel-mueller 2020-12-23
2886da0 jens-daniel-mueller 2020-12-19
02f0ee9 jens-daniel-mueller 2020-12-18
7bcb4eb jens-daniel-mueller 2020-12-18
7131186 jens-daniel-mueller 2020-12-17
5b48ef5 jens-daniel-mueller 2020-12-17
f3a708f jens-daniel-mueller 2020-12-17 

[[2]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
9480f94 jens-daniel-mueller 2021-06-29
d054c39 jens-daniel-mueller 2021-06-24
4fcd066 jens-daniel-mueller 2021-06-24
3b6387f jens-daniel-mueller 2021-06-18
c6b3da6 jens-daniel-mueller 2021-06-14
439ee80 jens-daniel-mueller 2021-06-11
33ffcab jens-daniel-mueller 2021-06-10
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
d3c3fde jens-daniel-mueller 2021-03-26
594d503 jens-daniel-mueller 2021-03-26

8.2 RMSE correlation between eras

RMSE was plotted to compare the agreement for one model applied to two adjecent eras (ie check whether the same predictor combination performs equal in both eras).

8.2.1 All models

# find max rmse to scale axis
max_rmse <-
  max(c(lm_all$rmse,
        lm_all$rmse_sum - lm_all$rmse))

lm_all %>%
  ggplot(aes(rmse, rmse_sum - rmse, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  coord_equal(xlim = c(0, max_rmse),
              ylim = c(0, max_rmse)) +
  geom_abline(slope = 1,
              col = 'red') +
  facet_grid(data_source ~ basin)

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
9480f94 jens-daniel-mueller 2021-06-29
d054c39 jens-daniel-mueller 2021-06-24
4fcd066 jens-daniel-mueller 2021-06-24
3b6387f jens-daniel-mueller 2021-06-18
c6b3da6 jens-daniel-mueller 2021-06-14
439ee80 jens-daniel-mueller 2021-06-11
33ffcab jens-daniel-mueller 2021-06-10
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
ac1a836 jens-daniel-mueller 2021-02-24
b03fbd3 jens-daniel-mueller 2021-02-24
3d3b4cc jens-daniel-mueller 2021-02-23
7b672f7 jens-daniel-mueller 2021-01-11
33ba23c jens-daniel-mueller 2021-01-07
318609d jens-daniel-mueller 2020-12-23
9d0b2d0 jens-daniel-mueller 2020-12-23
0aa2b50 jens-daniel-mueller 2020-12-23
2886da0 jens-daniel-mueller 2020-12-19
02f0ee9 jens-daniel-mueller 2020-12-18
7bcb4eb jens-daniel-mueller 2020-12-18
7131186 jens-daniel-mueller 2020-12-17
5b48ef5 jens-daniel-mueller 2020-12-17
e4ca289 jens-daniel-mueller 2020-12-16
158fe26 jens-daniel-mueller 2020-12-15
7a9a4cb jens-daniel-mueller 2020-12-15
61b263c jens-daniel-mueller 2020-12-15
984697e jens-daniel-mueller 2020-12-12
3ebff89 jens-daniel-mueller 2020-12-12
ba112d3 jens-daniel-mueller 2020-12-11
24a632f jens-daniel-mueller 2020-12-07
6a8004b jens-daniel-mueller 2020-12-07
70bf1a5 jens-daniel-mueller 2020-12-07
7555355 jens-daniel-mueller 2020-12-07
143d6fa jens-daniel-mueller 2020-12-07
090e4d5 jens-daniel-mueller 2020-12-02
7c25f7a jens-daniel-mueller 2020-12-02
b03ddb8 jens-daniel-mueller 2020-12-02
91435ae jens-daniel-mueller 2020-12-01
196be51 jens-daniel-mueller 2020-11-30 
rm(max_rmse)

8.2.2 Best models

# find max rmse to scale axis
max_rmse <-
  max(c(lm_best$rmse,
        lm_best$rmse_sum - lm_best$rmse))

lm_best %>%
  ggplot(aes(rmse, rmse_sum - rmse, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  coord_equal(xlim = c(0, max_rmse),
              ylim = c(0, max_rmse)) +
  geom_abline(slope = 1,
              col = 'red') +
  facet_grid(data_source ~ basin)

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
9480f94 jens-daniel-mueller 2021-06-29
d054c39 jens-daniel-mueller 2021-06-24
4fcd066 jens-daniel-mueller 2021-06-24
3b6387f jens-daniel-mueller 2021-06-18
c6b3da6 jens-daniel-mueller 2021-06-14
439ee80 jens-daniel-mueller 2021-06-11
33ffcab jens-daniel-mueller 2021-06-10
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
ac1a836 jens-daniel-mueller 2021-02-24
b03fbd3 jens-daniel-mueller 2021-02-24
3d3b4cc jens-daniel-mueller 2021-02-23
7b672f7 jens-daniel-mueller 2021-01-11
33ba23c jens-daniel-mueller 2021-01-07
318609d jens-daniel-mueller 2020-12-23
9d0b2d0 jens-daniel-mueller 2020-12-23
0aa2b50 jens-daniel-mueller 2020-12-23
2886da0 jens-daniel-mueller 2020-12-19
02f0ee9 jens-daniel-mueller 2020-12-18
7bcb4eb jens-daniel-mueller 2020-12-18
7131186 jens-daniel-mueller 2020-12-17
a84ff3c jens-daniel-mueller 2020-12-17 
rm(max_rmse)

8.3 Predictor counts

The number of models where a particular predictor was included were counted for each basin, density slab and compared eras

# calculate cases of predictor used
lm_all_stats <- lm_best_dcant %>% 
  pivot_longer(starts_with("delta_coeff_"),
               names_to = "term",
               names_prefix = "delta_coeff_",
               values_to = "delta_coeff") %>% 
  filter(term != "(Intercept)",
         delta_coeff != 0) %>% 
  group_by(data_source, basin, gamma_slab) %>% 
  count(term) %>% 
  ungroup() %>% 
  pivot_wider(values_from = n, names_from = term)

# print table
lm_all_stats %>%
  gt(rowname_col = "gamma_slab",
     groupname_col = c("basin", "data_source")) %>% 
  summary_rows(
    groups = TRUE,
    fns = list(total = "sum")
  )
aou nitrate oxygen sal silicate temp
Atlantic - mod
(-Inf,26] 5 6 5 10 8 8
(26,26.5] 7 10 6 7 7 6
(26.5,26.75] 6 10 6 4 5 8
(26.75,27] 7 8 5 6 9 6
(27,27.25] 8 9 5 6 8 6
(27.25,27.5] 6 10 5 5 6 8
(27.5,27.75] 8 10 6 6 7 5
(27.75,27.85] 8 10 6 6 7 5
(27.85,27.95] 7 10 7 6 8 4
(27.95,28.05] 8 10 6 6 7 5
(28.05,28.1] 7 10 7 6 6 6
(28.1,28.15] 7 10 7 6 6 6
(28.15,28.2] 7 10 7 6 6 6
(28.2, Inf] 6 10 6 10 5 4
total 97.00 133.00 84.00 90.00 95.00 83.00
Indo-Pacific - mod
(-Inf,26] 4 7 4 10 7 9
(26,26.5] 7 6 6 7 10 6
(26.5,26.75] 6 6 6 9 9 6
(26.75,27] 6 6 8 10 5 6
(27,27.25] 4 8 7 5 10 6
(27.25,27.5] 5 6 5 6 8 10
(27.5,27.75] 4 6 4 9 7 10
(27.75,27.85] 5 10 3 9 7 6
(27.85,27.95] 5 10 5 6 6 8
(27.95,28.05] 4 5 4 10 6 10
(28.05,28.1] 4 8 4 7 10 8
(28.1,28.2] 3 6 4 10 10 6
(28.2,28.3] 4 5 4 10 10 6
(28.3,28.4] 5 7 5 10 10 5
total 66.00 96.00 69.00 118.00 115.00 102.00
Atlantic - obs
(-Inf,26] 6 6 7 10 6 7
(26,26.5] 8 8 5 10 6 5
(26.5,26.75] 7 7 6 10 7 6
(26.75,27] 8 10 6 6 5 6
(27,27.25] 8 10 5 8 6 5
(27.25,27.5] 8 10 5 8 6 5
(27.5,27.75] 8 10 5 8 6 5
(27.75,27.85] 5 10 5 6 10 6
(27.85,27.95] 4 10 6 5 10 6
(27.95,28.05] 5 10 7 9 6 5
(28.05,28.1] 4 10 6 8 5 6
(28.1,28.15] 5 10 5 6 6 8
(28.15,28.2] 6 10 6 5 6 8
(28.2, Inf] 6 10 6 10 5 4
total 88.00 131.00 80.00 109.00 90.00 82.00
Indo-Pacific - obs
(-Inf,26] 6 10 6 6 5 8
(26,26.5] 5 10 6 6 10 5
(26.5,26.75] 5 10 5 5 9 6
(26.75,27] 4 10 4 5 10 6
(27,27.25] 4 8 4 7 10 8
(27.25,27.5] 5 8 4 7 10 7
(27.5,27.75] 5 10 6 7 6 7
(27.75,27.85] 6 10 6 5 4 8
(27.85,27.95] 6 10 6 5 5 8
(27.95,28.05] 6 10 6 5 4 8
(28.05,28.1] 5 10 5 8 6 6
(28.1,28.2] 5 10 5 6 7 7
(28.2,28.3] 6 8 6 9 5 10
(28.3,28.4] 4 7 5 10 8 8
(28.4, Inf] 7 8 3 7 6 9
total 79.00 139.00 77.00 98.00 105.00 111.00

8.4 RMSE alternatives

8.4.1 AIC

AIC is an alternative criterion to RMSE to judge model quality, but not (yet) taken into account.

lm_all %>% 
  ggplot(aes(rmse, AIC, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  facet_grid(data_source~basin)

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
9480f94 jens-daniel-mueller 2021-06-29
d054c39 jens-daniel-mueller 2021-06-24
4fcd066 jens-daniel-mueller 2021-06-24
3b6387f jens-daniel-mueller 2021-06-18
c6b3da6 jens-daniel-mueller 2021-06-14
439ee80 jens-daniel-mueller 2021-06-11
33ffcab jens-daniel-mueller 2021-06-10
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
ac1a836 jens-daniel-mueller 2021-02-24
b03fbd3 jens-daniel-mueller 2021-02-24
3d3b4cc jens-daniel-mueller 2021-02-23
7b672f7 jens-daniel-mueller 2021-01-11
33ba23c jens-daniel-mueller 2021-01-07
318609d jens-daniel-mueller 2020-12-23
9d0b2d0 jens-daniel-mueller 2020-12-23
0aa2b50 jens-daniel-mueller 2020-12-23
2886da0 jens-daniel-mueller 2020-12-19
02f0ee9 jens-daniel-mueller 2020-12-18
7bcb4eb jens-daniel-mueller 2020-12-18
7131186 jens-daniel-mueller 2020-12-17
5b48ef5 jens-daniel-mueller 2020-12-17 
lm_best %>% 
  ggplot(aes(rmse, AIC, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  facet_grid(data_source~basin)

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
9480f94 jens-daniel-mueller 2021-06-29
d054c39 jens-daniel-mueller 2021-06-24
4fcd066 jens-daniel-mueller 2021-06-24
3b6387f jens-daniel-mueller 2021-06-18
c6b3da6 jens-daniel-mueller 2021-06-14
439ee80 jens-daniel-mueller 2021-06-11
33ffcab jens-daniel-mueller 2021-06-10
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
ac1a836 jens-daniel-mueller 2021-02-24
b03fbd3 jens-daniel-mueller 2021-02-24
3d3b4cc jens-daniel-mueller 2021-02-23
7b672f7 jens-daniel-mueller 2021-01-11
33ba23c jens-daniel-mueller 2021-01-07
318609d jens-daniel-mueller 2020-12-23
9d0b2d0 jens-daniel-mueller 2020-12-23
0aa2b50 jens-daniel-mueller 2020-12-23
2886da0 jens-daniel-mueller 2020-12-19
02f0ee9 jens-daniel-mueller 2020-12-18
7bcb4eb jens-daniel-mueller 2020-12-18
7131186 jens-daniel-mueller 2020-12-17
5b48ef5 jens-daniel-mueller 2020-12-17

8.4.2 AIC vs BIC

BIC is an alternative criterion to RMSE to judge model quality, but not (yet) taken into account.

lm_all %>% 
  ggplot(aes(AIC, BIC, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  facet_grid(data_source~basin)

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
9480f94 jens-daniel-mueller 2021-06-29
d054c39 jens-daniel-mueller 2021-06-24
4fcd066 jens-daniel-mueller 2021-06-24
3b6387f jens-daniel-mueller 2021-06-18
c6b3da6 jens-daniel-mueller 2021-06-14
439ee80 jens-daniel-mueller 2021-06-11
33ffcab jens-daniel-mueller 2021-06-10
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20

8.5 RMSE vs VIF

GLODAP_glanced %>% 
  ggplot(aes(rmse, log10(vif_max),col=data_source)) +
  geom_hline(yintercept = 1) +
  geom_point() +
  scale_color_brewer(palette = "Set1") +
  facet_grid(gamma_slab~basin)

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
9480f94 jens-daniel-mueller 2021-06-29
d054c39 jens-daniel-mueller 2021-06-24
4fcd066 jens-daniel-mueller 2021-06-24
3b6387f jens-daniel-mueller 2021-06-18
c6b3da6 jens-daniel-mueller 2021-06-14
439ee80 jens-daniel-mueller 2021-06-11
33ffcab jens-daniel-mueller 2021-06-10
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20

8.6 Residual patterns

8.6.1 Fitted vs true

Plotted are fitted vs actual target variable values, here: rparams_local$MLR_target`

i_binwidth <- 2

p_prop_prop(df = GLODAP_augmented,
            var1 = params_local$MLR_target,
            var2 = ".fitted") +
  facet_grid(basin_AIP ~ data_source) +
  labs(title = "All models")

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13 
GLODAP_augmented %>%
  ggplot(aes(!!sym(params_local$MLR_target) - .fitted)) +
  geom_histogram(binwidth = i_binwidth) +
  labs(title = "All models") +
  scale_x_continuous(limits = c(-50,50)) +
  facet_wrap(~data_source)

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13 
p_prop_prop(df = GLODAP_augmented_best,
            var1 = params_local$MLR_target,
            var2 = ".fitted") +
  labs(title = "Selected models") +
  facet_grid(basin_AIP ~ data_source)

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13 
GLODAP_augmented_best %>%
  ggplot(aes(!!sym(params_local$MLR_target) - .fitted)) +
  geom_histogram(binwidth = i_binwidth) +
  labs(title = "Selected models") +
    scale_x_continuous(limits = c(-50,50)) +
  facet_wrap(~data_source)

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13

8.6.2 Pooled

In the following, we present residual patterns vs latitude across all domains.

i_ylim <- c(-30,30)

GLODAP_augmented_best %>%
  ggplot(aes(lat, .resid)) +
  geom_bin2d(binwidth = i_binwidth) +
  geom_hline(yintercept = 0, col = "black") +
  scale_fill_viridis_c(trans = "log10") +
  labs(
    title = paste(
      "Target variable:",
      params_local$MLR_target,
      "| Selected models",
      "| All domains"
    )
  ) +
  facet_wrap(~data_source)

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
9480f94 jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29
d054c39 jens-daniel-mueller 2021-06-24
4fcd066 jens-daniel-mueller 2021-06-24
3b6387f jens-daniel-mueller 2021-06-18
c6b3da6 jens-daniel-mueller 2021-06-14
439ee80 jens-daniel-mueller 2021-06-11
33ffcab jens-daniel-mueller 2021-06-10
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
27ae473 jens-daniel-mueller 2021-02-24
4bc00ea jens-daniel-mueller 2021-02-24

Due to the few large residuals, we limit the y axis range for the plots below.

GLODAP_augmented_best %>%
  ggplot(aes(lat, .resid)) +
  geom_bin2d(binwidth = i_binwidth) +
  geom_hline(yintercept = 0, col = "black") +
  scale_fill_viridis_c(trans = "log10") +
  coord_cartesian(ylim = i_ylim) +
  labs(
    title = paste(
      "Target variable:",
      params_local$MLR_target,
      "| Selected models",
      "| All domains"
    )
  ) +
  facet_wrap(~data_source)

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
9480f94 jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29
d054c39 jens-daniel-mueller 2021-06-24
4fcd066 jens-daniel-mueller 2021-06-24
3b6387f jens-daniel-mueller 2021-06-18
c6b3da6 jens-daniel-mueller 2021-06-14
439ee80 jens-daniel-mueller 2021-06-11
33ffcab jens-daniel-mueller 2021-06-10
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
d3c3fde jens-daniel-mueller 2021-03-26
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
27ae473 jens-daniel-mueller 2021-02-24
4bc00ea jens-daniel-mueller 2021-02-24

8.6.3 By model domain

In the following, we present residual patterns vs latitude for separate model domains, ie basins, density slabs and eras.

GLODAP_augmented_best %>% 
  group_split(basin, data_source) %>% 
  # head(1) %>% 
  map(~ ggplot(data = .x, aes(lat, .resid)) +
      geom_bin2d(binwidth = i_binwidth) +
      geom_hline(yintercept = 0, col = "black") +
      scale_fill_viridis_c(trans = "log10") +
      facet_grid(gamma_slab ~ era) +
      coord_cartesian(ylim = i_ylim) +
      labs(
        title = paste(
          "Target:",
          params_local$MLR_target,
          "| best models | basin:",
          unique(.x$basin),
          "| data_source:",
          unique(.x$data_source)
        )
      ))
[[1]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
9480f94 jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29
d054c39 jens-daniel-mueller 2021-06-24
4fcd066 jens-daniel-mueller 2021-06-24
c6b3da6 jens-daniel-mueller 2021-06-14
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
d3c3fde jens-daniel-mueller 2021-03-26
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
27ae473 jens-daniel-mueller 2021-02-24
a1ba577 jens-daniel-mueller 2021-02-24
071743d jens-daniel-mueller 2021-02-24 

[[2]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
9480f94 jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29
d054c39 jens-daniel-mueller 2021-06-24
4fcd066 jens-daniel-mueller 2021-06-24
3b6387f jens-daniel-mueller 2021-06-18
c6b3da6 jens-daniel-mueller 2021-06-14
439ee80 jens-daniel-mueller 2021-06-11
33ffcab jens-daniel-mueller 2021-06-10
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
d3c3fde jens-daniel-mueller 2021-03-26
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
27ae473 jens-daniel-mueller 2021-02-24
a1ba577 jens-daniel-mueller 2021-02-24 

[[3]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
9480f94 jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29
d054c39 jens-daniel-mueller 2021-06-24
4fcd066 jens-daniel-mueller 2021-06-24
c6b3da6 jens-daniel-mueller 2021-06-14
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
d3c3fde jens-daniel-mueller 2021-03-26
594d503 jens-daniel-mueller 2021-03-26 

[[4]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
9480f94 jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29
d054c39 jens-daniel-mueller 2021-06-24
4fcd066 jens-daniel-mueller 2021-06-24
3b6387f jens-daniel-mueller 2021-06-18
c6b3da6 jens-daniel-mueller 2021-06-14
439ee80 jens-daniel-mueller 2021-06-11
33ffcab jens-daniel-mueller 2021-06-10
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
d3c3fde jens-daniel-mueller 2021-03-26
594d503 jens-daniel-mueller 2021-03-26

8.6.4 Latitudinal mean

GLODAP_augmented_best <- GLODAP_augmented_best %>%
  mutate(lat_grid = as.numeric(as.character(cut(
    lat,
    seq(-90, 90, 10),
    seq(-85, 85, 10)
  ))))

lat_residual <- GLODAP_augmented_best %>%
  group_by(data_source, basin, gamma_slab, era, lat_grid) %>% 
  summarise(.resid_mean = mean(.resid)) %>% 
  ungroup()

lat_residual %>%
  group_split(basin) %>% 
  map(
  ~ggplot(data = .x,
          aes(lat_grid, .resid_mean, col=era)) +
  geom_line() +
  geom_point() +
  geom_hline(yintercept = 0, col = "black") +
  facet_grid(gamma_slab ~ data_source) +
    labs(title = paste("Basin:", unique(.x$basin)))
  )
[[1]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
d054c39 jens-daniel-mueller 2021-06-24
4fcd066 jens-daniel-mueller 2021-06-24
3b6387f jens-daniel-mueller 2021-06-18
c6b3da6 jens-daniel-mueller 2021-06-14
439ee80 jens-daniel-mueller 2021-06-11
33ffcab jens-daniel-mueller 2021-06-10
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
fec3558 jens-daniel-mueller 2021-02-24 

[[2]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
9480f94 jens-daniel-mueller 2021-06-29
d054c39 jens-daniel-mueller 2021-06-24
4fcd066 jens-daniel-mueller 2021-06-24
3b6387f jens-daniel-mueller 2021-06-18
c6b3da6 jens-daniel-mueller 2021-06-14
439ee80 jens-daniel-mueller 2021-06-11
33ffcab jens-daniel-mueller 2021-06-10
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26

8.6.5 Latitudinal offset

# calculate residual offset for adjacent eras
lat_residual_offset <- lat_residual  %>%
  select(data_source, basin, gamma_slab, era, lat_grid, .resid_mean) %>% 
  arrange(era) %>% 
  group_by(data_source, basin, gamma_slab, lat_grid) %>% 
  mutate(.resid_mean_offset = .resid_mean - lag(.resid_mean)) %>% 
  ungroup() %>% 
  drop_na()

lat_residual_offset %>%
  ggplot(aes(lat_grid, .resid_mean_offset, col=data_source)) +
  geom_line() +
  geom_point() +
  geom_hline(yintercept = 0, col = "black") +
  facet_grid(gamma_slab ~ basin)

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
9480f94 jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29
d054c39 jens-daniel-mueller 2021-06-24
4fcd066 jens-daniel-mueller 2021-06-24
3b6387f jens-daniel-mueller 2021-06-18
c6b3da6 jens-daniel-mueller 2021-06-14
439ee80 jens-daniel-mueller 2021-06-11
33ffcab jens-daniel-mueller 2021-06-10
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
fec3558 jens-daniel-mueller 2021-02-24

8.6.6 Spatial mean

GLODAP_augmented_best <- GLODAP_augmented_best %>%
  mutate(lon_grid = as.numeric(as.character(cut(
    lon,
    seq(-360, 500, 10),
    seq(-355, 500, 10)
  ))))

spatial_residual <- GLODAP_augmented_best %>%
  group_by(data_source, gamma_slab, era, lat_grid, lon_grid) %>% 
  summarise(.resid_mean = mean(.resid)) %>% 
  ungroup()

spatial_obs <- GLODAP_augmented_best %>%
  distinct(gamma_slab, era, lat, lon)

spatial_residual %>%
  rename(lat = lat_grid,
         lon = lon_grid) %>%
  group_split(data_source, gamma_slab) %>%
  # head(1) %>%
  map(~ p_map_dcant_slab(
    df = .x,
    var = ".resid_mean",
    col = "divergent",
    title_text = "Residual distribution",
    subtitle_text = paste(
      "Data source:",
      unique(.x$data_source),
      "| density slab:",
      unique(.x$gamma_slab)
      ),
    legend_title = paste(params_local$MLR_target, "\nbias")
  ) +
    geom_raster(data = spatial_obs,
                aes(lon, lat)) +
    facet_grid(era ~ .)
  )
[[1]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[2]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[3]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[4]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[5]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[6]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[7]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[8]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[9]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[10]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[11]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[12]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[13]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[14]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[15]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[16]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[17]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[18]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[19]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[20]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[21]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[22]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[23]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
9480f94 jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[24]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[25]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[26]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[27]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[28]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[29]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[30]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[31]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[32]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[33]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[34]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[35]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29

8.6.7 Spatial offset

# calculate residual offset for adjacent eras
spatial_residual_offset <- spatial_residual  %>%
  select(data_source, gamma_slab, era, lat_grid, lon_grid, .resid_mean) %>% 
  arrange(era) %>% 
  group_by(data_source, gamma_slab, lat_grid, lon_grid) %>% 
  mutate(.resid_mean_offset = .resid_mean - lag(.resid_mean)) %>% 
  ungroup() %>% 
  drop_na()

spatial_obs <- spatial_obs %>% 
  count(lon, lat, gamma_slab) %>% 
  filter(n == 2) %>% 
  select(-n)

spatial_residual_offset %>%
  rename(lat = lat_grid,
         lon = lon_grid) %>%
  group_split(data_source, gamma_slab) %>%
  # head(1) %>%
  map(~ p_map_dcant_slab(
    df = .x,
    var = ".resid_mean_offset",
    col = "divergent",
    title_text = "Residual offset distribution (era 2 - era 1)",
    subtitle_text = paste(
      "Data source:",
      unique(.x$data_source),
      "| density slab:",
      unique(.x$gamma_slab)
      ),
    legend_title = paste(params_local$MLR_target, "\nbias offset")
  ) +
    geom_raster(data = spatial_obs,
                aes(lon, lat))
  )
[[1]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[2]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[3]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
9480f94 jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[4]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[5]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[6]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[7]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[8]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[9]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[10]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[11]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[12]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[13]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[14]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[15]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[16]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[17]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[18]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[19]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[20]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[21]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[22]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[23]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[24]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[25]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[26]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[27]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[28]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[29]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[30]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[31]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[32]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
c18a8b1 jens-daniel-mueller 2021-07-09
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[33]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[34]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
971ce87 jens-daniel-mueller 2021-07-13
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29 

[[35]]

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
955f5dc jens-daniel-mueller 2021-06-29
147cb07 jens-daniel-mueller 2021-06-29

9 VIF

GLODAP_glanced %>% 
  ggplot(aes(as.factor(n),
             log10(vif_max),
             col=data_source)) + 
  geom_hline(yintercept = log10(params_local$vif_max)) +
  scale_color_brewer(palette = "Set1") +
  geom_violin() +
  facet_grid(gamma_slab~basin)

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
9480f94 jens-daniel-mueller 2021-06-29
d054c39 jens-daniel-mueller 2021-06-24
4fcd066 jens-daniel-mueller 2021-06-24
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
d37a85d jens-daniel-mueller 2021-05-31
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
330dcd0 jens-daniel-mueller 2021-03-20 
GLODAP_glanced_vif %>% 
  ggplot(aes(Variables,
             log10(VIF),
             col = data_source)) + 
  geom_hline(yintercept = log10(params_local$vif_max)) +
  geom_violin() +
  scale_color_brewer(palette = "Set1") +
  facet_grid(gamma_slab~basin) +
  coord_flip()

Version Author Date
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
9480f94 jens-daniel-mueller 2021-06-29
d054c39 jens-daniel-mueller 2021-06-24
4fcd066 jens-daniel-mueller 2021-06-24
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
d37a85d jens-daniel-mueller 2021-05-31
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
d3c3fde jens-daniel-mueller 2021-03-26
594d503 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
330dcd0 jens-daniel-mueller 2021-03-20 

sessionInfo()
R version 4.0.3 (2020-10-10)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: openSUSE Leap 15.2

Matrix products: default
BLAS:   /usr/local/R-4.0.3/lib64/R/lib/libRblas.so
LAPACK: /usr/local/R-4.0.3/lib64/R/lib/libRlapack.so

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
 [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8    
 [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8   
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C            
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
 [1] tictoc_1.0.1     gt_0.2.2         corrr_0.4.3      broom_0.7.5     
 [5] kableExtra_1.3.1 knitr_1.30       olsrr_0.5.3      GGally_2.0.0    
 [9] lubridate_1.7.9  ggforce_0.3.3    metR_0.9.0       scico_1.2.0     
[13] patchwork_1.1.1  collapse_1.5.0   forcats_0.5.0    stringr_1.4.0   
[17] dplyr_1.0.5      purrr_0.3.4      readr_1.4.0      tidyr_1.1.2     
[21] tibble_3.0.4     ggplot2_3.3.3    tidyverse_1.3.0  workflowr_1.6.2 

loaded via a namespace (and not attached):
 [1] fs_1.5.0                 webshot_0.5.2            RColorBrewer_1.1-2      
 [4] httr_1.4.2               rprojroot_2.0.2          tools_4.0.3             
 [7] backports_1.1.10         utf8_1.1.4               R6_2.5.0                
[10] nortest_1.0-4            DBI_1.1.0                colorspace_1.4-1        
[13] withr_2.3.0              gridExtra_2.3            tidyselect_1.1.0        
[16] curl_4.3                 compiler_4.0.3           git2r_0.27.1            
[19] cli_2.1.0                rvest_0.3.6              xml2_1.3.2              
[22] sass_0.2.0               labeling_0.4.2           scales_1.1.1            
[25] checkmate_2.0.0          goftest_1.2-2            digest_0.6.27           
[28] foreign_0.8-80           rmarkdown_2.5            rio_0.5.16              
[31] pkgconfig_2.0.3          htmltools_0.5.0          highr_0.8               
[34] dbplyr_1.4.4             rlang_0.4.10             readxl_1.3.1            
[37] rstudioapi_0.11          farver_2.0.3             generics_0.0.2          
[40] jsonlite_1.7.1           zip_2.1.1                car_3.0-10              
[43] magrittr_1.5             Matrix_1.2-18            Rcpp_1.0.5              
[46] munsell_0.5.0            fansi_0.4.1              abind_1.4-5             
[49] lifecycle_1.0.0          stringi_1.5.3            whisker_0.4             
[52] yaml_2.2.1               carData_3.0-4            MASS_7.3-53             
[55] plyr_1.8.6               grid_4.0.3               blob_1.2.1              
[58] parallel_4.0.3           promises_1.1.1           crayon_1.3.4            
[61] lattice_0.20-41          haven_2.3.1              hms_0.5.3               
[64] pillar_1.4.7             reprex_0.3.0             glue_1.4.2              
[67] evaluate_0.14            RcppArmadillo_0.10.1.2.0 data.table_1.13.2       
[70] modelr_0.1.8             vctrs_0.3.5              tweenr_1.0.2            
[73] httpuv_1.5.4             cellranger_1.1.0         gtable_0.3.0            
[76] polyclip_1.10-0          reshape_0.8.8            assertthat_0.2.1        
[79] openxlsx_4.2.3           xfun_0.18                RcppEigen_0.3.3.7.0     
[82] later_1.1.0.1            viridisLite_0.3.0        ellipsis_0.3.1          
[85] here_0.1