Last updated: 2021-06-10

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.


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.

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.

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.

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

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

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

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 ca07f9b. See the Past versions tab to see a history of the changes made to the R Markdown and HTML files.

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


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

Unstaged changes:
    Modified:   code/Workflowr_project_managment.R

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

1 Version ID

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

2 Required data

Required are:

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

3 Predictor combinations

Find all possible combinations of following considered predictor variables:

  • sal, temp, aou, 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 111.79255 125.31533 0.6311650 2000-2009 –> 2010-2019 4.681558 768.2401 800.9163
obs Atlantic (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 850.03516 869.88913 6.6769500 2000-2009 –> 2010-2019 14.634681 1276.2189 1310.6156
mod Atlantic (-Inf,26] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 61.26850 74.79128 0.3846114 2000-2009 –> 2010-2019 4.768848 735.7740 768.4502
mod Atlantic (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 83.57041 97.09319 0.4786069 2000-2009 –> 2010-2019 4.463756 736.3155 768.9916
obs Atlantic (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 842.61341 862.46738 6.4831722 2000-2009 –> 2010-2019 15.395101 1282.1603 1316.5570
mod Atlantic (-Inf,26] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 59.50004 73.02282 0.3780005 2000-2009 –> 2010-2019 4.250715 705.7199 738.3961
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate 852.25800 869.27569 6.7897804 2000-2009 –> 2010-2019 14.843101 1277.8507 1307.3336
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 850.64257 870.49655 6.6930634 2000-2009 –> 2010-2019 14.682434 1277.2945 1311.6913
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 850.68994 870.54391 6.6943216 2000-2009 –> 2010-2019 14.680514 1277.2949 1311.6917
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 56.01042 69.53320 0.3652871 2000-2009 –> 2010-2019 4.657767 725.6936 758.3697
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 844.98857 864.84255 6.5445665 2000-2009 –> 2010-2019 14.530605 1271.5913 1305.9881
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 81.87983 95.40261 0.4707397 2000-2009 –> 2010-2019 4.347945 728.3640 761.0401
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 82.34827 95.87105 0.4729066 2000-2009 –> 2010-2019 4.759587 751.7231 784.3993
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 111.22440 124.74718 0.6276591 2000-2009 –> 2010-2019 4.663799 766.8682 799.5444
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 850.90353 870.75751 6.6999980 2000-2009 –> 2010-2019 14.709771 1277.8564 1312.2532
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 55.84267 69.36545 0.3646868 2000-2009 –> 2010-2019 4.672811 726.3553 759.0314
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 81.69481 95.21759 0.4698866 2000-2009 –> 2010-2019 4.657279 745.7266 778.4028
obs Atlantic (-Inf,26] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 857.67229 877.52627 6.8823998 2000-2009 –> 2010-2019 15.373755 1291.5148 1325.9115
obs Atlantic (-Inf,26] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 856.75250 876.60647 6.8573248 2000-2009 –> 2010-2019 15.340516 1290.4815 1324.8782
obs Atlantic (-Inf,26] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 853.93169 873.78566 6.7809942 2000-2009 –> 2010-2019 15.228075 1287.1573 1321.5540
obs Atlantic (26,26.5] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 3998.91964 4025.46886 6.1223903 2000-2009 –> 2010-2019 11.591753 7253.7442 7305.8164
mod Atlantic (26,26.5] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 2097.65353 2125.14201 3.8933352 2000-2009 –> 2010-2019 8.363711 3459.0354 3510.6208
obs Atlantic (26,26.5] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 4000.91614 4031.89022 6.1223729 2000-2009 –> 2010-2019 11.569005 7253.4095 7314.1604
mod Atlantic (26,26.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 2052.02396 2079.51244 3.6635285 2000-2009 –> 2010-2019 7.669905 3362.7774 3414.3628
obs Atlantic (26,26.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 3867.15233 3898.12642 5.4934210 2000-2009 –> 2010-2019 10.951182 7121.7685 7182.5194
mod Atlantic (26,26.5] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2023.75666 2051.24514 3.5280207 2000-2009 –> 2010-2019 7.919193 3376.8797 3428.4651
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 2055.73443 2083.22291 3.6816979 2000-2009 –> 2010-2019 7.969326 3397.8330 3449.4184
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 3884.39674 3915.37082 5.5707272 2000-2009 –> 2010-2019 11.097407 7152.0636 7212.8145
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 2056.00902 2083.49750 3.6830461 2000-2009 –> 2010-2019 8.040453 3405.5658 3457.1512
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 3885.34256 3916.31665 5.5749986 2000-2009 –> 2010-2019 11.175524 7166.8135 7227.5643
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 2027.87079 2055.35928 3.5474269 2000-2009 –> 2010-2019 8.285851 3416.1559 3467.7414
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2028.56151 2056.04999 3.5506954 2000-2009 –> 2010-2019 8.406401 3428.1425 3479.7279
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 2057.01516 2084.50364 3.6879903 2000-2009 –> 2010-2019 8.117967 3414.2029 3465.7883
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 3895.58691 3926.56099 5.6214734 2000-2009 –> 2010-2019 11.282985 7188.3216 7249.0725
obs Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + phosphate 3997.46153 4024.01074 6.1151602 2000-2009 –> 2010-2019 11.676325 7269.5975 7321.6697
mod Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 2095.18887 2122.67735 3.8805619 2000-2009 –> 2010-2019 8.391056 3460.6983 3512.2837
obs Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 3999.46145 4030.43553 6.1151598 2000-2009 –> 2010-2019 11.597733 7258.7951 7319.5460
mod Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2091.50111 2118.98959 3.8615279 2000-2009 –> 2010-2019 8.408437 3460.7255 3512.3109
obs Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 3999.38713 4030.36121 6.1147915 2000-2009 –> 2010-2019 11.640200 7266.8146 7327.5655
obs Atlantic (26,26.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 3999.91682 4030.89091 6.1174169 2000-2009 –> 2010-2019 11.697842 7277.6485 7338.3994
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 2752.53775 2781.79787 4.1205501 2000-2009 –> 2010-2019 7.288724 4654.4547 4711.0524
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 4797.00490 4829.60199 5.2329338 2000-2009 –> 2010-2019 10.904280 8066.3097 8128.6295
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 2755.60407 2784.86419 4.1336505 2000-2009 –> 2010-2019 7.192359 4631.7116 4688.3093
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 4796.11774 4828.71483 5.2299511 2000-2009 –> 2010-2019 10.930477 8070.7186 8133.0384
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 2802.46494 2831.72506 4.3391182 2000-2009 –> 2010-2019 7.191744 4627.3739 4683.9715
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 4890.55083 4923.14792 5.5571843 2000-2009 –> 2010-2019 11.255054 8164.6706 8226.9905
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2808.65452 2837.91464 4.3670100 2000-2009 –> 2010-2019 7.165643 4619.5376 4676.1353
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 2762.21377 2791.47388 4.1620313 2000-2009 –> 2010-2019 7.087331 4605.5881 4662.1857
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 4793.22627 4825.82336 5.2202414 2000-2009 –> 2010-2019 10.919412 8067.5818 8129.9017
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 2809.63551 2838.89562 4.3714470 2000-2009 –> 2010-2019 7.166142 4619.4850 4676.0827
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + phosphate 4858.06057 4886.00093 5.4493497 2000-2009 –> 2010-2019 11.124851 8126.1212 8179.5382
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + phosphate_star 4855.11453 4883.05489 5.4390420 2000-2009 –> 2010-2019 11.139949 8127.7844 8181.2014
mod Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 2774.57390 2803.83401 4.2156274 2000-2009 –> 2010-2019 7.283764 4652.9404 4709.5381
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 4858.23236 4890.82944 5.4429508 2000-2009 –> 2010-2019 11.117730 8128.1616 8190.4814
mod Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2778.43789 2807.69800 4.2325236 2000-2009 –> 2010-2019 7.200564 4632.4587 4689.0563
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 4855.43817 4888.03526 5.4331854 2000-2009 –> 2010-2019 11.126275 8128.6919 8191.0118
mod Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 2787.65194 2812.73204 4.2819442 2000-2009 –> 2010-2019 7.346858 4663.2470 4711.7593
obs Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 4852.68738 4880.62774 5.4305644 2000-2009 –> 2010-2019 11.130788 8125.2335 8178.6505
mod Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 2782.21384 2811.47395 4.2491003 2000-2009 –> 2010-2019 7.169037 4624.2412 4680.8389
obs Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 4852.98340 4885.58049 5.4246207 2000-2009 –> 2010-2019 11.117623 8126.2213 8188.5411
mod Atlantic (26.75,27] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 4279.52569 4313.25050 2.4955272 2000-2009 –> 2010-2019 4.823243 7127.8356 7192.6359
mod Atlantic (26.75,27] cstar_tref ~ aou + silicate + phosphate + phosphate_star 4649.97600 4678.88298 3.0594876 2000-2009 –> 2010-2019 5.422155 7514.9458 7570.4889
mod Atlantic (26.75,27] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 4428.84026 4462.56508 2.7079231 2000-2009 –> 2010-2019 5.237372 7381.2093 7446.0096
mod Atlantic (26.75,27] cstar_tref ~ sal + aou + silicate + phosphate 4577.67404 4606.58103 2.9408392 2000-2009 –> 2010-2019 5.476729 7531.2269 7586.7700
mod Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate 4558.50583 4592.23064 2.9069807 2000-2009 –> 2010-2019 5.182370 7378.3500 7443.1503
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate 5688.70619 5722.88326 4.4422776 2000-2009 –> 2010-2019 10.820727 9714.0108 9779.1165
mod Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4578.57322 4612.29804 2.9390687 2000-2009 –> 2010-2019 5.288526 7438.5228 7503.3231
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5693.53293 5727.70999 4.4532869 2000-2009 –> 2010-2019 10.775605 9708.0010 9773.1066
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 5703.62424 5737.80130 4.4763925 2000-2009 –> 2010-2019 10.785002 9715.4309 9780.5366
mod Atlantic (26.75,27] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4585.84895 4619.57377 2.9507900 2000-2009 –> 2010-2019 5.333632 7463.4641 7528.2643
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5694.57579 5728.75286 4.4556692 2000-2009 –> 2010-2019 10.774861 9708.4374 9773.5431
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + silicate + phosphate_star 5713.16644 5742.46106 4.5029672 2000-2009 –> 2010-2019 10.908739 9741.7116 9797.5164
mod Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 4298.97034 4332.69515 2.5222141 2000-2009 –> 2010-2019 4.797921 7118.9891 7183.7894
obs Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 5694.26394 5728.44101 4.4549567 2000-2009 –> 2010-2019 10.923933 9736.8465 9801.9521
mod Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 4308.94620 4342.67101 2.5360161 2000-2009 –> 2010-2019 4.828661 7138.2493 7203.0496
obs Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5700.02706 5734.20412 4.4681425 2000-2009 –> 2010-2019 10.849642 9725.9179 9791.0236
mod Atlantic (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 4331.40826 4365.13307 2.5673704 2000-2009 –> 2010-2019 4.893221 7178.7147 7243.5150
obs Atlantic (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 5702.17323 5736.35029 4.4730629 2000-2009 –> 2010-2019 10.866115 9730.2816 9795.3873
obs Atlantic (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate_star 5707.01025 5736.30487 4.4887737 2000-2009 –> 2010-2019 10.912688 9739.0227 9794.8275
obs Atlantic (26.75,27] cstar_tref ~ temp + silicate + phosphate + phosphate_star 5702.37280 5731.66742 4.4781113 2000-2009 –> 2010-2019 10.929530 9739.6233 9795.4281
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 5425.31531 5461.29848 2.0647429 2000-2009 –> 2010-2019 4.383179 9412.0991 9481.5337
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 5572.87215 5608.85532 2.1890491 2000-2009 –> 2010-2019 4.667965 9677.3159 9746.7506
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate 5719.80926 5750.65198 2.3221085 2000-2009 –> 2010-2019 4.785560 9811.2515 9870.7669
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 5720.94014 5756.92331 2.3213090 2000-2009 –> 2010-2019 4.666181 9727.6560 9797.0906
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate 11086.21176 11124.32150 6.1629175 2000-2009 –> 2010-2019 12.060315 19058.2884 19132.2920
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 11060.38206 11098.49180 6.1165471 2000-2009 –> 2010-2019 11.899918 18983.8038 19057.8074
mod Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + silicate 5524.08632 5560.06950 2.1471437 2000-2009 –> 2010-2019 4.412109 9469.8497 9539.2843
mod Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate 5719.83849 5755.82166 2.3202960 2000-2009 –> 2010-2019 4.699869 9752.3799 9821.8146
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 11105.37749 11143.48723 6.1975515 2000-2009 –> 2010-2019 11.936637 19009.6438 19083.6474
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 11067.28206 11105.39180 6.1289000 2000-2009 –> 2010-2019 11.945348 19004.9157 19078.9193
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 11124.00078 11162.11052 6.2313918 2000-2009 –> 2010-2019 12.066062 19069.4296 19143.4332
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 11107.04111 11145.15085 6.2005670 2000-2009 –> 2010-2019 11.952942 19017.0711 19091.0747
obs Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + phosphate_star 11081.09727 11113.76276 6.1573077 2000-2009 –> 2010-2019 12.068437 19056.9698 19120.4015
mod Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate 5689.21924 5720.06196 2.2941352 2000-2009 –> 2010-2019 4.597101 9662.2332 9721.7486
mod Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 5690.67882 5726.66199 2.2936440 2000-2009 –> 2010-2019 4.554630 9633.3511 9702.7858
obs Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 11047.18520 11085.29494 6.0929905 2000-2009 –> 2010-2019 12.067911 19051.8066 19125.8102
mod Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5690.98777 5726.97094 2.2939248 2000-2009 –> 2010-2019 4.553735 9632.7452 9702.1799
obs Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 11006.34671 11044.45645 6.0206663 2000-2009 –> 2010-2019 11.863342 18955.1921 19029.1957
mod Atlantic (27,27.25] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 5606.54375 5642.52692 2.2184479 2000-2009 –> 2010-2019 4.480387 9549.9568 9619.3915
obs Atlantic (27,27.25] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 11017.54221 11055.65195 6.0404075 2000-2009 –> 2010-2019 11.917291 18980.9354 19054.9390
obs Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 13194.96510 13234.41223 5.8408099 2000-2009 –> 2010-2019 10.521615 21888.8228 21965.2972
mod Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate 11037.48008 11071.23072 3.5662526 2000-2009 –> 2010-2019 6.894390 20706.5385 20773.4008
mod Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 10896.45903 10935.83478 3.4439365 2000-2009 –> 2010-2019 6.734126 20525.2709 20603.2769
mod Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 10909.24250 10948.61825 3.4546965 2000-2009 –> 2010-2019 6.766724 20562.4253 20640.4313
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 13179.87557 13219.32270 5.8195599 2000-2009 –> 2010-2019 10.421769 21824.1170 21900.5914
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 11157.86355 11197.23930 3.6707774 2000-2009 –> 2010-2019 6.936125 20758.7537 20836.7597
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 13178.43971 13217.88684 5.8175419 2000-2009 –> 2010-2019 10.394743 21806.7163 21883.1907
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + silicate 11009.23901 11048.61476 3.5400325 2000-2009 –> 2010-2019 6.788235 20590.7353 20668.7414
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + phosphate 13246.02151 13279.83333 5.9161448 2000-2009 –> 2010-2019 10.533005 21897.5757 21963.1252
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + phosphate_star 13249.42486 13283.23669 5.9210103 2000-2009 –> 2010-2019 10.515232 21886.5767 21952.1262
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 13242.93846 13282.38559 5.9088855 2000-2009 –> 2010-2019 10.523794 21895.2540 21971.7284
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 10888.48721 10927.86296 3.4372436 2000-2009 –> 2010-2019 7.101934 20914.4292 20992.4352
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 13246.13405 13285.58117 5.9134482 2000-2009 –> 2010-2019 10.507561 21885.2167 21961.6911
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 11153.43861 11192.81436 3.6668159 2000-2009 –> 2010-2019 6.799423 20601.4405 20679.4465
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 13175.40231 13214.84944 5.8132753 2000-2009 –> 2010-2019 10.395609 21806.9626 21883.4370
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 10866.65552 10906.03127 3.4189807 2000-2009 –> 2010-2019 6.956576 20762.5651 20840.5711
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + phosphate + phosphate_star 13246.18760 13279.99942 5.9163821 2000-2009 –> 2010-2019 10.514419 21885.7720 21951.3214
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 10889.76604 10929.14179 3.4383164 2000-2009 –> 2010-2019 7.105437 20918.1506 20996.1566
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 13243.08854 13282.53567 5.9090997 2000-2009 –> 2010-2019 10.506660 21884.3689 21960.8433
mod Atlantic (27.25,27.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 11164.41152 11203.78727 3.6766475 2000-2009 –> 2010-2019 6.802581 20604.5570 20682.5630
obs Atlantic (27.5,27.75] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 14153.32599 14193.57125 5.0953584 2000-2009 –> 2010-2019 10.366142 25277.8361 25356.5618
obs Atlantic (27.5,27.75] cstar_tref ~ aou + silicate + phosphate + phosphate_star 14163.26072 14197.75665 5.1084812 2000-2009 –> 2010-2019 10.380409 25286.5532 25354.0324
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + nitrate + phosphate 18253.87260 18291.85778 2.1790248 2000-2009 –> 2010-2019 4.614400 33879.5952 33954.3360
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 17816.25038 17860.56643 2.0666127 2000-2009 –> 2010-2019 4.495252 33425.2451 33512.4427
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate 18565.63649 18603.62167 2.2624298 2000-2009 –> 2010-2019 4.704378 34209.5880 34284.3287
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 18532.82685 18577.14290 2.2529612 2000-2009 –> 2010-2019 4.687718 34158.8348 34246.0323
obs Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 14164.92916 14205.17442 5.1081162 2000-2009 –> 2010-2019 10.304252 25238.0041 25316.7298
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 18142.52596 18186.84200 2.1494698 2000-2009 –> 2010-2019 4.569169 33726.5850 33813.7825
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 18564.13903 18608.45508 2.2614766 2000-2009 –> 2010-2019 4.683767 34155.4334 34242.6310
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 14170.58753 14210.83279 5.1143492 2000-2009 –> 2010-2019 10.317845 25248.7665 25327.4922
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 18372.67580 18416.99185 2.2099063 2000-2009 –> 2010-2019 4.589834 33844.6674 33931.8649
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 14199.39596 14239.64122 5.1462015 2000-2009 –> 2010-2019 10.377317 25296.6644 25375.3901
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 18456.66054 18500.97658 2.2323811 2000-2009 –> 2010-2019 4.621863 33955.7428 34042.9403
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 18482.91224 18527.22829 2.2394529 2000-2009 –> 2010-2019 4.627015 33976.5604 34063.7579
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 18378.06876 18422.38481 2.2113427 2000-2009 –> 2010-2019 4.593036 33855.0745 33942.2720
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 14179.53324 14219.77850 5.1242190 2000-2009 –> 2010-2019 10.340302 25266.4246 25345.1503
obs Atlantic (27.5,27.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 14160.00203 14200.24729 5.1026948 2000-2009 –> 2010-2019 10.390464 25296.1138 25374.8395
obs Atlantic (27.5,27.75] cstar_tref ~ temp + aou + silicate + phosphate 14168.84335 14203.33929 5.1146312 2000-2009 –> 2010-2019 10.403006 25303.3681 25370.8473
obs Atlantic (27.5,27.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 14169.81681 14210.06207 5.1134998 2000-2009 –> 2010-2019 10.425263 25322.2543 25400.9800
obs Atlantic (27.5,27.75] cstar_tref ~ temp + silicate + phosphate + phosphate_star 14178.19629 14212.69222 5.1249512 2000-2009 –> 2010-2019 10.436884 25328.7488 25396.2280
obs Atlantic (27.75,27.85] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 9246.87643 9284.01745 5.3731962 2000-2009 –> 2010-2019 10.550341 16973.9946 17047.1021
obs Atlantic (27.75,27.85] cstar_tref ~ aou + silicate + phosphate + phosphate_star 9244.87751 9276.71267 5.3731982 2000-2009 –> 2010-2019 10.575164 16982.0387 17044.7023
mod Atlantic (27.75,27.85] cstar_tref ~ sal + aou + nitrate + silicate 5217.71566 5249.68629 1.3365542 2000-2009 –> 2010-2019 2.574047 9489.3023 9552.3167
obs Atlantic (27.75,27.85] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 9233.61663 9270.75765 5.3493248 2000-2009 –> 2010-2019 10.542566 16968.5516 17041.6591
mod Atlantic (27.75,27.85] cstar_tref ~ sal + aou + silicate + phosphate 5004.25167 5036.22230 1.2460951 2000-2009 –> 2010-2019 2.432525 9165.8581 9228.8725
obs Atlantic (27.75,27.85] cstar_tref ~ sal + aou + silicate + phosphate 9231.67711 9263.51227 5.3494334 2000-2009 –> 2010-2019 10.570502 16978.0684 17040.7320
obs Atlantic (27.75,27.85] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 9235.00875 9272.14977 5.3518260 2000-2009 –> 2010-2019 10.576448 16985.1130 17058.2205
mod Atlantic (27.75,27.85] cstar_tref ~ sal + nitrate + silicate + phosphate_star 5613.15984 5645.13046 1.5218382 2000-2009 –> 2010-2019 2.887483 10141.9357 10204.9500
mod Atlantic (27.75,27.85] cstar_tref ~ sal + silicate + phosphate + phosphate_star 5604.82687 5636.79749 1.5176806 2000-2009 –> 2010-2019 2.889200 10144.8051 10207.8195
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + silicate 5175.23214 5212.53120 1.3171772 2000-2009 –> 2010-2019 2.551227 9441.5491 9515.0658
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 4888.05129 4925.35035 1.1986664 2000-2009 –> 2010-2019 2.380380 9041.2598 9114.7766
obs Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 9232.58061 9269.72163 5.3474641 2000-2009 –> 2010-2019 10.557684 16975.7341 17048.8416
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5016.43796 5053.73702 1.2502692 2000-2009 –> 2010-2019 2.471106 9254.6588 9328.1756
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 5021.37629 5058.67536 1.2522979 2000-2009 –> 2010-2019 2.471753 9256.6404 9330.1572
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5010.40122 5047.70029 1.2477938 2000-2009 –> 2010-2019 2.464656 9240.1108 9313.6276
obs Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 9234.23148 9271.37250 5.3504293 2000-2009 –> 2010-2019 10.570165 16981.9796 17055.0871
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate_star 5097.71878 5129.68941 1.2849245 2000-2009 –> 2010-2019 2.589200 9506.4894 9569.5038
obs Atlantic (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 9246.92671 9284.06773 5.3732869 2000-2009 –> 2010-2019 10.557369 16977.4165 17050.5240
obs Atlantic (27.75,27.85] cstar_tref ~ temp + aou + silicate + phosphate 9244.92837 9276.76353 5.3732899 2000-2009 –> 2010-2019 10.583525 16986.0892 17048.7528
obs Atlantic (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 9247.27589 9284.41691 5.3739170 2000-2009 –> 2010-2019 10.565842 16981.5722 17054.6797
mod Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate 6319.08144 6350.96076 1.9806438 2000-2009 –> 2010-2019 4.899636 12825.4628 12888.3813
obs Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 19492.39980 19534.79034 5.3169996 2000-2009 –> 2010-2019 10.587780 36162.9183 36246.6237
mod Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 6319.59103 6356.78358 1.9796601 2000-2009 –> 2010-2019 4.821020 12757.6720 12831.0769
obs Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 18890.50061 18932.89116 4.8328203 2000-2009 –> 2010-2019 9.870016 35315.9706 35399.6760
mod Atlantic (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate 6095.43996 6127.31928 1.8383617 2000-2009 –> 2010-2019 4.572966 12431.6454 12494.5639
obs Atlantic (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate_star 18894.96175 18931.29651 4.8377762 2000-2009 –> 2010-2019 9.875246 35318.7264 35390.4739
obs Atlantic (27.85,27.95] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 19290.74840 19333.13895 5.1496116 2000-2009 –> 2010-2019 10.343369 35881.6844 35965.3898
obs Atlantic (27.85,27.95] cstar_tref ~ sal + silicate + phosphate + phosphate_star 19332.54446 19368.87921 5.1855123 2000-2009 –> 2010-2019 10.382794 35925.1476 35996.8951
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + silicate 6160.71854 6197.91108 1.8775498 2000-2009 –> 2010-2019 4.783787 12657.6793 12731.0842
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + phosphate 6554.12558 6586.00490 2.1420643 2000-2009 –> 2010-2019 4.893894 12906.7071 12969.6255
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 5797.85768 5835.05023 1.6636501 2000-2009 –> 2010-2019 4.398017 12135.8360 12209.2409
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 18728.82991 18771.22046 4.7104549 2000-2009 –> 2010-2019 9.716884 35121.1795 35204.8849
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5901.57928 5938.77182 1.7221748 2000-2009 –> 2010-2019 4.481362 12263.1246 12336.5295
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 18710.04999 18752.44054 4.6964431 2000-2009 –> 2010-2019 9.698111 35097.2556 35180.9610
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 18709.46811 18751.85865 4.6960096 2000-2009 –> 2010-2019 9.741353 35143.6753 35227.3807
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 5909.82263 5947.01517 1.7269134 2000-2009 –> 2010-2019 4.530353 12312.8636 12386.2685
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate 18733.58786 18769.92262 4.7155073 2000-2009 –> 2010-2019 9.762225 35167.2675 35239.0150
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5891.29916 5928.49170 1.7162835 2000-2009 –> 2010-2019 4.466938 12244.7662 12318.1711
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 18727.55711 18769.94766 4.7095039 2000-2009 –> 2010-2019 9.714778 35118.6589 35202.3643
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate_star 5991.79378 6023.67311 1.7759333 2000-2009 –> 2010-2019 4.808484 12597.7115 12660.6300
obs Atlantic (27.95,28.05] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 19947.78046 19990.29864 5.3979092 2000-2009 –> 2010-2019 12.105174 37715.2416 37798.9939
obs Atlantic (27.95,28.05] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 19532.94444 19575.46262 5.0601468 2000-2009 –> 2010-2019 11.622265 37183.4912 37267.2435
obs Atlantic (27.95,28.05] cstar_tref ~ sal + aou + silicate + phosphate_star 19619.10419 19655.54835 5.1301123 2000-2009 –> 2010-2019 11.706556 37279.3046 37351.0922
obs Atlantic (27.95,28.05] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 19669.03564 19711.55383 5.1685568 2000-2009 –> 2010-2019 11.775560 37356.0111 37439.7634
obs Atlantic (27.95,28.05] cstar_tref ~ sal + silicate + phosphate + phosphate_star 19797.79409 19834.23825 5.2749060 2000-2009 –> 2010-2019 11.895792 37493.9866 37565.7743
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate 9622.65320 9655.57288 3.5773782 2000-2009 –> 2010-2019 7.377761 17467.1436 17531.6222
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate_star 9652.48993 9685.40961 3.6074188 2000-2009 –> 2010-2019 7.406622 17496.0978 17560.5764
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 9486.18653 9524.59282 3.4412069 2000-2009 –> 2010-2019 7.240729 17332.0327 17407.2577
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 19651.32939 19693.84757 5.1543216 2000-2009 –> 2010-2019 11.753483 37331.9584 37415.7106
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 9497.95181 9536.35811 3.4525729 2000-2009 –> 2010-2019 7.249365 17341.7541 17416.9792
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 19596.75395 19639.27213 5.1106912 2000-2009 –> 2010-2019 11.686594 37258.5141 37342.2663
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate + phosphate_star 9717.12437 9750.04405 3.6733628 2000-2009 –> 2010-2019 7.459472 17550.9130 17615.3916
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 19680.70297 19723.22115 5.1779584 2000-2009 –> 2010-2019 11.893277 37454.5775 37538.3298
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 9472.18248 9510.58877 3.4277270 2000-2009 –> 2010-2019 7.177830 17280.7957 17356.0207
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + phosphate + phosphate_star 9647.68565 9680.60533 3.6025647 2000-2009 –> 2010-2019 7.401766 17491.2924 17555.7710
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + phosphate_star 9730.31137 9757.74444 3.6890316 2000-2009 –> 2010-2019 7.491122 17574.0797 17627.8119
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate 19783.77790 19820.22206 5.2634023 2000-2009 –> 2010-2019 11.983212 37559.2253 37631.0129
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 9498.19893 9536.60522 3.4528120 2000-2009 –> 2010-2019 7.249609 17342.0051 17417.2302
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 19635.57323 19678.09141 5.1416872 2000-2009 –> 2010-2019 11.732171 37309.1700 37392.9223
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate_star 9496.80661 9529.72629 3.4534001 2000-2009 –> 2010-2019 7.250205 17338.6187 17403.0973
obs Atlantic (28.05,28.1] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 14408.11038 14448.40374 5.2698011 2000-2009 –> 2010-2019 11.607201 26001.2442 26079.9006
obs Atlantic (28.05,28.1] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 14352.47349 14392.76686 5.2074175 2000-2009 –> 2010-2019 11.544615 25945.4939 26024.1503
obs Atlantic (28.05,28.1] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 14355.60130 14395.89467 5.2109049 2000-2009 –> 2010-2019 11.535357 25941.4827 26020.1391
obs Atlantic (28.05,28.1] cstar_tref ~ sal + silicate + phosphate + phosphate_star 14441.46649 14476.00366 5.3098325 2000-2009 –> 2010-2019 11.657792 26038.5038 26105.9235
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 14474.66161 14514.95498 5.3454051 2000-2009 –> 2010-2019 11.650603 26049.7310 26128.3873
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate 3563.04052 3592.71660 1.3363415 2000-2009 –> 2010-2019 2.330197 5730.9346 5788.4342
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate_star 3602.43304 3632.10912 1.3619161 2000-2009 –> 2010-2019 2.374152 5798.2916 5855.7912
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 3563.67700 3598.29909 1.3354649 2000-2009 –> 2010-2019 2.310712 5704.7291 5771.8120
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 3600.68877 3635.31087 1.3594643 2000-2009 –> 2010-2019 2.357800 5777.4469 5844.5298
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 14388.96406 14429.25743 5.2482492 2000-2009 –> 2010-2019 11.519671 25944.9869 26023.6433
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 14467.60022 14507.89359 5.3373320 2000-2009 –> 2010-2019 11.645811 26044.5144 26123.1708
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + phosphate_star 3648.90156 3678.57764 1.3927145 2000-2009 –> 2010-2019 2.415952 5861.2563 5918.7559
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + phosphate_star 14476.19064 14510.72781 5.3494443 2000-2009 –> 2010-2019 11.658351 26051.3452 26118.7650
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 3622.10518 3656.72728 1.3735477 2000-2009 –> 2010-2019 2.395962 5835.2322 5902.3151
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 14356.15895 14396.45231 5.2115269 2000-2009 –> 2010-2019 11.502607 25923.2801 26001.9365
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + phosphate + phosphate_star 3591.20186 3620.87794 1.3545750 2000-2009 –> 2010-2019 2.362408 5780.4082 5837.9079
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 3591.09990 3625.72200 1.3532055 2000-2009 –> 2010-2019 2.345058 5757.9167 5824.9996
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 14430.29566 14470.58903 5.2948846 2000-2009 –> 2010-2019 11.592305 26000.9883 26079.6447
mod Atlantic (28.05,28.1] cstar_tref ~ temp + aou + phosphate 3580.04281 3604.77288 1.3486177 2000-2009 –> 2010-2019 2.426445 5869.7127 5917.6291
mod Atlantic (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate 3582.02450 3611.70058 1.3486058 2000-2009 –> 2010-2019 2.401926 5838.5968 5896.0964
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 16565.62969 16606.63477 5.9375657 2000-2009 –> 2010-2019 15.026675 27624.4730 27702.7771
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + phosphate 4574.36068 4606.36274 1.0736404 2000-2009 –> 2010-2019 2.400369 8945.9994 9008.9338
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + phosphate_star 4622.25519 4654.25725 1.0905659 2000-2009 –> 2010-2019 2.425592 9009.8674 9072.8019
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + silicate + phosphate 4544.37672 4581.71245 1.0624842 2000-2009 –> 2010-2019 2.270661 8678.2031 8751.6266
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4601.20707 4638.54281 1.0823879 2000-2009 –> 2010-2019 2.304517 8764.4507 8837.8742
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 16600.09447 16641.09955 5.9772641 2000-2009 –> 2010-2019 15.082625 27664.3790 27742.6831
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 16548.81720 16589.82228 5.9182959 2000-2009 –> 2010-2019 15.001257 27605.5991 27683.9032
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + phosphate_star 4641.98037 4673.98243 1.0976139 2000-2009 –> 2010-2019 2.407151 8980.2053 9043.1398
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 4642.83673 4680.17246 1.0972040 2000-2009 –> 2010-2019 2.376534 8923.2711 8996.6946
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 16574.22931 16615.23438 5.9474464 2000-2009 –> 2010-2019 15.056555 27639.7672 27718.0713
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + phosphate + phosphate_star 4607.12094 4639.12300 1.0851890 2000-2009 –> 2010-2019 2.414556 8983.8507 9046.7851
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + phosphate + phosphate_star 16601.30522 16636.45243 5.9809759 2000-2009 –> 2010-2019 15.084773 27663.0666 27730.1844
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4580.49086 4617.82659 1.0750896 2000-2009 –> 2010-2019 2.286484 8721.1319 8794.5554
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 16600.04198 16641.04705 5.9772034 2000-2009 –> 2010-2019 15.079219 27663.2072 27741.5113
obs Atlantic (28.1,28.15] cstar_tref ~ temp + aou + nitrate + phosphate_star 16592.66559 16627.81279 5.9709932 2000-2009 –> 2010-2019 15.074461 27654.3167 27721.4346
obs Atlantic (28.1,28.15] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 16568.29139 16609.29646 5.9406221 2000-2009 –> 2010-2019 15.043344 27631.6929 27709.9971
mod Atlantic (28.1,28.15] cstar_tref ~ temp + aou + silicate + phosphate 4673.58506 4705.58712 1.1090017 2000-2009 –> 2010-2019 2.395358 8966.0521 9028.9865
obs Atlantic (28.1,28.15] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 16579.60071 16614.74792 5.9559291 2000-2009 –> 2010-2019 15.054158 27639.4984 27706.6163
obs Atlantic (28.1,28.15] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 16556.91458 16597.91965 5.9275689 2000-2009 –> 2010-2019 15.025201 27618.6124 27696.9165
mod Atlantic (28.1,28.15] cstar_tref ~ temp + silicate + phosphate + phosphate_star 4720.70230 4752.70436 1.1261987 2000-2009 –> 2010-2019 2.418062 9024.1137 9087.0482
mod Atlantic (28.15,28.2] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 7261.85653 7303.64687 0.8468053 2000-2009 –> 2010-2019 1.897770 13655.6608 13737.2350
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate 7321.50353 7357.32382 0.8558759 2000-2009 –> 2010-2019 1.898081 13676.9479 13746.8687
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate + silicate 7224.32626 7266.11660 0.8413304 2000-2009 –> 2010-2019 1.879699 13565.7511 13647.3252
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate 5567.37007 5601.07187 5.0988463 2000-2009 –> 2010-2019 12.844694 10191.2771 10256.4774
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5527.36129 5561.06309 4.9881025 2000-2009 –> 2010-2019 12.651857 10137.0972 10202.2975
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + nitrate + phosphate_star 7635.27741 7671.09771 0.9035714 2000-2009 –> 2010-2019 1.963973 14065.9129 14135.8337
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 7598.36760 7640.15794 0.8975155 2000-2009 –> 2010-2019 1.951499 14004.6319 14086.2060
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 5512.94465 5546.64645 4.9487898 2000-2009 –> 2010-2019 12.661672 10131.1791 10196.3794
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5543.75425 5577.45605 5.0331842 2000-2009 –> 2010-2019 12.708624 10155.5166 10220.7169
mod Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + phosphate 7327.79014 7363.61044 0.8568063 2000-2009 –> 2010-2019 1.898133 13679.5737 13749.4944
mod Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + phosphate_star 7361.15132 7396.97162 0.8617608 2000-2009 –> 2010-2019 1.908640 13736.0354 13805.9561
mod Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate 7252.13969 7287.95999 0.8456767 2000-2009 –> 2010-2019 1.896299 13642.5285 13712.4492
obs Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 5528.11872 5561.82052 4.9901766 2000-2009 –> 2010-2019 12.994848 10195.7405 10260.9408
obs Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5499.73324 5533.43504 4.9130356 2000-2009 –> 2010-2019 12.819017 10150.8555 10216.0558
obs Atlantic (28.15,28.2] cstar_tref ~ temp + aou + silicate + phosphate_star 5527.89451 5556.78177 4.9950426 2000-2009 –> 2010-2019 12.922129 10180.5625 10236.4485
mod Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 7412.04696 7447.86726 0.8693746 2000-2009 –> 2010-2019 1.920233 13803.4086 13873.3294
obs Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 5506.18327 5539.88507 4.9304590 2000-2009 –> 2010-2019 12.865842 10162.2426 10227.4429
mod Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + silicate + phosphate_star 7615.40588 7651.22618 0.9004735 2000-2009 –> 2010-2019 1.954604 14020.2747 14090.1954
obs Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + silicate + phosphate_star 5513.90221 5542.78947 4.9568294 2000-2009 –> 2010-2019 12.983007 10183.0926 10238.9785
obs Atlantic (28.15,28.2] cstar_tref ~ temp + silicate + phosphate + phosphate_star 5543.64806 5572.53532 5.0384185 2000-2009 –> 2010-2019 12.986745 10199.8750 10255.7609
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate 10052.47882 10084.87109 0.6872743 2000-2009 –> 2010-2019 1.438629 15897.7442 15959.4045
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + phosphate 9990.02714 10028.89786 0.6826851 2000-2009 –> 2010-2019 1.431933 15822.8310 15896.8233
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 10756.58577 10796.08008 3.1850159 2000-2009 –> 2010-2019 6.499538 18334.4439 18410.8692
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + phosphate_star 9905.41181 9944.28252 0.6767067 2000-2009 –> 2010-2019 1.425972 15738.3344 15812.3268
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate 9942.73034 9981.60106 0.6793369 2000-2009 –> 2010-2019 1.414445 15677.4171 15751.4094
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 10755.16827 10794.66258 3.1839329 2000-2009 –> 2010-2019 6.500724 18335.0041 18411.4295
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 9301.16501 9346.51418 0.6353770 2000-2009 –> 2010-2019 1.358540 14953.4740 15039.7984
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + phosphate + phosphate_star 10757.91153 10791.76380 3.1875583 2000-2009 –> 2010-2019 6.503495 18335.0027 18400.5102
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + silicate + phosphate 10753.76031 10787.61258 3.1843852 2000-2009 –> 2010-2019 6.501487 18331.8672 18397.3746
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 10735.72182 10775.21613 3.1691124 2000-2009 –> 2010-2019 6.482917 18312.9546 18389.3799
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate 9971.57593 10010.44664 0.6813770 2000-2009 –> 2010-2019 1.431256 15808.7196 15882.7119
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate 9907.41443 9952.76360 0.6767069 2000-2009 –> 2010-2019 1.425733 15740.6915 15827.0159
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate 10754.81559 10794.30990 3.1836635 2000-2009 –> 2010-2019 6.502138 18336.1177 18412.5431
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 9902.22649 9947.57565 0.6763420 2000-2009 –> 2010-2019 1.425582 15736.9725 15823.2969
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + silicate 9863.59462 9908.94378 0.6736314 2000-2009 –> 2010-2019 1.406402 15583.8787 15670.2031
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + silicate + phosphate 10726.16902 10765.66334 3.1618573 2000-2009 –> 2010-2019 6.478944 18306.2626 18382.6880
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 10769.32443 10808.81874 3.1947652 2000-2009 –> 2010-2019 6.469227 18312.0408 18388.4662
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 9947.49754 9992.84671 0.6795324 2000-2009 –> 2010-2019 1.431670 15802.1271 15888.4515
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 10796.89428 10836.38859 3.2159675 2000-2009 –> 2010-2019 6.497744 18346.0594 18422.4848
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 10742.15715 10781.65146 3.1740092 2000-2009 –> 2010-2019 6.455807 18291.3409 18367.7663
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 22563.23032 22606.59765 5.4311906 2000-2009 –> 2010-2019 10.970276 38074.5314 38158.5967
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 29326.59297 29371.74735 5.5516082 2000-2009 –> 2010-2019 11.849933 50316.9473 50404.6373
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + phosphate + phosphate_star 22894.48305 22931.65505 5.6867388 2000-2009 –> 2010-2019 11.581674 38711.9918 38784.0478
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 22673.87250 22717.23984 5.5147348 2000-2009 –> 2010-2019 11.168988 38287.0378 38371.1031
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 29644.96191 29690.11630 5.7437710 2000-2009 –> 2010-2019 12.265880 50860.0230 50947.7130
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate 30107.46736 30146.17112 6.0361346 2000-2009 –> 2010-2019 12.519479 51282.1619 51357.3248
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 22671.54490 22714.91223 5.5129641 2000-2009 –> 2010-2019 11.329341 38424.6443 38508.7096
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 29785.72275 29830.87713 5.8308393 2000-2009 –> 2010-2019 12.151500 50798.8607 50886.5507
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 22774.07227 22817.43961 5.5915025 2000-2009 –> 2010-2019 11.510661 38613.8801 38697.9454
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 29839.02203 29884.17641 5.8641512 2000-2009 –> 2010-2019 12.195472 50863.0056 50950.6956
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 23245.40852 23288.77586 5.9672004 2000-2009 –> 2010-2019 12.157979 39307.3059 39391.3712
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 30026.66347 30071.81785 5.9829484 2000-2009 –> 2010-2019 12.463050 51200.1378 51287.8278
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 23124.38721 23167.75455 5.8683921 2000-2009 –> 2010-2019 11.965243 39110.6067 39194.6720
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 30107.76147 30152.91585 6.0350341 2000-2009 –> 2010-2019 12.651639 51415.4023 51503.0922
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 23322.43681 23365.80415 6.0309552 2000-2009 –> 2010-2019 12.304377 39449.9765 39534.0418
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 30219.46209 30264.61647 6.1075178 2000-2009 –> 2010-2019 12.760607 51562.4937 51650.1837
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 22874.51920 22917.88653 5.6695322 2000-2009 –> 2010-2019 11.680688 38790.6696 38874.7350
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 29967.58748 30012.74186 5.9452896 2000-2009 –> 2010-2019 12.320889 51036.4253 51124.1153
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 23391.39430 23434.76163 6.0886075 2000-2009 –> 2010-2019 12.398439 39547.5797 39631.6450
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 30268.51393 30313.66831 6.1396225 2000-2009 –> 2010-2019 12.808064 51626.3796 51714.0696
mod Indo-Pacific (26,26.5] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 16275.97559 16317.72225 4.0928074 2000-2009 –> 2010-2019 8.541978 30504.6812 30587.0290
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 16566.47455 16608.22120 4.3048948 2000-2009 –> 2010-2019 8.574181 30593.6950 30676.0428
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 15842.99745 15884.74410 3.7959344 2000-2009 –> 2010-2019 7.862931 29633.2365 29715.5843
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 16255.29399 16297.04064 4.0781128 2000-2009 –> 2010-2019 8.358055 30294.6844 30377.0322
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 17497.72499 17539.66851 4.6524300 2000-2009 –> 2010-2019 9.226131 32184.5171 32267.2177
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 17517.20075 17559.14427 4.6677765 2000-2009 –> 2010-2019 9.258646 32222.6966 32305.3973
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + phosphate 17782.20341 17818.15499 4.8833448 2000-2009 –> 2010-2019 9.747186 32774.0329 32844.9192
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 15464.90958 15506.65623 3.5543641 2000-2009 –> 2010-2019 7.633327 29269.4921 29351.8399
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 17360.72309 17402.66660 4.5458920 2000-2009 –> 2010-2019 9.223113 32159.2435 32241.9441
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 15611.93516 15653.68181 3.6464199 2000-2009 –> 2010-2019 7.894369 29614.6955 29697.0433
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 17431.99101 17473.93453 4.6010047 2000-2009 –> 2010-2019 9.314520 32269.0991 32351.7997
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 17520.46072 17562.40423 4.6703502 2000-2009 –> 2010-2019 9.266656 32231.8641 32314.5647
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate_star 17574.84750 17610.79909 4.7150924 2000-2009 –> 2010-2019 9.483085 32467.3204 32538.2066
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 17342.57955 17384.52307 4.5319671 2000-2009 –> 2010-2019 9.245988 32180.2227 32262.9233
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 15654.38453 15696.13119 3.6734392 2000-2009 –> 2010-2019 7.953682 29694.1182 29776.4660
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 17431.40655 17473.35007 4.6005500 2000-2009 –> 2010-2019 9.315828 32270.3805 32353.0812
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate_star 17430.48139 17466.43297 4.6013862 2000-2009 –> 2010-2019 9.336258 32288.1558 32359.0421
mod Indo-Pacific (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 16024.93558 16066.68224 3.9179635 2000-2009 –> 2010-2019 8.196251 30062.4381 30144.7859
mod Indo-Pacific (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 16089.48659 16131.23325 3.9621954 2000-2009 –> 2010-2019 8.250511 30138.4188 30220.7666
mod Indo-Pacific (26,26.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 16136.21577 16177.96243 3.9945266 2000-2009 –> 2010-2019 8.300398 30205.0936 30287.4414
obs Indo-Pacific (26.5,26.75] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 15824.03239 15858.99105 5.6740114 2000-2009 –> 2010-2019 11.055651 32412.7546 32483.0568
mod Indo-Pacific (26.5,26.75] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 14208.57004 14249.72508 3.5515776 2000-2009 –> 2010-2019 7.888712 27037.0862 27118.1781
mod Indo-Pacific (26.5,26.75] cstar_tref ~ aou + silicate + phosphate + phosphate_star 14250.26232 14285.53807 3.5810666 2000-2009 –> 2010-2019 8.197619 27353.9874 27423.4947
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 15670.07047 15710.85557 5.5001680 2000-2009 –> 2010-2019 10.527988 31897.3658 31979.3850
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 15694.59086 15735.37596 5.5271426 2000-2009 –> 2010-2019 10.613086 31983.3107 32065.3299
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + phosphate 15875.63667 15910.59533 5.7327336 2000-2009 –> 2010-2019 11.108687 32458.7096 32529.0118
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 14031.18381 14072.33885 3.4343286 2000-2009 –> 2010-2019 7.801268 26890.1076 26971.1995
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 15822.59308 15863.37819 5.6701191 2000-2009 –> 2010-2019 11.037026 32398.6655 32480.6847
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 14201.34897 14242.50401 3.5467274 2000-2009 –> 2010-2019 8.026970 27174.0028 27255.0946
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 15697.08466 15737.86977 5.5298934 2000-2009 –> 2010-2019 10.625112 31995.5407 32077.5599
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate_star 15717.72507 15752.68373 5.5549298 2000-2009 –> 2010-2019 10.762650 32130.8920 32201.1942
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 14715.42309 14756.57813 3.9091272 2000-2009 –> 2010-2019 8.522772 27818.3510 27899.4429
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 15703.38911 15744.17421 5.5368537 2000-2009 –> 2010-2019 10.676353 32048.0880 32130.1072
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 14226.55109 14267.70613 3.5636840 2000-2009 –> 2010-2019 8.052473 27207.6662 27288.7581
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 15858.67865 15899.46375 5.7110903 2000-2009 –> 2010-2019 11.108133 32464.6750 32546.6942
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 14378.52266 14419.67770 3.6676662 2000-2009 –> 2010-2019 8.098821 27302.2617 27383.3536
obs Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 15805.71039 15846.49549 5.6510517 2000-2009 –> 2010-2019 11.070485 32433.8310 32515.8502
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 14595.04410 14636.19914 3.8210771 2000-2009 –> 2010-2019 8.348558 27614.2667 27695.3586
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + silicate + phosphate 14423.92704 14459.20279 3.7007180 2000-2009 –> 2010-2019 8.451766 27655.1560 27724.6633
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 14637.33554 14678.49058 3.8517825 2000-2009 –> 2010-2019 8.396024 27672.9651 27754.0570
mod Indo-Pacific (26.75,27] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 18480.99799 18517.70187 3.8037355 2000-2009 –> 2010-2019 7.731123 36842.0326 36915.3319
obs Indo-Pacific (26.75,27] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 24819.67159 24857.46283 5.3022360 2000-2009 –> 2010-2019 10.343536 49409.9319 49485.5575
mod Indo-Pacific (26.75,27] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 17495.06703 17537.88822 3.2825455 2000-2009 –> 2010-2019 7.202698 35845.9613 35931.4773
obs Indo-Pacific (26.75,27] cstar_tref ~ aou + phosphate + phosphate_star 25128.71418 25160.20688 5.5114887 2000-2009 –> 2010-2019 10.730823 49997.8845 50060.9059
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + phosphate 17817.70822 17860.52940 3.4443867 2000-2009 –> 2010-2019 7.145423 35789.9069 35875.4229
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + phosphate 24996.71825 25040.80803 5.4190007 2000-2009 –> 2010-2019 10.557538 49743.6100 49831.8399
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 17992.99435 18035.81554 3.5356331 2000-2009 –> 2010-2019 7.316263 36105.2849 36190.8009
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate 25243.97505 25288.06483 5.5883276 2000-2009 –> 2010-2019 10.753122 50032.1207 50120.3506
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 18015.04027 18057.86145 3.5472790 2000-2009 –> 2010-2019 7.331654 36133.8503 36219.3663
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + phosphate_star 18501.70569 18538.40957 3.8155028 2000-2009 –> 2010-2019 7.622801 36658.2776 36731.5769
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 17932.53807 17975.35925 3.5038923 2000-2009 –> 2010-2019 7.307094 36084.0208 36169.5367
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 25128.79302 25172.88280 5.5088000 2000-2009 –> 2010-2019 10.596130 49794.6211 49882.8510
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + phosphate 25063.12877 25100.92001 5.4653296 2000-2009 –> 2010-2019 10.665991 49905.2901 49980.9157
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 17478.54933 17521.37052 3.2744678 2000-2009 –> 2010-2019 7.203831 35844.8952 35930.4111
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 25034.29919 25078.38897 5.4444024 2000-2009 –> 2010-2019 10.530084 49697.5045 49785.7344
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 17662.47642 17705.29761 3.3655479 2000-2009 –> 2010-2019 7.394092 36192.9468 36278.4628
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 25176.59912 25220.68889 5.5416695 2000-2009 –> 2010-2019 10.706353 49964.5717 50052.8016
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 17690.56111 17733.38230 3.3796765 2000-2009 –> 2010-2019 7.413393 36229.4790 36314.9950
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 25177.12770 25221.21747 5.5420340 2000-2009 –> 2010-2019 10.705364 49962.9788 50051.2087
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate_star 25175.45895 25213.25019 5.5422624 2000-2009 –> 2010-2019 10.707037 49961.5743 50037.1999
mod Indo-Pacific (27,27.25] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 21220.76300 21259.13993 2.6507087 2000-2009 –> 2010-2019 5.832583 43168.3082 43244.8231
mod Indo-Pacific (27,27.25] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 21172.63174 21217.40482 2.6357529 2000-2009 –> 2010-2019 5.647203 42653.4921 42742.7595
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 31764.35787 31810.34770 4.9177833 2000-2009 –> 2010-2019 9.671041 62925.8179 63017.7429
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 31970.12292 32016.11275 5.0147145 2000-2009 –> 2010-2019 9.876072 63366.8042 63458.7292
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate 31765.40910 31804.82895 4.9192069 2000-2009 –> 2010-2019 9.673526 62927.2060 63005.9989
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate_star 31974.65167 32014.07152 5.0178211 2000-2009 –> 2010-2019 9.898371 63410.5474 63489.3403
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 31898.60899 31944.59881 4.9808112 2000-2009 –> 2010-2019 9.800819 63205.9386 63297.8636
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + silicate + phosphate + phosphate_star 31904.65299 31944.07284 4.9846132 2000-2009 –> 2010-2019 9.826481 63257.3148 63336.1076
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate 21163.86479 21208.63787 2.6331462 2000-2009 –> 2010-2019 5.798136 43068.1115 43157.3789
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 21562.48834 21607.26142 2.7543205 2000-2009 –> 2010-2019 6.015767 43722.3332 43811.6005
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate 31764.56927 31810.55910 4.9178820 2000-2009 –> 2010-2019 9.669849 62923.1907 63015.1157
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 31941.29515 31987.28498 5.0010201 2000-2009 –> 2010-2019 9.811972 63228.9536 63320.8785
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 21590.14629 21634.91937 2.7629320 2000-2009 –> 2010-2019 6.032217 43770.4304 43859.6977
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 22055.03799 22099.81108 2.9117762 2000-2009 –> 2010-2019 6.153363 44162.8794 44252.1467
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 32032.08335 32078.07317 5.0442752 2000-2009 –> 2010-2019 9.897042 63410.2625 63502.1875
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 31896.64606 31942.63589 4.9798839 2000-2009 –> 2010-2019 9.779639 63159.9319 63251.8569
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + phosphate 21550.64729 21589.02422 2.7512629 2000-2009 –> 2010-2019 5.979138 43620.3992 43696.9141
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 21551.26518 21596.03827 2.7508337 2000-2009 –> 2010-2019 5.870544 43332.8280 43422.0953
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 21910.91738 21955.69047 2.8647951 2000-2009 –> 2010-2019 6.103153 44010.2753 44099.5426
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 21931.56431 21976.33740 2.8714789 2000-2009 –> 2010-2019 6.113595 44040.7969 44130.0642
obs Indo-Pacific (27.25,27.5] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 33846.05936 33886.14441 4.2789618 2000-2009 –> 2010-2019 8.643485 66594.7796 66674.7107
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 33843.22439 33889.99028 4.2772056 2000-2009 –> 2010-2019 8.584856 66445.4931 66538.7460
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 29324.56917 29371.07842 3.1981761 2000-2009 –> 2010-2019 6.264805 57019.3071 57112.0399
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 33536.31232 33583.07821 4.1671895 2000-2009 –> 2010-2019 8.406268 65956.9651 66050.2180
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + silicate + phosphate 33541.25540 33581.34045 4.1696468 2000-2009 –> 2010-2019 8.419747 65989.2987 66069.2297
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 33783.53074 33830.29663 4.2555826 2000-2009 –> 2010-2019 8.589609 66454.8910 66548.1439
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + silicate + phosphate + phosphate_star 33784.10899 33824.19403 4.2565143 2000-2009 –> 2010-2019 8.633755 66565.7604 66645.6915
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 28898.54459 28945.05384 3.0803977 2000-2009 –> 2010-2019 6.265853 57007.6609 57100.3938
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 29159.37615 29205.88540 3.1519817 2000-2009 –> 2010-2019 6.406509 57502.3174 57595.0503
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 33518.04888 33564.81477 4.1607327 2000-2009 –> 2010-2019 8.355607 65820.0619 65913.3148
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 33731.67680 33778.44269 4.2368881 2000-2009 –> 2010-2019 8.483935 66173.5855 66266.8384
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 29178.73796 29225.24721 3.1573613 2000-2009 –> 2010-2019 6.417935 57541.9133 57634.6462
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 28786.35539 28832.86464 3.0501100 2000-2009 –> 2010-2019 6.101239 56425.8621 56518.5950
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 33674.41187 33721.17776 4.2163383 2000-2009 –> 2010-2019 8.468118 66128.9249 66222.1778
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 33678.45040 33725.21629 4.2177842 2000-2009 –> 2010-2019 8.455486 66095.4269 66188.6798
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + phosphate 28900.30646 28940.17153 3.0814184 2000-2009 –> 2010-2019 6.331608 57226.7121 57306.1974
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 28576.97249 28623.48174 2.9943774 2000-2009 –> 2010-2019 6.126879 56503.3696 56596.1025
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 28858.38188 28904.89112 3.0695206 2000-2009 –> 2010-2019 6.290158 57087.2278 57179.9606
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 28859.86591 28906.37516 3.0699218 2000-2009 –> 2010-2019 6.291675 57092.4856 57185.2184
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + nitrate + silicate + phosphate_star 28858.32268 28898.18775 3.0700453 2000-2009 –> 2010-2019 6.292751 57092.1650 57171.6503
obs Indo-Pacific (27.5,27.75] cstar_tref ~ aou + silicate + phosphate + phosphate_star 38445.88375 38486.87577 4.0006404 2000-2009 –> 2010-2019 7.670422 72669.2312 72750.7045
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 38974.58814 39022.41217 4.1574421 2000-2009 –> 2010-2019 7.798297 73100.3650 73195.4172
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 32819.85423 32867.15872 3.1903737 2000-2009 –> 2010-2019 6.103531 64981.0889 65075.8070
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 38212.11688 38259.94090 3.9323813 2000-2009 –> 2010-2019 7.770875 73002.9985 73098.0507
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 32993.24714 33040.55163 3.2341611 2000-2009 –> 2010-2019 6.147669 65156.0362 65250.7544
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate_star 33258.83569 33299.38239 3.3029183 2000-2009 –> 2010-2019 6.343046 65969.2620 66050.4490
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 32946.54900 32993.85349 3.2223095 2000-2009 –> 2010-2019 6.131486 65090.1144 65184.8326
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 33185.94790 33233.25239 3.2835298 2000-2009 –> 2010-2019 6.287558 65744.0387 65838.7568
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 38000.13768 38047.96171 3.8720042 2000-2009 –> 2010-2019 7.456049 71928.0443 72023.0964
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 33245.56266 33292.86715 3.2989548 2000-2009 –> 2010-2019 6.316799 65862.9424 65957.6606
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 38097.92787 38145.75189 3.8997413 2000-2009 –> 2010-2019 7.516661 72140.7184 72235.7706
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 33248.16512 33295.46960 3.2996299 2000-2009 –> 2010-2019 6.305267 65813.1767 65907.8948
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 32935.16777 32982.47226 3.2194276 2000-2009 –> 2010-2019 6.182476 65315.8001 65410.5183
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 38063.27337 38111.09739 3.8898893 2000-2009 –> 2010-2019 7.501226 72086.6279 72181.6800
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate 33285.02865 33325.57535 3.3097267 2000-2009 –> 2010-2019 6.346160 65979.7458 66060.9328
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 33216.52878 33263.83327 3.2914334 2000-2009 –> 2010-2019 6.301788 65801.8012 65896.5193
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 38083.71166 38131.53569 3.8956967 2000-2009 –> 2010-2019 7.511208 72121.6013 72216.6535
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate_star 38375.06459 38416.05661 3.9800133 2000-2009 –> 2010-2019 7.605168 72444.4682 72525.9414
obs Indo-Pacific (27.5,27.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 38519.65812 38567.48215 4.0216547 2000-2009 –> 2010-2019 7.737396 72901.6018 72996.6540
obs Indo-Pacific (27.5,27.75] cstar_tref ~ temp + aou + silicate + phosphate 38520.50395 38561.49597 4.0224902 2000-2009 –> 2010-2019 7.745516 72925.0899 73006.5632
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + phosphate 11509.57595 11550.30215 2.4450232 2000-2009 –> 2010-2019 4.678139 22635.2571 22716.7516
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 11550.32163 11591.04783 2.4651507 2000-2009 –> 2010-2019 4.706192 22693.7173 22775.2118
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + phosphate 11838.20834 11873.11651 2.6132123 2000-2009 –> 2010-2019 4.894998 23069.6921 23139.5445
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 11593.61759 11634.34379 2.4867196 2000-2009 –> 2010-2019 4.768098 22826.2087 22907.7032
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 21251.81100 21295.54154 3.9081514 2000-2009 –> 2010-2019 6.998613 37649.7809 37736.0427
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 11633.92235 11674.64855 2.5069679 2000-2009 –> 2010-2019 4.795751 22882.7169 22964.2114
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 21230.64967 21274.38021 3.8973331 2000-2009 –> 2010-2019 6.989298 37631.7471 37718.0089
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 11575.62452 11616.35072 2.4777331 2000-2009 –> 2010-2019 4.724381 22731.5128 22813.0073
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 21228.45479 21272.18534 3.8962127 2000-2009 –> 2010-2019 6.996942 37647.7577 37734.0195
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 11650.83682 11691.56302 2.5155144 2000-2009 –> 2010-2019 4.806526 22904.4961 22985.9907
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 21230.08918 21273.81972 3.8970470 2000-2009 –> 2010-2019 6.989416 37632.0275 37718.2893
mod Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 11502.53246 11543.25865 2.4415606 2000-2009 –> 2010-2019 4.838748 22982.7059 23064.2004
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 21256.50152 21300.23206 3.9105534 2000-2009 –> 2010-2019 7.000823 37654.0710 37740.3328
mod Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 11553.17829 11593.90449 2.4665680 2000-2009 –> 2010-2019 4.872442 23051.4378 23132.9323
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 21237.99052 21281.72106 3.9010826 2000-2009 –> 2010-2019 6.993039 37639.0712 37725.3330
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + silicate + phosphate_star 21243.34910 21280.83242 3.9048447 2000-2009 –> 2010-2019 6.996828 37642.4841 37716.4228
mod Indo-Pacific (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 11587.43439 11628.16058 2.4836278 2000-2009 –> 2010-2019 4.895442 23098.0236 23179.5181
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 21238.20869 21281.93923 3.9011941 2000-2009 –> 2010-2019 6.993539 37640.0977 37726.3595
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate_star 21236.41100 21273.89433 3.9012974 2000-2009 –> 2010-2019 7.002029 37653.7205 37727.6591
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + silicate + phosphate + phosphate_star 21242.37748 21279.86080 3.9043478 2000-2009 –> 2010-2019 6.996731 37642.3453 37716.2840
obs Indo-Pacific (27.85,27.95] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 27294.96577 27340.88351 3.3058164 2000-2009 –> 2010-2019 6.049764 48697.7500 48788.3997
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 13698.59450 13741.03153 2.0906480 2000-2009 –> 2010-2019 4.108106 26390.5990 26475.0549
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 27321.42077 27367.33852 3.3142088 2000-2009 –> 2010-2019 6.059277 48727.8019 48818.4515
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate 13092.94489 13129.31949 1.9009452 2000-2009 –> 2010-2019 3.853012 25585.9760 25658.3667
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + phosphate 13094.70843 13137.14546 1.9008744 2000-2009 –> 2010-2019 3.848055 25574.7580 25659.2139
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 13094.84547 13137.28250 1.9009154 2000-2009 –> 2010-2019 3.849471 25579.1139 25663.5698
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + silicate 13077.65023 13120.08726 1.8957716 2000-2009 –> 2010-2019 3.816281 25475.2499 25559.7057
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + phosphate 13722.22570 13758.60030 2.0991092 2000-2009 –> 2010-2019 4.152425 26517.5476 26589.9383
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 13573.89263 13616.32966 2.0499668 2000-2009 –> 2010-2019 4.097072 26353.1054 26437.5612
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 27275.45314 27321.37088 3.2996399 2000-2009 –> 2010-2019 6.048117 48692.7666 48783.4162
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 13591.65527 13634.09230 2.0557128 2000-2009 –> 2010-2019 4.108348 26386.9965 26471.4524
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 27266.37515 27312.29290 3.2967704 2000-2009 –> 2010-2019 6.056349 48719.1939 48809.8435
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 13138.94162 13181.37865 1.9141702 2000-2009 –> 2010-2019 3.871813 25651.0252 25735.4811
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 13607.70375 13650.14078 2.0609181 2000-2009 –> 2010-2019 4.115769 26409.4934 26493.9493
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 27266.26143 27312.17918 3.2967344 2000-2009 –> 2010-2019 6.057987 48724.4194 48815.0690
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 27258.08837 27304.00612 3.2941531 2000-2009 –> 2010-2019 6.042307 48674.3666 48765.0162
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 27254.88293 27300.80068 3.2931412 2000-2009 –> 2010-2019 6.052847 48708.1056 48798.7552
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + silicate + phosphate 27280.53256 27319.89063 3.3018795 2000-2009 –> 2010-2019 6.051026 48697.9894 48775.6891
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + silicate + phosphate_star 27273.52053 27312.87860 3.2996612 2000-2009 –> 2010-2019 6.059817 48726.1803 48803.8799
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 27254.89129 27300.80904 3.2931439 2000-2009 –> 2010-2019 6.054121 48712.1703 48802.8199
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 46540.70687 46583.18142 3.4341909 2000-2009 –> 2010-2019 6.285496 83247.5671 83331.5278
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 46541.03440 46590.58805 3.4338635 2000-2009 –> 2010-2019 6.272776 83185.0952 83283.0493
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + phosphate + phosphate_star 46697.52297 46732.91843 3.4654270 2000-2009 –> 2010-2019 6.320659 83422.8551 83492.8224
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + silicate + phosphate + phosphate_star 46696.91204 46739.38659 3.4649112 2000-2009 –> 2010-2019 6.312894 83386.4281 83470.3888
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 46542.68744 46592.24108 3.4341871 2000-2009 –> 2010-2019 6.285465 83251.4039 83349.3580
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + aou + phosphate + phosphate_star 46685.01552 46727.49008 3.4625619 2000-2009 –> 2010-2019 6.317641 83411.5479 83495.5086
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate 17706.38782 17744.84814 1.7343036 2000-2009 –> 2010-2019 3.460089 30496.2813 30571.2637
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate + silicate 17579.56591 17624.43629 1.7096128 2000-2009 –> 2010-2019 3.388556 30192.4854 30279.9649
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate 18006.95463 18045.41495 1.7933076 2000-2009 –> 2010-2019 3.620097 31166.7800 31241.7625
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate_star 18009.51847 18047.97879 1.7938195 2000-2009 –> 2010-2019 3.624760 31184.1091 31259.0916
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 18002.30708 18047.17745 1.7919812 2000-2009 –> 2010-2019 3.588822 31056.6226 31144.1022
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 46781.07214 46830.62579 3.4811795 2000-2009 –> 2010-2019 6.338969 83523.7203 83621.6744
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 18005.81989 18050.69026 1.7926820 2000-2009 –> 2010-2019 3.595354 31081.2079 31168.6875
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 18227.09863 18271.96900 1.8373846 2000-2009 –> 2010-2019 3.678838 31440.9270 31528.4066
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + phosphate + phosphate_star 18016.11292 18054.57324 1.7951367 2000-2009 –> 2010-2019 3.628096 31197.8689 31272.8514
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 18012.19502 18057.06539 1.7939546 2000-2009 –> 2010-2019 3.598159 31093.1106 31180.5902
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate_star 18283.22769 18321.68801 1.8493116 2000-2009 –> 2010-2019 3.707876 31555.2128 31630.1952
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 46629.34864 46678.90228 3.4511966 2000-2009 –> 2010-2019 6.303277 83342.2501 83440.2043
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 46743.62408 46793.17773 3.4737551 2000-2009 –> 2010-2019 6.353887 83602.1229 83700.0770
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 46728.21983 46777.77348 3.4707056 2000-2009 –> 2010-2019 6.348983 83577.1370 83675.0911
obs Indo-Pacific (28.05,28.1] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 35851.93567 35893.19075 2.9591420 2000-2009 –> 2010-2019 5.687445 58714.0705 58794.0779
obs Indo-Pacific (28.05,28.1] cstar_tref ~ aou + phosphate + phosphate_star 35926.04135 35960.42058 2.9749173 2000-2009 –> 2010-2019 5.704879 58791.9102 58858.5830
obs Indo-Pacific (28.05,28.1] cstar_tref ~ aou + silicate + phosphate + phosphate_star 35902.67495 35943.93003 2.9696500 2000-2009 –> 2010-2019 5.698634 58767.1641 58847.1715
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 35868.39223 35916.52315 2.9621321 2000-2009 –> 2010-2019 5.702685 58774.7816 58868.1236
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + phosphate 35866.97287 35915.10380 2.9618384 2000-2009 –> 2010-2019 5.716480 58821.7254 58915.0673
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + silicate 9290.46329 9331.46565 1.4554984 2000-2009 –> 2010-2019 3.127041 16737.5713 16817.5052
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 9477.71973 9518.72209 1.5091826 2000-2009 –> 2010-2019 3.087624 16704.4169 16784.3508
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 35866.15828 35914.28921 2.9616699 2000-2009 –> 2010-2019 5.674422 58676.3812 58769.7232
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 9482.97327 9523.97563 1.5107170 2000-2009 –> 2010-2019 3.093579 16720.4273 16800.3611
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 35986.23729 36034.36821 2.9866196 2000-2009 –> 2010-2019 5.722184 58875.4447 58968.7866
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 9432.95003 9473.95239 1.4961702 2000-2009 –> 2010-2019 3.079011 16670.3529 16750.2868
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 35955.51005 36003.64097 2.9802152 2000-2009 –> 2010-2019 5.709921 58824.4956 58917.8376
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 9482.33899 9523.34135 1.5105316 2000-2009 –> 2010-2019 3.093393 16719.7926 16799.7265
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 35987.71750 36035.84842 2.9869285 2000-2009 –> 2010-2019 5.722130 58875.6736 58969.0155
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate_star 9487.86994 9523.01483 1.5127336 2000-2009 –> 2010-2019 3.095597 16723.3275 16791.8422
obs Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 36048.36073 36096.49165 2.9996098 2000-2009 –> 2010-2019 5.722142 58892.5262 58985.8682
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate 9886.05501 9921.19990 1.6338459 2000-2009 –> 2010-2019 3.237057 17170.6377 17239.1524
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate_star 9894.30909 9929.45397 1.6364565 2000-2009 –> 2010-2019 3.245590 17193.0727 17261.5875
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + nitrate + silicate + phosphate_star 9884.24222 9919.38710 1.6332732 2000-2009 –> 2010-2019 3.241943 17181.8970 17250.4117
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + silicate + phosphate + phosphate_star 9892.23756 9927.38244 1.6358009 2000-2009 –> 2010-2019 3.244911 17190.9453 17259.4600
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 18973.45451 19019.82176 1.3299335 2000-2009 –> 2010-2019 2.771672 36449.8202 36541.6549
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 37477.94067 37526.60347 2.7369771 2000-2009 –> 2010-2019 5.271122 67006.8968 67102.7788
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + phosphate + phosphate_star 37542.77630 37584.48727 2.7488473 2000-2009 –> 2010-2019 5.288122 67095.1447 67177.3293
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + silicate + phosphate 37512.44010 37554.15107 2.7434531 2000-2009 –> 2010-2019 5.279457 67048.6099 67130.7945
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + nitrate + silicate + phosphate 19063.18350 19102.92686 1.3409435 2000-2009 –> 2010-2019 2.803853 36680.1750 36758.8905
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 19065.00677 19111.37402 1.3409222 2000-2009 –> 2010-2019 2.797740 36643.1728 36735.0075
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 37592.68844 37641.35124 2.7573883 2000-2009 –> 2010-2019 5.326212 67292.4402 67388.3222
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + nitrate + phosphate 37636.09454 37684.75734 2.7651490 2000-2009 –> 2010-2019 5.339617 67363.4262 67459.3082
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + nitrate + silicate 18063.36471 18109.73196 1.2254774 2000-2009 –> 2010-2019 2.747899 36072.4972 36164.3319
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate 17651.93210 17698.29935 1.1809876 2000-2009 –> 2010-2019 2.556123 34665.5430 34757.3778
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate 37497.14664 37545.80944 2.7403829 2000-2009 –> 2010-2019 5.257675 66942.2545 67038.1366
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 17670.98913 17717.35638 1.1830122 2000-2009 –> 2010-2019 2.560043 34698.0786 34789.9133
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 37599.00250 37647.66530 2.7585159 2000-2009 –> 2010-2019 5.295804 67143.5339 67239.4159
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 18515.61631 18561.98356 1.2763171 2000-2009 –> 2010-2019 2.715295 35973.2360 36065.0707
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 17629.66490 17676.03215 1.1786264 2000-2009 –> 2010-2019 2.565976 34729.7964 34821.6311
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 37675.54308 37724.20588 2.7722210 2000-2009 –> 2010-2019 5.326787 67305.3543 67401.2363
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 17667.43813 17713.80538 1.1826347 2000-2009 –> 2010-2019 2.559398 34692.6241 34784.4589
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 37582.72647 37631.38927 2.7556103 2000-2009 –> 2010-2019 5.288869 67107.2879 67203.1699
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate_star 17844.85951 17884.60287 1.2018608 2000-2009 –> 2010-2019 2.592329 34964.9583 35043.6738
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate_star 37705.77155 37747.48252 2.7780121 2000-2009 –> 2010-2019 5.340571 67372.8383 67455.0229
obs Indo-Pacific (28.2,28.3] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 11304.02585 11339.33105 2.0291561 2000-2009 –> 2010-2019 4.639219 19554.5016 19622.5490
obs Indo-Pacific (28.2,28.3] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 11223.49896 11264.68836 1.9978632 2000-2009 –> 2010-2019 4.605944 19473.3428 19552.7314
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 11246.01949 11287.20889 2.0063545 2000-2009 –> 2010-2019 4.613917 19495.1753 19574.5639
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 31795.42323 31848.74972 0.6963147 2000-2009 –> 2010-2019 1.350590 54011.3866 54115.9537
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 34922.48262 34975.80911 0.7726308 2000-2009 –> 2010-2019 1.541942 60753.5594 60858.1266
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + phosphate + phosphate_star 11377.82615 11413.13135 2.0575550 2000-2009 –> 2010-2019 4.666753 19627.1542 19695.2016
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 11367.94619 11409.13559 2.0529568 2000-2009 –> 2010-2019 4.662020 19619.0948 19698.4834
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + nitrate + silicate + phosphate 32249.09237 32294.80079 0.7069474 2000-2009 –> 2010-2019 1.453828 57417.7102 57507.3392
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 32071.04969 32124.37618 0.7027270 2000-2009 –> 2010-2019 1.362640 54478.4991 54583.0663
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + nitrate + phosphate 11228.46585 11269.65525 1.9997329 2000-2009 –> 2010-2019 4.657906 19544.2105 19623.5991
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 11360.91650 11402.10589 2.0502408 2000-2009 –> 2010-2019 4.658474 19610.9628 19690.3514
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + silicate + phosphate 34854.25187 34907.57836 0.7708795 2000-2009 –> 2010-2019 1.601782 62404.3366 62508.9038
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 34866.43673 34919.76322 0.7711919 2000-2009 –> 2010-2019 1.602397 62424.6473 62529.2144
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 32060.83678 32114.16327 0.7024883 2000-2009 –> 2010-2019 1.441526 56995.8503 57100.4175
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 34864.06369 34917.39018 0.7711311 2000-2009 –> 2010-2019 1.602298 62421.2517 62525.8188
obs Indo-Pacific (28.2,28.3] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 11260.23475 11301.42415 2.0117328 2000-2009 –> 2010-2019 4.621720 19512.6105 19591.9991
mod Indo-Pacific (28.2,28.3] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 36018.81701 36072.14350 0.8013222 2000-2009 –> 2010-2019 1.638815 63745.2449 63849.8121
obs Indo-Pacific (28.2,28.3] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 11231.53801 11272.72741 2.0008902 2000-2009 –> 2010-2019 4.648854 19533.9537 19613.3423
obs Indo-Pacific (28.2,28.3] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 11374.20798 11415.39738 2.0553792 2000-2009 –> 2010-2019 4.664929 19626.0024 19705.3910
mod Indo-Pacific (28.2,28.3] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 36018.28907 36071.61556 0.8013081 2000-2009 –> 2010-2019 1.640051 63778.0060 63882.5731
mod Indo-Pacific (28.3,28.4] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 20732.43039 20779.63400 1.2630285 2000-2009 –> 2010-2019 2.189973 33637.2800 33729.8167
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 20403.82340 20451.02701 1.2303560 2000-2009 –> 2010-2019 2.159120 33327.4944 33420.0312
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 4498.10060 4532.78975 2.0511139 2000-2009 –> 2010-2019 4.699357 7775.9274 7842.2918
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 4481.83788 4516.52703 2.0352761 2000-2009 –> 2010-2019 4.669781 7752.5702 7818.9345
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + silicate + phosphate_star 4502.07708 4531.81064 2.0569652 2000-2009 –> 2010-2019 4.707523 7779.0958 7835.9795
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 4529.61649 4564.30564 2.0821581 2000-2009 –> 2010-2019 4.715596 7799.7963 7866.1606
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + silicate + phosphate + phosphate_star 4530.87905 4560.61261 2.0853985 2000-2009 –> 2010-2019 4.720308 7799.8209 7856.7047
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + nitrate + phosphate 20381.09610 20428.29971 1.2281278 2000-2009 –> 2010-2019 2.154179 33276.6959 33369.2326
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 20463.53024 20510.73385 1.2362290 2000-2009 –> 2010-2019 2.167244 33410.4426 33502.9793
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + silicate + phosphate 4506.77432 4541.46346 2.0596113 2000-2009 –> 2010-2019 4.693139 7777.0008 7843.3652
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4480.82584 4515.51498 2.0342945 2000-2009 –> 2010-2019 4.684007 7759.4091 7825.7735
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 20379.00436 20426.20797 1.2279229 2000-2009 –> 2010-2019 2.155044 33285.6795 33378.2162
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 20353.86662 20401.07022 1.2254635 2000-2009 –> 2010-2019 2.191712 33657.3001 33749.8369
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 4506.14563 4540.83478 2.0589942 2000-2009 –> 2010-2019 4.708995 7784.8776 7851.2419
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4500.54548 4535.23463 2.0535055 2000-2009 –> 2010-2019 4.687868 7771.2044 7837.5688
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + aou + nitrate + phosphate 20712.54600 20753.00624 1.2612282 2000-2009 –> 2010-2019 2.187820 33611.7507 33691.0680
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 20714.18992 20761.39352 1.2611924 2000-2009 –> 2010-2019 2.185096 33587.4993 33680.0361
obs Indo-Pacific (28.3,28.4] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 4478.43487 4513.12402 2.0319775 2000-2009 –> 2010-2019 4.717373 7775.2644 7841.6288
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 20715.30729 20755.76753 1.2615060 2000-2009 –> 2010-2019 2.189225 33626.1730 33705.4902
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 20716.90849 20764.11210 1.2614659 2000-2009 –> 2010-2019 2.186555 33602.5287 33695.0655
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
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
d3c3fde jens-daniel-mueller 2021-03-26
594d503 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
ac1a836 jens-daniel-mueller 2021-02-24
b03fbd3 jens-daniel-mueller 2021-02-24
3d3b4cc jens-daniel-mueller 2021-02-23
7b672f7 jens-daniel-mueller 2021-01-11
33ba23c jens-daniel-mueller 2021-01-07
318609d jens-daniel-mueller 2020-12-23
9d0b2d0 jens-daniel-mueller 2020-12-23
0aa2b50 jens-daniel-mueller 2020-12-23
2886da0 jens-daniel-mueller 2020-12-19
02f0ee9 jens-daniel-mueller 2020-12-18
7bcb4eb jens-daniel-mueller 2020-12-18
7131186 jens-daniel-mueller 2020-12-17
5b48ef5 jens-daniel-mueller 2020-12-17
f3a708f jens-daniel-mueller 2020-12-17

[[2]]

Version Author Date
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
d3c3fde jens-daniel-mueller 2021-03-26
594d503 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26

8.1.2 Best models

left_join(lm_best_target %>% select(basin, gamma_slab, era, model),
          GLODAP_glanced) %>%
  group_split(era) %>%
  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
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
d3c3fde jens-daniel-mueller 2021-03-26
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
ac1a836 jens-daniel-mueller 2021-02-24
b03fbd3 jens-daniel-mueller 2021-02-24
3d3b4cc jens-daniel-mueller 2021-02-23
7b672f7 jens-daniel-mueller 2021-01-11
33ba23c jens-daniel-mueller 2021-01-07
318609d jens-daniel-mueller 2020-12-23
9d0b2d0 jens-daniel-mueller 2020-12-23
0aa2b50 jens-daniel-mueller 2020-12-23
2886da0 jens-daniel-mueller 2020-12-19
02f0ee9 jens-daniel-mueller 2020-12-18
7bcb4eb jens-daniel-mueller 2020-12-18
7131186 jens-daniel-mueller 2020-12-17
5b48ef5 jens-daniel-mueller 2020-12-17
f3a708f jens-daniel-mueller 2020-12-17

[[2]]

Version Author Date
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
d3c3fde jens-daniel-mueller 2021-03-26
594d503 jens-daniel-mueller 2021-03-26

8.2 RMSE correlation between eras

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

8.2.1 All models

# find max rmse to scale axis
max_rmse <-
  max(c(lm_all_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
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
ac1a836 jens-daniel-mueller 2021-02-24
b03fbd3 jens-daniel-mueller 2021-02-24
3d3b4cc jens-daniel-mueller 2021-02-23
7b672f7 jens-daniel-mueller 2021-01-11
33ba23c jens-daniel-mueller 2021-01-07
318609d jens-daniel-mueller 2020-12-23
9d0b2d0 jens-daniel-mueller 2020-12-23
0aa2b50 jens-daniel-mueller 2020-12-23
2886da0 jens-daniel-mueller 2020-12-19
02f0ee9 jens-daniel-mueller 2020-12-18
7bcb4eb jens-daniel-mueller 2020-12-18
7131186 jens-daniel-mueller 2020-12-17
5b48ef5 jens-daniel-mueller 2020-12-17
e4ca289 jens-daniel-mueller 2020-12-16
158fe26 jens-daniel-mueller 2020-12-15
7a9a4cb jens-daniel-mueller 2020-12-15
61b263c jens-daniel-mueller 2020-12-15
984697e jens-daniel-mueller 2020-12-12
3ebff89 jens-daniel-mueller 2020-12-12
ba112d3 jens-daniel-mueller 2020-12-11
24a632f jens-daniel-mueller 2020-12-07
6a8004b jens-daniel-mueller 2020-12-07
70bf1a5 jens-daniel-mueller 2020-12-07
7555355 jens-daniel-mueller 2020-12-07
143d6fa jens-daniel-mueller 2020-12-07
090e4d5 jens-daniel-mueller 2020-12-02
7c25f7a jens-daniel-mueller 2020-12-02
b03ddb8 jens-daniel-mueller 2020-12-02
91435ae jens-daniel-mueller 2020-12-01
196be51 jens-daniel-mueller 2020-11-30
rm(max_rmse)

8.2.2 Best models

# find max rmse to scale axis
max_rmse <-
  max(c(lm_best_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
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
ac1a836 jens-daniel-mueller 2021-02-24
b03fbd3 jens-daniel-mueller 2021-02-24
3d3b4cc jens-daniel-mueller 2021-02-23
7b672f7 jens-daniel-mueller 2021-01-11
33ba23c jens-daniel-mueller 2021-01-07
318609d jens-daniel-mueller 2020-12-23
9d0b2d0 jens-daniel-mueller 2020-12-23
0aa2b50 jens-daniel-mueller 2020-12-23
2886da0 jens-daniel-mueller 2020-12-19
02f0ee9 jens-daniel-mueller 2020-12-18
7bcb4eb jens-daniel-mueller 2020-12-18
7131186 jens-daniel-mueller 2020-12-17
a84ff3c jens-daniel-mueller 2020-12-17
rm(max_rmse)

8.3 Predictor counts

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

# calculate cases of predictor used
lm_all_stats <- lm_best_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
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
ac1a836 jens-daniel-mueller 2021-02-24
b03fbd3 jens-daniel-mueller 2021-02-24
3d3b4cc jens-daniel-mueller 2021-02-23
7b672f7 jens-daniel-mueller 2021-01-11
33ba23c jens-daniel-mueller 2021-01-07
318609d jens-daniel-mueller 2020-12-23
9d0b2d0 jens-daniel-mueller 2020-12-23
0aa2b50 jens-daniel-mueller 2020-12-23
2886da0 jens-daniel-mueller 2020-12-19
02f0ee9 jens-daniel-mueller 2020-12-18
7bcb4eb jens-daniel-mueller 2020-12-18
7131186 jens-daniel-mueller 2020-12-17
5b48ef5 jens-daniel-mueller 2020-12-17
lm_best_eras %>% 
  ggplot(aes(rmse, AIC, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  facet_grid(data_source~basin)

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

8.4.2 AIC vs BIC

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

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

Version Author Date
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20

8.5 RMSE vs VIF

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

Version Author Date
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20

8.6 Residual patterns

8.6.1 Fitted vs 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
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
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
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
27ae473 jens-daniel-mueller 2021-02-24
4bc00ea jens-daniel-mueller 2021-02-24
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
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
27ae473 jens-daniel-mueller 2021-02-24
4bc00ea jens-daniel-mueller 2021-02-24

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

GLODAP_augmented_best %>%
  ggplot(aes(lat, .resid)) +
  geom_bin2d(binwidth = i_binwidth) +
  geom_hline(yintercept = 0, col = "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
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
d3c3fde jens-daniel-mueller 2021-03-26
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
27ae473 jens-daniel-mueller 2021-02-24
4bc00ea jens-daniel-mueller 2021-02-24

8.6.3 By model domain

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

GLODAP_augmented_best %>% 
  group_split(basin, 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
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
d3c3fde jens-daniel-mueller 2021-03-26
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
27ae473 jens-daniel-mueller 2021-02-24
a1ba577 jens-daniel-mueller 2021-02-24
071743d jens-daniel-mueller 2021-02-24

[[2]]

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

[[3]]

Version Author Date
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
d3c3fde jens-daniel-mueller 2021-03-26
594d503 jens-daniel-mueller 2021-03-26

[[4]]

Version Author Date
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
d3c3fde jens-daniel-mueller 2021-03-26
594d503 jens-daniel-mueller 2021-03-26

8.6.4 Latitudinal mean

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

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

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

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

[[2]]

Version Author Date
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26

8.6.5 Latitudinal offset

# calculate residual offset for adjacent eras
lat_residual_offset <- lat_residual  %>%
  select(data_source, basin, gamma_slab, era, lat_grid, .resid_mean) %>% 
  arrange(era) %>% 
  group_by(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
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
fec3558 jens-daniel-mueller 2021-02-24

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
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
d37a85d jens-daniel-mueller 2021-05-31
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
330dcd0 jens-daniel-mueller 2021-03-20
GLODAP_glanced_vif %>% 
  ggplot(aes(Variables,
             log10(VIF),
             col = data_source)) + 
  geom_hline(yintercept = log10(params_local$vif_max)) +
  geom_violin() +
  scale_color_brewer(palette = "Set1") +
  facet_grid(gamma_slab~basin) +
  coord_flip()

Version Author Date
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
d37a85d jens-daniel-mueller 2021-05-31
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
d3c3fde jens-daniel-mueller 2021-03-26
594d503 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
330dcd0 jens-daniel-mueller 2021-03-20

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

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

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

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

other attached packages:
 [1] tictoc_1.0.1     gt_0.2.2         corrr_0.4.3      broom_0.7.5     
 [5] kableExtra_1.3.1 knitr_1.30       olsrr_0.5.3      GGally_2.0.0    
 [9] lubridate_1.7.9  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