Last updated: 2021-06-04

Checks: 7 0

Knit directory: emlr_obs_v_XXX/

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


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

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

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

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

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

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

Great! You are using Git for version control. Tracking code development and connecting the code version to the results is critical for reproducibility.

The results in this page were generated with repository version db7df0e. 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:  Cant_eMLR.Rproj
    Untracked:  README.md
    Untracked:  _workflowr.yml
    Untracked:  analysis/_site.yml
    Untracked:  code/
    Untracked:  data/

Unstaged changes:
    Modified:   analysis/analysis_anomalous_changes.Rmd

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 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 112.29226 125.81504 0.6342647 2000-2009 –> 2010-2019 5.153268 793.7007 826.3769
obs Atlantic (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 850.02112 869.87509 6.6765781 2000-2009 –> 2010-2019 14.634309 1276.2049 1310.6016
mod Atlantic (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 84.13195 97.65473 0.4812490 2000-2009 –> 2010-2019 4.931078 762.0233 794.6995
obs Atlantic (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 842.59746 862.45143 6.4827617 2000-2009 –> 2010-2019 15.394691 1282.1443 1316.5410
mod Atlantic (-Inf,26] cstar_tref ~ sal + nitrate + phosphate + phosphate_star 118.30580 129.89675 0.6861049 2000-2009 –> 2010-2019 5.209118 797.9165 825.9246
mod Atlantic (-Inf,26] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 59.45820 72.98098 0.3778455 2000-2009 –> 2010-2019 4.789329 735.3763 768.0525
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate 852.24214 869.25983 6.7893532 2000-2009 –> 2010-2019 14.842674 1277.8349 1307.3178
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 111.84622 125.36900 0.6314971 2000-2009 –> 2010-2019 5.361170 803.6434 836.3196
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 850.62668 870.48065 6.6926412 2000-2009 –> 2010-2019 14.682012 1277.2786 1311.6754
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 850.67393 870.52790 6.6938962 2000-2009 –> 2010-2019 14.680089 1277.2789 1311.6757
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 56.06135 69.58413 0.3654695 2000-2009 –> 2010-2019 5.165715 751.2354 783.9116
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 844.97735 864.83132 6.5442751 2000-2009 –> 2010-2019 14.530314 1271.5801 1305.9768
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 82.47591 95.99869 0.4734987 2000-2009 –> 2010-2019 4.725077 749.9761 782.6522
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 111.74389 125.26667 0.6308639 2000-2009 –> 2010-2019 5.119055 791.5924 824.2686
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 850.88683 870.74080 6.6995538 2000-2009 –> 2010-2019 14.709327 1277.8397 1312.2365
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 55.96190 69.48468 0.3651134 2000-2009 –> 2010-2019 5.337238 759.1571 791.8332
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 82.29310 95.81588 0.4726508 2000-2009 –> 2010-2019 5.359487 781.5434 814.2195
obs Atlantic (-Inf,26] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 857.62009 877.47406 6.8809742 2000-2009 –> 2010-2019 15.372329 1291.4626 1325.8593
obs Atlantic (-Inf,26] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 856.70107 876.55504 6.8559255 2000-2009 –> 2010-2019 15.339116 1290.4300 1324.8268
obs Atlantic (-Inf,26] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 853.88419 873.73817 6.7797162 2000-2009 –> 2010-2019 15.226797 1287.1098 1321.5065
obs Atlantic (26,26.5] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 3998.03734 4024.58655 6.1180143 2000-2009 –> 2010-2019 11.587377 7252.8619 7304.9341
mod Atlantic (26,26.5] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 2119.48172 2146.97021 4.0083129 2000-2009 –> 2010-2019 8.425221 3475.3046 3526.8900
obs Atlantic (26,26.5] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 4000.03103 4031.00511 6.1179831 2000-2009 –> 2010-2019 11.564615 7252.5244 7313.2753
mod Atlantic (26,26.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 2109.81869 2137.30717 3.9570009 2000-2009 –> 2010-2019 7.937653 3417.5962 3469.1816
obs Atlantic (26,26.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 3866.60012 3897.57420 5.4909633 2000-2009 –> 2010-2019 10.948724 7121.2163 7181.9672
mod Atlantic (26,26.5] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2075.36390 2102.85238 3.7793297 2000-2009 –> 2010-2019 8.063200 3417.0574 3468.6428
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 2112.95021 2140.43869 3.9735573 2000-2009 –> 2010-2019 8.178543 3446.0569 3497.6423
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 3883.83709 3914.81117 5.5682013 2000-2009 –> 2010-2019 11.094882 7151.5040 7212.2549
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 2113.23972 2140.72820 3.9750915 2000-2009 –> 2010-2019 8.236472 3452.5014 3504.0868
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 3884.78886 3915.76294 5.5724976 2000-2009 –> 2010-2019 11.173023 7166.2597 7227.0106
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 2079.14826 2106.63675 3.7984477 2000-2009 –> 2010-2019 8.387477 3452.6327 3504.2181
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 2114.10194 2141.59042 3.9796640 2000-2009 –> 2010-2019 8.304232 3460.1638 3511.7492
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 3895.03534 3926.00942 5.6189613 2000-2009 –> 2010-2019 11.280473 7187.7700 7248.5209
obs Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + phosphate 3996.58675 4023.13596 6.1108267 2000-2009 –> 2010-2019 11.671992 7268.7227 7320.7949
mod Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 2117.53298 2145.02146 3.9979115 2000-2009 –> 2010-2019 8.435856 3475.5509 3527.1364
obs Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 3998.58589 4029.55997 6.1108225 2000-2009 –> 2010-2019 11.593396 7257.9195 7318.6704
mod Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2113.75516 2141.24364 3.9778243 2000-2009 –> 2010-2019 8.440997 3474.3920 3525.9774
obs Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 3998.51880 4029.49289 6.1104903 2000-2009 –> 2010-2019 11.635898 7265.9463 7326.6972
mod Atlantic (26,26.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 2115.93499 2143.42347 3.9894024 2000-2009 –> 2010-2019 8.479962 3479.3982 3530.9836
obs Atlantic (26,26.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 3999.05791 4030.03199 6.1131604 2000-2009 –> 2010-2019 11.693586 7276.7896 7337.5405
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 2751.64574 2780.90585 4.1167469 2000-2009 –> 2010-2019 7.521296 4706.3791 4762.9768
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 4796.38740 4828.98448 5.2308575 2000-2009 –> 2010-2019 10.902204 8065.6922 8128.0120
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 2754.71307 2783.97318 4.1298395 2000-2009 –> 2010-2019 7.413858 4682.9897 4739.5874
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 4795.49514 4828.09222 5.2278588 2000-2009 –> 2010-2019 10.928384 8070.0960 8132.4158
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 2802.06719 2831.32730 4.3373319 2000-2009 –> 2010-2019 7.372535 4672.5123 4729.1099
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 4889.77404 4922.37112 5.5544107 2000-2009 –> 2010-2019 11.252280 8163.8938 8226.2137
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2808.25973 2837.51985 4.3652257 2000-2009 –> 2010-2019 7.308854 4656.2186 4712.8163
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 2761.34070 2790.60081 4.1582714 2000-2009 –> 2010-2019 7.283549 4653.2514 4709.8491
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 4792.60438 4825.20146 5.2181554 2000-2009 –> 2010-2019 10.917326 8066.9599 8129.2798
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 2830.50668 2859.76679 4.4669232 2000-2009 –> 2010-2019 7.503786 4701.3530 4757.9507
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 2809.24186 2838.50197 4.3696660 2000-2009 –> 2010-2019 7.298476 4653.4964 4710.0941
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + phosphate 4857.34142 4885.28178 5.4468317 2000-2009 –> 2010-2019 11.122333 8125.4020 8178.8190
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + phosphate_star 4854.39016 4882.33051 5.4365105 2000-2009 –> 2010-2019 11.137417 8127.0600 8180.4770
mod Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 2773.80767 2803.06778 4.2122849 2000-2009 –> 2010-2019 7.484564 4699.4558 4756.0534
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 4857.50477 4890.10186 5.4404063 2000-2009 –> 2010-2019 11.115186 8127.4340 8189.7538
mod Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2777.66870 2806.92882 4.2291548 2000-2009 –> 2010-2019 7.379878 4675.5315 4732.1292
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 4854.70576 4887.30284 5.4306286 2000-2009 –> 2010-2019 11.123718 8127.9595 8190.2794
obs Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 4851.96416 4879.90452 5.4280409 2000-2009 –> 2010-2019 11.128265 8124.5103 8177.9273
mod Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 2781.45978 2810.71990 4.2457847 2000-2009 –> 2010-2019 7.341943 4666.4996 4723.0972
obs Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 4852.25205 4884.84914 5.4220716 2000-2009 –> 2010-2019 11.115074 8125.4899 8187.8098
mod Atlantic (26.75,27] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 4277.07275 4310.79756 2.4921808 2000-2009 –> 2010-2019 5.000943 7219.1601 7283.9603
mod Atlantic (26.75,27] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 4408.85742 4442.58223 2.6784826 2000-2009 –> 2010-2019 5.214940 7364.6905 7429.4908
mod Atlantic (26.75,27] cstar_tref ~ sal + aou + silicate + phosphate 4561.11696 4590.02394 2.9143228 2000-2009 –> 2010-2019 5.450786 7514.9530 7570.4961
mod Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate 4546.12606 4579.85087 2.8873603 2000-2009 –> 2010-2019 5.208852 7391.0841 7455.8844
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate 5687.58692 5721.76399 4.4397285 2000-2009 –> 2010-2019 10.818177 9712.8916 9777.9972
mod Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4565.38721 4599.11203 2.9179445 2000-2009 –> 2010-2019 5.266662 7424.9431 7489.7433
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5692.37611 5726.55317 4.4506458 2000-2009 –> 2010-2019 10.772964 9706.8442 9771.9498
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 5702.39453 5736.57159 4.4735705 2000-2009 –> 2010-2019 10.782180 9714.2012 9779.3069
mod Atlantic (26.75,27] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4572.47149 4606.19630 2.9292747 2000-2009 –> 2010-2019 5.295379 7441.2611 7506.0614
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5693.41408 5727.59114 4.4530155 2000-2009 –> 2010-2019 10.772208 9707.2757 9772.3813
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + silicate + phosphate_star 5711.91022 5741.20485 4.5000673 2000-2009 –> 2010-2019 10.905839 9740.4554 9796.2602
mod Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 4294.21475 4327.93956 2.5156611 2000-2009 –> 2010-2019 4.894133 7169.5315 7234.3317
obs Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 5693.10635 5727.28341 4.4523128 2000-2009 –> 2010-2019 10.921289 9735.6889 9800.7945
mod Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 4303.46710 4337.19192 2.5284263 2000-2009 –> 2010-2019 4.897186 7173.6610 7238.4613
obs Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5698.83918 5733.01624 4.4654215 2000-2009 –> 2010-2019 10.846921 9724.7301 9789.8357
mod Atlantic (26.75,27] cstar_tref ~ temp + aou + silicate + phosphate 4691.53961 4720.44660 3.1298487 2000-2009 –> 2010-2019 5.569642 7596.7259 7652.2690
mod Atlantic (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 4324.96342 4358.68823 2.5583347 2000-2009 –> 2010-2019 4.930600 7197.0089 7261.8091
obs Atlantic (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 5700.96956 5735.14663 4.4703027 2000-2009 –> 2010-2019 10.863355 9729.0780 9794.1836
obs Atlantic (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate_star 5705.73104 5735.02566 4.4858300 2000-2009 –> 2010-2019 10.909744 9737.7435 9793.5483
obs Atlantic (26.75,27] cstar_tref ~ temp + silicate + phosphate + phosphate_star 5701.13118 5730.42581 4.4752608 2000-2009 –> 2010-2019 10.926679 9738.3817 9794.1865
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 5596.17651 5632.15968 2.2093544 2000-2009 –> 2010-2019 4.784496 9767.5704 9837.0051
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 5677.14374 5713.12691 2.2813771 2000-2009 –> 2010-2019 5.008754 9949.5097 10018.9444
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate 5794.50538 5825.34810 2.3918567 2000-2009 –> 2010-2019 5.020402 9999.9841 10059.4996
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 5787.50158 5823.48475 2.3833395 2000-2009 –> 2010-2019 4.921896 9933.7408 10003.1754
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate 11084.91399 11123.02373 6.1605793 2000-2009 –> 2010-2019 12.057977 19056.9907 19130.9943
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 11059.06345 11097.17319 6.1141893 2000-2009 –> 2010-2019 11.897560 18982.4852 19056.4888
mod Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + silicate 5781.28854 5817.27171 2.3774799 2000-2009 –> 2010-2019 4.989611 9977.7551 10047.1898
mod Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate 5786.95409 5822.93726 2.3828226 2000-2009 –> 2010-2019 4.933042 9941.2517 10010.6863
mod Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5807.80221 5843.78538 2.4025861 2000-2009 –> 2010-2019 4.974434 9976.9460 10046.3807
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 11104.39075 11142.50049 6.1957637 2000-2009 –> 2010-2019 11.934850 19008.6570 19082.6606
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 11065.97065 11104.08039 6.1265503 2000-2009 –> 2010-2019 11.942998 19003.6043 19077.6079
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 11123.05373 11161.16347 6.2296665 2000-2009 –> 2010-2019 12.064337 19068.4826 19142.4862
mod Atlantic (27,27.25] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5816.09449 5852.07766 2.4104925 2000-2009 –> 2010-2019 4.990838 9991.0378 10060.4725
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 11106.05636 11144.16610 6.1987819 2000-2009 –> 2010-2019 11.951157 19016.0863 19090.0899
obs Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + phosphate_star 11079.90152 11112.56701 6.1551552 2000-2009 –> 2010-2019 12.066285 19055.7741 19119.2057
obs Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 11045.99071 11084.10045 6.0908628 2000-2009 –> 2010-2019 12.065784 19050.6121 19124.6157
mod Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5917.70267 5953.68584 2.5095109 2000-2009 –> 2010-2019 5.090338 10092.9736 10162.4083
obs Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 11005.11332 11043.22306 6.0184954 2000-2009 –> 2010-2019 11.861171 18953.9587 19027.9623
mod Atlantic (27,27.25] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 5781.28455 5817.26772 2.3774762 2000-2009 –> 2010-2019 4.951916 9952.1993 10021.6339
obs Atlantic (27,27.25] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 11016.32500 11054.43474 6.0382580 2000-2009 –> 2010-2019 11.915142 18979.7182 19053.7218
obs Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 13193.32810 13232.77523 5.8385008 2000-2009 –> 2010-2019 10.519306 21887.1858 21963.6602
mod Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate 11281.53610 11315.28674 3.7850923 2000-2009 –> 2010-2019 7.121819 20960.0895 21026.9518
mod Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 11109.87218 11149.24793 3.6280400 2000-2009 –> 2010-2019 6.863480 20676.8655 20754.8715
mod Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 11130.98424 11170.35999 3.6467791 2000-2009 –> 2010-2019 6.968727 20795.1853 20873.1913
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 13178.22233 13217.66946 5.8172365 2000-2009 –> 2010-2019 10.419445 21822.4638 21898.9382
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 11352.00621 11391.38196 3.8488656 2000-2009 –> 2010-2019 6.976101 20793.6858 20871.6918
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 13176.76737 13216.21450 5.8151924 2000-2009 –> 2010-2019 10.392393 21805.0440 21881.5184
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + silicate 11260.84060 11300.21635 3.7641877 2000-2009 –> 2010-2019 6.952605 20773.8990 20851.9050
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + phosphate 13244.44225 13278.25407 5.9138884 2000-2009 –> 2010-2019 10.530748 21895.9965 21961.5460
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + phosphate_star 13247.83591 13281.64774 5.9187382 2000-2009 –> 2010-2019 10.512960 21884.9878 21950.5373
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 13241.38475 13280.83187 5.9066683 2000-2009 –> 2010-2019 10.521577 21893.7003 21970.1747
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 11077.17048 11116.54623 3.5992037 2000-2009 –> 2010-2019 7.114699 20949.9931 21027.9991
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 13244.57132 13284.01845 5.9112165 2000-2009 –> 2010-2019 10.505329 21883.6540 21960.1284
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 11352.05360 11391.42935 3.8489101 2000-2009 –> 2010-2019 6.903681 20707.3634 20785.3694
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 13173.73628 13213.18341 5.8109364 2000-2009 –> 2010-2019 10.393270 21805.2966 21881.7710
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 11062.12180 11101.49755 3.5860110 2000-2009 –> 2010-2019 6.980510 20805.9147 20883.9207
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + phosphate + phosphate_star 13244.60065 13278.41247 5.9141147 2000-2009 –> 2010-2019 10.512152 21884.1850 21949.7345
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 11077.96081 11117.33656 3.5998979 2000-2009 –> 2010-2019 7.111980 20947.2044 21025.2104
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 13241.52728 13280.97441 5.9068717 2000-2009 –> 2010-2019 10.504432 21882.8076 21959.2820
mod Atlantic (27.25,27.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 11431.88681 11471.26256 3.9246259 2000-2009 –> 2010-2019 7.000374 20812.4073 20890.4133
obs Atlantic (27.5,27.75] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 14151.73814 14191.98340 5.0936150 2000-2009 –> 2010-2019 10.364399 25276.2482 25354.9739
obs Atlantic (27.5,27.75] cstar_tref ~ aou + silicate + phosphate + phosphate_star 14161.71092 14196.20685 5.1067752 2000-2009 –> 2010-2019 10.378703 25285.0034 25352.4826
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + nitrate + phosphate 18450.52123 18488.50641 2.2312680 2000-2009 –> 2010-2019 4.664316 34069.7820 34144.5228
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 17999.09591 18043.41196 2.1126445 2000-2009 –> 2010-2019 4.538851 33601.3140 33688.5116
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate 18658.51556 18696.50075 2.2878891 2000-2009 –> 2010-2019 4.732929 34311.0206 34385.7613
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 18629.29507 18673.61112 2.2792994 2000-2009 –> 2010-2019 4.718648 34268.0441 34355.2416
obs Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 14163.38695 14203.63220 5.1064187 2000-2009 –> 2010-2019 10.302555 25236.4619 25315.1876
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 18347.58096 18391.89701 2.2032348 2000-2009 –> 2010-2019 4.623351 33932.8049 34020.0024
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 18656.60205 18700.91809 2.2868106 2000-2009 –> 2010-2019 4.714101 34261.8405 34349.0380
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 14169.03180 14209.27705 5.1126347 2000-2009 –> 2010-2019 10.316131 25247.2108 25325.9365
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 18498.44218 18542.75822 2.2436471 2000-2009 –> 2010-2019 4.624860 33974.0835 34061.2810
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 14197.82632 14238.07158 5.1444609 2000-2009 –> 2010-2019 10.375576 25295.0947 25373.8204
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 18662.37815 18706.69420 2.2884026 2000-2009 –> 2010-2019 4.677892 34161.4802 34248.6777
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 18597.88826 18642.20430 2.2706909 2000-2009 –> 2010-2019 4.662151 34102.5666 34189.7641
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 18501.61749 18545.93353 2.2445056 2000-2009 –> 2010-2019 4.627907 33983.4709 34070.6685
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 14177.97164 14218.21690 5.1224947 2000-2009 –> 2010-2019 10.338578 25264.8630 25343.5887
obs Atlantic (27.5,27.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 14158.40362 14198.64888 5.1009373 2000-2009 –> 2010-2019 10.388706 25294.5154 25373.2411
obs Atlantic (27.5,27.75] cstar_tref ~ temp + aou + silicate + phosphate 14167.28203 14201.77797 5.1129104 2000-2009 –> 2010-2019 10.401285 25301.8068 25369.2859
obs Atlantic (27.5,27.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 14168.21060 14208.45586 5.1117300 2000-2009 –> 2010-2019 10.423493 25320.6481 25399.3738
obs Atlantic (27.5,27.75] cstar_tref ~ temp + silicate + phosphate + phosphate_star 14176.62572 14211.12166 5.1232168 2000-2009 –> 2010-2019 10.435149 25327.1782 25394.6574
obs Atlantic (27.75,27.85] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 9246.19236 9283.33338 5.3719621 2000-2009 –> 2010-2019 10.549107 16973.3105 17046.4181
obs Atlantic (27.75,27.85] cstar_tref ~ aou + silicate + phosphate + phosphate_star 9244.19323 9276.02839 5.3719637 2000-2009 –> 2010-2019 10.573929 16981.3544 17044.0180
mod Atlantic (27.75,27.85] cstar_tref ~ sal + aou + nitrate + silicate 5223.07992 5255.05055 1.3389101 2000-2009 –> 2010-2019 2.573897 9489.3765 9552.3908
obs Atlantic (27.75,27.85] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 9232.90194 9270.04296 5.3480411 2000-2009 –> 2010-2019 10.541283 16967.8369 17040.9444
mod Atlantic (27.75,27.85] cstar_tref ~ sal + aou + silicate + phosphate 5009.86162 5041.83224 1.2483922 2000-2009 –> 2010-2019 2.426688 9153.5126 9216.5270
obs Atlantic (27.75,27.85] cstar_tref ~ sal + aou + silicate + phosphate 9230.96086 9262.79602 5.3481470 2000-2009 –> 2010-2019 10.569216 16977.3522 17040.0158
obs Atlantic (27.75,27.85] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 9234.29568 9271.43670 5.3505447 2000-2009 –> 2010-2019 10.575167 16984.3999 17057.5075
mod Atlantic (27.75,27.85] cstar_tref ~ sal + nitrate + silicate + phosphate_star 5615.64947 5647.62009 1.5230826 2000-2009 –> 2010-2019 2.885221 10137.7144 10200.7288
mod Atlantic (27.75,27.85] cstar_tref ~ sal + silicate + phosphate + phosphate_star 5607.34615 5639.31677 1.5189363 2000-2009 –> 2010-2019 2.886270 10139.3462 10202.3606
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + silicate 5181.77196 5219.07103 1.3200082 2000-2009 –> 2010-2019 2.552107 9443.9573 9517.4741
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 4895.94927 4933.24833 1.2017785 2000-2009 –> 2010-2019 2.374573 9029.3862 9102.9029
obs Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 9231.86087 9269.00189 5.3461719 2000-2009 –> 2010-2019 10.556392 16975.0144 17048.1219
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5023.97147 5061.27053 1.2533653 2000-2009 –> 2010-2019 2.465168 9242.8062 9316.3230
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 5028.84807 5066.14713 1.2553735 2000-2009 –> 2010-2019 2.466459 9246.1377 9319.6545
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5017.98177 5055.28083 1.2509030 2000-2009 –> 2010-2019 2.458787 9228.3615 9301.8782
obs Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 9233.51531 9270.65633 5.3491428 2000-2009 –> 2010-2019 10.568878 16981.2634 17054.3709
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate_star 5104.02180 5135.99242 1.2875861 2000-2009 –> 2010-2019 2.581531 9492.0365 9555.0509
obs Atlantic (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 9246.24246 9283.38348 5.3720525 2000-2009 –> 2010-2019 10.556134 16976.7323 17049.8398
obs Atlantic (27.75,27.85] cstar_tref ~ temp + aou + silicate + phosphate 9244.24385 9276.07901 5.3720550 2000-2009 –> 2010-2019 10.582290 16985.4047 17048.0683
obs Atlantic (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 9246.59213 9283.73315 5.3726833 2000-2009 –> 2010-2019 10.564608 16980.8884 17053.9960
mod Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate 6335.30776 6367.18708 1.9913857 2000-2009 –> 2010-2019 4.910848 12842.1089 12905.0274
obs Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 19491.48588 19533.87642 5.3162288 2000-2009 –> 2010-2019 10.587009 36162.0044 36245.7097
mod Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 6337.25137 6374.44391 1.9913483 2000-2009 –> 2010-2019 4.833981 12776.5001 12849.9050
obs Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 18889.37946 18931.77001 4.8319609 2000-2009 –> 2010-2019 9.869156 35314.8495 35398.5548
mod Atlantic (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate 6067.00470 6098.88402 1.8210192 2000-2009 –> 2010-2019 4.556299 12403.8543 12466.7727
obs Atlantic (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate_star 18893.80456 18930.13932 4.8368882 2000-2009 –> 2010-2019 9.874358 35317.5692 35389.3167
obs Atlantic (27.85,27.95] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 19289.80928 19332.19983 5.1488445 2000-2009 –> 2010-2019 10.342602 35880.7453 35964.4506
obs Atlantic (27.85,27.95] cstar_tref ~ sal + silicate + phosphate + phosphate_star 19331.52239 19367.85714 5.1846717 2000-2009 –> 2010-2019 10.381953 35924.1255 35995.8730
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + silicate 6176.78798 6213.98052 1.8876339 2000-2009 –> 2010-2019 4.794736 12674.5241 12747.9290
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + phosphate 6526.38326 6558.26259 2.1223470 2000-2009 –> 2010-2019 4.875153 12879.8902 12942.8087
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 5757.25381 5794.44635 1.6412849 2000-2009 –> 2010-2019 4.376268 12095.8200 12169.2249
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 18727.70586 18770.09641 4.7096151 2000-2009 –> 2010-2019 9.716044 35120.0555 35203.7609
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5861.65180 5898.84434 1.6994059 2000-2009 –> 2010-2019 4.459201 12223.7717 12297.1765
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 18708.87253 18751.26307 4.6955660 2000-2009 –> 2010-2019 9.697233 35096.0782 35179.7835
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 18708.42002 18750.81057 4.6952290 2000-2009 –> 2010-2019 9.740572 35142.6272 35226.3326
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 5875.42951 5912.62205 1.7072285 2000-2009 –> 2010-2019 4.511439 12279.1878 12352.5927
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate 18732.46864 18768.80339 4.7146701 2000-2009 –> 2010-2019 9.761388 35166.1483 35237.8958
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5851.27596 5888.46850 1.6935385 2000-2009 –> 2010-2019 4.444789 12205.3082 12278.7131
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 18726.43013 18768.82068 4.7086621 2000-2009 –> 2010-2019 9.713936 35117.5319 35201.2373
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate_star 5954.64421 5986.52353 1.7540771 2000-2009 –> 2010-2019 4.788139 12561.8612 12624.7797
obs Atlantic (27.95,28.05] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 19947.07164 19989.58982 5.3973133 2000-2009 –> 2010-2019 12.104578 37714.5328 37798.2850
obs Atlantic (27.95,28.05] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 19532.13881 19574.65699 5.0595118 2000-2009 –> 2010-2019 11.621630 37182.6856 37266.4379
obs Atlantic (27.95,28.05] cstar_tref ~ sal + aou + silicate + phosphate_star 19618.25855 19654.70271 5.1294366 2000-2009 –> 2010-2019 11.705880 37278.4589 37350.2466
obs Atlantic (27.95,28.05] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 19668.25049 19710.76867 5.1679247 2000-2009 –> 2010-2019 11.774928 37355.2259 37438.9782
obs Atlantic (27.95,28.05] cstar_tref ~ sal + silicate + phosphate + phosphate_star 19796.95276 19833.39692 5.2742148 2000-2009 –> 2010-2019 11.895101 37493.1453 37564.9329
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate 9622.43132 9655.35100 3.5771558 2000-2009 –> 2010-2019 7.370873 17461.9294 17526.4080
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate_star 9652.26384 9685.18352 3.6071902 2000-2009 –> 2010-2019 7.399792 17490.9259 17555.4045
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 9485.98573 9524.39202 3.4410133 2000-2009 –> 2010-2019 7.233640 17326.6665 17401.8916
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 19650.50167 19693.01985 5.1536571 2000-2009 –> 2010-2019 11.752819 37331.1307 37414.8829
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 9497.74698 9536.15327 3.4523747 2000-2009 –> 2010-2019 7.242187 17336.3164 17411.5415
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 19595.91703 19638.43521 5.1100250 2000-2009 –> 2010-2019 11.685927 37257.6771 37341.4294
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate + phosphate_star 9716.85928 9749.77896 3.6730899 2000-2009 –> 2010-2019 7.451946 17545.1947 17609.6733
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 19680.08439 19722.60258 5.1774595 2000-2009 –> 2010-2019 11.892778 37453.9589 37537.7112
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 9471.95391 9510.36020 3.4275074 2000-2009 –> 2010-2019 7.169174 17274.1617 17349.3867
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + phosphate + phosphate_star 9647.46254 9680.38222 3.6023394 2000-2009 –> 2010-2019 7.394938 17486.1219 17550.6005
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + phosphate_star 9730.01669 9757.44975 3.6887269 2000-2009 –> 2010-2019 7.484020 17568.6963 17622.4285
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate 19783.08936 19819.53352 5.2628379 2000-2009 –> 2010-2019 11.982648 37558.5367 37630.3244
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 9497.99478 9536.40108 3.4526144 2000-2009 –> 2010-2019 7.242426 17336.5633 17411.7883
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 19634.73754 19677.25573 5.1410180 2000-2009 –> 2010-2019 11.731502 37308.3344 37392.0866
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate_star 9496.60669 9529.52637 3.4532066 2000-2009 –> 2010-2019 7.243019 17333.1762 17397.6548
obs Atlantic (28.05,28.1] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 14407.92395 14448.21732 5.2695909 2000-2009 –> 2010-2019 11.606991 26001.0578 26079.7141
obs Atlantic (28.05,28.1] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 14352.25288 14392.54625 5.2071716 2000-2009 –> 2010-2019 11.544369 25945.2733 26023.9297
obs Atlantic (28.05,28.1] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 14355.39070 14395.68407 5.2106700 2000-2009 –> 2010-2019 11.535122 25941.2721 26019.9285
obs Atlantic (28.05,28.1] cstar_tref ~ sal + silicate + phosphate + phosphate_star 14441.24722 14475.78439 5.3095833 2000-2009 –> 2010-2019 11.657543 26038.2845 26105.7043
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 14474.53023 14514.82360 5.3452547 2000-2009 –> 2010-2019 11.650452 26049.5996 26128.2560
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate 3562.36572 3592.04181 1.3359076 2000-2009 –> 2010-2019 2.327088 5726.1471 5783.6467
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate_star 3601.10429 3630.78038 1.3610455 2000-2009 –> 2010-2019 2.370657 5793.0016 5850.5012
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 3563.87765 3598.49975 1.3355939 2000-2009 –> 2010-2019 2.308748 5701.6505 5768.7334
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 3600.92724 3635.54934 1.3596203 2000-2009 –> 2010-2019 2.355833 5774.4372 5841.5201
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 14388.75587 14429.04924 5.2480153 2000-2009 –> 2010-2019 11.519437 25944.7787 26023.4351
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 14467.47291 14507.76628 5.3371865 2000-2009 –> 2010-2019 11.645665 26044.3871 26123.0435
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + phosphate_star 3647.49651 3677.17260 1.3917731 2000-2009 –> 2010-2019 2.412706 5856.4112 5913.9108
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + phosphate_star 14476.05487 14510.59204 5.3492888 2000-2009 –> 2010-2019 11.658195 26051.2094 26118.6292
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 3625.08955 3659.71165 1.3755218 2000-2009 –> 2010-2019 2.395743 5834.9398 5902.0227
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 14355.99668 14396.29005 5.2113459 2000-2009 –> 2010-2019 11.502427 25923.1178 26001.7742
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + phosphate + phosphate_star 3590.21036 3619.88645 1.3539289 2000-2009 –> 2010-2019 2.359139 5775.4401 5832.9397
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 3591.22882 3625.85092 1.3532895 2000-2009 –> 2010-2019 2.343063 5754.8428 5821.9257
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 14430.09414 14470.38751 5.2946562 2000-2009 –> 2010-2019 11.592076 26000.7868 26079.4432
mod Atlantic (28.05,28.1] cstar_tref ~ temp + aou + phosphate 3577.36094 3602.09101 1.3468783 2000-2009 –> 2010-2019 2.421644 5862.6894 5910.6058
mod Atlantic (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate 3579.09797 3608.77406 1.3467079 2000-2009 –> 2010-2019 2.397622 5832.1805 5889.6801
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 16565.60044 16606.60551 5.9375321 2000-2009 –> 2010-2019 15.026642 27624.4437 27702.7478
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + phosphate 4580.32385 4612.32591 1.0757334 2000-2009 –> 2010-2019 2.401785 8950.6549 9013.5893
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + phosphate_star 4628.03720 4660.03926 1.0926272 2000-2009 –> 2010-2019 2.426969 9014.3358 9077.2703
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + silicate + phosphate 4547.07786 4584.41359 1.0634218 2000-2009 –> 2010-2019 2.270791 8679.1921 8752.6156
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4604.27546 4641.61120 1.0834731 2000-2009 –> 2010-2019 2.304791 8765.8184 8839.2419
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 16600.03889 16641.04397 5.9771998 2000-2009 –> 2010-2019 15.082561 27664.3234 27742.6275
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 16548.78789 16589.79296 5.9182623 2000-2009 –> 2010-2019 15.001223 27605.5698 27683.8739
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + phosphate_star 4648.96382 4680.96587 1.1001201 2000-2009 –> 2010-2019 2.409019 8985.9389 9048.8733
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 4650.47799 4687.81372 1.0999456 2000-2009 –> 2010-2019 2.378600 8929.5597 9002.9832
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 16574.15941 16615.16448 5.9473660 2000-2009 –> 2010-2019 15.056474 27639.6973 27718.0015
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + phosphate + phosphate_star 4613.47686 4645.47891 1.0874439 2000-2009 –> 2010-2019 2.416132 8988.8980 9051.8324
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + phosphate + phosphate_star 16601.26615 16636.41335 5.9809307 2000-2009 –> 2010-2019 15.084728 27663.0275 27730.1454
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4583.89126 4621.22699 1.0762842 2000-2009 –> 2010-2019 2.286871 8722.8240 8796.2475
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 16599.99154 16640.99661 5.9771451 2000-2009 –> 2010-2019 15.079161 27663.1568 27741.4609
obs Atlantic (28.1,28.15] cstar_tref ~ temp + aou + nitrate + phosphate_star 16592.64465 16627.79185 5.9709690 2000-2009 –> 2010-2019 15.074437 27654.2958 27721.4136
obs Atlantic (28.1,28.15] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 16568.23698 16609.24206 5.9405596 2000-2009 –> 2010-2019 15.043281 27631.6385 27709.9426
mod Atlantic (28.1,28.15] cstar_tref ~ temp + aou + silicate + phosphate 4661.65665 4693.65871 1.1046899 2000-2009 –> 2010-2019 2.390335 8952.7067 9015.6411
obs Atlantic (28.1,28.15] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 16579.57996 16614.72717 5.9559052 2000-2009 –> 2010-2019 15.054134 27639.4777 27706.5955
obs Atlantic (28.1,28.15] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 16556.86242 16597.86750 5.9275092 2000-2009 –> 2010-2019 15.025141 27618.5602 27696.8644
mod Atlantic (28.1,28.15] cstar_tref ~ temp + silicate + phosphate + phosphate_star 4709.25714 4741.25920 1.1219970 2000-2009 –> 2010-2019 2.413151 9011.2609 9074.1954
mod Atlantic (28.15,28.2] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 7255.79282 7297.58317 0.8459183 2000-2009 –> 2010-2019 1.897094 13650.4664 13732.0406
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate 7314.86494 7350.68524 0.8548944 2000-2009 –> 2010-2019 1.897297 13671.1340 13741.0547
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate + silicate 7217.16492 7258.95526 0.8402897 2000-2009 –> 2010-2019 1.878862 13559.4416 13641.0157
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate 5567.39033 5601.09213 5.0989030 2000-2009 –> 2010-2019 12.844751 10191.2973 10256.4976
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5527.37771 5561.07951 4.9881475 2000-2009 –> 2010-2019 12.651902 10137.1136 10202.3140
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + nitrate + phosphate_star 7630.14454 7665.96483 0.9027702 2000-2009 –> 2010-2019 1.963345 14061.4890 14131.4097
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 7592.82172 7634.61207 0.8966556 2000-2009 –> 2010-2019 1.950818 13999.8196 14081.3938
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 5512.98660 5546.68840 4.9489037 2000-2009 –> 2010-2019 12.661785 10131.2211 10196.4214
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5543.77540 5577.47720 5.0332426 2000-2009 –> 2010-2019 12.708683 10155.5377 10220.7380
mod Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + phosphate 7320.76901 7356.58930 0.8557672 2000-2009 –> 2010-2019 1.897285 13673.3466 13743.2673
mod Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + phosphate_star 7354.22975 7390.05004 0.8607305 2000-2009 –> 2010-2019 1.907796 13729.8868 13799.8075
mod Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate 7245.04408 7280.86438 0.8446402 2000-2009 –> 2010-2019 1.895449 13636.2015 13706.1222
obs Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 5528.14467 5561.84647 4.9902477 2000-2009 –> 2010-2019 12.994919 10195.7665 10260.9668
obs Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5499.75341 5533.45521 4.9130900 2000-2009 –> 2010-2019 12.819071 10150.8757 10216.0760
obs Atlantic (28.15,28.2] cstar_tref ~ temp + aou + silicate + phosphate_star 5527.90676 5556.79402 4.9950762 2000-2009 –> 2010-2019 12.922162 10180.5748 10236.4607
mod Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 7405.37130 7441.19160 0.8683721 2000-2009 –> 2010-2019 1.919415 13797.4951 13867.4158
obs Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 5506.20818 5539.90998 4.9305264 2000-2009 –> 2010-2019 12.865910 10162.2675 10227.4678
mod Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + silicate + phosphate_star 7610.10760 7645.92790 0.8996493 2000-2009 –> 2010-2019 1.953957 14015.7075 14085.6282
obs Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + silicate + phosphate_star 5513.93926 5542.82652 4.9569302 2000-2009 –> 2010-2019 12.983107 10183.1296 10239.0156
obs Atlantic (28.15,28.2] cstar_tref ~ temp + silicate + phosphate + phosphate_star 5543.66553 5572.55279 5.0384668 2000-2009 –> 2010-2019 12.986793 10199.8924 10255.7784
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate 9976.68015 10009.07241 0.6818803 2000-2009 –> 2010-2019 1.432832 15819.1804 15880.8407
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + phosphate 9911.28660 9950.15732 0.6771201 2000-2009 –> 2010-2019 1.425834 15740.4170 15814.4094
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 10756.61904 10796.11335 3.1850413 2000-2009 –> 2010-2019 6.499563 18334.4772 18410.9025
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + phosphate_star 9830.13305 9869.00376 0.6714320 2000-2009 –> 2010-2019 1.420186 15659.5420 15733.5344
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate 9868.17724 9907.04796 0.6740925 2000-2009 –> 2010-2019 1.408261 15596.2718 15670.2641
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 10755.19494 10794.68925 3.1839533 2000-2009 –> 2010-2019 6.500744 18335.0308 18411.4561
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 9231.00629 9276.35546 0.6307600 2000-2009 –> 2010-2019 1.352751 14874.9648 14961.2892
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + phosphate + phosphate_star 10757.95854 10791.81081 3.1875943 2000-2009 –> 2010-2019 6.503531 18335.0497 18400.5572
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + silicate + phosphate 10753.79212 10787.64438 3.1844095 2000-2009 –> 2010-2019 6.501512 18331.8990 18397.4064
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 10735.74732 10775.24163 3.1691318 2000-2009 –> 2010-2019 6.482937 18312.9800 18389.4054
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate 9900.84642 9939.71714 0.6763856 2000-2009 –> 2010-2019 1.425810 15734.8605 15808.8528
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate 9833.76949 9879.11866 0.6715462 2000-2009 –> 2010-2019 1.420042 15663.4033 15749.7277
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate 10754.86238 10794.35669 3.1836992 2000-2009 –> 2010-2019 6.502173 18336.1645 18412.5898
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 9828.37701 9873.72618 0.6711699 2000-2009 –> 2010-2019 1.419890 15659.5494 15745.8738
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + silicate 9794.11265 9839.46182 0.6687836 2000-2009 –> 2010-2019 1.400542 15507.2772 15593.6016
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + silicate + phosphate 10726.20824 10765.70255 3.1618870 2000-2009 –> 2010-2019 6.478974 18306.3018 18382.7272
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 10769.42383 10808.91814 3.1948414 2000-2009 –> 2010-2019 6.469304 18312.1402 18388.5656
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 9875.03474 9920.38390 0.6744330 2000-2009 –> 2010-2019 1.426039 15726.0216 15812.3460
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 10796.97655 10836.47086 3.2160309 2000-2009 –> 2010-2019 6.497808 18346.1417 18422.5671
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 10742.23745 10781.73177 3.1740704 2000-2009 –> 2010-2019 6.455869 18291.4212 18367.8466
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 22625.96829 22669.33563 5.4784065 2000-2009 –> 2010-2019 10.986523 38109.5164 38193.5817
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 29326.05469 29371.20907 5.5512888 2000-2009 –> 2010-2019 11.849614 50316.4090 50404.0990
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + phosphate + phosphate_star 22936.47629 22973.64829 5.7197821 2000-2009 –> 2010-2019 11.567971 38714.5755 38786.6315
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 22719.14138 22762.50872 5.5492860 2000-2009 –> 2010-2019 11.168521 38301.5552 38385.6205
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 29645.24600 29690.40039 5.7439454 2000-2009 –> 2010-2019 12.266054 50860.3071 50947.9971
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate 30106.88923 30145.59299 6.0357616 2000-2009 –> 2010-2019 12.519106 51281.5838 51356.7467
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 22692.23293 22735.60026 5.5287223 2000-2009 –> 2010-2019 11.269609 38380.6665 38464.7318
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 29785.22061 29830.37499 5.8305263 2000-2009 –> 2010-2019 12.151187 50798.3585 50886.0485
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 22787.09332 22830.46066 5.6015567 2000-2009 –> 2010-2019 11.440968 38559.7578 38643.8231
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 29838.50708 29883.66146 5.8638285 2000-2009 –> 2010-2019 12.195149 50862.4907 50950.1807
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 23273.78330 23317.15063 5.9906068 2000-2009 –> 2010-2019 12.145553 39306.9460 39391.0113
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 30026.27169 30071.42607 5.9826978 2000-2009 –> 2010-2019 12.462799 51199.7460 51287.4360
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 23121.34980 23164.71713 5.8659333 2000-2009 –> 2010-2019 11.882858 39042.2481 39126.3135
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 30107.97736 30153.13174 6.0351734 2000-2009 –> 2010-2019 12.651778 51415.6181 51503.3081
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 23310.96282 23354.33016 6.0214154 2000-2009 –> 2010-2019 12.208728 39370.0877 39454.1530
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 30219.64587 30264.80025 6.1076377 2000-2009 –> 2010-2019 12.760727 51562.6775 51650.3675
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 22880.26790 22923.63523 5.6740307 2000-2009 –> 2010-2019 11.600313 38726.0290 38810.0943
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 29967.06653 30012.22092 5.9449586 2000-2009 –> 2010-2019 12.320558 51035.9043 51123.5943
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 23378.11648 23421.48382 6.0774638 2000-2009 –> 2010-2019 12.300020 39465.3571 39549.4224
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 30268.66172 30313.81610 6.1397195 2000-2009 –> 2010-2019 12.808161 51626.5274 51714.2174
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 15946.49359 15988.24025 3.8648772 2000-2009 –> 2010-2019 8.051398 29878.1407 29960.4885
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 16291.48026 16333.22692 4.1038584 2000-2009 –> 2010-2019 8.396213 30345.0093 30427.3571
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 17499.40257 17541.34609 4.6537499 2000-2009 –> 2010-2019 9.227451 32186.1946 32268.8953
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 16423.25924 16465.00590 4.1989970 2000-2009 –> 2010-2019 8.609157 30608.9707 30691.3185
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 17518.90551 17560.84903 4.6691222 2000-2009 –> 2010-2019 9.259992 32224.4014 32307.1020
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + phosphate 17783.56068 17819.51227 4.8844656 2000-2009 –> 2010-2019 9.748307 32775.3902 32846.2764
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 15498.84850 15540.59515 3.5754055 2000-2009 –> 2010-2019 7.667946 29319.6549 29402.0027
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 17362.22915 17404.17266 4.5470498 2000-2009 –> 2010-2019 9.224271 32160.7495 32243.4502
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 15607.68616 15649.43282 3.6437264 2000-2009 –> 2010-2019 7.858836 29572.5581 29654.9059
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 17433.49844 17475.44196 4.6021776 2000-2009 –> 2010-2019 9.315693 32270.6065 32353.3071
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 16466.61649 16508.36314 4.2307787 2000-2009 –> 2010-2019 8.673007 30687.6986 30770.0464
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 17522.16439 17564.10790 4.6716958 2000-2009 –> 2010-2019 9.268001 32233.5677 32316.2684
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate_star 17576.37760 17612.32918 4.7163125 2000-2009 –> 2010-2019 9.484305 32468.8505 32539.7367
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 17344.09741 17386.04093 4.5331304 2000-2009 –> 2010-2019 9.247151 32181.7405 32264.4412
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 15643.60284 15685.34950 3.6665577 2000-2009 –> 2010-2019 7.907049 29637.7840 29720.1318
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 17432.90919 17474.85270 4.6017191 2000-2009 –> 2010-2019 9.316997 32271.8832 32354.5838
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate_star 17431.99837 17467.94996 4.6025667 2000-2009 –> 2010-2019 9.337438 32289.6728 32360.5591
mod Indo-Pacific (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 16041.05343 16082.80008 3.9289614 2000-2009 –> 2010-2019 8.279125 30159.8930 30242.2408
mod Indo-Pacific (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 16076.24857 16117.99523 3.9530839 2000-2009 –> 2010-2019 8.273950 30162.0983 30244.4461
mod Indo-Pacific (26,26.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 16113.53334 16155.28000 3.9788002 2000-2009 –> 2010-2019 8.305664 30206.1545 30288.5023
obs Indo-Pacific (26.5,26.75] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 15825.13259 15860.09125 5.6752570 2000-2009 –> 2010-2019 11.056896 32413.8548 32484.1570
mod Indo-Pacific (26.5,26.75] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 14268.92941 14310.08445 3.5923801 2000-2009 –> 2010-2019 7.927286 27095.1644 27176.2563
mod Indo-Pacific (26.5,26.75] cstar_tref ~ aou + silicate + phosphate + phosphate_star 14308.04136 14343.31711 3.6204394 2000-2009 –> 2010-2019 8.234995 27409.8451 27479.3525
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 15670.84882 15711.63392 5.5010222 2000-2009 –> 2010-2019 10.528842 31898.1442 31980.1633
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 15695.42502 15736.21012 5.5280626 2000-2009 –> 2010-2019 10.614006 31984.1449 32066.1641
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + phosphate 15876.44722 15911.40588 5.7336608 2000-2009 –> 2010-2019 11.109614 32459.5202 32529.8224
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 14073.95644 14115.11148 3.4622415 2000-2009 –> 2010-2019 7.828106 26931.7872 27012.8790
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 15823.43389 15864.21899 5.6710704 2000-2009 –> 2010-2019 11.037977 32399.5063 32481.5255
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 14211.84045 14252.99549 3.5537765 2000-2009 –> 2010-2019 8.033125 27183.6078 27264.6996
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 15697.92274 15738.70784 5.5308182 2000-2009 –> 2010-2019 10.626036 31996.3788 32078.3980
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate_star 15718.55107 15753.50973 5.5558453 2000-2009 –> 2010-2019 10.763565 32131.7180 32202.0202
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 14729.67672 14770.83176 3.9196864 2000-2009 –> 2010-2019 8.532219 27831.5336 27912.6255
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 15704.22881 15745.01391 5.5377814 2000-2009 –> 2010-2019 10.677281 32048.9277 32130.9469
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 14233.52705 14274.68209 3.5683919 2000-2009 –> 2010-2019 8.056294 27213.7649 27294.8567
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 15859.56068 15900.34578 5.7120955 2000-2009 –> 2010-2019 11.109139 32465.5571 32547.5763
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 14372.53221 14413.68725 3.6635105 2000-2009 –> 2010-2019 8.092736 27294.3375 27375.4294
obs Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 15806.92494 15847.71004 5.6524213 2000-2009 –> 2010-2019 11.071854 32435.0455 32517.0647
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 14561.61592 14602.77096 3.7969801 2000-2009 –> 2010-2019 8.322725 27579.1355 27660.2274
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + silicate + phosphate 14418.86889 14454.14464 3.6971771 2000-2009 –> 2010-2019 8.446481 27648.4675 27717.9748
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 14599.63768 14640.79272 3.8244003 2000-2009 –> 2010-2019 8.366932 27633.5961 27714.6880
mod Indo-Pacific (26.75,27] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 18487.14886 18523.85274 3.8072270 2000-2009 –> 2010-2019 7.747617 36869.9453 36943.2447
obs Indo-Pacific (26.75,27] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 24821.34008 24859.13131 5.3033370 2000-2009 –> 2010-2019 10.344637 49411.6003 49487.2260
mod Indo-Pacific (26.75,27] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 17497.76106 17540.58225 3.2838649 2000-2009 –> 2010-2019 7.216756 35870.0161 35955.5321
obs Indo-Pacific (26.75,27] cstar_tref ~ aou + phosphate + phosphate_star 25130.50149 25161.99419 5.5127146 2000-2009 –> 2010-2019 10.732049 49999.6718 50062.6932
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + phosphate 17821.88216 17864.70335 3.4465319 2000-2009 –> 2010-2019 7.156352 35809.6879 35895.2039
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + phosphate 24998.28867 25042.37845 5.4200598 2000-2009 –> 2010-2019 10.558597 49745.1804 49833.4103
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 17996.67489 18039.49608 3.5375747 2000-2009 –> 2010-2019 7.322058 36115.6733 36201.1893
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate 25245.62063 25289.71040 5.5894721 2000-2009 –> 2010-2019 10.754266 50033.7663 50121.9962
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 18018.71558 18061.53677 3.5492243 2000-2009 –> 2010-2019 7.337237 36143.8520 36229.3679
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + phosphate_star 18506.78895 18543.49283 3.8183970 2000-2009 –> 2010-2019 7.628462 36668.1423 36741.4417
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 17935.70648 17978.52767 3.5055486 2000-2009 –> 2010-2019 7.311061 36091.1894 36176.7054
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 25130.36650 25174.45628 5.5098787 2000-2009 –> 2010-2019 10.597209 49796.1946 49884.4245
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + phosphate 25064.87356 25102.66480 5.4665164 2000-2009 –> 2010-2019 10.667178 49907.0349 49982.6605
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 17480.37921 17523.20040 3.2753617 2000-2009 –> 2010-2019 7.211383 35857.8719 35943.3878
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 25036.00113 25080.09090 5.4455556 2000-2009 –> 2010-2019 10.531237 49699.2064 49787.4363
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 17663.83128 17706.65247 3.3662281 2000-2009 –> 2010-2019 7.396583 36197.2599 36282.7758
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 25178.32135 25222.41112 5.5428572 2000-2009 –> 2010-2019 10.707541 49966.2939 50054.5238
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 17691.92378 17734.74497 3.3803636 2000-2009 –> 2010-2019 7.415674 36233.4437 36318.9597
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 25178.83915 25222.92892 5.5432144 2000-2009 –> 2010-2019 10.706544 49964.6902 50052.9201
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate_star 25177.17491 25214.96615 5.5434460 2000-2009 –> 2010-2019 10.708221 49963.2903 50038.9159
mod Indo-Pacific (27,27.25] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 21219.69659 21258.07352 2.6503897 2000-2009 –> 2010-2019 5.831415 43164.9695 43241.4844
mod Indo-Pacific (27,27.25] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 21172.30987 21217.08295 2.6356572 2000-2009 –> 2010-2019 5.647711 42654.8785 42744.1459
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 31766.09969 31812.08952 4.9185960 2000-2009 –> 2010-2019 9.671853 62927.5598 63019.4847
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 31971.87322 32017.86304 5.0155471 2000-2009 –> 2010-2019 9.876905 63368.5545 63460.4795
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate 31767.15437 31806.57422 4.9200214 2000-2009 –> 2010-2019 9.674341 62928.9513 63007.7442
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate_star 31976.39133 32015.81119 5.0186492 2000-2009 –> 2010-2019 9.899200 63412.2871 63491.0799
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 31900.34561 31946.33544 4.9816318 2000-2009 –> 2010-2019 9.801640 63207.6753 63299.6002
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + silicate + phosphate + phosphate_star 31906.37759 31945.79744 4.9854287 2000-2009 –> 2010-2019 9.827297 63259.0394 63337.8322
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate 21162.86065 21207.63373 2.6328477 2000-2009 –> 2010-2019 5.798122 43067.8729 43157.1403
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 21561.73631 21606.50939 2.7540867 2000-2009 –> 2010-2019 6.016211 43723.3502 43812.6175
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate 31766.31737 31812.30720 4.9186975 2000-2009 –> 2010-2019 9.670665 62924.9388 63016.8638
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 31943.04602 31989.03585 5.0018508 2000-2009 –> 2010-2019 9.812803 63230.7044 63322.6294
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 21589.48250 21634.25558 2.7627250 2000-2009 –> 2010-2019 6.032678 43771.5052 43860.7725
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 22057.99586 22102.76895 2.9127484 2000-2009 –> 2010-2019 6.155442 44168.7446 44258.0119
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 32033.84159 32079.83142 5.0451165 2000-2009 –> 2010-2019 9.897883 63412.0207 63503.9457
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 31898.37976 31944.36958 4.9807029 2000-2009 –> 2010-2019 9.780458 63161.6656 63253.5906
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + phosphate 21549.66452 21588.04145 2.7509577 2000-2009 –> 2010-2019 5.978853 43619.4714 43695.9863
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 21550.40902 21595.18211 2.7505679 2000-2009 –> 2010-2019 5.871431 43335.1176 43424.3849
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 21909.89015 21954.66323 2.8644630 2000-2009 –> 2010-2019 6.104141 44012.7177 44101.9850
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 21930.63474 21975.40782 2.8711776 2000-2009 –> 2010-2019 6.114601 44043.2998 44132.5671
obs Indo-Pacific (27.25,27.5] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 33847.30684 33887.39189 4.2794151 2000-2009 –> 2010-2019 8.643938 66596.0271 66675.9582
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 33844.50071 33891.26660 4.2776692 2000-2009 –> 2010-2019 8.585319 66446.7694 66540.0223
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 29328.93038 29375.43963 3.1994048 2000-2009 –> 2010-2019 6.266425 57025.0612 57117.7940
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 33537.94755 33584.71344 4.1677681 2000-2009 –> 2010-2019 8.406846 65958.6003 66051.8532
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + silicate + phosphate 33542.90307 33582.98812 4.1702301 2000-2009 –> 2010-2019 8.420331 65990.9464 66070.8774
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 33785.14035 33831.90624 4.2561642 2000-2009 –> 2010-2019 8.590191 66456.5006 66549.7536
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + silicate + phosphate + phosphate_star 33785.70797 33825.79302 4.2570922 2000-2009 –> 2010-2019 8.634333 66567.3594 66647.2904
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 28902.83770 28949.34694 3.0815626 2000-2009 –> 2010-2019 6.267116 57012.2903 57105.0232
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 29163.75621 29210.26546 3.1531979 2000-2009 –> 2010-2019 6.407790 57506.9172 57599.6500
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 33519.69736 33566.46325 4.1613151 2000-2009 –> 2010-2019 8.356190 65821.7103 65914.9632
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 33733.35020 33780.11609 4.2374901 2000-2009 –> 2010-2019 8.484537 66175.2589 66268.5118
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 29183.13894 29229.64819 3.1585854 2000-2009 –> 2010-2019 6.419221 57546.5201 57639.2529
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 28789.92395 28836.43320 3.0510688 2000-2009 –> 2010-2019 6.102420 56430.2265 56522.9593
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 33676.03377 33722.79966 4.2169189 2000-2009 –> 2010-2019 8.468698 66130.5468 66223.7997
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 33680.10283 33726.86872 4.2183760 2000-2009 –> 2010-2019 8.456078 66097.0793 66190.3323
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + phosphate 28904.36489 28944.22996 3.0825201 2000-2009 –> 2010-2019 6.332736 57230.8573 57310.3426
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 28579.59344 28626.10269 2.9950687 2000-2009 –> 2010-2019 6.127676 56506.3590 56599.0918
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 28860.99571 28907.50495 3.0702273 2000-2009 –> 2010-2019 6.290925 57090.0448 57182.7777
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 28862.43550 28908.94475 3.0706166 2000-2009 –> 2010-2019 6.292429 57095.2549 57187.9877
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + nitrate + silicate + phosphate_star 28860.85164 28900.71671 3.0707292 2000-2009 –> 2010-2019 6.293493 57094.8898 57174.3751
obs Indo-Pacific (27.5,27.75] cstar_tref ~ aou + silicate + phosphate + phosphate_star 38447.08563 38488.07766 4.0009914 2000-2009 –> 2010-2019 7.670773 72670.4331 72751.9064
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 38975.66467 39023.48870 4.1577688 2000-2009 –> 2010-2019 7.798623 73101.4415 73196.4937
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 32813.83492 32861.13941 3.1888643 2000-2009 –> 2010-2019 6.097364 64954.3933 65049.1115
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 38213.27537 38261.09940 3.9327139 2000-2009 –> 2010-2019 7.771208 73004.1570 73099.2092
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 32989.40947 33036.71396 3.2331855 2000-2009 –> 2010-2019 6.139879 65121.9431 65216.6612
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate_star 33256.91360 33297.46031 3.3024192 2000-2009 –> 2010-2019 6.335964 65939.3358 66020.5228
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 32942.04772 32989.35221 3.2211694 2000-2009 –> 2010-2019 6.123887 65056.8956 65151.6137
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 33184.62742 33231.93190 3.2831890 2000-2009 –> 2010-2019 6.278978 65707.2364 65801.9546
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 38001.36301 38049.18704 3.8723505 2000-2009 –> 2010-2019 7.456395 71929.2696 72024.3217
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 33244.30524 33291.60972 3.2986287 2000-2009 –> 2010-2019 6.308168 65826.0774 65920.7956
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 38099.13268 38146.95671 3.9000843 2000-2009 –> 2010-2019 7.517004 72141.9233 72236.9754
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 33247.34581 33294.65030 3.2994173 2000-2009 –> 2010-2019 6.294928 65768.7545 65863.4727
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 32927.07057 32974.37505 3.2173789 2000-2009 –> 2010-2019 6.176488 65290.5105 65385.2287
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 38064.48036 38112.30439 3.8902320 2000-2009 –> 2010-2019 7.501569 72087.8349 72182.8870
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate 33284.94412 33325.49082 3.3097047 2000-2009 –> 2010-2019 6.336182 65937.2298 66018.4168
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 33215.32944 33262.63392 3.2911231 2000-2009 –> 2010-2019 6.293077 65764.4984 65859.2166
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 38084.91693 38132.74096 3.8960395 2000-2009 –> 2010-2019 7.511551 72122.8066 72217.8588
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate_star 38376.32106 38417.31309 3.9803783 2000-2009 –> 2010-2019 7.605533 72445.7246 72527.1979
obs Indo-Pacific (27.5,27.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 38520.84036 38568.66439 4.0220018 2000-2009 –> 2010-2019 7.737743 72902.7840 72997.8362
obs Indo-Pacific (27.5,27.75] cstar_tref ~ temp + aou + silicate + phosphate 38521.68691 38562.67893 4.0228376 2000-2009 –> 2010-2019 7.745863 72926.2729 73007.7461
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + phosphate 11556.77237 11597.49856 2.4683524 2000-2009 –> 2010-2019 4.709137 22699.5958 22781.0903
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 11596.89368 11637.61988 2.4883593 2000-2009 –> 2010-2019 4.737141 22757.5299 22839.0244
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + phosphate 11872.40929 11907.31746 2.6312571 2000-2009 –> 2010-2019 4.919449 23117.9125 23187.7650
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 11632.21374 11672.93993 2.5061062 2000-2009 –> 2010-2019 4.793865 22878.7691 22960.2636
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 21252.40742 21296.13796 3.9084568 2000-2009 –> 2010-2019 6.998918 37650.3773 37736.6391
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 11672.59525 11713.32145 2.5265514 2000-2009 –> 2010-2019 4.821805 22935.5053 23016.9999
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 21231.23696 21274.96750 3.8976329 2000-2009 –> 2010-2019 6.989598 37632.3344 37718.5962
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 11622.45629 11663.18248 2.5011908 2000-2009 –> 2010-2019 4.755754 22795.9303 22877.4248
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 21229.03813 21272.76867 3.8965104 2000-2009 –> 2010-2019 6.997239 37648.3411 37734.6029
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 11689.97553 11730.70173 2.5354023 2000-2009 –> 2010-2019 4.832979 22957.9420 23039.4366
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 21230.67566 21274.40620 3.8973464 2000-2009 –> 2010-2019 6.989715 37632.6139 37718.8757
mod Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 11541.34925 11582.07545 2.4607044 2000-2009 –> 2010-2019 4.861932 23029.9417 23111.4362
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 21257.07364 21300.80418 3.9108465 2000-2009 –> 2010-2019 7.001116 37654.6431 37740.9049
mod Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 11591.69907 11632.42526 2.4857598 2000-2009 –> 2010-2019 4.895679 23098.3581 23179.8527
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 21238.55439 21282.28493 3.9013707 2000-2009 –> 2010-2019 6.993327 37639.6351 37725.8969
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + silicate + phosphate_star 21243.91471 21281.39803 3.9051340 2000-2009 –> 2010-2019 6.997117 37643.0498 37716.9884
mod Indo-Pacific (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 11626.44982 11667.17602 2.5032014 2000-2009 –> 2010-2019 4.919206 23145.7188 23227.2133
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 21238.77263 21282.50317 3.9014823 2000-2009 –> 2010-2019 6.993828 37640.6616 37726.9234
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate_star 21236.97473 21274.45805 3.9015855 2000-2009 –> 2010-2019 7.002318 37654.2842 37728.2229
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + silicate + phosphate + phosphate_star 21242.94236 21280.42568 3.9046367 2000-2009 –> 2010-2019 6.997020 37642.9102 37716.8488
obs Indo-Pacific (27.85,27.95] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 27295.31392 27341.23166 3.3059267 2000-2009 –> 2010-2019 6.049874 48698.0982 48788.7478
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 13742.82261 13785.25964 2.1052696 2000-2009 –> 2010-2019 4.137604 26478.7446 26563.2005
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 27321.79626 27367.71401 3.3143280 2000-2009 –> 2010-2019 6.059396 48728.1773 48818.8270
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate 13413.83464 13450.20924 1.9995395 2000-2009 –> 2010-2019 3.980118 25993.5498 26065.9406
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + phosphate 13399.86548 13442.30250 1.9945142 2000-2009 –> 2010-2019 3.973377 25976.3989 26060.8548
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 13393.25634 13435.69337 1.9924380 2000-2009 –> 2010-2019 3.972090 25972.1728 26056.6287
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + silicate 13410.05493 13452.49196 1.9977192 2000-2009 –> 2010-2019 3.946574 25895.2422 25979.6981
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 13848.88385 13891.32088 2.1407508 2000-2009 –> 2010-2019 4.215953 26709.5874 26794.0433
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 27275.83131 27321.74905 3.2997595 2000-2009 –> 2010-2019 6.048237 48693.1448 48783.7944
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 13858.48979 13900.92682 2.1439937 2000-2009 –> 2010-2019 4.224132 26733.3968 26817.8527
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 27266.75726 27312.67501 3.2968911 2000-2009 –> 2010-2019 6.056470 48719.5760 48810.2256
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 13431.19269 13473.62972 2.0043845 2000-2009 –> 2010-2019 3.993259 26037.8952 26122.3511
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 13781.48968 13823.92671 2.1181365 2000-2009 –> 2010-2019 4.201160 26664.6830 26749.1389
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 13871.91844 13914.35547 2.1485354 2000-2009 –> 2010-2019 4.230835 26753.0344 26837.4903
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 27266.64375 27312.56150 3.2968552 2000-2009 –> 2010-2019 6.058108 48724.8017 48815.4513
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 27258.44169 27304.35944 3.2942646 2000-2009 –> 2010-2019 6.042419 48674.7199 48765.3695
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 27255.23767 27301.15542 3.2932532 2000-2009 –> 2010-2019 6.052959 48708.4603 48799.1099
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + silicate + phosphate 27280.89110 27320.24916 3.3019929 2000-2009 –> 2010-2019 6.051139 48698.3480 48776.0476
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + silicate + phosphate_star 27273.88017 27313.23824 3.2997750 2000-2009 –> 2010-2019 6.059931 48726.5399 48804.2396
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 27255.24605 27301.16380 3.2932558 2000-2009 –> 2010-2019 6.054233 48712.5250 48803.1746
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 46541.13744 46583.61199 3.4342752 2000-2009 –> 2010-2019 6.285581 83247.9977 83331.9584
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 46541.45588 46591.00953 3.4339460 2000-2009 –> 2010-2019 6.272858 83185.5167 83283.4708
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + phosphate + phosphate_star 46697.95028 46733.34574 3.4655115 2000-2009 –> 2010-2019 6.320744 83423.2824 83493.2497
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + silicate + phosphate + phosphate_star 46697.32812 46739.80267 3.4649934 2000-2009 –> 2010-2019 6.312977 83386.8441 83470.8048
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 46543.11816 46592.67181 3.4342714 2000-2009 –> 2010-2019 6.285550 83251.8346 83349.7888
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + aou + phosphate + phosphate_star 46685.43900 46727.91355 3.4626455 2000-2009 –> 2010-2019 6.317724 83411.9714 83495.9321
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate 17741.87900 17780.33932 1.7411685 2000-2009 –> 2010-2019 3.484695 30598.2932 30673.2757
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate + silicate 17594.69204 17639.56241 1.7124937 2000-2009 –> 2010-2019 3.406533 30265.8306 30353.3102
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate 18023.61864 18062.07896 1.7966371 2000-2009 –> 2010-2019 3.641985 31249.1870 31324.1695
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate_star 18025.71318 18064.17350 1.7970560 2000-2009 –> 2010-2019 3.646347 31265.1638 31340.1462
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 18013.54358 18058.41396 1.7942239 2000-2009 –> 2010-2019 3.607792 31128.1261 31215.6057
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 46781.52470 46831.07835 3.4812693 2000-2009 –> 2010-2019 6.339058 83524.1728 83622.1270
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 18016.86923 18061.73961 1.7948882 2000-2009 –> 2010-2019 3.614081 31151.5937 31239.0732
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 18220.32800 18265.19838 1.8360004 2000-2009 –> 2010-2019 3.692246 31486.1911 31573.6706
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + phosphate + phosphate_star 18031.90339 18070.36371 1.7982946 2000-2009 –> 2010-2019 3.649561 31278.2979 31353.2804
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 18022.77504 18067.64541 1.7960685 2000-2009 –> 2010-2019 3.616733 31162.7590 31250.2386
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate_star 18263.38496 18301.84528 1.8452315 2000-2009 –> 2010-2019 3.716280 31578.9110 31653.8935
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 46629.77355 46679.32719 3.4512803 2000-2009 –> 2010-2019 6.303361 83342.6750 83440.6292
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 46744.04255 46793.59620 3.4738379 2000-2009 –> 2010-2019 6.353970 83602.5414 83700.4955
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 46728.63993 46778.19357 3.4707887 2000-2009 –> 2010-2019 6.349066 83577.5571 83675.5112
obs Indo-Pacific (28.05,28.1] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 35852.07589 35893.33096 2.9591710 2000-2009 –> 2010-2019 5.687474 58714.2107 58794.2181
obs Indo-Pacific (28.05,28.1] cstar_tref ~ aou + phosphate + phosphate_star 35926.16993 35960.54916 2.9749440 2000-2009 –> 2010-2019 5.704906 58792.0388 58858.7116
obs Indo-Pacific (28.05,28.1] cstar_tref ~ aou + silicate + phosphate + phosphate_star 35902.78533 35944.04041 2.9696729 2000-2009 –> 2010-2019 5.698657 58767.2745 58847.2819
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 35868.59593 35916.72685 2.9621743 2000-2009 –> 2010-2019 5.702727 58774.9853 58868.3273
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + phosphate 35867.10877 35915.23969 2.9618665 2000-2009 –> 2010-2019 5.716508 58821.8613 58915.2032
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + silicate 9265.75300 9306.75536 1.4485583 2000-2009 –> 2010-2019 3.113485 16697.6084 16777.5423
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 9446.71790 9487.72027 1.5001599 2000-2009 –> 2010-2019 3.076453 16668.1776 16748.1115
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 35866.30790 35914.43882 2.9617008 2000-2009 –> 2010-2019 5.674453 58676.5308 58769.8728
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 9452.27545 9493.27782 1.5017734 2000-2009 –> 2010-2019 3.082556 16684.6735 16764.6074
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 35986.37987 36034.51079 2.9866493 2000-2009 –> 2010-2019 5.722214 58875.5872 58968.9292
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 9405.64960 9446.65197 1.4882905 2000-2009 –> 2010-2019 3.069064 16638.0263 16717.9602
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 35955.65067 36003.78159 2.9802445 2000-2009 –> 2010-2019 5.709951 58824.6362 58917.9782
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 9451.67556 9492.67793 1.5015991 2000-2009 –> 2010-2019 3.082389 16684.0920 16764.0258
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 35987.86093 36035.99185 2.9869584 2000-2009 –> 2010-2019 5.722160 58875.8170 58969.1590
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate_star 9456.12732 9491.27221 1.5034742 2000-2009 –> 2010-2019 3.084264 16686.5439 16755.0586
obs Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 36048.47613 36096.60705 2.9996340 2000-2009 –> 2010-2019 5.722166 58892.6416 58985.9836
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate 9809.20560 9844.35049 1.6097393 2000-2009 –> 2010-2019 3.212356 17092.3618 17160.8765
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate_star 9817.81883 9852.96372 1.6124233 2000-2009 –> 2010-2019 3.221063 17115.4021 17183.9168
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + nitrate + silicate + phosphate_star 9806.84309 9841.98798 1.6090038 2000-2009 –> 2010-2019 3.217185 17103.3295 17171.8442
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + silicate + phosphate + phosphate_star 9815.89842 9851.04330 1.6118245 2000-2009 –> 2010-2019 3.220450 17113.4475 17181.9622
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 18902.82673 18949.19398 1.3215178 2000-2009 –> 2010-2019 2.761150 36364.8930 36456.7278
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 37477.89906 37526.56186 2.7369697 2000-2009 –> 2010-2019 5.271115 67006.8551 67102.7372
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + phosphate + phosphate_star 37542.70513 37584.41610 2.7488346 2000-2009 –> 2010-2019 5.288110 67095.0735 67177.2581
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + silicate + phosphate 37512.40826 37554.11924 2.7434475 2000-2009 –> 2010-2019 5.279451 67048.5781 67130.7627
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + nitrate + silicate + phosphate 18987.18106 19026.92442 1.3318147 2000-2009 –> 2010-2019 2.793148 36593.6304 36672.3459
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 18989.17916 19035.54641 1.3318144 2000-2009 –> 2010-2019 2.786718 36554.4819 36646.3166
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 37592.64991 37641.31271 2.7573815 2000-2009 –> 2010-2019 5.326205 67292.4016 67388.2837
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + nitrate + phosphate 37635.96887 37684.63167 2.7651265 2000-2009 –> 2010-2019 5.339595 67363.3005 67459.1826
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + nitrate + silicate 18115.30832 18161.67557 1.2312121 2000-2009 –> 2010-2019 2.751798 36112.6327 36204.4674
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate 17796.10026 17842.46751 1.1963902 2000-2009 –> 2010-2019 2.569909 34798.2023 34890.0371
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate 37497.09463 37545.75744 2.7403736 2000-2009 –> 2010-2019 5.257666 66942.2025 67038.0846
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 17813.93446 17860.30171 1.1983094 2000-2009 –> 2010-2019 2.573727 34829.5488 34921.3835
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 37598.93120 37647.59401 2.7585032 2000-2009 –> 2010-2019 5.295791 67143.4626 67239.3446
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 18524.66650 18571.03375 1.2773557 2000-2009 –> 2010-2019 2.715568 35977.0733 36068.9081
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 17758.61188 17804.97913 1.1923658 2000-2009 –> 2010-2019 2.578367 34849.2279 34941.0626
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 37675.46172 37724.12453 2.7722064 2000-2009 –> 2010-2019 5.326772 67305.2729 67401.1549
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 17809.78784 17856.15509 1.1978629 2000-2009 –> 2010-2019 2.573006 34823.4545 34915.2892
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 37582.65888 37631.32168 2.7555982 2000-2009 –> 2010-2019 5.288857 67107.2203 67203.1023
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate_star 18059.40588 18099.14924 1.2252616 2000-2009 –> 2010-2019 2.614733 35172.4856 35251.2011
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate_star 37705.67625 37747.38722 2.7779950 2000-2009 –> 2010-2019 5.340554 67372.7430 67454.9276
obs Indo-Pacific (28.2,28.3] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 11303.87755 11339.18275 2.0290995 2000-2009 –> 2010-2019 4.639162 19554.3533 19622.4007
obs Indo-Pacific (28.2,28.3] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 11223.37906 11264.56846 1.9978181 2000-2009 –> 2010-2019 4.605899 19473.2229 19552.6115
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 11245.89137 11287.08077 2.0063061 2000-2009 –> 2010-2019 4.613869 19495.0472 19574.4358
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 32204.51369 32257.84018 0.7058531 2000-2009 –> 2010-2019 1.359917 54413.2638 54517.8310
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 35013.00008 35066.32657 0.7749602 2000-2009 –> 2010-2019 1.544393 60847.5849 60952.1520
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + phosphate + phosphate_star 11377.70701 11413.01221 2.0575088 2000-2009 –> 2010-2019 4.666707 19627.0351 19695.0825
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 11367.83406 11409.02346 2.0529134 2000-2009 –> 2010-2019 4.661977 19618.9827 19698.3713
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + nitrate + silicate + phosphate 32519.31596 32565.02437 0.7133295 2000-2009 –> 2010-2019 1.460256 57689.3088 57778.9378
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 32404.57717 32457.90366 0.7105653 2000-2009 –> 2010-2019 1.370325 54806.8658 54911.4329
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + nitrate + phosphate 11228.35327 11269.54267 1.9996905 2000-2009 –> 2010-2019 4.657863 19544.0979 19623.4866
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 11360.81121 11402.00061 2.0502001 2000-2009 –> 2010-2019 4.658433 19610.8575 19690.2461
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + silicate + phosphate 35268.98639 35322.31288 0.7815861 2000-2009 –> 2010-2019 1.612586 62821.6956 62926.2627
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 35279.94222 35333.26870 0.7818709 2000-2009 –> 2010-2019 1.613166 62840.5753 62945.1424
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 32363.41756 32416.74405 0.7095933 2000-2009 –> 2010-2019 1.448710 57300.7971 57405.3642
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 35277.63284 35330.95933 0.7818109 2000-2009 –> 2010-2019 1.613069 62837.2773 62941.8445
obs Indo-Pacific (28.2,28.3] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 11260.11328 11301.30268 2.0116868 2000-2009 –> 2010-2019 4.621674 19512.4890 19591.8777
mod Indo-Pacific (28.2,28.3] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 36557.02203 36610.34852 0.8157947 2000-2009 –> 2010-2019 1.651219 64228.2416 64332.8088
obs Indo-Pacific (28.2,28.3] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 11231.40512 11272.59452 2.0008401 2000-2009 –> 2010-2019 4.648804 19533.8208 19613.2095
obs Indo-Pacific (28.2,28.3] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 11374.10102 11415.29042 2.0553378 2000-2009 –> 2010-2019 4.664887 19625.8955 19705.2841
mod Indo-Pacific (28.2,28.3] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 36547.44385 36600.77034 0.8155348 2000-2009 –> 2010-2019 1.652222 64252.3766 64356.9437
mod Indo-Pacific (28.3,28.4] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 20687.63117 20728.09140 1.2587245 2000-2009 –> 2010-2019 2.189990 33635.1223 33714.4396
mod Indo-Pacific (28.3,28.4] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 20689.09349 20736.29709 1.2586705 2000-2009 –> 2010-2019 2.187613 33614.6156 33707.1524
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 20397.26210 20444.46570 1.2297123 2000-2009 –> 2010-2019 2.160542 33342.2550 33434.7917
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 4498.24789 4532.93704 2.0512579 2000-2009 –> 2010-2019 4.699501 7776.0747 7842.4391
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 4481.96413 4516.65327 2.0353986 2000-2009 –> 2010-2019 4.669903 7752.6964 7819.0608
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + silicate + phosphate_star 4502.23262 4531.96617 2.0571176 2000-2009 –> 2010-2019 4.707676 7779.2513 7836.1350
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 4529.77374 4564.46289 2.0823141 2000-2009 –> 2010-2019 4.715752 7799.9535 7866.3179
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + silicate + phosphate + phosphate_star 4531.04156 4560.77512 2.0855600 2000-2009 –> 2010-2019 4.720469 7799.9834 7856.8672
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + nitrate + phosphate 20373.74749 20420.95110 1.2274082 2000-2009 –> 2010-2019 2.155517 33290.6427 33383.1795
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 20454.80557 20502.00917 1.2353691 2000-2009 –> 2010-2019 2.168423 33422.7076 33515.2443
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + silicate + phosphate 4506.89542 4541.58457 2.0597302 2000-2009 –> 2010-2019 4.693258 7777.1219 7843.4863
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4480.94608 4515.63523 2.0344111 2000-2009 –> 2010-2019 4.684123 7759.5294 7825.8937
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 20372.04688 20419.25049 1.2272417 2000-2009 –> 2010-2019 2.156434 33300.1491 33392.6858
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 4506.23558 4540.92473 2.0590825 2000-2009 –> 2010-2019 4.709083 7784.9675 7851.3319
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4500.67160 4535.36075 2.0536290 2000-2009 –> 2010-2019 4.687992 7771.3305 7837.6949
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + aou + nitrate + phosphate 20667.94702 20708.40725 1.2567499 2000-2009 –> 2010-2019 2.185293 33587.3348 33666.6520
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 20669.78774 20716.99135 1.2567339 2000-2009 –> 2010-2019 2.182628 33563.7626 33656.2993
obs Indo-Pacific (28.3,28.4] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 4478.54998 4513.23913 2.0320890 2000-2009 –> 2010-2019 4.717484 7775.3795 7841.7439
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 20670.81548 20711.27571 1.2570374 2000-2009 –> 2010-2019 2.186716 33601.9390 33681.2562
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 20672.62700 20719.83061 1.2570185 2000-2009 –> 2010-2019 2.184108 33578.9812 33671.5179
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
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
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
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
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
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
db7df0e jens-daniel-mueller 2021-06-04
2edc791 jens-daniel-mueller 2021-06-04
207339d jens-daniel-mueller 2021-06-03
315710b jens-daniel-mueller 2021-06-03
be90356 jens-daniel-mueller 2021-06-02
d37a85d jens-daniel-mueller 2021-05-31
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
ac1a836 jens-daniel-mueller 2021-02-24
b03fbd3 jens-daniel-mueller 2021-02-24
3d3b4cc jens-daniel-mueller 2021-02-23
7b672f7 jens-daniel-mueller 2021-01-11
33ba23c jens-daniel-mueller 2021-01-07
318609d jens-daniel-mueller 2020-12-23
9d0b2d0 jens-daniel-mueller 2020-12-23
0aa2b50 jens-daniel-mueller 2020-12-23
2886da0 jens-daniel-mueller 2020-12-19
02f0ee9 jens-daniel-mueller 2020-12-18
7bcb4eb jens-daniel-mueller 2020-12-18
7131186 jens-daniel-mueller 2020-12-17
a84ff3c jens-daniel-mueller 2020-12-17
rm(max_rmse)

8.3 Predictor counts

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

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

# print table
lm_all_stats %>%
  gt(rowname_col = "gamma_slab",
     groupname_col = c("basin", "data_source")) %>% 
  summary_rows(
    groups = TRUE,
    fns = list(total = "sum")
  )
aou nitrate phosphate phosphate_star sal silicate temp
Atlantic - mod
(-Inf,26] 5 7 8 7 10 6 6
(26,26.5] 8 8 8 7 6 6 7
(26.5,26.75] 6 7 6 7 7 7 10
(26.75,27] 8 5 8 5 5 10 7
(27,27.25] 8 5 6 6 8 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 90.00 79.00 77.00 86.00 116.00 97.00 106.00
Indo-Pacific - mod
(-Inf,26] 8 5 7 7 10 5 7
(26,26.5] 7 6 7 7 7 7 9
(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] 7 7 5 5 10 6 9
(27.95,28.05] 6 3 4 6 10 6 10
(28.05,28.1] 5 3 4 7 6 10 10
(28.1,28.2] 4 6 6 5 10 10 7
(28.2,28.3] 5 7 8 5 8 10 6
(28.3,28.4] 7 10 9 7 4 3 7
total 87.00 93.00 89.00 89.00 100.00 99.00 114.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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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