Last updated: 2021-06-08

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: 2a0c856

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 2a0c856. 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/

Untracked files:
    Untracked:  emlr_obs_v_XXX.Rproj

Unstaged changes:
    Deleted:    Cant_eMLR.Rproj
    Modified:   code/Workflowr_project_managment.R

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 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, nitrate, silicate, phosphate, phosphate_star
# 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:

  • 112

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.

library(furrr)
library(purrr)

map(c("hello", "world"), ~.x)
future_map(c("hello", "world"), ~.x)



library(tictoc)
plan(sequential)

tic()
nothingness <- future_map(c(2, 2, 2), ~Sys.sleep(.x))
toc()

plan(multisession, workers = 3)

tic()
nothingness <- future_map(c(2, 2, 2), ~Sys.sleep(.x))
toc()


future_map_dbl(1:4, function(x){
  Sys.sleep(1)
  x^2
}, .progress = TRUE)
# 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)
  )


# test parallel processing

# expand with model definitions


# test multidplyr
 
# library(parallel)
# library(multidplyr)
# 
# n_cores <- detectCores()
# system.time(
# cluster <- new_cluster(n_cores - 2)
# )
# 
# cluster
# 
# system.time(
#   GLODAP_nested_lm_part <-
#     GLODAP_nested_lm %>%
#     #group_by(gamma_slab, era, basin, data_source, model, predictors, n) %>% 
#     partition(cluster)
# )
# 
# 
# GLODAP_nested_lm
# GLODAP_nested_lm_part
# 
# system.time(cluster %>%
#               cluster_library(c("tidyverse",
#                                 "broom",
#                                 "olsrr",
#                                 "purrr")))
# 
# cluster_call(cluster, search())
# 
# system.time(
# GLODAP_nested_lm_fit <- GLODAP_nested_lm_part %>%
#   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)
#   )
# )
# 
# system.time(
#   # 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)
# )
# 
# 
# system.time(
# GLODAP_nested_lm_fit <- GLODAP_nested_lm_fit %>%
#   collect()
# )


### test furrr

# library(furrr)
# 
# plan(multisession, workers = 10)
# 
# tic()
# GLODAP_nested_lm_fit_furrr <- GLODAP_nested_lm %>%
#   mutate(
#     fit = future_map2(.x = data, .y = model,
#                ~ lm(as.formula(.y), data = .x)),
#     tidied = future_map(fit, tidy),
#     glanced = future_map(fit, glance),
#     augmented = future_map(fit, augment),
#     vif = future_map(fit, ols_vif_tol)
#   )
# toc()
# 
# plan(sequential)

# print(object.size(GLODAP_nested), units = "MB")
# print(object.size(GLODAP_nested_lm), units = "MB")
# print(object.size(GLODAP_nested_lm_fit), units = "MB")

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:

  • 1000

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_eras <- 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(eras = paste(lag(era), era, sep = " --> "),
         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_eras <- lm_all_eras %>%
    group_by(data_source, basin, gamma_slab, eras) %>%
    slice_min(order_by = aic_sum,
              with_ties = FALSE,
              n = params_local$MLR_number) %>%
    ungroup() %>%
    arrange(basin, gamma_slab, eras, model)
} 

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

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


# print table
lm_best_eras %>% 
  kable() %>%
  add_header_above() %>%
  kable_styling() %>%
  scroll_box(width = "100%", height = "400px")
data_source basin gamma_slab model AIC BIC rmse eras rmse_sum aic_sum bic_sum
mod Atlantic (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 114.56657 128.08935 0.6485658 2000-2009 –> 2010-2019 4.698958 771.0141 803.6903
obs Atlantic (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 849.96994 869.82392 6.6752223 2000-2009 –> 2010-2019 14.632954 1276.1537 1310.5504
mod Atlantic (-Inf,26] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 62.03165 75.55443 0.3874998 2000-2009 –> 2010-2019 4.771737 736.5372 769.2134
mod Atlantic (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 86.80798 100.33076 0.4940420 2000-2009 –> 2010-2019 4.479191 739.5530 772.2292
obs Atlantic (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 842.53866 862.39264 6.4812495 2000-2009 –> 2010-2019 15.393179 1282.0855 1316.4823
mod Atlantic (-Inf,26] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 59.79918 73.32196 0.3791107 2000-2009 –> 2010-2019 4.251825 706.0191 738.6952
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate 852.18294 869.20064 6.7877585 2000-2009 –> 2010-2019 14.841080 1277.7757 1307.2586
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 850.56756 870.42153 6.6910713 2000-2009 –> 2010-2019 14.680442 1277.2195 1311.6163
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 850.61431 870.46829 6.6923129 2000-2009 –> 2010-2019 14.678505 1277.2193 1311.6161
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 56.87756 70.40034 0.3684058 2000-2009 –> 2010-2019 4.660886 726.5607 759.2369
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 844.93848 864.79246 6.5432658 2000-2009 –> 2010-2019 14.529305 1271.5412 1305.9380
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 85.30201 98.82479 0.4868013 2000-2009 –> 2010-2019 4.364007 731.7861 764.4623
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 85.69833 99.22111 0.4886964 2000-2009 –> 2010-2019 4.775377 755.0732 787.7494
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 114.09905 127.62183 0.6455999 2000-2009 –> 2010-2019 4.681740 769.7429 802.4190
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 850.82425 870.67822 6.6978903 2000-2009 –> 2010-2019 14.707663 1277.7772 1312.1739
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 57.12564 70.64842 0.3693029 2000-2009 –> 2010-2019 4.677427 727.6382 760.3144
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 85.12985 98.65263 0.4859804 2000-2009 –> 2010-2019 4.673373 749.1617 781.8379
obs Atlantic (-Inf,26] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 857.40643 877.26040 6.8751426 2000-2009 –> 2010-2019 15.366498 1291.2489 1325.6456
obs Atlantic (-Inf,26] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 856.49076 876.34473 6.8502062 2000-2009 –> 2010-2019 15.333397 1290.2197 1324.6165
obs Atlantic (-Inf,26] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 853.69077 873.54475 6.7745145 2000-2009 –> 2010-2019 15.221595 1286.9164 1321.3131
obs Atlantic (26,26.5] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 3994.37020 4020.91942 6.0998602 2000-2009 –> 2010-2019 11.569223 7249.1948 7301.2670
mod Atlantic (26,26.5] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 2095.48713 2122.97561 3.8821054 2000-2009 –> 2010-2019 8.352481 3456.8690 3508.4544
obs Atlantic (26,26.5] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 3996.34274 4027.31682 6.0997244 2000-2009 –> 2010-2019 11.546356 7248.8361 7309.5870
mod Atlantic (26,26.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 2051.10054 2078.58903 3.6590206 2000-2009 –> 2010-2019 7.665397 3361.8540 3413.4394
obs Atlantic (26,26.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 3864.35746 3895.33155 5.4809931 2000-2009 –> 2010-2019 10.938754 7118.9737 7179.7246
mod Atlantic (26,26.5] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2021.73227 2049.22076 3.5185108 2000-2009 –> 2010-2019 7.909683 3374.8553 3426.4407
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 2054.83375 2082.32223 3.6772792 2000-2009 –> 2010-2019 7.964907 3396.9323 3448.5177
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 3881.55976 3912.53385 5.5579348 2000-2009 –> 2010-2019 11.084615 7149.2267 7209.9775
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 2055.10880 2082.59728 3.6786280 2000-2009 –> 2010-2019 8.036035 3404.6656 3456.2510
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 3882.53669 3913.51078 5.5623366 2000-2009 –> 2010-2019 11.162862 7164.0076 7224.7585
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 2025.89652 2053.38501 3.5381011 2000-2009 –> 2010-2019 8.276525 3414.1817 3465.7671
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2026.59715 2054.08563 3.5414078 2000-2009 –> 2010-2019 8.397113 3426.1781 3477.7635
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 2056.12326 2083.61174 3.6836072 2000-2009 –> 2010-2019 8.113584 3413.3110 3464.8964
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 3892.79056 3923.76464 5.6087491 2000-2009 –> 2010-2019 11.270261 7185.5253 7246.2761
obs Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + phosphate 3992.95222 4019.50144 6.0928549 2000-2009 –> 2010-2019 11.654020 7265.0882 7317.1604
mod Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 2092.99597 2120.48445 3.8692322 2000-2009 –> 2010-2019 8.379726 3458.5054 3510.0908
obs Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 3994.93888 4025.91297 6.0927890 2000-2009 –> 2010-2019 11.575362 7254.2725 7315.0234
mod Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2089.24838 2116.73686 3.8499467 2000-2009 –> 2010-2019 8.396856 3458.4728 3510.0582
obs Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 3994.90258 4025.87666 6.0926098 2000-2009 –> 2010-2019 11.618018 7262.3301 7323.0809
obs Atlantic (26,26.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 3995.48221 4026.45629 6.0954722 2000-2009 –> 2010-2019 11.675898 7273.2139 7333.9648
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 2747.93325 2777.19337 4.1009560 2000-2009 –> 2010-2019 7.269130 4649.8502 4706.4479
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 4793.86995 4826.46703 5.2224013 2000-2009 –> 2010-2019 10.893748 8063.1747 8125.4946
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 2751.00431 2780.26442 4.1140143 2000-2009 –> 2010-2019 7.172723 4627.1119 4683.7095
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 4792.95603 4825.55312 5.2193348 2000-2009 –> 2010-2019 10.919861 8067.5569 8129.8767
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 2800.46774 2829.72786 4.3301563 2000-2009 –> 2010-2019 7.182782 4625.3767 4681.9743
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 4886.56631 4919.16340 5.5429720 2000-2009 –> 2010-2019 11.240841 8160.6861 8223.0059
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2806.67172 2835.93184 4.3580555 2000-2009 –> 2010-2019 7.156689 4617.5548 4674.1525
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 2757.70756 2786.96768 4.1426615 2000-2009 –> 2010-2019 7.067962 4601.0819 4657.6795
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 4790.06873 4822.66582 5.2096589 2000-2009 –> 2010-2019 10.908830 8064.4243 8126.7441
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 2807.65852 2836.91864 4.3625097 2000-2009 –> 2010-2019 7.157205 4617.5081 4674.1057
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + phosphate 4854.38310 4882.32346 5.4364858 2000-2009 –> 2010-2019 11.111987 8122.4437 8175.8607
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + phosphate_star 4851.40998 4879.35034 5.4261080 2000-2009 –> 2010-2019 11.127015 8124.0798 8177.4968
mod Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 2770.62481 2799.88493 4.1984287 2000-2009 –> 2010-2019 7.266565 4648.9913 4705.5890
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 4854.51057 4887.10766 5.4299474 2000-2009 –> 2010-2019 11.104727 8124.4398 8186.7596
mod Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2774.47259 2803.73271 4.2151853 2000-2009 –> 2010-2019 7.183225 4628.4934 4685.0910
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 4851.69134 4884.28843 5.4201181 2000-2009 –> 2010-2019 11.113208 8124.9451 8187.2649
mod Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 2783.38534 2808.46544 4.2630735 2000-2009 –> 2010-2019 7.327987 4658.9804 4707.4927
obs Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 4848.98920 4876.92956 5.4176728 2000-2009 –> 2010-2019 11.117896 8121.5353 8174.9523
mod Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 2778.32764 2807.58776 4.2320406 2000-2009 –> 2010-2019 7.151978 4620.3550 4676.9527
obs Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 4849.24247 4881.83956 5.4115945 2000-2009 –> 2010-2019 11.104596 8122.4804 8184.8002
mod Atlantic (26.75,27] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 4271.57383 4305.29865 2.4846952 2000-2009 –> 2010-2019 4.812411 7119.8838 7184.6840
mod Atlantic (26.75,27] cstar_tref ~ aou + silicate + phosphate + phosphate_star 4647.05376 4675.96075 3.0546006 2000-2009 –> 2010-2019 5.417268 7512.0235 7567.5666
mod Atlantic (26.75,27] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 4421.58545 4455.31026 2.6971974 2000-2009 –> 2010-2019 5.226647 7373.9545 7438.7548
mod Atlantic (26.75,27] cstar_tref ~ sal + aou + silicate + phosphate 4573.23449 4602.14148 2.9337056 2000-2009 –> 2010-2019 5.469595 7526.7873 7582.3304
mod Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate 4554.27530 4588.00011 2.9002609 2000-2009 –> 2010-2019 5.175650 7374.1195 7438.9197
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate 5682.95401 5717.13108 4.4291929 2000-2009 –> 2010-2019 10.807642 9708.2587 9773.3643
mod Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4574.30323 4608.02805 2.9322114 2000-2009 –> 2010-2019 5.281668 7434.2529 7499.0531
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5687.58149 5721.75856 4.4397161 2000-2009 –> 2010-2019 10.762035 9702.0495 9767.1552
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 5697.28995 5731.46701 4.4618752 2000-2009 –> 2010-2019 10.770484 9709.0967 9774.2023
mod Atlantic (26.75,27] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4581.57191 4615.29672 2.9438940 2000-2009 –> 2010-2019 5.326736 7459.1870 7523.9873
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5688.59857 5722.77563 4.4420324 2000-2009 –> 2010-2019 10.761224 9702.4602 9767.5658
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + silicate + phosphate_star 5706.69133 5735.98596 4.4880396 2000-2009 –> 2010-2019 10.893812 9735.2365 9791.0413
mod Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 4290.63188 4324.35670 2.5107352 2000-2009 –> 2010-2019 4.786442 7110.6507 7175.4509
obs Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 5688.30887 5722.48593 4.4413725 2000-2009 –> 2010-2019 10.910349 9730.8914 9795.9971
mod Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 4300.53917 4334.26398 2.5243797 2000-2009 –> 2010-2019 4.817025 7129.8423 7194.6426
obs Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5693.91081 5728.08787 4.4541500 2000-2009 –> 2010-2019 10.835650 9719.8017 9784.9073
mod Atlantic (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 4323.03559 4356.76040 2.5556381 2000-2009 –> 2010-2019 4.881489 7170.3420 7235.1423
obs Atlantic (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 5695.97408 5730.15114 4.4588654 2000-2009 –> 2010-2019 10.851918 9724.0825 9789.1881
obs Atlantic (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate_star 5700.41712 5729.71174 4.4736224 2000-2009 –> 2010-2019 10.897536 9732.4296 9788.2344
obs Atlantic (26.75,27] cstar_tref ~ temp + silicate + phosphate + phosphate_star 5695.97725 5725.27188 4.4634482 2000-2009 –> 2010-2019 10.914867 9733.2278 9789.0326
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 5420.29491 5456.27808 2.0606401 2000-2009 –> 2010-2019 4.379076 9407.0787 9476.5133
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 5567.75184 5603.73501 2.1846127 2000-2009 –> 2010-2019 4.663529 9672.1956 9741.6303
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate 5714.41322 5745.25594 2.3171494 2000-2009 –> 2010-2019 4.780601 9805.8554 9865.3709
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 5715.71060 5751.69377 2.3165044 2000-2009 –> 2010-2019 4.661377 9722.4264 9791.8611
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate 11079.47655 11117.58629 6.1507924 2000-2009 –> 2010-2019 12.048190 19051.5532 19125.5568
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 11053.53847 11091.64821 6.1043198 2000-2009 –> 2010-2019 11.887691 18976.9602 19050.9638
mod Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + silicate 5523.15822 5559.14139 2.1463543 2000-2009 –> 2010-2019 4.411320 9468.9216 9538.3562
mod Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate 5714.27018 5750.25336 2.3151828 2000-2009 –> 2010-2019 4.694756 9746.8116 9816.2462
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 11100.28113 11138.39087 6.1883230 2000-2009 –> 2010-2019 11.927409 19004.5474 19078.5510
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 11060.47604 11098.58578 6.1167152 2000-2009 –> 2010-2019 11.933163 18998.1097 19072.1133
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 11119.11260 11157.22234 6.2224917 2000-2009 –> 2010-2019 12.057162 19064.5415 19138.5450
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 11101.95511 11140.06485 6.1913528 2000-2009 –> 2010-2019 11.943728 19011.9851 19085.9887
obs Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + phosphate_star 11074.90245 11107.56794 6.1461647 2000-2009 –> 2010-2019 12.057294 19050.7750 19114.2066
mod Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate 5688.89637 5719.73909 2.2938417 2000-2009 –> 2010-2019 4.596808 9661.9103 9721.4257
mod Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 5690.43627 5726.41944 2.2934236 2000-2009 –> 2010-2019 4.554410 9633.1085 9702.5432
obs Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 11040.99902 11079.10876 6.0819793 2000-2009 –> 2010-2019 12.056900 19045.6204 19119.6240
mod Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5690.59150 5726.57467 2.2935647 2000-2009 –> 2010-2019 4.553375 9632.3490 9701.7836
obs Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 10999.95793 11038.06767 6.0094298 2000-2009 –> 2010-2019 11.852105 18948.8034 19022.8070
mod Atlantic (27,27.25] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 5604.52389 5640.50706 2.2166733 2000-2009 –> 2010-2019 4.478612 9547.9369 9617.3716
obs Atlantic (27,27.25] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 11011.23805 11049.34779 6.0292833 2000-2009 –> 2010-2019 11.906167 18974.6312 19048.6348
obs Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 13186.43099 13225.87812 5.8287821 2000-2009 –> 2010-2019 10.509588 21880.2887 21956.7631
mod Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate 11037.87289 11071.62353 3.5665945 2000-2009 –> 2010-2019 6.894732 20706.9313 20773.7936
mod Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 10896.26184 10935.63759 3.4437708 2000-2009 –> 2010-2019 6.733961 20525.0737 20603.0797
mod Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 10908.64410 10948.01985 3.4541920 2000-2009 –> 2010-2019 6.766220 20561.8269 20639.8329
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 13171.25657 13210.70370 5.8074569 2000-2009 –> 2010-2019 10.409666 21815.4980 21891.9724
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 11157.61842 11196.99417 3.6705579 2000-2009 –> 2010-2019 6.935905 20758.5086 20836.5146
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 13169.72008 13209.16721 5.8053020 2000-2009 –> 2010-2019 10.382503 21797.9967 21874.4711
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + silicate 11009.95406 11049.32981 3.5406502 2000-2009 –> 2010-2019 6.788853 20591.4504 20669.4564
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + phosphate 13237.78979 13271.60162 5.9043932 2000-2009 –> 2010-2019 10.521253 21889.3440 21954.8935
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + phosphate_star 13241.14197 13274.95379 5.9091759 2000-2009 –> 2010-2019 10.503397 21878.2938 21943.8433
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 13234.84058 13274.28771 5.8973389 2000-2009 –> 2010-2019 10.512247 21887.1562 21963.6306
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 10887.24084 10926.61659 3.4361983 2000-2009 –> 2010-2019 7.100889 20913.1828 20991.1888
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 13237.98853 13277.43566 5.9018248 2000-2009 –> 2010-2019 10.495938 21877.0712 21953.5456
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 11153.29597 11192.67172 3.6666883 2000-2009 –> 2010-2019 6.799295 20601.2978 20679.3038
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 13166.71605 13206.16318 5.8010911 2000-2009 –> 2010-2019 10.383424 21798.2763 21874.7508
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 10865.57138 10904.94713 3.4180763 2000-2009 –> 2010-2019 6.955672 20761.4810 20839.4870
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + phosphate + phosphate_star 13237.91533 13271.72715 5.9045722 2000-2009 –> 2010-2019 10.502609 21877.4997 21943.0492
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 10888.53310 10927.90885 3.4372821 2000-2009 –> 2010-2019 7.104403 20916.9177 20994.9237
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 13234.95085 13274.39798 5.8974960 2000-2009 –> 2010-2019 10.495056 21876.2312 21952.7056
mod Atlantic (27.25,27.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 11166.27339 11205.64913 3.6783183 2000-2009 –> 2010-2019 6.804252 20606.4189 20684.4249
obs Atlantic (27.5,27.75] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 14145.05777 14185.30302 5.0862868 2000-2009 –> 2010-2019 10.357071 25269.5679 25348.2936
obs Atlantic (27.5,27.75] cstar_tref ~ aou + silicate + phosphate + phosphate_star 14155.19288 14189.68881 5.0996065 2000-2009 –> 2010-2019 10.371534 25278.4854 25345.9646
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + nitrate + phosphate 18252.01169 18289.99688 2.1785363 2000-2009 –> 2010-2019 4.613911 33877.7343 33952.4751
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 17810.75482 17855.07087 2.0652449 2000-2009 –> 2010-2019 4.493884 33419.7496 33506.9471
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate 18558.69946 18596.68465 2.2605397 2000-2009 –> 2010-2019 4.702488 34202.6509 34277.3917
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 18525.94842 18570.26447 2.2510949 2000-2009 –> 2010-2019 4.685852 34151.9564 34239.1539
obs Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 14156.90023 14197.14548 5.0992849 2000-2009 –> 2010-2019 10.295421 25229.9752 25308.7009
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 18140.64965 18184.96569 2.1489840 2000-2009 –> 2010-2019 4.568684 33724.7087 33811.9062
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 18557.17406 18601.49011 2.2595797 2000-2009 –> 2010-2019 4.681870 34148.4685 34235.6660
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 14162.48736 14202.73262 5.1054287 2000-2009 –> 2010-2019 10.308925 25240.6663 25319.3920
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 18366.65132 18410.96737 2.2083028 2000-2009 –> 2010-2019 4.588231 33838.6429 33925.8404
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 14191.22146 14231.46672 5.1371432 2000-2009 –> 2010-2019 10.368259 25288.4899 25367.2156
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 18455.81360 18500.12965 2.2321533 2000-2009 –> 2010-2019 4.621635 33954.8958 34042.0933
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 18476.74637 18521.06242 2.2377899 2000-2009 –> 2010-2019 4.625352 33970.3945 34057.5921
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 18371.97047 18416.28652 2.2097185 2000-2009 –> 2010-2019 4.591412 33848.9762 33936.1738
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 14171.40179 14211.64705 5.1152468 2000-2009 –> 2010-2019 10.331330 25258.2932 25337.0189
obs Atlantic (27.5,27.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 14151.67776 14191.92302 5.0935487 2000-2009 –> 2010-2019 10.381318 25287.7895 25366.5152
obs Atlantic (27.5,27.75] cstar_tref ~ temp + aou + silicate + phosphate 14160.71450 14195.21043 5.1056787 2000-2009 –> 2010-2019 10.394053 25295.2392 25362.7184
obs Atlantic (27.5,27.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 14161.45106 14201.69632 5.1042886 2000-2009 –> 2010-2019 10.416052 25313.8886 25392.6143
obs Atlantic (27.5,27.75] cstar_tref ~ temp + silicate + phosphate + phosphate_star 14170.01834 14204.51428 5.1159265 2000-2009 –> 2010-2019 10.427859 25320.5708 25388.0500
obs Atlantic (27.75,27.85] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 9243.32439 9280.46541 5.3667911 2000-2009 –> 2010-2019 10.543936 16970.4426 17043.5501
obs Atlantic (27.75,27.85] cstar_tref ~ aou + silicate + phosphate + phosphate_star 9241.32462 9273.15978 5.3667915 2000-2009 –> 2010-2019 10.568757 16978.4858 17041.1494
mod Atlantic (27.75,27.85] cstar_tref ~ sal + aou + nitrate + silicate 5209.53261 5241.50323 1.3329684 2000-2009 –> 2010-2019 2.570461 9481.1193 9544.1337
obs Atlantic (27.75,27.85] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 9229.90361 9267.04463 5.3426593 2000-2009 –> 2010-2019 10.535901 16964.8385 17037.9461
mod Atlantic (27.75,27.85] cstar_tref ~ sal + aou + silicate + phosphate 4994.59016 5026.56079 1.2421489 2000-2009 –> 2010-2019 2.428579 9156.1966 9219.2110
obs Atlantic (27.75,27.85] cstar_tref ~ sal + aou + silicate + phosphate 9227.95610 9259.79126 5.3427535 2000-2009 –> 2010-2019 10.563823 16974.3474 17037.0110
obs Atlantic (27.75,27.85] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 9231.30421 9268.44523 5.3451726 2000-2009 –> 2010-2019 10.569795 16981.4085 17054.5160
mod Atlantic (27.75,27.85] cstar_tref ~ sal + nitrate + silicate + phosphate_star 5606.50157 5638.47220 1.5185152 2000-2009 –> 2010-2019 2.884160 10135.2774 10198.2918
mod Atlantic (27.75,27.85] cstar_tref ~ sal + silicate + phosphate + phosphate_star 5598.11784 5630.08847 1.5143414 2000-2009 –> 2010-2019 2.885861 10138.0961 10201.1105
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + silicate 5166.84891 5204.14798 1.3135570 2000-2009 –> 2010-2019 2.547607 9433.1659 9506.6826
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 4877.67076 4914.96982 1.1945884 2000-2009 –> 2010-2019 2.376301 9030.8793 9104.3960
obs Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 9228.84118 9265.98220 5.3407536 2000-2009 –> 2010-2019 10.550973 16971.9947 17045.1022
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5006.55683 5043.85589 1.2462199 2000-2009 –> 2010-2019 2.467057 9244.7777 9318.2944
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 5011.55609 5048.85515 1.2482670 2000-2009 –> 2010-2019 2.467722 9246.8202 9320.3370
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5000.49129 5037.79035 1.2437408 2000-2009 –> 2010-2019 2.460603 9230.2008 9303.7176
obs Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 9230.51088 9267.65190 5.3437489 2000-2009 –> 2010-2019 10.563484 16978.2590 17051.3665
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate_star 5088.54517 5120.51579 1.2810605 2000-2009 –> 2010-2019 2.585336 9497.3158 9560.3302
obs Atlantic (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 9243.37370 9280.51472 5.3668800 2000-2009 –> 2010-2019 10.550962 16973.8635 17046.9710
obs Atlantic (27.75,27.85] cstar_tref ~ temp + aou + silicate + phosphate 9241.37422 9273.20938 5.3668809 2000-2009 –> 2010-2019 10.577116 16982.5351 17045.1987
obs Atlantic (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 9243.72548 9280.86650 5.3675140 2000-2009 –> 2010-2019 10.559439 16978.0218 17051.1293
mod Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate 6311.85701 6343.73633 1.9758799 2000-2009 –> 2010-2019 4.894872 12818.2384 12881.1568
obs Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 19487.67947 19530.07002 5.3130198 2000-2009 –> 2010-2019 10.583800 36158.1980 36241.9033
mod Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 6312.40507 6349.59761 1.9749238 2000-2009 –> 2010-2019 4.816283 12750.4861 12823.8910
obs Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 18884.70789 18927.09843 4.8283815 2000-2009 –> 2010-2019 9.865577 35310.1779 35393.8833
mod Atlantic (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate 6086.81798 6118.69731 1.8330858 2000-2009 –> 2010-2019 4.567691 12423.0234 12485.9419
obs Atlantic (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate_star 18888.98056 18925.31531 4.8331883 2000-2009 –> 2010-2019 9.870658 35312.7452 35384.4927
obs Atlantic (27.85,27.95] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 19285.90101 19328.29156 5.1456534 2000-2009 –> 2010-2019 10.339410 35876.8370 35960.5424
obs Atlantic (27.85,27.95] cstar_tref ~ sal + silicate + phosphate + phosphate_star 19327.26084 19363.59559 5.1811680 2000-2009 –> 2010-2019 10.378449 35919.8640 35991.6114
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + silicate 6153.07632 6190.26887 1.8727730 2000-2009 –> 2010-2019 4.779011 12650.0371 12723.4420
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + phosphate 6547.65369 6579.53301 2.1374482 2000-2009 –> 2010-2019 4.889278 12900.2352 12963.1536
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 5787.93869 5825.13123 1.6581586 2000-2009 –> 2010-2019 4.392525 12125.9170 12199.3219
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 18723.02847 18765.41901 4.7061220 2000-2009 –> 2010-2019 9.712551 35115.3781 35199.0835
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5892.08701 5929.27956 1.7167343 2000-2009 –> 2010-2019 4.475922 12253.6324 12327.0373
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 18703.96824 18746.35879 4.6919144 2000-2009 –> 2010-2019 9.693582 35091.1739 35174.8793
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 18704.06548 18746.45603 4.6919868 2000-2009 –> 2010-2019 9.737330 35138.2727 35221.9781
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 5900.41757 5937.61011 1.7215080 2000-2009 –> 2010-2019 4.524948 12303.4585 12376.8634
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate 18727.81158 18764.14633 4.7111885 2000-2009 –> 2010-2019 9.757906 35161.4912 35233.2387
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5881.76475 5918.95729 1.7108376 2000-2009 –> 2010-2019 4.461492 12235.2318 12308.6367
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 18721.74029 18764.13083 4.7051604 2000-2009 –> 2010-2019 9.710434 35112.8421 35196.5475
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate_star 5982.89396 6014.77328 1.7706726 2000-2009 –> 2010-2019 4.803223 12588.8117 12651.7302
obs Atlantic (27.95,28.05] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 19944.08159 19986.59977 5.3948001 2000-2009 –> 2010-2019 12.102064 37711.5427 37795.2950
obs Atlantic (27.95,28.05] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 19528.74031 19571.25849 5.0568342 2000-2009 –> 2010-2019 11.618953 37179.2871 37263.0394
obs Atlantic (27.95,28.05] cstar_tref ~ sal + aou + silicate + phosphate_star 19614.68916 19651.13332 5.1265855 2000-2009 –> 2010-2019 11.703029 37274.8895 37346.6772
obs Atlantic (27.95,28.05] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 19664.93788 19707.45607 5.1652589 2000-2009 –> 2010-2019 11.772262 37351.9133 37435.6656
obs Atlantic (27.95,28.05] cstar_tref ~ sal + silicate + phosphate + phosphate_star 19793.40021 19829.84437 5.2712971 2000-2009 –> 2010-2019 11.892183 37489.5927 37561.3804
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate 9620.02828 9652.94796 3.5747474 2000-2009 –> 2010-2019 7.375130 17464.5187 17528.9973
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate_star 9649.87330 9682.79298 3.6047742 2000-2009 –> 2010-2019 7.403977 17493.4811 17557.9597
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 9483.61843 9522.02472 3.4387310 2000-2009 –> 2010-2019 7.238253 17329.4646 17404.6896
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 19647.00826 19689.52644 5.1508536 2000-2009 –> 2010-2019 11.750015 37327.6373 37411.3895
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 9495.37859 9533.78489 3.4500838 2000-2009 –> 2010-2019 7.246876 17339.1809 17414.4060
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 19592.38492 19634.90311 5.1072144 2000-2009 –> 2010-2019 11.683117 37254.1450 37337.8973
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate + phosphate_star 9714.34956 9747.26924 3.6705071 2000-2009 –> 2010-2019 7.456616 17548.1382 17612.6168
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 19677.48104 19719.99923 5.1753604 2000-2009 –> 2010-2019 11.890679 37451.3556 37535.1078
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 9469.49310 9507.89939 3.4251443 2000-2009 –> 2010-2019 7.175248 17278.1063 17353.3313
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + phosphate + phosphate_star 9645.08061 9678.00029 3.5999354 2000-2009 –> 2010-2019 7.399137 17488.6873 17553.1659
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + phosphate_star 9727.40734 9754.84040 3.6860303 2000-2009 –> 2010-2019 7.488121 17571.1757 17624.9078
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate 19780.18771 19816.63187 5.2604598 2000-2009 –> 2010-2019 11.980270 37555.6351 37627.4227
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 9495.62897 9534.03527 3.4503259 2000-2009 –> 2010-2019 7.247123 17339.4352 17414.6602
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 19631.21038 19673.72856 5.1381943 2000-2009 –> 2010-2019 11.728678 37304.8072 37388.5595
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate_star 9494.25623 9527.17591 3.4509325 2000-2009 –> 2010-2019 7.247738 17336.0683 17400.5469
obs Atlantic (28.05,28.1] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 14407.13707 14447.43044 5.2687034 2000-2009 –> 2010-2019 11.606104 26000.2709 26078.9273
obs Atlantic (28.05,28.1] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 14351.32074 14391.61411 5.2061328 2000-2009 –> 2010-2019 11.543330 25944.3412 26022.9976
obs Atlantic (28.05,28.1] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 14354.50116 14394.79453 5.2096780 2000-2009 –> 2010-2019 11.534130 25940.3826 26019.0389
obs Atlantic (28.05,28.1] cstar_tref ~ sal + silicate + phosphate + phosphate_star 14440.32057 14474.85774 5.3085303 2000-2009 –> 2010-2019 11.656490 26037.3579 26104.7776
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 14473.97727 14514.27063 5.3446221 2000-2009 –> 2010-2019 11.649820 26049.0466 26127.7030
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate 3558.77561 3588.45170 1.3336016 2000-2009 –> 2010-2019 2.327457 5726.6697 5784.1693
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate_star 3598.21990 3627.89598 1.3591576 2000-2009 –> 2010-2019 2.371393 5794.0785 5851.5781
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 3559.47456 3594.09665 1.3327668 2000-2009 –> 2010-2019 2.308014 5700.5267 5767.6096
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 3596.57205 3631.19415 1.3567737 2000-2009 –> 2010-2019 2.355109 5773.3302 5840.4131
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 14387.87636 14428.16973 5.2470274 2000-2009 –> 2010-2019 11.518449 25943.8992 26022.5556
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 14466.93731 14507.23068 5.3365747 2000-2009 –> 2010-2019 11.645054 26043.8515 26122.5079
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + phosphate_star 3644.68557 3674.36166 1.3898917 2000-2009 –> 2010-2019 2.413129 5857.0403 5914.5399
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + phosphate_star 14475.48329 14510.02046 5.3486344 2000-2009 –> 2010-2019 11.657541 26050.6379 26118.0576
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 3618.07273 3652.69483 1.3708849 2000-2009 –> 2010-2019 2.393299 5831.1998 5898.2827
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 14355.31272 14395.60609 5.2105830 2000-2009 –> 2010-2019 11.501664 25922.4339 26001.0902
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + phosphate + phosphate_star 3586.96778 3616.64386 1.3518178 2000-2009 –> 2010-2019 2.359650 5776.1742 5833.6738
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 3586.94005 3621.56215 1.3504993 2000-2009 –> 2010-2019 2.342352 5753.7568 5820.8397
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 14429.24291 14469.53628 5.2936916 2000-2009 –> 2010-2019 11.591112 25999.9356 26078.5919
mod Atlantic (28.05,28.1] cstar_tref ~ temp + aou + phosphate 3576.19509 3600.92516 1.3461229 2000-2009 –> 2010-2019 2.423950 5865.8650 5913.7814
mod Atlantic (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate 3578.16459 3607.84067 1.3461031 2000-2009 –> 2010-2019 2.399424 5834.7369 5892.2365
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 16565.47810 16606.48318 5.9373916 2000-2009 –> 2010-2019 15.026502 27624.3214 27702.6255
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + phosphate 4571.14518 4603.14723 1.0725136 2000-2009 –> 2010-2019 2.399242 8942.7839 9005.7183
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + phosphate_star 4619.01954 4651.02159 1.0894141 2000-2009 –> 2010-2019 2.424440 9006.6318 9069.5662
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + silicate + phosphate 4540.57815 4577.91388 1.0611669 2000-2009 –> 2010-2019 2.269344 8674.4045 8747.8280
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4597.49244 4634.82817 1.0810756 2000-2009 –> 2010-2019 2.303205 8760.7360 8834.1595
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 16599.80439 16640.80947 5.9769288 2000-2009 –> 2010-2019 15.082290 27664.0889 27742.3930
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 16548.66533 16589.67040 5.9181221 2000-2009 –> 2010-2019 15.001083 27605.4472 27683.7513
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + phosphate_star 4638.28190 4670.28396 1.0962890 2000-2009 –> 2010-2019 2.405826 8976.5069 9039.4413
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 4639.23073 4676.56646 1.0959127 2000-2009 –> 2010-2019 2.375243 8919.6651 8993.0886
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 16573.86412 16614.86919 5.9470265 2000-2009 –> 2010-2019 15.056135 27639.4020 27717.7062
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + phosphate + phosphate_star 4603.89707 4635.89913 1.0840470 2000-2009 –> 2010-2019 2.413414 8980.6268 9043.5613
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + phosphate + phosphate_star 16601.10204 16636.24925 5.9807409 2000-2009 –> 2010-2019 15.084538 27662.8634 27729.9813
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4576.72809 4614.06382 1.0737693 2000-2009 –> 2010-2019 2.285164 8717.3691 8790.7926
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 16599.77887 16640.78395 5.9768993 2000-2009 –> 2010-2019 15.078915 27662.9441 27741.2482
obs Atlantic (28.1,28.15] cstar_tref ~ temp + aou + nitrate + phosphate_star 16592.55768 16627.70489 5.9708686 2000-2009 –> 2010-2019 15.074336 27654.2088 27721.3267
obs Atlantic (28.1,28.15] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 16568.00750 16609.01257 5.9402960 2000-2009 –> 2010-2019 15.043018 27631.4090 27709.7132
mod Atlantic (28.1,28.15] cstar_tref ~ temp + aou + silicate + phosphate 4670.36611 4702.36817 1.1078365 2000-2009 –> 2010-2019 2.394193 8962.8331 9025.7676
obs Atlantic (28.1,28.15] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 16579.49381 16614.64101 5.9558060 2000-2009 –> 2010-2019 15.054035 27639.3915 27706.5094
obs Atlantic (28.1,28.15] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 16556.64252 16597.64759 5.9272571 2000-2009 –> 2010-2019 15.024889 27618.3403 27696.6444
mod Atlantic (28.1,28.15] cstar_tref ~ temp + silicate + phosphate + phosphate_star 4717.54144 4749.54350 1.1250367 2000-2009 –> 2010-2019 2.416900 9020.9529 9083.8873
mod Atlantic (28.15,28.2] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 7256.96291 7298.75325 0.8460894 2000-2009 –> 2010-2019 1.897055 13650.7672 13732.3413
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate 7316.46458 7352.28488 0.8551308 2000-2009 –> 2010-2019 1.897336 13671.9090 13741.8297
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate + silicate 7219.54827 7261.33861 0.8406359 2000-2009 –> 2010-2019 1.879005 13560.9731 13642.5473
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate 5567.47742 5601.17922 5.0991468 2000-2009 –> 2010-2019 12.844995 10191.3844 10256.5847
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5527.44854 5561.15034 4.9883414 2000-2009 –> 2010-2019 12.652096 10137.1845 10202.3848
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + nitrate + phosphate_star 7631.45173 7667.27203 0.9029742 2000-2009 –> 2010-2019 1.963376 14062.0873 14132.0080
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 7594.17745 7635.96779 0.8968657 2000-2009 –> 2010-2019 1.950850 14000.4417 14082.0159
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 5513.16601 5546.86781 4.9493910 2000-2009 –> 2010-2019 12.662273 10131.4005 10196.6008
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5543.86633 5577.56813 5.0334938 2000-2009 –> 2010-2019 12.708934 10155.6286 10220.8290
mod Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + phosphate 7323.10701 7358.92730 0.8561131 2000-2009 –> 2010-2019 1.897440 13674.8905 13744.8112
mod Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + phosphate_star 7356.50303 7392.32332 0.8610687 2000-2009 –> 2010-2019 1.907948 13731.3871 13801.3078
mod Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate 7247.65747 7283.47777 0.8450218 2000-2009 –> 2010-2019 1.895644 13638.0463 13707.9670
obs Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 5528.25602 5561.95782 4.9905527 2000-2009 –> 2010-2019 12.995224 10195.8778 10261.0781
obs Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5499.84021 5533.54201 4.9133241 2000-2009 –> 2010-2019 12.819305 10150.9625 10216.1628
obs Atlantic (28.15,28.2] cstar_tref ~ temp + aou + silicate + phosphate_star 5527.95986 5556.84712 4.9952218 2000-2009 –> 2010-2019 12.922308 10180.6279 10236.5138
mod Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 7407.47426 7443.29455 0.8686878 2000-2009 –> 2010-2019 1.919546 13798.8359 13868.7567
obs Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 5506.31516 5540.01696 4.9308159 2000-2009 –> 2010-2019 12.866199 10162.3745 10227.5748
mod Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + silicate + phosphate_star 7611.08827 7646.90856 0.8998018 2000-2009 –> 2010-2019 1.953932 14015.9571 14085.8778
obs Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + silicate + phosphate_star 5514.09792 5542.98518 4.9573619 2000-2009 –> 2010-2019 12.983539 10183.2883 10239.1742
obs Atlantic (28.15,28.2] cstar_tref ~ temp + silicate + phosphate + phosphate_star 5543.74083 5572.62809 5.0386750 2000-2009 –> 2010-2019 12.987001 10199.9677 10255.8537
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate 10039.62543 10072.01769 0.6863566 2000-2009 –> 2010-2019 1.437711 15884.8908 15946.5511
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + phosphate 9977.38371 10016.25442 0.6817884 2000-2009 –> 2010-2019 1.431036 15810.1875 15884.1799
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 10756.76371 10796.25802 3.1851519 2000-2009 –> 2010-2019 6.499674 18334.6218 18411.0472
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + phosphate_star 9892.93610 9931.80681 0.6758297 2000-2009 –> 2010-2019 1.425095 15725.8587 15799.8510
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate 9929.50155 9968.37226 0.6784034 2000-2009 –> 2010-2019 1.413512 15664.1883 15738.1806
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 10755.31150 10794.80581 3.1840423 2000-2009 –> 2010-2019 6.500833 18335.1473 18411.5727
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 9288.18169 9333.53085 0.6345201 2000-2009 –> 2010-2019 1.357683 14940.4906 15026.8150
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + phosphate + phosphate_star 10758.16192 10792.01418 3.1877498 2000-2009 –> 2010-2019 6.503686 18335.2531 18400.7606
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + silicate + phosphate 10753.93079 10787.78306 3.1845154 2000-2009 –> 2010-2019 6.501618 18332.0376 18397.5451
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 10735.85917 10775.35348 3.1692168 2000-2009 –> 2010-2019 6.483022 18313.0919 18389.5173
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate 9958.87192 9997.74264 0.6804777 2000-2009 –> 2010-2019 1.430357 15796.0156 15870.0079
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate 9894.92764 9940.27680 0.6758291 2000-2009 –> 2010-2019 1.424855 15728.2047 15814.5291
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate 10755.06453 10794.55884 3.1838537 2000-2009 –> 2010-2019 6.502328 18336.3666 18412.7920
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 9889.74553 9935.09470 0.6754651 2000-2009 –> 2010-2019 1.424705 15724.4915 15810.8159
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + silicate 9850.51512 9895.86428 0.6727162 2000-2009 –> 2010-2019 1.405487 15570.7992 15657.1236
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + silicate + phosphate 10726.37846 10765.87277 3.1620162 2000-2009 –> 2010-2019 6.479103 18306.4721 18382.8974
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 10769.84987 10809.34418 3.1951679 2000-2009 –> 2010-2019 6.469630 18312.5663 18388.9916
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 9935.08462 9980.43379 0.6786561 2000-2009 –> 2010-2019 1.430794 15789.7142 15876.0386
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 10797.32945 10836.82376 3.2163032 2000-2009 –> 2010-2019 6.498080 18346.4946 18422.9200
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 10742.58234 10782.07665 3.1743330 2000-2009 –> 2010-2019 6.456131 18291.7661 18368.1915
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 22555.27957 22598.64690 5.4252360 2000-2009 –> 2010-2019 10.964321 38066.5807 38150.6460
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 29324.13446 29369.28884 5.5501496 2000-2009 –> 2010-2019 11.848474 50314.4888 50402.1788
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + phosphate + phosphate_star 22888.95434 22926.12634 5.6824027 2000-2009 –> 2010-2019 11.577338 38706.4631 38778.5191
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 22669.97139 22713.33872 5.5117674 2000-2009 –> 2010-2019 11.166020 38283.1366 38367.2019
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 29646.79685 29691.95123 5.7448976 2000-2009 –> 2010-2019 12.267006 50861.8579 50949.5479
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate 30104.74215 30143.44591 6.0343766 2000-2009 –> 2010-2019 12.517721 51279.4367 51354.5996
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 22662.21117 22705.57850 5.5058693 2000-2009 –> 2010-2019 11.322246 38415.3105 38499.3758
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 29783.41887 29828.57325 5.8294036 2000-2009 –> 2010-2019 12.150064 50796.5568 50884.2468
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 22764.84930 22808.21663 5.5843919 2000-2009 –> 2010-2019 11.503551 38604.6571 38688.7224
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 29836.64710 29881.80148 5.8626628 2000-2009 –> 2010-2019 12.193984 50860.6307 50948.3207
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 23241.09913 23284.46647 5.9636536 2000-2009 –> 2010-2019 12.154433 39302.9965 39387.0618
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 30024.92026 30070.07465 5.9818337 2000-2009 –> 2010-2019 12.461935 51198.3946 51286.0846
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 23121.31706 23164.68439 5.8659068 2000-2009 –> 2010-2019 11.962758 39107.5366 39191.6019
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 30109.20519 30154.35957 6.0359654 2000-2009 –> 2010-2019 12.652570 51416.8460 51504.5360
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 23320.00921 23363.37655 6.0289355 2000-2009 –> 2010-2019 12.302358 39447.5489 39531.6142
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 30220.72991 30265.88429 6.1083454 2000-2009 –> 2010-2019 12.761435 51563.7616 51651.4515
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 22865.50495 22908.87228 5.6624855 2000-2009 –> 2010-2019 11.673642 38781.6554 38865.7207
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 29965.17212 30010.32650 5.9437550 2000-2009 –> 2010-2019 12.319355 51034.0099 51121.6999
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 23388.70863 23432.07597 6.0863518 2000-2009 –> 2010-2019 12.396184 39544.8940 39628.9593
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 30269.58971 30314.74410 6.1403285 2000-2009 –> 2010-2019 12.808770 51627.4554 51715.1454
mod Indo-Pacific (26,26.5] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 16280.36951 16322.11616 4.0959361 2000-2009 –> 2010-2019 8.545107 30509.0751 30591.4229
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 16565.91500 16607.66166 4.3044759 2000-2009 –> 2010-2019 8.573762 30593.1355 30675.4833
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 15841.73982 15883.48648 3.7951042 2000-2009 –> 2010-2019 7.862101 29631.9788 29714.3266
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 16255.91535 16297.66201 4.0785535 2000-2009 –> 2010-2019 8.358496 30295.3058 30377.6536
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 17506.76733 17548.71084 4.6595489 2000-2009 –> 2010-2019 9.233250 32193.5594 32276.2600
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 17526.38377 17568.32729 4.6750301 2000-2009 –> 2010-2019 9.265900 32231.8797 32314.5803
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + phosphate 17789.54549 17825.49707 4.8894111 2000-2009 –> 2010-2019 9.753253 32781.3750 32852.2613
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 15464.84658 15506.59324 3.5543252 2000-2009 –> 2010-2019 7.633288 29269.4292 29351.7769
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 17368.87885 17410.82237 4.5521654 2000-2009 –> 2010-2019 9.229387 32167.3992 32250.0999
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 15612.63625 15654.38290 3.6468645 2000-2009 –> 2010-2019 7.894814 29615.3966 29697.7444
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 17440.14780 17482.09131 4.6073549 2000-2009 –> 2010-2019 9.320871 32277.2558 32359.9565
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 17529.63762 17571.58113 4.6776029 2000-2009 –> 2010-2019 9.273908 32241.0410 32323.7416
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate_star 17583.11248 17619.06406 4.7216865 2000-2009 –> 2010-2019 9.489679 32475.5853 32546.4716
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 17350.79960 17392.74312 4.5382706 2000-2009 –> 2010-2019 9.252292 32188.4427 32271.1433
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 15655.21200 15696.95866 3.6739679 2000-2009 –> 2010-2019 7.954211 29694.9457 29777.2935
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 17439.53819 17481.48170 4.6068800 2000-2009 –> 2010-2019 9.322158 32278.5122 32361.2128
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate_star 17438.68928 17474.64086 4.6077768 2000-2009 –> 2010-2019 9.342648 32296.3637 32367.2500
mod Indo-Pacific (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 16031.42294 16073.16960 3.9223864 2000-2009 –> 2010-2019 8.200674 30068.9255 30151.2733
mod Indo-Pacific (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 16096.43667 16138.18333 3.9669874 2000-2009 –> 2010-2019 8.255303 30145.3689 30227.7167
mod Indo-Pacific (26,26.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 16143.27791 16185.02457 3.9994357 2000-2009 –> 2010-2019 8.305307 30212.1557 30294.5035
obs Indo-Pacific (26.5,26.75] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 15829.91440 15864.87306 5.6806742 2000-2009 –> 2010-2019 11.062314 32418.6366 32488.9388
mod Indo-Pacific (26.5,26.75] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 14216.12700 14257.28204 3.5566606 2000-2009 –> 2010-2019 7.893795 27044.6432 27125.7350
mod Indo-Pacific (26.5,26.75] cstar_tref ~ aou + silicate + phosphate + phosphate_star 14257.94894 14293.22468 3.5862798 2000-2009 –> 2010-2019 8.202832 27361.6740 27431.1813
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 15674.26417 15715.04927 5.5047721 2000-2009 –> 2010-2019 10.532592 31901.5595 31983.5787
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 15699.07628 15739.86138 5.5320913 2000-2009 –> 2010-2019 10.618035 31987.7961 32069.8153
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + phosphate 15879.99086 15914.94952 5.7377161 2000-2009 –> 2010-2019 11.113669 32463.0638 32533.3660
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 14038.23946 14079.39450 3.4389175 2000-2009 –> 2010-2019 7.805857 26897.1632 26978.2551
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 15827.10813 15867.89323 5.6752293 2000-2009 –> 2010-2019 11.042136 32403.1805 32485.1997
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 14208.47680 14249.63184 3.5515150 2000-2009 –> 2010-2019 8.031757 27181.1306 27262.2225
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 15701.59051 15742.37561 5.5348671 2000-2009 –> 2010-2019 10.630085 32000.0466 32082.0658
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate_star 15722.16661 15757.12527 5.5598546 2000-2009 –> 2010-2019 10.767575 32135.3336 32205.6357
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 14721.63297 14762.78801 3.9137240 2000-2009 –> 2010-2019 8.527369 27824.5609 27905.6528
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 15707.90318 15748.68828 5.5418427 2000-2009 –> 2010-2019 10.681342 32052.6021 32134.6213
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 14233.66112 14274.81617 3.5684824 2000-2009 –> 2010-2019 8.057271 27214.7762 27295.8681
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 15863.40842 15904.19352 5.7164824 2000-2009 –> 2010-2019 11.113525 32469.4048 32551.4240
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 14386.57448 14427.72952 3.6732592 2000-2009 –> 2010-2019 8.104415 27310.3135 27391.4054
obs Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 15812.19158 15852.97669 5.6583640 2000-2009 –> 2010-2019 11.077797 32440.3122 32522.3314
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 14603.12240 14644.27744 3.8269233 2000-2009 –> 2010-2019 8.354404 27622.3450 27703.4369
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + silicate + phosphate 14432.04608 14467.32183 3.7064086 2000-2009 –> 2010-2019 8.457457 27663.2750 27732.7823
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 14645.40446 14686.55951 3.8576688 2000-2009 –> 2010-2019 8.401911 27681.0341 27762.1259
mod Indo-Pacific (26.75,27] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 18484.28674 18520.99061 3.8056019 2000-2009 –> 2010-2019 7.732989 36845.3213 36918.6207
obs Indo-Pacific (26.75,27] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 24828.55953 24866.35077 5.3081036 2000-2009 –> 2010-2019 10.349403 49418.8198 49494.4454
mod Indo-Pacific (26.75,27] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 17499.74502 17542.56621 3.2848369 2000-2009 –> 2010-2019 7.204989 35850.6393 35936.1552
obs Indo-Pacific (26.75,27] cstar_tref ~ aou + phosphate + phosphate_star 25138.21621 25169.70891 5.5180095 2000-2009 –> 2010-2019 10.737344 50007.3865 50070.4079
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + phosphate 17817.47507 17860.29626 3.4442669 2000-2009 –> 2010-2019 7.145303 35789.6738 35875.1897
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + phosphate 25005.08252 25049.17230 5.4246440 2000-2009 –> 2010-2019 10.563181 49751.9742 49840.2041
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 17993.53722 18036.35840 3.5359194 2000-2009 –> 2010-2019 7.316549 36105.8278 36191.3437
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate 25252.72656 25296.81633 5.5944168 2000-2009 –> 2010-2019 10.759211 50040.8722 50129.1021
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 18015.68447 18058.50566 3.5476199 2000-2009 –> 2010-2019 7.331995 36134.4945 36220.0105
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + phosphate_star 18504.22917 18540.93305 3.8169393 2000-2009 –> 2010-2019 7.624238 36660.8010 36734.1004
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 17937.18822 17980.00941 3.5063235 2000-2009 –> 2010-2019 7.309525 36088.6710 36174.1869
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 25137.16917 25181.25895 5.5145450 2000-2009 –> 2010-2019 10.601875 49802.9973 49891.2272
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + phosphate 25072.40842 25110.19965 5.4716444 2000-2009 –> 2010-2019 10.672306 49914.5697 49990.1954
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 17484.87212 17527.69330 3.2775575 2000-2009 –> 2010-2019 7.206920 35851.2180 35936.7339
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 25043.35456 25087.44434 5.4505409 2000-2009 –> 2010-2019 10.536222 49706.5598 49794.7897
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 17669.76556 17712.58675 3.3692092 2000-2009 –> 2010-2019 7.397753 36200.2360 36285.7519
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 25185.75645 25229.84623 5.5479880 2000-2009 –> 2010-2019 10.712672 49973.7291 50061.9590
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 17697.95223 17740.77342 3.3834047 2000-2009 –> 2010-2019 7.417121 36236.8702 36322.3861
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 25186.22836 25230.31814 5.5483138 2000-2009 –> 2010-2019 10.711644 49972.0794 50060.3093
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate_star 25184.58362 25222.37485 5.5485591 2000-2009 –> 2010-2019 10.713334 49970.6990 50046.3246
mod Indo-Pacific (27,27.25] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 21228.26707 21266.64400 2.6529547 2000-2009 –> 2010-2019 5.834829 43175.8123 43252.3271
mod Indo-Pacific (27,27.25] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 21179.55042 21224.32351 2.6378120 2000-2009 –> 2010-2019 5.649262 42660.4108 42749.6781
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 31773.62076 31819.61059 4.9221063 2000-2009 –> 2010-2019 9.675364 62935.0808 63027.0058
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 31979.42556 32025.41539 5.0191416 2000-2009 –> 2010-2019 9.880500 63376.1068 63468.0318
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate 31774.69004 31814.10989 4.9235396 2000-2009 –> 2010-2019 9.677859 62936.4870 63015.2798
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate_star 31983.89848 32023.31834 5.0222244 2000-2009 –> 2010-2019 9.902775 63419.7942 63498.5871
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 31907.84143 31953.83126 4.9851752 2000-2009 –> 2010-2019 9.805183 63215.1711 63307.0961
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + silicate + phosphate + phosphate_star 31913.82228 31953.24213 4.9889507 2000-2009 –> 2010-2019 9.830819 63266.4841 63345.2769
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate 21164.97129 21209.74437 2.6334750 2000-2009 –> 2010-2019 5.798465 43069.2180 43158.4854
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 21563.38632 21608.15941 2.7545996 2000-2009 –> 2010-2019 6.016046 43723.2312 43812.4985
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate 31773.86508 31819.85491 4.9222204 2000-2009 –> 2010-2019 9.674188 62932.4865 63024.4115
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 31950.60152 31996.59135 5.0054369 2000-2009 –> 2010-2019 9.816389 63238.2599 63330.1849
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 21591.10181 21635.87489 2.7632300 2000-2009 –> 2010-2019 6.032515 43771.3859 43860.6532
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 22058.02633 22102.79941 2.9127585 2000-2009 –> 2010-2019 6.154345 44165.8677 44255.1350
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 32041.42637 32087.41620 5.0487477 2000-2009 –> 2010-2019 9.901514 63419.6055 63511.5305
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 31905.86322 31951.85304 4.9842398 2000-2009 –> 2010-2019 9.783995 63169.1491 63261.0741
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + phosphate 21556.95332 21595.33025 2.7532218 2000-2009 –> 2010-2019 5.981097 43626.7052 43703.2201
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 21557.47588 21602.24897 2.7527627 2000-2009 –> 2010-2019 5.872473 43339.0387 43428.3060
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 21916.83292 21961.60601 2.8667085 2000-2009 –> 2010-2019 6.105066 44016.1908 44105.4582
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 21937.46261 21982.23569 2.8733911 2000-2009 –> 2010-2019 6.115507 44046.6952 44135.9625
obs Indo-Pacific (27.25,27.5] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 33852.66981 33892.75486 4.2813641 2000-2009 –> 2010-2019 8.645887 66601.3901 66681.3212
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 33849.98533 33896.75123 4.2796616 2000-2009 –> 2010-2019 8.587312 66452.2540 66545.5070
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 29341.56009 29388.06934 3.2029657 2000-2009 –> 2010-2019 6.269594 57036.2980 57129.0308
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 33544.96001 33591.72590 4.1702503 2000-2009 –> 2010-2019 8.409329 65965.6128 66058.8657
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + silicate + phosphate 33549.96846 33590.05351 4.1727325 2000-2009 –> 2010-2019 8.422833 65998.0117 66077.9428
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 33792.04120 33838.80709 4.2586587 2000-2009 –> 2010-2019 8.592685 66463.4015 66556.6544
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + silicate + phosphate + phosphate_star 33792.56393 33832.64897 4.2595710 2000-2009 –> 2010-2019 8.636812 66574.2153 66654.1464
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 28918.12764 28964.63689 3.0857152 2000-2009 –> 2010-2019 6.271170 57027.2440 57119.9768
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 29178.52604 29225.03529 3.1573024 2000-2009 –> 2010-2019 6.411829 57521.4673 57614.2002
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 33526.76645 33573.53234 4.1638134 2000-2009 –> 2010-2019 8.358688 65828.7794 65922.0323
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 33740.52276 33787.28865 4.2400715 2000-2009 –> 2010-2019 8.487118 66182.4314 66275.6843
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 29197.84906 29244.35831 3.1626802 2000-2009 –> 2010-2019 6.423254 57561.0244 57653.7573
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 28804.81135 28851.32059 3.0550720 2000-2009 –> 2010-2019 6.106201 56444.3181 56537.0510
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 33682.98797 33729.75386 4.2194095 2000-2009 –> 2010-2019 8.471189 66137.5010 66230.7539
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 33687.18679 33733.95268 4.2209140 2000-2009 –> 2010-2019 8.458616 66104.1633 66197.4162
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + phosphate 28919.70002 28959.56509 3.0866863 2000-2009 –> 2010-2019 6.336876 57246.1057 57325.5910
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 28595.00918 28641.51843 2.9991380 2000-2009 –> 2010-2019 6.131639 56521.4063 56614.1392
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 28875.87488 28922.38413 3.0742534 2000-2009 –> 2010-2019 6.294891 57104.7208 57197.4536
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 28877.26734 28923.77659 3.0746304 2000-2009 –> 2010-2019 6.296384 57109.8870 57202.6198
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + nitrate + silicate + phosphate_star 28875.64302 28915.50809 3.0747322 2000-2009 –> 2010-2019 6.297438 57109.4853 57188.9706
obs Indo-Pacific (27.5,27.75] cstar_tref ~ aou + silicate + phosphate + phosphate_star 38452.21569 38493.20772 4.0024898 2000-2009 –> 2010-2019 7.672271 72675.5632 72757.0364
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 38980.26054 39028.08457 4.1591639 2000-2009 –> 2010-2019 7.800018 73106.0374 73201.0896
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 32843.91798 32891.22247 3.1964150 2000-2009 –> 2010-2019 6.109572 65005.1526 65099.8708
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 38218.22189 38266.04591 3.9341341 2000-2009 –> 2010-2019 7.772628 73009.1035 73104.1557
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 33016.86953 33064.17401 3.2401729 2000-2009 –> 2010-2019 6.153681 65179.6586 65274.3768
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate_star 33282.55827 33323.10497 3.3090839 2000-2009 –> 2010-2019 6.349211 65992.9846 66074.1716
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 32970.24085 33017.54534 3.2283169 2000-2009 –> 2010-2019 6.137493 65113.8063 65208.5244
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 33209.74049 33257.04497 3.2896774 2000-2009 –> 2010-2019 6.293705 65767.8313 65862.5494
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 38006.59414 38054.41817 3.8738294 2000-2009 –> 2010-2019 7.457874 71934.5007 72029.5529
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 33269.18750 33316.49199 3.3050877 2000-2009 –> 2010-2019 6.322932 65886.5672 65981.2854
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 38104.27634 38152.10037 3.9015488 2000-2009 –> 2010-2019 7.518469 72147.0669 72242.1191
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 33270.60112 33317.90561 3.3054550 2000-2009 –> 2010-2019 6.311092 65835.6127 65930.3309
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 32959.12083 33006.42532 3.2254958 2000-2009 –> 2010-2019 6.188545 65339.7532 65434.4714
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 38069.63340 38117.45743 3.8916955 2000-2009 –> 2010-2019 7.503032 72092.9879 72188.0401
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate 33307.84714 33348.39384 3.3156693 2000-2009 –> 2010-2019 6.352102 66002.5643 66083.7513
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 33240.15672 33287.46120 3.2975531 2000-2009 –> 2010-2019 6.307907 65825.4291 65920.1473
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 38090.06258 38137.88661 3.8975031 2000-2009 –> 2010-2019 7.513014 72127.9523 72223.0044
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate_star 38381.68340 38422.67542 3.9819366 2000-2009 –> 2010-2019 7.607091 72451.0870 72532.5602
obs Indo-Pacific (27.5,27.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 38525.88675 38573.71078 4.0234836 2000-2009 –> 2010-2019 7.739225 72907.8304 73002.8826
obs Indo-Pacific (27.5,27.75] cstar_tref ~ temp + aou + silicate + phosphate 38526.73634 38567.72837 4.0243206 2000-2009 –> 2010-2019 7.747346 72931.3223 73012.7956
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + phosphate 11520.78014 11561.50633 2.4505414 2000-2009 –> 2010-2019 4.683657 22646.4613 22727.9558
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 11561.36187 11602.08807 2.4706328 2000-2009 –> 2010-2019 4.711674 22704.7575 22786.2521
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + phosphate 11848.68595 11883.59412 2.6187273 2000-2009 –> 2010-2019 4.900513 23080.1697 23150.0221
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 11604.42839 11645.15458 2.4921346 2000-2009 –> 2010-2019 4.773513 22837.0195 22918.5140
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 21254.95049 21298.68103 3.9097590 2000-2009 –> 2010-2019 7.000221 37652.9204 37739.1822
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 11644.60498 11685.33117 2.5123622 2000-2009 –> 2010-2019 4.801145 22893.3995 22974.8940
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 21233.74127 21277.47181 3.8989117 2000-2009 –> 2010-2019 6.990876 37634.8387 37721.1005
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 11586.61566 11627.34186 2.4832187 2000-2009 –> 2010-2019 4.729866 22742.5039 22823.9984
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 21231.52563 21275.25617 3.8977803 2000-2009 –> 2010-2019 6.998509 37650.8286 37737.0904
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 11661.49199 11702.21818 2.5209132 2000-2009 –> 2010-2019 4.811924 22915.1513 22996.6458
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 21233.17656 21276.90710 3.8986233 2000-2009 –> 2010-2019 6.990992 37635.1148 37721.3766
mod Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 11513.44008 11554.16628 2.4469249 2000-2009 –> 2010-2019 4.844112 22993.6135 23075.1080
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 21259.51364 21303.24418 3.9120967 2000-2009 –> 2010-2019 7.002366 37657.0831 37743.3449
mod Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 11563.91968 11604.64588 2.4719047 2000-2009 –> 2010-2019 4.877779 23062.1792 23143.6737
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 21240.95933 21284.68987 3.9026000 2000-2009 –> 2010-2019 6.994557 37642.0400 37728.3018
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + silicate + phosphate_star 21246.32707 21283.81039 3.9063683 2000-2009 –> 2010-2019 6.998351 37645.4621 37719.4008
mod Indo-Pacific (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 11598.12174 11638.84793 2.4889742 2000-2009 –> 2010-2019 4.900789 23108.7109 23190.2054
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 21241.17790 21284.90844 3.9027117 2000-2009 –> 2010-2019 6.995057 37643.0669 37729.3287
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate_star 21239.37911 21276.86243 3.9028146 2000-2009 –> 2010-2019 7.003546 37656.6886 37730.6272
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + silicate + phosphate + phosphate_star 21245.35163 21282.83495 3.9058692 2000-2009 –> 2010-2019 6.998252 37645.3194 37719.2581
obs Indo-Pacific (27.85,27.95] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 27296.80670 27342.72445 3.3063997 2000-2009 –> 2010-2019 6.050347 48699.5910 48790.2406
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 13704.99792 13747.43495 2.0927587 2000-2009 –> 2010-2019 4.110217 26397.0024 26481.4583
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 27323.40490 27369.32264 3.3148390 2000-2009 –> 2010-2019 6.059907 48729.7860 48820.4356
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate 13106.11907 13142.49367 1.9048956 2000-2009 –> 2010-2019 3.856962 25599.1502 25671.5409
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + phosphate 13107.97410 13150.41113 1.9048521 2000-2009 –> 2010-2019 3.852033 25588.0237 25672.4796
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 13107.94468 13150.38171 1.9048433 2000-2009 –> 2010-2019 3.853399 25592.2131 25676.6690
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + silicate 13091.07750 13133.51453 1.8997871 2000-2009 –> 2010-2019 3.820296 25488.6771 25573.1330
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + phosphate 13733.26276 13769.63736 2.1027631 2000-2009 –> 2010-2019 4.156079 26528.5846 26600.9754
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 13586.23131 13628.66834 2.0539565 2000-2009 –> 2010-2019 4.101062 26365.4440 26449.8999
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 27277.45145 27323.36920 3.3002719 2000-2009 –> 2010-2019 6.048749 48694.7649 48785.4145
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 13603.75618 13646.19321 2.0596365 2000-2009 –> 2010-2019 4.112272 26399.0975 26483.5533
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 27268.39418 27314.31192 3.2974084 2000-2009 –> 2010-2019 6.056987 48721.2129 48811.8625
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 13151.75839 13194.19542 1.9180401 2000-2009 –> 2010-2019 3.875683 25663.8420 25748.2978
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 13619.69210 13662.12913 2.0648151 2000-2009 –> 2010-2019 4.119666 26421.4817 26505.9376
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 27268.28157 27314.19932 3.2973728 2000-2009 –> 2010-2019 6.058625 48726.4395 48817.0891
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 27259.95654 27305.87429 3.2947429 2000-2009 –> 2010-2019 6.042897 48676.2348 48766.8844
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 27256.75857 27302.67631 3.2937333 2000-2009 –> 2010-2019 6.053439 48709.9812 48800.6308
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + silicate + phosphate 27282.42805 27321.78612 3.3024794 2000-2009 –> 2010-2019 6.051626 48699.8849 48777.5846
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + silicate + phosphate_star 27275.42185 27314.77992 3.3002626 2000-2009 –> 2010-2019 6.060418 48728.0816 48805.7812
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 27256.76703 27302.68478 3.2937359 2000-2009 –> 2010-2019 6.054713 48714.0460 48804.6956
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 46542.98458 46585.45913 3.4346369 2000-2009 –> 2010-2019 6.285942 83249.8448 83333.8055
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 46543.26415 46592.81779 3.4343000 2000-2009 –> 2010-2019 6.273212 83187.3249 83285.2791
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + phosphate + phosphate_star 46699.78326 46735.17872 3.4658736 2000-2009 –> 2010-2019 6.321106 83425.1154 83495.0827
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + silicate + phosphate + phosphate_star 46699.11311 46741.58766 3.4653461 2000-2009 –> 2010-2019 6.313329 83388.6291 83472.5898
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 46544.96598 46594.51962 3.4346333 2000-2009 –> 2010-2019 6.285911 83253.6824 83351.6366
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + aou + phosphate + phosphate_star 46687.25574 46729.73030 3.4630042 2000-2009 –> 2010-2019 6.318083 83413.7881 83497.7488
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate 17710.89336 17749.35368 1.7351736 2000-2009 –> 2010-2019 3.460959 30500.7868 30575.7693
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate + silicate 17582.77230 17627.64268 1.7102231 2000-2009 –> 2010-2019 3.389167 30195.6918 30283.1713
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate 18010.17871 18048.63903 1.7939513 2000-2009 –> 2010-2019 3.620740 31170.0041 31244.9865
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate_star 18012.70225 18051.16257 1.7944553 2000-2009 –> 2010-2019 3.625396 31187.2929 31262.2753
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 18005.23154 18050.10191 1.7925646 2000-2009 –> 2010-2019 3.589405 31059.5471 31147.0267
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 46783.46456 46833.01820 3.4816543 2000-2009 –> 2010-2019 6.339443 83526.1127 83624.0668
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 18008.72165 18053.59202 1.7932611 2000-2009 –> 2010-2019 3.595933 31084.1097 31171.5892
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 18228.52283 18273.39320 1.8376759 2000-2009 –> 2010-2019 3.679130 31442.3512 31529.8308
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + phosphate + phosphate_star 18019.25373 18057.71405 1.7957644 2000-2009 –> 2010-2019 3.628723 31201.0097 31275.9922
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 18015.04993 18059.92030 1.7945247 2000-2009 –> 2010-2019 3.598729 31095.9655 31183.4451
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate_star 18283.47112 18321.93144 1.8493617 2000-2009 –> 2010-2019 3.707926 31555.4562 31630.4387
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 46631.59623 46681.14987 3.4516389 2000-2009 –> 2010-2019 6.303719 83344.4977 83442.4519
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 46745.83762 46795.39127 3.4741935 2000-2009 –> 2010-2019 6.354325 83604.3364 83702.2906
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 46730.44196 46779.99560 3.4711453 2000-2009 –> 2010-2019 6.349423 83579.3591 83677.3133
obs Indo-Pacific (28.05,28.1] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 35852.68046 35893.93554 2.9592960 2000-2009 –> 2010-2019 5.687599 58714.8153 58794.8227
obs Indo-Pacific (28.05,28.1] cstar_tref ~ aou + phosphate + phosphate_star 35926.72497 35961.10420 2.9750594 2000-2009 –> 2010-2019 5.705021 58792.5939 58859.2667
obs Indo-Pacific (28.05,28.1] cstar_tref ~ aou + silicate + phosphate + phosphate_star 35903.26293 35944.51801 2.9697720 2000-2009 –> 2010-2019 5.698756 58767.7521 58847.7595
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 35869.47045 35917.60137 2.9623553 2000-2009 –> 2010-2019 5.702908 58775.8598 58869.2018
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + phosphate 35867.69481 35915.82573 2.9619878 2000-2009 –> 2010-2019 5.716629 58822.4473 58915.7893
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + silicate 9292.08095 9333.08331 1.4559539 2000-2009 –> 2010-2019 3.127496 16739.1890 16819.1228
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 9478.43771 9519.44008 1.5093922 2000-2009 –> 2010-2019 3.087833 16705.1349 16785.0688
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 35866.95224 35915.08316 2.9618342 2000-2009 –> 2010-2019 5.674587 58677.1752 58770.5171
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 9483.67700 9524.67936 1.5109226 2000-2009 –> 2010-2019 3.093785 16721.1310 16801.0649
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 35986.99416 36035.12508 2.9867775 2000-2009 –> 2010-2019 5.722342 58876.2015 58969.5435
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 9433.81914 9474.82151 1.4964218 2000-2009 –> 2010-2019 3.079263 16671.2220 16751.1559
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 35956.25666 36004.38758 2.9803707 2000-2009 –> 2010-2019 5.710077 58825.2422 58918.5842
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 9483.04454 9524.04691 1.5107378 2000-2009 –> 2010-2019 3.093600 16720.4982 16800.4320
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 35988.47885 36036.60978 2.9870873 2000-2009 –> 2010-2019 5.722289 58876.4350 58969.7769
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate_star 9488.61114 9523.75603 1.5129505 2000-2009 –> 2010-2019 3.095814 16724.0687 16792.5834
obs Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 36048.97483 36097.10575 2.9997385 2000-2009 –> 2010-2019 5.722271 58893.1403 58986.4823
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate 9886.60913 9921.75402 1.6340211 2000-2009 –> 2010-2019 3.237232 17171.1918 17239.7065
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate_star 9894.86812 9930.01300 1.6366335 2000-2009 –> 2010-2019 3.245767 17193.6318 17262.1465
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + nitrate + silicate + phosphate_star 9884.93656 9920.08144 1.6334925 2000-2009 –> 2010-2019 3.242162 17182.5913 17251.1060
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + silicate + phosphate + phosphate_star 9892.79925 9927.94413 1.6359787 2000-2009 –> 2010-2019 3.245089 17191.5070 17260.0217
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 18975.79797 19022.16522 1.3302136 2000-2009 –> 2010-2019 2.771952 36452.1636 36543.9984
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 37477.73512 37526.39792 2.7369406 2000-2009 –> 2010-2019 5.271086 67006.6912 67102.5732
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + phosphate + phosphate_star 37542.41545 37584.12643 2.7487831 2000-2009 –> 2010-2019 5.288058 67094.7839 67176.9684
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + silicate + phosphate 37512.28583 37553.99680 2.7434257 2000-2009 –> 2010-2019 5.279429 67048.4557 67130.6402
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + nitrate + silicate + phosphate 19064.84747 19104.59083 1.3411441 2000-2009 –> 2010-2019 2.804053 36681.8390 36760.5545
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 19066.66164 19113.02889 1.3411217 2000-2009 –> 2010-2019 2.797939 36644.8276 36736.6624
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 37592.49887 37641.16167 2.7573545 2000-2009 –> 2010-2019 5.326179 67292.2506 67388.1326
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + nitrate + phosphate 37635.44728 37684.11008 2.7650331 2000-2009 –> 2010-2019 5.339501 67362.7789 67458.6610
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + nitrate + silicate 18072.66784 18119.03509 1.2265025 2000-2009 –> 2010-2019 2.748924 36081.8003 36173.6351
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate 17661.43749 17707.80474 1.1819970 2000-2009 –> 2010-2019 2.557133 34675.0484 34766.8832
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate 37496.88634 37545.54914 2.7403367 2000-2009 –> 2010-2019 5.257629 66941.9942 67037.8763
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 17680.41773 17726.78498 1.1840152 2000-2009 –> 2010-2019 2.561046 34707.5072 34799.3419
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 37598.64071 37647.30351 2.7584513 2000-2009 –> 2010-2019 5.295739 67143.1721 67239.0541
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 18519.73338 18566.10063 1.2767894 2000-2009 –> 2010-2019 2.715768 35977.3530 36069.1878
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 17639.36493 17685.73218 1.1796544 2000-2009 –> 2010-2019 2.567004 34739.4964 34831.3312
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 37675.12833 37723.79113 2.7721466 2000-2009 –> 2010-2019 5.326713 67304.9395 67400.8215
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 17676.85825 17723.22550 1.1836364 2000-2009 –> 2010-2019 2.560400 34702.0443 34793.8790
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 37582.38415 37631.04695 2.7555492 2000-2009 –> 2010-2019 5.288808 67106.9456 67202.8276
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate_star 17855.92901 17895.67237 1.2030572 2000-2009 –> 2010-2019 2.593526 34976.0278 35054.7433
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate_star 37705.28355 37746.99452 2.7779244 2000-2009 –> 2010-2019 5.340483 67372.3503 67454.5349
obs Indo-Pacific (28.2,28.3] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 11303.26545 11338.57065 2.0288656 2000-2009 –> 2010-2019 4.638928 19553.7412 19621.7886
obs Indo-Pacific (28.2,28.3] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 11222.88822 11264.07762 1.9976334 2000-2009 –> 2010-2019 4.605714 19472.7321 19552.1207
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 11245.36549 11286.55488 2.0061074 2000-2009 –> 2010-2019 4.613670 19494.5213 19573.9099
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 31796.58556 31849.91205 0.6963416 2000-2009 –> 2010-2019 1.350617 54012.5489 54117.1160
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 34920.12127 34973.44776 0.7725701 2000-2009 –> 2010-2019 1.541881 60751.1981 60855.7652
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + phosphate + phosphate_star 11377.21838 11412.52358 2.0573195 2000-2009 –> 2010-2019 4.666518 19626.5465 19694.5939
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 11367.37528 11408.56468 2.0527361 2000-2009 –> 2010-2019 4.661799 19618.5239 19697.9125
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + nitrate + silicate + phosphate 32248.48400 32294.19242 0.7069331 2000-2009 –> 2010-2019 1.453814 57417.1019 57506.7308
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 32071.28062 32124.60710 0.7027324 2000-2009 –> 2010-2019 1.362645 54478.7301 54583.2972
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + nitrate + phosphate 11227.89375 11269.08315 1.9995174 2000-2009 –> 2010-2019 4.657690 19543.6384 19623.0270
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 11360.38156 11401.57096 2.0500342 2000-2009 –> 2010-2019 4.658267 19610.4278 19689.8165
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + silicate + phosphate 34855.19451 34908.52099 0.7709037 2000-2009 –> 2010-2019 1.601806 62405.2793 62509.8464
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 34867.37022 34920.69670 0.7712159 2000-2009 –> 2010-2019 1.602421 62425.5808 62530.1479
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 32060.56113 32113.88761 0.7024819 2000-2009 –> 2010-2019 1.441520 56995.5747 57100.1418
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 34864.99596 34918.32244 0.7711550 2000-2009 –> 2010-2019 1.602322 62422.1840 62526.7511
obs Indo-Pacific (28.2,28.3] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 11259.61554 11300.80494 2.0114982 2000-2009 –> 2010-2019 4.621486 19511.9913 19591.3799
mod Indo-Pacific (28.2,28.3] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 36026.75069 36080.07717 0.8015337 2000-2009 –> 2010-2019 1.639027 63753.1786 63857.7457
obs Indo-Pacific (28.2,28.3] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 11230.85911 11272.04851 2.0006344 2000-2009 –> 2010-2019 4.648599 19533.2748 19612.6634
obs Indo-Pacific (28.2,28.3] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 11373.66416 11414.85356 2.0551687 2000-2009 –> 2010-2019 4.664718 19625.4586 19704.8472
mod Indo-Pacific (28.2,28.3] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 36026.08018 36079.40667 0.8015158 2000-2009 –> 2010-2019 1.640259 63785.7971 63890.3643
mod Indo-Pacific (28.3,28.4] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 20728.40269 20775.60629 1.2626229 2000-2009 –> 2010-2019 2.189567 33633.2523 33725.7890
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 20401.59419 20448.79779 1.2301373 2000-2009 –> 2010-2019 2.158901 33325.2652 33417.8020
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 4498.88323 4533.57237 2.0518792 2000-2009 –> 2010-2019 4.700122 7776.7101 7843.0744
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 4482.50999 4517.19914 2.0359282 2000-2009 –> 2010-2019 4.670433 7753.2423 7819.6066
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + silicate + phosphate_star 4502.90300 4532.63656 2.0577751 2000-2009 –> 2010-2019 4.708333 7779.9217 7836.8054
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 4530.45108 4565.14023 2.0829865 2000-2009 –> 2010-2019 4.716424 7800.6308 7866.9952
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + silicate + phosphate + phosphate_star 4531.74127 4561.47482 2.0862557 2000-2009 –> 2010-2019 4.721165 7800.6831 7857.5669
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + nitrate + phosphate 20378.81926 20426.02286 1.2279048 2000-2009 –> 2010-2019 2.153956 33274.4190 33366.9558
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 20461.27020 20508.47380 1.2360062 2000-2009 –> 2010-2019 2.167022 33408.1825 33500.7193
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + silicate + phosphate 4507.41912 4542.10827 2.0602444 2000-2009 –> 2010-2019 4.693772 7777.6456 7844.0100
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4481.46639 4516.15554 2.0349157 2000-2009 –> 2010-2019 4.684628 7760.0497 7826.4141
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 20376.74701 20423.95062 1.2277019 2000-2009 –> 2010-2019 2.154822 33283.4221 33375.9589
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 20351.28586 20398.48947 1.2252113 2000-2009 –> 2010-2019 2.191460 33654.7194 33747.2561
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 4506.62716 4541.31630 2.0594669 2000-2009 –> 2010-2019 4.709468 7785.3591 7851.7234
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4501.21670 4535.90585 2.0541626 2000-2009 –> 2010-2019 4.688525 7771.8756 7838.2400
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + aou + nitrate + phosphate 20708.44295 20748.90318 1.2608156 2000-2009 –> 2010-2019 2.187408 33607.6477 33686.9649
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 20710.09364 20757.29724 1.2607804 2000-2009 –> 2010-2019 2.184684 33583.4031 33675.9398
obs Indo-Pacific (28.3,28.4] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 4479.04843 4513.73758 2.0325718 2000-2009 –> 2010-2019 4.717967 7775.8780 7842.2423
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 20711.21431 20751.67454 1.2610943 2000-2009 –> 2010-2019 2.188813 33622.0800 33701.3972
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 20712.82263 20760.02624 1.2610549 2000-2009 –> 2010-2019 2.186144 33598.4429 33690.9796
obs Indo-Pacific (28.4, Inf] cstar_tref ~ sal + temp + aou + phosphate -Inf -Inf 0.0000000 2000-2009 –> 2010-2019 1.656889 -Inf -Inf
obs Indo-Pacific (28.4, Inf] cstar_tref ~ sal + temp + nitrate + phosphate_star -Inf -Inf 0.0000000 2000-2009 –> 2010-2019 1.652663 -Inf -Inf
obs Indo-Pacific (28.4, Inf] cstar_tref ~ sal + temp + silicate + phosphate -Inf -Inf 0.0000000 2000-2009 –> 2010-2019 1.519547 -Inf -Inf
obs Indo-Pacific (28.4, Inf] cstar_tref ~ sal + temp + silicate + phosphate_star -Inf -Inf 0.0000000 2000-2009 –> 2010-2019 1.538645 -Inf -Inf
obs Indo-Pacific (28.4, Inf] cstar_tref ~ temp + aou + silicate + phosphate -Inf -Inf 0.0000000 2000-2009 –> 2010-2019 1.581647 -Inf -Inf
obs Indo-Pacific (28.4, Inf] cstar_tref ~ temp + aou + silicate + phosphate_star -Inf -Inf 0.0000000 2000-2009 –> 2010-2019 1.545706 -Inf -Inf
obs Indo-Pacific (28.4, Inf] cstar_tref ~ temp + nitrate + phosphate + phosphate_star -Inf -Inf 0.0000000 2000-2009 –> 2010-2019 1.853069 -Inf -Inf
obs Indo-Pacific (28.4, Inf] cstar_tref ~ temp + nitrate + silicate + phosphate -Inf -Inf 0.0000000 2000-2009 –> 2010-2019 1.567579 -Inf -Inf
obs Indo-Pacific (28.4, Inf] cstar_tref ~ temp + nitrate + silicate + phosphate_star -Inf -Inf 0.0000000 2000-2009 –> 2010-2019 1.674364 -Inf -Inf
obs Indo-Pacific (28.4, Inf] cstar_tref ~ temp + silicate + phosphate + phosphate_star -Inf -Inf 0.0000000 2000-2009 –> 2010-2019 1.565887 -Inf -Inf

7.3 Target variable coefficients

A data frame to map the target variable is prepared.

# create table with two era belonging to one eras
eras_forward <- GLODAP_glanced %>%
  arrange(era) %>% 
  group_by(basin, gamma_slab, model) %>% 
  mutate(eras = paste(era, lead(era), sep = " --> ")) %>% 
  ungroup() %>% 
  select(era, eras) %>% 
  unique()

eras_backward <- GLODAP_glanced %>%
  arrange(era) %>% 
  group_by(basin, gamma_slab, model) %>% 
  mutate(eras = paste(lag(era), era, sep = " --> ")) %>% 
  ungroup() %>% 
  select(era, eras) %>% 
  unique()

eras_era <- full_join(eras_backward, eras_forward) %>% 
  filter(str_detect(eras, "NA") == FALSE)

# extend best model selection from eras to era
lm_best <- full_join(
  lm_best_eras %>% select(data_source, basin, gamma_slab, model, eras),
  eras_era) %>% 
  drop_na() %>% 
  arrange(era, basin, gamma_slab, data_source)

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

rm(eras_era, eras_forward, eras_backward)

7.4 Cant coeffcients

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

# subtract coefficients of adjacent era  
lm_best_cant <- lm_best %>%
  arrange(era) %>%
  group_by(data_source, basin, gamma_slab, eras, model, term) %>%
  mutate(delta_coeff = estimate - lag(estimate)) %>%
  ungroup() %>%
  arrange(basin, gamma_slab, model, term, eras) %>%
  drop_na() %>%
  select(-c(era,estimate))

# pivot to wide format
lm_best_cant <- lm_best_cant %>%
  pivot_wider(values_from = delta_coeff,
              names_from = term,
              names_prefix = "delta_coeff_",
              values_fill = 0)

7.5 Write files

# create table of target varaible coefficients in wide format
lm_best_target <- lm_best %>% 
  pivot_wider(names_from = "term",
              names_prefix = "coeff_",
              values_from = "estimate",
              values_fill = 0
              )


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

lm_best_cant %>%
  write_csv(paste(path_version_data,
                  "lm_best_cant.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) %>%
  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
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
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) %>%
  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
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
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_eras$rmse,
        lm_all_eras$rmse_sum - lm_all_eras$rmse))

lm_all_eras %>%
  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
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_eras$rmse,
        lm_best_eras$rmse_sum - lm_best_eras$rmse))

lm_best_eras %>%
  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
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_cant %>% 
  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 phosphate phosphate_star sal silicate temp
Atlantic - mod
(-Inf,26] 6 6 6 8 10 8 6
(26,26.5] 9 7 7 7 7 6 7
(26.5,26.75] 6 7 7 7 6 6 10
(26.75,27] 8 5 8 6 5 10 6
(27,27.25] 9 7 6 4 6 10 6
(27.25,27.5] 6 8 4 7 9 8 7
(27.5,27.75] 7 5 8 5 10 7 6
(27.75,27.85] 5 4 4 6 10 10 6
(27.85,27.95] 7 4 4 5 10 9 7
(27.95,28.05] 4 2 4 8 10 5 10
(28.05,28.1] 6 2 6 6 8 5 10
(28.1,28.15] 5 2 6 7 8 6 10
(28.15,28.2] 6 10 2 6 5 5 9
(28.2, Inf] 9 10 3 4 10 3 5
total 93.00 79.00 75.00 86.00 114.00 98.00 105.00
Indo-Pacific - mod
(-Inf,26] 8 5 7 7 10 5 7
(26,26.5] 8 6 8 7 6 8 7
(26.5,26.75] 7 5 7 7 4 10 8
(26.75,27] 6 10 6 8 5 5 8
(27,27.25] 7 10 7 7 4 5 8
(27.25,27.5] 6 10 6 6 5 6 9
(27.5,27.75] 5 5 6 6 10 10 6
(27.75,27.85] 7 6 7 6 7 6 10
(27.85,27.95] 8 6 5 5 10 5 9
(27.95,28.05] 6 3 4 6 10 6 10
(28.05,28.1] 5 3 4 7 6 10 10
(28.1,28.2] 4 6 6 5 10 10 7
(28.2,28.3] 5 7 8 5 8 10 6
(28.3,28.4] 6 10 9 6 5 4 8
total 88.00 92.00 90.00 88.00 100.00 100.00 113.00
Atlantic - obs
(-Inf,26] 8 9 6 6 7 5 8
(26,26.5] 8 10 8 7 4 4 7
(26.5,26.75] 7 9 7 6 4 4 10
(26.75,27] 4 5 5 8 5 10 10
(27,27.25] 6 8 5 8 6 6 10
(27.25,27.5] 7 4 7 7 10 3 9
(27.5,27.75] 7 3 9 7 4 10 7
(27.75,27.85] 7 5 10 5 5 10 5
(27.85,27.95] 5 4 7 6 10 10 5
(27.95,28.05] 5 4 7 6 10 10 5
(28.05,28.1] 4 7 5 9 10 7 6
(28.1,28.15] 4 7 5 10 6 5 10
(28.15,28.2] 5 5 5 8 4 10 10
(28.2, Inf] 8 4 9 6 10 6 5
total 85.00 84.00 95.00 99.00 95.00 100.00 107.00
Indo-Pacific - obs
(-Inf,26] 8 6 6 6 10 5 8
(26,26.5] 5 5 5 7 10 5 10
(26.5,26.75] 6 7 7 6 8 4 9
(26.75,27] 7 8 7 6 3 6 8
(27,27.25] 6 4 6 7 10 10 4
(27.25,27.5] 6 5 8 7 9 8 4
(27.5,27.75] 7 4 7 6 7 9 7
(27.75,27.85] 5 5 5 8 4 10 10
(27.85,27.95] 8 4 7 7 4 10 8
(27.95,28.05] 9 6 9 8 3 6 4
(28.05,28.1] 8 5 8 6 6 7 6
(28.1,28.2] 6 4 7 6 10 8 6
(28.2,28.3] 10 6 10 8 5 4 5
(28.3,28.4] 6 3 7 8 9 10 5
(28.4, Inf] 3 4 6 6 4 7 10
total 100.00 76.00 105.00 102.00 102.00 109.00 104.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_eras %>% 
  ggplot(aes(rmse, AIC, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  facet_grid(data_source~basin)

Version Author Date
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_eras %>% 
  ggplot(aes(rmse, AIC, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  facet_grid(data_source~basin)

Version Author Date
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_eras %>% 
  ggplot(aes(AIC, BIC, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  facet_grid(data_source~basin)

Version Author Date
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
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 actual models

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

GLODAP_augmented_best <- left_join(
  lm_best_target %>% select(data_source, basin, gamma_slab, era, model),
  GLODAP_augmented
)
# calculate equal axis limits and binwidth
axis_lims <- GLODAP_augmented %>%
  summarise(
    max_value = max(
      c(max(.fitted, max(!!sym(params_local$MLR_target))))
      ),
    min_value = min(
      c(min(.fitted, min(!!sym(params_local$MLR_target))))
    )
  )

i_binwidth <- 1
# binwidth_value <- (axis_lims$max_value - axis_lims$min_value) / 40
axis_lims <- c(axis_lims$min_value, axis_lims$max_value)

GLODAP_augmented %>%
  ggplot(aes(!!sym(params_local$MLR_target), .fitted)) +
  geom_bin2d(binwidth = i_binwidth) +
  scale_fill_viridis_c() +
  geom_abline(slope = 1,
              col = 'red') +
  coord_equal(xlim = axis_lims,
              ylim = axis_lims) +
  labs(title = "All models") +
  facet_wrap(~data_source)

Version Author Date
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
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 
GLODAP_augmented_best %>%
  ggplot(aes(!!sym(params_local$MLR_target), .fitted)) +
  geom_bin2d(binwidth = i_binwidth) +
  scale_fill_viridis_c() +
  geom_abline(slope = 1,
              col = 'red') +
  coord_equal(xlim = axis_lims,
              ylim = axis_lims) +
  labs(title = "Selected models") +
  facet_wrap(~data_source)

Version Author Date
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 
rm(binwidth_value, axis_lims)

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 = "white") +
  scale_fill_viridis_c() +
  labs(
    title = paste(
      "Target variable:",
      params_local$MLR_target,
      "| Selected models",
      "| All domains"
    )
  ) +
  facet_wrap(~data_source)

Version Author Date
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 = "white") +
  scale_fill_viridis_c() +
  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
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, era) %>% 
  map(~ ggplot(data = .x, aes(lat, .resid)) +
      geom_bin2d(binwidth = i_binwidth) +
      geom_hline(yintercept = 0, col = "black") +
      scale_fill_viridis_c() +
      facet_grid(gamma_slab ~ data_source) +
      coord_cartesian(ylim = i_ylim) +
      labs(
        title = paste(
          "Target:",
          params_local$MLR_target,
          "| best models | basin:",
          unique(.x$basin),
          "| era:",
          unique(.x$era)
        )
      ))
[[1]]

Version Author Date
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
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
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
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
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
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(basin, gamma_slab, lat_grid) %>% 
  mutate(eras = paste(lag(era), era, sep = " --> "),
         .resid_mean_offset = .resid_mean - lag(.resid_mean)
         ) %>% 
  ungroup() %>% 
  select(-era) %>% 
  drop_na() %>% 
  filter(eras != paste(unique(lat_residual$era)[1],
                       unique(lat_residual$era)[3],
                       sep = " --> "))

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

Version Author Date
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

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
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
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  metR_0.9.0       scico_1.2.0      patchwork_1.1.1 
[13] collapse_1.5.0   forcats_0.5.0    stringr_1.4.0    dplyr_1.0.5     
[17] purrr_0.3.4      readr_1.4.0      tidyr_1.1.2      tibble_3.0.4    
[21] 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            plyr_1.8.6              
[55] grid_4.0.3               blob_1.2.1               parallel_4.0.3          
[58] promises_1.1.1           crayon_1.3.4             lattice_0.20-41         
[61] haven_2.3.1              hms_0.5.3                pillar_1.4.7            
[64] reprex_0.3.0             glue_1.4.2               evaluate_0.14           
[67] RcppArmadillo_0.10.1.2.0 data.table_1.13.2        modelr_0.1.8            
[70] vctrs_0.3.5              httpuv_1.5.4             cellranger_1.1.0        
[73] gtable_0.3.0             reshape_0.8.8            assertthat_0.2.1        
[76] xfun_0.18                openxlsx_4.2.3           RcppEigen_0.3.3.7.0     
[79] later_1.1.0.1            viridisLite_0.3.0        ellipsis_0.3.1          
[82] here_0.1