Last updated: 2021-06-14

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: aeddb8d

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 aeddb8d. 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 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 132.95178 146.47456 0.7766662 2000-2009 –> 2010-2019 4.827059 789.3993 822.0755
obs Atlantic (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 850.02359 869.87756 6.6766435 2000-2009 –> 2010-2019 14.634375 1276.2073 1310.6041
mod Atlantic (-Inf,26] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 80.22693 93.74971 0.4631729 2000-2009 –> 2010-2019 4.847410 754.7325 787.4086
mod Atlantic (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 110.14181 123.66459 0.6210326 2000-2009 –> 2010-2019 4.606182 762.8869 795.5630
obs Atlantic (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 842.77014 862.62411 6.4872056 2000-2009 –> 2010-2019 15.399135 1282.3170 1316.7137
mod Atlantic (-Inf,26] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 77.59713 91.11991 0.4513838 2000-2009 –> 2010-2019 4.324098 723.8170 756.4932
mod Atlantic (-Inf,26] cstar_tref ~ sal + nitrate + silicate + phosphate_star 84.83213 96.42308 0.4941589 2000-2009 –> 2010-2019 4.891141 757.9995 786.0077
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate 852.24553 869.26322 6.7894444 2000-2009 –> 2010-2019 14.842766 1277.8383 1307.3212
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 850.56896 870.42293 6.6911085 2000-2009 –> 2010-2019 14.680480 1277.2209 1311.6177
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 850.61284 870.46681 6.6922737 2000-2009 –> 2010-2019 14.678466 1277.2178 1311.6146
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 78.70244 92.22522 0.4563018 2000-2009 –> 2010-2019 4.748782 748.3856 781.0617
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 845.37941 865.23338 6.5547246 2000-2009 –> 2010-2019 14.540764 1271.9822 1306.3789
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 109.46911 122.99188 0.6169502 2000-2009 –> 2010-2019 4.494156 755.9532 788.6294
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 132.85392 146.37670 0.7759214 2000-2009 –> 2010-2019 4.812061 788.4977 821.1739
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 850.80320 870.65718 6.6973310 2000-2009 –> 2010-2019 14.707104 1277.7561 1312.1529
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 81.74180 95.26458 0.4701031 2000-2009 –> 2010-2019 4.778228 752.2544 784.9306
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 109.37850 122.90128 0.6164025 2000-2009 –> 2010-2019 4.803795 773.4103 806.0865
obs Atlantic (-Inf,26] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 856.51571 876.36968 6.8508845 2000-2009 –> 2010-2019 15.342240 1290.3582 1324.7549
obs Atlantic (-Inf,26] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 855.62688 875.48085 6.8267633 2000-2009 –> 2010-2019 15.309954 1289.3559 1323.7526
obs Atlantic (-Inf,26] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 852.96084 872.81481 6.7549201 2000-2009 –> 2010-2019 15.202001 1286.1865 1320.5832
obs Atlantic (26,26.5] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 3973.14361 3999.69283 5.9958310 2000-2009 –> 2010-2019 11.465193 7227.9682 7280.0404
mod Atlantic (26,26.5] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 2087.85084 2115.33933 3.8427795 2000-2009 –> 2010-2019 8.313155 3449.2327 3500.8181
obs Atlantic (26,26.5] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 3974.55108 4005.52516 5.9929527 2000-2009 –> 2010-2019 11.439584 7227.0444 7287.7953
mod Atlantic (26,26.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 2051.71785 2079.20634 3.6620335 2000-2009 –> 2010-2019 7.668410 3362.4713 3414.0567
obs Atlantic (26,26.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 3853.86255 3884.83664 5.4345761 2000-2009 –> 2010-2019 10.892337 7108.4788 7169.2297
mod Atlantic (26,26.5] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2016.39028 2043.87877 3.4935387 2000-2009 –> 2010-2019 7.884711 3369.5134 3421.0988
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 2055.54710 2083.03558 3.6807784 2000-2009 –> 2010-2019 7.968407 3397.6457 3449.2311
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 3870.80054 3901.77462 5.5096859 2000-2009 –> 2010-2019 11.036366 7138.4674 7199.2183
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 2055.82076 2083.30925 3.6821218 2000-2009 –> 2010-2019 8.039528 3405.3776 3456.9630
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 3871.95137 3902.92545 5.5148266 2000-2009 –> 2010-2019 11.115352 7153.4223 7214.1731
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 2020.79873 2048.28721 3.5141339 2000-2009 –> 2010-2019 8.252558 3409.0839 3460.6693
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2021.54686 2049.03534 3.5176411 2000-2009 –> 2010-2019 8.373347 3421.1278 3472.7132
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 2056.87335 2084.36183 3.6872930 2000-2009 –> 2010-2019 8.117270 3414.0611 3465.6465
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 3882.25381 3913.22790 5.5610616 2000-2009 –> 2010-2019 11.222574 7174.9885 7235.7394
obs Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + phosphate 3971.99505 3998.54426 5.9902529 2000-2009 –> 2010-2019 11.551418 7244.1310 7296.2032
mod Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 2085.26031 2112.74880 3.8295293 2000-2009 –> 2010-2019 8.340023 3450.7698 3502.3552
obs Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 3973.47957 4004.45365 5.9877511 2000-2009 –> 2010-2019 11.470324 7232.8132 7293.5641
mod Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2081.23749 2108.72597 3.8090435 2000-2009 –> 2010-2019 8.355953 3450.4619 3502.0473
obs Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 3973.67194 4004.64602 5.9886846 2000-2009 –> 2010-2019 11.514093 7241.0994 7301.8503
obs Atlantic (26,26.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 3974.58217 4005.55626 5.9931037 2000-2009 –> 2010-2019 11.573529 7252.3139 7313.0648
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 2726.45851 2755.71863 4.0107952 2000-2009 –> 2010-2019 7.178969 4628.3755 4684.9732
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 4781.08845 4813.68553 5.1796785 2000-2009 –> 2010-2019 10.851025 8050.3932 8112.7131
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 2729.52408 2758.78420 4.0235436 2000-2009 –> 2010-2019 7.082252 4605.6317 4662.2293
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 4780.02262 4812.61971 5.1761317 2000-2009 –> 2010-2019 10.876657 8054.6235 8116.9433
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 2794.15487 2823.41499 4.3019507 2000-2009 –> 2010-2019 7.154576 4619.0638 4675.6614
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 4868.55787 4901.15496 5.4791899 2000-2009 –> 2010-2019 11.177059 8142.6777 8204.9975
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2800.36688 2829.62699 4.3297042 2000-2009 –> 2010-2019 7.128337 4611.2500 4667.8476
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 2736.69312 2765.95323 4.0535148 2000-2009 –> 2010-2019 6.978815 4580.0674 4636.6651
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 4777.16617 4809.76325 5.1666383 2000-2009 –> 2010-2019 10.865809 8051.5217 8113.8415
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 2801.37501 2830.63513 4.3342251 2000-2009 –> 2010-2019 7.128920 4611.2245 4667.8222
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + phosphate 4837.47167 4865.41203 5.3777192 2000-2009 –> 2010-2019 11.053221 8105.5322 8158.9492
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + phosphate_star 4834.36936 4862.30972 5.3670079 2000-2009 –> 2010-2019 11.067915 8107.0392 8160.4562
mod Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 2752.33082 2781.59094 4.1196675 2000-2009 –> 2010-2019 7.187804 4630.6973 4687.2950
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 4837.25390 4869.85099 5.3700598 2000-2009 –> 2010-2019 11.044839 8107.1831 8169.5030
mod Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2756.04217 2785.30229 4.1355256 2000-2009 –> 2010-2019 7.103566 4610.0629 4666.6606
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 4834.31961 4866.91670 5.3599425 2000-2009 –> 2010-2019 11.053032 8107.5734 8169.8932
mod Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 2763.22631 2788.30641 4.1750312 2000-2009 –> 2010-2019 7.239945 4638.8214 4687.3337
obs Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 4831.97713 4859.91749 5.3587629 2000-2009 –> 2010-2019 11.058986 8104.5233 8157.9403
mod Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 2760.32020 2789.58032 4.1538809 2000-2009 –> 2010-2019 7.073818 4602.3476 4658.9453
obs Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 4831.89617 4864.49325 5.3516010 2000-2009 –> 2010-2019 11.044603 8105.1341 8167.4539
mod Atlantic (26.75,27] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 4240.27178 4273.99660 2.4425103 2000-2009 –> 2010-2019 4.770226 7088.5817 7153.3820
mod Atlantic (26.75,27] cstar_tref ~ aou + silicate + phosphate + phosphate_star 4641.90307 4670.81005 3.0460059 2000-2009 –> 2010-2019 5.408673 7506.8728 7562.4159
mod Atlantic (26.75,27] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 4390.47272 4424.19754 2.6516793 2000-2009 –> 2010-2019 5.181129 7342.8418 7407.6421
mod Atlantic (26.75,27] cstar_tref ~ sal + aou + silicate + phosphate 4558.21119 4587.11818 2.9096939 2000-2009 –> 2010-2019 5.445584 7511.7640 7567.3071
mod Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate 4541.25047 4574.97528 2.8796694 2000-2009 –> 2010-2019 5.155058 7361.0946 7425.8949
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate 5656.49344 5690.67051 4.3694968 2000-2009 –> 2010-2019 10.747946 9681.7981 9746.9037
mod Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4560.77146 4594.49627 2.9105859 2000-2009 –> 2010-2019 5.260043 7420.7211 7485.5213
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5660.07838 5694.25545 4.3775372 2000-2009 –> 2010-2019 10.699856 9674.5464 9739.6521
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 5668.40994 5702.58700 4.3962806 2000-2009 –> 2010-2019 10.704890 9680.2167 9745.3223
mod Atlantic (26.75,27] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4567.89619 4601.62100 2.9219522 2000-2009 –> 2010-2019 5.304794 7445.5113 7510.3116
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5660.99451 5695.17158 4.3795943 2000-2009 –> 2010-2019 10.698786 9674.8562 9739.9618
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + silicate + phosphate_star 5677.62041 5706.91504 4.4216275 2000-2009 –> 2010-2019 10.827400 9706.1656 9761.9704
mod Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 4256.68746 4290.41228 2.4645432 2000-2009 –> 2010-2019 4.740250 7076.7062 7141.5065
obs Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 5660.55536 5694.73242 4.3786081 2000-2009 –> 2010-2019 10.847584 9703.1379 9768.2436
mod Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 4265.99116 4299.71597 2.4771186 2000-2009 –> 2010-2019 4.769763 7095.2943 7160.0946
obs Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5665.11355 5699.29062 4.3888552 2000-2009 –> 2010-2019 10.770355 9691.0044 9756.1100
mod Atlantic (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 4288.30980 4322.03461 2.5075479 2000-2009 –> 2010-2019 4.833398 7135.6162 7200.4165
obs Atlantic (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 5666.86320 5701.04026 4.3927949 2000-2009 –> 2010-2019 10.785847 9694.9716 9760.0772
obs Atlantic (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate_star 5670.13688 5699.43150 4.4046911 2000-2009 –> 2010-2019 10.828605 9702.1493 9757.9541
obs Atlantic (26.75,27] cstar_tref ~ temp + silicate + phosphate + phosphate_star 5666.26719 5695.56181 4.3959588 2000-2009 –> 2010-2019 10.847377 9703.5177 9759.3225
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 5403.90244 5439.88561 2.0473003 2000-2009 –> 2010-2019 4.365737 9390.6862 9460.1209
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 5547.81109 5583.79426 2.1674213 2000-2009 –> 2010-2019 4.646337 9652.2549 9721.6895
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate 5694.40925 5725.25197 2.2988574 2000-2009 –> 2010-2019 4.762309 9785.8515 9845.3669
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 5696.38150 5732.36467 2.2988321 2000-2009 –> 2010-2019 4.643704 9703.0973 9772.5320
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate 11056.43901 11094.54875 6.1094992 2000-2009 –> 2010-2019 12.006897 19028.5157 19102.5193
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 11029.77887 11067.88861 6.0620586 2000-2009 –> 2010-2019 11.845430 18953.2006 19027.2042
mod Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + silicate 5533.21222 5569.19539 2.1549211 2000-2009 –> 2010-2019 4.419886 9478.9756 9548.4102
mod Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate 5692.73169 5728.71486 2.2955103 2000-2009 –> 2010-2019 4.675084 9725.2731 9794.7077
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 11081.53915 11119.64889 6.1545031 2000-2009 –> 2010-2019 11.893589 18985.8054 19059.8090
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 11036.98050 11075.09024 6.0748371 2000-2009 –> 2010-2019 11.891285 18974.6142 19048.6177
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 11101.08725 11139.19699 6.1897818 2000-2009 –> 2010-2019 12.024453 19046.5161 19120.5197
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 11083.26215 11121.37189 6.1576045 2000-2009 –> 2010-2019 11.909980 18993.2921 19067.2957
obs Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + phosphate_star 11053.97426 11086.63975 6.1086690 2000-2009 –> 2010-2019 12.019799 19029.8468 19093.2785
mod Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate 5702.20474 5733.04745 2.3059685 2000-2009 –> 2010-2019 4.608935 9675.2187 9734.7341
mod Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 5704.09286 5740.07604 2.3058663 2000-2009 –> 2010-2019 4.566853 9646.7651 9716.1998
obs Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 11020.08413 11058.19387 6.0448987 2000-2009 –> 2010-2019 12.019820 19024.7055 19098.7091
mod Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5703.24179 5739.22497 2.3050889 2000-2009 –> 2010-2019 4.564899 9644.9992 9714.4339
obs Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 10978.01041 11016.12015 5.9709883 2000-2009 –> 2010-2019 11.813664 18926.8558 19000.8594
mod Atlantic (27,27.25] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 5607.64067 5643.62384 2.2194122 2000-2009 –> 2010-2019 4.481351 9551.0537 9620.4884
obs Atlantic (27,27.25] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 10989.63952 11027.74926 5.9913262 2000-2009 –> 2010-2019 11.868210 18953.0327 19027.0363
obs Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 13160.06856 13199.51569 5.7917839 2000-2009 –> 2010-2019 10.472589 21853.9263 21930.4007
mod Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate 11047.31894 11081.06958 3.5748251 2000-2009 –> 2010-2019 6.902963 20716.3773 20783.2396
mod Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 10902.58531 10941.96106 3.4490889 2000-2009 –> 2010-2019 6.739279 20531.3972 20609.4032
mod Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 10912.11643 10951.49218 3.4571201 2000-2009 –> 2010-2019 6.769148 20565.2992 20643.3052
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 13144.56163 13184.00875 5.7701306 2000-2009 –> 2010-2019 10.372339 21788.8031 21865.2775
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 11160.19972 11199.57547 3.6728707 2000-2009 –> 2010-2019 6.938218 20761.0899 20839.0959
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 13142.47457 13181.92170 5.7672225 2000-2009 –> 2010-2019 10.344423 21770.7512 21847.2256
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + silicate 11021.54018 11060.91593 3.5506747 2000-2009 –> 2010-2019 6.798877 20603.0365 20681.0425
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + phosphate 13211.88473 13245.69656 5.8675632 2000-2009 –> 2010-2019 10.484423 21863.4390 21928.9884
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + phosphate_star 13214.88928 13248.70111 5.8718231 2000-2009 –> 2010-2019 10.466045 21852.0411 21917.5906
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 13209.31604 13248.76316 5.8610916 2000-2009 –> 2010-2019 10.476000 21861.6316 21938.1060
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 10887.45031 10926.82606 3.4363740 2000-2009 –> 2010-2019 7.101064 20913.3923 20991.3983
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 13212.13059 13251.57772 5.8650776 2000-2009 –> 2010-2019 10.459191 21851.2133 21927.6877
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 11156.38306 11195.75881 3.6694515 2000-2009 –> 2010-2019 6.802058 20604.3849 20682.3909
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 13139.70144 13179.14857 5.7633607 2000-2009 –> 2010-2019 10.345694 21771.2617 21847.7361
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 10866.78034 10906.15609 3.4190848 2000-2009 –> 2010-2019 6.956681 20762.6900 20840.6960
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + phosphate + phosphate_star 13211.76150 13245.57332 5.8673886 2000-2009 –> 2010-2019 10.465426 21851.3459 21916.8953
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 10888.81762 10928.19337 3.4375207 2000-2009 –> 2010-2019 7.104642 20917.2022 20995.2082
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 13209.18141 13248.62854 5.8609010 2000-2009 –> 2010-2019 10.458461 21850.4617 21926.9361
mod Atlantic (27.25,27.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 11184.37255 11223.74830 3.6945998 2000-2009 –> 2010-2019 6.820533 20624.5180 20702.5240
obs Atlantic (27.5,27.75] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 14123.06674 14163.31199 5.0622376 2000-2009 –> 2010-2019 10.333021 25247.5768 25326.3025
obs Atlantic (27.5,27.75] cstar_tref ~ aou + silicate + phosphate + phosphate_star 14133.97107 14168.46700 5.0763359 2000-2009 –> 2010-2019 10.348263 25257.2636 25324.7428
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + nitrate + phosphate 18254.93278 18292.91796 2.1793032 2000-2009 –> 2010-2019 4.614678 33880.6554 33955.3962
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 17796.33340 17840.64945 2.0616596 2000-2009 –> 2010-2019 4.490299 33405.3281 33492.5257
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate 18532.59316 18570.57834 2.2534407 2000-2009 –> 2010-2019 4.695389 34176.5446 34251.2854
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 18500.86236 18545.17840 2.2443014 2000-2009 –> 2010-2019 4.679059 34126.8703 34214.0678
obs Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 14135.81855 14176.06381 5.0761690 2000-2009 –> 2010-2019 10.272305 25208.8935 25287.6192
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 18145.18546 18189.50151 2.1501587 2000-2009 –> 2010-2019 4.569858 33729.2445 33816.4420
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 18531.11210 18575.42815 2.2524958 2000-2009 –> 2010-2019 4.674786 34122.4065 34209.6041
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 14141.02773 14181.27299 5.0818711 2000-2009 –> 2010-2019 10.285367 25219.2067 25297.9324
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 18346.55989 18390.87594 2.2029638 2000-2009 –> 2010-2019 4.582892 33818.5515 33905.7490
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 14169.16511 14209.41036 5.1127816 2000-2009 –> 2010-2019 10.343897 25266.4335 25345.1592
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 18464.65283 18508.96888 2.2345317 2000-2009 –> 2010-2019 4.624014 33963.7350 34050.9326
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 18454.93003 18499.24607 2.2319157 2000-2009 –> 2010-2019 4.619478 33948.5782 34035.7757
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 18351.41366 18395.72971 2.2042524 2000-2009 –> 2010-2019 4.585946 33828.4194 33915.6170
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 14149.70083 14189.94609 5.0913790 2000-2009 –> 2010-2019 10.307462 25236.5922 25315.3179
obs Atlantic (27.5,27.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 14129.35475 14169.60001 5.0691025 2000-2009 –> 2010-2019 10.356871 25265.4665 25344.1922
obs Atlantic (27.5,27.75] cstar_tref ~ temp + aou + silicate + phosphate 14139.15247 14173.64840 5.0820077 2000-2009 –> 2010-2019 10.370382 25273.6772 25341.1564
obs Atlantic (27.5,27.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 14138.83884 14179.08410 5.0794743 2000-2009 –> 2010-2019 10.391237 25291.2764 25370.0021
obs Atlantic (27.5,27.75] cstar_tref ~ temp + silicate + phosphate + phosphate_star 14148.13754 14182.63348 5.0918582 2000-2009 –> 2010-2019 10.403790 25298.6900 25366.1692
obs Atlantic (27.75,27.85] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 9233.49201 9270.63303 5.3491009 2000-2009 –> 2010-2019 10.526246 16960.6102 17033.7177
obs Atlantic (27.75,27.85] cstar_tref ~ aou + silicate + phosphate + phosphate_star 9231.49608 9263.33124 5.3491082 2000-2009 –> 2010-2019 10.551074 16968.6573 17031.3209
mod Atlantic (27.75,27.85] cstar_tref ~ sal + aou + nitrate + silicate 5173.73538 5205.70600 1.3173948 2000-2009 –> 2010-2019 2.554887 9445.3221 9508.3364
obs Atlantic (27.75,27.85] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 9219.15570 9256.29672 5.3234119 2000-2009 –> 2010-2019 10.516654 16954.0906 17027.1982
mod Atlantic (27.75,27.85] cstar_tref ~ sal + aou + silicate + phosphate 4952.21222 4984.18284 1.2249870 2000-2009 –> 2010-2019 2.411417 9113.8187 9176.8330
obs Atlantic (27.75,27.85] cstar_tref ~ sal + aou + silicate + phosphate 9217.17947 9249.01463 5.3234543 2000-2009 –> 2010-2019 10.544523 16963.5708 17026.2344
obs Atlantic (27.75,27.85] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 9220.38728 9257.52830 5.3256138 2000-2009 –> 2010-2019 10.550236 16970.4915 17043.5991
mod Atlantic (27.75,27.85] cstar_tref ~ sal + nitrate + silicate + phosphate_star 5577.28913 5609.25975 1.5040216 2000-2009 –> 2010-2019 2.869667 10106.0650 10169.0793
mod Atlantic (27.75,27.85] cstar_tref ~ sal + silicate + phosphate + phosphate_star 5568.68143 5600.65205 1.4997774 2000-2009 –> 2010-2019 2.871297 10108.6597 10171.6741
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + silicate 5129.95023 5167.24929 1.2977408 2000-2009 –> 2010-2019 2.531791 9396.2672 9469.7839
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 4831.72900 4869.02806 1.1767060 2000-2009 –> 2010-2019 2.358419 8984.9375 9058.4543
obs Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 9217.89443 9255.03545 5.3211577 2000-2009 –> 2010-2019 10.531377 16961.0479 17034.1554
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4962.49905 4999.79811 1.2283242 2000-2009 –> 2010-2019 2.449161 9200.7199 9274.2366
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 4967.78126 5005.08032 1.2304561 2000-2009 –> 2010-2019 2.449911 9203.0454 9276.5621
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4956.33627 4993.63534 1.2258415 2000-2009 –> 2010-2019 2.442704 9186.0458 9259.5626
obs Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 9219.61953 9256.76055 5.3242410 2000-2009 –> 2010-2019 10.543976 16967.3676 17040.4751
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate_star 5047.16381 5079.13444 1.2637744 2000-2009 –> 2010-2019 2.568050 9455.9345 9518.9488
obs Atlantic (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 9233.50098 9270.64200 5.3491170 2000-2009 –> 2010-2019 10.533199 16963.9908 17037.0983
obs Atlantic (27.75,27.85] cstar_tref ~ temp + aou + silicate + phosphate 9231.50412 9263.33928 5.3491227 2000-2009 –> 2010-2019 10.559358 16972.6650 17035.3286
obs Atlantic (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 9233.83358 9270.97460 5.3497145 2000-2009 –> 2010-2019 10.541639 16968.1299 17041.2374
mod Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate 6272.84959 6304.72891 1.9503548 2000-2009 –> 2010-2019 4.869347 12779.2310 12842.1494
obs Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 19473.23877 19515.62932 5.3008631 2000-2009 –> 2010-2019 10.571644 36143.7573 36227.4626
mod Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 6273.69009 6310.88264 1.9496012 2000-2009 –> 2010-2019 4.790961 12711.7711 12785.1760
obs Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 18865.02567 18907.41622 4.8133300 2000-2009 –> 2010-2019 9.850525 35290.4957 35374.2010
mod Atlantic (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate 6039.14869 6071.02801 1.8041887 2000-2009 –> 2010-2019 4.538793 12375.3541 12438.2726
obs Atlantic (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate_star 18868.38656 18904.72132 4.8174249 2000-2009 –> 2010-2019 9.854895 35292.1512 35363.8987
obs Atlantic (27.85,27.95] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 19270.25518 19312.64572 5.1328983 2000-2009 –> 2010-2019 10.326655 35861.1912 35944.8965
obs Atlantic (27.85,27.95] cstar_tref ~ sal + silicate + phosphate + phosphate_star 19309.46175 19345.79651 5.1665598 2000-2009 –> 2010-2019 10.363841 35902.0649 35973.8124
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + silicate 6112.43748 6149.63002 1.8475750 2000-2009 –> 2010-2019 4.753813 12609.3983 12682.8032
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + phosphate 6515.42325 6547.30257 2.1146075 2000-2009 –> 2010-2019 4.866437 12868.0047 12930.9232
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 5733.93669 5771.12923 1.6285777 2000-2009 –> 2010-2019 4.362944 12071.9150 12145.3199
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 18703.14289 18745.53344 4.6913002 2000-2009 –> 2010-2019 9.697729 35095.4925 35179.1979
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5840.07759 5877.27014 1.6872287 2000-2009 –> 2010-2019 4.446416 12201.6230 12275.0278
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 18682.65392 18725.04447 4.6760775 2000-2009 –> 2010-2019 9.677745 35069.8596 35153.5649
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 18686.57177 18728.96231 4.6789845 2000-2009 –> 2010-2019 9.724328 35120.7790 35204.4843
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 5849.02844 5886.22098 1.6922702 2000-2009 –> 2010-2019 4.495710 12252.0694 12325.4743
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate 18708.33204 18744.66679 4.6966532 2000-2009 –> 2010-2019 9.743371 35142.0117 35213.7592
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5829.55620 5866.74874 1.6813217 2000-2009 –> 2010-2019 4.431976 12183.0232 12256.4281
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 18701.74999 18744.14053 4.6902637 2000-2009 –> 2010-2019 9.695537 35092.8518 35176.5572
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate_star 5934.11537 5965.99469 1.7421151 2000-2009 –> 2010-2019 4.774666 12540.0331 12602.9516
obs Atlantic (27.95,28.05] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 19920.09774 19962.61593 5.3746838 2000-2009 –> 2010-2019 12.081948 37687.5589 37771.3112
obs Atlantic (27.95,28.05] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 19499.97690 19542.49509 5.0342289 2000-2009 –> 2010-2019 11.596347 37150.5237 37234.2760
obs Atlantic (27.95,28.05] cstar_tref ~ sal + aou + silicate + phosphate_star 19584.24414 19620.68829 5.1023317 2000-2009 –> 2010-2019 11.678776 37244.4445 37316.2322
obs Atlantic (27.95,28.05] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 19637.82822 19680.34640 5.1434936 2000-2009 –> 2010-2019 11.750497 37324.8037 37408.5559
obs Atlantic (27.95,28.05] cstar_tref ~ sal + silicate + phosphate + phosphate_star 19764.23854 19800.68270 5.2474075 2000-2009 –> 2010-2019 11.868294 37460.4311 37532.2187
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate 9605.11085 9638.03053 3.5598329 2000-2009 –> 2010-2019 7.360215 17449.6013 17514.0798
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate_star 9634.96172 9667.88140 3.5897404 2000-2009 –> 2010-2019 7.388944 17478.5695 17543.0481
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 9468.62186 9507.02815 3.4243081 2000-2009 –> 2010-2019 7.223830 17314.4680 17389.6930
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 19617.19352 19659.71171 5.1269883 2000-2009 –> 2010-2019 11.726150 37297.8225 37381.5748
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 9480.28832 9518.69462 3.4355230 2000-2009 –> 2010-2019 7.232315 17324.0906 17399.3157
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 19561.88847 19604.40665 5.0830115 2000-2009 –> 2010-2019 11.658914 37223.6486 37307.4009
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate + phosphate_star 9698.71714 9731.63682 3.6544608 2000-2009 –> 2010-2019 7.440569 17532.5058 17596.9844
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 19653.70185 19696.22004 5.1562267 2000-2009 –> 2010-2019 11.871545 37427.5764 37511.3287
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 9453.77314 9492.17944 3.4100870 2000-2009 –> 2010-2019 7.160190 17262.3864 17337.6114
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + phosphate + phosphate_star 9630.23201 9663.15169 3.5849850 2000-2009 –> 2010-2019 7.384187 17473.8387 17538.3173
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + phosphate_star 9711.16674 9738.59981 3.6692906 2000-2009 –> 2010-2019 7.471381 17554.9351 17608.6672
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate 19754.11109 19790.55525 5.2391363 2000-2009 –> 2010-2019 11.958946 37529.5585 37601.3461
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 9480.55493 9518.96122 3.4357798 2000-2009 –> 2010-2019 7.232577 17324.3611 17399.5862
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 19601.08730 19643.60549 5.1141420 2000-2009 –> 2010-2019 11.704626 37274.6841 37358.4364
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate_star 9479.29520 9512.21488 3.4364927 2000-2009 –> 2010-2019 7.233298 17321.1073 17385.5859
obs Atlantic (28.05,28.1] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 14397.87547 14438.16883 5.2582693 2000-2009 –> 2010-2019 11.595670 25991.0093 26069.6657
obs Atlantic (28.05,28.1] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 14340.80802 14381.10139 5.1944313 2000-2009 –> 2010-2019 11.531629 25933.8285 26012.4848
obs Atlantic (28.05,28.1] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 14344.31895 14384.61232 5.1983363 2000-2009 –> 2010-2019 11.522788 25930.2003 26008.8567
obs Atlantic (28.05,28.1] cstar_tref ~ sal + silicate + phosphate + phosphate_star 14429.72026 14464.25743 5.2964994 2000-2009 –> 2010-2019 11.644459 26026.7576 26094.1773
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 14466.74453 14507.03790 5.3363545 2000-2009 –> 2010-2019 11.641552 26041.8139 26120.4703
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate 3550.30794 3579.98402 1.3281783 2000-2009 –> 2010-2019 2.322033 5718.2020 5775.7016
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate_star 3589.40608 3619.08216 1.3534049 2000-2009 –> 2010-2019 2.365641 5785.2647 5842.7643
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 3551.49221 3586.11431 1.3276570 2000-2009 –> 2010-2019 2.302904 5692.5443 5759.6272
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 3588.59054 3623.21263 1.3515724 2000-2009 –> 2010-2019 2.349908 5765.3486 5832.4316
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 14377.47729 14417.77066 5.2353615 2000-2009 –> 2010-2019 11.506783 25933.5001 26012.1565
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 14459.86768 14500.16104 5.3285056 2000-2009 –> 2010-2019 11.636984 26036.7819 26115.4383
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + phosphate_star 3635.40586 3665.08195 1.3836987 2000-2009 –> 2010-2019 2.406936 5847.7606 5905.2602
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + phosphate_star 14468.05893 14502.59610 5.3401416 2000-2009 –> 2010-2019 11.649048 26043.2135 26110.6332
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 3610.65269 3645.27479 1.3659985 2000-2009 –> 2010-2019 2.388413 5823.7797 5890.8626
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 14346.87803 14387.17139 5.2011845 2000-2009 –> 2010-2019 11.492265 25913.9992 25992.6555
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + phosphate + phosphate_star 3578.07643 3607.75252 1.3460460 2000-2009 –> 2010-2019 2.353879 5767.2828 5824.7825
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 3578.66624 3613.28833 1.3451328 2000-2009 –> 2010-2019 2.336986 5745.4830 5812.5659
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 14419.17831 14459.47168 5.2823001 2000-2009 –> 2010-2019 11.579720 25989.8710 26068.5273
mod Atlantic (28.05,28.1] cstar_tref ~ temp + aou + phosphate 3569.79879 3594.52886 1.3419857 2000-2009 –> 2010-2019 2.419813 5859.4687 5907.3851
mod Atlantic (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate 3571.58217 3601.25825 1.3418458 2000-2009 –> 2010-2019 2.395166 5828.1544 5885.6541
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 16563.27184 16604.27691 5.9348594 2000-2009 –> 2010-2019 15.023969 27622.1151 27700.4192
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + phosphate 4576.07324 4608.07530 1.0742411 2000-2009 –> 2010-2019 2.400969 8947.7119 9010.6464
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + phosphate_star 4623.21453 4655.21659 1.0909076 2000-2009 –> 2010-2019 2.425934 9010.8268 9073.7612
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + silicate + phosphate 4541.52262 4578.85835 1.0614943 2000-2009 –> 2010-2019 2.269671 8675.3490 8748.7725
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4598.30703 4635.64277 1.0813632 2000-2009 –> 2010-2019 2.303492 8761.5506 8834.9741
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 16596.61665 16637.62172 5.9732461 2000-2009 –> 2010-2019 15.078607 27660.9011 27739.2053
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 16546.34993 16587.35500 5.9154733 2000-2009 –> 2010-2019 14.998434 27603.1318 27681.4359
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + phosphate_star 4639.83303 4671.83508 1.0968444 2000-2009 –> 2010-2019 2.406382 8978.0580 9040.9924
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 4641.28089 4678.61662 1.0966467 2000-2009 –> 2010-2019 2.375977 8921.7153 8995.1388
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 16570.19116 16611.19623 5.9428046 2000-2009 –> 2010-2019 15.051913 27635.7291 27714.0332
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + phosphate + phosphate_star 4608.15461 4640.15667 1.0855554 2000-2009 –> 2010-2019 2.414923 8984.8844 9047.8188
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + phosphate + phosphate_star 16598.50891 16633.65611 5.9777431 2000-2009 –> 2010-2019 15.081540 27660.2703 27727.3881
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4577.19302 4614.52876 1.0739323 2000-2009 –> 2010-2019 2.285327 8717.8340 8791.2575
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 16596.76398 16637.76905 5.9734163 2000-2009 –> 2010-2019 15.075432 27659.9292 27738.2333
obs Atlantic (28.1,28.15] cstar_tref ~ temp + aou + nitrate + phosphate_star 16590.80452 16625.95173 5.9688450 2000-2009 –> 2010-2019 15.072313 27652.4557 27719.5735
obs Atlantic (28.1,28.15] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 16564.95564 16605.96072 5.9367919 2000-2009 –> 2010-2019 15.039514 27628.3572 27706.6613
mod Atlantic (28.1,28.15] cstar_tref ~ temp + aou + silicate + phosphate 4672.98644 4704.98850 1.1087849 2000-2009 –> 2010-2019 2.395141 8965.4535 9028.3879
obs Atlantic (28.1,28.15] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 16577.64806 16612.79527 5.9536809 2000-2009 –> 2010-2019 15.051910 27637.5458 27704.6636
obs Atlantic (28.1,28.15] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 16553.60041 16594.60548 5.9237718 2000-2009 –> 2010-2019 15.021403 27615.2982 27693.6023
mod Atlantic (28.1,28.15] cstar_tref ~ temp + silicate + phosphate + phosphate_star 4719.57828 4751.58034 1.1257854 2000-2009 –> 2010-2019 2.417648 9022.9897 9085.9241
mod Atlantic (28.15,28.2] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 7251.41773 7293.20807 0.8452789 2000-2009 –> 2010-2019 1.896244 13645.2220 13726.7962
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate 7309.65193 7345.47222 0.8541245 2000-2009 –> 2010-2019 1.896330 13665.0963 13735.0170
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate + silicate 7215.31024 7257.10059 0.8400204 2000-2009 –> 2010-2019 1.878389 13556.7351 13638.3092
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate 5567.92808 5601.62988 5.1004081 2000-2009 –> 2010-2019 12.846256 10191.8351 10257.0354
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5527.65810 5561.35990 4.9889152 2000-2009 –> 2010-2019 12.652670 10137.3940 10202.5943
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + nitrate + phosphate_star 7625.11460 7660.93489 0.9019858 2000-2009 –> 2010-2019 1.962388 14055.7501 14125.6708
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 7586.97765 7628.76800 0.8957504 2000-2009 –> 2010-2019 1.949734 13993.2419 14074.8161
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 5515.21434 5548.91614 4.9549584 2000-2009 –> 2010-2019 12.667840 10133.4488 10198.6491
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5544.28144 5577.98324 5.0346407 2000-2009 –> 2010-2019 12.710081 10156.0438 10221.2441
mod Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + phosphate 7317.66502 7353.48531 0.8553082 2000-2009 –> 2010-2019 1.896635 13669.4485 13739.3693
mod Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + phosphate_star 7350.91497 7386.73526 0.8602375 2000-2009 –> 2010-2019 1.907117 13725.7990 13795.7197
mod Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate 7244.22005 7280.04035 0.8445200 2000-2009 –> 2010-2019 1.895143 13634.6089 13704.5296
obs Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 5529.51769 5563.21949 4.9940096 2000-2009 –> 2010-2019 12.998681 10197.1395 10262.3398
obs Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5500.71160 5534.41340 4.9156745 2000-2009 –> 2010-2019 12.821656 10151.8339 10217.0342
obs Atlantic (28.15,28.2] cstar_tref ~ temp + aou + silicate + phosphate_star 5527.96408 5556.85134 4.9952333 2000-2009 –> 2010-2019 12.922320 10180.6321 10236.5181
mod Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 7401.69776 7437.51806 0.8678209 2000-2009 –> 2010-2019 1.918679 13793.0594 13862.9802
obs Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 5507.45097 5541.15277 4.9338907 2000-2009 –> 2010-2019 12.869274 10163.5103 10228.7106
mod Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + silicate + phosphate_star 7603.04155 7638.86185 0.8985513 2000-2009 –> 2010-2019 1.952682 14007.9104 14077.8311
obs Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + silicate + phosphate_star 5515.89837 5544.78563 4.9622630 2000-2009 –> 2010-2019 12.988440 10185.0887 10240.9747
obs Atlantic (28.15,28.2] cstar_tref ~ temp + silicate + phosphate + phosphate_star 5543.97677 5572.86403 5.0393275 2000-2009 –> 2010-2019 12.987654 10200.2037 10256.0896
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate 9943.48230 9975.87457 0.6795313 2000-2009 –> 2010-2019 1.430886 15788.7477 15850.4080
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + phosphate 9882.84828 9921.71899 0.6751213 2000-2009 –> 2010-2019 1.424369 15715.6521 15789.6444
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 10759.89314 10799.38745 3.1875443 2000-2009 –> 2010-2019 6.502066 18337.7513 18414.1766
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + phosphate_star 9799.01743 9837.88814 0.6692637 2000-2009 –> 2010-2019 1.418529 15631.9400 15705.9324
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate 9830.21228 9869.08300 0.6714375 2000-2009 –> 2010-2019 1.406546 15564.8990 15638.8914
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 10758.80137 10798.29568 3.1867094 2000-2009 –> 2010-2019 6.503500 18338.6372 18415.0626
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 9188.27071 9233.61987 0.6279642 2000-2009 –> 2010-2019 1.351127 14840.5797 14926.9041
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + phosphate + phosphate_star 10762.45594 10796.30821 3.1910357 2000-2009 –> 2010-2019 6.506972 18339.5471 18405.0546
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + silicate + phosphate 10757.92822 10791.78049 3.1875711 2000-2009 –> 2010-2019 6.504673 18336.0351 18401.5425
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 10738.90209 10778.39640 3.1715314 2000-2009 –> 2010-2019 6.485336 18316.1348 18392.5602
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate 9863.05765 9901.92837 0.6737339 2000-2009 –> 2010-2019 1.423613 15700.2013 15774.1937
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate 9800.76234 9846.11151 0.6692460 2000-2009 –> 2010-2019 1.418272 15634.0394 15720.3638
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate 10759.06438 10798.55870 3.1869105 2000-2009 –> 2010-2019 6.505385 18340.3665 18416.7918
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 9795.61966 9840.96883 0.6688883 2000-2009 –> 2010-2019 1.418128 15630.3656 15716.6900
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + silicate 9751.55292 9796.90208 0.6658313 2000-2009 –> 2010-2019 1.398602 15471.8370 15558.1614
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + silicate + phosphate 10730.30538 10769.79969 3.1649967 2000-2009 –> 2010-2019 6.482083 18310.3990 18386.8243
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 10778.31299 10817.80730 3.2016623 2000-2009 –> 2010-2019 6.476125 18321.0294 18397.4547
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 9841.51454 9886.86371 0.6720871 2000-2009 –> 2010-2019 1.424225 15696.1441 15782.4685
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 10803.87011 10843.36442 3.2213544 2000-2009 –> 2010-2019 6.503131 18353.0353 18429.4606
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 10749.63881 10789.13312 3.1797118 2000-2009 –> 2010-2019 6.461510 18298.8226 18375.2480
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 22523.38582 22566.75315 5.4014155 2000-2009 –> 2010-2019 10.940501 38034.6869 38118.7522
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 29320.63877 29365.79315 5.5480762 2000-2009 –> 2010-2019 11.846401 50310.9931 50398.6831
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + phosphate + phosphate_star 22869.06300 22906.23501 5.6668293 2000-2009 –> 2010-2019 11.561765 38686.5718 38758.6277
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 22659.16412 22702.53145 5.5035551 2000-2009 –> 2010-2019 11.157808 38272.3294 38356.3947
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 29666.51640 29711.67078 5.7570189 2000-2009 –> 2010-2019 12.279128 50881.5775 50969.2675
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate 30097.80033 30136.50408 6.0299010 2000-2009 –> 2010-2019 12.513246 51272.4949 51347.6578
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 22622.62847 22665.99580 5.4758826 2000-2009 –> 2010-2019 11.292259 38375.7278 38459.7931
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 29778.78953 29823.94391 5.8265200 2000-2009 –> 2010-2019 12.147181 50791.9274 50879.6174
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 22725.74668 22769.11402 5.5543455 2000-2009 –> 2010-2019 11.473504 38565.5545 38649.6198
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 29831.48249 29876.63687 5.8594275 2000-2009 –> 2010-2019 12.190748 50855.4661 50943.1561
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 23228.42173 23271.78907 5.9532317 2000-2009 –> 2010-2019 12.144011 39290.3191 39374.3844
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 30023.26168 30068.41606 5.9807734 2000-2009 –> 2010-2019 12.460875 51196.7360 51284.4260
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 23114.00693 23157.37426 5.8599936 2000-2009 –> 2010-2019 11.956845 39100.2264 39184.2917
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 30125.05430 30170.20868 6.0461991 2000-2009 –> 2010-2019 12.662804 51432.6951 51520.3851
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 23315.89970 23359.26704 6.0255182 2000-2009 –> 2010-2019 12.298940 39443.4393 39527.5047
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 30235.34126 30280.49564 6.1178924 2000-2009 –> 2010-2019 12.770982 51578.3729 51666.0629
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 22827.33802 22870.70535 5.6327460 2000-2009 –> 2010-2019 11.643902 38743.4885 38827.5538
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 29959.55048 30004.70486 5.9401847 2000-2009 –> 2010-2019 12.315784 51028.3883 51116.0783
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 23383.12989 23426.49723 6.0816690 2000-2009 –> 2010-2019 12.391501 39539.3153 39623.3806
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 30283.15067 30328.30505 6.1492350 2000-2009 –> 2010-2019 12.817676 51641.0164 51728.7063
mod Indo-Pacific (26,26.5] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 16317.75800 16359.50466 4.1226561 2000-2009 –> 2010-2019 8.571827 30546.4636 30628.8114
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 16574.92971 16616.67637 4.3112296 2000-2009 –> 2010-2019 8.580516 30602.1502 30684.4980
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 15846.46193 15888.20859 3.7982222 2000-2009 –> 2010-2019 7.865219 29636.7009 29719.0487
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 16273.10124 16314.84789 4.0907619 2000-2009 –> 2010-2019 8.370704 30312.4917 30394.8395
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 17570.49225 17612.43576 4.7100283 2000-2009 –> 2010-2019 9.283729 32257.2843 32339.9849
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 17591.02115 17632.96467 4.7264063 2000-2009 –> 2010-2019 9.317276 32296.5170 32379.2177
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + phosphate 17841.80954 17877.76112 4.9328120 2000-2009 –> 2010-2019 9.796654 32833.6390 32904.5253
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 15477.32849 15519.07514 3.5620492 2000-2009 –> 2010-2019 7.641012 29281.9111 29364.2589
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 17426.14699 17468.09051 4.5964604 2000-2009 –> 2010-2019 9.273682 32224.6674 32307.3680
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 15628.96446 15670.71111 3.6572352 2000-2009 –> 2010-2019 7.905185 29631.7248 29714.0726
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 17497.47131 17539.41482 4.6522305 2000-2009 –> 2010-2019 9.365746 32334.5794 32417.2800
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 17594.22723 17636.17074 4.7289693 2000-2009 –> 2010-2019 9.325275 32305.6306 32388.3312
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate_star 17641.91445 17677.86603 4.7688677 2000-2009 –> 2010-2019 9.536860 32534.3873 32605.2736
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 17409.08752 17451.03104 4.5832206 2000-2009 –> 2010-2019 9.297242 32246.7306 32329.4313
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 15672.16979 15713.91644 3.6848191 2000-2009 –> 2010-2019 7.965062 29711.9035 29794.2513
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 17496.71116 17538.65467 4.6516325 2000-2009 –> 2010-2019 9.366910 32335.6851 32418.3858
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate_star 17496.36433 17532.31591 4.6529330 2000-2009 –> 2010-2019 9.387804 32354.0388 32424.9250
mod Indo-Pacific (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 16081.59670 16123.34336 3.9567624 2000-2009 –> 2010-2019 8.235050 30119.0992 30201.4470
mod Indo-Pacific (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 16149.11357 16190.86023 4.0034968 2000-2009 –> 2010-2019 8.291812 30198.0458 30280.3936
mod Indo-Pacific (26,26.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 16196.45790 16238.20456 4.0365968 2000-2009 –> 2010-2019 8.342468 30265.3357 30347.6835
obs Indo-Pacific (26.5,26.75] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 15864.65691 15899.61557 5.7201887 2000-2009 –> 2010-2019 11.101828 32453.3791 32523.6813
mod Indo-Pacific (26.5,26.75] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 14259.05688 14300.21192 3.5856744 2000-2009 –> 2010-2019 7.922808 27087.5730 27168.6649
mod Indo-Pacific (26.5,26.75] cstar_tref ~ aou + silicate + phosphate + phosphate_star 14301.41521 14336.69096 3.6159022 2000-2009 –> 2010-2019 8.232455 27405.1403 27474.6476
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 15700.67680 15741.46191 5.5338582 2000-2009 –> 2010-2019 10.561678 31927.9721 32009.9913
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 15726.81967 15767.60477 5.5627985 2000-2009 –> 2010-2019 10.648742 32015.5395 32097.5587
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + phosphate 15906.96220 15941.92086 5.7686760 2000-2009 –> 2010-2019 11.144629 32490.0352 32560.3373
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 14077.79066 14118.94570 3.4647547 2000-2009 –> 2010-2019 7.831694 26936.7144 27017.8063
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 15854.53286 15895.31796 5.7063682 2000-2009 –> 2010-2019 11.073275 32430.6053 32512.6244
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 14248.87230 14290.02734 3.5787699 2000-2009 –> 2010-2019 8.059012 27221.5261 27302.6180
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 15729.40628 15770.19138 5.5656701 2000-2009 –> 2010-2019 10.660888 32027.8624 32109.8815
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate_star 15749.31611 15784.27477 5.5900535 2000-2009 –> 2010-2019 10.797773 32162.4831 32232.7852
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 14756.04822 14797.20326 3.9392977 2000-2009 –> 2010-2019 8.552943 27858.9761 27940.0680
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 15735.10623 15775.89134 5.5720033 2000-2009 –> 2010-2019 10.711503 32079.8052 32161.8243
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 14273.98271 14315.13775 3.5958172 2000-2009 –> 2010-2019 8.084606 27255.0978 27336.1897
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 15891.81883 15932.60393 5.7489782 2000-2009 –> 2010-2019 11.146021 32497.8152 32579.8344
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 14433.05028 14474.20532 3.7057102 2000-2009 –> 2010-2019 8.136866 27356.7893 27437.8812
obs Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 15849.32561 15890.11071 5.7004427 2000-2009 –> 2010-2019 11.119876 32477.4462 32559.4654
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 14650.09557 14691.25061 3.8610952 2000-2009 –> 2010-2019 8.388576 27669.3182 27750.4101
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + silicate + phosphate 14478.65195 14513.92770 3.7392444 2000-2009 –> 2010-2019 8.490293 27709.8809 27779.3882
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 14692.33527 14733.49031 3.8920840 2000-2009 –> 2010-2019 8.436326 27727.9649 27809.0567
mod Indo-Pacific (26.75,27] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 18504.22789 18540.93177 3.8169386 2000-2009 –> 2010-2019 7.744326 36865.2625 36938.5618
obs Indo-Pacific (26.75,27] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 24864.78053 24902.57177 5.3320830 2000-2009 –> 2010-2019 10.373383 49455.0408 49530.6664
mod Indo-Pacific (26.75,27] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 17528.27279 17571.09397 3.2988447 2000-2009 –> 2010-2019 7.218997 35879.1671 35964.6830
obs Indo-Pacific (26.75,27] cstar_tref ~ aou + phosphate + phosphate_star 25177.52594 25209.01864 5.5450681 2000-2009 –> 2010-2019 10.764402 50046.6963 50109.7176
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + phosphate 17816.33380 17859.15499 3.4436806 2000-2009 –> 2010-2019 7.144717 35788.5325 35874.0484
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + phosphate 25039.87441 25083.96419 5.4481809 2000-2009 –> 2010-2019 10.586718 49786.7661 49874.9960
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 17996.96373 18039.78492 3.5377271 2000-2009 –> 2010-2019 7.318357 36109.2543 36194.7702
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate 25292.44368 25336.53346 5.6221351 2000-2009 –> 2010-2019 10.786929 50080.5893 50168.8192
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 18019.69779 18062.51898 3.5497443 2000-2009 –> 2010-2019 7.334119 36138.5079 36224.0238
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + phosphate_star 18518.82562 18555.52950 3.8252589 2000-2009 –> 2010-2019 7.632557 36675.3975 36748.6969
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 17964.30660 18007.12779 3.5205357 2000-2009 –> 2010-2019 7.323737 36115.7893 36201.3053
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 25174.72457 25218.81434 5.5403769 2000-2009 –> 2010-2019 10.627707 49840.5527 49928.7826
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + phosphate 25110.65334 25148.44457 5.4977471 2000-2009 –> 2010-2019 10.698409 49952.8147 50028.4403
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 17523.19695 17566.01814 3.2963480 2000-2009 –> 2010-2019 7.225711 35889.5428 35975.0588
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 25082.27319 25126.36297 5.4770021 2000-2009 –> 2010-2019 10.562684 49745.4785 49833.7084
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 17713.70496 17756.52615 3.3913642 2000-2009 –> 2010-2019 7.419909 36244.1754 36329.6913
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 25226.04320 25270.13298 5.5758715 2000-2009 –> 2010-2019 10.740555 50014.0158 50102.2457
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 17742.46669 17785.28788 3.4059452 2000-2009 –> 2010-2019 7.439661 36281.3846 36366.9005
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 25226.54356 25270.63334 5.5762187 2000-2009 –> 2010-2019 10.739549 50012.3946 50100.6245
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate_star 25224.89812 25262.68936 5.5764647 2000-2009 –> 2010-2019 10.741239 50011.0135 50086.6391
mod Indo-Pacific (27,27.25] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 21273.64533 21312.02226 2.6665772 2000-2009 –> 2010-2019 5.848451 43221.1905 43297.7054
mod Indo-Pacific (27,27.25] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 21221.59337 21266.36646 2.6503588 2000-2009 –> 2010-2019 5.661809 42702.4538 42791.7211
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 31802.23763 31848.22746 4.9354858 2000-2009 –> 2010-2019 9.688743 62963.6977 63055.6227
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 32011.38811 32057.37794 5.0343823 2000-2009 –> 2010-2019 9.895740 63408.0694 63499.9944
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate 31803.10961 31842.52946 4.9368306 2000-2009 –> 2010-2019 9.691150 62964.9065 63043.6994
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate_star 32016.03941 32055.45926 5.0375597 2000-2009 –> 2010-2019 9.918110 63451.9351 63530.7280
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 31938.65695 31984.64678 4.9997688 2000-2009 –> 2010-2019 9.819777 63245.9866 63337.9116
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + silicate + phosphate + phosphate_star 31944.95777 31984.37763 5.0037072 2000-2009 –> 2010-2019 9.845575 63297.6196 63376.4124
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate 21171.51227 21216.28536 2.6354199 2000-2009 –> 2010-2019 5.800410 43075.7590 43165.0263
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 21568.94261 21613.71569 2.7563277 2000-2009 –> 2010-2019 6.017774 43728.7875 43818.0548
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate 31802.93395 31848.92378 4.9358118 2000-2009 –> 2010-2019 9.687779 62961.5554 63053.4804
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 31980.39762 32026.38744 5.0196044 2000-2009 –> 2010-2019 9.830556 63268.0560 63359.9810
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 21597.02925 21641.80233 2.7650792 2000-2009 –> 2010-2019 6.034365 43777.3133 43866.5806
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 22076.89464 22121.66773 2.9189681 2000-2009 –> 2010-2019 6.160554 44184.7360 44274.0033
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 32068.93361 32114.92344 5.0619386 2000-2009 –> 2010-2019 9.914705 63447.1127 63539.0377
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 31935.74464 31981.73447 4.9983877 2000-2009 –> 2010-2019 9.798143 63199.0305 63290.9555
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + phosphate 21595.50853 21633.88546 2.7652288 2000-2009 –> 2010-2019 5.993104 43665.2604 43741.7753
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 21595.43711 21640.21019 2.7645824 2000-2009 –> 2010-2019 5.884293 43376.9999 43466.2673
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 21953.11704 21997.89013 2.8784725 2000-2009 –> 2010-2019 6.116830 44052.4749 44141.7423
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 21973.66830 22018.44138 2.8851570 2000-2009 –> 2010-2019 6.127273 44082.9009 44172.1682
obs Indo-Pacific (27.25,27.5] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 33863.08425 33903.16930 4.2851515 2000-2009 –> 2010-2019 8.649674 66611.8045 66691.7356
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 33860.76015 33907.52604 4.2835785 2000-2009 –> 2010-2019 8.591228 66463.0289 66556.2818
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 29439.73017 29486.23942 3.2307795 2000-2009 –> 2010-2019 6.297408 57134.4681 57227.2009
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 33569.90425 33616.67014 4.1790917 2000-2009 –> 2010-2019 8.418170 65990.5570 66083.8099
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + silicate + phosphate 33574.66643 33614.75147 4.1814917 2000-2009 –> 2010-2019 8.431592 66022.7097 66102.6408
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 33818.25693 33865.02282 4.2681483 2000-2009 –> 2010-2019 8.602175 66489.6172 66582.8701
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + silicate + phosphate + phosphate_star 33818.88987 33858.97491 4.2691025 2000-2009 –> 2010-2019 8.646343 66600.5413 66680.4723
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 29025.16348 29071.67273 3.1149423 2000-2009 –> 2010-2019 6.300397 57134.2798 57227.0126
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 29283.36814 29329.87739 3.1865917 2000-2009 –> 2010-2019 6.441119 57626.3094 57719.0423
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 33549.70110 33596.46699 4.1719293 2000-2009 –> 2010-2019 8.366804 65851.7141 65944.9670
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 33763.77279 33810.53868 4.2484497 2000-2009 –> 2010-2019 8.495497 66205.6815 66298.9344
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 29302.52658 29349.03583 3.1919732 2000-2009 –> 2010-2019 6.452547 57665.7019 57758.4348
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 28906.30713 28952.81638 3.0825044 2000-2009 –> 2010-2019 6.133633 56545.8139 56638.5467
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 33704.77396 33751.53985 4.2272214 2000-2009 –> 2010-2019 8.479001 66159.2870 66252.5399
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 33710.49455 33757.26044 4.2292751 2000-2009 –> 2010-2019 8.466977 66127.4711 66220.7240
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + phosphate 29025.51510 29065.38017 3.1155875 2000-2009 –> 2010-2019 6.365777 57351.9208 57431.4061
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 28693.05584 28739.56509 3.0251490 2000-2009 –> 2010-2019 6.157650 56619.4530 56712.1858
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 28970.93610 29017.44535 3.1001006 2000-2009 –> 2010-2019 6.320739 57199.7820 57292.5149
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 28971.85358 29018.36283 3.1003511 2000-2009 –> 2010-2019 6.322104 57204.4732 57297.2061
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + nitrate + silicate + phosphate_star 28969.93550 29009.80057 3.1003735 2000-2009 –> 2010-2019 6.323079 57203.7778 57283.2631
obs Indo-Pacific (27.5,27.75] cstar_tref ~ aou + silicate + phosphate + phosphate_star 38456.27039 38497.26242 4.0036746 2000-2009 –> 2010-2019 7.673456 72679.6179 72761.0911
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 38984.32804 39032.15207 4.1603989 2000-2009 –> 2010-2019 7.801253 73110.1049 73205.1571
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 32979.00396 33026.30845 3.2305417 2000-2009 –> 2010-2019 6.143699 65140.2386 65234.9568
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 38229.49970 38277.32373 3.9373740 2000-2009 –> 2010-2019 7.775868 73020.3813 73115.4335
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 33149.94527 33197.24975 3.2742493 2000-2009 –> 2010-2019 6.187757 65312.7344 65407.4525
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate_star 33417.12790 33457.67460 3.3442778 2000-2009 –> 2010-2019 6.384405 66127.5542 66208.7412
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 33103.58046 33150.88494 3.2623362 2000-2009 –> 2010-2019 6.171512 65247.1459 65341.8641
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 33344.69321 33391.99770 3.3247650 2000-2009 –> 2010-2019 6.328793 65902.7840 65997.5022
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 38018.02195 38065.84598 3.8770621 2000-2009 –> 2010-2019 7.461107 71945.9285 72040.9807
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 33403.22370 33450.52819 3.3400990 2000-2009 –> 2010-2019 6.357944 66020.6034 66115.3216
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 38116.21694 38164.04097 3.9049508 2000-2009 –> 2010-2019 7.521871 72159.0075 72254.0597
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 33398.01390 33445.31839 3.3387312 2000-2009 –> 2010-2019 6.344368 65963.0255 66057.7436
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 33093.38703 33140.69152 3.2597230 2000-2009 –> 2010-2019 6.222772 65474.0194 65568.7376
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 38080.43764 38128.26167 3.8947658 2000-2009 –> 2010-2019 7.506103 72103.7921 72198.8443
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate 33437.66615 33478.21285 3.3496820 2000-2009 –> 2010-2019 6.386115 66132.3833 66213.5703
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 33374.20629 33421.51077 3.3324881 2000-2009 –> 2010-2019 6.342842 65959.4787 66054.1968
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 38101.93561 38149.75963 3.9008823 2000-2009 –> 2010-2019 7.516394 72139.8253 72234.8775
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate_star 38394.39233 38435.38435 3.9856322 2000-2009 –> 2010-2019 7.610787 72463.7959 72545.2692
obs Indo-Pacific (27.5,27.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 38529.98412 38577.80815 4.0246871 2000-2009 –> 2010-2019 7.740428 72911.9278 73006.9800
obs Indo-Pacific (27.5,27.75] cstar_tref ~ temp + aou + silicate + phosphate 38531.15773 38572.14975 4.0256195 2000-2009 –> 2010-2019 7.748645 72935.7437 73017.2169
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + phosphate 11581.76073 11622.48692 2.4807941 2000-2009 –> 2010-2019 4.713910 22707.4419 22788.9364
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 11621.48901 11662.21521 2.5007041 2000-2009 –> 2010-2019 4.741745 22764.8847 22846.3792
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + phosphate 11905.89182 11940.79999 2.6490436 2000-2009 –> 2010-2019 4.930829 23137.3755 23207.2280
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 11663.86677 11704.59296 2.5221181 2000-2009 –> 2010-2019 4.803496 22896.4578 22977.9524
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 21266.71779 21310.44833 3.9157903 2000-2009 –> 2010-2019 7.006252 37664.6877 37750.9495
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 11703.36522 11744.09142 2.5422421 2000-2009 –> 2010-2019 4.831025 22952.1597 23033.6542
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 21246.17932 21289.90986 3.9052694 2000-2009 –> 2010-2019 6.997234 37647.2768 37733.5386
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 11646.48338 11687.20957 2.5133119 2000-2009 –> 2010-2019 4.759959 22802.3716 22883.8661
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 21244.31841 21288.04895 3.9043175 2000-2009 –> 2010-2019 7.005046 37663.6214 37749.8831
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 11720.09273 11760.81892 2.5508130 2000-2009 –> 2010-2019 4.841824 22973.7521 23055.2466
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 21245.62876 21289.35930 3.9049878 2000-2009 –> 2010-2019 6.997357 37647.5670 37733.8288
mod Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 11572.94476 11613.67095 2.4763975 2000-2009 –> 2010-2019 4.873585 23053.1182 23134.6127
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 21270.74315 21314.47369 3.9178556 2000-2009 –> 2010-2019 7.008125 37668.3126 37754.5744
mod Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 11622.54730 11663.27350 2.5012367 2000-2009 –> 2010-2019 4.907111 23120.8068 23202.3013
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 21252.68647 21296.41701 3.9085996 2000-2009 –> 2010-2019 7.000556 37653.7672 37740.0290
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + silicate + phosphate_star 21256.94411 21294.42743 3.9118049 2000-2009 –> 2010-2019 7.003788 37656.0792 37730.0178
mod Indo-Pacific (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 11656.45351 11697.17970 2.5183589 2000-2009 –> 2010-2019 4.930173 23167.0427 23248.5372
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 21252.83211 21296.56265 3.9086742 2000-2009 –> 2010-2019 7.001019 37654.7211 37740.9829
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate_star 21250.91097 21288.39429 3.9087146 2000-2009 –> 2010-2019 7.009447 37668.2204 37742.1591
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + silicate + phosphate + phosphate_star 21256.02291 21293.50623 3.9113328 2000-2009 –> 2010-2019 7.003716 37655.9907 37729.9294
obs Indo-Pacific (27.85,27.95] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 27307.00515 27352.92290 3.3096330 2000-2009 –> 2010-2019 6.053580 48709.7894 48800.4391
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 13742.55766 13784.99469 2.1051817 2000-2009 –> 2010-2019 4.122640 26434.5621 26519.0180
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 27332.88104 27378.79878 3.3178509 2000-2009 –> 2010-2019 6.062919 48739.2621 48829.9117
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate 13179.82546 13216.20006 1.9271493 2000-2009 –> 2010-2019 3.879216 25672.8565 25745.2473
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + phosphate 13181.80201 13224.23904 1.9271421 2000-2009 –> 2010-2019 3.874323 25661.8516 25746.3075
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 13180.95484 13223.39187 1.9268849 2000-2009 –> 2010-2019 3.875440 25665.2233 25749.6791
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + silicate 13166.24848 13208.68551 1.9224247 2000-2009 –> 2010-2019 3.842934 25563.8481 25648.3040
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + phosphate 13794.45230 13830.82690 2.1231365 2000-2009 –> 2010-2019 4.176452 26589.7742 26662.1649
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 13654.68762 13697.12465 2.0762331 2000-2009 –> 2010-2019 4.123338 26433.9003 26518.3562
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 27289.53495 27335.45270 3.3040962 2000-2009 –> 2010-2019 6.052574 48706.8484 48797.4980
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 13670.97518 13713.41221 2.0815688 2000-2009 –> 2010-2019 4.134204 26466.3165 26550.7723
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 27282.08115 27327.99890 3.3017366 2000-2009 –> 2010-2019 6.061316 48734.8999 48825.5495
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 13223.18371 13265.62074 1.9397499 2000-2009 –> 2010-2019 3.897393 25735.2673 25819.7232
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 13686.28544 13728.72247 2.0865968 2000-2009 –> 2010-2019 4.141448 26488.0751 26572.5310
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 27282.09480 27328.01254 3.3017409 2000-2009 –> 2010-2019 6.062993 48740.2527 48830.9024
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 27272.55394 27318.47168 3.2987232 2000-2009 –> 2010-2019 6.046878 48688.8322 48779.4818
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 27270.72144 27316.63918 3.2981439 2000-2009 –> 2010-2019 6.057849 48723.9441 48814.5937
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + silicate + phosphate 27293.42647 27332.78453 3.3059623 2000-2009 –> 2010-2019 6.055109 48710.8833 48788.5830
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + silicate + phosphate_star 27287.84973 27327.20780 3.3041958 2000-2009 –> 2010-2019 6.064352 48740.5095 48818.2091
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 27270.74882 27316.66657 3.2981526 2000-2009 –> 2010-2019 6.059129 48728.0278 48818.6774
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 46570.69817 46613.17272 3.4400680 2000-2009 –> 2010-2019 6.291374 83277.5584 83361.5191
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 46570.78261 46620.33625 3.4396923 2000-2009 –> 2010-2019 6.278605 83214.8434 83312.7975
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + phosphate + phosphate_star 46726.36028 46761.75574 3.4711292 2000-2009 –> 2010-2019 6.326361 83451.6924 83521.6597
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + silicate + phosphate + phosphate_star 46725.45436 46767.92891 3.4705542 2000-2009 –> 2010-2019 6.318537 83414.9704 83498.9311
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 46572.67821 46622.23185 3.4400641 2000-2009 –> 2010-2019 6.291342 83281.3947 83379.3488
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + aou + phosphate + phosphate_star 46713.97446 46756.44901 3.4682834 2000-2009 –> 2010-2019 6.323362 83440.5069 83524.4676
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate 17728.97137 17767.43169 1.7386687 2000-2009 –> 2010-2019 3.464454 30518.8648 30593.8473
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate + silicate 17593.50339 17638.37376 1.7122671 2000-2009 –> 2010-2019 3.391211 30206.4229 30293.9024
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate 18027.84769 18066.30801 1.7974830 2000-2009 –> 2010-2019 3.624272 31187.6731 31262.6555
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate_star 18030.31513 18068.77545 1.7979767 2000-2009 –> 2010-2019 3.628918 31204.9057 31279.8882
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 18021.26961 18066.13998 1.7957675 2000-2009 –> 2010-2019 3.592608 31075.5852 31163.0647
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 46812.90932 46862.46296 3.4875040 2000-2009 –> 2010-2019 6.345293 83555.5574 83653.5116
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 18024.80759 18069.67796 1.7964749 2000-2009 –> 2010-2019 3.599147 31100.1956 31187.6752
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 18237.58399 18282.45437 1.8395303 2000-2009 –> 2010-2019 3.680984 31451.4124 31538.8919
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + phosphate + phosphate_star 18036.65695 18075.11727 1.7992464 2000-2009 –> 2010-2019 3.632205 31218.4129 31293.3954
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 18030.90817 18075.77854 1.7976952 2000-2009 –> 2010-2019 3.601900 31111.8237 31199.3033
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate_star 18286.81110 18325.27142 1.8500493 2000-2009 –> 2010-2019 3.708614 31558.7962 31633.7787
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 46661.30567 46710.85932 3.4574903 2000-2009 –> 2010-2019 6.309571 83374.2072 83472.1613
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 46776.21987 46825.77351 3.4802166 2000-2009 –> 2010-2019 6.360348 83634.7187 83732.6728
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 46760.73947 46810.29311 3.4771464 2000-2009 –> 2010-2019 6.355424 83609.6566 83707.6108
obs Indo-Pacific (28.05,28.1] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 35866.85515 35908.11022 2.9622279 2000-2009 –> 2010-2019 5.690531 58728.9900 58808.9974
obs Indo-Pacific (28.05,28.1] cstar_tref ~ aou + phosphate + phosphate_star 35940.31596 35974.69519 2.9778855 2000-2009 –> 2010-2019 5.707847 58806.1848 58872.8577
obs Indo-Pacific (28.05,28.1] cstar_tref ~ aou + silicate + phosphate + phosphate_star 35915.54551 35956.80059 2.9723214 2000-2009 –> 2010-2019 5.701305 58780.0347 58860.0421
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 35883.85327 35931.98419 2.9653334 2000-2009 –> 2010-2019 5.705886 58790.2427 58883.5846
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + phosphate 35879.16015 35927.29108 2.9643613 2000-2009 –> 2010-2019 5.719002 58833.9127 58927.2546
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + silicate 9299.40343 9340.40580 1.4580175 2000-2009 –> 2010-2019 3.129560 16746.5115 16826.4453
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 9480.25842 9521.26078 1.5099239 2000-2009 –> 2010-2019 3.088365 16706.9556 16786.8895
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 35878.06056 35926.19148 2.9641336 2000-2009 –> 2010-2019 5.676886 58688.2835 58781.6254
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 9485.43633 9526.43869 1.5114369 2000-2009 –> 2010-2019 3.094299 16722.8903 16802.8242
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 35997.76321 36045.89413 2.9890254 2000-2009 –> 2010-2019 5.724590 58886.9706 58980.3125
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 9437.56812 9478.57048 1.4975073 2000-2009 –> 2010-2019 3.080348 16674.9710 16754.9049
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 35967.57235 36015.70327 2.9827277 2000-2009 –> 2010-2019 5.712434 58836.5579 58929.8999
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 9484.84505 9525.84742 1.5112640 2000-2009 –> 2010-2019 3.094126 16722.2987 16802.2325
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 35999.16760 36047.29852 2.9893187 2000-2009 –> 2010-2019 5.724521 58887.1237 58980.4656
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate_star 9490.08008 9525.22496 1.5133804 2000-2009 –> 2010-2019 3.096244 16725.5376 16794.0524
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate_star 35997.26847 36038.52355 2.9893398 2000-2009 –> 2010-2019 5.724973 58884.7123 58964.7197
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate 9892.15976 9927.30465 1.6357763 2000-2009 –> 2010-2019 3.238987 17176.7424 17245.2571
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate_star 9900.57352 9935.71841 1.6384406 2000-2009 –> 2010-2019 3.247574 17199.3372 17267.8519
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + nitrate + silicate + phosphate_star 9891.85601 9927.00089 1.6356802 2000-2009 –> 2010-2019 3.244350 17189.5108 17258.0255
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + silicate + phosphate + phosphate_star 9898.56498 9933.70987 1.6378042 2000-2009 –> 2010-2019 3.246914 17197.2727 17265.7874
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 18999.44436 19045.81161 1.3330438 2000-2009 –> 2010-2019 2.774782 36475.8100 36567.6448
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 37470.80466 37519.46746 2.7357127 2000-2009 –> 2010-2019 5.269858 66999.7607 67095.6428
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + phosphate + phosphate_star 37532.94241 37574.65338 2.7470975 2000-2009 –> 2010-2019 5.286373 67085.3108 67167.4954
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + silicate + phosphate 37504.70987 37546.42084 2.7420803 2000-2009 –> 2010-2019 5.278084 67040.8797 67123.0643
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + nitrate + silicate + phosphate 19087.00622 19126.74958 1.3438178 2000-2009 –> 2010-2019 2.806727 36703.9978 36782.7133
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 19088.87805 19135.24530 1.3438023 2000-2009 –> 2010-2019 2.800620 36667.0440 36758.8788
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 37586.35302 37635.01582 2.7562574 2000-2009 –> 2010-2019 5.325081 67286.1048 67381.9868
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + nitrate + phosphate 37625.57214 37674.23494 2.7632657 2000-2009 –> 2010-2019 5.337734 67352.9038 67448.7858
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + nitrate + silicate 18148.26731 18194.63456 1.2348648 2000-2009 –> 2010-2019 2.757287 36157.3998 36249.2345
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate 17727.97981 17774.34706 1.1890875 2000-2009 –> 2010-2019 2.564223 34741.5908 34833.4255
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate 37488.46808 37537.13088 2.7388434 2000-2009 –> 2010-2019 5.256135 66933.5760 67029.4580
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 17746.51702 17792.88427 1.1910703 2000-2009 –> 2010-2019 2.568102 34773.6065 34865.4412
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 37589.88235 37638.54515 2.7568874 2000-2009 –> 2010-2019 5.294175 67134.4137 67230.2958
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 18561.56389 18607.93114 1.2815988 2000-2009 –> 2010-2019 2.720577 36019.1835 36111.0183
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 17708.49792 17754.86517 1.1870072 2000-2009 –> 2010-2019 2.574357 34808.6294 34900.4642
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 37666.69049 37715.35329 2.7706324 2000-2009 –> 2010-2019 5.325198 67296.5017 67392.3837
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 17742.93637 17789.30362 1.1906870 2000-2009 –> 2010-2019 2.567451 34768.1224 34859.9571
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 37573.57350 37622.23631 2.7539776 2000-2009 –> 2010-2019 5.287237 67098.1350 67194.0170
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate_star 17930.45480 17970.19815 1.2111427 2000-2009 –> 2010-2019 2.601611 35050.5536 35129.2691
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate_star 37697.87798 37739.58895 2.7765926 2000-2009 –> 2010-2019 5.339151 67364.9448 67447.1293
obs Indo-Pacific (28.2,28.3] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 11277.68197 11312.98717 2.0191140 2000-2009 –> 2010-2019 4.629177 19528.1577 19596.2051
obs Indo-Pacific (28.2,28.3] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 11199.59632 11240.78572 1.9888902 2000-2009 –> 2010-2019 4.596971 19449.4402 19528.8288
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 11219.15842 11260.34782 1.9962308 2000-2009 –> 2010-2019 4.603793 19468.3142 19547.7028
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 31806.76549 31860.09198 0.6965774 2000-2009 –> 2010-2019 1.350853 54022.7289 54127.2960
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 34910.84312 34964.16961 0.7723317 2000-2009 –> 2010-2019 1.541643 60741.9199 60846.4871
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + phosphate + phosphate_star 11349.62648 11384.93168 2.0466569 2000-2009 –> 2010-2019 4.655855 19598.9546 19667.0020
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 11340.91667 11382.10607 2.0425332 2000-2009 –> 2010-2019 4.651596 19592.0653 19671.4539
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + nitrate + silicate + phosphate 32248.86745 32294.57587 0.7069421 2000-2009 –> 2010-2019 1.453823 57417.4853 57507.1143
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 32076.38310 32129.70958 0.7028516 2000-2009 –> 2010-2019 1.362764 54483.8325 54588.3997
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + nitrate + phosphate 11203.82281 11245.01221 1.9904739 2000-2009 –> 2010-2019 4.648647 19519.5675 19598.9561
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 11334.91796 11376.10736 2.0402270 2000-2009 –> 2010-2019 4.648460 19584.9642 19664.3529
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + silicate + phosphate 34862.11215 34915.43864 0.7710810 2000-2009 –> 2010-2019 1.601983 62412.1969 62516.7641
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 34874.25314 34927.57962 0.7713924 2000-2009 –> 2010-2019 1.602597 62432.4637 62537.0308
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 32062.54602 32115.87251 0.7025283 2000-2009 –> 2010-2019 1.441566 56997.5596 57102.1267
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 34871.86790 34925.19439 0.7713313 2000-2009 –> 2010-2019 1.602498 62429.0559 62533.6230
obs Indo-Pacific (28.2,28.3] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 11237.01149 11278.20089 2.0029537 2000-2009 –> 2010-2019 4.612941 19489.3872 19568.7759
mod Indo-Pacific (28.2,28.3] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 36079.01827 36132.34476 0.8029282 2000-2009 –> 2010-2019 1.640421 63805.4462 63910.0133
obs Indo-Pacific (28.2,28.3] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 11207.86922 11249.05862 1.9919912 2000-2009 –> 2010-2019 4.639956 19510.2849 19589.6736
obs Indo-Pacific (28.2,28.3] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 11349.60783 11390.79723 2.0458790 2000-2009 –> 2010-2019 4.655429 19601.4023 19680.7909
mod Indo-Pacific (28.2,28.3] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 36077.52613 36130.85262 0.8028883 2000-2009 –> 2010-2019 1.641631 63837.2431 63941.8102
mod Indo-Pacific (28.3,28.4] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 20703.92726 20744.38749 1.2603615 2000-2009 –> 2010-2019 2.189667 33631.1980 33710.5152
mod Indo-Pacific (28.3,28.4] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 20705.11409 20752.31769 1.2602798 2000-2009 –> 2010-2019 2.187224 33609.9637 33702.5004
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 20391.47441 20438.67802 1.2291448 2000-2009 –> 2010-2019 2.157909 33315.1454 33407.6822
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 4493.00233 4527.69148 2.0461356 2000-2009 –> 2010-2019 4.694379 7770.8292 7837.1935
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 4476.55291 4511.24206 2.0301556 2000-2009 –> 2010-2019 4.664660 7747.2852 7813.6496
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + silicate + phosphate_star 4496.83808 4526.57163 2.0518350 2000-2009 –> 2010-2019 4.702393 7773.8567 7830.7405
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 4524.62273 4559.31188 2.0772079 2000-2009 –> 2010-2019 4.710645 7794.8025 7861.1669
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + silicate + phosphate + phosphate_star 4525.77911 4555.51266 2.0803353 2000-2009 –> 2010-2019 4.715245 7794.7210 7851.6047
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + nitrate + phosphate 20368.38467 20415.58827 1.2268833 2000-2009 –> 2010-2019 2.152935 33263.9844 33356.5212
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 20450.84761 20498.05121 1.2349792 2000-2009 –> 2010-2019 2.165995 33397.7599 33490.2967
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + silicate + phosphate 4502.82522 4537.51437 2.0557381 2000-2009 –> 2010-2019 4.689266 7773.0517 7839.4161
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4476.36674 4511.05588 2.0299754 2000-2009 –> 2010-2019 4.679688 7754.9500 7821.3144
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 20366.46049 20413.66409 1.2266950 2000-2009 –> 2010-2019 2.153816 33273.1356 33365.6724
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 4502.83383 4537.52298 2.0557466 2000-2009 –> 2010-2019 4.705747 7781.5658 7847.9301
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4496.38438 4531.07352 2.0494367 2000-2009 –> 2010-2019 4.683800 7767.0433 7833.4077
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + aou + nitrate + phosphate 20684.67634 20725.13657 1.2584279 2000-2009 –> 2010-2019 2.185020 33583.8811 33663.1983
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 20686.34932 20733.55292 1.2583950 2000-2009 –> 2010-2019 2.182298 33559.6587 33652.1955
obs Indo-Pacific (28.3,28.4] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 4473.72820 4508.41735 2.0274240 2000-2009 –> 2010-2019 4.712820 7770.5577 7836.9221
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 20687.51931 20727.97954 1.2587132 2000-2009 –> 2010-2019 2.186432 33598.3850 33677.7022
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 20689.15096 20736.35457 1.2586763 2000-2009 –> 2010-2019 2.183765 33574.7712 33667.3079
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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