Last updated: 2021-06-10

Checks: 7 0

Knit directory: emlr_obs_v_XXX/

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

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

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

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


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

Untracked files:
    Untracked:  .Rprofile
    Untracked:  .gitattributes
    Untracked:  .gitignore
    Untracked:  README.md
    Untracked:  _workflowr.yml
    Untracked:  analysis/_site.yml
    Untracked:  code/
    Untracked:  data/
    Untracked:  emlr_obs_v_XXX.Rproj

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

1 Version ID

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

2 Required data

Required are:

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

3 Predictor combinations

Find all possible combinations of following considered predictor variables:

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

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

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

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

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

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

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

# convert to tibble
lm_all <- as_tibble(lm_all)

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

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

4 Apply predictor threshold

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

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

This results in a total number of MLR models of:

  • 112

5 Fit all models

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

  • cstar_tref

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

library(furrr)
library(purrr)

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



library(tictoc)
plan(sequential)

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

plan(multisession, workers = 3)

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


future_map_dbl(1:4, function(x){
  Sys.sleep(1)
  x^2
}, .progress = TRUE)
# prepare nested data frame

GLODAP %>%
  filter_all(any_vars(is.na(.)))

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

# expand with model definitions

GLODAP_nested_lm <- expand_grid(
  GLODAP_nested,
  lm_all#[1:3,]
)

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


# test parallel processing

# expand with model definitions


# test multidplyr
 
# library(parallel)
# library(multidplyr)
# 
# n_cores <- detectCores()
# system.time(
# cluster <- new_cluster(n_cores - 2)
# )
# 
# cluster
# 
# system.time(
#   GLODAP_nested_lm_part <-
#     GLODAP_nested_lm %>%
#     #group_by(gamma_slab, era, basin, data_source, model, predictors, n) %>% 
#     partition(cluster)
# )
# 
# 
# GLODAP_nested_lm
# GLODAP_nested_lm_part
# 
# system.time(cluster %>%
#               cluster_library(c("tidyverse",
#                                 "broom",
#                                 "olsrr",
#                                 "purrr")))
# 
# cluster_call(cluster, search())
# 
# system.time(
# GLODAP_nested_lm_fit <- GLODAP_nested_lm_part %>%
#   mutate(
#     fit = map2(.x = data, .y = model,
#                ~ lm(as.formula(.y), data = .x)),
#     tidied = map(fit, tidy),
#     glanced = map(fit, glance),
#     augmented = map(fit, augment),
#     vif = map(fit, ols_vif_tol)
#   )
# )
# 
# system.time(
#   # extract glanced model output (model diagnostics, such as AIC)
#   GLODAP_glanced <- GLODAP_nested_lm_fit %>%
#     select(-c(data, fit, tidied, augmented, vif))
#     # unnest(glanced) %>%
#     # rename(n_predictors = n)
# )
# 
# 
# system.time(
# GLODAP_nested_lm_fit <- GLODAP_nested_lm_fit %>%
#   collect()
# )


### test furrr

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

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

6 Tidy models

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

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

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

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

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

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

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

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


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

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

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

7 Prepare coeffcients

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

7.1 VIF threshold

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

  • 1000

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

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

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

7.2 Predictor selection

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

  • 10

The criterion used to select the best models was:

  • rmse

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

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

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

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

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

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


# print table
lm_best_eras %>% 
  kable() %>%
  add_header_above() %>%
  kable_styling() %>%
  scroll_box(width = "100%", height = "400px")
data_source basin gamma_slab model AIC BIC rmse eras rmse_sum aic_sum bic_sum
mod Atlantic (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 111.45352 124.97630 0.6290706 2000-2009 –> 2010-2019 4.679463 767.9011 800.5772
obs Atlantic (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 850.04536 869.89934 6.6772204 2000-2009 –> 2010-2019 14.634952 1276.2291 1310.6258
mod Atlantic (-Inf,26] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 61.25482 74.77760 0.3845598 2000-2009 –> 2010-2019 4.768797 735.7604 768.4365
mod Atlantic (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 83.19634 96.71912 0.4768549 2000-2009 –> 2010-2019 4.462004 735.9414 768.6176
obs Atlantic (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 842.62498 862.47896 6.4834699 2000-2009 –> 2010-2019 15.395399 1282.1718 1316.5686
mod Atlantic (-Inf,26] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 59.55653 73.07931 0.3782099 2000-2009 –> 2010-2019 4.250925 705.7764 738.4526
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate 852.26946 869.28715 6.7900892 2000-2009 –> 2010-2019 14.843410 1277.8622 1307.3451
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 850.65407 870.50805 6.6933688 2000-2009 –> 2010-2019 14.682740 1277.3060 1311.7028
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 850.70152 870.55549 6.6946292 2000-2009 –> 2010-2019 14.680822 1277.3065 1311.7033
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 56.00703 69.52981 0.3652750 2000-2009 –> 2010-2019 4.657755 725.6902 758.3663
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 844.99684 864.85082 6.5447813 2000-2009 –> 2010-2019 14.530820 1271.5996 1305.9963
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 81.48201 95.00479 0.4689073 2000-2009 –> 2010-2019 4.346113 727.9661 760.6423
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 81.96023 95.48301 0.4711109 2000-2009 –> 2010-2019 4.757791 751.3351 784.0113
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 110.87147 124.39425 0.6254911 2000-2009 –> 2010-2019 4.661631 766.5153 799.1915
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 850.91560 870.76957 6.7003188 2000-2009 –> 2010-2019 14.710092 1277.8685 1312.2652
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 55.79514 69.31792 0.3645169 2000-2009 –> 2010-2019 4.672642 726.3077 758.9839
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 81.29552 94.81830 0.4680508 2000-2009 –> 2010-2019 4.655443 745.3274 778.0035
obs Atlantic (-Inf,26] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 857.70903 877.56300 6.8834031 2000-2009 –> 2010-2019 15.374758 1291.5515 1325.9482
obs Atlantic (-Inf,26] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 856.78870 876.64267 6.8583100 2000-2009 –> 2010-2019 15.341501 1290.5177 1324.9144
obs Atlantic (-Inf,26] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 853.96517 873.81914 6.7818951 2000-2009 –> 2010-2019 15.228976 1287.1908 1321.5875
obs Atlantic (26,26.5] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 3999.53755 4026.08676 6.1254567 2000-2009 –> 2010-2019 11.594819 7254.3622 7306.4343
mod Atlantic (26,26.5] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 2097.95276 2125.44124 3.8948889 2000-2009 –> 2010-2019 8.365265 3459.3346 3510.9200
obs Atlantic (26,26.5] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 4001.53551 4032.50960 6.1254466 2000-2009 –> 2010-2019 11.572078 7254.0289 7314.7798
mod Atlantic (26,26.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 2052.16339 2079.65188 3.6642096 2000-2009 –> 2010-2019 7.670586 3362.9168 3414.5023
obs Atlantic (26,26.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 3867.54184 3898.51592 5.4951553 2000-2009 –> 2010-2019 10.952916 7122.1581 7182.9089
mod Atlantic (26,26.5] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2024.04141 2051.52989 3.5293605 2000-2009 –> 2010-2019 7.920533 3377.1645 3428.7499
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 2055.87068 2083.35916 3.6823668 2000-2009 –> 2010-2019 7.969995 3397.9693 3449.5547
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 3884.79124 3915.76533 5.5725084 2000-2009 –> 2010-2019 11.099189 7152.4581 7213.2090
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 2056.14520 2083.63369 3.6837149 2000-2009 –> 2010-2019 8.041122 3405.7020 3457.2874
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 3885.73293 3916.70702 5.5767625 2000-2009 –> 2010-2019 11.177288 7167.2038 7227.9547
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 2028.14876 2055.63725 3.5487419 2000-2009 –> 2010-2019 8.287166 3416.4339 3468.0193
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2028.83814 2056.32663 3.5520053 2000-2009 –> 2010-2019 8.407711 3428.4191 3480.0045
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 2057.15020 2084.63868 3.6886544 2000-2009 –> 2010-2019 8.118631 3414.3380 3465.9234
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 3895.97570 3926.94978 5.6232449 2000-2009 –> 2010-2019 11.284757 7188.7104 7249.4613
obs Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + phosphate 3998.07423 4024.62345 6.1181973 2000-2009 –> 2010-2019 11.679362 7270.2102 7322.2824
mod Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 2095.49168 2122.98016 3.8821289 2000-2009 –> 2010-2019 8.392623 3461.0011 3512.5865
obs Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 4000.07421 4031.04829 6.1181972 2000-2009 –> 2010-2019 11.600770 7259.4078 7320.1587
mod Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2091.81192 2119.30040 3.8631285 2000-2009 –> 2010-2019 8.410037 3461.0363 3512.6217
obs Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 3999.99486 4030.96894 6.1178038 2000-2009 –> 2010-2019 11.643212 7267.4223 7328.1732
obs Atlantic (26,26.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 4000.51805 4031.49213 6.1203981 2000-2009 –> 2010-2019 11.700824 7278.2497 7339.0006
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 2753.16219 2782.42230 4.1232145 2000-2009 –> 2010-2019 7.291389 4655.0792 4711.6768
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 4797.43995 4830.03703 5.2343971 2000-2009 –> 2010-2019 10.905744 8066.7447 8129.0646
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 2756.22777 2785.48789 4.1363203 2000-2009 –> 2010-2019 7.195029 4632.3353 4688.9330
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 4796.55634 4829.15343 5.2314255 2000-2009 –> 2010-2019 10.931951 8071.1572 8133.4770
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 2802.74633 2832.00645 4.3403823 2000-2009 –> 2010-2019 7.193008 4627.6553 4684.2529
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 4891.09595 4923.69304 5.5591315 2000-2009 –> 2010-2019 11.257001 8165.2158 8227.5356
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2808.93379 2838.19391 4.3682727 2000-2009 –> 2010-2019 7.166906 4619.8169 4676.4145
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 2762.82496 2792.08508 4.1646655 2000-2009 –> 2010-2019 7.089965 4606.1993 4662.7969
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 4793.66440 4826.26148 5.2217115 2000-2009 –> 2010-2019 10.920882 8068.0199 8130.3398
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 2809.91398 2839.17409 4.3727074 2000-2009 –> 2010-2019 7.167402 4619.7635 4676.3612
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + phosphate 4858.56585 4886.50621 5.4511195 2000-2009 –> 2010-2019 11.126621 8126.6264 8180.0434
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + phosphate_star 4855.62345 4883.56381 5.4408212 2000-2009 –> 2010-2019 11.141728 8128.2933 8181.7103
mod Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 2775.11060 2804.37072 4.2179702 2000-2009 –> 2010-2019 7.286107 4653.4771 4710.0748
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 4858.74349 4891.34057 5.4447390 2000-2009 –> 2010-2019 11.119518 8128.6727 8190.9926
mod Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2778.97661 2808.23673 4.2348847 2000-2009 –> 2010-2019 7.202925 4632.9974 4689.5950
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 4855.95267 4888.54976 5.4349822 2000-2009 –> 2010-2019 11.128072 8129.2064 8191.5263
mod Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 2788.23076 2813.31086 4.2845107 2000-2009 –> 2010-2019 7.349424 4663.8258 4712.3381
obs Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 4853.19552 4881.13588 5.4323382 2000-2009 –> 2010-2019 11.132562 8125.7417 8179.1587
mod Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 2782.74202 2812.00214 4.2514242 2000-2009 –> 2010-2019 7.171361 4624.7694 4681.3671
obs Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 4853.49719 4886.09427 5.4264121 2000-2009 –> 2010-2019 11.119414 8126.7351 8189.0549
mod Atlantic (26.75,27] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 4280.62621 4314.35102 2.4970301 2000-2009 –> 2010-2019 4.824745 7128.9361 7193.7364
mod Atlantic (26.75,27] cstar_tref ~ aou + silicate + phosphate + phosphate_star 4650.39808 4679.30506 3.0601941 2000-2009 –> 2010-2019 5.422862 7515.3679 7570.9109
mod Atlantic (26.75,27] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 4429.84146 4463.56627 2.7094066 2000-2009 –> 2010-2019 5.238856 7382.2105 7447.0108
mod Atlantic (26.75,27] cstar_tref ~ sal + aou + silicate + phosphate 4578.29914 4607.20612 2.9418450 2000-2009 –> 2010-2019 5.477735 7531.8520 7587.3951
mod Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate 4559.10414 4592.82895 2.9079323 2000-2009 –> 2010-2019 5.183321 7378.9483 7443.7486
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate 5689.49106 5723.66813 4.4440659 2000-2009 –> 2010-2019 10.822515 9714.7957 9779.9013
mod Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4579.17588 4612.90069 2.9400378 2000-2009 –> 2010-2019 5.289495 7439.1255 7503.9258
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5694.34380 5728.52087 4.4551391 2000-2009 –> 2010-2019 10.777458 9708.8119 9773.9175
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 5704.48578 5738.66284 4.4783707 2000-2009 –> 2010-2019 10.786980 9716.2925 9781.3981
mod Atlantic (26.75,27] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4586.45223 4620.17705 2.9517640 2000-2009 –> 2010-2019 5.334606 7464.0674 7528.8676
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5695.39007 5729.56713 4.4575301 2000-2009 –> 2010-2019 10.776722 9709.2517 9774.3573
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + silicate + phosphate_star 5714.04633 5743.34096 4.5049996 2000-2009 –> 2010-2019 10.910772 9742.5915 9798.3963
mod Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 4300.12051 4333.84532 2.5238016 2000-2009 –> 2010-2019 4.799508 7120.1393 7184.9396
obs Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 5695.07538 5729.25244 4.4568109 2000-2009 –> 2010-2019 10.925787 9737.6580 9802.7636
mod Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 4310.10485 4343.82966 2.5376241 2000-2009 –> 2010-2019 4.830269 7139.4080 7204.2082
obs Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5700.85943 5735.03649 4.4700502 2000-2009 –> 2010-2019 10.851550 9726.7503 9791.8559
mod Atlantic (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 4332.56108 4366.28589 2.5689900 2000-2009 –> 2010-2019 4.894840 7179.8675 7244.6678
obs Atlantic (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 5703.01658 5737.19364 4.4749978 2000-2009 –> 2010-2019 10.868051 9731.1250 9796.2306
obs Atlantic (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate_star 5707.90627 5737.20089 4.4908367 2000-2009 –> 2010-2019 10.914751 9739.9187 9795.7235
obs Atlantic (26.75,27] cstar_tref ~ temp + silicate + phosphate + phosphate_star 5703.24269 5732.53732 4.4801094 2000-2009 –> 2010-2019 10.931528 9740.4932 9796.2980
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 5426.02689 5462.01007 2.0653251 2000-2009 –> 2010-2019 4.383761 9412.8106 9482.2453
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 5573.59031 5609.57348 2.1896720 2000-2009 –> 2010-2019 4.668588 9678.0341 9747.4688
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate 5720.56050 5751.40322 2.3227997 2000-2009 –> 2010-2019 4.786251 9812.0027 9871.5181
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 5721.66809 5757.65126 2.3219786 2000-2009 –> 2010-2019 4.666851 9728.3839 9797.8186
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate 11087.11837 11125.22811 6.1645514 2000-2009 –> 2010-2019 12.061949 19059.1951 19133.1987
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 11061.30322 11099.41296 6.1181948 2000-2009 –> 2010-2019 11.901566 18984.7250 19058.7286
mod Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + silicate 5524.25433 5560.23750 2.1472866 2000-2009 –> 2010-2019 4.412252 9470.0177 9539.4523
mod Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate 5720.61201 5756.59518 2.3210072 2000-2009 –> 2010-2019 4.700581 9753.1534 9822.5881
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 11106.06811 11144.17785 6.1988032 2000-2009 –> 2010-2019 11.937889 19010.3344 19084.3380
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 11068.19819 11106.30793 6.1305420 2000-2009 –> 2010-2019 11.946990 19005.8318 19079.8354
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 11124.66378 11162.77352 6.2325999 2000-2009 –> 2010-2019 12.067271 19070.0926 19144.0962
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 11107.73035 11145.84009 6.2018167 2000-2009 –> 2010-2019 11.954192 19017.7603 19091.7639
obs Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + phosphate_star 11081.93318 11114.59867 6.1588128 2000-2009 –> 2010-2019 12.069942 19057.8057 19121.2374
mod Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate 5689.30161 5720.14433 2.2942100 2000-2009 –> 2010-2019 4.597176 9662.3156 9721.8310
mod Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 5690.75009 5726.73326 2.2937088 2000-2009 –> 2010-2019 4.554695 9633.4224 9702.8570
obs Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 11048.02035 11086.13009 6.0944785 2000-2009 –> 2010-2019 12.069399 19052.6418 19126.6454
mod Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5691.07911 5727.06228 2.2940078 2000-2009 –> 2010-2019 4.553818 9632.8366 9702.2712
obs Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 11007.20899 11045.31873 6.0221845 2000-2009 –> 2010-2019 11.864860 18956.0544 19030.0580
mod Atlantic (27,27.25] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 5606.85277 5642.83594 2.2187195 2000-2009 –> 2010-2019 4.480659 9550.2658 9619.7005
obs Atlantic (27,27.25] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 11018.39322 11056.50296 6.0419107 2000-2009 –> 2010-2019 11.918794 18981.7864 19055.7900
obs Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 13196.10668 13235.55381 5.8424207 2000-2009 –> 2010-2019 10.523226 21889.9644 21966.4388
mod Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate 11037.44840 11071.19904 3.5662251 2000-2009 –> 2010-2019 6.894362 20706.5068 20773.3691
mod Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 10896.50685 10935.88260 3.4439767 2000-2009 –> 2010-2019 6.734166 20525.3187 20603.3247
mod Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 10909.34310 10948.71885 3.4547813 2000-2009 –> 2010-2019 6.766809 20562.5259 20640.5319
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 13181.02845 13220.47558 5.8211808 2000-2009 –> 2010-2019 10.423390 21825.2699 21901.7443
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 11157.91492 11197.29067 3.6708234 2000-2009 –> 2010-2019 6.936171 20758.8051 20836.8111
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 13179.60585 13219.05298 5.8191808 2000-2009 –> 2010-2019 10.396382 21807.8824 21884.3568
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + silicate 11009.16484 11048.54059 3.5399684 2000-2009 –> 2010-2019 6.788171 20590.6612 20668.6672
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + phosphate 13247.12289 13280.93471 5.9177189 2000-2009 –> 2010-2019 10.534579 21898.6771 21964.2266
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + phosphate_star 13250.53296 13284.34479 5.9225953 2000-2009 –> 2010-2019 10.516817 21887.6848 21953.2343
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 13244.02206 13283.46919 5.9104323 2000-2009 –> 2010-2019 10.525341 21896.3376 21972.8120
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 10888.67413 10928.04988 3.4374004 2000-2009 –> 2010-2019 7.102091 20914.6161 20992.6221
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 13247.22390 13286.67102 5.9150051 2000-2009 –> 2010-2019 10.509118 21886.3066 21962.7810
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 11153.47642 11192.85217 3.6668498 2000-2009 –> 2010-2019 6.799456 20601.4783 20679.4843
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 13176.56408 13216.01121 5.8149069 2000-2009 –> 2010-2019 10.397240 21808.1244 21884.5988
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 10866.82116 10906.19691 3.4191189 2000-2009 –> 2010-2019 6.956715 20762.7308 20840.7368
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + phosphate + phosphate_star 13247.29432 13281.10614 5.9179639 2000-2009 –> 2010-2019 10.516001 21886.8787 21952.4282
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 10889.95116 10929.32691 3.4384717 2000-2009 –> 2010-2019 7.105593 20918.3358 20996.3418
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 13244.17738 13283.62451 5.9106540 2000-2009 –> 2010-2019 10.508214 21885.4577 21961.9321
mod Atlantic (27.25,27.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 11164.18404 11203.55979 3.6764434 2000-2009 –> 2010-2019 6.802377 20604.3295 20682.3355
obs Atlantic (27.5,27.75] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 14154.43379 14194.67905 5.0965750 2000-2009 –> 2010-2019 10.367359 25278.9439 25357.6696
obs Atlantic (27.5,27.75] cstar_tref ~ aou + silicate + phosphate + phosphate_star 14164.34209 14198.83803 5.1096719 2000-2009 –> 2010-2019 10.381599 25287.6346 25355.1138
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + nitrate + phosphate 18254.16992 18292.15510 2.1791029 2000-2009 –> 2010-2019 4.614477 33879.8926 33954.6333
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 17817.03319 17861.34924 2.0668077 2000-2009 –> 2010-2019 4.495447 33426.0279 33513.2255
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate 18566.60089 18604.58607 2.2626927 2000-2009 –> 2010-2019 4.704641 34210.5524 34285.2931
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 18533.78389 18578.09994 2.2532210 2000-2009 –> 2010-2019 4.687978 34159.7918 34246.9894
obs Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 14166.00522 14206.25047 5.1093009 2000-2009 –> 2010-2019 10.305437 25239.0802 25317.8059
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 18142.82673 18187.14278 2.1495477 2000-2009 –> 2010-2019 4.569247 33726.8858 33814.0833
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 18565.10716 18609.42321 2.2617404 2000-2009 –> 2010-2019 4.684031 34156.4016 34243.5991
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 14171.67298 14211.91824 5.1155458 2000-2009 –> 2010-2019 10.319042 25249.8519 25328.5776
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 18373.52176 18417.83780 2.2101316 2000-2009 –> 2010-2019 4.590060 33845.5133 33932.7109
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 14200.49101 14240.73627 5.1474162 2000-2009 –> 2010-2019 10.378532 25297.7594 25376.4851
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 18456.82122 18501.13727 2.2324243 2000-2009 –> 2010-2019 4.621906 33955.9034 34043.1010
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 18483.77561 18528.09166 2.2396859 2000-2009 –> 2010-2019 4.627248 33977.4238 34064.6213
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 18378.92441 18423.24046 2.2115706 2000-2009 –> 2010-2019 4.593264 33855.9302 33943.1277
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 14180.62274 14220.86800 5.1254223 2000-2009 –> 2010-2019 10.341506 25267.5141 25346.2398
obs Atlantic (27.5,27.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 14161.11715 14201.36240 5.1039213 2000-2009 –> 2010-2019 10.391690 25297.2289 25375.9546
obs Atlantic (27.5,27.75] cstar_tref ~ temp + aou + silicate + phosphate 14169.93271 14204.42865 5.1158321 2000-2009 –> 2010-2019 10.404207 25304.4574 25371.9366
obs Atlantic (27.5,27.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 14170.93733 14211.18259 5.1147348 2000-2009 –> 2010-2019 10.426498 25323.3749 25402.1006
obs Atlantic (27.5,27.75] cstar_tref ~ temp + silicate + phosphate + phosphate_star 14179.29205 14213.78799 5.1261617 2000-2009 –> 2010-2019 10.438094 25329.8446 25397.3237
obs Atlantic (27.75,27.85] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 9247.35421 9284.49523 5.3740583 2000-2009 –> 2010-2019 10.551203 16974.4724 17047.5799
obs Atlantic (27.75,27.85] cstar_tref ~ aou + silicate + phosphate + phosphate_star 9245.35546 9277.19062 5.3740606 2000-2009 –> 2010-2019 10.576026 16982.5166 17045.1802
mod Atlantic (27.75,27.85] cstar_tref ~ sal + aou + nitrate + silicate 5218.81645 5250.78707 1.3370373 2000-2009 –> 2010-2019 2.574530 9490.4031 9553.4175
obs Atlantic (27.75,27.85] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 9234.11570 9271.25672 5.3502213 2000-2009 –> 2010-2019 10.543463 16969.0506 17042.1582
mod Atlantic (27.75,27.85] cstar_tref ~ sal + aou + silicate + phosphate 5005.55001 5037.52063 1.2466264 2000-2009 –> 2010-2019 2.433057 9167.1564 9230.1708
obs Atlantic (27.75,27.85] cstar_tref ~ sal + aou + silicate + phosphate 9232.17728 9264.01244 5.3503319 2000-2009 –> 2010-2019 10.571401 16978.5686 17041.2322
obs Atlantic (27.75,27.85] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 9235.50669 9272.64771 5.3527209 2000-2009 –> 2010-2019 10.577343 16985.6109 17058.7185
mod Atlantic (27.75,27.85] cstar_tref ~ sal + nitrate + silicate + phosphate_star 5614.05513 5646.02576 1.5222855 2000-2009 –> 2010-2019 2.887931 10142.8310 10205.8453
mod Atlantic (27.75,27.85] cstar_tref ~ sal + silicate + phosphate + phosphate_star 5605.72893 5637.69955 1.5181301 2000-2009 –> 2010-2019 2.889649 10145.7072 10208.7216
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + silicate 5176.35987 5213.65893 1.3176649 2000-2009 –> 2010-2019 2.551715 9442.6768 9516.1936
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 4889.44603 4926.74509 1.1992154 2000-2009 –> 2010-2019 2.380928 9042.6545 9116.1713
obs Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 9233.08321 9270.22423 5.3483667 2000-2009 –> 2010-2019 10.558586 16976.2367 17049.3442
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5017.76500 5055.06406 1.2508140 2000-2009 –> 2010-2019 2.471651 9255.9858 9329.5026
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 5022.69525 5059.99431 1.2528402 2000-2009 –> 2010-2019 2.472295 9257.9594 9331.4761
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5011.73214 5049.03121 1.2483391 2000-2009 –> 2010-2019 2.465202 9241.4417 9314.9585
obs Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 9234.73158 9271.87260 5.3513279 2000-2009 –> 2010-2019 10.571063 16982.4797 17055.5872
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate_star 5098.95157 5130.92219 1.2854446 2000-2009 –> 2010-2019 2.589720 9507.7222 9570.7366
obs Atlantic (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 9247.40463 9284.54565 5.3741493 2000-2009 –> 2010-2019 10.558231 16977.8944 17051.0019
obs Atlantic (27.75,27.85] cstar_tref ~ temp + aou + silicate + phosphate 9245.40648 9277.24164 5.3741527 2000-2009 –> 2010-2019 10.584388 16986.5673 17049.2309
obs Atlantic (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 9247.75346 9284.89448 5.3747789 2000-2009 –> 2010-2019 10.566703 16982.0498 17055.1573
mod Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate 6320.04286 6351.92218 1.9812787 2000-2009 –> 2010-2019 4.900271 12826.4242 12889.3427
obs Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 19493.03946 19535.43001 5.3175391 2000-2009 –> 2010-2019 10.588320 36163.5580 36247.2633
mod Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 6320.54734 6357.73988 1.9802912 2000-2009 –> 2010-2019 4.821651 12758.6283 12832.0332
obs Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 18891.28521 18933.67576 4.8334219 2000-2009 –> 2010-2019 9.870617 35316.7552 35400.4606
mod Atlantic (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate 6096.58651 6128.46583 1.8390644 2000-2009 –> 2010-2019 4.573669 12432.7919 12495.7104
obs Atlantic (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate_star 18895.77145 18932.10620 4.8383976 2000-2009 –> 2010-2019 9.875868 35319.5361 35391.2836
obs Atlantic (27.85,27.95] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 19291.40586 19333.79641 5.1501487 2000-2009 –> 2010-2019 10.343906 35882.3419 35966.0472
obs Atlantic (27.85,27.95] cstar_tref ~ sal + silicate + phosphate + phosphate_star 19333.25956 19369.59432 5.1861006 2000-2009 –> 2010-2019 10.383382 35925.8627 35997.6102
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + silicate 6161.73577 6198.92832 1.8781866 2000-2009 –> 2010-2019 4.784424 12658.6966 12732.1015
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + phosphate 6554.98677 6586.86609 2.1426793 2000-2009 –> 2010-2019 4.894509 12907.5683 12970.4867
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 5799.17665 5836.36919 1.6643817 2000-2009 –> 2010-2019 4.398748 12137.1550 12210.5599
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 18729.61686 18772.00741 4.7110430 2000-2009 –> 2010-2019 9.717472 35121.9665 35205.6719
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5902.84144 5940.03399 1.7228995 2000-2009 –> 2010-2019 4.482087 12264.3868 12337.7917
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 18710.87410 18753.26464 4.6970571 2000-2009 –> 2010-2019 9.698725 35098.0797 35181.7851
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 18710.20223 18752.59278 4.6965565 2000-2009 –> 2010-2019 9.741900 35144.4094 35228.1148
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 5911.07328 5948.26582 1.7276335 2000-2009 –> 2010-2019 4.531073 12314.1142 12387.5191
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate 18734.37145 18770.70620 4.7160934 2000-2009 –> 2010-2019 9.762811 35168.0511 35239.7986
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5892.56694 5929.75948 1.7170089 2000-2009 –> 2010-2019 4.467663 12246.0340 12319.4389
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 18728.34610 18770.73665 4.7100934 2000-2009 –> 2010-2019 9.715367 35119.4479 35203.1533
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate_star 5992.97747 6024.85679 1.7766341 2000-2009 –> 2010-2019 4.809185 12598.8952 12661.8137
obs Atlantic (27.95,28.05] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 19948.27457 19990.79276 5.3983247 2000-2009 –> 2010-2019 12.105589 37715.7357 37799.4880
obs Atlantic (27.95,28.05] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 19533.50603 19576.02422 5.0605894 2000-2009 –> 2010-2019 11.622708 37184.0528 37267.8051
obs Atlantic (27.95,28.05] cstar_tref ~ sal + aou + silicate + phosphate_star 19619.69357 19656.13773 5.1305832 2000-2009 –> 2010-2019 11.707027 37279.8939 37351.6816
obs Atlantic (27.95,28.05] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 19669.58294 19712.10112 5.1689974 2000-2009 –> 2010-2019 11.776001 37356.5584 37440.3107
obs Atlantic (27.95,28.05] cstar_tref ~ sal + silicate + phosphate + phosphate_star 19798.38040 19834.82455 5.2753877 2000-2009 –> 2010-2019 11.896274 37494.5729 37566.3606
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate 9623.00203 9655.92171 3.5777280 2000-2009 –> 2010-2019 7.378110 17467.4924 17531.9710
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate_star 9652.83764 9685.75732 3.6077704 2000-2009 –> 2010-2019 7.406974 17496.4455 17560.9241
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 9486.52795 9524.93425 3.4415363 2000-2009 –> 2010-2019 7.241058 17332.3741 17407.5991
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 19651.90629 19694.42447 5.1547848 2000-2009 –> 2010-2019 11.753947 37332.5353 37416.2875
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 9498.29391 9536.70020 3.4529039 2000-2009 –> 2010-2019 7.249696 17342.0962 17417.3213
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 19597.33727 19639.85546 5.1111556 2000-2009 –> 2010-2019 11.687058 37259.0974 37342.8497
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate + phosphate_star 9717.49295 9750.41263 3.6737423 2000-2009 –> 2010-2019 7.459851 17551.2816 17615.7602
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 19681.13449 19723.65267 5.1783064 2000-2009 –> 2010-2019 11.893625 37455.0090 37538.7613
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 9472.53996 9510.94625 3.4280705 2000-2009 –> 2010-2019 7.178174 17281.1532 17356.3782
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + phosphate + phosphate_star 9648.03183 9680.95151 3.6029143 2000-2009 –> 2010-2019 7.402116 17491.6385 17556.1171
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + phosphate_star 9730.69708 9758.13014 3.6894304 2000-2009 –> 2010-2019 7.491521 17574.4654 17628.1976
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate 19784.25804 19820.70220 5.2637960 2000-2009 –> 2010-2019 11.983606 37559.7054 37631.4931
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 9498.54059 9536.94688 3.4531426 2000-2009 –> 2010-2019 7.249940 17342.3468 17417.5718
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 19636.15567 19678.67386 5.1421537 2000-2009 –> 2010-2019 11.732637 37309.7525 37393.5048
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate_star 9497.14571 9530.06539 3.4537283 2000-2009 –> 2010-2019 7.250533 17338.9578 17403.4364
obs Atlantic (28.05,28.1] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 14408.24031 14448.53368 5.2699477 2000-2009 –> 2010-2019 11.607348 26001.3741 26080.0305
obs Atlantic (28.05,28.1] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 14352.62719 14392.92056 5.2075888 2000-2009 –> 2010-2019 11.544787 25945.6476 26024.3040
obs Atlantic (28.05,28.1] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 14355.74805 14396.04141 5.2110686 2000-2009 –> 2010-2019 11.535521 25941.6294 26020.2858
obs Atlantic (28.05,28.1] cstar_tref ~ sal + silicate + phosphate + phosphate_star 14441.61926 14476.15643 5.3100061 2000-2009 –> 2010-2019 11.657966 26038.6566 26106.0763
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 14474.75327 14515.04664 5.3455099 2000-2009 –> 2010-2019 11.650708 26049.8226 26128.4790
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate 3563.60761 3593.28370 1.3367062 2000-2009 –> 2010-2019 2.330561 5731.5017 5789.0013
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate_star 3602.99318 3632.66926 1.3622832 2000-2009 –> 2010-2019 2.374519 5798.8518 5856.3514
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 3564.23574 3598.85784 1.3358240 2000-2009 –> 2010-2019 2.311071 5705.2879 5772.3708
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 3601.23609 3635.85819 1.3598224 2000-2009 –> 2010-2019 2.358158 5777.9942 5845.0771
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 14389.10912 14429.40249 5.2484121 2000-2009 –> 2010-2019 11.519834 25945.1320 26023.7883
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 14467.68904 14507.98241 5.3374334 2000-2009 –> 2010-2019 11.645912 26044.6033 26123.2596
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + phosphate_star 3649.46189 3679.13797 1.3930901 2000-2009 –> 2010-2019 2.416327 5861.8166 5919.3162
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + phosphate_star 14476.28534 14510.82252 5.3495527 2000-2009 –> 2010-2019 11.658459 26051.4399 26118.8597
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 3622.64130 3657.26340 1.3739021 2000-2009 –> 2010-2019 2.396317 5835.7683 5902.8512
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 14356.27209 14396.56545 5.2116531 2000-2009 –> 2010-2019 11.502734 25923.3932 26002.0496
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + phosphate + phosphate_star 3591.76480 3621.44088 1.3549420 2000-2009 –> 2010-2019 2.362775 5780.9712 5838.4708
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 3591.65295 3626.27504 1.3535657 2000-2009 –> 2010-2019 2.345419 5758.4697 5825.5526
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 14430.43608 14470.72944 5.2950438 2000-2009 –> 2010-2019 11.592464 26001.1287 26079.7851
mod Atlantic (28.05,28.1] cstar_tref ~ temp + aou + phosphate 3580.55508 3605.28515 1.3489502 2000-2009 –> 2010-2019 2.426778 5870.2250 5918.1414
mod Atlantic (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate 3582.53815 3612.21423 1.3489392 2000-2009 –> 2010-2019 2.402260 5839.1104 5896.6100
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 16565.65015 16606.65522 5.9375891 2000-2009 –> 2010-2019 15.026699 27624.4934 27702.7975
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + phosphate 4574.78913 4606.79119 1.0737907 2000-2009 –> 2010-2019 2.400519 8946.4278 9009.3623
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + phosphate_star 4622.68620 4654.68826 1.0907194 2000-2009 –> 2010-2019 2.425745 9010.2984 9073.2329
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + silicate + phosphate 4544.88193 4582.21766 1.0626595 2000-2009 –> 2010-2019 2.270836 8678.7083 8752.1318
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4601.70110 4639.03683 1.0825625 2000-2009 –> 2010-2019 2.304692 8764.9447 8838.3682
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 16600.13321 16641.13829 5.9773088 2000-2009 –> 2010-2019 15.082670 27664.4177 27742.7218
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 16548.83770 16589.84277 5.9183193 2000-2009 –> 2010-2019 15.001280 27605.6196 27683.9237
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + phosphate_star 4642.47248 4674.47454 1.0977903 2000-2009 –> 2010-2019 2.407328 8980.6974 9043.6319
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 4643.31633 4680.65206 1.0973759 2000-2009 –> 2010-2019 2.376706 8923.7507 8997.1742
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 16574.27801 16615.28308 5.9475024 2000-2009 –> 2010-2019 15.056611 27639.8159 27718.1201
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + phosphate + phosphate_star 4607.55043 4639.55249 1.0853412 2000-2009 –> 2010-2019 2.414708 8984.2802 9047.2146
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + phosphate + phosphate_star 16601.33250 16636.47971 5.9810074 2000-2009 –> 2010-2019 15.084805 27663.0939 27730.2117
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4580.99128 4618.32701 1.0752653 2000-2009 –> 2010-2019 2.286660 8721.6323 8795.0558
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 16600.07715 16641.08223 5.9772441 2000-2009 –> 2010-2019 15.079260 27663.2424 27741.5465
obs Atlantic (28.1,28.15] cstar_tref ~ temp + aou + nitrate + phosphate_star 16592.68025 16627.82746 5.9710101 2000-2009 –> 2010-2019 15.074478 27654.3314 27721.4492
obs Atlantic (28.1,28.15] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 16568.32932 16609.33439 5.9406657 2000-2009 –> 2010-2019 15.043387 27631.7308 27710.0350
mod Atlantic (28.1,28.15] cstar_tref ~ temp + aou + silicate + phosphate 4674.01351 4706.01557 1.1091569 2000-2009 –> 2010-2019 2.395513 8966.4805 9029.4150
obs Atlantic (28.1,28.15] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 16579.61525 16614.76245 5.9559458 2000-2009 –> 2010-2019 15.054175 27639.5130 27706.6308
obs Atlantic (28.1,28.15] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 16556.95095 16597.95602 5.9276106 2000-2009 –> 2010-2019 15.025242 27618.6487 27696.9529
mod Atlantic (28.1,28.15] cstar_tref ~ temp + silicate + phosphate + phosphate_star 4721.12299 4753.12505 1.1263534 2000-2009 –> 2010-2019 2.418217 9024.5344 9087.4689
mod Atlantic (28.15,28.2] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 7262.50579 7304.29613 0.8469004 2000-2009 –> 2010-2019 1.897865 13656.3101 13737.8842
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate 7322.17199 7357.99229 0.8559747 2000-2009 –> 2010-2019 1.898180 13677.6164 13747.5371
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate + silicate 7224.96024 7266.75058 0.8414226 2000-2009 –> 2010-2019 1.879791 13566.3851 13647.9592
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate 5567.35604 5601.05784 5.0988071 2000-2009 –> 2010-2019 12.844655 10191.2630 10256.4633
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5527.34992 5561.05172 4.9880714 2000-2009 –> 2010-2019 12.651826 10137.0858 10202.2862
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + nitrate + phosphate_star 7635.78534 7671.60564 0.9036508 2000-2009 –> 2010-2019 1.964053 14066.4209 14136.3416
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 7598.92364 7640.71399 0.8976017 2000-2009 –> 2010-2019 1.951586 14005.1879 14086.7621
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 5512.91554 5546.61734 4.9487107 2000-2009 –> 2010-2019 12.661592 10131.1500 10196.3503
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5543.73959 5577.44139 5.0331437 2000-2009 –> 2010-2019 12.708584 10155.5019 10220.7022
mod Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + phosphate 7328.41159 7364.23189 0.8568983 2000-2009 –> 2010-2019 1.898225 13680.1951 13750.1158
mod Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + phosphate_star 7361.76815 7397.58845 0.8618526 2000-2009 –> 2010-2019 1.908732 13736.6522 13806.5729
mod Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate 7252.73464 7288.55494 0.8457636 2000-2009 –> 2010-2019 1.896386 13643.1235 13713.0442
obs Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 5528.10073 5561.80253 4.9901274 2000-2009 –> 2010-2019 12.994798 10195.7225 10260.9229
obs Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5499.71927 5533.42107 4.9129980 2000-2009 –> 2010-2019 12.818979 10150.8415 10216.0418
obs Atlantic (28.15,28.2] cstar_tref ~ temp + aou + silicate + phosphate_star 5527.88605 5556.77330 4.9950194 2000-2009 –> 2010-2019 12.922106 10180.5540 10236.4400
mod Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 7412.65376 7448.47406 0.8694657 2000-2009 –> 2010-2019 1.920324 13804.0154 13873.9362
obs Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 5506.16600 5539.86780 4.9304123 2000-2009 –> 2010-2019 12.865795 10162.2253 10227.4256
mod Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + silicate + phosphate_star 7615.97881 7651.79910 0.9005627 2000-2009 –> 2010-2019 1.954693 14020.8476 14090.7683
obs Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + silicate + phosphate_star 5513.87650 5542.76376 4.9567595 2000-2009 –> 2010-2019 12.982937 10183.0668 10238.9528
obs Atlantic (28.15,28.2] cstar_tref ~ temp + silicate + phosphate + phosphate_star 5543.63597 5572.52322 5.0383850 2000-2009 –> 2010-2019 12.986711 10199.8629 10255.7488
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate 10054.18240 10086.57466 0.6873960 2000-2009 –> 2010-2019 1.438751 15899.4478 15961.1081
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + phosphate 9991.70298 10030.57369 0.6828040 2000-2009 –> 2010-2019 1.432052 15824.5068 15898.4991
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 10756.56280 10796.05711 3.1849984 2000-2009 –> 2010-2019 6.499520 18334.4209 18410.8463
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + phosphate_star 9907.06551 9945.93623 0.6768230 2000-2009 –> 2010-2019 1.426088 15739.9881 15813.9805
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate 9944.48357 9983.35429 0.6794607 2000-2009 –> 2010-2019 1.414569 15679.1703 15753.1627
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 10755.14990 10794.64421 3.1839189 2000-2009 –> 2010-2019 6.500710 18334.9857 18411.4111
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 9302.88624 9348.23540 0.6354907 2000-2009 –> 2010-2019 1.358653 14955.1952 15041.5196
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + phosphate + phosphate_star 10757.87903 10791.73130 3.1875335 2000-2009 –> 2010-2019 6.503470 18334.9702 18400.4777
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + silicate + phosphate 10753.73837 10787.59064 3.1843684 2000-2009 –> 2010-2019 6.501470 18331.8452 18397.3527
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 10735.70428 10775.19859 3.1690990 2000-2009 –> 2010-2019 6.482904 18312.9370 18389.3624
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate 9973.25980 10012.13051 0.6814962 2000-2009 –> 2010-2019 1.431376 15810.4035 15884.3958
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate 9909.06960 9954.41877 0.6768233 2000-2009 –> 2010-2019 1.425849 15742.3467 15828.6711
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate 10754.78323 10794.27754 3.1836388 2000-2009 –> 2010-2019 6.502113 18336.0853 18412.5107
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 9903.88089 9949.23006 0.6764584 2000-2009 –> 2010-2019 1.425698 15738.6269 15824.9513
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + silicate 9865.32815 9910.67732 0.6737528 2000-2009 –> 2010-2019 1.406524 15585.6123 15671.9367
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + silicate + phosphate 10726.14194 10765.63625 3.1618367 2000-2009 –> 2010-2019 6.478924 18306.2355 18382.6609
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 10769.25549 10808.74980 3.1947123 2000-2009 –> 2010-2019 6.469175 18311.9719 18388.3972
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 9949.14292 9994.49209 0.6796486 2000-2009 –> 2010-2019 1.431786 15803.7725 15890.0969
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 10796.83724 10836.33155 3.2159234 2000-2009 –> 2010-2019 6.497700 18346.0024 18422.4277
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 10742.10149 10781.59580 3.1739668 2000-2009 –> 2010-2019 6.455765 18291.2853 18367.7106
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 22564.33777 22607.70510 5.4320205 2000-2009 –> 2010-2019 10.971106 38075.6389 38159.7042
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 29326.98669 29372.14107 5.5518418 2000-2009 –> 2010-2019 11.850167 50317.3410 50405.0310
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + phosphate + phosphate_star 22895.26593 22932.43793 5.6873531 2000-2009 –> 2010-2019 11.582288 38712.7747 38784.8307
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 22674.44283 22717.81016 5.5151688 2000-2009 –> 2010-2019 11.169422 38287.6081 38371.6734
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 29644.78249 29689.93687 5.7436609 2000-2009 –> 2010-2019 12.265770 50859.8436 50947.5335
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate 30107.88575 30146.58951 6.0364045 2000-2009 –> 2010-2019 12.519749 51282.5803 51357.7432
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 22672.83280 22716.20013 5.5139438 2000-2009 –> 2010-2019 11.330320 38425.9322 38509.9975
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 29786.08949 29831.24387 5.8310678 2000-2009 –> 2010-2019 12.151728 50799.2274 50886.9174
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 22775.34402 22818.71136 5.5924837 2000-2009 –> 2010-2019 11.511643 38615.1518 38699.2171
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 29839.39747 29884.55186 5.8643865 2000-2009 –> 2010-2019 12.195707 50863.3811 50951.0711
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 23246.02431 23289.39164 5.9677074 2000-2009 –> 2010-2019 12.158486 39307.9217 39391.9870
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 30026.95247 30072.10686 5.9831332 2000-2009 –> 2010-2019 12.463235 51200.4268 51288.1168
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 23124.84124 23168.20858 5.8687597 2000-2009 –> 2010-2019 11.965611 39111.0607 39195.1260
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 30107.62772 30152.78210 6.0349478 2000-2009 –> 2010-2019 12.651553 51415.2685 51502.9585
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 23322.80334 23366.17067 6.0312601 2000-2009 –> 2010-2019 12.304682 39450.3430 39534.4083
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 30219.35028 30264.50466 6.1074448 2000-2009 –> 2010-2019 12.760534 51562.3819 51650.0719
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 22875.76193 22919.12927 5.6705044 2000-2009 –> 2010-2019 11.681661 38791.9124 38875.9777
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 29967.96662 30013.12100 5.9455306 2000-2009 –> 2010-2019 12.321130 51036.8044 51124.4944
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 23391.79408 23435.16142 6.0889433 2000-2009 –> 2010-2019 12.398775 39547.9795 39632.0448
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 30268.42699 30313.58137 6.1395654 2000-2009 –> 2010-2019 12.808007 51626.2927 51713.9827
mod Indo-Pacific (26,26.5] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 16275.44971 16317.19636 4.0924331 2000-2009 –> 2010-2019 8.541604 30504.1553 30586.5031
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 16566.59765 16608.34430 4.3049869 2000-2009 –> 2010-2019 8.574273 30593.8181 30676.1659
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 15843.22697 15884.97362 3.7960859 2000-2009 –> 2010-2019 7.863082 29633.4660 29715.8138
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 16255.26645 16297.01311 4.0780933 2000-2009 –> 2010-2019 8.358035 30294.6569 30377.0047
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 17496.57081 17538.51432 4.6515221 2000-2009 –> 2010-2019 9.225223 32183.3629 32266.0635
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 17516.02756 17557.97107 4.6668506 2000-2009 –> 2010-2019 9.257720 32221.5235 32304.2241
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + phosphate 17781.27098 17817.22257 4.8825749 2000-2009 –> 2010-2019 9.746417 32773.1005 32843.9867
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 15464.98997 15506.73662 3.5544138 2000-2009 –> 2010-2019 7.633377 29269.5725 29351.9203
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 17359.68892 17401.63243 4.5450972 2000-2009 –> 2010-2019 9.222318 32158.2093 32240.9099
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 15611.91081 15653.65747 3.6464045 2000-2009 –> 2010-2019 7.894354 29614.6712 29697.0190
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 17430.95556 17472.89907 4.6001992 2000-2009 –> 2010-2019 9.313715 32268.0636 32350.7642
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 17519.28826 17561.23178 4.6694244 2000-2009 –> 2010-2019 9.265730 32230.6916 32313.3922
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate_star 17573.79571 17609.74730 4.7142539 2000-2009 –> 2010-2019 9.482247 32466.2686 32537.1548
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 17341.53730 17383.48081 4.5311685 2000-2009 –> 2010-2019 9.245190 32179.1804 32261.8810
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 15654.34245 15696.08910 3.6734123 2000-2009 –> 2010-2019 7.953655 29694.0761 29776.4239
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 17430.37444 17472.31795 4.5997472 2000-2009 –> 2010-2019 9.315025 32269.3484 32352.0490
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate_star 17429.43934 17465.39092 4.6005755 2000-2009 –> 2010-2019 9.335447 32287.1138 32358.0000
mod Indo-Pacific (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 16024.13689 16065.88355 3.9174193 2000-2009 –> 2010-2019 8.195707 30061.6394 30143.9872
mod Indo-Pacific (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 16088.62512 16130.37178 3.9616018 2000-2009 –> 2010-2019 8.249917 30137.5574 30219.9052
mod Indo-Pacific (26,26.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 16135.33848 16177.08513 3.9939172 2000-2009 –> 2010-2019 8.299789 30204.2163 30286.5641
obs Indo-Pacific (26.5,26.75] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 15823.27289 15858.23155 5.6731516 2000-2009 –> 2010-2019 11.054791 32411.9951 32482.2973
mod Indo-Pacific (26.5,26.75] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 14207.61089 14248.76594 3.5509330 2000-2009 –> 2010-2019 7.888067 27036.1271 27117.2189
mod Indo-Pacific (26.5,26.75] cstar_tref ~ aou + silicate + phosphate + phosphate_star 14249.28531 14284.56106 3.5804045 2000-2009 –> 2010-2019 8.196957 27353.0103 27422.5177
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 15669.53485 15710.31996 5.4995803 2000-2009 –> 2010-2019 10.527400 31896.8302 31978.8494
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 15694.01637 15734.80147 5.5265091 2000-2009 –> 2010-2019 10.612453 31982.7362 32064.7554
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + phosphate 15875.07821 15910.03687 5.7320949 2000-2009 –> 2010-2019 11.108048 32458.1512 32528.4533
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 14030.29375 14071.44879 3.4337502 2000-2009 –> 2010-2019 7.800689 26889.2175 26970.3094
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 15822.01370 15862.79880 5.6694637 2000-2009 –> 2010-2019 11.036370 32398.0861 32480.1053
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 14200.44652 14241.60156 3.5461217 2000-2009 –> 2010-2019 8.026364 27173.1003 27254.1922
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 15696.50745 15737.29255 5.5292566 2000-2009 –> 2010-2019 10.624475 31994.9635 32076.9827
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate_star 15717.15620 15752.11486 5.5542993 2000-2009 –> 2010-2019 10.762019 32130.3232 32200.6253
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 14714.63496 14755.79000 3.9085442 2000-2009 –> 2010-2019 8.522189 27817.5629 27898.6548
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 15702.81074 15743.59585 5.5362148 2000-2009 –> 2010-2019 10.675714 32047.5097 32129.5289
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 14225.65060 14266.80564 3.5630767 2000-2009 –> 2010-2019 8.051865 27206.7657 27287.8576
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 15858.07050 15898.85560 5.7103974 2000-2009 –> 2010-2019 11.107440 32464.0669 32546.0861
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 14377.49507 14418.65011 3.6669530 2000-2009 –> 2010-2019 8.098108 27301.2341 27382.3260
obs Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 15804.87130 15845.65640 5.6501057 2000-2009 –> 2010-2019 11.069539 32432.9919 32515.0111
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 14594.00988 14635.16492 3.8203293 2000-2009 –> 2010-2019 8.347810 27613.2325 27694.3244
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + silicate + phosphate 14422.88984 14458.16559 3.6999916 2000-2009 –> 2010-2019 8.451040 27654.1188 27723.6261
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 14636.30198 14677.45702 3.8510291 2000-2009 –> 2010-2019 8.395271 27671.9316 27753.0234
mod Indo-Pacific (26.75,27] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 18480.59192 18517.29580 3.8035051 2000-2009 –> 2010-2019 7.730892 36841.6265 36914.9259
obs Indo-Pacific (26.75,27] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 24818.51808 24856.30932 5.3014749 2000-2009 –> 2010-2019 10.342775 49408.7783 49484.4040
mod Indo-Pacific (26.75,27] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 17494.48692 17537.30811 3.2822615 2000-2009 –> 2010-2019 7.202414 35845.3812 35930.8971
obs Indo-Pacific (26.75,27] cstar_tref ~ aou + phosphate + phosphate_star 25127.47751 25158.97021 5.5106405 2000-2009 –> 2010-2019 10.729975 49996.6478 50059.6692
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + phosphate 17817.77401 17860.59520 3.4444205 2000-2009 –> 2010-2019 7.145457 35789.9727 35875.4886
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + phosphate 24995.63247 25039.72224 5.4182685 2000-2009 –> 2010-2019 10.556805 49742.5242 49830.7541
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 17992.95572 18035.77691 3.5356127 2000-2009 –> 2010-2019 7.316242 36105.2463 36190.7622
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate 25242.83662 25286.92640 5.5875360 2000-2009 –> 2010-2019 10.752330 50030.9823 50119.2122
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 18014.98802 18057.80920 3.5472514 2000-2009 –> 2010-2019 7.331626 36133.7981 36219.3140
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + phosphate_star 18501.40033 18538.10421 3.8153290 2000-2009 –> 2010-2019 7.622628 36657.9722 36731.2716
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 17931.95600 17974.77719 3.5035881 2000-2009 –> 2010-2019 7.306789 36083.4387 36168.9547
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 25127.70489 25171.79466 5.5080541 2000-2009 –> 2010-2019 10.595384 49793.5330 49881.7629
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + phosphate 25061.92172 25099.71295 5.4645087 2000-2009 –> 2010-2019 10.665170 49904.0830 49979.7087
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 17477.75139 17520.57257 3.2740780 2000-2009 –> 2010-2019 7.203441 35844.0972 35929.6132
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 25033.12198 25077.21176 5.4436049 2000-2009 –> 2010-2019 10.529286 49696.3273 49784.5572
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 17661.54828 17704.36947 3.3650820 2000-2009 –> 2010-2019 7.393626 36192.0187 36277.5346
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 25175.40755 25219.49732 5.5408478 2000-2009 –> 2010-2019 10.705531 49963.3801 50051.6100
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 17689.61913 17732.44032 3.3792017 2000-2009 –> 2010-2019 7.412918 36228.5371 36314.0530
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 25175.94362 25220.03340 5.5412174 2000-2009 –> 2010-2019 10.704548 49961.7947 50050.0246
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate_star 25174.27176 25212.06299 5.5414437 2000-2009 –> 2010-2019 10.706218 49960.3871 50036.0127
mod Indo-Pacific (27,27.25] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 21219.81726 21258.19419 2.6504258 2000-2009 –> 2010-2019 5.832300 43167.3625 43243.8773
mod Indo-Pacific (27,27.25] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 21171.76387 21216.53695 2.6354948 2000-2009 –> 2010-2019 5.646945 42652.6243 42741.8916
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 31763.15282 31809.14264 4.9172212 2000-2009 –> 2010-2019 9.670478 62924.6129 63016.5379
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 31968.91173 32014.90156 5.0141384 2000-2009 –> 2010-2019 9.875496 63365.5930 63457.5180
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate 31764.20164 31803.62149 4.9186435 2000-2009 –> 2010-2019 9.672963 62925.9986 63004.7914
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate_star 31973.44786 32012.86772 5.0172481 2000-2009 –> 2010-2019 9.897799 63409.3436 63488.1365
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 31897.40738 31943.39721 4.9802435 2000-2009 –> 2010-2019 9.800251 63204.7370 63296.6620
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + silicate + phosphate + phosphate_star 31903.45974 31942.87959 4.9840490 2000-2009 –> 2010-2019 9.825917 63256.1215 63334.9144
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate 21163.76450 21208.53758 2.6331164 2000-2009 –> 2010-2019 5.798106 43068.0112 43157.2786
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 21562.41158 21607.18467 2.7542966 2000-2009 –> 2010-2019 6.015743 43722.2564 43811.5238
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate 31763.35986 31809.34969 4.9173178 2000-2009 –> 2010-2019 9.669285 62921.9813 63013.9063
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 31940.08359 31986.07342 5.0004454 2000-2009 –> 2010-2019 9.811397 63227.7420 63319.6670
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 21590.06164 21634.83472 2.7629056 2000-2009 –> 2010-2019 6.032191 43770.3457 43859.6130
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 22054.68004 22099.45312 2.9116586 2000-2009 –> 2010-2019 6.153245 44162.5214 44251.7887
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 32030.86655 32076.85638 5.0436930 2000-2009 –> 2010-2019 9.896460 63409.0457 63500.9707
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 31895.44650 31941.43633 4.9793172 2000-2009 –> 2010-2019 9.779072 63158.7324 63250.6574
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + phosphate 21549.85686 21588.23379 2.7510174 2000-2009 –> 2010-2019 5.978892 43619.6088 43696.1236
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 21550.48717 21595.26025 2.7505922 2000-2009 –> 2010-2019 5.870303 43332.0500 43421.3173
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 21910.17498 21954.94807 2.8645551 2000-2009 –> 2010-2019 6.102913 44009.5329 44098.8002
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 21930.82401 21975.59709 2.8712390 2000-2009 –> 2010-2019 6.113355 44040.0566 44129.3239
obs Indo-Pacific (27.25,27.5] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 33845.19506 33885.28011 4.2786478 2000-2009 –> 2010-2019 8.643171 66593.9153 66673.8464
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 33842.33998 33889.10587 4.2768845 2000-2009 –> 2010-2019 8.584535 66444.6087 66537.8616
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 29322.33252 29368.84177 3.1975462 2000-2009 –> 2010-2019 6.264175 57017.0704 57109.8033
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 33535.17846 33581.94435 4.1667884 2000-2009 –> 2010-2019 8.405867 65955.8312 66049.0841
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + silicate + phosphate 33540.11288 33580.19793 4.1692423 2000-2009 –> 2010-2019 8.419343 65988.1562 66068.0872
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 33782.41455 33829.18044 4.2551793 2000-2009 –> 2010-2019 8.589206 66453.7748 66547.0278
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + silicate + phosphate + phosphate_star 33783.00019 33823.08524 4.2561136 2000-2009 –> 2010-2019 8.633354 66564.6516 66644.5827
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 28895.96523 28942.47448 3.0796979 2000-2009 –> 2010-2019 6.265153 57005.0816 57097.8144
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 29156.85356 29203.36281 3.1512815 2000-2009 –> 2010-2019 6.405808 57499.7948 57592.5277
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 33516.90581 33563.67170 4.1603289 2000-2009 –> 2010-2019 8.355204 65818.9188 65912.1717
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 33730.51627 33777.28216 4.2364707 2000-2009 –> 2010-2019 8.483518 66172.4249 66265.6779
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 29176.22047 29222.72971 3.1566613 2000-2009 –> 2010-2019 6.417235 57539.3958 57632.1287
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 28783.92600 28830.43525 3.0494575 2000-2009 –> 2010-2019 6.100586 56423.4328 56516.1656
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 33673.28718 33720.05307 4.2159357 2000-2009 –> 2010-2019 8.467715 66127.8002 66221.0531
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 33677.30449 33724.07038 4.2173739 2000-2009 –> 2010-2019 8.455075 66094.2810 66187.5339
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + phosphate 28897.75259 28937.61766 3.0807254 2000-2009 –> 2010-2019 6.330915 57224.1583 57303.6436
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 28574.59934 28621.10859 2.9937516 2000-2009 –> 2010-2019 6.126253 56500.9965 56593.7293
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 28856.08005 28902.58930 3.0688983 2000-2009 –> 2010-2019 6.289536 57084.9260 57177.6588
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 28857.57624 28904.08549 3.0693028 2000-2009 –> 2010-2019 6.291056 57090.1959 57182.9287
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + nitrate + silicate + phosphate_star 28856.04436 28895.90942 3.0694293 2000-2009 –> 2010-2019 6.292135 57089.8867 57169.3720
obs Indo-Pacific (27.5,27.75] cstar_tref ~ aou + silicate + phosphate + phosphate_star 38445.04909 38486.04111 4.0003966 2000-2009 –> 2010-2019 7.670178 72668.3966 72749.8698
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 38973.84058 39021.66461 4.1572153 2000-2009 –> 2010-2019 7.798070 73099.6174 73194.6696
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 32816.67256 32863.97704 3.1895758 2000-2009 –> 2010-2019 6.102733 64977.9072 65072.6254
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 38211.31244 38259.13646 3.9321504 2000-2009 –> 2010-2019 7.770644 73002.1941 73097.2462
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 32990.12385 33037.42834 3.2333671 2000-2009 –> 2010-2019 6.146875 65152.9130 65247.6311
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate_star 33255.69880 33296.24550 3.3021039 2000-2009 –> 2010-2019 6.342231 65966.1251 66047.3121
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 32943.41655 32990.72104 3.2215161 2000-2009 –> 2010-2019 6.130692 65086.9820 65181.7001
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 33182.80178 33230.10627 3.2827178 2000-2009 –> 2010-2019 6.286745 65740.8925 65835.6107
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 37999.28679 38047.11082 3.8717637 2000-2009 –> 2010-2019 7.455808 71927.1934 72022.2455
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 33242.43873 33289.74322 3.2981447 2000-2009 –> 2010-2019 6.315989 65859.8185 65954.5366
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 38097.09123 38144.91525 3.8995031 2000-2009 –> 2010-2019 7.516423 72139.8818 72234.9340
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 33245.19912 33292.50360 3.2988605 2000-2009 –> 2010-2019 6.304498 65810.2107 65904.9288
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 32932.00065 32979.30513 3.2186261 2000-2009 –> 2010-2019 6.181675 65312.6330 65407.3512
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 38062.43522 38110.25924 3.8896513 2000-2009 –> 2010-2019 7.500988 72085.7897 72180.8419
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate 33282.01200 33322.55870 3.3089418 2000-2009 –> 2010-2019 6.345375 65976.7291 66057.9161
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 33213.40447 33260.70896 3.2906251 2000-2009 –> 2010-2019 6.300979 65798.6769 65893.3950
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 38082.87470 38130.69873 3.8954587 2000-2009 –> 2010-2019 7.510970 72120.7644 72215.8165
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate_star 38374.19199 38415.18401 3.9797598 2000-2009 –> 2010-2019 7.604914 72443.5955 72525.0688
obs Indo-Pacific (27.5,27.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 38518.83711 38566.66114 4.0214137 2000-2009 –> 2010-2019 7.737155 72900.7808 72995.8329
obs Indo-Pacific (27.5,27.75] cstar_tref ~ temp + aou + silicate + phosphate 38519.68244 38560.67446 4.0222490 2000-2009 –> 2010-2019 7.745275 72924.2684 73005.7417
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + phosphate 11508.09371 11548.81990 2.4442941 2000-2009 –> 2010-2019 4.677410 22633.7749 22715.2694
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 11548.86110 11589.58729 2.4644264 2000-2009 –> 2010-2019 4.705468 22692.2568 22773.7513
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + phosphate 11836.82221 11871.73038 2.6124836 2000-2009 –> 2010-2019 4.894269 23068.3059 23138.1584
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 11592.18751 11632.91370 2.4860041 2000-2009 –> 2010-2019 4.767382 22824.7786 22906.2731
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 21251.39668 21295.12722 3.9079393 2000-2009 –> 2010-2019 6.998401 37649.3666 37735.6284
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 11632.50922 11673.23542 2.5062552 2000-2009 –> 2010-2019 4.795039 22881.3037 22962.7982
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 21230.24171 21273.97225 3.8971248 2000-2009 –> 2010-2019 6.989090 37631.3392 37717.6010
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 11574.17047 11614.89667 2.4770083 2000-2009 –> 2010-2019 4.723656 22730.0587 22811.5532
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 21228.04957 21271.78011 3.8960059 2000-2009 –> 2010-2019 6.996735 37647.3525 37733.6143
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 11649.42731 11690.15351 2.5148011 2000-2009 –> 2010-2019 4.805812 22903.0866 22984.5812
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 21229.68177 21273.41231 3.8968390 2000-2009 –> 2010-2019 6.989208 37631.6200 37717.8818
mod Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 11501.08964 11541.81583 2.4408519 2000-2009 –> 2010-2019 4.838039 22981.2631 23062.7576
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 21256.10410 21299.83464 3.9103498 2000-2009 –> 2010-2019 7.000619 37653.6736 37739.9354
mod Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 11551.75747 11592.48367 2.4658630 2000-2009 –> 2010-2019 4.871737 23050.0170 23131.5115
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 21237.59885 21281.32939 3.9008824 2000-2009 –> 2010-2019 6.992839 37638.6795 37724.9413
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + silicate + phosphate_star 21242.95621 21280.43953 3.9046438 2000-2009 –> 2010-2019 6.996626 37642.0913 37716.0299
mod Indo-Pacific (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 11586.02070 11626.74689 2.4829214 2000-2009 –> 2010-2019 4.894736 23096.6099 23178.1044
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 21237.81696 21281.54750 3.9009939 2000-2009 –> 2010-2019 6.993339 37639.7060 37725.9677
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate_star 21236.01942 21273.50274 3.9010973 2000-2009 –> 2010-2019 7.001829 37653.3289 37727.2676
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + silicate + phosphate + phosphate_star 21241.98509 21279.46841 3.9041471 2000-2009 –> 2010-2019 6.996530 37641.9529 37715.8916
obs Indo-Pacific (27.85,27.95] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 27294.72435 27340.64209 3.3057399 2000-2009 –> 2010-2019 6.049687 48697.5086 48788.1583
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 13697.76797 13740.20500 2.0903758 2000-2009 –> 2010-2019 4.107834 26389.7724 26474.2283
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 27321.16032 27367.07807 3.3141260 2000-2009 –> 2010-2019 6.059194 48727.5414 48818.1910
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate 13091.22350 13127.59810 1.9004296 2000-2009 –> 2010-2019 3.852496 25584.2546 25656.6453
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + phosphate 13092.97322 13135.41025 1.9003547 2000-2009 –> 2010-2019 3.847536 25573.0228 25657.4787
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 13093.13245 13135.56948 1.9004024 2000-2009 –> 2010-2019 3.848958 25577.4009 25661.8568
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + silicate 13075.89524 13118.33227 1.8952474 2000-2009 –> 2010-2019 3.815757 25473.4949 25557.9508
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + phosphate 13720.78229 13757.15688 2.0986318 2000-2009 –> 2010-2019 4.151947 26516.1041 26588.4949
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 13572.27702 13614.71405 2.0494450 2000-2009 –> 2010-2019 4.096550 26351.4897 26435.9456
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 27275.19083 27321.10858 3.2995570 2000-2009 –> 2010-2019 6.048034 48692.5043 48783.1539
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 13590.07118 13632.50821 2.0551997 2000-2009 –> 2010-2019 4.107835 26385.4125 26469.8683
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 27266.11011 27312.02785 3.2966866 2000-2009 –> 2010-2019 6.056266 48718.9288 48809.5785
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 13137.26582 13179.70285 1.9136648 2000-2009 –> 2010-2019 3.871308 25649.3494 25733.8053
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 13606.13452 13648.57155 2.0604085 2000-2009 –> 2010-2019 4.115259 26407.9242 26492.3800
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 27265.99624 27311.91398 3.2966506 2000-2009 –> 2010-2019 6.057903 48724.1542 48814.8038
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 27257.84336 27303.76111 3.2940757 2000-2009 –> 2010-2019 6.042230 48674.1216 48764.7712
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 27254.63694 27300.55468 3.2930636 2000-2009 –> 2010-2019 6.052769 48707.8596 48798.5092
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + silicate + phosphate 27280.28393 27319.64199 3.3018008 2000-2009 –> 2010-2019 6.050947 48697.7408 48775.4405
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + silicate + phosphate_star 27273.27113 27312.62919 3.2995824 2000-2009 –> 2010-2019 6.059738 48725.9308 48803.6305
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 27254.64528 27300.56303 3.2930662 2000-2009 –> 2010-2019 6.054043 48711.9242 48802.5739
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 46540.40836 46582.88291 3.4341325 2000-2009 –> 2010-2019 6.285438 83247.2686 83331.2293
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 46540.74219 46590.29583 3.4338063 2000-2009 –> 2010-2019 6.272719 83184.8030 83282.7571
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + phosphate + phosphate_star 46697.22671 46732.62217 3.4653685 2000-2009 –> 2010-2019 6.320601 83422.5589 83492.5261
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + silicate + phosphate + phosphate_star 46696.62357 46739.09812 3.4648542 2000-2009 –> 2010-2019 6.312837 83386.1396 83470.1003
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 46542.38881 46591.94245 3.4341286 2000-2009 –> 2010-2019 6.285407 83251.1053 83349.0594
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + aou + phosphate + phosphate_star 46684.72192 46727.19647 3.4625040 2000-2009 –> 2010-2019 6.317583 83411.2543 83495.2150
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate 17705.86068 17744.32100 1.7342019 2000-2009 –> 2010-2019 3.459987 30495.7541 30570.7366
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate + silicate 17579.21239 17624.08276 1.7095456 2000-2009 –> 2010-2019 3.388489 30192.1319 30279.6114
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate 18006.59273 18045.05305 1.7932354 2000-2009 –> 2010-2019 3.620024 31166.4181 31241.4006
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate_star 18009.16190 18047.62222 1.7937483 2000-2009 –> 2010-2019 3.624689 31183.7525 31258.7350
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 18001.98447 18046.85484 1.7919169 2000-2009 –> 2010-2019 3.588758 31056.3000 31143.7796
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 46780.75829 46830.31193 3.4811172 2000-2009 –> 2010-2019 6.338906 83523.4064 83621.3605
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 18005.50024 18050.37061 1.7926182 2000-2009 –> 2010-2019 3.595290 31080.8883 31168.3678
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 18226.97157 18271.84194 1.8373586 2000-2009 –> 2010-2019 3.678812 31440.8000 31528.2795
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + phosphate + phosphate_star 18015.76195 18054.22227 1.7950665 2000-2009 –> 2010-2019 3.628026 31197.5179 31272.5004
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 18011.88150 18056.75187 1.7938920 2000-2009 –> 2010-2019 3.598097 31092.7971 31180.2766
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate_star 18283.25709 18321.71740 1.8493176 2000-2009 –> 2010-2019 3.707882 31555.2422 31630.2246
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 46629.05403 46678.60767 3.4511387 2000-2009 –> 2010-2019 6.303219 83341.9555 83439.9097
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 46743.33394 46792.88759 3.4736976 2000-2009 –> 2010-2019 6.353829 83601.8328 83699.7869
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 46727.92856 46777.48221 3.4706480 2000-2009 –> 2010-2019 6.348925 83576.8457 83674.7999
obs Indo-Pacific (28.05,28.1] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 35851.83862 35893.09370 2.9591219 2000-2009 –> 2010-2019 5.687425 58713.9735 58793.9808
obs Indo-Pacific (28.05,28.1] cstar_tref ~ aou + phosphate + phosphate_star 35925.95238 35960.33161 2.9748988 2000-2009 –> 2010-2019 5.704860 58791.8213 58858.4941
obs Indo-Pacific (28.05,28.1] cstar_tref ~ aou + silicate + phosphate + phosphate_star 35902.59864 35943.85372 2.9696342 2000-2009 –> 2010-2019 5.698618 58767.0878 58847.0952
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 35868.25103 35916.38195 2.9621029 2000-2009 –> 2010-2019 5.702655 58774.6404 58867.9824
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + phosphate 35866.87882 35915.00974 2.9618190 2000-2009 –> 2010-2019 5.716460 58821.6313 58914.9733
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + silicate 9290.27584 9331.27821 1.4554456 2000-2009 –> 2010-2019 3.126988 16737.3839 16817.3177
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 9477.64929 9518.65165 1.5091620 2000-2009 –> 2010-2019 3.087603 16704.3465 16784.2803
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 35866.05469 35914.18561 2.9616484 2000-2009 –> 2010-2019 5.674401 58676.2776 58769.6196
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 9482.90466 9523.90702 1.5106969 2000-2009 –> 2010-2019 3.093559 16720.3587 16800.2925
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 35986.13858 36034.26950 2.9865990 2000-2009 –> 2010-2019 5.722164 58875.3460 58968.6879
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 9432.85999 9473.86236 1.4961442 2000-2009 –> 2010-2019 3.078985 16670.2629 16750.1967
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 35955.41270 36003.54363 2.9801950 2000-2009 –> 2010-2019 5.709901 58824.3983 58917.7402
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 9482.27014 9523.27251 1.5105115 2000-2009 –> 2010-2019 3.093373 16719.7238 16799.6576
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 35987.61819 36035.74911 2.9869077 2000-2009 –> 2010-2019 5.722110 58875.5743 58968.9162
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate_star 9487.79631 9522.94119 1.5127120 2000-2009 –> 2010-2019 3.095575 16723.2539 16791.7686
obs Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 36048.28091 36096.41183 2.9995931 2000-2009 –> 2010-2019 5.722125 58892.4464 58985.7883
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate 9886.00266 9921.14754 1.6338294 2000-2009 –> 2010-2019 3.237041 17170.5853 17239.1000
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate_star 9894.25600 9929.40089 1.6364397 2000-2009 –> 2010-2019 3.245573 17193.0196 17261.5344
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + nitrate + silicate + phosphate_star 9884.17130 9919.31619 1.6332508 2000-2009 –> 2010-2019 3.241920 17181.8261 17250.3408
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + silicate + phosphate + phosphate_star 9892.18414 9927.32902 1.6357840 2000-2009 –> 2010-2019 3.244894 17190.8918 17259.4066
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 18973.17233 19019.53958 1.3298997 2000-2009 –> 2010-2019 2.771638 36449.5380 36541.3727
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 37477.97028 37526.63308 2.7369823 2000-2009 –> 2010-2019 5.271127 67006.9264 67102.8084
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + phosphate + phosphate_star 37542.82648 37584.53745 2.7488562 2000-2009 –> 2010-2019 5.288131 67095.1949 67177.3795
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + silicate + phosphate 37512.46292 37554.17389 2.7434572 2000-2009 –> 2010-2019 5.279461 67048.6327 67130.8173
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + nitrate + silicate + phosphate 19062.99102 19102.73438 1.3409203 2000-2009 –> 2010-2019 2.803829 36679.9826 36758.6981
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 19064.81548 19111.18273 1.3408992 2000-2009 –> 2010-2019 2.797717 36642.9815 36734.8162
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 37592.71590 37641.37870 2.7573932 2000-2009 –> 2010-2019 5.326217 67292.4676 67388.3497
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + nitrate + phosphate 37636.18259 37684.84539 2.7651648 2000-2009 –> 2010-2019 5.339633 67363.5143 67459.3963
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + nitrate + silicate 18062.16389 18108.53114 1.2253451 2000-2009 –> 2010-2019 2.747767 36071.2964 36163.1311
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate 17650.70708 17697.07433 1.1808576 2000-2009 –> 2010-2019 2.555993 34664.3180 34756.1528
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate 37497.18347 37545.84627 2.7403894 2000-2009 –> 2010-2019 5.257681 66942.2914 67038.1734
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 17669.77417 17716.14142 1.1828830 2000-2009 –> 2010-2019 2.559914 34696.8636 34788.6984
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 37599.05275 37647.71555 2.7585249 2000-2009 –> 2010-2019 5.295813 67143.5841 67239.4662
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 18515.10142 18561.46867 1.2762580 2000-2009 –> 2010-2019 2.715237 35972.7211 36064.5558
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 17628.41448 17674.78173 1.1784940 2000-2009 –> 2010-2019 2.565844 34728.5460 34820.3807
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 37675.60031 37724.26311 2.7722313 2000-2009 –> 2010-2019 5.326797 67305.4115 67401.2935
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 17666.22433 17712.59158 1.1825057 2000-2009 –> 2010-2019 2.559269 34691.4103 34783.2451
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 37582.77414 37631.43694 2.7556188 2000-2009 –> 2010-2019 5.288878 67107.3356 67203.2176
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate_star 17843.42659 17883.16994 1.2017061 2000-2009 –> 2010-2019 2.592175 34963.5254 35042.2409
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate_star 37705.83847 37747.54945 2.7780242 2000-2009 –> 2010-2019 5.340583 67372.9053 67455.0898
obs Indo-Pacific (28.2,28.3] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 11304.12994 11339.43514 2.0291959 2000-2009 –> 2010-2019 4.639259 19554.6057 19622.6531
obs Indo-Pacific (28.2,28.3] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 11223.58334 11264.77274 1.9978950 2000-2009 –> 2010-2019 4.605975 19473.4272 19552.8158
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 11246.10956 11287.29896 2.0063885 2000-2009 –> 2010-2019 4.613951 19495.2654 19574.6540
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 31795.30928 31848.63577 0.6963120 2000-2009 –> 2010-2019 1.350587 54011.2726 54115.8398
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 34922.82781 34976.15430 0.7726396 2000-2009 –> 2010-2019 1.541951 60753.9046 60858.4718
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + phosphate + phosphate_star 11377.90993 11413.21513 2.0575874 2000-2009 –> 2010-2019 4.666786 19627.2380 19695.2854
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 11368.02511 11409.21450 2.0529873 2000-2009 –> 2010-2019 4.662050 19619.1737 19698.5623
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + nitrate + silicate + phosphate 32249.21193 32294.92035 0.7069502 2000-2009 –> 2010-2019 1.453831 57417.8298 57507.4588
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 32071.05845 32124.38494 0.7027272 2000-2009 –> 2010-2019 1.362640 54478.5079 54583.0750
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + nitrate + phosphate 11228.54514 11269.73454 1.9997627 2000-2009 –> 2010-2019 4.657936 19544.2898 19623.6784
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 11360.99065 11402.18005 2.0502694 2000-2009 –> 2010-2019 4.658502 19611.0369 19690.4256
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + silicate + phosphate 34854.15961 34907.48609 0.7708771 2000-2009 –> 2010-2019 1.601779 62404.2444 62508.8115
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 34866.34565 34919.67213 0.7711896 2000-2009 –> 2010-2019 1.602394 62424.5562 62529.1234
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 32060.91276 32114.23925 0.7024901 2000-2009 –> 2010-2019 1.441528 56995.9263 57100.4934
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 34863.97277 34917.29926 0.7711288 2000-2009 –> 2010-2019 1.602296 62421.1608 62525.7279
obs Indo-Pacific (28.2,28.3] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 11260.32019 11301.50959 2.0117652 2000-2009 –> 2010-2019 4.621753 19512.6959 19592.0846
mod Indo-Pacific (28.2,28.3] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 36017.79678 36071.12327 0.8012950 2000-2009 –> 2010-2019 1.638788 63744.2247 63848.7918
obs Indo-Pacific (28.2,28.3] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 11231.63142 11272.82082 2.0009254 2000-2009 –> 2010-2019 4.648890 19534.0471 19613.4358
obs Indo-Pacific (28.2,28.3] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 11374.28330 11415.47270 2.0554083 2000-2009 –> 2010-2019 4.664958 19626.0777 19705.4664
mod Indo-Pacific (28.2,28.3] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 36017.28772 36070.61421 0.8012814 2000-2009 –> 2010-2019 1.640024 63777.0047 63881.5718
mod Indo-Pacific (28.3,28.4] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 20732.98334 20780.18695 1.2630843 2000-2009 –> 2010-2019 2.190028 33637.8329 33730.3697
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 20404.13899 20451.34259 1.2303870 2000-2009 –> 2010-2019 2.159151 33327.8100 33420.3468
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 4497.99867 4532.68782 2.0510142 2000-2009 –> 2010-2019 4.699257 7775.8255 7842.1899
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 4481.75059 4516.43973 2.0351914 2000-2009 –> 2010-2019 4.669696 7752.4829 7818.8472
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + silicate + phosphate_star 4501.96942 4531.70297 2.0568596 2000-2009 –> 2010-2019 4.707418 7778.9881 7835.8718
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 4529.50761 4564.19676 2.0820500 2000-2009 –> 2010-2019 4.715488 7799.6874 7866.0517
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + silicate + phosphate + phosphate_star 4530.76652 4560.50008 2.0852867 2000-2009 –> 2010-2019 4.720196 7799.7084 7856.5921
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + nitrate + phosphate 20381.41806 20428.62167 1.2281593 2000-2009 –> 2010-2019 2.154211 33277.0178 33369.5546
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 20463.84970 20511.05331 1.2362605 2000-2009 –> 2010-2019 2.167276 33410.7620 33503.2988
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + silicate + phosphate 4506.69057 4541.37972 2.0595291 2000-2009 –> 2010-2019 4.693057 7776.9171 7843.2814
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4480.74270 4515.43185 2.0342139 2000-2009 –> 2010-2019 4.683926 7759.3260 7825.6904
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 20379.32375 20426.52736 1.2279542 2000-2009 –> 2010-2019 2.155075 33285.9989 33378.5356
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 20354.22921 20401.43281 1.2254989 2000-2009 –> 2010-2019 2.191748 33657.6627 33750.1995
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 4506.08357 4540.77272 2.0589333 2000-2009 –> 2010-2019 4.708934 7784.8155 7851.1799
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4500.45825 4535.14740 2.0534202 2000-2009 –> 2010-2019 4.687783 7771.1172 7837.4815
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + aou + nitrate + phosphate 20713.10898 20753.56921 1.2612849 2000-2009 –> 2010-2019 2.187877 33612.3137 33691.6309
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 20714.75200 20761.95560 1.2612490 2000-2009 –> 2010-2019 2.185152 33588.0614 33680.5982
obs Indo-Pacific (28.3,28.4] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 4478.35531 4513.04446 2.0319004 2000-2009 –> 2010-2019 4.717296 7775.1849 7841.5492
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 20715.86893 20756.32916 1.2615625 2000-2009 –> 2010-2019 2.189281 33626.7346 33706.0518
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 20717.46918 20764.67279 1.2615223 2000-2009 –> 2010-2019 2.186611 33603.0894 33695.6262
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
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
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
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
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
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
2cbe18c jens-daniel-mueller 2021-06-10
69c79d0 jens-daniel-mueller 2021-06-08
1772903 jens-daniel-mueller 2021-06-07
594ed9a jens-daniel-mueller 2021-06-04
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
ac1a836 jens-daniel-mueller 2021-02-24
b03fbd3 jens-daniel-mueller 2021-02-24
3d3b4cc jens-daniel-mueller 2021-02-23
7b672f7 jens-daniel-mueller 2021-01-11
33ba23c jens-daniel-mueller 2021-01-07
318609d jens-daniel-mueller 2020-12-23
9d0b2d0 jens-daniel-mueller 2020-12-23
0aa2b50 jens-daniel-mueller 2020-12-23
2886da0 jens-daniel-mueller 2020-12-19
02f0ee9 jens-daniel-mueller 2020-12-18
7bcb4eb jens-daniel-mueller 2020-12-18
7131186 jens-daniel-mueller 2020-12-17
a84ff3c jens-daniel-mueller 2020-12-17 
rm(max_rmse)

8.3 Predictor counts

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

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

# print table
lm_all_stats %>%
  gt(rowname_col = "gamma_slab",
     groupname_col = c("basin", "data_source")) %>% 
  summary_rows(
    groups = TRUE,
    fns = list(total = "sum")
  )
aou nitrate phosphate phosphate_star sal silicate temp
Atlantic - mod
(-Inf,26] 6 6 6 8 10 8 6
(26,26.5] 9 7 7 7 7 6 7
(26.5,26.75] 6 7 7 7 6 6 10
(26.75,27] 8 5 8 6 5 10 6
(27,27.25] 9 7 6 4 6 10 6
(27.25,27.5] 6 8 4 7 9 8 7
(27.5,27.75] 7 5 8 5 10 7 6
(27.75,27.85] 5 4 4 6 10 10 6
(27.85,27.95] 7 4 4 5 10 9 7
(27.95,28.05] 4 2 4 8 10 5 10
(28.05,28.1] 6 2 6 6 8 5 10
(28.1,28.15] 5 2 6 7 8 6 10
(28.15,28.2] 6 10 2 6 5 5 9
(28.2, Inf] 9 10 3 4 10 3 5
total 93.00 79.00 75.00 86.00 114.00 98.00 105.00
Indo-Pacific - mod
(-Inf,26] 8 5 7 7 10 5 7
(26,26.5] 8 6 8 7 6 8 7
(26.5,26.75] 7 5 7 7 4 10 8
(26.75,27] 6 10 6 8 5 5 8
(27,27.25] 7 10 7 7 4 5 8
(27.25,27.5] 6 10 6 6 5 6 9
(27.5,27.75] 5 5 6 6 10 10 6
(27.75,27.85] 7 6 7 6 7 6 10
(27.85,27.95] 8 6 5 5 10 5 9
(27.95,28.05] 6 3 4 6 10 6 10
(28.05,28.1] 5 3 4 7 6 10 10
(28.1,28.2] 4 6 6 5 10 10 7
(28.2,28.3] 5 7 8 5 8 10 6
(28.3,28.4] 6 10 9 6 5 4 8
total 88.00 92.00 90.00 88.00 100.00 100.00 113.00
Atlantic - obs
(-Inf,26] 8 9 6 6 7 5 8
(26,26.5] 8 10 8 7 4 4 7
(26.5,26.75] 7 9 7 6 4 4 10
(26.75,27] 4 5 5 8 5 10 10
(27,27.25] 6 8 5 8 6 6 10
(27.25,27.5] 7 4 7 7 10 3 9
(27.5,27.75] 7 3 9 7 4 10 7
(27.75,27.85] 7 5 10 5 5 10 5
(27.85,27.95] 5 4 7 6 10 10 5
(27.95,28.05] 5 4 7 6 10 10 5
(28.05,28.1] 4 7 5 9 10 7 6
(28.1,28.15] 4 7 5 10 6 5 10
(28.15,28.2] 5 5 5 8 4 10 10
(28.2, Inf] 8 4 9 6 10 6 5
total 85.00 84.00 95.00 99.00 95.00 100.00 107.00
Indo-Pacific - obs
(-Inf,26] 8 6 6 6 10 5 8
(26,26.5] 5 5 5 7 10 5 10
(26.5,26.75] 6 7 7 6 8 4 9
(26.75,27] 7 8 7 6 3 6 8
(27,27.25] 6 4 6 7 10 10 4
(27.25,27.5] 6 5 8 7 9 8 4
(27.5,27.75] 7 4 7 6 7 9 7
(27.75,27.85] 5 5 5 8 4 10 10
(27.85,27.95] 8 4 7 7 4 10 8
(27.95,28.05] 9 6 9 8 3 6 4
(28.05,28.1] 8 5 8 6 6 7 6
(28.1,28.2] 6 4 7 6 10 8 6
(28.2,28.3] 10 6 10 8 5 4 5
(28.3,28.4] 6 3 7 8 9 10 5
(28.4, Inf] 3 4 6 6 4 7 10
total 100.00 76.00 105.00 102.00 102.00 109.00 104.00

8.4 RMSE alternatives

8.4.1 AIC

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

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

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