Last updated: 2021-06-18

Checks: 7 0

Knit directory: emlr_obs_v_XXX/

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

R Markdown file: up-to-date

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

Environment: empty

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

Seed: set.seed(20200707)

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

Session information: recorded

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

Cache: none

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

File paths: relative

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

Repository version: 9a6e1b8

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

1 Version ID

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

2 Required data

Required are:

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

3 Predictor combinations

Find all possible combinations of following considered predictor variables:

  • sal, temp, aou, 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 133.20428 146.72706 0.7785912 2000-2009 –> 2010-2019 4.828984 789.6518 822.3280
obs Atlantic (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 850.01040 869.86437 6.6762940 2000-2009 –> 2010-2019 14.634025 1276.1941 1310.5909
mod Atlantic (-Inf,26] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 80.56216 94.08494 0.4646976 2000-2009 –> 2010-2019 4.848935 755.0677 787.7439
mod Atlantic (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 110.47653 123.99931 0.6230739 2000-2009 –> 2010-2019 4.608223 763.2216 795.8977
obs Atlantic (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 842.75911 862.61309 6.4869217 2000-2009 –> 2010-2019 15.398851 1282.3060 1316.7027
mod Atlantic (-Inf,26] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 77.94453 91.46731 0.4529238 2000-2009 –> 2010-2019 4.325639 724.1644 756.8406
mod Atlantic (-Inf,26] cstar_tref ~ sal + nitrate + silicate + phosphate_star 85.22918 96.82014 0.4960863 2000-2009 –> 2010-2019 4.893068 758.3966 786.4047
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate 852.23682 869.25451 6.7892097 2000-2009 –> 2010-2019 14.842531 1277.8295 1307.3125
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 850.55882 870.41279 6.6908393 2000-2009 –> 2010-2019 14.680210 1277.2108 1311.6075
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 850.60290 870.45687 6.6920097 2000-2009 –> 2010-2019 14.678202 1277.2079 1311.6046
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 79.09507 92.61785 0.4580616 2000-2009 –> 2010-2019 4.750542 748.7782 781.4544
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 845.36249 865.21646 6.5542846 2000-2009 –> 2010-2019 14.540323 1271.9652 1306.3620
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 109.81195 123.33473 0.6190274 2000-2009 –> 2010-2019 4.496233 756.2961 788.9722
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 133.10955 146.63233 0.7778684 2000-2009 –> 2010-2019 4.814008 788.7534 821.4295
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 850.79437 870.64834 6.6970962 2000-2009 –> 2010-2019 14.706869 1277.7473 1312.1440
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 82.16312 95.68590 0.4720489 2000-2009 –> 2010-2019 4.780173 752.6757 785.3519
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 109.72228 123.24506 0.6184835 2000-2009 –> 2010-2019 4.805876 773.7541 806.4303
obs Atlantic (-Inf,26] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 856.56356 876.41753 6.8521855 2000-2009 –> 2010-2019 15.343541 1290.4060 1324.8028
obs Atlantic (-Inf,26] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 855.67327 875.52725 6.8280203 2000-2009 –> 2010-2019 15.311211 1289.4022 1323.7990
obs Atlantic (-Inf,26] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 853.00049 872.85446 6.7559830 2000-2009 –> 2010-2019 15.203064 1286.2261 1320.6228
obs Atlantic (26,26.5] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 3974.10444 4000.65365 6.0005013 2000-2009 –> 2010-2019 11.469864 7228.9290 7281.0012
mod Atlantic (26,26.5] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 2087.77673 2115.26521 3.8423998 2000-2009 –> 2010-2019 8.312775 3449.1586 3500.7440
obs Atlantic (26,26.5] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 3975.56134 4006.53542 5.9978610 2000-2009 –> 2010-2019 11.444493 7228.0547 7288.8056
mod Atlantic (26,26.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 2051.74656 2079.23504 3.6621737 2000-2009 –> 2010-2019 7.668550 3362.5000 3414.0854
obs Atlantic (26,26.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 3854.21802 3885.19210 5.4361418 2000-2009 –> 2010-2019 10.893903 7108.8342 7169.5851
mod Atlantic (26,26.5] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2016.34783 2043.83631 3.4933409 2000-2009 –> 2010-2019 7.884513 3369.4709 3421.0563
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 2055.57634 2083.06483 3.6809220 2000-2009 –> 2010-2019 7.968550 3397.6749 3449.2603
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 3871.18290 3902.15699 5.5113933 2000-2009 –> 2010-2019 11.038074 7138.8498 7199.6007
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 2055.84998 2083.33846 3.6822652 2000-2009 –> 2010-2019 8.039672 3405.4068 3456.9922
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 3872.32519 3903.29928 5.5164975 2000-2009 –> 2010-2019 11.117023 7153.7961 7214.5470
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 2020.75861 2048.24709 3.5139460 2000-2009 –> 2010-2019 8.252370 3409.0438 3460.6292
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2021.50725 2048.99573 3.5174553 2000-2009 –> 2010-2019 8.373161 3421.0882 3472.6736
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 2056.90285 2084.39133 3.6874381 2000-2009 –> 2010-2019 8.117415 3414.0906 3465.6760
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 3882.63304 3913.60712 5.5627709 2000-2009 –> 2010-2019 11.224283 7175.3677 7236.1186
obs Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + phosphate 3972.94136 3999.49057 5.9948483 2000-2009 –> 2010-2019 11.556013 7245.0773 7297.1495
mod Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 2085.18485 2112.67333 3.8291440 2000-2009 –> 2010-2019 8.339638 3450.6943 3502.2797
obs Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 3974.47183 4005.44591 5.9925678 2000-2009 –> 2010-2019 11.475141 7233.8055 7294.5563
mod Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2081.15868 2108.64716 3.8086433 2000-2009 –> 2010-2019 8.355552 3450.3831 3501.9685
obs Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 3974.65317 4005.62726 5.9934485 2000-2009 –> 2010-2019 11.518857 7242.0807 7302.8315
obs Atlantic (26,26.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 3975.54638 4006.52046 5.9977883 2000-2009 –> 2010-2019 11.578214 7253.2781 7314.0290
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 2726.24594 2755.50606 4.0099127 2000-2009 –> 2010-2019 7.178087 4628.1629 4684.7606
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 4781.54070 4814.13778 5.1811842 2000-2009 –> 2010-2019 10.852531 8050.8455 8113.1653
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 2729.31122 2758.57133 4.0226571 2000-2009 –> 2010-2019 7.081366 4605.4188 4662.0164
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 4780.48251 4813.07960 5.1776618 2000-2009 –> 2010-2019 10.878187 8055.0833 8117.4032
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 2794.11765 2823.37777 4.3017850 2000-2009 –> 2010-2019 7.154411 4619.0266 4675.6242
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 4869.31685 4901.91394 5.4818632 2000-2009 –> 2010-2019 11.179732 8143.4367 8205.7565
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2800.32942 2829.58954 4.3295363 2000-2009 –> 2010-2019 7.128169 4611.2125 4667.8102
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 2736.48546 2765.74558 4.0526436 2000-2009 –> 2010-2019 6.977944 4579.8598 4636.4574
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 4777.62281 4810.21990 5.1681548 2000-2009 –> 2010-2019 10.867326 8051.9784 8114.2982
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 2801.33775 2830.59787 4.3340579 2000-2009 –> 2010-2019 7.128753 4611.1873 4667.7849
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + phosphate 4838.11488 4866.05524 5.3799427 2000-2009 –> 2010-2019 11.055444 8106.1755 8159.5925
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + phosphate_star 4835.01908 4862.95944 5.3692494 2000-2009 –> 2010-2019 11.070156 8107.6889 8161.1059
mod Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 2752.15393 2781.41405 4.1189132 2000-2009 –> 2010-2019 7.187050 4630.5204 4687.1181
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 4837.91003 4870.50711 5.3723247 2000-2009 –> 2010-2019 11.047104 8107.8393 8170.1591
mod Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2755.86354 2785.12366 4.1347610 2000-2009 –> 2010-2019 7.102801 4609.8843 4666.4820
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 4834.98173 4867.57882 5.3622238 2000-2009 –> 2010-2019 11.055313 8108.2355 8170.5553
mod Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 2763.03028 2788.11038 4.1741840 2000-2009 –> 2010-2019 7.239098 4638.6254 4687.1376
obs Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 4832.62335 4860.56371 5.3609889 2000-2009 –> 2010-2019 11.061212 8105.1695 8158.5865
mod Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 2760.14618 2789.40630 4.1531326 2000-2009 –> 2010-2019 7.073069 4602.1736 4658.7712
obs Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 4832.55491 4865.15200 5.3538671 2000-2009 –> 2010-2019 11.046869 8105.7928 8168.1126
mod Atlantic (26.75,27] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 4240.04990 4273.77471 2.4422139 2000-2009 –> 2010-2019 4.769929 7088.3598 7153.1601
mod Atlantic (26.75,27] cstar_tref ~ aou + silicate + phosphate + phosphate_star 4641.91544 4670.82242 3.0460265 2000-2009 –> 2010-2019 5.408694 7506.8852 7562.4283
mod Atlantic (26.75,27] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 4390.25346 4423.97827 2.6513613 2000-2009 –> 2010-2019 5.180811 7342.6225 7407.4228
mod Atlantic (26.75,27] cstar_tref ~ sal + aou + silicate + phosphate 4558.14441 4587.05140 2.9095876 2000-2009 –> 2010-2019 5.445477 7511.6973 7567.2403
mod Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate 4541.20077 4574.92559 2.8795912 2000-2009 –> 2010-2019 5.154980 7361.0449 7425.8452
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate 5657.73200 5691.90906 4.3722730 2000-2009 –> 2010-2019 10.750722 9683.0367 9748.1423
mod Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4560.71430 4594.43911 2.9104949 2000-2009 –> 2010-2019 5.259952 7420.6639 7485.4642
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5661.38471 5695.56177 4.3804707 2000-2009 –> 2010-2019 10.702789 9675.8528 9740.9584
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 5669.83656 5704.01363 4.3994981 2000-2009 –> 2010-2019 10.708107 9681.6433 9746.7489
mod Atlantic (26.75,27] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4567.83682 4601.56163 2.9218573 2000-2009 –> 2010-2019 5.304699 7445.4519 7510.2522
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5662.30907 5696.48614 4.3825477 2000-2009 –> 2010-2019 10.701740 9676.1707 9741.2763
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + silicate + phosphate_star 5679.09957 5708.39420 4.4249827 2000-2009 –> 2010-2019 10.830755 9707.6448 9763.4496
mod Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 4256.43285 4290.15766 2.4641999 2000-2009 –> 2010-2019 4.739906 7076.4516 7141.2519
obs Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 5661.85333 5696.03039 4.3815236 2000-2009 –> 2010-2019 10.850500 9704.4359 9769.5415
mod Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 4265.72883 4299.45364 2.4767631 2000-2009 –> 2010-2019 4.769408 7095.0320 7159.8322
obs Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5666.46731 5700.64438 4.3919032 2000-2009 –> 2010-2019 10.773403 9692.3582 9757.4638
mod Atlantic (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 4288.04186 4321.76668 2.5071804 2000-2009 –> 2010-2019 4.833031 7135.3483 7200.1486
obs Atlantic (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 5668.24226 5702.41933 4.3959026 2000-2009 –> 2010-2019 10.788955 9696.3507 9761.4563
obs Atlantic (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate_star 5671.63292 5700.92755 4.4080716 2000-2009 –> 2010-2019 10.831986 9703.6454 9759.4502
obs Atlantic (26.75,27] cstar_tref ~ temp + silicate + phosphate + phosphate_star 5667.69915 5696.99377 4.3991881 2000-2009 –> 2010-2019 10.850607 9704.9496 9760.7545
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 5403.86805 5439.85122 2.0472724 2000-2009 –> 2010-2019 4.365709 9390.6518 9460.0865
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 5547.73125 5583.71442 2.1673528 2000-2009 –> 2010-2019 4.646269 9652.1750 9721.6097
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate 5694.28779 5725.13051 2.2987467 2000-2009 –> 2010-2019 4.762198 9785.7300 9845.2454
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 5696.26202 5732.24519 2.2987233 2000-2009 –> 2010-2019 4.643595 9702.9779 9772.4125
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate 11058.27890 11096.38864 6.1127868 2000-2009 –> 2010-2019 12.010184 19030.3556 19104.3592
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 11031.64720 11069.75694 6.0653712 2000-2009 –> 2010-2019 11.848742 18955.0690 19029.0726
mod Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + silicate 5533.41566 5569.39883 2.1550948 2000-2009 –> 2010-2019 4.420060 9479.1790 9548.6137
mod Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate 5692.59621 5728.57938 2.2953871 2000-2009 –> 2010-2019 4.674960 9725.1376 9794.5723
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 11082.84697 11120.95671 6.1568570 2000-2009 –> 2010-2019 11.895943 18987.1133 19061.1169
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 11038.83863 11076.94837 6.0781386 2000-2009 –> 2010-2019 11.894586 18976.4723 19050.4759
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 11102.33117 11140.44091 6.1920336 2000-2009 –> 2010-2019 12.026704 19047.7600 19121.7636
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 11084.56702 11122.67676 6.1599544 2000-2009 –> 2010-2019 11.912329 18994.5970 19068.6006
obs Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + phosphate_star 11055.63158 11088.29707 6.1116300 2000-2009 –> 2010-2019 12.022760 19031.5041 19094.9358
mod Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate 5702.41173 5733.25445 2.3061576 2000-2009 –> 2010-2019 4.609124 9675.4257 9734.9411
mod Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 5704.30230 5740.28547 2.3060576 2000-2009 –> 2010-2019 4.567044 9646.9746 9716.4092
obs Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 11021.73202 11059.84176 6.0478120 2000-2009 –> 2010-2019 12.022733 19026.3534 19100.3570
mod Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5703.44005 5739.42322 2.3052699 2000-2009 –> 2010-2019 4.565080 9645.1975 9714.6322
obs Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 10979.71384 11017.82359 5.9739631 2000-2009 –> 2010-2019 11.816639 18928.5593 19002.5629
mod Atlantic (27,27.25] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 5607.75503 5643.73820 2.2195128 2000-2009 –> 2010-2019 4.481452 9551.1681 9620.6028
obs Atlantic (27,27.25] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 10991.31808 11029.42783 5.9942675 2000-2009 –> 2010-2019 11.871151 18954.7113 19028.7149
obs Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 13162.57559 13202.02272 5.7952923 2000-2009 –> 2010-2019 10.476098 21856.4333 21932.9077
mod Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate 11047.45464 11081.20528 3.5749435 2000-2009 –> 2010-2019 6.903081 20716.5130 20783.3753
mod Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 10902.69104 10942.06679 3.4491779 2000-2009 –> 2010-2019 6.739368 20531.5029 20609.5089
mod Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 10912.19521 10951.57096 3.4571866 2000-2009 –> 2010-2019 6.769214 20565.3780 20643.3840
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 13147.09382 13186.54095 5.7736609 2000-2009 –> 2010-2019 10.375870 21791.3353 21867.8097
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 11160.29079 11199.66654 3.6729523 2000-2009 –> 2010-2019 6.938299 20761.1810 20839.1870
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 13145.03800 13184.48513 5.7707946 2000-2009 –> 2010-2019 10.347995 21773.3146 21849.7890
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + silicate 11021.69698 11061.07273 3.5508106 2000-2009 –> 2010-2019 6.799013 20603.1933 20681.1993
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + phosphate 13214.29659 13248.10842 5.8709825 2000-2009 –> 2010-2019 10.487842 21865.8508 21931.4003
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + phosphate_star 13217.31708 13251.12890 5.8752675 2000-2009 –> 2010-2019 10.469489 21854.4689 21920.0184
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 13211.68655 13251.13367 5.8644486 2000-2009 –> 2010-2019 10.479357 21864.0021 21940.4765
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 10887.49001 10926.86576 3.4364073 2000-2009 –> 2010-2019 7.101098 20913.4320 20991.4380
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 13214.51595 13253.96307 5.8684579 2000-2009 –> 2010-2019 10.462571 21853.5986 21930.0730
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 11156.48095 11195.85670 3.6695392 2000-2009 –> 2010-2019 6.802146 20604.4828 20682.4888
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 13142.25447 13181.70160 5.7669159 2000-2009 –> 2010-2019 10.349249 21773.8148 21850.2892
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 10866.83154 10906.20729 3.4191275 2000-2009 –> 2010-2019 6.956723 20762.7412 20840.7472
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + phosphate + phosphate_star 13214.18596 13247.99778 5.8708256 2000-2009 –> 2010-2019 10.468863 21853.7703 21919.3198
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 10888.85808 10928.23383 3.4375547 2000-2009 –> 2010-2019 7.104676 20917.2427 20995.2487
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 13211.56430 13251.01143 5.8642754 2000-2009 –> 2010-2019 10.461835 21852.8446 21929.3190
mod Atlantic (27.25,27.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 11184.59422 11223.96997 3.6947997 2000-2009 –> 2010-2019 6.820733 20624.7397 20702.7457
obs Atlantic (27.5,27.75] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 14125.48140 14165.72666 5.0648727 2000-2009 –> 2010-2019 10.335657 25249.9915 25328.7172
obs Atlantic (27.5,27.75] cstar_tref ~ aou + silicate + phosphate + phosphate_star 14136.32193 14170.81786 5.0789084 2000-2009 –> 2010-2019 10.350836 25259.6145 25327.0936
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + nitrate + phosphate 18255.07010 18293.05528 2.1793392 2000-2009 –> 2010-2019 4.614714 33880.7927 33955.5335
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 17796.25665 17840.57270 2.0616405 2000-2009 –> 2010-2019 4.490279 33405.2514 33492.4489
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate 18532.37847 18570.36365 2.2533824 2000-2009 –> 2010-2019 4.695331 34176.3299 34251.0707
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 18500.65257 18544.96861 2.2442446 2000-2009 –> 2010-2019 4.679002 34126.6605 34213.8580
obs Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 14138.16083 14178.40609 5.0787321 2000-2009 –> 2010-2019 10.274868 25211.2358 25289.9615
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 18145.32693 18189.64298 2.1501953 2000-2009 –> 2010-2019 4.569895 33729.3860 33816.5835
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 18530.89553 18575.21158 2.2524370 2000-2009 –> 2010-2019 4.674727 34122.1899 34209.3875
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 14143.39186 14183.63712 5.0844610 2000-2009 –> 2010-2019 10.287957 25221.5708 25300.2965
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 18346.41559 18390.73164 2.2029255 2000-2009 –> 2010-2019 4.582854 33818.4072 33905.6047
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 14171.55417 14211.79943 5.1154148 2000-2009 –> 2010-2019 10.346530 25268.8226 25347.5483
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 18464.84193 18509.15797 2.2345826 2000-2009 –> 2010-2019 4.624065 33963.9241 34051.1217
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 18454.77088 18499.08693 2.2318729 2000-2009 –> 2010-2019 4.619435 33948.4190 34035.6166
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 18351.26430 18395.58035 2.2042128 2000-2009 –> 2010-2019 4.585906 33828.2701 33915.4676
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 14152.07539 14192.32065 5.0939852 2000-2009 –> 2010-2019 10.310069 25238.9668 25317.6925
obs Atlantic (27.5,27.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 14131.78783 14172.03309 5.0717613 2000-2009 –> 2010-2019 10.359530 25267.8996 25346.6253
obs Atlantic (27.5,27.75] cstar_tref ~ temp + aou + silicate + phosphate 14141.52311 14176.01904 5.0846048 2000-2009 –> 2010-2019 10.372979 25276.0478 25343.5270
obs Atlantic (27.5,27.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 14141.28549 14181.53075 5.0821534 2000-2009 –> 2010-2019 10.393916 25293.7230 25372.4487
obs Atlantic (27.5,27.75] cstar_tref ~ temp + silicate + phosphate + phosphate_star 14150.52416 14185.02010 5.0944779 2000-2009 –> 2010-2019 10.406410 25301.0767 25368.5558
obs Atlantic (27.75,27.85] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 9234.50317 9271.64419 5.3509175 2000-2009 –> 2010-2019 10.528063 16961.6214 17034.7289
obs Atlantic (27.75,27.85] cstar_tref ~ aou + silicate + phosphate + phosphate_star 9232.50656 9264.34172 5.3509236 2000-2009 –> 2010-2019 10.552889 16969.6677 17032.3313
mod Atlantic (27.75,27.85] cstar_tref ~ sal + aou + nitrate + silicate 5173.27778 5205.24840 1.3171969 2000-2009 –> 2010-2019 2.554689 9444.8645 9507.8788
obs Atlantic (27.75,27.85] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 9220.21597 9257.35699 5.3253075 2000-2009 –> 2010-2019 10.518549 16955.1509 17028.2584
mod Atlantic (27.75,27.85] cstar_tref ~ sal + aou + silicate + phosphate 4951.66626 4983.63688 1.2247674 2000-2009 –> 2010-2019 2.411198 9113.2727 9176.2871
obs Atlantic (27.75,27.85] cstar_tref ~ sal + aou + silicate + phosphate 9218.24136 9250.07652 5.3253529 2000-2009 –> 2010-2019 10.546422 16964.6327 17027.2963
obs Atlantic (27.75,27.85] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 9221.44404 9258.58506 5.3275040 2000-2009 –> 2010-2019 10.552126 16971.5483 17044.6558
mod Atlantic (27.75,27.85] cstar_tref ~ sal + nitrate + silicate + phosphate_star 5576.91446 5608.88508 1.5038366 2000-2009 –> 2010-2019 2.869482 10105.6903 10168.7047
mod Atlantic (27.75,27.85] cstar_tref ~ sal + silicate + phosphate + phosphate_star 5568.30367 5600.27429 1.4995914 2000-2009 –> 2010-2019 2.871111 10108.2819 10171.2963
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + silicate 5129.48191 5166.78098 1.2975413 2000-2009 –> 2010-2019 2.531592 9395.7989 9469.3156
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 4831.14217 4868.44123 1.1764794 2000-2009 –> 2010-2019 2.358192 8984.3507 9057.8675
obs Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 9218.96179 9256.10281 5.3230652 2000-2009 –> 2010-2019 10.533285 16962.1153 17035.2228
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4961.93752 4999.23658 1.2280978 2000-2009 –> 2010-2019 2.448935 9200.1583 9273.6751
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 4967.22366 5004.52272 1.2302309 2000-2009 –> 2010-2019 2.449686 9202.4878 9276.0045
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4955.77316 4993.07222 1.2256149 2000-2009 –> 2010-2019 2.442477 9185.4827 9258.9995
obs Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 9220.68093 9257.82195 5.3261390 2000-2009 –> 2010-2019 10.545875 16968.4290 17041.5365
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate_star 5046.64648 5078.61710 1.2635597 2000-2009 –> 2010-2019 2.567836 9455.4171 9518.4315
obs Atlantic (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 9234.51230 9271.65332 5.3509339 2000-2009 –> 2010-2019 10.535016 16965.0021 17038.1096
obs Atlantic (27.75,27.85] cstar_tref ~ temp + aou + silicate + phosphate 9232.51485 9264.35001 5.3509385 2000-2009 –> 2010-2019 10.561173 16973.6757 17036.3393
obs Atlantic (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 9234.84397 9271.98499 5.3515299 2000-2009 –> 2010-2019 10.543455 16969.1403 17042.2478
mod Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate 6272.40472 6304.28405 1.9500657 2000-2009 –> 2010-2019 4.869058 12778.7861 12841.7046
obs Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 19474.46730 19516.85785 5.3018962 2000-2009 –> 2010-2019 10.572677 36144.9858 36228.6912
mod Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 6273.24743 6310.43997 1.9493135 2000-2009 –> 2010-2019 4.790673 12711.3284 12784.7333
obs Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 18866.52719 18908.91773 4.8144766 2000-2009 –> 2010-2019 9.851672 35291.9972 35375.7026
mod Atlantic (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate 6038.61453 6070.49385 1.8038675 2000-2009 –> 2010-2019 4.538472 12374.8200 12437.7384
obs Atlantic (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate_star 18869.94130 18906.27606 4.8186132 2000-2009 –> 2010-2019 9.856083 35293.7060 35365.4534
obs Atlantic (27.85,27.95] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 19271.50147 19313.89201 5.1339132 2000-2009 –> 2010-2019 10.327670 35862.4375 35946.1428
obs Atlantic (27.85,27.95] cstar_tref ~ sal + silicate + phosphate + phosphate_star 19310.84064 19347.17540 5.1676900 2000-2009 –> 2010-2019 10.364971 35903.4438 35975.1913
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + silicate 6111.96731 6149.15985 1.8472854 2000-2009 –> 2010-2019 4.753523 12608.9281 12682.3330
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + phosphate 6515.02096 6546.90028 2.1143240 2000-2009 –> 2010-2019 4.866154 12867.6024 12930.5209
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 5733.32090 5770.51344 1.6282434 2000-2009 –> 2010-2019 4.362610 12071.2993 12144.7041
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 18704.62372 18747.01427 4.6924023 2000-2009 –> 2010-2019 9.698831 35096.9733 35180.6787
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5839.48834 5876.68089 1.6868973 2000-2009 –> 2010-2019 4.446085 12201.0337 12274.4386
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 18684.21845 18726.60900 4.6772381 2000-2009 –> 2010-2019 9.678906 35071.4241 35155.1295
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 18687.93686 18730.32740 4.6799978 2000-2009 –> 2010-2019 9.725341 35122.1441 35205.8494
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 5848.44487 5885.63741 1.6919411 2000-2009 –> 2010-2019 4.495381 12251.4858 12324.8907
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate 18709.80754 18746.14230 4.6977527 2000-2009 –> 2010-2019 9.744470 35143.4872 35215.2347
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5828.96432 5866.15687 1.6809900 2000-2009 –> 2010-2019 4.431644 12182.4314 12255.8362
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 18703.23534 18745.62588 4.6913690 2000-2009 –> 2010-2019 9.696642 35094.3371 35178.0425
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate_star 5933.56445 5965.44378 1.7417952 2000-2009 –> 2010-2019 4.774346 12539.4822 12602.4007
obs Atlantic (27.95,28.05] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 19921.08757 19963.60575 5.3755125 2000-2009 –> 2010-2019 12.082777 37688.5487 37772.3010
obs Atlantic (27.95,28.05] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 19501.09300 19543.61119 5.0351041 2000-2009 –> 2010-2019 11.597223 37151.6398 37235.3921
obs Atlantic (27.95,28.05] cstar_tref ~ sal + aou + silicate + phosphate_star 19585.42251 19621.86667 5.1032683 2000-2009 –> 2010-2019 11.679712 37245.6229 37317.4105
obs Atlantic (27.95,28.05] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 19638.92387 19681.44205 5.1443714 2000-2009 –> 2010-2019 11.751375 37325.8993 37409.6516
obs Atlantic (27.95,28.05] cstar_tref ~ sal + silicate + phosphate + phosphate_star 19765.42378 19801.86793 5.2483763 2000-2009 –> 2010-2019 11.869262 37461.6163 37533.4039
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate 9604.94833 9637.86801 3.5596708 2000-2009 –> 2010-2019 7.360053 17449.4387 17513.9173
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate_star 9634.79964 9667.71932 3.5895773 2000-2009 –> 2010-2019 7.388781 17478.4075 17542.8861
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 9468.46427 9506.87056 3.4241568 2000-2009 –> 2010-2019 7.223678 17314.3104 17389.5355
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 19618.34528 19660.86346 5.1279081 2000-2009 –> 2010-2019 11.727070 37298.9743 37382.7265
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 9480.13041 9518.53671 3.4353710 2000-2009 –> 2010-2019 7.232163 17323.9327 17399.1578
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 19563.05023 19605.56841 5.0839314 2000-2009 –> 2010-2019 11.659834 37224.8103 37308.5626
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate + phosphate_star 9698.54418 9731.46386 3.6542836 2000-2009 –> 2010-2019 7.440392 17532.3328 17596.8114
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 19654.52209 19697.04027 5.1568855 2000-2009 –> 2010-2019 11.872204 37428.3966 37512.1489
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 9453.60771 9492.01400 3.4099289 2000-2009 –> 2010-2019 7.160032 17262.2209 17337.4459
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + phosphate + phosphate_star 9630.07071 9662.99039 3.5848230 2000-2009 –> 2010-2019 7.384025 17473.6774 17538.1560
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + phosphate_star 9710.98529 9738.41836 3.6691040 2000-2009 –> 2010-2019 7.471194 17554.7536 17608.4858
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate 19755.04214 19791.48630 5.2398961 2000-2009 –> 2010-2019 11.959706 37530.4895 37602.2772
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 9480.39723 9518.80352 3.4356279 2000-2009 –> 2010-2019 7.232425 17324.2034 17399.4285
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 19602.25066 19644.76884 5.1150688 2000-2009 –> 2010-2019 11.705552 37275.8475 37359.5997
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate_star 9479.13891 9512.05859 3.4363421 2000-2009 –> 2010-2019 7.233147 17320.9510 17385.4296
obs Atlantic (28.05,28.1] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 14398.11319 14438.40656 5.2585368 2000-2009 –> 2010-2019 11.595937 25991.2470 26069.9034
obs Atlantic (28.05,28.1] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 14341.09742 14381.39079 5.1947531 2000-2009 –> 2010-2019 11.531951 25934.1179 26012.7742
obs Atlantic (28.05,28.1] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 14344.59262 14384.88599 5.1986408 2000-2009 –> 2010-2019 11.523093 25930.4740 26009.1304
obs Atlantic (28.05,28.1] cstar_tref ~ sal + silicate + phosphate + phosphate_star 14430.00752 14464.54469 5.2968251 2000-2009 –> 2010-2019 11.644785 26027.0448 26094.4646
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 14466.89935 14507.19272 5.3365314 2000-2009 –> 2010-2019 11.641729 26041.9687 26120.6251
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate 3550.02695 3579.70304 1.3279987 2000-2009 –> 2010-2019 2.321854 5717.9210 5775.4207
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate_star 3589.12856 3618.80465 1.3532242 2000-2009 –> 2010-2019 2.365460 5784.9871 5842.4868
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 3551.21426 3585.83636 1.3274795 2000-2009 –> 2010-2019 2.302727 5692.2664 5759.3493
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 3588.31822 3622.94032 1.3513953 2000-2009 –> 2010-2019 2.349731 5765.0763 5832.1592
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 14377.74632 14418.03969 5.2356629 2000-2009 –> 2010-2019 11.507084 25933.7692 26012.4255
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 14460.01635 14500.30972 5.3286751 2000-2009 –> 2010-2019 11.637154 26036.9306 26115.5869
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + phosphate_star 3635.12750 3664.80358 1.3835134 2000-2009 –> 2010-2019 2.406751 5847.4822 5904.9818
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + phosphate_star 14468.22003 14502.75720 5.3403257 2000-2009 –> 2010-2019 11.649232 26043.3746 26110.7943
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 3610.38538 3645.00748 1.3658228 2000-2009 –> 2010-2019 2.388237 5823.5124 5890.5953
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 14347.07778 14387.37114 5.2014069 2000-2009 –> 2010-2019 11.492487 25914.1989 25992.8553
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + phosphate + phosphate_star 3577.79772 3607.47380 1.3458655 2000-2009 –> 2010-2019 2.353698 5767.0041 5824.5037
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 3578.39131 3613.01341 1.3449549 2000-2009 –> 2010-2019 2.336808 5745.2081 5812.2910
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 14419.43834 14459.73171 5.2825941 2000-2009 –> 2010-2019 11.580014 25990.1310 26068.7874
mod Atlantic (28.05,28.1] cstar_tref ~ temp + aou + phosphate 3569.54698 3594.27705 1.3418231 2000-2009 –> 2010-2019 2.419650 5859.2169 5907.1333
mod Atlantic (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate 3571.32803 3601.00411 1.3416817 2000-2009 –> 2010-2019 2.395002 5827.9003 5885.3999
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 16563.29999 16604.30507 5.9348917 2000-2009 –> 2010-2019 15.024002 27622.1433 27700.4474
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + phosphate 4575.85016 4607.85221 1.0741628 2000-2009 –> 2010-2019 2.400891 8947.4889 9010.4233
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + phosphate_star 4622.99040 4654.99246 1.0908278 2000-2009 –> 2010-2019 2.425854 9010.6026 9073.5371
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + silicate + phosphate 4541.26026 4578.59599 1.0614033 2000-2009 –> 2010-2019 2.269580 8675.0866 8748.5101
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4598.05039 4635.38612 1.0812726 2000-2009 –> 2010-2019 2.303402 8761.2940 8834.7175
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 16596.68605 16637.69112 5.9733263 2000-2009 –> 2010-2019 15.078687 27660.9705 27739.2747
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 16546.37739 16587.38246 5.9155047 2000-2009 –> 2010-2019 14.998466 27603.1593 27681.4634
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + phosphate_star 4639.57807 4671.58013 1.0967531 2000-2009 –> 2010-2019 2.406290 8977.8030 9040.7375
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 4641.03015 4678.36588 1.0965569 2000-2009 –> 2010-2019 2.375887 8921.4646 8994.8881
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 16570.28199 16611.28706 5.9429090 2000-2009 –> 2010-2019 15.052017 27635.8199 27714.1240
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + phosphate + phosphate_star 4607.93145 4639.93351 1.0854762 2000-2009 –> 2010-2019 2.414843 8984.6612 9047.5956
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + phosphate + phosphate_star 16598.55183 16633.69903 5.9777927 2000-2009 –> 2010-2019 15.081590 27660.3132 27727.4310
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4576.93355 4614.26928 1.0738413 2000-2009 –> 2010-2019 2.285236 8717.5745 8790.9981
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 16596.82538 16637.83046 5.9734872 2000-2009 –> 2010-2019 15.075503 27659.9906 27738.2947
obs Atlantic (28.1,28.15] cstar_tref ~ temp + aou + nitrate + phosphate_star 16590.82104 16625.96825 5.9688641 2000-2009 –> 2010-2019 15.072332 27652.4722 27719.5900
obs Atlantic (28.1,28.15] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 16565.02380 16606.02887 5.9368701 2000-2009 –> 2010-2019 15.039592 27628.4253 27706.7295
mod Atlantic (28.1,28.15] cstar_tref ~ temp + aou + silicate + phosphate 4672.76304 4704.76510 1.1087040 2000-2009 –> 2010-2019 2.395060 8965.2301 9028.1645
obs Atlantic (28.1,28.15] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 16577.66358 16612.81079 5.9536987 2000-2009 –> 2010-2019 15.051928 27637.5613 27704.6791
obs Atlantic (28.1,28.15] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 16553.66452 16594.66959 5.9238453 2000-2009 –> 2010-2019 15.021477 27615.3623 27693.6665
mod Atlantic (28.1,28.15] cstar_tref ~ temp + silicate + phosphate + phosphate_star 4719.35944 4751.36150 1.1257049 2000-2009 –> 2010-2019 2.417568 9022.7709 9085.7053
mod Atlantic (28.15,28.2] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 7251.07917 7292.86952 0.8452295 2000-2009 –> 2010-2019 1.896195 13644.8834 13726.4576
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate 7309.30443 7345.12472 0.8540732 2000-2009 –> 2010-2019 1.896278 13664.7488 13734.6695
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate + silicate 7214.97836 7256.76870 0.8399722 2000-2009 –> 2010-2019 1.878341 13556.4032 13637.9773
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate 5567.88857 5601.59037 5.1002975 2000-2009 –> 2010-2019 12.846146 10191.7956 10256.9959
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5527.62342 5561.32522 4.9888202 2000-2009 –> 2010-2019 12.652575 10137.3594 10202.5597
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + nitrate + phosphate_star 7624.85397 7660.67426 0.9019451 2000-2009 –> 2010-2019 1.962347 14055.4895 14125.4102
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 7586.69137 7628.48172 0.8957061 2000-2009 –> 2010-2019 1.949690 13992.9556 14074.5298
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 5515.14642 5548.84822 4.9547737 2000-2009 –> 2010-2019 12.667655 10133.3809 10198.5812
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5544.23989 5577.94169 5.0345259 2000-2009 –> 2010-2019 12.709966 10156.0022 10221.2025
mod Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + phosphate 7317.34166 7353.16196 0.8552604 2000-2009 –> 2010-2019 1.896588 13669.1252 13739.0459
mod Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + phosphate_star 7350.59424 7386.41453 0.8601898 2000-2009 –> 2010-2019 1.907069 13725.4783 13795.3990
mod Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate 7243.90888 7279.72917 0.8444745 2000-2009 –> 2010-2019 1.895097 13634.2977 13704.2184
obs Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 5529.47311 5563.17491 4.9938875 2000-2009 –> 2010-2019 12.998558 10197.0949 10262.2952
obs Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5500.67440 5534.37620 4.9155741 2000-2009 –> 2010-2019 12.821555 10151.7967 10216.9970
obs Atlantic (28.15,28.2] cstar_tref ~ temp + aou + silicate + phosphate_star 5527.93510 5556.82235 4.9951539 2000-2009 –> 2010-2019 12.922240 10180.6031 10236.4891
mod Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 7401.38263 7437.20292 0.8677737 2000-2009 –> 2010-2019 1.918632 13792.7443 13862.6650
obs Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 5507.40714 5541.10894 4.9337720 2000-2009 –> 2010-2019 12.869155 10163.4665 10228.6668
mod Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + silicate + phosphate_star 7602.74729 7638.56759 0.8985056 2000-2009 –> 2010-2019 1.952636 14007.6161 14077.5368
obs Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + silicate + phosphate_star 5515.83705 5544.72431 4.9620960 2000-2009 –> 2010-2019 12.988273 10185.0274 10240.9134
obs Atlantic (28.15,28.2] cstar_tref ~ temp + silicate + phosphate + phosphate_star 5543.94021 5572.82747 5.0392264 2000-2009 –> 2010-2019 12.987553 10200.1671 10256.0531
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate 9942.69937 9975.09163 0.6794760 2000-2009 –> 2010-2019 1.430831 15787.9647 15849.6250
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + phosphate 9882.07876 9920.94947 0.6750673 2000-2009 –> 2010-2019 1.424315 15714.8826 15788.8749
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 10759.82339 10799.31770 3.1874909 2000-2009 –> 2010-2019 6.502013 18337.6815 18414.1069
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + phosphate_star 9798.25985 9837.13056 0.6692110 2000-2009 –> 2010-2019 1.418476 15631.1825 15705.1748
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate 9829.40615 9868.27687 0.6713812 2000-2009 –> 2010-2019 1.406490 15564.0929 15638.0853
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 10758.74571 10798.24002 3.1866669 2000-2009 –> 2010-2019 6.503458 18338.5815 18415.0069
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 9187.48710 9232.83627 0.6279130 2000-2009 –> 2010-2019 1.351076 14839.7961 14926.1205
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + phosphate + phosphate_star 10762.36640 10796.21866 3.1909671 2000-2009 –> 2010-2019 6.506904 18339.4576 18404.9650
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + silicate + phosphate 10757.86405 10791.71631 3.1875220 2000-2009 –> 2010-2019 6.504624 18335.9709 18401.4783
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 10738.84248 10778.33680 3.1714860 2000-2009 –> 2010-2019 6.485291 18316.0752 18392.5006
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate 9862.28569 9901.15640 0.6736798 2000-2009 –> 2010-2019 1.423559 15699.4294 15773.4217
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate 9800.00428 9845.35344 0.6691933 2000-2009 –> 2010-2019 1.418219 15633.2813 15719.6057
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate 10758.97651 10798.47082 3.1868433 2000-2009 –> 2010-2019 6.505318 18340.2786 18416.7040
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 9794.86201 9840.21117 0.6688357 2000-2009 –> 2010-2019 1.418075 15629.6080 15715.9324
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + silicate 9750.75775 9796.10692 0.6657763 2000-2009 –> 2010-2019 1.398547 15471.0419 15557.3663
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + silicate + phosphate 10730.22721 10769.72152 3.1649373 2000-2009 –> 2010-2019 6.482024 18310.3208 18386.7462
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 10778.16000 10817.65431 3.2015448 2000-2009 –> 2010-2019 6.476007 18320.8764 18397.3018
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 9840.76093 9886.11010 0.6720344 2000-2009 –> 2010-2019 1.424172 15695.3905 15781.7149
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 10803.73772 10843.23203 3.2212521 2000-2009 –> 2010-2019 6.503029 18352.9029 18429.3282
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 10749.50963 10789.00394 3.1796132 2000-2009 –> 2010-2019 6.461411 18298.6934 18375.1188
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 22523.15867 22566.52600 5.4012463 2000-2009 –> 2010-2019 10.940331 38034.4598 38118.5251
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 29319.78443 29364.93881 5.5475696 2000-2009 –> 2010-2019 11.845894 50310.1387 50397.8287
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + phosphate + phosphate_star 22868.97156 22906.14356 5.6667578 2000-2009 –> 2010-2019 11.561693 38686.4803 38758.5363
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 22659.19343 22702.56076 5.5035774 2000-2009 –> 2010-2019 11.157830 38272.3587 38356.4240
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 29664.50369 29709.65807 5.7557805 2000-2009 –> 2010-2019 12.277889 50879.5648 50967.2548
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate 30097.26990 30135.97365 6.0295591 2000-2009 –> 2010-2019 12.512904 51271.9645 51347.1273
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 22622.29787 22665.66520 5.4756329 2000-2009 –> 2010-2019 11.292010 38375.3972 38459.4625
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 29778.03433 29823.18872 5.8260497 2000-2009 –> 2010-2019 12.146710 50791.1722 50878.8622
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 22725.41507 22768.78241 5.5540914 2000-2009 –> 2010-2019 11.473250 38565.2229 38649.2882
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 29830.76423 29875.91861 5.8589777 2000-2009 –> 2010-2019 12.190299 50854.7478 50942.4378
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 23228.37698 23271.74432 5.9531950 2000-2009 –> 2010-2019 12.143974 39290.2743 39374.3396
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 30022.42489 30067.57927 5.9802385 2000-2009 –> 2010-2019 12.460340 51195.8992 51283.5892
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 23114.05656 23157.42389 5.8600337 2000-2009 –> 2010-2019 11.956885 39100.2761 39184.3414
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 30123.29336 30168.44774 6.0450612 2000-2009 –> 2010-2019 12.661666 51430.9341 51518.6241
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 23315.97773 23359.34507 6.0255830 2000-2009 –> 2010-2019 12.299005 39443.5174 39527.5827
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 30233.66296 30278.81734 6.1167950 2000-2009 –> 2010-2019 12.769884 51576.6946 51664.3846
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 22827.01241 22870.37975 5.6324929 2000-2009 –> 2010-2019 11.643649 38743.1629 38827.2282
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 29958.88286 30004.03724 5.9397608 2000-2009 –> 2010-2019 12.315361 51027.7206 51115.4106
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 23383.18639 23426.55373 6.0817164 2000-2009 –> 2010-2019 12.391548 39539.3718 39623.4371
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 30281.54263 30326.69701 6.1481782 2000-2009 –> 2010-2019 12.816619 51639.4083 51727.0983
mod Indo-Pacific (26,26.5] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 16318.33448 16360.08113 4.1230694 2000-2009 –> 2010-2019 8.572240 30547.0401 30629.3879
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 16575.15216 16616.89881 4.3113964 2000-2009 –> 2010-2019 8.580683 30602.3726 30684.7204
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 15846.71932 15888.46597 3.7983922 2000-2009 –> 2010-2019 7.865389 29636.9583 29719.3061
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 16273.42889 16315.17554 4.0909951 2000-2009 –> 2010-2019 8.370937 30312.8193 30395.1671
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 17566.84088 17608.78439 4.7071211 2000-2009 –> 2010-2019 9.280822 32253.6329 32336.3336
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 17587.33887 17629.28238 4.7234644 2000-2009 –> 2010-2019 9.314334 32292.8348 32375.5354
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + phosphate 17838.72527 17874.67686 4.9302401 2000-2009 –> 2010-2019 9.794082 32830.5548 32901.4410
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 15477.71124 15519.45790 3.5622863 2000-2009 –> 2010-2019 7.641249 29282.2938 29364.6416
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 17422.67886 17464.62237 4.5937657 2000-2009 –> 2010-2019 9.270987 32221.1992 32303.8999
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 15629.37645 15671.12311 3.6574973 2000-2009 –> 2010-2019 7.905447 29632.1368 29714.4846
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 17494.03172 17535.97523 4.6495255 2000-2009 –> 2010-2019 9.363041 32331.1398 32413.8404
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 17590.54848 17632.49199 4.7260286 2000-2009 –> 2010-2019 9.322334 32301.9518 32384.6525
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate_star 17638.49868 17674.45026 4.7661141 2000-2009 –> 2010-2019 9.534107 32530.9715 32601.8578
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 17405.59482 17447.53834 4.5805146 2000-2009 –> 2010-2019 9.294536 32243.2379 32325.9386
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 15672.58431 15714.33097 3.6850847 2000-2009 –> 2010-2019 7.965328 29712.3180 29794.6658
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 17493.27872 17535.22224 4.6489335 2000-2009 –> 2010-2019 9.364211 32332.2527 32414.9533
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate_star 17492.90568 17528.85727 4.6502127 2000-2009 –> 2010-2019 9.385084 32350.5801 32421.4664
mod Indo-Pacific (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 16082.32916 16124.07582 3.9572664 2000-2009 –> 2010-2019 8.235554 30119.8317 30202.1795
mod Indo-Pacific (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 16149.86727 16191.61393 4.0040216 2000-2009 –> 2010-2019 8.292337 30198.7995 30281.1473
mod Indo-Pacific (26,26.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 16197.21348 16238.96014 4.0371273 2000-2009 –> 2010-2019 8.342999 30266.0913 30348.4391
obs Indo-Pacific (26.5,26.75] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 15862.48298 15897.44164 5.7177081 2000-2009 –> 2010-2019 11.099347 32451.2052 32521.5073
mod Indo-Pacific (26.5,26.75] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 14259.76826 14300.92331 3.5861572 2000-2009 –> 2010-2019 7.923291 27088.2844 27169.3763
mod Indo-Pacific (26.5,26.75] cstar_tref ~ aou + silicate + phosphate + phosphate_star 14302.13142 14337.40717 3.6163924 2000-2009 –> 2010-2019 8.232945 27405.8565 27475.3638
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 15699.00128 15739.78638 5.5320085 2000-2009 –> 2010-2019 10.559829 31926.2966 32008.3158
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 15725.06151 15765.84661 5.5608474 2000-2009 –> 2010-2019 10.646791 32013.7814 32095.8006
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + phosphate 15905.27330 15940.23196 5.7667324 2000-2009 –> 2010-2019 11.142686 32488.3463 32558.6484
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 14078.48505 14119.64009 3.4652100 2000-2009 –> 2010-2019 7.832149 26937.4088 27018.5007
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 15852.78461 15893.56971 5.7043781 2000-2009 –> 2010-2019 11.071285 32428.8570 32510.8762
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 14249.55541 14290.71045 3.5792326 2000-2009 –> 2010-2019 8.059475 27222.2092 27303.3011
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 15727.64243 15768.42753 5.5637117 2000-2009 –> 2010-2019 10.658930 32026.0985 32108.1177
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate_star 15747.57266 15782.53132 5.5881093 2000-2009 –> 2010-2019 10.795829 32160.7396 32231.0418
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 14756.63488 14797.78992 3.9397351 2000-2009 –> 2010-2019 8.553380 27859.5628 27940.6547
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 15733.33618 15774.12128 5.5700358 2000-2009 –> 2010-2019 10.709535 32078.0351 32160.0543
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 14274.66248 14315.81752 3.5962799 2000-2009 –> 2010-2019 8.085069 27255.7776 27336.8695
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 15890.01307 15930.79817 5.7469072 2000-2009 –> 2010-2019 11.143950 32496.0095 32578.0286
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 14433.77627 14474.93131 3.7062193 2000-2009 –> 2010-2019 8.137375 27357.5153 27438.6072
obs Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 15846.97299 15887.75810 5.6977675 2000-2009 –> 2010-2019 11.117201 32475.0936 32557.1128
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 14650.80587 14691.96091 3.8616142 2000-2009 –> 2010-2019 8.389095 27670.0285 27751.1204
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + silicate + phosphate 14479.37883 14514.65458 3.7397588 2000-2009 –> 2010-2019 8.490807 27710.6077 27780.1151
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 14693.04186 14734.19690 3.8926045 2000-2009 –> 2010-2019 8.436846 27728.6714 27809.7633
mod Indo-Pacific (26.75,27] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 18504.60473 18541.30860 3.8171531 2000-2009 –> 2010-2019 7.744540 36865.6393 36938.9387
obs Indo-Pacific (26.75,27] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 24861.56265 24899.35388 5.3299483 2000-2009 –> 2010-2019 10.371248 49451.8229 49527.4485
mod Indo-Pacific (26.75,27] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 17528.79366 17571.61485 3.2991010 2000-2009 –> 2010-2019 7.219253 35879.6879 35965.2039
obs Indo-Pacific (26.75,27] cstar_tref ~ aou + phosphate + phosphate_star 25174.16169 25205.65438 5.5427471 2000-2009 –> 2010-2019 10.762082 50043.3320 50106.3534
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + phosphate 17816.51710 17859.33828 3.4437748 2000-2009 –> 2010-2019 7.144811 35788.7158 35874.2317
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + phosphate 25036.85284 25080.94261 5.4461328 2000-2009 –> 2010-2019 10.584670 49783.7446 49871.9745
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 17997.18713 18040.00832 3.5378450 2000-2009 –> 2010-2019 7.318474 36109.4777 36194.9936
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate 25289.35394 25333.44371 5.6199739 2000-2009 –> 2010-2019 10.784768 50077.4996 50165.7295
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 18019.92654 18062.74773 3.5498654 2000-2009 –> 2010-2019 7.334240 36138.7366 36224.2525
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + phosphate_star 18519.15108 18555.85496 3.8254446 2000-2009 –> 2010-2019 7.632743 36675.7229 36749.0223
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 17964.79555 18007.61673 3.5207925 2000-2009 –> 2010-2019 7.323994 36116.2783 36201.7942
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 25171.73348 25215.82326 5.5383151 2000-2009 –> 2010-2019 10.625645 49837.5616 49925.7915
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + phosphate 25107.35150 25145.14274 5.4954886 2000-2009 –> 2010-2019 10.696150 49949.5128 50025.1384
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 17523.82229 17566.64348 3.2966555 2000-2009 –> 2010-2019 7.226018 35890.1681 35975.6841
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 25079.04646 25123.13624 5.4748033 2000-2009 –> 2010-2019 10.560485 49742.2517 49830.4816
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 17714.37950 17757.20068 3.3917054 2000-2009 –> 2010-2019 7.420250 36244.8499 36330.3659
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 25222.80944 25266.89922 5.5736282 2000-2009 –> 2010-2019 10.738312 50010.7820 50099.0119
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 17743.14594 17785.96713 3.4062903 2000-2009 –> 2010-2019 7.440007 36282.0639 36367.5798
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 25223.32912 25267.41890 5.5739886 2000-2009 –> 2010-2019 10.737319 50009.1802 50097.4101
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate_star 25221.67564 25259.46687 5.5742290 2000-2009 –> 2010-2019 10.739004 50007.7910 50083.4166
mod Indo-Pacific (27,27.25] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 21274.41083 21312.78776 2.6668076 2000-2009 –> 2010-2019 5.848682 43221.9560 43298.4709
mod Indo-Pacific (27,27.25] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 21222.32351 21267.09660 2.6505772 2000-2009 –> 2010-2019 5.662027 42703.1839 42792.4512
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 31798.95389 31844.94372 4.9339487 2000-2009 –> 2010-2019 9.687206 62960.4140 63052.3389
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 32008.13030 32054.12012 5.0328268 2000-2009 –> 2010-2019 9.894185 63404.8116 63496.7365
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate 31799.81890 31839.23875 4.9352898 2000-2009 –> 2010-2019 9.689609 62961.6158 63040.4087
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate_star 32012.80215 32052.22200 5.0360130 2000-2009 –> 2010-2019 9.916563 63448.6979 63527.4907
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 31935.40964 31981.39947 4.9982289 2000-2009 –> 2010-2019 9.818237 63242.7393 63334.6643
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + silicate + phosphate + phosphate_star 31941.73475 31981.15460 5.0021776 2000-2009 –> 2010-2019 9.844046 63294.3965 63373.1894
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate 21171.86094 21216.63403 2.6355236 2000-2009 –> 2010-2019 5.800514 43076.1077 43165.3750
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 21569.25313 21614.02622 2.7564243 2000-2009 –> 2010-2019 6.017871 43729.0980 43818.3653
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate 31799.64295 31845.63277 4.9342712 2000-2009 –> 2010-2019 9.686239 62958.2644 63050.1894
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 31977.12070 32023.11052 5.0180443 2000-2009 –> 2010-2019 9.828996 63264.7791 63356.7041
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 21597.34174 21642.11483 2.7651767 2000-2009 –> 2010-2019 6.034462 43777.6258 43866.8931
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 22077.31098 22122.08407 2.9191053 2000-2009 –> 2010-2019 6.160692 44185.1524 44274.4197
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 32065.63555 32111.62538 5.0603552 2000-2009 –> 2010-2019 9.913122 63443.8147 63535.7397
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 31932.49494 31978.48477 4.9968472 2000-2009 –> 2010-2019 9.796602 63195.7808 63287.7058
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + phosphate 21596.17738 21634.55431 2.7654375 2000-2009 –> 2010-2019 5.993312 43665.9293 43742.4441
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 21596.09871 21640.87179 2.7647888 2000-2009 –> 2010-2019 5.884499 43377.6615 43466.9289
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 21953.73906 21998.51214 2.8786746 2000-2009 –> 2010-2019 6.117032 44053.0970 44142.3643
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 21974.28804 22019.06113 2.8853589 2000-2009 –> 2010-2019 6.127475 44083.5206 44172.7880
obs Indo-Pacific (27.25,27.5] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 33860.73695 33900.82200 4.2842975 2000-2009 –> 2010-2019 8.648821 66609.4572 66689.3883
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 33858.36703 33905.13292 4.2827083 2000-2009 –> 2010-2019 8.590358 66460.6357 66553.8887
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 29440.90513 29487.41438 3.2311139 2000-2009 –> 2010-2019 6.297742 57135.6430 57228.3759
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 33566.96910 33613.73499 4.1780504 2000-2009 –> 2010-2019 8.417129 65987.6219 66080.8748
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + silicate + phosphate 33571.70837 33611.79341 4.1804417 2000-2009 –> 2010-2019 8.430542 66019.7516 66099.6827
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 33815.38627 33862.15216 4.2671081 2000-2009 –> 2010-2019 8.601135 66486.7466 66579.9995
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + silicate + phosphate + phosphate_star 33816.03893 33856.12398 4.2680693 2000-2009 –> 2010-2019 8.645310 66597.6903 66677.6214
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 29026.51844 29073.02769 3.1153140 2000-2009 –> 2010-2019 6.300769 57135.6348 57228.3676
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 29284.69116 29331.20041 3.1869630 2000-2009 –> 2010-2019 6.441490 57627.6325 57720.3653
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 33546.72814 33593.49403 4.1708764 2000-2009 –> 2010-2019 8.365751 65848.7411 65941.9940
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 33760.76670 33807.53259 4.2473655 2000-2009 –> 2010-2019 8.494413 66202.6754 66295.9283
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 29303.84676 29350.35601 3.1923444 2000-2009 –> 2010-2019 6.452918 57667.0221 57759.7550
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 28907.59202 28954.10127 3.0828533 2000-2009 –> 2010-2019 6.133982 56547.0988 56639.8316
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 33701.84624 33748.61213 4.2261708 2000-2009 –> 2010-2019 8.477950 66156.3592 66249.6122
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 33707.52357 33754.28946 4.2282084 2000-2009 –> 2010-2019 8.465910 66124.5001 66217.7530
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + phosphate 29026.86079 29066.72586 3.1159568 2000-2009 –> 2010-2019 6.366147 57353.2665 57432.7518
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 28694.32043 28740.82968 3.0254859 2000-2009 –> 2010-2019 6.157987 56620.7176 56713.4504
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 28972.16135 29018.67060 3.1004352 2000-2009 –> 2010-2019 6.321073 57201.0073 57293.7401
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 28973.07299 29019.58224 3.1006841 2000-2009 –> 2010-2019 6.322437 57205.6926 57298.4255
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + nitrate + silicate + phosphate_star 28971.15262 29011.01769 3.1007058 2000-2009 –> 2010-2019 6.323412 57204.9949 57284.4802
obs Indo-Pacific (27.5,27.75] cstar_tref ~ aou + silicate + phosphate + phosphate_star 38454.14084 38495.13287 4.0030523 2000-2009 –> 2010-2019 7.672834 72677.4883 72758.9616
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 38982.42151 39030.24554 4.1598200 2000-2009 –> 2010-2019 7.800675 73108.1984 73203.2505
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 32980.58838 33027.89286 3.2309441 2000-2009 –> 2010-2019 6.144101 65141.8230 65236.5412
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 38227.49441 38275.31843 3.9367977 2000-2009 –> 2010-2019 7.775291 73018.3760 73113.4282
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 33151.50001 33198.80450 3.2746495 2000-2009 –> 2010-2019 6.188157 65314.2891 65409.0073
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate_star 33418.68574 33459.23244 3.3446874 2000-2009 –> 2010-2019 6.384815 66129.1121 66210.2991
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 33105.14009 33152.44458 3.2627363 2000-2009 –> 2010-2019 6.171912 65248.7055 65343.4237
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 33346.25584 33393.56032 3.3251734 2000-2009 –> 2010-2019 6.329201 65904.3466 65999.0648
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 38015.89884 38063.72287 3.8764613 2000-2009 –> 2010-2019 7.460506 71943.8054 72038.8576
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 33404.77532 33452.07981 3.3405064 2000-2009 –> 2010-2019 6.358351 66022.1551 66116.8732
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 38114.13409 38161.95811 3.9043572 2000-2009 –> 2010-2019 7.521277 72156.9247 72251.9768
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 33399.49128 33446.79576 3.3391190 2000-2009 –> 2010-2019 6.344756 65964.5028 66059.2210
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 33094.96401 33142.26849 3.2601271 2000-2009 –> 2010-2019 6.223176 65475.5964 65570.3145
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 38078.34311 38126.16713 3.8941704 2000-2009 –> 2010-2019 7.505507 72101.6976 72196.7498
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate 33439.16789 33479.71460 3.3500775 2000-2009 –> 2010-2019 6.386511 66133.8850 66215.0720
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 33375.75830 33423.06278 3.3328947 2000-2009 –> 2010-2019 6.343249 65961.0307 66055.7489
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 38099.85140 38147.67542 3.9002889 2000-2009 –> 2010-2019 7.515800 72137.7411 72232.7932
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate_star 38392.22574 38433.21776 3.9850019 2000-2009 –> 2010-2019 7.610156 72461.6293 72543.1026
obs Indo-Pacific (27.5,27.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 38527.88984 38575.71386 4.0240719 2000-2009 –> 2010-2019 7.739813 72909.8335 73004.8857
obs Indo-Pacific (27.5,27.75] cstar_tref ~ temp + aou + silicate + phosphate 38529.06440 38570.05642 4.0250045 2000-2009 –> 2010-2019 7.748030 72933.6503 73015.1236
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + phosphate 11582.49612 11623.22231 2.4811612 2000-2009 –> 2010-2019 4.714277 22708.1773 22789.6718
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 11622.21369 11662.93989 2.5010688 2000-2009 –> 2010-2019 4.742110 22765.6094 22847.1039
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + phosphate 11906.57917 11941.48733 2.6494099 2000-2009 –> 2010-2019 4.931195 23138.0629 23207.9153
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 11664.57598 11705.30218 2.5224780 2000-2009 –> 2010-2019 4.803856 22897.1671 22978.6616
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 21265.73264 21309.46318 3.9152850 2000-2009 –> 2010-2019 7.005747 37663.7026 37749.9644
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 11704.06600 11744.79220 2.5426006 2000-2009 –> 2010-2019 4.831384 22952.8605 23034.3550
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 21245.21474 21288.94528 3.9047760 2000-2009 –> 2010-2019 6.996741 37646.3122 37732.5740
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 11647.20477 11687.93097 2.5136768 2000-2009 –> 2010-2019 4.760324 22803.0930 22884.5875
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 21243.36271 21287.09325 3.9038288 2000-2009 –> 2010-2019 7.004558 37662.6656 37748.9274
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 11720.79167 11761.51787 2.5511717 2000-2009 –> 2010-2019 4.842183 22974.4510 23055.9455
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 21244.66567 21288.39621 3.9044951 2000-2009 –> 2010-2019 6.996864 37646.6039 37732.8657
mod Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 11573.66152 11614.38772 2.4767547 2000-2009 –> 2010-2019 4.873942 23053.8350 23135.3295
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 21269.79440 21313.52494 3.9173687 2000-2009 –> 2010-2019 7.007638 37667.3639 37753.6257
mod Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 11623.25316 11663.97936 2.5015919 2000-2009 –> 2010-2019 4.907466 23121.5126 23203.0072
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 21251.75549 21295.48603 3.9081230 2000-2009 –> 2010-2019 7.000080 37652.8362 37739.0980
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + silicate + phosphate_star 21256.00311 21293.48643 3.9113227 2000-2009 –> 2010-2019 7.003305 37655.1382 37729.0768
mod Indo-Pacific (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 11657.15572 11697.88191 2.5187147 2000-2009 –> 2010-2019 4.930529 23167.7449 23249.2394
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 21251.90058 21295.63112 3.9081973 2000-2009 –> 2010-2019 7.000542 37653.7896 37740.0514
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate_star 21249.97904 21287.46236 3.9082374 2000-2009 –> 2010-2019 7.008969 37667.2885 37741.2272
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + silicate + phosphate + phosphate_star 21255.08352 21292.56684 3.9108516 2000-2009 –> 2010-2019 7.003235 37655.0513 37728.9900
obs Indo-Pacific (27.85,27.95] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 27306.42126 27352.33901 3.3094478 2000-2009 –> 2010-2019 6.053395 48709.2055 48799.8552
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 13743.08420 13785.52123 2.1053564 2000-2009 –> 2010-2019 4.122815 26435.0887 26519.5446
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 27332.24803 27378.16577 3.3176497 2000-2009 –> 2010-2019 6.062718 48738.6291 48829.2787
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate 13180.79924 13217.17384 1.9274450 2000-2009 –> 2010-2019 3.879512 25673.8303 25746.2211
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + phosphate 13182.77368 13225.21071 1.9274372 2000-2009 –> 2010-2019 3.874618 25662.8233 25747.2792
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 13181.91640 13224.35343 1.9271769 2000-2009 –> 2010-2019 3.875732 25666.1848 25750.6407
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + silicate 13167.23798 13209.67501 1.9227244 2000-2009 –> 2010-2019 3.843234 25564.8376 25649.2935
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + phosphate 13795.26325 13831.63784 2.1234078 2000-2009 –> 2010-2019 4.176723 26590.5851 26662.9759
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 13655.58322 13698.02025 2.0765262 2000-2009 –> 2010-2019 4.123631 26434.7959 26519.2518
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 27288.91418 27334.83192 3.3038996 2000-2009 –> 2010-2019 6.052377 48706.2276 48796.8773
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 13671.85548 13714.29251 2.0818576 2000-2009 –> 2010-2019 4.134493 26467.1968 26551.6526
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 27281.46433 27327.38207 3.3015414 2000-2009 –> 2010-2019 6.061120 48734.2831 48824.9327
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 13224.12618 13266.56321 1.9400380 2000-2009 –> 2010-2019 3.897681 25736.2097 25820.6656
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 13687.15827 13729.59530 2.0868839 2000-2009 –> 2010-2019 4.141735 26488.9479 26573.4038
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 27281.47835 27327.39610 3.3015459 2000-2009 –> 2010-2019 6.062798 48739.6363 48830.2859
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 27271.97695 27317.89470 3.2985408 2000-2009 –> 2010-2019 6.046695 48688.2552 48778.9048
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 27270.15117 27316.06892 3.2979637 2000-2009 –> 2010-2019 6.057669 48723.3738 48814.0234
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + silicate + phosphate 27292.83026 27332.18833 3.3057734 2000-2009 –> 2010-2019 6.054920 48710.2871 48787.9868
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + silicate + phosphate_star 27287.26124 27326.61931 3.3040095 2000-2009 –> 2010-2019 6.064165 48739.9210 48817.6206
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 27270.17866 27316.09641 3.2979724 2000-2009 –> 2010-2019 6.058949 48727.4576 48818.1072
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 46570.02549 46612.50004 3.4399361 2000-2009 –> 2010-2019 6.291242 83276.8857 83360.8464
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 46570.12556 46619.67921 3.4395635 2000-2009 –> 2010-2019 6.278476 83214.1864 83312.1405
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + phosphate + phosphate_star 46725.68751 46761.08297 3.4709960 2000-2009 –> 2010-2019 6.326228 83451.0197 83520.9869
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + silicate + phosphate + phosphate_star 46724.80068 46767.27523 3.4704248 2000-2009 –> 2010-2019 6.318408 83414.3167 83498.2774
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 46572.00523 46621.55887 3.4399321 2000-2009 –> 2010-2019 6.291210 83280.7217 83378.6758
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + aou + phosphate + phosphate_star 46713.30916 46755.78371 3.4681519 2000-2009 –> 2010-2019 6.323231 83439.8416 83523.8023
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate 17729.68183 17768.14215 1.7388062 2000-2009 –> 2010-2019 3.464592 30519.5753 30594.5577
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate + silicate 17594.14007 17639.01045 1.7123885 2000-2009 –> 2010-2019 3.391332 30207.0596 30294.5391
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate 18028.43935 18066.89967 1.7976014 2000-2009 –> 2010-2019 3.624391 31188.2647 31263.2472
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate_star 18030.90379 18069.36411 1.7980945 2000-2009 –> 2010-2019 3.629036 31205.4944 31280.4769
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 18021.84001 18066.71039 1.7958816 2000-2009 –> 2010-2019 3.592722 31076.1556 31163.6351
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 46812.21359 46861.76724 3.4873656 2000-2009 –> 2010-2019 6.345155 83554.8617 83652.8159
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 18025.37597 18070.24635 1.7965885 2000-2009 –> 2010-2019 3.599260 31100.7640 31188.2436
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 18238.03747 18282.90785 1.8396231 2000-2009 –> 2010-2019 3.681077 31451.8659 31539.3454
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + phosphate + phosphate_star 18037.24229 18075.70260 1.7993636 2000-2009 –> 2010-2019 3.632323 31218.9983 31293.9807
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 18031.47302 18076.34339 1.7978082 2000-2009 –> 2010-2019 3.602013 31112.3886 31199.8681
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate_star 18287.18738 18325.64770 1.8501268 2000-2009 –> 2010-2019 3.708692 31559.1725 31634.1549
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 46660.65814 46710.21179 3.4573627 2000-2009 –> 2010-2019 6.309443 83373.5596 83471.5138
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 46775.58952 46825.14317 3.4800915 2000-2009 –> 2010-2019 6.360223 83634.0883 83732.0425
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 46760.10551 46809.65915 3.4770207 2000-2009 –> 2010-2019 6.355298 83609.0227 83706.9768
obs Indo-Pacific (28.05,28.1] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 35866.63840 35907.89348 2.9621831 2000-2009 –> 2010-2019 5.690486 58728.7732 58808.7806
obs Indo-Pacific (28.05,28.1] cstar_tref ~ aou + phosphate + phosphate_star 35940.11636 35974.49559 2.9778440 2000-2009 –> 2010-2019 5.707806 58805.9852 58872.6581
obs Indo-Pacific (28.05,28.1] cstar_tref ~ aou + silicate + phosphate + phosphate_star 35915.36968 35956.62476 2.9722849 2000-2009 –> 2010-2019 5.701269 58779.8589 58859.8662
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 35883.52586 35931.65678 2.9652655 2000-2009 –> 2010-2019 5.705818 58789.9152 58883.2572
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + phosphate 35878.93401 35927.06494 2.9643144 2000-2009 –> 2010-2019 5.718955 58833.6865 58927.0285
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + silicate 9299.69593 9340.69830 1.4580999 2000-2009 –> 2010-2019 3.129643 16746.8040 16826.7378
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 9480.47802 9521.48039 1.5099880 2000-2009 –> 2010-2019 3.088429 16707.1752 16787.1091
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 35877.80798 35925.93891 2.9640813 2000-2009 –> 2010-2019 5.676834 58688.0309 58781.3729
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 9485.65453 9526.65689 1.5115006 2000-2009 –> 2010-2019 3.094363 16723.1085 16803.0424
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 35997.52177 36045.65269 2.9889750 2000-2009 –> 2010-2019 5.724540 58886.7292 58980.0711
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 9437.79841 9478.80078 1.4975740 2000-2009 –> 2010-2019 3.080415 16675.2013 16755.1352
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 35967.33774 36015.46866 2.9826788 2000-2009 –> 2010-2019 5.712385 58836.3233 58929.6653
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 9485.06336 9526.06573 1.5113278 2000-2009 –> 2010-2019 3.094190 16722.5170 16802.4509
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 35998.92413 36047.05505 2.9892679 2000-2009 –> 2010-2019 5.724470 58886.8802 58980.2222
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate_star 9490.30094 9525.44582 1.5134450 2000-2009 –> 2010-2019 3.096309 16725.7585 16794.2732
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate_star 35997.02427 36038.27935 2.9892888 2000-2009 –> 2010-2019 5.724922 58884.4681 58964.4755
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate 9892.33679 9927.48167 1.6358323 2000-2009 –> 2010-2019 3.239044 17176.9194 17245.4342
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate_star 9900.75011 9935.89500 1.6384966 2000-2009 –> 2010-2019 3.247630 17199.5137 17268.0285
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + nitrate + silicate + phosphate_star 9892.04100 9927.18589 1.6357388 2000-2009 –> 2010-2019 3.244408 17189.6958 17258.2105
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + silicate + phosphate + phosphate_star 9898.74180 9933.88668 1.6378602 2000-2009 –> 2010-2019 3.246970 17197.4495 17265.9642
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 18999.77354 19046.14079 1.3330832 2000-2009 –> 2010-2019 2.774821 36476.1392 36567.9739
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 37470.70856 37519.37136 2.7356957 2000-2009 –> 2010-2019 5.269841 66999.6646 67095.5467
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + phosphate + phosphate_star 37532.88195 37574.59292 2.7470868 2000-2009 –> 2010-2019 5.286362 67085.2504 67167.4349
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + silicate + phosphate 37504.59301 37546.30399 2.7420595 2000-2009 –> 2010-2019 5.278063 67040.7628 67122.9474
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + nitrate + silicate + phosphate 19087.28553 19127.02889 1.3438515 2000-2009 –> 2010-2019 2.806761 36704.2771 36782.9926
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 19089.15673 19135.52397 1.3438360 2000-2009 –> 2010-2019 2.800653 36667.3227 36759.1575
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 37586.25861 37634.92141 2.7562406 2000-2009 –> 2010-2019 5.325065 67286.0103 67381.8924
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + nitrate + phosphate 37625.60684 37674.26964 2.7632719 2000-2009 –> 2010-2019 5.337740 67352.9385 67448.8205
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + nitrate + silicate 18149.04592 18195.41317 1.2349512 2000-2009 –> 2010-2019 2.757373 36158.1784 36250.0131
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate 17728.80573 17775.17298 1.1891758 2000-2009 –> 2010-2019 2.564311 34742.4167 34834.2514
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate 37488.38169 37537.04449 2.7388280 2000-2009 –> 2010-2019 5.256120 66933.4896 67029.3716
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 17747.33733 17793.70458 1.1911581 2000-2009 –> 2010-2019 2.568189 34774.4268 34866.2615
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 37589.82924 37638.49204 2.7568779 2000-2009 –> 2010-2019 5.294166 67134.3606 67230.2426
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 18562.00834 18608.37559 1.2816500 2000-2009 –> 2010-2019 2.720629 36019.6280 36111.4627
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 17709.33657 17755.70382 1.1870967 2000-2009 –> 2010-2019 2.574447 34809.4681 34901.3028
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 37666.65824 37715.32104 2.7706267 2000-2009 –> 2010-2019 5.325193 67296.4694 67392.3515
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 17743.75629 17790.12354 1.1907748 2000-2009 –> 2010-2019 2.567538 34768.9423 34860.7770
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 37573.51327 37622.17607 2.7539669 2000-2009 –> 2010-2019 5.287226 67098.0747 67193.9567
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate_star 17931.37647 17971.11982 1.2112431 2000-2009 –> 2010-2019 2.601711 35051.4753 35130.1908
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate_star 37697.87807 37739.58904 2.7765927 2000-2009 –> 2010-2019 5.339151 67364.9448 67447.1294
obs Indo-Pacific (28.2,28.3] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 11277.66930 11312.97449 2.0191092 2000-2009 –> 2010-2019 4.629172 19528.1450 19596.1924
obs Indo-Pacific (28.2,28.3] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 11199.54616 11240.73556 1.9888714 2000-2009 –> 2010-2019 4.596952 19449.3900 19528.7786
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 11219.10336 11260.29275 1.9962101 2000-2009 –> 2010-2019 4.603773 19468.2591 19547.6478
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 31807.11665 31860.44313 0.6965855 2000-2009 –> 2010-2019 1.350861 54023.0800 54127.6471
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 34910.91367 34964.24016 0.7723335 2000-2009 –> 2010-2019 1.541645 60741.9905 60846.5576
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + phosphate + phosphate_star 11349.55162 11384.85682 2.0466281 2000-2009 –> 2010-2019 4.655826 19598.8797 19666.9271
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 11340.83700 11382.02640 2.0425025 2000-2009 –> 2010-2019 4.651566 19591.9856 19671.3742
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + nitrate + silicate + phosphate 32249.09693 32294.80534 0.7069475 2000-2009 –> 2010-2019 1.453828 57417.7148 57507.3437
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 32076.66931 32129.99579 0.7028583 2000-2009 –> 2010-2019 1.362771 54484.1188 54588.6859
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + nitrate + phosphate 11203.75420 11244.94360 1.9904481 2000-2009 –> 2010-2019 4.648621 19519.4989 19598.8875
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 11334.83319 11376.02259 2.0401944 2000-2009 –> 2010-2019 4.648427 19584.8795 19664.2681
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + silicate + phosphate 34862.40003 34915.72652 0.7710884 2000-2009 –> 2010-2019 1.601991 62412.4848 62517.0519
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 34874.54021 34927.86670 0.7713998 2000-2009 –> 2010-2019 1.602605 62432.7508 62537.3179
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 32062.80080 32116.12729 0.7025342 2000-2009 –> 2010-2019 1.441572 56997.8143 57102.3815
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 34872.15494 34925.48142 0.7713386 2000-2009 –> 2010-2019 1.602505 62429.3430 62533.9101
obs Indo-Pacific (28.2,28.3] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 11236.97082 11278.16022 2.0029384 2000-2009 –> 2010-2019 4.612926 19489.3466 19568.7352
mod Indo-Pacific (28.2,28.3] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 36079.73289 36133.05938 0.8029473 2000-2009 –> 2010-2019 1.640440 63806.1608 63910.7279
obs Indo-Pacific (28.2,28.3] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 11207.84438 11249.03378 1.9919819 2000-2009 –> 2010-2019 4.639946 19510.2601 19589.6487
obs Indo-Pacific (28.2,28.3] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 11349.53631 11390.72571 2.0458515 2000-2009 –> 2010-2019 4.655401 19601.3307 19680.7194
mod Indo-Pacific (28.2,28.3] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 36078.23156 36131.55805 0.8029072 2000-2009 –> 2010-2019 1.641650 63837.9485 63942.5156
mod Indo-Pacific (28.3,28.4] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 20703.78864 20744.24887 1.2603476 2000-2009 –> 2010-2019 2.189653 33631.0593 33710.3766
mod Indo-Pacific (28.3,28.4] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 20704.97612 20752.17973 1.2602659 2000-2009 –> 2010-2019 2.187210 33609.8257 33702.3625
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 20391.45252 20438.65612 1.2291426 2000-2009 –> 2010-2019 2.157906 33315.1235 33407.6603
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 4492.68896 4527.37810 2.0458300 2000-2009 –> 2010-2019 4.694073 7770.5158 7836.8801
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 4476.27609 4510.96524 2.0298877 2000-2009 –> 2010-2019 4.664392 7747.0084 7813.3727
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + silicate + phosphate_star 4496.50920 4526.24276 2.0515134 2000-2009 –> 2010-2019 4.702071 7773.5279 7830.4116
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 4524.29254 4558.98169 2.0768810 2000-2009 –> 2010-2019 4.710318 7794.4723 7860.8367
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + silicate + phosphate + phosphate_star 4525.43895 4555.17251 2.0799981 2000-2009 –> 2010-2019 4.714907 7794.3808 7851.2646
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + nitrate + phosphate 20368.36020 20415.56380 1.2268809 2000-2009 –> 2010-2019 2.152932 33263.9600 33356.4967
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 20450.82267 20498.02628 1.2349767 2000-2009 –> 2010-2019 2.165992 33397.7350 33490.2717
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + silicate + phosphate 4502.56600 4537.25515 2.0554842 2000-2009 –> 2010-2019 4.689012 7772.7925 7839.1569
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4476.10508 4510.79423 2.0297223 2000-2009 –> 2010-2019 4.679434 7754.6884 7821.0527
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 20366.43728 20413.64088 1.2266928 2000-2009 –> 2010-2019 2.153813 33273.1124 33365.6491
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 4502.63595 4537.32509 2.0555527 2000-2009 –> 2010-2019 4.705553 7781.3679 7847.7322
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4496.11451 4530.80366 2.0491731 2000-2009 –> 2010-2019 4.683536 7766.7734 7833.1378
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + aou + nitrate + phosphate 20684.53392 20724.99415 1.2584136 2000-2009 –> 2010-2019 2.185006 33583.7387 33663.0559
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 20686.20735 20733.41096 1.2583808 2000-2009 –> 2010-2019 2.182284 33559.5168 33652.0535
obs Indo-Pacific (28.3,28.4] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 4473.47431 4508.16345 2.0271787 2000-2009 –> 2010-2019 4.712574 7770.3038 7836.6682
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 20687.37747 20727.83770 1.2586990 2000-2009 –> 2010-2019 2.186418 33598.2431 33677.5603
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 20689.00960 20736.21320 1.2586621 2000-2009 –> 2010-2019 2.183751 33574.6298 33667.1666
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
c6b3da6 jens-daniel-mueller 2021-06-14
439ee80 jens-daniel-mueller 2021-06-11
33ffcab jens-daniel-mueller 2021-06-10
7e1f407 jens-daniel-mueller 2021-06-10
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
d3c3fde jens-daniel-mueller 2021-03-26
594d503 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26

8.1.2 Best models

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

[[2]]

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

8.2 RMSE correlation between eras

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

8.2.1 All models

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

8.2.2 Best models

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

8.3 Predictor counts

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

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

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

8.4.2 AIC vs BIC

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

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

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

8.5 RMSE vs VIF

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

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

8.6 Residual patterns

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

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

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

8.6.3 By model domain

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

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

[[2]]

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

[[3]]

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

[[4]]

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

8.6.4 Latitudinal mean

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

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

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

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

[[2]]

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

8.6.5 Latitudinal offset

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

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