Last updated: 2021-05-28

Checks: 7 0

Knit directory: emlr_obs_v_XXX/

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

R Markdown file: up-to-date

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

Environment: empty

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

Seed: set.seed(20200707)

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

Session information: recorded

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

Cache: none

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

File paths: relative

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

Repository version: a61014f

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 a61014f. 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:  code/local_job.R

Unstaged changes:
    Modified:   code/Workflowr_project_managment.R
    Modified:   data/auxillary/params_local.rds

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

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 ~ aou + nitrate + silicate + phosphate + phosphate_star 240.4224 253.9452 2.2275212 2000-2009 –> 2010-2019 6.297827 960.0866 993.4076
mod Atlantic (-Inf,26] cstar ~ aou + silicate + phosphate + phosphate_star 253.9937 265.5847 2.5949025 2000-2009 –> 2010-2019 6.680270 972.5813 1001.1421
mod Atlantic (-Inf,26] cstar ~ nitrate + silicate + phosphate + phosphate_star 241.5783 253.1692 2.2975172 2000-2009 –> 2010-2019 6.571520 971.4507 1000.0115
obs Atlantic (-Inf,26] cstar ~ sal + aou + nitrate + phosphate + phosphate_star 887.6208 907.5850 7.3422307 2000-2009 –> 2010-2019 15.271280 1278.4837 1312.3708
mod Atlantic (-Inf,26] cstar ~ sal + aou + silicate + phosphate + phosphate_star 253.3712 266.8940 2.5290357 2000-2009 –> 2010-2019 6.612622 973.8498 1007.1707
obs Atlantic (-Inf,26] cstar ~ sal + aou + silicate + phosphate + phosphate_star 887.6341 907.5983 7.3426114 2000-2009 –> 2010-2019 15.467755 1281.1355 1315.0226
mod Atlantic (-Inf,26] cstar ~ sal + nitrate + silicate + phosphate + phosphate_star 238.9933 252.5161 2.1965288 2000-2009 –> 2010-2019 6.336539 962.9025 996.2235
obs Atlantic (-Inf,26] cstar ~ sal + temp + aou + nitrate 887.4993 904.6115 7.3963057 2000-2009 –> 2010-2019 15.357501 1276.7992 1305.8453
obs Atlantic (-Inf,26] cstar ~ sal + temp + aou + nitrate + phosphate 887.9021 907.8663 7.3503026 2000-2009 –> 2010-2019 15.311491 1279.2019 1313.0890
obs Atlantic (-Inf,26] cstar ~ sal + temp + aou + nitrate + phosphate_star 887.9215 907.8858 7.3508615 2000-2009 –> 2010-2019 15.312049 1279.2214 1313.1085
obs Atlantic (-Inf,26] cstar ~ sal + temp + aou + nitrate + silicate 889.3926 909.3568 7.3932221 2000-2009 –> 2010-2019 14.961944 1275.2325 1309.1196
mod Atlantic (-Inf,26] cstar ~ sal + temp + aou + silicate + phosphate 253.5623 267.0851 2.5337780 2000-2009 –> 2010-2019 6.398451 960.2662 993.5871
obs Atlantic (-Inf,26] cstar ~ sal + temp + aou + silicate + phosphate 887.9241 907.8883 7.3509342 2000-2009 –> 2010-2019 15.546304 1282.3549 1316.2420
mod Atlantic (-Inf,26] cstar ~ sal + temp + nitrate + phosphate + phosphate_star 256.9240 270.4468 2.6186775 2000-2009 –> 2010-2019 6.824262 984.7619 1018.0829
obs Atlantic (-Inf,26] cstar ~ sal + temp + nitrate + phosphate + phosphate_star 887.9770 907.9412 7.3524545 2000-2009 –> 2010-2019 15.333074 1279.5401 1313.4272
mod Atlantic (-Inf,26] cstar ~ temp + aou + nitrate + silicate + phosphate 240.9119 254.4346 2.2382348 2000-2009 –> 2010-2019 6.000146 940.8782 974.1992
obs Atlantic (-Inf,26] cstar ~ temp + aou + nitrate + silicate + phosphate_star 897.6780 917.6422 7.6364178 2000-2009 –> 2010-2019 15.533933 1288.1106 1321.9977
mod Atlantic (-Inf,26] cstar ~ temp + aou + silicate + phosphate 253.5113 265.1023 2.5826598 2000-2009 –> 2010-2019 6.470364 959.7007 988.2615
mod Atlantic (-Inf,26] cstar ~ temp + nitrate + silicate + phosphate + phosphate_star 241.1507 254.6735 2.2434814 2000-2009 –> 2010-2019 6.515000 972.8776 1006.1986
obs Atlantic (-Inf,26] cstar ~ temp + nitrate + silicate + phosphate + phosphate_star 896.8062 916.7704 7.6104557 2000-2009 –> 2010-2019 15.491777 1287.0170 1320.9041
obs Atlantic (26,26.5] cstar ~ aou + nitrate + phosphate + phosphate_star 4049.9148 4076.4250 6.5184055 2000-2009 –> 2010-2019 13.190514 7617.6042 7669.8192
mod Atlantic (26,26.5] cstar ~ aou + nitrate + silicate + phosphate + phosphate_star 2483.5253 2511.6042 4.9900230 2000-2009 –> 2010-2019 10.119256 4226.0297 4279.6267
obs Atlantic (26,26.5] cstar ~ aou + nitrate + silicate + phosphate + phosphate_star 3994.3232 4025.2517 6.2192828 2000-2009 –> 2010-2019 12.565071 7510.2575 7571.1750
mod Atlantic (26,26.5] cstar ~ sal + aou + nitrate + phosphate + phosphate_star 2502.3042 2530.3831 5.1061916 2000-2009 –> 2010-2019 9.753065 4188.9096 4242.5066
obs Atlantic (26,26.5] cstar ~ sal + aou + nitrate + phosphate + phosphate_star 4031.6961 4062.6246 6.4117882 2000-2009 –> 2010-2019 13.082846 7601.2166 7662.1341
obs Atlantic (26,26.5] cstar ~ sal + aou + nitrate + silicate + phosphate 4043.9596 4074.8881 6.4762462 2000-2009 –> 2010-2019 13.135225 7611.5374 7672.4548
mod Atlantic (26,26.5] cstar ~ sal + aou + silicate + phosphate + phosphate_star 2462.5732 2490.6521 4.8635272 2000-2009 –> 2010-2019 9.734330 4175.8169 4229.4139
obs Atlantic (26,26.5] cstar ~ sal + aou + silicate + phosphate + phosphate_star 4049.6489 4080.5775 6.5063696 2000-2009 –> 2010-2019 13.099355 7606.5497 7667.4672
mod Atlantic (26,26.5] cstar ~ sal + temp + aou + nitrate + phosphate 2501.3627 2529.4415 5.1003034 2000-2009 –> 2010-2019 10.017836 4220.0106 4273.6076
mod Atlantic (26,26.5] cstar ~ sal + temp + aou + nitrate + phosphate_star 2501.3506 2529.4295 5.1002282 2000-2009 –> 2010-2019 10.047873 4223.4539 4277.0509
mod Atlantic (26,26.5] cstar ~ sal + temp + aou + silicate + phosphate 2462.0834 2490.1623 4.8606086 2000-2009 –> 2010-2019 10.061382 4212.4276 4266.0246
obs Atlantic (26,26.5] cstar ~ sal + temp + aou + silicate + phosphate 4051.7399 4082.6685 6.5174761 2000-2009 –> 2010-2019 13.141636 7613.6976 7674.6151
mod Atlantic (26,26.5] cstar ~ sal + temp + aou + silicate + phosphate_star 2462.0851 2490.1639 4.8606186 2000-2009 –> 2010-2019 10.134110 4220.2883 4273.8853
obs Atlantic (26,26.5] cstar ~ sal + temp + aou + silicate + phosphate_star 4054.7820 4085.7106 6.5336681 2000-2009 –> 2010-2019 13.231903 7628.6610 7689.5785
mod Atlantic (26,26.5] cstar ~ sal + temp + nitrate + phosphate + phosphate_star 2501.2960 2529.3748 5.0998864 2000-2009 –> 2010-2019 10.105259 4229.9649 4283.5619
mod Atlantic (26,26.5] cstar ~ temp + aou + nitrate + silicate + phosphate 2481.4845 2509.5634 4.9775587 2000-2009 –> 2010-2019 10.141849 4227.8442 4281.4412
obs Atlantic (26,26.5] cstar ~ temp + aou + nitrate + silicate + phosphate 3999.8795 4030.8081 6.2475331 2000-2009 –> 2010-2019 12.627507 7521.5734 7582.4908
mod Atlantic (26,26.5] cstar ~ temp + aou + nitrate + silicate + phosphate_star 2477.6336 2505.7124 4.9541234 2000-2009 –> 2010-2019 10.139446 4226.2937 4279.8907
obs Atlantic (26,26.5] cstar ~ temp + aou + nitrate + silicate + phosphate_star 4001.9876 4032.9161 6.2582847 2000-2009 –> 2010-2019 12.679229 7530.5435 7591.4610
obs Atlantic (26,26.5] cstar ~ temp + nitrate + silicate + phosphate + phosphate_star 4008.3357 4039.2643 6.2907739 2000-2009 –> 2010-2019 12.763841 7545.5588 7606.4763
mod Atlantic (26.5,26.75] cstar ~ aou + nitrate + silicate + phosphate + phosphate_star 3128.8563 3158.6331 4.8361338 2000-2009 –> 2010-2019 9.406399 5734.6218 5793.0220
mod Atlantic (26.5,26.75] cstar ~ sal + aou + nitrate + silicate + phosphate 3177.8978 3207.6746 5.0696456 2000-2009 –> 2010-2019 9.538699 5763.9113 5822.3115
mod Atlantic (26.5,26.75] cstar ~ sal + temp + aou + nitrate + phosphate 3130.1486 3159.9254 4.8421467 2000-2009 –> 2010-2019 9.470734 5747.0983 5805.4984
obs Atlantic (26.5,26.75] cstar ~ sal + temp + aou + nitrate + phosphate 5763.6678 5797.2612 5.9655692 2000-2009 –> 2010-2019 12.067885 10328.5935 10394.0943
mod Atlantic (26.5,26.75] cstar ~ sal + temp + aou + nitrate + phosphate_star 3132.9066 3162.6834 4.8550050 2000-2009 –> 2010-2019 9.476074 5748.4225 5806.8226
obs Atlantic (26.5,26.75] cstar ~ sal + temp + aou + nitrate + phosphate_star 5760.8055 5794.3989 5.9560590 2000-2009 –> 2010-2019 12.070713 10328.5792 10394.0800
obs Atlantic (26.5,26.75] cstar ~ sal + temp + aou + silicate + phosphate 5784.7145 5818.3079 6.0359679 2000-2009 –> 2010-2019 12.139014 10349.8090 10415.3098
obs Atlantic (26.5,26.75] cstar ~ sal + temp + aou + silicate + phosphate_star 5780.7786 5814.3719 6.0227398 2000-2009 –> 2010-2019 12.129598 10346.7534 10412.2542
mod Atlantic (26.5,26.75] cstar ~ sal + temp + nitrate + phosphate + phosphate_star 3139.1067 3168.8835 4.8840353 2000-2009 –> 2010-2019 9.511965 5755.9312 5814.3313
obs Atlantic (26.5,26.75] cstar ~ sal + temp + nitrate + phosphate + phosphate_star 5758.5567 5792.1501 5.9485976 2000-2009 –> 2010-2019 12.068883 10327.6284 10393.1292
obs Atlantic (26.5,26.75] cstar ~ sal + temp + silicate + phosphate + phosphate_star 5780.1455 5813.7389 6.0206149 2000-2009 –> 2010-2019 12.129225 10346.5248 10412.0256
mod Atlantic (26.5,26.75] cstar ~ temp + aou + nitrate + phosphate 3135.2977 3160.8207 4.8755473 2000-2009 –> 2010-2019 9.558816 5760.6062 5810.6635
mod Atlantic (26.5,26.75] cstar ~ temp + aou + nitrate + phosphate_star 3137.9541 3163.4770 4.8880161 2000-2009 –> 2010-2019 9.561335 5761.3867 5811.4439
mod Atlantic (26.5,26.75] cstar ~ temp + aou + nitrate + silicate + phosphate 3136.9586 3166.7354 4.8739578 2000-2009 –> 2010-2019 9.404515 5735.0276 5793.4277
obs Atlantic (26.5,26.75] cstar ~ temp + aou + nitrate + silicate + phosphate 5788.1700 5821.7633 6.0476051 2000-2009 –> 2010-2019 12.164112 10356.3707 10421.8715
mod Atlantic (26.5,26.75] cstar ~ temp + aou + nitrate + silicate + phosphate_star 3139.7384 3169.5152 4.8870025 2000-2009 –> 2010-2019 9.423799 5739.0212 5797.4213
obs Atlantic (26.5,26.75] cstar ~ temp + aou + nitrate + silicate + phosphate_star 5784.0510 5817.6444 6.0337358 2000-2009 –> 2010-2019 12.154724 10353.2846 10418.7854
mod Atlantic (26.5,26.75] cstar ~ temp + nitrate + silicate + phosphate + phosphate_star 3144.5954 3174.3722 4.9098795 2000-2009 –> 2010-2019 9.473746 5749.1253 5807.5254
obs Atlantic (26.5,26.75] cstar ~ temp + nitrate + silicate + phosphate + phosphate_star 5782.5251 5816.1185 6.0286060 2000-2009 –> 2010-2019 12.151169 10352.1213 10417.6221
obs Atlantic (26.5,26.75] cstar ~ temp + silicate + phosphate + phosphate_star 5789.3396 5818.1339 6.0582993 2000-2009 –> 2010-2019 12.181198 10357.0132 10413.1567
mod Atlantic (26.75,27] cstar ~ aou + nitrate + silicate + phosphate + phosphate_star 5380.9318 5415.5066 3.2586699 2000-2009 –> 2010-2019 6.321165 9456.0012 9523.3683
mod Atlantic (26.75,27] cstar ~ aou + silicate + phosphate + phosphate_star 5659.7639 5689.3995 3.7336310 2000-2009 –> 2010-2019 6.809893 9740.0099 9797.7531
mod Atlantic (26.75,27] cstar ~ sal + aou + nitrate + silicate + phosphate 5440.1256 5474.7004 3.3534789 2000-2009 –> 2010-2019 6.631355 9623.9404 9691.3075
mod Atlantic (26.75,27] cstar ~ sal + temp + aou + silicate + phosphate 5566.5575 5601.1322 3.5653208 2000-2009 –> 2010-2019 6.587308 9620.3225 9687.6896
obs Atlantic (26.75,27] cstar ~ sal + temp + aou + silicate + phosphate 7008.8704 7044.1050 5.2868336 2000-2009 –> 2010-2019 11.109799 12468.3115 12536.8119
mod Atlantic (26.75,27] cstar ~ sal + temp + aou + silicate + phosphate_star 5576.6746 5611.2494 3.5828399 2000-2009 –> 2010-2019 6.639125 9648.4965 9715.8635
obs Atlantic (26.75,27] cstar ~ sal + temp + aou + silicate + phosphate_star 7001.9971 7037.2316 5.2708356 2000-2009 –> 2010-2019 11.059764 12451.4015 12519.9020
obs Atlantic (26.75,27] cstar ~ sal + temp + nitrate + silicate + phosphate_star 6997.9079 7033.1425 5.2613411 2000-2009 –> 2010-2019 11.030842 12441.5573 12510.0577
mod Atlantic (26.75,27] cstar ~ sal + temp + silicate + phosphate + phosphate_star 5580.7160 5615.2908 3.5898621 2000-2009 –> 2010-2019 6.668508 9664.2017 9731.5687
obs Atlantic (26.75,27] cstar ~ sal + temp + silicate + phosphate + phosphate_star 7001.8407 7037.0753 5.2704724 2000-2009 –> 2010-2019 11.056955 12450.5215 12519.0220
mod Atlantic (26.75,27] cstar ~ temp + aou + nitrate + silicate + phosphate 5365.0598 5399.6345 3.2337071 2000-2009 –> 2010-2019 6.243269 9412.2329 9479.5999
obs Atlantic (26.75,27] cstar ~ temp + aou + nitrate + silicate + phosphate 7005.3185 7040.5531 5.2785604 2000-2009 –> 2010-2019 11.119485 12470.0313 12538.5318
mod Atlantic (26.75,27] cstar ~ temp + aou + nitrate + silicate + phosphate_star 5365.1621 5399.7369 3.2338675 2000-2009 –> 2010-2019 6.250061 9415.8566 9483.2237
obs Atlantic (26.75,27] cstar ~ temp + aou + nitrate + silicate + phosphate_star 6996.1772 7031.4117 5.2573275 2000-2009 –> 2010-2019 11.044803 12445.1519 12513.6523
mod Atlantic (26.75,27] cstar ~ temp + aou + silicate + phosphate 5676.3230 5705.9585 3.7637056 2000-2009 –> 2010-2019 6.828761 9750.7298 9808.4730
obs Atlantic (26.75,27] cstar ~ temp + aou + silicate + phosphate_star 7000.3844 7030.5855 5.2717360 2000-2009 –> 2010-2019 11.132466 12468.8925 12527.6072
mod Atlantic (26.75,27] cstar ~ temp + nitrate + silicate + phosphate + phosphate_star 5369.5302 5404.1049 3.2407185 2000-2009 –> 2010-2019 6.275398 9430.0009 9497.3679
obs Atlantic (26.75,27] cstar ~ temp + nitrate + silicate + phosphate + phosphate_star 6996.6545 7031.8890 5.2584341 2000-2009 –> 2010-2019 11.053017 12447.7301 12516.2306
obs Atlantic (26.75,27] cstar ~ temp + nitrate + silicate + phosphate_star 6996.7680 7026.9691 5.2633367 2000-2009 –> 2010-2019 11.087583 12454.5853 12513.3000
obs Atlantic (26.75,27] cstar ~ temp + silicate + phosphate + phosphate_star 7000.2005 7030.4015 5.2713084 2000-2009 –> 2010-2019 11.129615 12468.0007 12526.7154
mod Atlantic (27,27.25] cstar ~ sal + aou + nitrate + silicate + phosphate 6485.3552 6521.9488 2.5366990 2000-2009 –> 2010-2019 5.205336 11895.6938 11967.4600
mod Atlantic (27,27.25] cstar ~ sal + aou + silicate + phosphate 6627.6852 6659.0512 2.6731861 2000-2009 –> 2010-2019 5.400889 12085.2365 12146.7504
mod Atlantic (27,27.25] cstar ~ sal + aou + silicate + phosphate + phosphate_star 6626.3309 6662.9246 2.6699322 2000-2009 –> 2010-2019 5.218099 11932.8261 12004.5923
obs Atlantic (27,27.25] cstar ~ sal + temp + aou + nitrate + phosphate 12006.3624 12045.0945 5.9852080 2000-2009 –> 2010-2019 11.678197 22322.3260 22398.8368
obs Atlantic (27,27.25] cstar ~ sal + temp + aou + nitrate + phosphate_star 11984.2649 12022.9970 5.9499304 2000-2009 –> 2010-2019 11.556396 22250.2713 22326.7820
mod Atlantic (27,27.25] cstar ~ sal + temp + aou + nitrate + silicate 6670.5580 6707.1517 2.7131554 2000-2009 –> 2010-2019 5.203622 11925.5648 11997.3310
mod Atlantic (27,27.25] cstar ~ sal + temp + aou + silicate + phosphate 6625.1147 6661.7083 2.6687534 2000-2009 –> 2010-2019 5.281897 11988.2139 12059.9801
mod Atlantic (27,27.25] cstar ~ sal + temp + aou + silicate + phosphate_star 6651.3252 6687.9189 2.6942739 2000-2009 –> 2010-2019 5.391124 12085.3050 12157.0712
obs Atlantic (27,27.25] cstar ~ sal + temp + aou + silicate + phosphate_star 12034.2418 12072.9739 6.0300150 2000-2009 –> 2010-2019 11.611802 22285.8574 22362.3682
obs Atlantic (27,27.25] cstar ~ sal + temp + nitrate + phosphate + phosphate_star 11990.5690 12029.3012 5.9599735 2000-2009 –> 2010-2019 11.591451 22271.0955 22347.6063
obs Atlantic (27,27.25] cstar ~ sal + temp + nitrate + silicate + phosphate_star 12049.9698 12088.7019 6.0554402 2000-2009 –> 2010-2019 11.686428 22330.2126 22406.7234
obs Atlantic (27,27.25] cstar ~ sal + temp + silicate + phosphate + phosphate_star 12034.4075 12073.1396 6.0302823 2000-2009 –> 2010-2019 11.615593 22288.0820 22364.5927
obs Atlantic (27,27.25] cstar ~ temp + aou + nitrate + phosphate_star 12036.7712 12069.9702 6.0373260 2000-2009 –> 2010-2019 11.667704 22314.6608 22380.2415
mod Atlantic (27,27.25] cstar ~ temp + aou + nitrate + silicate 6762.6753 6794.0413 2.8074794 2000-2009 –> 2010-2019 5.332570 12046.7196 12108.2335
mod Atlantic (27,27.25] cstar ~ temp + aou + nitrate + silicate + phosphate 6757.9821 6794.5757 2.8006645 2000-2009 –> 2010-2019 5.321878 12040.5718 12112.3380
obs Atlantic (27,27.25] cstar ~ temp + aou + nitrate + silicate + phosphate 11987.4723 12026.2045 5.9550380 2000-2009 –> 2010-2019 11.583322 22266.1484 22342.6591
mod Atlantic (27,27.25] cstar ~ temp + aou + nitrate + silicate + phosphate_star 6748.3558 6784.9494 2.7908922 2000-2009 –> 2010-2019 5.313823 12032.4769 12104.2430
obs Atlantic (27,27.25] cstar ~ temp + aou + nitrate + silicate + phosphate_star 11961.8726 12000.6047 5.9143942 2000-2009 –> 2010-2019 11.448045 22185.2355 22261.7463
mod Atlantic (27,27.25] cstar ~ temp + nitrate + silicate + phosphate + phosphate_star 6652.6118 6689.2055 2.6955328 2000-2009 –> 2010-2019 5.302906 12010.7416 12082.5078
obs Atlantic (27,27.25] cstar ~ temp + nitrate + silicate + phosphate + phosphate_star 11967.3064 12006.0386 5.9229981 2000-2009 –> 2010-2019 11.478993 22203.8144 22280.3252
mod Atlantic (27.25,27.5] cstar ~ sal + aou + nitrate + silicate 13147.8178 13182.4036 3.9352816 2000-2009 –> 2010-2019 7.414532 24979.7768 25048.5860
mod Atlantic (27.25,27.5] cstar ~ sal + aou + nitrate + silicate + phosphate 12873.6677 12914.0178 3.7111873 2000-2009 –> 2010-2019 7.137588 24639.7577 24720.0352
mod Atlantic (27.25,27.5] cstar ~ sal + aou + nitrate + silicate + phosphate_star 12904.8980 12945.2481 3.7358766 2000-2009 –> 2010-2019 7.215056 24738.7670 24819.0444
obs Atlantic (27.25,27.5] cstar ~ sal + temp + aou + nitrate + phosphate 14602.0085 14642.2114 5.7208433 2000-2009 –> 2010-2019 10.369934 26154.7108 26233.9500
mod Atlantic (27.25,27.5] cstar ~ sal + temp + aou + nitrate + phosphate_star 12992.3472 13032.6972 3.8058874 2000-2009 –> 2010-2019 7.232295 24758.4470 24838.7244
obs Atlantic (27.25,27.5] cstar ~ sal + temp + aou + nitrate + phosphate_star 14608.3018 14648.5047 5.7286549 2000-2009 –> 2010-2019 10.348064 26135.9993 26215.2385
mod Atlantic (27.25,27.5] cstar ~ sal + temp + aou + nitrate + silicate 13109.9879 13150.3379 3.9021433 2000-2009 –> 2010-2019 7.201891 24709.0746 24789.3520
obs Atlantic (27.25,27.5] cstar ~ sal + temp + aou + phosphate_star 14641.0376 14675.4972 5.7719639 2000-2009 –> 2010-2019 10.476897 26238.3535 26306.2728
obs Atlantic (27.25,27.5] cstar ~ sal + temp + aou + silicate + phosphate 14635.9236 14676.1265 5.7630675 2000-2009 –> 2010-2019 10.474005 26240.2185 26319.4577
mod Atlantic (27.25,27.5] cstar ~ sal + temp + aou + silicate + phosphate_star 12757.7421 12798.0922 3.6209601 2000-2009 –> 2010-2019 7.381440 24936.3550 25016.6325
obs Atlantic (27.25,27.5] cstar ~ sal + temp + aou + silicate + phosphate_star 14640.3458 14680.5487 5.7685960 2000-2009 –> 2010-2019 10.472005 26238.3964 26317.6356
mod Atlantic (27.25,27.5] cstar ~ sal + temp + nitrate + phosphate + phosphate_star 12993.8200 13034.1701 3.8070777 2000-2009 –> 2010-2019 7.194091 24708.6451 24788.9225
obs Atlantic (27.25,27.5] cstar ~ sal + temp + nitrate + phosphate + phosphate_star 14600.6347 14640.8375 5.7191394 2000-2009 –> 2010-2019 10.347314 26135.7330 26214.9722
mod Atlantic (27.25,27.5] cstar ~ sal + temp + nitrate + silicate + phosphate_star 12748.4881 12788.8382 3.6138527 2000-2009 –> 2010-2019 7.317068 24859.0604 24939.3378
obs Atlantic (27.25,27.5] cstar ~ sal + temp + phosphate + phosphate_star 14636.9500 14671.4096 5.7668505 2000-2009 –> 2010-2019 10.473034 26235.3029 26303.2222
mod Atlantic (27.25,27.5] cstar ~ sal + temp + silicate + phosphate + phosphate_star 12759.9196 12800.2696 3.6226344 2000-2009 –> 2010-2019 7.380287 24935.1978 25015.4753
obs Atlantic (27.25,27.5] cstar ~ sal + temp + silicate + phosphate + phosphate_star 14636.0568 14676.2596 5.7632339 2000-2009 –> 2010-2019 10.466571 26234.0475 26313.2867
obs Atlantic (27.25,27.5] cstar ~ temp + aou + nitrate + silicate + phosphate_star 14619.4089 14659.6118 5.7424680 2000-2009 –> 2010-2019 10.480989 26246.4954 26325.7345
mod Atlantic (27.25,27.5] cstar ~ temp + nitrate + silicate + phosphate + phosphate_star 13193.6020 13233.9521 3.9720346 2000-2009 –> 2010-2019 7.223186 24726.9021 24807.1795
obs Atlantic (27.25,27.5] cstar ~ temp + nitrate + silicate + phosphate + phosphate_star 14610.1559 14650.3588 5.7309584 2000-2009 –> 2010-2019 10.472550 26239.7715 26319.0106
obs Atlantic (27.5,27.75] cstar ~ aou + nitrate + silicate + phosphate + phosphate_star 15882.6090 15923.6680 5.0819293 2000-2009 –> 2010-2019 10.155324 29869.8295 29951.0610
obs Atlantic (27.5,27.75] cstar ~ aou + silicate + phosphate + phosphate_star 15917.5043 15952.6977 5.1180315 2000-2009 –> 2010-2019 10.195741 29906.6285 29976.2555
mod Atlantic (27.5,27.75] cstar ~ sal + aou + nitrate + phosphate 20599.5045 20638.0074 2.3547315 2000-2009 –> 2010-2019 4.841386 39805.0277 39881.4708
mod Atlantic (27.5,27.75] cstar ~ sal + aou + nitrate + silicate + phosphate 20178.1527 20223.0728 2.2470925 2000-2009 –> 2010-2019 4.692927 39249.3237 39338.5073
mod Atlantic (27.5,27.75] cstar ~ sal + aou + silicate + phosphate 20643.6988 20682.2017 2.3662611 2000-2009 –> 2010-2019 4.837994 39799.6400 39876.0831
mod Atlantic (27.5,27.75] cstar ~ sal + aou + silicate + phosphate + phosphate_star 20629.2443 20674.1644 2.3619618 2000-2009 –> 2010-2019 4.831010 39778.2328 39867.4164
obs Atlantic (27.5,27.75] cstar ~ sal + aou + silicate + phosphate + phosphate_star 15915.6247 15956.6837 5.1142232 2000-2009 –> 2010-2019 10.146283 29865.2791 29946.5105
mod Atlantic (27.5,27.75] cstar ~ sal + temp + aou + nitrate + phosphate 20521.9715 20566.8915 2.3341238 2000-2009 –> 2010-2019 4.819158 39724.1252 39813.3088
mod Atlantic (27.5,27.75] cstar ~ sal + temp + aou + silicate + phosphate 20644.5396 20689.4596 2.3659579 2000-2009 –> 2010-2019 4.830956 39780.0048 39869.1884
obs Atlantic (27.5,27.75] cstar ~ sal + temp + aou + silicate + phosphate 15914.0648 15955.1238 5.1126928 2000-2009 –> 2010-2019 10.136712 29856.3756 29937.6071
mod Atlantic (27.5,27.75] cstar ~ sal + temp + aou + silicate + phosphate_star 20531.9177 20576.8378 2.3366910 2000-2009 –> 2010-2019 4.789172 39625.4433 39714.6269
obs Atlantic (27.5,27.75] cstar ~ sal + temp + aou + silicate + phosphate_star 15926.1713 15967.2303 5.1245825 2000-2009 –> 2010-2019 10.168919 29887.0147 29968.2462
mod Atlantic (27.5,27.75] cstar ~ sal + temp + nitrate + phosphate + phosphate_star 20817.2941 20862.2142 2.4115655 2000-2009 –> 2010-2019 4.906041 40050.6879 40139.8715
mod Atlantic (27.5,27.75] cstar ~ sal + temp + nitrate + silicate + phosphate_star 20597.9802 20642.9002 2.3538144 2000-2009 –> 2010-2019 4.812100 39710.9835 39800.1671
mod Atlantic (27.5,27.75] cstar ~ sal + temp + silicate + phosphate + phosphate_star 20531.3866 20576.3067 2.3365539 2000-2009 –> 2010-2019 4.788265 39622.3255 39711.5091
obs Atlantic (27.5,27.75] cstar ~ sal + temp + silicate + phosphate + phosphate_star 15915.1849 15956.2439 5.1137917 2000-2009 –> 2010-2019 10.144033 29863.1796 29944.4110
obs Atlantic (27.5,27.75] cstar ~ temp + aou + nitrate + silicate + phosphate 15881.7348 15922.7938 5.0810770 2000-2009 –> 2010-2019 10.154692 29869.1544 29950.3859
obs Atlantic (27.5,27.75] cstar ~ temp + aou + silicate + phosphate 15915.0600 15950.2534 5.1156318 2000-2009 –> 2010-2019 10.193916 29904.7035 29974.3305
obs Atlantic (27.5,27.75] cstar ~ temp + nitrate + silicate + phosphate + phosphate_star 15885.4155 15926.4745 5.0846665 2000-2009 –> 2010-2019 10.167394 29881.0756 29962.3071
obs Atlantic (27.5,27.75] cstar ~ temp + silicate + phosphate + phosphate_star 15916.4612 15951.6546 5.1170073 2000-2009 –> 2010-2019 10.204948 29914.8286 29984.4556
mod Atlantic (27.75,27.85] cstar ~ sal + aou + nitrate + silicate 5681.0497 5713.5082 1.3456478 2000-2009 –> 2010-2019 2.626581 11019.2701 11083.9876
obs Atlantic (27.75,27.85] cstar ~ sal + aou + nitrate + silicate + phosphate 10045.3712 10082.9804 5.6497864 2000-2009 –> 2010-2019 10.625660 19123.9426 19198.7397
mod Atlantic (27.75,27.85] cstar ~ sal + aou + silicate + phosphate 5413.3735 5445.8320 1.2409283 2000-2009 –> 2010-2019 2.448393 10562.8209 10627.5384
obs Atlantic (27.75,27.85] cstar ~ sal + aou + silicate + phosphate 10046.5918 10078.8283 5.6555042 2000-2009 –> 2010-2019 10.662887 19142.0878 19206.1996
obs Atlantic (27.75,27.85] cstar ~ sal + nitrate + silicate + phosphate + phosphate_star 10036.0821 10073.6913 5.6333275 2000-2009 –> 2010-2019 10.617438 19119.6120 19194.4091
mod Atlantic (27.75,27.85] cstar ~ sal + nitrate + silicate + phosphate_star 6111.4206 6143.8790 1.5328566 2000-2009 –> 2010-2019 2.947335 11766.5966 11831.3141
mod Atlantic (27.75,27.85] cstar ~ sal + silicate + phosphate + phosphate_star 6100.8862 6133.3447 1.5279771 2000-2009 –> 2010-2019 2.943227 11757.8049 11822.5224
obs Atlantic (27.75,27.85] cstar ~ sal + silicate + phosphate + phosphate_star 10036.8439 10069.0804 5.6382160 2000-2009 –> 2010-2019 10.658449 19140.0235 19204.1353
mod Atlantic (27.75,27.85] cstar ~ sal + temp + aou + nitrate + silicate 5625.1624 5663.0306 1.3222769 2000-2009 –> 2010-2019 2.603209 10965.3807 11040.8845
mod Atlantic (27.75,27.85] cstar ~ sal + temp + aou + silicate + phosphate 5268.3231 5306.1913 1.1869096 2000-2009 –> 2010-2019 2.381927 10386.6518 10462.1556
obs Atlantic (27.75,27.85] cstar ~ sal + temp + aou + silicate + phosphate 10041.9249 10079.5342 5.6436746 2000-2009 –> 2010-2019 10.651049 19139.4163 19214.2134
mod Atlantic (27.75,27.85] cstar ~ sal + temp + aou + silicate + phosphate_star 5402.0063 5439.8745 1.2359179 2000-2009 –> 2010-2019 2.469862 10622.7817 10698.2855
obs Atlantic (27.75,27.85] cstar ~ sal + temp + aou + silicate + phosphate_star 10031.1986 10068.8079 5.6246941 2000-2009 –> 2010-2019 10.646406 19137.2617 19212.0588
mod Atlantic (27.75,27.85] cstar ~ sal + temp + nitrate + silicate + phosphate_star 5409.0296 5446.8978 1.2385479 2000-2009 –> 2010-2019 2.473707 10632.9524 10708.4561
mod Atlantic (27.75,27.85] cstar ~ sal + temp + silicate + phosphate + phosphate_star 5396.2307 5434.0989 1.2337594 2000-2009 –> 2010-2019 2.465184 10610.4747 10685.9784
obs Atlantic (27.75,27.85] cstar ~ sal + temp + silicate + phosphate + phosphate_star 10038.4288 10076.0380 5.6374811 2000-2009 –> 2010-2019 10.648883 19138.3303 19213.1274
mod Atlantic (27.75,27.85] cstar ~ sal + temp + silicate + phosphate_star 5480.7926 5513.2510 1.2665099 2000-2009 –> 2010-2019 2.555165 10838.2220 10902.9395
obs Atlantic (27.75,27.85] cstar ~ temp + aou + nitrate + silicate + phosphate 10064.8357 10102.4450 5.6844307 2000-2009 –> 2010-2019 10.672513 19150.7540 19225.5511
obs Atlantic (27.75,27.85] cstar ~ temp + aou + nitrate + silicate + phosphate_star 10055.0332 10092.6424 5.6669571 2000-2009 –> 2010-2019 10.663579 19146.0795 19220.8766
obs Atlantic (27.75,27.85] cstar ~ temp + nitrate + silicate + phosphate + phosphate_star 10061.6209 10099.2302 5.6786942 2000-2009 –> 2010-2019 10.668224 19148.4091 19223.2062
mod Atlantic (27.85,27.95] cstar ~ sal + aou + nitrate + silicate 6516.7494 6549.0533 1.8242519 2000-2009 –> 2010-2019 4.920004 14833.2919 14897.9701
obs Atlantic (27.85,27.95] cstar ~ sal + aou + nitrate + silicate + phosphate 21131.2942 21174.2846 5.2372604 2000-2009 –> 2010-2019 10.706381 41524.9106 41610.5445
mod Atlantic (27.85,27.95] cstar ~ sal + aou + nitrate + silicate + phosphate_star 6516.0711 6553.7590 1.8227351 2000-2009 –> 2010-2019 4.915850 14831.8363 14907.2942
obs Atlantic (27.85,27.95] cstar ~ sal + aou + nitrate + silicate + phosphate_star 20600.3474 20643.3378 4.8476362 2000-2009 –> 2010-2019 10.010391 40617.1798 40702.8137
mod Atlantic (27.85,27.95] cstar ~ sal + aou + silicate + phosphate 6236.0480 6268.3519 1.6719584 2000-2009 –> 2010-2019 4.623774 14397.4760 14462.1542
obs Atlantic (27.85,27.95] cstar ~ sal + aou + silicate + phosphate_star 20604.7712 20641.6201 4.8521724 2000-2009 –> 2010-2019 10.014928 40619.6046 40693.0051
obs Atlantic (27.85,27.95] cstar ~ sal + nitrate + silicate + phosphate + phosphate_star 20934.8741 20977.8644 5.0895999 2000-2009 –> 2010-2019 10.440820 41186.0502 41271.6841
obs Atlantic (27.85,27.95] cstar ~ sal + silicate + phosphate + phosphate_star 20965.9911 21002.8400 5.1142009 2000-2009 –> 2010-2019 10.471493 41222.5796 41295.9801
mod Atlantic (27.85,27.95] cstar ~ sal + temp + aou + nitrate + silicate 6311.2225 6348.9104 1.7103886 2000-2009 –> 2010-2019 4.783832 14606.2016 14681.6596
mod Atlantic (27.85,27.95] cstar ~ sal + temp + aou + phosphate 6775.6684 6807.9724 1.9769982 2000-2009 –> 2010-2019 4.948825 14959.1088 15023.7870
mod Atlantic (27.85,27.95] cstar ~ sal + temp + aou + silicate + phosphate 5832.2145 5869.9025 1.4739714 2000-2009 –> 2010-2019 4.370765 13934.3401 14009.7980
obs Atlantic (27.85,27.95] cstar ~ sal + temp + aou + silicate + phosphate 20549.9436 20592.9340 4.8121899 2000-2009 –> 2010-2019 9.924443 40502.5261 40588.1600
mod Atlantic (27.85,27.95] cstar ~ sal + temp + aou + silicate + phosphate_star 5957.1357 5994.8236 1.5322784 2000-2009 –> 2010-2019 4.465872 14100.3890 14175.8469
obs Atlantic (27.85,27.95] cstar ~ sal + temp + aou + silicate + phosphate_star 20508.1814 20551.1717 4.7830172 2000-2009 –> 2010-2019 9.887030 40450.2206 40535.8545
obs Atlantic (27.85,27.95] cstar ~ sal + temp + nitrate + silicate + phosphate 20590.8830 20633.8733 4.8409605 2000-2009 –> 2010-2019 10.002462 40606.1280 40691.7619
mod Atlantic (27.85,27.95] cstar ~ sal + temp + nitrate + silicate + phosphate_star 5970.4662 6008.1542 1.5386351 2000-2009 –> 2010-2019 4.494614 14138.4853 14213.9433
obs Atlantic (27.85,27.95] cstar ~ sal + temp + silicate + phosphate 20597.4345 20634.2834 4.8469919 2000-2009 –> 2010-2019 10.013141 40616.5619 40689.9624
mod Atlantic (27.85,27.95] cstar ~ sal + temp + silicate + phosphate + phosphate_star 5945.1964 5982.8843 1.5266074 2000-2009 –> 2010-2019 4.453281 14080.7550 14156.2130
obs Atlantic (27.85,27.95] cstar ~ sal + temp + silicate + phosphate + phosphate_star 20543.6594 20586.6498 4.8077888 2000-2009 –> 2010-2019 9.918124 40493.7903 40579.4242
mod Atlantic (27.85,27.95] cstar ~ sal + temp + silicate + phosphate_star 6074.9360 6107.2399 1.5903607 2000-2009 –> 2010-2019 4.625836 14327.4165 14392.0947
obs Atlantic (27.95,28.05] cstar ~ sal + aou + nitrate + silicate + phosphate 22944.7798 22988.1972 5.5971871 2000-2009 –> 2010-2019 11.975510 46195.4142 46282.0564
obs Atlantic (27.95,28.05] cstar ~ sal + aou + nitrate + silicate + phosphate_star 22517.3942 22560.8116 5.2789025 2000-2009 –> 2010-2019 11.525464 45619.7809 45706.4232
obs Atlantic (27.95,28.05] cstar ~ sal + aou + silicate + phosphate_star 22596.1399 22633.3548 5.3376168 2000-2009 –> 2010-2019 11.602989 45717.8808 45792.1456
obs Atlantic (27.95,28.05] cstar ~ sal + nitrate + silicate + phosphate + phosphate_star 22647.4373 22690.8547 5.3737841 2000-2009 –> 2010-2019 11.645842 45778.7535 45865.3958
obs Atlantic (27.95,28.05] cstar ~ sal + silicate + phosphate + phosphate_star 22749.7763 22786.9912 5.4511432 2000-2009 –> 2010-2019 11.738212 45896.0686 45970.3334
mod Atlantic (27.95,28.05] cstar ~ sal + temp + aou + phosphate 10153.2795 10186.6431 3.3893728 2000-2009 –> 2010-2019 7.183448 20244.4926 20310.9319
mod Atlantic (27.95,28.05] cstar ~ sal + temp + aou + phosphate_star 10185.9010 10219.2646 3.4182736 2000-2009 –> 2010-2019 7.223558 20287.9168 20354.3561
mod Atlantic (27.95,28.05] cstar ~ sal + temp + aou + silicate + phosphate 10017.9201 10056.8443 3.2703363 2000-2009 –> 2010-2019 7.027954 20075.7740 20153.2866
obs Atlantic (27.95,28.05] cstar ~ sal + temp + aou + silicate + phosphate 22665.7774 22709.1948 5.3873018 2000-2009 –> 2010-2019 11.677915 45818.0725 45904.7148
mod Atlantic (27.95,28.05] cstar ~ sal + temp + aou + silicate + phosphate_star 10031.1733 10070.0975 3.2816370 2000-2009 –> 2010-2019 7.041501 20091.2158 20168.7284
obs Atlantic (27.95,28.05] cstar ~ sal + temp + aou + silicate + phosphate_star 22587.8357 22631.2531 5.3300879 2000-2009 –> 2010-2019 11.588862 45704.0939 45790.7362
mod Atlantic (27.95,28.05] cstar ~ sal + temp + nitrate + phosphate_star 10248.5724 10281.9360 3.4744904 2000-2009 –> 2010-2019 7.344959 20412.7871 20479.2264
obs Atlantic (27.95,28.05] cstar ~ sal + temp + nitrate + silicate + phosphate 22787.9662 22831.3836 5.4782345 2000-2009 –> 2010-2019 11.994181 46190.2086 46276.8509
mod Atlantic (27.95,28.05] cstar ~ sal + temp + nitrate + silicate + phosphate_star 10007.0121 10045.9363 3.2610645 2000-2009 –> 2010-2019 7.019482 20065.6458 20143.1584
obs Atlantic (27.95,28.05] cstar ~ sal + temp + nitrate + silicate + phosphate_star 22840.9461 22884.3635 5.5181375 2000-2009 –> 2010-2019 11.990337 46195.3460 46281.9883
mod Atlantic (27.95,28.05] cstar ~ sal + temp + phosphate + phosphate_star 10181.2312 10214.5949 3.4141215 2000-2009 –> 2010-2019 7.217478 20281.3914 20347.8307
mod Atlantic (27.95,28.05] cstar ~ sal + temp + phosphate_star 10260.3228 10288.1258 3.4869478 2000-2009 –> 2010-2019 7.375544 20439.6482 20495.0143
mod Atlantic (27.95,28.05] cstar ~ sal + temp + silicate + phosphate + phosphate_star 10031.3576 10070.2818 3.2817945 2000-2009 –> 2010-2019 7.041492 20091.2380 20168.7506
obs Atlantic (27.95,28.05] cstar ~ sal + temp + silicate + phosphate + phosphate_star 22637.9478 22681.3652 5.3668031 2000-2009 –> 2010-2019 11.642604 45773.5007 45860.1430
mod Atlantic (27.95,28.05] cstar ~ sal + temp + silicate + phosphate_star 10029.6675 10063.0311 3.2820592 2000-2009 –> 2010-2019 7.047352 20092.9940 20159.4333
obs Atlantic (28.05,28.1] cstar ~ sal + aou + nitrate + silicate + phosphate_star 18842.1482 18884.2254 5.4986911 2000-2009 –> 2010-2019 12.211342 36547.4406 36630.7256
obs Atlantic (28.05,28.1] cstar ~ sal + aou + silicate + phosphate_star 18907.4820 18943.5482 5.5604567 2000-2009 –> 2010-2019 12.281318 36617.2822 36688.6693
obs Atlantic (28.05,28.1] cstar ~ sal + nitrate + silicate + phosphate + phosphate_star 18858.8153 18900.8925 5.5139157 2000-2009 –> 2010-2019 12.232525 36568.8313 36652.1163
obs Atlantic (28.05,28.1] cstar ~ sal + temp + aou + nitrate + phosphate_star 18907.4362 18949.5134 5.5585699 2000-2009 –> 2010-2019 12.244577 36591.5549 36674.8399
mod Atlantic (28.05,28.1] cstar ~ sal + temp + aou + phosphate 3618.5704 3648.5997 1.2428864 2000-2009 –> 2010-2019 2.313359 6628.4397 6687.9634
mod Atlantic (28.05,28.1] cstar ~ sal + temp + aou + phosphate_star 3661.7943 3691.8236 1.2675019 2000-2009 –> 2010-2019 2.362004 6716.4188 6775.9424
obs Atlantic (28.05,28.1] cstar ~ sal + temp + aou + phosphate_star 18936.7037 18972.7698 5.5874774 2000-2009 –> 2010-2019 12.275249 36620.2271 36691.6142
mod Atlantic (28.05,28.1] cstar ~ sal + temp + aou + silicate + phosphate 3619.4120 3654.4462 1.2422333 2000-2009 –> 2010-2019 2.312696 6631.2627 6700.7070
mod Atlantic (28.05,28.1] cstar ~ sal + temp + aou + silicate + phosphate_star 3659.9425 3694.9767 1.2652888 2000-2009 –> 2010-2019 2.359524 6716.0747 6785.5189
obs Atlantic (28.05,28.1] cstar ~ sal + temp + aou + silicate + phosphate_star 18876.8715 18918.9487 5.5304568 2000-2009 –> 2010-2019 12.199857 36547.7495 36631.0345
obs Atlantic (28.05,28.1] cstar ~ sal + temp + nitrate + phosphate + phosphate_star 18888.1775 18930.2546 5.5408393 2000-2009 –> 2010-2019 12.230909 36575.5303 36658.8153
mod Atlantic (28.05,28.1] cstar ~ sal + temp + nitrate + phosphate_star 3707.4623 3737.4916 1.2940392 2000-2009 –> 2010-2019 2.423692 6825.8137 6885.3374
obs Atlantic (28.05,28.1] cstar ~ sal + temp + nitrate + phosphate_star 18912.0131 18948.0793 5.5646379 2000-2009 –> 2010-2019 12.254804 36597.4426 36668.8297
mod Atlantic (28.05,28.1] cstar ~ sal + temp + nitrate + silicate + phosphate_star 3679.6070 3714.6411 1.2766284 2000-2009 –> 2010-2019 2.394138 6778.1700 6847.6142
obs Atlantic (28.05,28.1] cstar ~ sal + temp + nitrate + silicate + phosphate_star 18836.4066 18878.4837 5.4934561 2000-2009 –> 2010-2019 12.176882 36518.4692 36601.7542
mod Atlantic (28.05,28.1] cstar ~ sal + temp + phosphate + phosphate_star 3650.1261 3680.1554 1.2608094 2000-2009 –> 2010-2019 2.348982 6693.0574 6752.5810
mod Atlantic (28.05,28.1] cstar ~ sal + temp + silicate + phosphate + phosphate_star 3649.9721 3685.0063 1.2595778 2000-2009 –> 2010-2019 2.347700 6694.8103 6764.2545
obs Atlantic (28.05,28.1] cstar ~ sal + temp + silicate + phosphate + phosphate_star 18926.0902 18968.1674 5.5757979 2000-2009 –> 2010-2019 12.268861 36615.8245 36699.1095
mod Atlantic (28.05,28.1] cstar ~ temp + aou + phosphate 3662.6997 3687.7241 1.2691739 2000-2009 –> 2010-2019 2.404336 6788.8593 6838.4623
mod Atlantic (28.05,28.1] cstar ~ temp + aou + silicate + phosphate 3663.8197 3693.8490 1.2686673 2000-2009 –> 2010-2019 2.401988 6788.7064 6848.2300
obs Atlantic (28.1,28.15] cstar ~ sal + temp + aou + nitrate + phosphate_star 15175.3187 15215.6836 6.0086810 2000-2009 –> 2010-2019 13.444738 28993.1944 29072.8180
mod Atlantic (28.1,28.15] cstar ~ sal + temp + aou + phosphate 4814.7717 4847.2808 1.0227239 2000-2009 –> 2010-2019 2.322530 10361.1289 10426.0746
mod Atlantic (28.1,28.15] cstar ~ sal + temp + aou + phosphate_star 4864.0541 4896.5632 1.0379630 2000-2009 –> 2010-2019 2.346889 10433.4293 10498.3751
mod Atlantic (28.1,28.15] cstar ~ sal + temp + aou + silicate + phosphate 4758.7756 4796.7028 1.0050766 2000-2009 –> 2010-2019 2.229304 10109.9246 10185.6946
mod Atlantic (28.1,28.15] cstar ~ sal + temp + aou + silicate + phosphate_star 4822.0279 4859.9552 1.0243385 2000-2009 –> 2010-2019 2.263054 10211.9094 10287.6794
obs Atlantic (28.1,28.15] cstar ~ sal + temp + aou + silicate + phosphate_star 15181.6654 15222.0303 6.0167660 2000-2009 –> 2010-2019 13.446412 28996.0647 29075.6883
obs Atlantic (28.1,28.15] cstar ~ sal + temp + nitrate + phosphate + phosphate_star 15155.4473 15195.8122 5.9834373 2000-2009 –> 2010-2019 13.414872 28970.8166 29050.4401
mod Atlantic (28.1,28.15] cstar ~ sal + temp + nitrate + phosphate_star 4866.3577 4898.8668 1.0386809 2000-2009 –> 2010-2019 2.347368 10435.1309 10500.0766
mod Atlantic (28.1,28.15] cstar ~ sal + temp + nitrate + silicate + phosphate_star 4867.7814 4905.7087 1.0385013 2000-2009 –> 2010-2019 2.335508 10409.0413 10484.8113
obs Atlantic (28.1,28.15] cstar ~ sal + temp + nitrate + silicate + phosphate_star 15170.2232 15210.5881 6.0021978 2000-2009 –> 2010-2019 13.427964 28982.5176 29062.1411
mod Atlantic (28.1,28.15] cstar ~ sal + temp + phosphate + phosphate_star 4849.8595 4882.3686 1.0335506 2000-2009 –> 2010-2019 2.337062 10405.5871 10470.5328
obs Atlantic (28.1,28.15] cstar ~ sal + temp + phosphate + phosphate_star 15183.4490 15218.0475 6.0215910 2000-2009 –> 2010-2019 13.457889 28999.4552 29067.7040
mod Atlantic (28.1,28.15] cstar ~ sal + temp + silicate + phosphate + phosphate_star 4800.0509 4837.9781 1.0176044 2000-2009 –> 2010-2019 2.246741 10164.3745 10240.1445
obs Atlantic (28.1,28.15] cstar ~ sal + temp + silicate + phosphate + phosphate_star 15181.1895 15221.5545 6.0161594 2000-2009 –> 2010-2019 13.447015 28996.2450 29075.8685
obs Atlantic (28.1,28.15] cstar ~ sal + temp + silicate + phosphate_star 15179.9208 15214.5193 6.0170916 2000-2009 –> 2010-2019 13.448162 28993.0930 29061.3418
obs Atlantic (28.1,28.15] cstar ~ temp + aou + nitrate + silicate + phosphate_star 15172.5732 15212.9381 6.0051869 2000-2009 –> 2010-2019 13.452687 28996.6454 29076.2690
mod Atlantic (28.1,28.15] cstar ~ temp + aou + silicate + phosphate 4918.6434 4951.1525 1.0551084 2000-2009 –> 2010-2019 2.342361 10433.0535 10497.9993
obs Atlantic (28.1,28.15] cstar ~ temp + nitrate + phosphate + phosphate_star 15168.3484 15202.9470 6.0023571 2000-2009 –> 2010-2019 13.459524 28995.6485 29063.8972
obs Atlantic (28.1,28.15] cstar ~ temp + nitrate + silicate + phosphate + phosphate_star 15159.3819 15199.7468 5.9884273 2000-2009 –> 2010-2019 13.434920 28982.9093 29062.5328
mod Atlantic (28.1,28.15] cstar ~ temp + silicate + phosphate + phosphate_star 4973.3226 5005.8317 1.0725659 2000-2009 –> 2010-2019 2.367205 10506.5677 10571.5134
mod Atlantic (28.15,28.2] cstar ~ sal + aou + nitrate + silicate + phosphate_star 7387.5577 7429.7363 0.8078969 2000-2009 –> 2010-2019 1.820270 15427.9289 15511.6791
mod Atlantic (28.15,28.2] cstar ~ sal + temp + aou + nitrate 7416.7507 7452.9038 0.8120279 2000-2009 –> 2010-2019 1.810954 15380.1299 15451.9159
obs Atlantic (28.15,28.2] cstar ~ sal + temp + aou + nitrate + phosphate 4819.3564 4852.8396 3.6770151 2000-2009 –> 2010-2019 7.322842 9797.2710 9864.4867
mod Atlantic (28.15,28.2] cstar ~ sal + temp + aou + nitrate + silicate 7313.2217 7355.4003 0.7981368 2000-2009 –> 2010-2019 1.790716 15242.8560 15326.6063
obs Atlantic (28.15,28.2] cstar ~ sal + temp + aou + silicate + phosphate 4746.4253 4779.9086 3.5282570 2000-2009 –> 2010-2019 7.166982 9720.7720 9787.9877
obs Atlantic (28.15,28.2] cstar ~ sal + temp + aou + silicate + phosphate_star 4739.4902 4772.9734 3.5144286 2000-2009 –> 2010-2019 7.186424 9730.4932 9797.7089
mod Atlantic (28.15,28.2] cstar ~ sal + temp + nitrate + phosphate_star 7828.0650 7864.2181 0.8685167 2000-2009 –> 2010-2019 1.894535 15941.5141 16013.3000
mod Atlantic (28.15,28.2] cstar ~ sal + temp + nitrate + silicate + phosphate_star 7761.9732 7804.1518 0.8589008 2000-2009 –> 2010-2019 1.878772 15843.7277 15927.4779
obs Atlantic (28.15,28.2] cstar ~ sal + temp + silicate + phosphate + phosphate_star 4740.0384 4773.5217 3.5155198 2000-2009 –> 2010-2019 7.178246 9726.4162 9793.6319
mod Atlantic (28.15,28.2] cstar ~ temp + aou + nitrate + phosphate 7429.2680 7465.4211 0.8136915 2000-2009 –> 2010-2019 1.810464 15380.5482 15452.3341
mod Atlantic (28.15,28.2] cstar ~ temp + aou + nitrate + phosphate_star 7466.5851 7502.7382 0.8186715 2000-2009 –> 2010-2019 1.820403 15445.6926 15517.4785
mod Atlantic (28.15,28.2] cstar ~ temp + aou + nitrate + silicate 7359.8928 7396.0459 0.8045138 2000-2009 –> 2010-2019 1.804204 15327.5611 15399.3470
obs Atlantic (28.15,28.2] cstar ~ temp + aou + nitrate + silicate + phosphate 4759.3630 4792.8463 3.5541999 2000-2009 –> 2010-2019 7.198671 9736.5970 9803.8128
obs Atlantic (28.15,28.2] cstar ~ temp + aou + nitrate + silicate + phosphate_star 4761.0059 4794.4892 3.5575079 2000-2009 –> 2010-2019 7.240404 9757.4332 9824.6490
obs Atlantic (28.15,28.2] cstar ~ temp + aou + silicate + phosphate 4761.5959 4790.2958 3.5627291 2000-2009 –> 2010-2019 7.207294 9736.8769 9794.4904
obs Atlantic (28.15,28.2] cstar ~ temp + aou + silicate + phosphate_star 4763.9552 4792.6551 3.5674920 2000-2009 –> 2010-2019 7.250429 9758.4028 9816.0163
mod Atlantic (28.15,28.2] cstar ~ temp + nitrate + phosphate + phosphate_star 7525.5453 7561.6984 0.8266019 2000-2009 –> 2010-2019 1.834227 15537.5531 15609.3390
obs Atlantic (28.15,28.2] cstar ~ temp + nitrate + silicate + phosphate + phosphate_star 4758.3376 4791.8209 3.5521368 2000-2009 –> 2010-2019 7.224422 9749.4848 9816.7006
mod Atlantic (28.15,28.2] cstar ~ temp + nitrate + silicate + phosphate_star 7779.5121 7815.6652 0.8616491 2000-2009 –> 2010-2019 1.882483 15864.5569 15936.3428
obs Atlantic (28.15,28.2] cstar ~ temp + silicate + phosphate + phosphate_star 4760.6039 4789.3038 3.5607284 2000-2009 –> 2010-2019 7.233404 9749.9457 9807.5592
mod Atlantic (28.2, Inf] cstar ~ sal + aou + nitrate 9802.2820 9835.1640 0.6089561 2000-2009 –> 2010-2019 1.297371 19117.8937 19182.7790
mod Atlantic (28.2, Inf] cstar ~ sal + aou + nitrate + phosphate 9720.5262 9759.9847 0.6041679 2000-2009 –> 2010-2019 1.289535 18998.6524 19076.5148
obs Atlantic (28.2, Inf] cstar ~ sal + aou + nitrate + phosphate + phosphate_star 11026.5194 11066.0506 3.3512980 2000-2009 –> 2010-2019 6.583887 22109.3065 22188.5013
mod Atlantic (28.2, Inf] cstar ~ sal + aou + nitrate + phosphate_star 9640.6785 9680.1370 0.5996382 2000-2009 –> 2010-2019 1.285333 18923.0592 19000.9216
mod Atlantic (28.2, Inf] cstar ~ sal + aou + nitrate + silicate 9639.1439 9678.6023 0.5995515 2000-2009 –> 2010-2019 1.267143 18683.3967 18761.2591
mod Atlantic (28.2, Inf] cstar ~ sal + aou + nitrate + silicate + phosphate_star 8876.1085 8922.1434 0.5578425 2000-2009 –> 2010-2019 1.208781 17697.5388 17788.3783
obs Atlantic (28.2, Inf] cstar ~ sal + aou + phosphate + phosphate_star 11024.7221 11058.6059 3.3514600 2000-2009 –> 2010-2019 6.592584 22116.7693 22184.6504
obs Atlantic (28.2, Inf] cstar ~ sal + aou + silicate + phosphate + phosphate_star 10982.7820 11022.3132 3.3164974 2000-2009 –> 2010-2019 6.557457 22076.6125 22155.8072
obs Atlantic (28.2, Inf] cstar ~ sal + nitrate + silicate + phosphate + phosphate_star 11047.9103 11087.4415 3.3684508 2000-2009 –> 2010-2019 6.601239 22130.9605 22210.1552
obs Atlantic (28.2, Inf] cstar ~ sal + silicate + phosphate + phosphate_star 11052.1730 11086.0569 3.3734894 2000-2009 –> 2010-2019 6.614841 22144.5196 22212.4008
mod Atlantic (28.2, Inf] cstar ~ sal + temp + aou + nitrate 9727.6618 9767.1202 0.6045744 2000-2009 –> 2010-2019 1.291290 19023.2852 19101.1476
mod Atlantic (28.2, Inf] cstar ~ sal + temp + aou + nitrate + phosphate 9646.7354 9692.7702 0.5998675 2000-2009 –> 2010-2019 1.285089 18924.9735 19015.8130
mod Atlantic (28.2, Inf] cstar ~ sal + temp + aou + nitrate + phosphate_star 9641.3472 9687.3820 0.5995630 2000-2009 –> 2010-2019 1.285116 18923.8823 19014.7217
mod Atlantic (28.2, Inf] cstar ~ sal + temp + aou + nitrate + silicate 9564.2174 9610.2523 0.5952202 2000-2009 –> 2010-2019 1.259352 18564.2311 18655.0705
obs Atlantic (28.2, Inf] cstar ~ sal + temp + aou + silicate + phosphate 10985.7636 11025.2947 3.3188581 2000-2009 –> 2010-2019 6.593545 22123.8000 22202.9947
obs Atlantic (28.2, Inf] cstar ~ sal + temp + aou + silicate + phosphate_star 10995.6546 11035.1858 3.3267020 2000-2009 –> 2010-2019 6.545566 22060.2736 22139.4683
mod Atlantic (28.2, Inf] cstar ~ sal + temp + nitrate + phosphate + phosphate_star 9690.3014 9736.3362 0.6023357 2000-2009 –> 2010-2019 1.291451 19018.9698 19109.8092
obs Atlantic (28.2, Inf] cstar ~ sal + temp + nitrate + phosphate + phosphate_star 11053.5698 11093.1010 3.3730037 2000-2009 –> 2010-2019 6.600259 22129.3057 22208.5004
obs Atlantic (28.2, Inf] cstar ~ sal + temp + phosphate + phosphate_star 11054.7757 11088.6595 3.3755855 2000-2009 –> 2010-2019 6.610218 22138.2617 22206.1429
obs Atlantic (28.2, Inf] cstar ~ sal + temp + silicate + phosphate + phosphate_star 10977.1658 11016.6969 3.3120549 2000-2009 –> 2010-2019 6.546640 22062.5887 22141.7834
mod Indo-Pacific (-Inf,26] cstar ~ sal + aou + nitrate + phosphate + phosphate_star 23935.4904 23979.2993 5.3639472 2000-2009 –> 2010-2019 11.130849 41542.6062 41627.9115
obs Indo-Pacific (-Inf,26] cstar ~ sal + aou + nitrate + phosphate + phosphate_star 31747.9561 31793.4891 6.0152701 2000-2009 –> 2010-2019 12.816255 55761.1830 55850.0310
mod Indo-Pacific (-Inf,26] cstar ~ sal + aou + phosphate + phosphate_star 24999.9587 25037.5092 6.1585643 2000-2009 –> 2010-2019 12.246967 42906.0580 42979.1768
mod Indo-Pacific (-Inf,26] cstar ~ sal + aou + silicate + phosphate + phosphate_star 24793.6201 24837.4291 5.9945865 2000-2009 –> 2010-2019 11.751389 42391.0112 42476.3165
obs Indo-Pacific (-Inf,26] cstar ~ sal + aou + silicate + phosphate + phosphate_star 32617.6059 32663.1389 6.5689771 2000-2009 –> 2010-2019 13.511879 56779.4062 56868.2542
mod Indo-Pacific (-Inf,26] cstar ~ sal + temp + aou + nitrate 23815.6775 23853.2280 5.2827108 2000-2009 –> 2010-2019 11.978078 42249.0036 42322.1224
obs Indo-Pacific (-Inf,26] cstar ~ sal + temp + aou + nitrate 31936.0533 31975.0816 6.1321759 2000-2009 –> 2010-2019 13.070481 56091.0886 56167.2440
mod Indo-Pacific (-Inf,26] cstar ~ sal + temp + aou + nitrate + phosphate 23748.0287 23791.8376 5.2352653 2000-2009 –> 2010-2019 11.258468 41596.3943 41681.6995
obs Indo-Pacific (-Inf,26] cstar ~ sal + temp + aou + nitrate + phosphate 31820.6490 31866.1820 6.0597092 2000-2009 –> 2010-2019 12.868941 55842.5936 55931.4416
mod Indo-Pacific (-Inf,26] cstar ~ sal + temp + aou + nitrate + phosphate_star 23749.2069 23793.0159 5.2360643 2000-2009 –> 2010-2019 11.346237 41677.1087 41762.4140
obs Indo-Pacific (-Inf,26] cstar ~ sal + temp + aou + nitrate + phosphate_star 31832.5780 31878.1110 6.0670330 2000-2009 –> 2010-2019 12.883771 55862.4496 55951.2976
mod Indo-Pacific (-Inf,26] cstar ~ sal + temp + aou + nitrate + silicate 23649.3703 23693.1792 5.1687864 2000-2009 –> 2010-2019 11.328812 41622.3539 41707.6592
obs Indo-Pacific (-Inf,26] cstar ~ sal + temp + aou + nitrate + silicate 31912.5602 31958.0932 6.1163672 2000-2009 –> 2010-2019 13.032362 56046.4255 56135.2735
mod Indo-Pacific (-Inf,26] cstar ~ sal + temp + aou + silicate + phosphate 25182.2113 25226.0203 6.3040511 2000-2009 –> 2010-2019 12.469620 43160.1856 43245.4909
obs Indo-Pacific (-Inf,26] cstar ~ sal + temp + aou + silicate + phosphate 32739.8974 32785.4304 6.6508245 2000-2009 –> 2010-2019 13.681366 56991.9380 57080.7860
mod Indo-Pacific (-Inf,26] cstar ~ sal + temp + aou + silicate + phosphate_star 25284.3975 25328.2065 6.3880498 2000-2009 –> 2010-2019 12.703933 43396.0077 43481.3129
obs Indo-Pacific (-Inf,26] cstar ~ sal + temp + aou + silicate + phosphate_star 32782.7509 32828.2839 6.6797462 2000-2009 –> 2010-2019 13.743080 57068.2685 57157.1165
mod Indo-Pacific (-Inf,26] cstar ~ sal + temp + nitrate + phosphate + phosphate_star 23781.3797 23825.1886 5.2579309 2000-2009 –> 2010-2019 11.456148 41788.6548 41873.9601
obs Indo-Pacific (-Inf,26] cstar ~ sal + temp + nitrate + phosphate + phosphate_star 31905.1149 31950.6479 6.1117580 2000-2009 –> 2010-2019 12.966714 55975.2068 56064.0548
obs Indo-Pacific (-Inf,26] cstar ~ sal + temp + silicate + phosphate + phosphate_star 32804.1550 32849.6880 6.6942388 2000-2009 –> 2010-2019 13.775240 57107.6443 57196.4923
mod Indo-Pacific (26,26.5] cstar ~ aou + nitrate + silicate + phosphate + phosphate_star 16506.4999 16548.1682 4.4000093 2000-2009 –> 2010-2019 9.507247 33012.7220 33095.7077
mod Indo-Pacific (26,26.5] cstar ~ sal + aou + nitrate + phosphate + phosphate_star 16603.0132 16644.6815 4.4753318 2000-2009 –> 2010-2019 9.503477 33024.8291 33107.8149
mod Indo-Pacific (26,26.5] cstar ~ sal + aou + silicate + phosphate + phosphate_star 16192.0123 16233.6806 4.1632563 2000-2009 –> 2010-2019 8.865669 32251.6250 32334.6107
mod Indo-Pacific (26,26.5] cstar ~ sal + temp + aou + nitrate + phosphate 16838.2080 16879.8763 4.6643307 2000-2009 –> 2010-2019 9.635312 33198.1898 33281.1755
obs Indo-Pacific (26,26.5] cstar ~ sal + temp + aou + nitrate + phosphate 17101.4450 17143.0365 5.0499059 2000-2009 –> 2010-2019 10.320166 33956.2913 34039.2750
obs Indo-Pacific (26,26.5] cstar ~ sal + temp + aou + nitrate + phosphate_star 17111.0313 17152.6229 5.0585210 2000-2009 –> 2010-2019 10.325433 33962.4036 34045.3872
mod Indo-Pacific (26,26.5] cstar ~ sal + temp + aou + silicate + phosphate 16437.6270 16479.2953 4.3470347 2000-2009 –> 2010-2019 8.958590 32391.7273 32474.7130
obs Indo-Pacific (26,26.5] cstar ~ sal + temp + aou + silicate + phosphate 17256.9646 17298.5562 5.1914989 2000-2009 –> 2010-2019 10.472046 34122.4695 34205.4531
mod Indo-Pacific (26,26.5] cstar ~ sal + temp + aou + silicate + phosphate_star 16566.3088 16607.9771 4.4465356 2000-2009 –> 2010-2019 9.154089 32631.8301 32714.8158
obs Indo-Pacific (26,26.5] cstar ~ sal + temp + aou + silicate + phosphate_star 17276.0943 17317.6858 5.2091875 2000-2009 –> 2010-2019 10.511903 34164.4976 34247.4813
obs Indo-Pacific (26,26.5] cstar ~ sal + temp + nitrate + phosphate + phosphate_star 17114.0239 17155.6154 5.0612134 2000-2009 –> 2010-2019 10.329051 33966.3567 34049.3403
obs Indo-Pacific (26,26.5] cstar ~ sal + temp + nitrate + phosphate_star 17232.2548 17267.9047 5.1705778 2000-2009 –> 2010-2019 10.616510 34264.3261 34335.4550
obs Indo-Pacific (26,26.5] cstar ~ sal + temp + nitrate + silicate + phosphate_star 17211.5675 17253.1590 5.1497616 2000-2009 –> 2010-2019 10.489981 34138.4933 34221.4770
obs Indo-Pacific (26,26.5] cstar ~ sal + temp + phosphate + phosphate_star 17370.5288 17406.1787 5.2992794 2000-2009 –> 2010-2019 10.852135 34508.8774 34580.0063
mod Indo-Pacific (26,26.5] cstar ~ sal + temp + silicate + phosphate + phosphate_star 16596.1374 16637.8057 4.4699233 2000-2009 –> 2010-2019 9.200439 32687.9744 32770.9602
obs Indo-Pacific (26,26.5] cstar ~ sal + temp + silicate + phosphate + phosphate_star 17276.1310 17317.7226 5.2092215 2000-2009 –> 2010-2019 10.513725 34166.3768 34249.3605
obs Indo-Pacific (26,26.5] cstar ~ sal + temp + silicate + phosphate_star 17274.1393 17309.7892 5.2092292 2000-2009 –> 2010-2019 10.556247 34206.0194 34277.1482
mod Indo-Pacific (26,26.5] cstar ~ temp + aou + nitrate + silicate + phosphate 16569.5291 16611.1975 4.4490547 2000-2009 –> 2010-2019 9.380665 32886.5080 32969.4937
mod Indo-Pacific (26,26.5] cstar ~ temp + aou + nitrate + silicate + phosphate_star 16633.6709 16675.3392 4.4995270 2000-2009 –> 2010-2019 9.409142 32926.4764 33009.4622
mod Indo-Pacific (26,26.5] cstar ~ temp + nitrate + silicate + phosphate + phosphate_star 16670.7873 16712.4557 4.5289946 2000-2009 –> 2010-2019 9.444098 32969.6344 33052.6201
mod Indo-Pacific (26.5,26.75] cstar ~ aou + nitrate + silicate + phosphate + phosphate_star 15129.0822 15170.1170 4.4422716 2000-2009 –> 2010-2019 9.364191 29772.4477 29854.0463
mod Indo-Pacific (26.5,26.75] cstar ~ aou + silicate + phosphate + phosphate_star 15228.4086 15263.5813 4.5297838 2000-2009 –> 2010-2019 9.623768 30036.6352 30106.5768
mod Indo-Pacific (26.5,26.75] cstar ~ sal + temp + aou + nitrate + phosphate 15389.6013 15430.6361 4.6707682 2000-2009 –> 2010-2019 9.734464 30170.8684 30252.4670
obs Indo-Pacific (26.5,26.75] cstar ~ sal + temp + aou + nitrate + phosphate 16988.4144 17029.8372 5.3278098 2000-2009 –> 2010-2019 10.873881 34973.9924 35057.1472
obs Indo-Pacific (26.5,26.75] cstar ~ sal + temp + aou + nitrate + phosphate_star 17025.4142 17066.8369 5.3638377 2000-2009 –> 2010-2019 10.955930 35058.4101 35141.5649
obs Indo-Pacific (26.5,26.75] cstar ~ sal + temp + aou + phosphate 17294.2447 17329.7499 5.6350799 2000-2009 –> 2010-2019 11.448819 35548.2786 35619.5542
mod Indo-Pacific (26.5,26.75] cstar ~ sal + temp + aou + silicate + phosphate 15037.3047 15078.3394 4.3644664 2000-2009 –> 2010-2019 9.187789 29582.4063 29664.0049
obs Indo-Pacific (26.5,26.75] cstar ~ sal + temp + aou + silicate + phosphate 17290.0762 17331.4989 5.6287519 2000-2009 –> 2010-2019 11.429783 35533.5547 35616.7095
mod Indo-Pacific (26.5,26.75] cstar ~ sal + temp + aou + silicate + phosphate_star 15158.6072 15199.6419 4.4675953 2000-2009 –> 2010-2019 9.386077 29798.5802 29880.1787
obs Indo-Pacific (26.5,26.75] cstar ~ sal + temp + aou + silicate + phosphate_star 17326.7083 17368.1311 5.6664354 2000-2009 –> 2010-2019 11.494999 35597.3549 35680.5096
obs Indo-Pacific (26.5,26.75] cstar ~ sal + temp + nitrate + phosphate + phosphate_star 17031.2941 17072.7168 5.3695856 2000-2009 –> 2010-2019 10.967484 35070.2447 35153.3995
obs Indo-Pacific (26.5,26.75] cstar ~ sal + temp + nitrate + phosphate_star 17092.5773 17128.0825 5.4318389 2000-2009 –> 2010-2019 11.083485 35184.3580 35255.6335
mod Indo-Pacific (26.5,26.75] cstar ~ sal + temp + nitrate + silicate + phosphate_star 15313.1606 15354.1954 4.6025314 2000-2009 –> 2010-2019 9.671050 30099.0500 30180.6485
obs Indo-Pacific (26.5,26.75] cstar ~ sal + temp + nitrate + silicate + phosphate_star 17055.6022 17097.0250 5.3934133 2000-2009 –> 2010-2019 10.989538 35092.7340 35175.8888
obs Indo-Pacific (26.5,26.75] cstar ~ sal + temp + phosphate + phosphate_star 17329.2132 17364.7184 5.6710870 2000-2009 –> 2010-2019 11.518309 35616.1993 35687.4749
mod Indo-Pacific (26.5,26.75] cstar ~ sal + temp + silicate + phosphate + phosphate_star 15168.4129 15209.4477 4.4760376 2000-2009 –> 2010-2019 9.402604 29816.3610 29897.9595
obs Indo-Pacific (26.5,26.75] cstar ~ sal + temp + silicate + phosphate + phosphate_star 17319.7801 17361.2028 5.6592890 2000-2009 –> 2010-2019 11.483613 35586.2520 35669.4068
mod Indo-Pacific (26.5,26.75] cstar ~ sal + temp + silicate + phosphate_star 15313.9278 15349.1005 4.6049841 2000-2009 –> 2010-2019 9.676911 30101.0818 30171.0234
mod Indo-Pacific (26.5,26.75] cstar ~ temp + aou + nitrate + silicate + phosphate 15282.5805 15323.6153 4.5755131 2000-2009 –> 2010-2019 9.585397 30011.9660 30093.5646
mod Indo-Pacific (26.5,26.75] cstar ~ temp + aou + nitrate + silicate + phosphate_star 15386.8518 15427.8866 4.6682963 2000-2009 –> 2010-2019 9.763877 30198.5994 30280.1979
obs Indo-Pacific (26.75,27] cstar ~ aou + nitrate + phosphate + phosphate_star 25312.7051 25350.6234 5.2777826 2000-2009 –> 2010-2019 10.563443 52161.1164 52237.3037
obs Indo-Pacific (26.75,27] cstar ~ aou + nitrate + silicate + phosphate 25526.2906 25564.2089 5.4169215 2000-2009 –> 2010-2019 10.886868 52672.9309 52749.1182
obs Indo-Pacific (26.75,27] cstar ~ sal + aou + nitrate + silicate + phosphate 25496.3054 25540.5434 5.3958537 2000-2009 –> 2010-2019 10.826607 52582.3706 52671.2557
mod Indo-Pacific (26.75,27] cstar ~ sal + temp + aou + nitrate + phosphate 19205.7438 19248.7117 3.9924811 2000-2009 –> 2010-2019 8.150866 39708.3902 39794.6828
obs Indo-Pacific (26.75,27] cstar ~ sal + temp + aou + nitrate + phosphate 25446.9860 25491.2240 5.3635289 2000-2009 –> 2010-2019 10.758454 52475.4507 52564.3358
mod Indo-Pacific (26.75,27] cstar ~ sal + temp + aou + nitrate + phosphate_star 19365.1420 19408.1099 4.0865303 2000-2009 –> 2010-2019 8.315354 39988.7939 40075.0866
mod Indo-Pacific (26.75,27] cstar ~ sal + temp + nitrate + phosphate + phosphate_star 19376.5458 19419.5137 4.0933432 2000-2009 –> 2010-2019 8.324676 40004.4717 40090.7643
mod Indo-Pacific (26.75,27] cstar ~ sal + temp + nitrate + phosphate_star 19477.6965 19514.5262 4.1554858 2000-2009 –> 2010-2019 8.397761 40122.2276 40196.1927
mod Indo-Pacific (26.75,27] cstar ~ sal + temp + nitrate + silicate + phosphate_star 19384.2499 19427.2178 4.0979522 2000-2009 –> 2010-2019 8.319758 39995.9375 40082.2301
obs Indo-Pacific (26.75,27] cstar ~ sal + temp + nitrate + silicate + phosphate_star 25407.1337 25451.3718 5.3375505 2000-2009 –> 2010-2019 10.691435 52369.1473 52458.0324
mod Indo-Pacific (26.75,27] cstar ~ temp + aou + nitrate + phosphate 19658.5468 19695.3765 4.2667238 2000-2009 –> 2010-2019 8.662215 40558.6726 40632.6377
obs Indo-Pacific (26.75,27] cstar ~ temp + aou + nitrate + phosphate 25461.8357 25499.7540 5.3745507 2000-2009 –> 2010-2019 10.792090 52524.7007 52600.8880
mod Indo-Pacific (26.75,27] cstar ~ temp + aou + nitrate + silicate + phosphate 19427.5155 19470.4834 4.1239327 2000-2009 –> 2010-2019 8.466286 40242.0205 40328.3132
obs Indo-Pacific (26.75,27] cstar ~ temp + aou + nitrate + silicate + phosphate 25347.1538 25391.3918 5.2986885 2000-2009 –> 2010-2019 10.615356 52248.4647 52337.3499
mod Indo-Pacific (26.75,27] cstar ~ temp + aou + nitrate + silicate + phosphate_star 19566.5435 19609.5114 4.2085373 2000-2009 –> 2010-2019 8.638873 40525.5532 40611.8459
obs Indo-Pacific (26.75,27] cstar ~ temp + aou + nitrate + silicate + phosphate_star 25421.8465 25466.0845 5.3471266 2000-2009 –> 2010-2019 10.730726 52432.0244 52520.9095
mod Indo-Pacific (26.75,27] cstar ~ temp + nitrate + silicate + phosphate + phosphate_star 19575.7656 19618.7335 4.2142103 2000-2009 –> 2010-2019 8.649105 40542.1854 40628.4780
obs Indo-Pacific (26.75,27] cstar ~ temp + nitrate + silicate + phosphate + phosphate_star 25424.0298 25468.2678 5.3485491 2000-2009 –> 2010-2019 10.731469 52433.1082 52521.9933
mod Indo-Pacific (26.75,27] cstar ~ temp + nitrate + silicate + phosphate_star 19643.9328 19680.7624 4.2576254 2000-2009 –> 2010-2019 8.711919 40639.7940 40713.7591
obs Indo-Pacific (26.75,27] cstar ~ temp + nitrate + silicate + phosphate_star 25422.4594 25460.3777 5.3488291 2000-2009 –> 2010-2019 10.732429 52430.6374 52506.8247
mod Indo-Pacific (27,27.25] cstar ~ aou + nitrate + phosphate + phosphate_star 23659.7649 23698.1877 3.4204792 2000-2009 –> 2010-2019 6.651096 47929.2483 48006.3775
mod Indo-Pacific (27,27.25] cstar ~ aou + nitrate + silicate + phosphate + phosphate_star 23439.4573 23484.2839 3.3363607 2000-2009 –> 2010-2019 6.450639 47367.6575 47457.6414
obs Indo-Pacific (27,27.25] cstar ~ sal + aou + nitrate + silicate + phosphate 31977.7010 32023.8146 4.7586840 2000-2009 –> 2010-2019 9.521699 67207.2066 67300.1098
obs Indo-Pacific (27,27.25] cstar ~ sal + aou + nitrate + silicate + phosphate_star 32031.7691 32077.8826 4.7827233 2000-2009 –> 2010-2019 9.637994 67487.9909 67580.8941
obs Indo-Pacific (27,27.25] cstar ~ sal + aou + silicate + phosphate 31976.0039 32015.5298 4.7588183 2000-2009 –> 2010-2019 9.521845 67203.5380 67283.1693
obs Indo-Pacific (27,27.25] cstar ~ sal + aou + silicate + phosphate_star 32067.2383 32106.7642 4.7994538 2000-2009 –> 2010-2019 9.665534 67547.7432 67627.3745
obs Indo-Pacific (27,27.25] cstar ~ sal + nitrate + silicate + phosphate + phosphate_star 32000.9904 32047.1039 4.7690239 2000-2009 –> 2010-2019 9.590132 67373.7648 67466.6680
obs Indo-Pacific (27,27.25] cstar ~ sal + silicate + phosphate + phosphate_star 32029.2508 32068.7767 4.7824923 2000-2009 –> 2010-2019 9.616174 67430.8063 67510.4376
mod Indo-Pacific (27,27.25] cstar ~ sal + temp + aou + nitrate + phosphate 23677.9436 23722.7702 3.4266832 2000-2009 –> 2010-2019 6.665067 47971.9057 48061.8897
mod Indo-Pacific (27,27.25] cstar ~ sal + temp + aou + nitrate + phosphate_star 23852.1245 23896.9511 3.4941925 2000-2009 –> 2010-2019 6.833119 48432.2680 48522.2520
obs Indo-Pacific (27,27.25] cstar ~ sal + temp + aou + silicate + phosphate 31966.5716 32012.6852 4.7537508 2000-2009 –> 2010-2019 9.516773 67196.0960 67288.9992
obs Indo-Pacific (27,27.25] cstar ~ sal + temp + aou + silicate + phosphate_star 32031.5333 32077.6468 4.7826182 2000-2009 –> 2010-2019 9.598574 67391.6708 67484.5739
mod Indo-Pacific (27,27.25] cstar ~ sal + temp + nitrate + phosphate + phosphate_star 23862.3945 23907.2211 3.4982142 2000-2009 –> 2010-2019 6.844869 48464.1766 48554.1606
mod Indo-Pacific (27,27.25] cstar ~ sal + temp + nitrate + silicate + phosphate_star 23845.7177 23890.5443 3.4916859 2000-2009 –> 2010-2019 6.845709 48468.0860 48558.0700
obs Indo-Pacific (27,27.25] cstar ~ sal + temp + nitrate + silicate + phosphate_star 32164.5275 32210.6410 4.8422658 2000-2009 –> 2010-2019 9.693259 67610.3328 67703.2360
obs Indo-Pacific (27,27.25] cstar ~ sal + temp + silicate + phosphate + phosphate_star 32017.4076 32063.5211 4.7763262 2000-2009 –> 2010-2019 9.582086 67352.4982 67445.4014
mod Indo-Pacific (27,27.25] cstar ~ temp + aou + nitrate + phosphate 23759.9383 23798.3611 3.4590735 2000-2009 –> 2010-2019 6.752342 48209.2058 48286.3349
mod Indo-Pacific (27,27.25] cstar ~ temp + aou + nitrate + silicate + phosphate 23646.8278 23691.6544 3.4147613 2000-2009 –> 2010-2019 6.639594 47901.5389 47991.5229
mod Indo-Pacific (27,27.25] cstar ~ temp + aou + nitrate + silicate + phosphate_star 23845.4248 23890.2514 3.4915713 2000-2009 –> 2010-2019 6.830979 48426.9173 48516.9013
mod Indo-Pacific (27,27.25] cstar ~ temp + nitrate + silicate + phosphate + phosphate_star 23849.1506 23893.9772 3.4930287 2000-2009 –> 2010-2019 6.836731 48442.6723 48532.6563
mod Indo-Pacific (27.25,27.5] cstar ~ sal + aou + nitrate + silicate + phosphate 30679.4551 30726.1277 3.3859381 2000-2009 –> 2010-2019 6.495305 60889.4621 60982.9290
obs Indo-Pacific (27.25,27.5] cstar ~ sal + aou + nitrate + silicate + phosphate 32915.4449 32962.2735 3.8559260 2000-2009 –> 2010-2019 8.162799 67945.6745 68039.4936
mod Indo-Pacific (27.25,27.5] cstar ~ sal + aou + nitrate + silicate + phosphate_star 31020.4639 31067.1364 3.4867590 2000-2009 –> 2010-2019 6.809320 62014.7341 62108.2010
obs Indo-Pacific (27.25,27.5] cstar ~ sal + aou + nitrate + silicate + phosphate_star 33149.4884 33196.3170 3.9326173 2000-2009 –> 2010-2019 8.424969 68692.5190 68786.3381
obs Indo-Pacific (27.25,27.5] cstar ~ sal + aou + silicate + phosphate 33003.0685 33043.2073 3.8851155 2000-2009 –> 2010-2019 8.201440 68057.9586 68138.3750
obs Indo-Pacific (27.25,27.5] cstar ~ sal + aou + silicate + phosphate_star 33306.4095 33346.5483 3.9855602 2000-2009 –> 2010-2019 8.501361 68910.7585 68991.1748
mod Indo-Pacific (27.25,27.5] cstar ~ sal + nitrate + silicate + phosphate + phosphate_star 30884.0447 30930.7172 3.4460707 2000-2009 –> 2010-2019 6.705056 61649.8368 61743.3037
obs Indo-Pacific (27.25,27.5] cstar ~ sal + nitrate + silicate + phosphate + phosphate_star 33005.8272 33052.6558 3.8853636 2000-2009 –> 2010-2019 8.283231 68290.3373 68384.1564
obs Indo-Pacific (27.25,27.5] cstar ~ sal + silicate + phosphate + phosphate_star 33182.3686 33222.5074 3.9441767 2000-2009 –> 2010-2019 8.382933 68577.4296 68657.8460
mod Indo-Pacific (27.25,27.5] cstar ~ sal + temp + aou + nitrate + phosphate 31359.7517 31406.4243 3.5900504 2000-2009 –> 2010-2019 6.794226 61924.9616 62018.4285
mod Indo-Pacific (27.25,27.5] cstar ~ sal + temp + aou + nitrate + phosphate_star 31520.4160 31567.0885 3.6400244 2000-2009 –> 2010-2019 6.908722 62321.3981 62414.8650
obs Indo-Pacific (27.25,27.5] cstar ~ sal + temp + aou + silicate + phosphate 32953.4188 33000.2474 3.8682669 2000-2009 –> 2010-2019 8.137392 67876.5855 67970.4045
obs Indo-Pacific (27.25,27.5] cstar ~ sal + temp + aou + silicate + phosphate_star 32988.1261 33034.9547 3.8795807 2000-2009 –> 2010-2019 8.196091 68045.5408 68139.3599
mod Indo-Pacific (27.25,27.5] cstar ~ sal + temp + nitrate + silicate + phosphate_star 31055.3627 31102.0353 3.4972449 2000-2009 –> 2010-2019 6.568616 61119.9707 61213.4376
obs Indo-Pacific (27.25,27.5] cstar ~ sal + temp + nitrate + silicate + phosphate_star 33155.0261 33201.8547 3.9344503 2000-2009 –> 2010-2019 8.256703 68228.6102 68322.4293
obs Indo-Pacific (27.25,27.5] cstar ~ sal + temp + silicate + phosphate + phosphate_star 32980.9821 33027.8107 3.8772492 2000-2009 –> 2010-2019 8.184925 68013.4792 68107.2983
mod Indo-Pacific (27.25,27.5] cstar ~ temp + aou + nitrate + silicate + phosphate 31274.7267 31321.3992 3.5638819 2000-2009 –> 2010-2019 6.701625 61592.1704 61685.6373
mod Indo-Pacific (27.25,27.5] cstar ~ temp + aou + nitrate + silicate + phosphate_star 31483.7744 31530.4469 3.6285663 2000-2009 –> 2010-2019 6.848048 62105.3435 62198.8104
mod Indo-Pacific (27.25,27.5] cstar ~ temp + nitrate + silicate + phosphate + phosphate_star 31478.5017 31525.1743 3.6269205 2000-2009 –> 2010-2019 6.847867 62105.4498 62198.9167
mod Indo-Pacific (27.25,27.5] cstar ~ temp + nitrate + silicate + phosphate_star 31495.1328 31535.1378 3.6327394 2000-2009 –> 2010-2019 6.855873 62128.1055 62208.2200
obs Indo-Pacific (27.5,27.75] cstar ~ sal + aou + nitrate + phosphate + phosphate_star 37831.8361 37879.5260 4.0345249 2000-2009 –> 2010-2019 7.814059 75429.3639 75524.8645
mod Indo-Pacific (27.5,27.75] cstar ~ sal + aou + nitrate + silicate + phosphate 32964.6319 33011.9978 3.1547522 2000-2009 –> 2010-2019 6.090693 67837.3015 67932.6260
obs Indo-Pacific (27.5,27.75] cstar ~ sal + aou + nitrate + silicate + phosphate 36300.1231 36347.8130 3.5999567 2000-2009 –> 2010-2019 7.494234 74306.6024 74402.1031
mod Indo-Pacific (27.5,27.75] cstar ~ sal + aou + nitrate + silicate + phosphate_star 32837.7726 32885.1385 3.1237175 2000-2009 –> 2010-2019 6.068431 67752.1083 67847.4328
obs Indo-Pacific (27.5,27.75] cstar ~ sal + aou + silicate + phosphate 36576.2326 36617.1096 3.6752255 2000-2009 –> 2010-2019 7.600616 74689.5252 74771.3829
mod Indo-Pacific (27.5,27.75] cstar ~ sal + aou + silicate + phosphate_star 33406.5748 33447.1741 3.2658057 2000-2009 –> 2010-2019 6.339404 68917.1794 68998.8861
mod Indo-Pacific (27.5,27.75] cstar ~ sal + nitrate + silicate + phosphate + phosphate_star 32862.0599 32909.4257 3.1296354 2000-2009 –> 2010-2019 6.068768 67749.9015 67845.2259
obs Indo-Pacific (27.5,27.75] cstar ~ sal + nitrate + silicate + phosphate + phosphate_star 36579.0153 36626.7053 3.6754395 2000-2009 –> 2010-2019 7.774554 75286.4646 75381.9653
obs Indo-Pacific (27.5,27.75] cstar ~ sal + temp + aou + nitrate + phosphate 37696.2085 37743.8984 3.9940159 2000-2009 –> 2010-2019 7.805878 75410.1960 75505.6966
mod Indo-Pacific (27.5,27.75] cstar ~ sal + temp + aou + silicate + phosphate 32961.8011 33009.1670 3.1540564 2000-2009 –> 2010-2019 6.197617 68337.2449 68432.5694
obs Indo-Pacific (27.5,27.75] cstar ~ sal + temp + aou + silicate + phosphate 36262.0455 36309.7354 3.5897719 2000-2009 –> 2010-2019 7.278464 73526.8993 73622.4000
mod Indo-Pacific (27.5,27.75] cstar ~ sal + temp + aou + silicate + phosphate_star 32972.9339 33020.2998 3.1567940 2000-2009 –> 2010-2019 6.214293 68412.1933 68507.5178
obs Indo-Pacific (27.5,27.75] cstar ~ sal + temp + aou + silicate + phosphate_star 36197.9635 36245.6534 3.5726966 2000-2009 –> 2010-2019 7.290699 73571.0428 73666.5434
mod Indo-Pacific (27.5,27.75] cstar ~ sal + temp + nitrate + silicate + phosphate 32756.7029 32804.0687 3.1040447 2000-2009 –> 2010-2019 6.152639 68155.2276 68250.5521
mod Indo-Pacific (27.5,27.75] cstar ~ sal + temp + nitrate + silicate + phosphate_star 33392.7404 33440.1062 3.2617783 2000-2009 –> 2010-2019 6.242819 68478.3151 68573.6396
obs Indo-Pacific (27.5,27.75] cstar ~ sal + temp + nitrate + silicate + phosphate_star 36357.4725 36405.1624 3.6153508 2000-2009 –> 2010-2019 7.325921 73703.1907 73798.6913
mod Indo-Pacific (27.5,27.75] cstar ~ sal + temp + silicate + phosphate 33028.7850 33069.3843 3.1710580 2000-2009 –> 2010-2019 6.252686 68575.8267 68657.5334
mod Indo-Pacific (27.5,27.75] cstar ~ sal + temp + silicate + phosphate + phosphate_star 32948.7937 32996.1596 3.1508608 2000-2009 –> 2010-2019 6.201233 68355.4624 68450.7869
obs Indo-Pacific (27.5,27.75] cstar ~ sal + temp + silicate + phosphate + phosphate_star 36206.8464 36254.5363 3.5750587 2000-2009 –> 2010-2019 7.291025 73572.4332 73667.9339
obs Indo-Pacific (27.5,27.75] cstar ~ sal + temp + silicate + phosphate_star 36544.1448 36585.0219 3.6664614 2000-2009 –> 2010-2019 7.406779 73997.0483 74078.9060
mod Indo-Pacific (27.75,27.85] cstar ~ sal + aou + nitrate + silicate + phosphate 10943.9421 10984.6231 2.2133263 2000-2009 –> 2010-2019 4.604850 23456.9247 23538.9849
mod Indo-Pacific (27.75,27.85] cstar ~ sal + aou + nitrate + silicate + phosphate_star 10962.0426 11002.7236 2.2214542 2000-2009 –> 2010-2019 4.630210 23514.1970 23596.2572
mod Indo-Pacific (27.75,27.85] cstar ~ sal + temp + aou + nitrate + phosphate 10928.9522 10969.6331 2.2066176 2000-2009 –> 2010-2019 4.510686 23238.6744 23320.7346
mod Indo-Pacific (27.75,27.85] cstar ~ sal + temp + aou + nitrate + phosphate_star 10932.3431 10973.0241 2.2081334 2000-2009 –> 2010-2019 4.521258 23263.4688 23345.5290
mod Indo-Pacific (27.75,27.85] cstar ~ sal + temp + aou + nitrate + silicate 10994.5018 11035.1827 2.2361047 2000-2009 –> 2010-2019 4.627946 23508.2090 23590.2692
mod Indo-Pacific (27.75,27.85] cstar ~ sal + temp + aou + phosphate 11037.1259 11071.9953 2.2564037 2000-2009 –> 2010-2019 4.623015 23490.9758 23561.3131
mod Indo-Pacific (27.75,27.85] cstar ~ sal + temp + aou + silicate + phosphate 11035.8347 11076.5157 2.2549003 2000-2009 –> 2010-2019 4.619367 23486.7382 23568.7985
mod Indo-Pacific (27.75,27.85] cstar ~ sal + temp + aou + silicate + phosphate_star 11038.6086 11079.2895 2.2561673 2000-2009 –> 2010-2019 4.629316 23509.5091 23591.5694
obs Indo-Pacific (27.75,27.85] cstar ~ sal + temp + aou + silicate + phosphate_star 19205.4458 19248.7549 3.4937998 2000-2009 –> 2010-2019 6.763560 37132.5923 37218.9040
mod Indo-Pacific (27.75,27.85] cstar ~ sal + temp + nitrate + phosphate + phosphate_star 10935.0860 10975.7670 2.2093603 2000-2009 –> 2010-2019 4.529340 23282.3566 23364.4168
obs Indo-Pacific (27.75,27.85] cstar ~ sal + temp + nitrate + silicate + phosphate_star 19180.5893 19223.8984 3.4817390 2000-2009 –> 2010-2019 6.759168 37123.8546 37210.1663
mod Indo-Pacific (27.75,27.85] cstar ~ sal + temp + silicate + phosphate + phosphate_star 11039.2474 11079.9284 2.2564592 2000-2009 –> 2010-2019 4.632449 23516.6747 23598.7349
obs Indo-Pacific (27.75,27.85] cstar ~ sal + temp + silicate + phosphate + phosphate_star 19203.3851 19246.6943 3.4927984 2000-2009 –> 2010-2019 6.762709 37130.8493 37217.1610
obs Indo-Pacific (27.75,27.85] cstar ~ sal + temp + silicate + phosphate_star 19211.4611 19248.5832 3.4976980 2000-2009 –> 2010-2019 6.787498 37178.6461 37252.6276
obs Indo-Pacific (27.75,27.85] cstar ~ temp + aou + nitrate + silicate + phosphate 19249.4011 19292.7103 3.5152303 2000-2009 –> 2010-2019 6.784509 37175.5353 37261.8471
obs Indo-Pacific (27.75,27.85] cstar ~ temp + aou + nitrate + silicate + phosphate_star 19195.3487 19238.6578 3.4888955 2000-2009 –> 2010-2019 6.758920 37123.0527 37209.3645
obs Indo-Pacific (27.75,27.85] cstar ~ temp + aou + silicate + phosphate_star 19261.1825 19298.3047 3.5219764 2000-2009 –> 2010-2019 6.792593 37188.1331 37262.1146
obs Indo-Pacific (27.75,27.85] cstar ~ temp + nitrate + silicate + phosphate + phosphate_star 19201.0158 19244.3249 3.4916473 2000-2009 –> 2010-2019 6.761914 37129.2284 37215.5402
obs Indo-Pacific (27.75,27.85] cstar ~ temp + nitrate + silicate + phosphate_star 19218.7732 19255.8953 3.5012578 2000-2009 –> 2010-2019 6.780368 37163.5654 37237.5468
obs Indo-Pacific (27.75,27.85] cstar ~ temp + silicate + phosphate + phosphate_star 19257.2638 19294.3859 3.5200568 2000-2009 –> 2010-2019 6.790874 37184.6358 37258.6173
mod Indo-Pacific (27.85,27.95] cstar ~ sal + aou + nitrate + silicate + phosphate 12788.2304 12830.3860 1.9671223 2000-2009 –> 2010-2019 4.081717 26832.2438 26916.9741
mod Indo-Pacific (27.85,27.95] cstar ~ sal + aou + nitrate + silicate + phosphate_star 12771.7141 12813.8698 1.9617998 2000-2009 –> 2010-2019 4.038337 26698.1861 26782.9164
mod Indo-Pacific (27.85,27.95] cstar ~ sal + temp + aou + nitrate 12702.1437 12738.2771 1.9401746 2000-2009 –> 2010-2019 3.945207 26399.8256 26472.4517
mod Indo-Pacific (27.85,27.95] cstar ~ sal + temp + aou + nitrate + phosphate 12585.4427 12627.5984 1.9027610 2000-2009 –> 2010-2019 3.906457 26280.8114 26365.5417
mod Indo-Pacific (27.85,27.95] cstar ~ sal + temp + aou + nitrate + phosphate_star 12583.9697 12626.1254 1.9023013 2000-2009 –> 2010-2019 3.906707 26281.6309 26366.3612
mod Indo-Pacific (27.85,27.95] cstar ~ sal + temp + aou + nitrate + silicate 12588.9220 12631.0776 1.9038473 2000-2009 –> 2010-2019 3.885475 26212.6145 26297.3449
mod Indo-Pacific (27.85,27.95] cstar ~ sal + temp + aou + silicate + phosphate 12762.5574 12804.7131 1.9588553 2000-2009 –> 2010-2019 4.081125 26830.0191 26914.7495
obs Indo-Pacific (27.85,27.95] cstar ~ sal + temp + aou + silicate + phosphate 26684.7282 26730.5337 3.2495689 2000-2009 –> 2010-2019 6.158403 50327.2644 50418.3327
mod Indo-Pacific (27.85,27.95] cstar ~ sal + temp + aou + silicate + phosphate_star 12763.5710 12805.7267 1.9591810 2000-2009 –> 2010-2019 4.086962 26847.8163 26932.5467
obs Indo-Pacific (27.85,27.95] cstar ~ sal + temp + aou + silicate + phosphate_star 26645.5975 26691.4030 3.2372085 2000-2009 –> 2010-2019 6.154834 50316.8074 50407.8757
mod Indo-Pacific (27.85,27.95] cstar ~ sal + temp + nitrate + phosphate + phosphate_star 12582.1494 12624.3051 1.9017334 2000-2009 –> 2010-2019 3.916307 26312.5577 26397.2880
obs Indo-Pacific (27.85,27.95] cstar ~ sal + temp + nitrate + silicate + phosphate_star 26624.1402 26669.9456 3.2304507 2000-2009 –> 2010-2019 6.157118 50324.7524 50415.8207
mod Indo-Pacific (27.85,27.95] cstar ~ sal + temp + silicate + phosphate + phosphate_star 12763.8761 12806.0318 1.9592791 2000-2009 –> 2010-2019 4.089624 26855.9157 26940.6460
obs Indo-Pacific (27.85,27.95] cstar ~ sal + temp + silicate + phosphate + phosphate_star 26644.2410 26690.0465 3.2367809 2000-2009 –> 2010-2019 6.155551 50319.1782 50410.2465
obs Indo-Pacific (27.85,27.95] cstar ~ temp + aou + nitrate + silicate + phosphate 26657.8549 26703.6604 3.2410753 2000-2009 –> 2010-2019 6.149866 50300.2496 50391.3178
obs Indo-Pacific (27.85,27.95] cstar ~ temp + aou + nitrate + silicate + phosphate_star 26630.7209 26676.5264 3.2325218 2000-2009 –> 2010-2019 6.150144 50301.9204 50392.9887
obs Indo-Pacific (27.85,27.95] cstar ~ temp + aou + silicate + phosphate_star 26666.3250 26705.5868 3.2443818 2000-2009 –> 2010-2019 6.162013 50335.5532 50413.6117
obs Indo-Pacific (27.85,27.95] cstar ~ temp + nitrate + silicate + phosphate + phosphate_star 26631.6548 26677.4603 3.2328158 2000-2009 –> 2010-2019 6.151528 50306.4032 50397.4715
obs Indo-Pacific (27.85,27.95] cstar ~ temp + nitrate + silicate + phosphate_star 26634.6896 26673.9515 3.2344014 2000-2009 –> 2010-2019 6.161162 50333.6049 50411.6634
obs Indo-Pacific (27.85,27.95] cstar ~ temp + silicate + phosphate + phosphate_star 26664.6541 26703.9159 3.2438539 2000-2009 –> 2010-2019 6.162624 50337.5920 50415.6505
obs Indo-Pacific (27.95,28.05] cstar ~ aou + nitrate + phosphate + phosphate_star 44167.5728 44209.9229 3.1620941 2000-2009 –> 2010-2019 6.013512 84644.2184 84728.6414
obs Indo-Pacific (27.95,28.05] cstar ~ aou + nitrate + silicate + phosphate + phosphate_star 44129.6690 44179.0775 3.1547581 2000-2009 –> 2010-2019 5.996881 84554.7586 84653.2520
obs Indo-Pacific (27.95,28.05] cstar ~ aou + phosphate + phosphate_star 44203.9972 44239.2890 3.1691743 2000-2009 –> 2010-2019 6.021223 84682.2708 84752.6232
obs Indo-Pacific (27.95,28.05] cstar ~ aou + silicate + phosphate + phosphate_star 44163.3527 44205.7029 3.1613175 2000-2009 –> 2010-2019 6.006583 84604.5691 84688.9921
obs Indo-Pacific (27.95,28.05] cstar ~ sal + aou + nitrate + phosphate + phosphate_star 44169.2004 44218.6089 3.1620256 2000-2009 –> 2010-2019 6.013323 84647.1563 84745.6498
obs Indo-Pacific (27.95,28.05] cstar ~ sal + aou + phosphate + phosphate_star 44205.8318 44248.1820 3.1691438 2000-2009 –> 2010-2019 6.020882 84684.3224 84768.7453
mod Indo-Pacific (27.95,28.05] cstar ~ sal + temp + aou + nitrate 16644.5967 16682.5833 1.7941045 2000-2009 –> 2010-2019 3.529143 30607.9103 30682.9299
mod Indo-Pacific (27.95,28.05] cstar ~ sal + temp + aou + nitrate + silicate 16358.8661 16403.1839 1.7329897 2000-2009 –> 2010-2019 3.425438 30148.1681 30235.6910
mod Indo-Pacific (27.95,28.05] cstar ~ sal + temp + aou + phosphate 16787.1456 16825.1322 1.8251761 2000-2009 –> 2010-2019 3.685693 31244.9572 31319.9768
mod Indo-Pacific (27.95,28.05] cstar ~ sal + temp + aou + phosphate_star 16784.7001 16822.6868 1.8246385 2000-2009 –> 2010-2019 3.689634 31259.5409 31334.5605
mod Indo-Pacific (27.95,28.05] cstar ~ sal + temp + aou + silicate + phosphate 16524.5224 16568.8401 1.7679166 2000-2009 –> 2010-2019 3.604437 30892.4013 30979.9242
obs Indo-Pacific (27.95,28.05] cstar ~ sal + temp + aou + silicate + phosphate 44094.6839 44144.0923 3.1483403 2000-2009 –> 2010-2019 5.998646 84566.9311 84665.4245
mod Indo-Pacific (27.95,28.05] cstar ~ sal + temp + aou + silicate + phosphate_star 16525.1284 16569.4461 1.7680457 2000-2009 –> 2010-2019 3.610531 30915.9697 31003.4926
mod Indo-Pacific (27.95,28.05] cstar ~ sal + temp + nitrate + silicate + phosphate_star 16552.2811 16596.5988 1.7738378 2000-2009 –> 2010-2019 3.660085 31109.3662 31196.8891
obs Indo-Pacific (27.95,28.05] cstar ~ sal + temp + nitrate + silicate + phosphate_star 44167.9293 44217.3377 3.1617916 2000-2009 –> 2010-2019 6.032945 84759.7199 84858.2133
mod Indo-Pacific (27.95,28.05] cstar ~ sal + temp + phosphate + phosphate_star 16783.1594 16821.1461 1.8242999 2000-2009 –> 2010-2019 3.691553 31266.5667 31341.5864
mod Indo-Pacific (27.95,28.05] cstar ~ sal + temp + silicate + phosphate + phosphate_star 16525.2304 16569.5482 1.7680674 2000-2009 –> 2010-2019 3.612282 30922.7167 31010.2396
mod Indo-Pacific (27.95,28.05] cstar ~ sal + temp + silicate + phosphate_star 16559.5556 16597.5422 1.7758205 2000-2009 –> 2010-2019 3.678135 31174.7082 31249.7278
obs Indo-Pacific (27.95,28.05] cstar ~ temp + aou + nitrate + silicate + phosphate 44167.4467 44216.8552 3.1617028 2000-2009 –> 2010-2019 6.016810 84667.3063 84765.7997
obs Indo-Pacific (27.95,28.05] cstar ~ temp + aou + silicate + phosphate 44208.7743 44251.1244 3.1696866 2000-2009 –> 2010-2019 6.036583 84774.2258 84858.6487
obs Indo-Pacific (28.05,28.1] cstar ~ aou + nitrate + phosphate + phosphate_star 39686.7016 39728.5249 3.0111777 2000-2009 –> 2010-2019 5.681985 69066.8346 69148.9690
obs Indo-Pacific (28.05,28.1] cstar ~ aou + silicate + phosphate + phosphate_star 39781.6939 39823.5173 3.0294100 2000-2009 –> 2010-2019 5.704411 69181.0194 69263.1538
obs Indo-Pacific (28.05,28.1] cstar ~ sal + aou + nitrate + silicate + phosphate 39483.2835 39532.0774 2.9721252 2000-2009 –> 2010-2019 5.631463 68812.7844 68908.6078
obs Indo-Pacific (28.05,28.1] cstar ~ sal + aou + silicate + phosphate 39524.0981 39565.9214 2.9802228 2000-2009 –> 2010-2019 5.646661 68884.2109 68966.3453
obs Indo-Pacific (28.05,28.1] cstar ~ sal + nitrate + silicate + phosphate + phosphate_star 39571.9018 39620.6958 2.9889101 2000-2009 –> 2010-2019 5.698125 69128.6579 69224.4813
mod Indo-Pacific (28.05,28.1] cstar ~ sal + temp + aou + nitrate + silicate 8796.0101 8836.2704 1.5994595 2000-2009 –> 2010-2019 3.225834 16808.4264 16888.2312
mod Indo-Pacific (28.05,28.1] cstar ~ sal + temp + aou + silicate + phosphate 8722.8374 8763.0977 1.5744873 2000-2009 –> 2010-2019 3.137767 16569.1486 16648.9534
obs Indo-Pacific (28.05,28.1] cstar ~ sal + temp + aou + silicate + phosphate 39522.8704 39571.6643 2.9796115 2000-2009 –> 2010-2019 5.613791 68736.1223 68831.9457
mod Indo-Pacific (28.05,28.1] cstar ~ sal + temp + aou + silicate + phosphate_star 8722.7424 8763.0028 1.5744552 2000-2009 –> 2010-2019 3.142932 16582.9873 16662.7921
obs Indo-Pacific (28.05,28.1] cstar ~ sal + temp + aou + silicate + phosphate_star 39608.5537 39657.3476 2.9958799 2000-2009 –> 2010-2019 5.653845 68931.7407 69027.5641
mod Indo-Pacific (28.05,28.1] cstar ~ sal + temp + nitrate + silicate + phosphate_star 8729.8435 8770.1038 1.5768614 2000-2009 –> 2010-2019 3.145932 16591.6778 16671.4826
obs Indo-Pacific (28.05,28.1] cstar ~ sal + temp + nitrate + silicate + phosphate_star 39769.1988 39817.9928 3.0266208 2000-2009 –> 2010-2019 5.692335 69127.9884 69223.8118
mod Indo-Pacific (28.05,28.1] cstar ~ sal + temp + silicate + phosphate + phosphate_star 8722.7561 8763.0164 1.5744598 2000-2009 –> 2010-2019 3.143084 16583.3957 16663.2006
obs Indo-Pacific (28.05,28.1] cstar ~ sal + temp + silicate + phosphate + phosphate_star 39588.1757 39636.9696 2.9920028 2000-2009 –> 2010-2019 5.647606 68900.4906 68996.3140
mod Indo-Pacific (28.05,28.1] cstar ~ sal + temp + silicate + phosphate_star 8759.7059 8794.2147 1.5877033 2000-2009 –> 2010-2019 3.156901 16619.8804 16688.2845
obs Indo-Pacific (28.05,28.1] cstar ~ temp + aou + nitrate + silicate + phosphate 39850.3042 39899.0981 3.0422607 2000-2009 –> 2010-2019 5.701545 69179.5584 69275.3819
mod Indo-Pacific (28.05,28.1] cstar ~ temp + aou + silicate + phosphate 9029.8609 9064.3697 1.6826777 2000-2009 –> 2010-2019 3.275636 16953.1246 17021.5287
mod Indo-Pacific (28.05,28.1] cstar ~ temp + aou + silicate + phosphate_star 9035.3407 9069.8496 1.6846618 2000-2009 –> 2010-2019 3.284429 16976.5081 17044.9123
mod Indo-Pacific (28.05,28.1] cstar ~ temp + nitrate + silicate + phosphate_star 9037.2102 9071.7191 1.6853393 2000-2009 –> 2010-2019 3.284745 16977.4295 17045.8336
mod Indo-Pacific (28.05,28.1] cstar ~ temp + silicate + phosphate + phosphate_star 9034.0132 9068.5220 1.6841809 2000-2009 –> 2010-2019 3.284111 16975.6097 17044.0138
obs Indo-Pacific (28.1,28.2] cstar ~ sal + aou + nitrate + phosphate + phosphate_star 26769.9138 26816.2255 2.7320261 2000-2009 –> 2010-2019 5.363394 53683.8913 53776.6615
obs Indo-Pacific (28.1,28.2] cstar ~ sal + aou + nitrate + silicate + phosphate 26630.0387 26676.3504 2.6976240 2000-2009 –> 2010-2019 5.347440 53622.7664 53715.5365
obs Indo-Pacific (28.1,28.2] cstar ~ sal + aou + phosphate + phosphate_star 26826.7064 26866.4021 2.7466168 2000-2009 –> 2010-2019 5.381952 53755.6667 53835.1840
mod Indo-Pacific (28.1,28.2] cstar ~ sal + aou + silicate + phosphate 18172.1863 18211.4621 1.4127764 2000-2009 –> 2010-2019 2.975100 37428.4208 37506.9877
obs Indo-Pacific (28.1,28.2] cstar ~ sal + aou + silicate + phosphate 26643.4761 26683.1719 2.7013994 2000-2009 –> 2010-2019 5.354942 53650.0448 53729.5621
mod Indo-Pacific (28.1,28.2] cstar ~ sal + nitrate + silicate + phosphate 18135.7374 18175.0133 1.4077819 2000-2009 –> 2010-2019 2.899052 36911.7115 36990.2783
obs Indo-Pacific (28.1,28.2] cstar ~ sal + nitrate + silicate + phosphate + phosphate_star 26574.4090 26620.7206 2.6840626 2000-2009 –> 2010-2019 5.380728 53764.6871 53857.4573
obs Indo-Pacific (28.1,28.2] cstar ~ sal + silicate + phosphate + phosphate_star 26592.2956 26631.9913 2.6889027 2000-2009 –> 2010-2019 5.392587 53809.8772 53889.3945
mod Indo-Pacific (28.1,28.2] cstar ~ sal + temp + aou + nitrate + silicate 16999.0546 17044.8764 1.2603352 2000-2009 –> 2010-2019 2.779837 35970.5319 36062.1932
obs Indo-Pacific (28.1,28.2] cstar ~ sal + temp + aou + phosphate + phosphate_star 26776.5208 26822.8324 2.7336619 2000-2009 –> 2010-2019 5.367064 53699.2094 53791.9796
mod Indo-Pacific (28.1,28.2] cstar ~ sal + temp + aou + silicate + phosphate 16935.7385 16981.5604 1.2526054 2000-2009 –> 2010-2019 2.634032 34924.2676 35015.9289
obs Indo-Pacific (28.1,28.2] cstar ~ sal + temp + aou + silicate + phosphate 26622.9517 26669.2633 2.6958925 2000-2009 –> 2010-2019 5.324109 53523.4232 53616.1934
mod Indo-Pacific (28.1,28.2] cstar ~ sal + temp + aou + silicate + phosphate_star 16946.5887 16992.4106 1.2539267 2000-2009 –> 2010-2019 2.637522 34951.3036 35042.9649
obs Indo-Pacific (28.1,28.2] cstar ~ sal + temp + aou + silicate + phosphate_star 26596.2676 26642.5792 2.6893831 2000-2009 –> 2010-2019 5.347480 53624.1657 53716.9359
mod Indo-Pacific (28.1,28.2] cstar ~ sal + temp + nitrate + silicate 17454.8225 17494.0984 1.3176577 2000-2009 –> 2010-2019 2.933229 37056.8565 37135.4233
mod Indo-Pacific (28.1,28.2] cstar ~ sal + temp + nitrate + silicate + phosphate_star 16925.0404 16970.8622 1.2513041 2000-2009 –> 2010-2019 2.644002 34997.4130 35089.0744
mod Indo-Pacific (28.1,28.2] cstar ~ sal + temp + silicate + phosphate 17442.5273 17481.8031 1.3160846 2000-2009 –> 2010-2019 2.882248 36724.0874 36802.6543
mod Indo-Pacific (28.1,28.2] cstar ~ sal + temp + silicate + phosphate + phosphate_star 16944.7814 16990.6033 1.2537065 2000-2009 –> 2010-2019 2.637001 34947.2507 35038.9120
obs Indo-Pacific (28.1,28.2] cstar ~ sal + temp + silicate + phosphate + phosphate_star 26594.2725 26640.5842 2.6888970 2000-2009 –> 2010-2019 5.341227 53597.6918 53690.4619
mod Indo-Pacific (28.1,28.2] cstar ~ sal + temp + silicate + phosphate_star 17024.6251 17063.9010 1.2637159 2000-2009 –> 2010-2019 2.660601 35125.9708 35204.5376
obs Indo-Pacific (28.2,28.3] cstar ~ sal + aou + nitrate + phosphate + phosphate_star 7960.6681 7999.5631 1.9310787 2000-2009 –> 2010-2019 4.580626 15116.3567 15192.4068
mod Indo-Pacific (28.2,28.3] cstar ~ sal + aou + nitrate + silicate + phosphate 30151.6098 30204.4893 0.7043133 2000-2009 –> 2010-2019 1.377867 54477.6649 54582.2191
mod Indo-Pacific (28.2,28.3] cstar ~ sal + aou + nitrate + silicate + phosphate_star 33441.0402 33493.9197 0.7914262 2000-2009 –> 2010-2019 1.575212 61366.5133 61471.0675
obs Indo-Pacific (28.2,28.3] cstar ~ sal + aou + phosphate + phosphate_star 7958.8219 7992.1605 1.9311563 2000-2009 –> 2010-2019 4.580960 15112.7995 15177.9854
obs Indo-Pacific (28.2,28.3] cstar ~ sal + aou + silicate + phosphate + phosphate_star 7960.4781 7999.3731 1.9309828 2000-2009 –> 2010-2019 4.564923 15098.5372 15174.5873
mod Indo-Pacific (28.2,28.3] cstar ~ sal + nitrate + silicate + phosphate 30912.7296 30958.0549 0.7236273 2000-2009 –> 2010-2019 1.480916 58019.4771 58109.0950
mod Indo-Pacific (28.2,28.3] cstar ~ sal + nitrate + silicate + phosphate + phosphate_star 30397.1016 30449.9811 0.7104696 2000-2009 –> 2010-2019 1.389424 54912.7955 55017.3498
obs Indo-Pacific (28.2,28.3] cstar ~ sal + temp + aou + nitrate + phosphate 7977.2830 8016.1780 1.9394829 2000-2009 –> 2010-2019 4.593965 15138.5239 15214.5740
obs Indo-Pacific (28.2,28.3] cstar ~ sal + temp + aou + phosphate 7978.2731 8011.6116 1.9409992 2000-2009 –> 2010-2019 4.595536 15137.5764 15202.7622
obs Indo-Pacific (28.2,28.3] cstar ~ sal + temp + aou + phosphate + phosphate_star 7960.7468 7999.6418 1.9311184 2000-2009 –> 2010-2019 4.572405 15107.1173 15183.1674
mod Indo-Pacific (28.2,28.3] cstar ~ sal + temp + aou + silicate + phosphate 31802.0674 31854.9469 0.7467523 2000-2009 –> 2010-2019 1.590105 61467.0637 61571.6179
obs Indo-Pacific (28.2,28.3] cstar ~ sal + temp + aou + silicate + phosphate 7973.9145 8012.8095 1.9377760 2000-2009 –> 2010-2019 4.592159 15135.0443 15211.0944
mod Indo-Pacific (28.2,28.3] cstar ~ sal + temp + aou + silicate + phosphate_star 31834.5636 31887.4431 0.7476130 2000-2009 –> 2010-2019 1.591300 61508.9750 61613.5293
obs Indo-Pacific (28.2,28.3] cstar ~ sal + temp + nitrate + phosphate + phosphate_star 8021.8849 8060.7799 1.9622249 2000-2009 –> 2010-2019 4.635986 15204.7198 15280.7699
mod Indo-Pacific (28.2,28.3] cstar ~ sal + temp + nitrate + silicate + phosphate 30913.3059 30966.1854 0.7235908 2000-2009 –> 2010-2019 1.474075 57807.7167 57912.2709
mod Indo-Pacific (28.2,28.3] cstar ~ sal + temp + nitrate + silicate + phosphate_star 31956.9458 32009.8253 0.7508636 2000-2009 –> 2010-2019 1.611191 62095.1868 62199.7411
obs Indo-Pacific (28.2,28.3] cstar ~ sal + temp + phosphate + phosphate_star 8033.7835 8067.1221 1.9693659 2000-2009 –> 2010-2019 4.643241 15214.7463 15279.9321
mod Indo-Pacific (28.2,28.3] cstar ~ sal + temp + silicate + phosphate + phosphate_star 31831.7916 31884.6711 0.7475396 2000-2009 –> 2010-2019 1.591185 61505.0303 61609.5845
obs Indo-Pacific (28.2,28.3] cstar ~ sal + temp + silicate + phosphate + phosphate_star 8006.2217 8045.1166 1.9542081 2000-2009 –> 2010-2019 4.624477 15185.1573 15261.2074
mod Indo-Pacific (28.2,28.3] cstar ~ sal + temp + silicate + phosphate_star 32040.1148 32085.4401 0.7531342 2000-2009 –> 2010-2019 1.618322 62310.1545 62399.7724
mod Indo-Pacific (28.3,28.4] cstar ~ aou + nitrate + phosphate + phosphate_star 16327.2576 16366.5323 1.1812892 2000-2009 –> 2010-2019 2.110745 29393.3449 29471.5399
mod Indo-Pacific (28.3,28.4] cstar ~ aou + nitrate + silicate + phosphate + phosphate_star 16324.0338 16369.8543 1.1806897 2000-2009 –> 2010-2019 2.107692 29366.4828 29457.7104
mod Indo-Pacific (28.3,28.4] cstar ~ sal + aou + nitrate + phosphate + phosphate_star 16239.5855 16285.4060 1.1710396 2000-2009 –> 2010-2019 2.100117 29303.7253 29394.9529
obs Indo-Pacific (28.3,28.4] cstar ~ sal + aou + nitrate + phosphate + phosphate_star 2457.1767 2487.4473 2.1603882 2000-2009 –> 2010-2019 4.667075 5075.5977 5136.1262
obs Indo-Pacific (28.3,28.4] cstar ~ sal + aou + nitrate + silicate + phosphate 2460.4238 2490.6943 2.1666831 2000-2009 –> 2010-2019 4.585688 5039.1804 5099.7088
obs Indo-Pacific (28.3,28.4] cstar ~ sal + aou + nitrate + silicate + phosphate_star 2451.3737 2481.6442 2.1491836 2000-2009 –> 2010-2019 4.564580 5028.4671 5088.9956
obs Indo-Pacific (28.3,28.4] cstar ~ sal + aou + silicate + phosphate + phosphate_star 2505.4119 2535.6824 2.2558106 2000-2009 –> 2010-2019 4.684556 5088.6449 5149.1733
obs Indo-Pacific (28.3,28.4] cstar ~ sal + nitrate + silicate + phosphate + phosphate_star 2452.9253 2483.1958 2.1521737 2000-2009 –> 2010-2019 4.566793 5029.6599 5090.1884
obs Indo-Pacific (28.3,28.4] cstar ~ sal + silicate + phosphate + phosphate_star 2503.7143 2529.6604 2.2564219 2000-2009 –> 2010-2019 4.685706 5085.1944 5137.0759
mod Indo-Pacific (28.3,28.4] cstar ~ sal + temp + aou + nitrate + phosphate 16220.7757 16266.5962 1.1689010 2000-2009 –> 2010-2019 2.095285 29256.7529 29347.9805
mod Indo-Pacific (28.3,28.4] cstar ~ sal + temp + aou + nitrate + phosphate_star 16278.4799 16324.3004 1.1754743 2000-2009 –> 2010-2019 2.106746 29365.5021 29456.7297
obs Indo-Pacific (28.3,28.4] cstar ~ sal + temp + aou + nitrate + phosphate_star 2460.0386 2490.3091 2.1659353 2000-2009 –> 2010-2019 4.668329 5076.5500 5137.0785
obs Indo-Pacific (28.3,28.4] cstar ~ sal + temp + aou + silicate + phosphate 2489.9130 2520.1835 2.2246988 2000-2009 –> 2010-2019 4.656134 5074.3793 5134.9077
obs Indo-Pacific (28.3,28.4] cstar ~ sal + temp + aou + silicate + phosphate_star 2496.5458 2526.8163 2.2379604 2000-2009 –> 2010-2019 4.670537 5081.5349 5142.0634
mod Indo-Pacific (28.3,28.4] cstar ~ sal + temp + nitrate + phosphate + phosphate_star 16222.3127 16268.1332 1.1690756 2000-2009 –> 2010-2019 2.096552 29269.7263 29360.9539
obs Indo-Pacific (28.3,28.4] cstar ~ sal + temp + silicate + phosphate + phosphate_star 2492.0718 2522.3424 2.2290065 2000-2009 –> 2010-2019 4.658136 5075.4813 5136.0097
mod Indo-Pacific (28.3,28.4] cstar ~ temp + aou + nitrate + phosphate 16310.2663 16349.5410 1.1793403 2000-2009 –> 2010-2019 2.106100 29348.1759 29426.3710
mod Indo-Pacific (28.3,28.4] cstar ~ temp + aou + nitrate + silicate + phosphate 16305.7517 16351.5721 1.1785939 2000-2009 –> 2010-2019 2.102565 29316.4386 29407.6662
mod Indo-Pacific (28.3,28.4] cstar ~ temp + nitrate + phosphate + phosphate_star 16313.7560 16353.0307 1.1797403 2000-2009 –> 2010-2019 2.107640 29363.5973 29441.7924
mod Indo-Pacific (28.3,28.4] cstar ~ temp + nitrate + silicate + phosphate + phosphate_star 16309.3779 16355.1984 1.1790093 2000-2009 –> 2010-2019 2.104182 29332.6655 29423.8931

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

8.3 Predictor counts

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

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

# print table
lm_all_stats %>%
  gt(rowname_col = "gamma_slab",
     groupname_col = c("basin", "data_source")) %>% 
  summary_rows(
    groups = TRUE,
    fns = list(total = "sum")
  )
aou nitrate phosphate phosphate_star sal silicate temp
Atlantic - mod
(-Inf,26] 6 6 10 7 4 9 5
(26,26.5] 9 7 7 7 7 6 7
(26.5,26.75] 8 10 7 6 4 5 8
(26.75,27] 8 5 8 6 4 10 7
(27,27.25] 9 6 6 4 6 10 7
(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 95.00 81.00 79.00 84.00 105.00 98.00 104.00
Indo-Pacific - mod
(-Inf,26] 9 6 6 6 10 4 7
(26,26.5] 8 6 8 7 6 8 7
(26.5,26.75] 7 5 6 7 6 9 8
(26.75,27] 5 10 5 7 5 5 10
(27,27.25] 7 10 7 7 4 5 8
(27.25,27.5] 6 10 5 7 6 8 7
(27.5,27.75] 5 5 6 6 10 10 6
(27.75,27.85] 8 6 6 5 10 6 8
(27.85,27.95] 8 7 5 5 10 6 8
(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 4 5 4 10 10 8
(28.2,28.3] 4 6 6 6 10 10 6
(28.3,28.4] 7 10 9 7 4 3 7
total 89.00 91.00 82.00 87.00 107.00 100.00 110.00
Atlantic - obs
(-Inf,26] 8 8 6 6 8 5 8
(26,26.5] 9 7 8 7 5 8 5
(26.5,26.75] 6 6 7 7 6 7 10
(26.75,27] 5 5 5 8 4 10 10
(27,27.25] 6 8 5 8 6 6 10
(27.25,27.5] 6 5 6 8 8 5 10
(27.5,27.75] 7 3 9 7 4 10 7
(27.75,27.85] 6 5 8 6 7 10 6
(27.85,27.95] 5 4 7 6 10 10 5
(27.95,28.05] 5 5 6 7 10 10 5
(28.05,28.1] 5 6 3 10 10 6 7
(28.1,28.15] 3 6 5 10 7 6 10
(28.15,28.2] 7 4 7 6 4 9 10
(28.2, Inf] 5 3 9 9 10 6 5
total 83.00 75.00 91.00 105.00 99.00 108.00 108.00
Indo-Pacific - obs
(-Inf,26] 8 6 6 6 10 5 8
(26,26.5] 4 5 5 8 10 5 10
(26.5,26.75] 5 5 6 7 10 4 10
(26.75,27] 7 10 7 5 3 7 7
(27,27.25] 6 4 6 7 10 10 4
(27.25,27.5] 6 4 6 7 10 10 4
(27.5,27.75] 6 5 7 6 10 8 6
(27.75,27.85] 4 5 4 9 4 10 10
(27.85,27.95] 5 5 5 8 4 10 10
(27.95,28.05] 9 5 9 7 4 6 4
(28.05,28.1] 7 5 8 6 7 9 5
(28.1,28.2] 7 3 9 7 10 7 4
(28.2,28.3] 7 3 10 7 10 3 7
(28.3,28.4] 7 5 7 8 10 8 4
total 88.00 70.00 95.00 98.00 112.00 102.00 93.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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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