Last updated: 2021-05-31

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: 1c2aff8

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

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


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

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

1 Version ID

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

2 Required data

Required are:

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

3 Predictor combinations

Find all possible combinations of following considered predictor variables:

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

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

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

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

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

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

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

# convert to tibble
lm_all <- as_tibble(lm_all)

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

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

4 Apply predictor threshold

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

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

This results in a total number of MLR models of:

  • 112

5 Fit all models

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

  • cstar_tref

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

library(furrr)
library(purrr)

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



library(tictoc)
plan(sequential)

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

plan(multisession, workers = 3)

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


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

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

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

# expand with model definitions

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

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


# test parallel processing

# expand with model definitions


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


### test furrr

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

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

6 Tidy models

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

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

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

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

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

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

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

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


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

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

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

7 Prepare coeffcients

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

7.1 VIF threshold

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

  • 1000

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

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

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

7.2 Predictor selection

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

  • 10

The criterion used to select the best models was:

  • rmse

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

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

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

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

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

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


# print table
lm_best_eras %>% 
  kable() %>%
  add_header_above() %>%
  kable_styling() %>%
  scroll_box(width = "100%", height = "400px")
data_source basin gamma_slab model AIC BIC rmse eras rmse_sum aic_sum bic_sum
mod Atlantic (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 134.20222 147.72500 0.7862462 2000-2009 –> 2010-2019 4.708972 844.6335 877.9545
obs Atlantic (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 949.55407 970.24486 6.5225536 2000-2009 –> 2010-2019 14.393513 1416.2295 1452.1410
mod Atlantic (-Inf,26] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 82.44135 95.96413 0.4733383 2000-2009 –> 2010-2019 4.721359 812.7893 846.1103
mod Atlantic (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 112.33837 125.86115 0.6345515 2000-2009 –> 2010-2019 4.467732 816.9967 850.3177
obs Atlantic (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 938.02658 958.71737 6.2631058 2000-2009 –> 2010-2019 14.867857 1416.2895 1452.2010
mod Atlantic (-Inf,26] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 79.93748 93.46026 0.4618604 2000-2009 –> 2010-2019 4.192559 777.8210 811.1419
mod Atlantic (-Inf,26] cstar_tref ~ sal + nitrate + silicate + phosphate_star 87.37399 98.96494 0.5066282 2000-2009 –> 2010-2019 4.759928 816.0324 844.5932
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate 953.24023 970.97519 6.6544622 2000-2009 –> 2010-2019 14.562664 1418.5293 1449.3106
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 950.16567 970.85646 6.5366151 2000-2009 –> 2010-2019 14.437280 1417.3308 1453.2423
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 950.22295 970.91374 6.5379337 2000-2009 –> 2010-2019 14.437606 1417.3717 1453.2832
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 81.32743 94.85021 0.4681972 2000-2009 –> 2010-2019 4.790806 816.0268 849.3478
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 940.63643 961.32722 6.3209265 2000-2009 –> 2010-2019 14.145662 1406.5461 1442.4576
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 111.71440 125.23718 0.6306816 2000-2009 –> 2010-2019 4.375467 810.5401 843.8611
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 939.26983 959.96062 6.2905835 2000-2009 –> 2010-2019 14.964265 1418.5700 1454.4815
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 134.12183 147.64461 0.7856267 2000-2009 –> 2010-2019 4.686882 843.1811 876.5020
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 950.47513 971.16592 6.5437416 2000-2009 –> 2010-2019 14.464332 1417.9677 1453.8792
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 84.51038 98.03316 0.4830378 2000-2009 –> 2010-2019 4.695626 812.7643 846.0853
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 111.62890 125.15168 0.6301531 2000-2009 –> 2010-2019 4.632465 827.0816 860.4026
obs Atlantic (-Inf,26] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 956.68052 977.37131 6.6882956 2000-2009 –> 2010-2019 14.945348 1429.5813 1465.4928
obs Atlantic (-Inf,26] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 953.92448 974.61527 6.6237038 2000-2009 –> 2010-2019 14.851347 1426.3614 1462.2729
obs Atlantic (26,26.5] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 4296.23330 4323.28591 5.8913530 2000-2009 –> 2010-2019 11.616291 8187.0156 8240.5784
obs Atlantic (26,26.5] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 4298.23270 4329.79409 5.8913504 2000-2009 –> 2010-2019 11.582341 8183.7235 8246.2135
mod Atlantic (26,26.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 2199.44233 2227.46954 3.5935423 2000-2009 –> 2010-2019 7.492506 3786.7239 3840.2692
obs Atlantic (26,26.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 4193.45466 4225.01605 5.4488752 2000-2009 –> 2010-2019 11.169105 8085.2283 8147.7183
mod Atlantic (26,26.5] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2167.19620 2195.22341 3.4532934 2000-2009 –> 2010-2019 7.755772 3810.2177 3863.7630
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 2203.47295 2231.50016 3.6114686 2000-2009 –> 2010-2019 7.871640 3840.9013 3894.4467
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 4210.96892 4242.53030 5.5204538 2000-2009 –> 2010-2019 11.377774 8131.7783 8194.2682
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 2203.82485 2231.85206 3.6130379 2000-2009 –> 2010-2019 7.928880 3848.6016 3902.1470
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 4212.83420 4244.39559 5.5281322 2000-2009 –> 2010-2019 11.472910 8151.8140 8214.3040
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + silicate 2225.99637 2254.02358 3.7133008 2000-2009 –> 2010-2019 8.344760 3910.7211 3964.2664
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 2171.23707 2199.26428 3.4705640 2000-2009 –> 2010-2019 8.215823 3869.7009 3923.2462
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2172.01245 2200.03966 3.4738878 2000-2009 –> 2010-2019 8.321335 3882.5356 3936.0809
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 2204.99686 2233.02407 3.6182695 2000-2009 –> 2010-2019 8.002621 3858.6876 3912.2330
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 4223.25793 4254.81931 5.5712380 2000-2009 –> 2010-2019 11.608542 8181.1725 8243.6625
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 2165.05710 2193.08431 3.4441857 2000-2009 –> 2010-2019 8.401023 3888.2109 3941.7562
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 2172.42427 2200.45148 3.4756545 2000-2009 –> 2010-2019 8.345684 3885.5781 3939.1234
obs Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + phosphate 4296.28730 4323.33991 5.8915901 2000-2009 –> 2010-2019 11.750187 8215.3637 8268.9265
obs Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 4298.28303 4329.84441 5.8915713 2000-2009 –> 2010-2019 11.651776 8198.5945 8261.0845
obs Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 4299.25056 4330.81194 5.8958204 2000-2009 –> 2010-2019 11.713850 8211.8083 8274.2982
obs Atlantic (26,26.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 4300.92771 4332.48909 5.9031933 2000-2009 –> 2010-2019 11.783572 8226.5540 8289.0439
mod Atlantic (26.5,26.75] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 2948.72915 2978.57294 3.9591985 2000-2009 –> 2010-2019 7.528083 5336.3601 5394.8272
mod Atlantic (26.5,26.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 3033.90129 3063.74508 4.2937391 2000-2009 –> 2010-2019 7.582105 5349.3298 5407.7969
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 5252.77984 5286.01295 5.2355023 2000-2009 –> 2010-2019 10.946422 9319.5904 9384.0647
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 5250.78316 5284.01627 5.2293711 2000-2009 –> 2010-2019 10.956878 9321.3117 9385.7860
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 3001.07451 3030.91829 4.1615780 2000-2009 –> 2010-2019 7.410401 5305.8326 5364.2997
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 3008.22426 3038.06805 4.1900121 2000-2009 –> 2010-2019 7.439722 5313.2231 5371.6902
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5321.86430 5355.09741 5.4521243 2000-2009 –> 2010-2019 11.174421 9391.2261 9455.7005
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 5247.85818 5281.09129 5.2204024 2000-2009 –> 2010-2019 10.950748 9319.0220 9383.4963
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 3009.38809 3039.23188 4.1946589 2000-2009 –> 2010-2019 7.455493 5317.4008 5375.8679
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + phosphate 5294.03024 5322.51576 5.3700889 2000-2009 –> 2010-2019 11.087459 9360.2879 9415.5516
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + phosphate_star 5290.55753 5319.04305 5.3591560 2000-2009 –> 2010-2019 11.099072 9361.8607 9417.1244
mod Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 2961.05611 2990.89990 4.0059533 2000-2009 –> 2010-2019 7.391533 5302.1812 5360.6483
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 5294.82128 5328.05438 5.3662803 2000-2009 –> 2010-2019 11.083534 9363.0528 9427.5271
mod Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2965.90135 2995.74514 4.0244815 2000-2009 –> 2010-2019 7.376628 5298.2731 5356.7402
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5291.44554 5324.67865 5.3556599 2000-2009 –> 2010-2019 11.087443 9362.9310 9427.4053
mod Atlantic (26.5,26.75] cstar_tref ~ temp + aou + silicate + phosphate 3003.76658 3029.34697 4.1802162 2000-2009 –> 2010-2019 7.583734 5347.5523 5397.6669
mod Atlantic (26.5,26.75] cstar_tref ~ temp + aou + silicate + phosphate_star 3010.35960 3035.93999 4.2065467 2000-2009 –> 2010-2019 7.596592 5350.6472 5400.7619
obs Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 5288.56734 5317.05286 5.3529004 2000-2009 –> 2010-2019 11.093591 9360.0436 9415.3073
mod Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 2970.57190 3000.41569 4.0424229 2000-2009 –> 2010-2019 7.389621 5301.6406 5360.1077
obs Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 5289.42948 5322.66258 5.3493271 2000-2009 –> 2010-2019 11.081773 9361.0631 9425.5374
mod Atlantic (26.75,27] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 4778.87860 4813.44659 2.4396774 2000-2009 –> 2010-2019 4.642171 8326.5246 8393.8849
mod Atlantic (26.75,27] cstar_tref ~ aou + silicate + phosphate + phosphate_star 5189.73060 5219.36031 2.9804820 2000-2009 –> 2010-2019 5.225393 8765.8978 8823.6352
mod Atlantic (26.75,27] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 4929.31153 4963.87952 2.6243171 2000-2009 –> 2010-2019 5.077154 8649.2058 8716.5661
mod Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate 5084.83503 5119.40302 2.8299083 2000-2009 –> 2010-2019 4.990717 8601.9086 8669.2689
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate 6486.58522 6521.51057 4.7770330 2000-2009 –> 2010-2019 10.559174 11610.2625 11678.0238
mod Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5110.98826 5145.55625 2.8660298 2000-2009 –> 2010-2019 5.094784 8677.5983 8744.9586
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 6489.28922 6524.21457 4.7829892 2000-2009 –> 2010-2019 10.511165 11597.8697 11665.6309
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 6493.72472 6528.65007 4.7927757 2000-2009 –> 2010-2019 10.512453 11599.9148 11667.6760
mod Atlantic (26.75,27] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5120.92250 5155.49049 2.8798710 2000-2009 –> 2010-2019 5.142029 8711.3353 8778.6955
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 6490.21152 6525.13686 4.7850225 2000-2009 –> 2010-2019 10.510612 11598.0649 11665.8262
mod Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 4786.89558 4821.46357 2.4491813 2000-2009 –> 2010-2019 4.585500 8285.7318 8353.0921
obs Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 6489.52460 6524.44994 4.7835081 2000-2009 –> 2010-2019 10.627591 11630.3573 11698.1185
mod Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 4796.18591 4830.75391 2.4602409 2000-2009 –> 2010-2019 4.610067 8305.1069 8372.4672
obs Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 6496.53023 6531.45558 4.7989761 2000-2009 –> 2010-2019 10.568013 11616.5545 11684.3158
mod Atlantic (26.75,27] cstar_tref ~ temp + aou + silicate + phosphate 5241.08889 5270.71860 3.0556492 2000-2009 –> 2010-2019 5.314963 8827.4883 8885.2257
obs Atlantic (26.75,27] cstar_tref ~ temp + aou + silicate + phosphate_star 6494.78274 6524.71875 4.7995346 2000-2009 –> 2010-2019 10.624090 11628.2267 11686.3078
mod Atlantic (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 4819.46607 4854.03406 2.4881747 2000-2009 –> 2010-2019 4.670280 8352.2317 8419.5920
obs Atlantic (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 6497.15414 6532.07949 4.8003561 2000-2009 –> 2010-2019 10.578628 11619.7536 11687.5148
obs Atlantic (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate_star 6495.55866 6525.49467 4.8012510 2000-2009 –> 2010-2019 10.604401 11623.0752 11681.1563
obs Atlantic (26.75,27] cstar_tref ~ temp + silicate + phosphate + phosphate_star 6495.83099 6525.76700 4.8018536 2000-2009 –> 2010-2019 10.624192 11628.6621 11686.7432
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 5933.54928 5970.14800 2.0728792 2000-2009 –> 2010-2019 4.251917 10888.2533 10960.0246
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 6095.47075 6132.06947 2.1983145 2000-2009 –> 2010-2019 4.578574 11248.7314 11320.5027
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate 6236.25224 6267.62257 2.3152054 2000-2009 –> 2010-2019 4.624137 11319.1188 11380.6370
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 6237.72086 6274.31958 2.3147591 2000-2009 –> 2010-2019 4.603713 11303.0516 11374.8228
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 12034.77723 12073.51309 6.0205163 2000-2009 –> 2010-2019 11.698318 21900.2564 21976.4639
mod Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + silicate 6035.05819 6071.65691 2.1506509 2000-2009 –> 2010-2019 4.365093 11025.9936 11097.7648
mod Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate 6236.33277 6272.93149 2.3135935 2000-2009 –> 2010-2019 4.587284 11286.6237 11358.3950
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 12106.60768 12145.34354 6.1372639 2000-2009 –> 2010-2019 11.744446 21933.0128 22009.2202
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 12043.07140 12081.80726 6.0338828 2000-2009 –> 2010-2019 11.736927 21922.3998 21998.6073
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 12128.73163 12167.46749 6.1736765 2000-2009 –> 2010-2019 11.831323 21983.1089 22059.3163
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 12108.43863 12147.17449 6.1402692 2000-2009 –> 2010-2019 11.756562 21939.9119 22016.1194
obs Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + phosphate_star 12055.60855 12088.81071 6.0573818 2000-2009 –> 2010-2019 11.784558 21946.1198 22011.4404
mod Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate 6198.35081 6229.72114 2.2835839 2000-2009 –> 2010-2019 4.563533 11252.8211 11314.3393
mod Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 6199.93639 6236.53511 2.2832405 2000-2009 –> 2010-2019 4.539918 11233.3426 11305.1139
obs Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 12030.56051 12069.29636 6.0137322 2000-2009 –> 2010-2019 11.794596 21952.2015 22028.4089
mod Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 6199.95941 6236.55813 2.2832596 2000-2009 –> 2010-2019 4.519170 11212.5832 11284.3545
obs Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 11984.03827 12022.77413 5.9393901 2000-2009 –> 2010-2019 11.602539 21841.4501 21917.6575
obs Atlantic (27,27.25] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 12065.67258 12098.87474 6.0737036 2000-2009 –> 2010-2019 11.829825 21971.9222 22037.2429
mod Atlantic (27,27.25] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 6107.23205 6143.83077 2.2077159 2000-2009 –> 2010-2019 4.350747 11024.4789 11096.2501
obs Atlantic (27,27.25] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 11997.07133 12035.80719 5.9601237 2000-2009 –> 2010-2019 11.647431 21867.7729 21943.9803
obs Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 14552.63685 14592.81541 5.7223835 2000-2009 –> 2010-2019 10.364240 25715.2327 25794.2104
mod Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate 12792.40947 12827.03081 3.5912077 2000-2009 –> 2010-2019 6.901767 24398.7686 24467.6107
mod Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 12563.10304 12603.49460 3.4200977 2000-2009 –> 2010-2019 6.699726 24129.8585 24210.1743
mod Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 12587.37288 12627.76444 3.4376617 2000-2009 –> 2010-2019 6.741717 24187.0164 24267.3322
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 14534.17708 14574.35564 5.6994457 2000-2009 –> 2010-2019 10.271180 25639.3260 25718.3038
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 12885.67808 12926.06965 3.6610561 2000-2009 –> 2010-2019 6.954804 24471.4739 24551.7897
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 14531.58953 14571.76809 5.6962378 2000-2009 –> 2010-2019 10.238143 25612.0334 25691.0111
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + silicate 12739.37714 12779.76870 3.5497367 2000-2009 –> 2010-2019 6.757556 24208.0141 24288.3299
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + phosphate_star 14605.92641 14640.36518 5.7916391 2000-2009 –> 2010-2019 10.347537 25695.9791 25763.6743
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 14602.99887 14643.17743 5.7854330 2000-2009 –> 2010-2019 10.363923 25713.7207 25792.6984
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 12454.41759 12494.80915 3.3425367 2000-2009 –> 2010-2019 7.042740 24555.9282 24636.2440
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 14605.69006 14645.86862 5.7888216 2000-2009 –> 2010-2019 10.343389 25696.6405 25775.6182
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 12887.53693 12927.92849 3.6624927 2000-2009 –> 2010-2019 6.836128 24308.6904 24389.0062
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 14528.85044 14569.02900 5.6928440 2000-2009 –> 2010-2019 10.241098 25614.5657 25693.5435
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 12437.77454 12478.16610 3.3308161 2000-2009 –> 2010-2019 6.956794 24449.4763 24529.7921
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + phosphate + phosphate_star 14602.75140 14637.19017 5.7876395 2000-2009 –> 2010-2019 10.348312 25696.7574 25764.4526
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 12457.88734 12498.27890 3.3449854 2000-2009 –> 2010-2019 7.044990 24559.1599 24639.4757
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 14602.62261 14642.80117 5.7849594 2000-2009 –> 2010-2019 10.343504 25696.8678 25775.8455
mod Atlantic (27.25,27.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 12809.62050 12850.01206 3.6027555 2000-2009 –> 2010-2019 6.698085 24120.0484 24200.3642
obs Atlantic (27.25,27.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 14539.47159 14579.65015 5.7060152 2000-2009 –> 2010-2019 10.365446 25716.3287 25795.3065
obs Atlantic (27.5,27.75] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 15647.66583 15688.57552 5.1824840 2000-2009 –> 2010-2019 10.114806 29215.9745 29296.9087
obs Atlantic (27.5,27.75] cstar_tref ~ aou + silicate + phosphate + phosphate_star 15655.93335 15690.99879 5.1929240 2000-2009 –> 2010-2019 10.125829 29222.7732 29292.1454
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + nitrate + phosphate 20430.73524 20469.29228 2.2648434 2000-2009 –> 2010-2019 4.656857 39197.4663 39273.9255
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 19872.19011 19917.17332 2.1299738 2000-2009 –> 2010-2019 4.494961 38547.9056 38637.1081
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate 20585.33760 20623.89464 2.3035216 2000-2009 –> 2010-2019 4.685402 39317.3167 39393.7759
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 20544.26215 20589.24537 2.2926791 2000-2009 –> 2010-2019 4.671017 39266.0586 39355.2610
obs Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 15657.71632 15698.62600 5.1927031 2000-2009 –> 2010-2019 10.064412 29170.4317 29251.3659
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 20275.46291 20320.44613 2.2261636 2000-2009 –> 2010-2019 4.614257 39030.7680 39119.9705
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 20582.18468 20627.16789 2.3022219 2000-2009 –> 2010-2019 4.674913 39284.5227 39373.7251
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 15661.61037 15702.52006 5.1966679 2000-2009 –> 2010-2019 10.077818 29183.0305 29263.9648
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 20354.35173 20399.33495 2.2454824 2000-2009 –> 2010-2019 4.598419 38988.2483 39077.4508
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 15687.08526 15727.99495 5.2226804 2000-2009 –> 2010-2019 10.134993 29237.1176 29318.0519
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 20555.38008 20600.36330 2.2954727 2000-2009 –> 2010-2019 4.681754 39304.4705 39393.6730
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 20467.55509 20512.53831 2.2734976 2000-2009 –> 2010-2019 4.630596 39115.9189 39205.1213
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 20362.02204 20407.00525 2.2473696 2000-2009 –> 2010-2019 4.599867 38994.3915 39083.5940
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 15669.45002 15710.35970 5.2046592 2000-2009 –> 2010-2019 10.100210 29204.1146 29285.0488
obs Atlantic (27.5,27.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 15652.25344 15693.16312 5.1871461 2000-2009 –> 2010-2019 10.137678 29237.1307 29318.0649
obs Atlantic (27.5,27.75] cstar_tref ~ temp + aou + silicate + phosphate 15659.96909 15695.03454 5.1970333 2000-2009 –> 2010-2019 10.147768 29243.0309 29312.4032
obs Atlantic (27.5,27.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 15660.78811 15701.69779 5.1958305 2000-2009 –> 2010-2019 10.171946 29268.8403 29349.7746
obs Atlantic (27.5,27.75] cstar_tref ~ temp + silicate + phosphate + phosphate_star 15668.31808 15703.38352 5.2055448 2000-2009 –> 2010-2019 10.181664 29274.3731 29343.7454
obs Atlantic (27.75,27.85] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 10172.21322 10210.00901 5.4062121 2000-2009 –> 2010-2019 10.256868 19534.1409 19609.3993
obs Atlantic (27.75,27.85] cstar_tref ~ aou + silicate + phosphate + phosphate_star 10170.27487 10202.67126 5.4063141 2000-2009 –> 2010-2019 10.276006 19542.4149 19606.9221
mod Atlantic (27.75,27.85] cstar_tref ~ sal + aou + nitrate + silicate 5736.35454 5768.88161 1.3416691 2000-2009 –> 2010-2019 2.566404 10918.2171 10982.9882
obs Atlantic (27.75,27.85] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 10158.19945 10195.99524 5.3830931 2000-2009 –> 2010-2019 10.257136 19535.1238 19610.3822
mod Atlantic (27.75,27.85] cstar_tref ~ sal + aou + silicate + phosphate 5494.07183 5526.59890 1.2478450 2000-2009 –> 2010-2019 2.425867 10551.9591 10616.7302
obs Atlantic (27.75,27.85] cstar_tref ~ sal + aou + silicate + phosphate 10156.20056 10188.59695 5.3830949 2000-2009 –> 2010-2019 10.279216 19545.2173 19609.7245
mod Atlantic (27.75,27.85] cstar_tref ~ sal + nitrate + silicate + phosphate_star 6160.98585 6193.51291 1.5234438 2000-2009 –> 2010-2019 2.918891 11758.8505 11823.6216
mod Atlantic (27.75,27.85] cstar_tref ~ sal + silicate + phosphate + phosphate_star 6151.63109 6184.15815 1.5191854 2000-2009 –> 2010-2019 2.918001 11757.1819 11821.9530
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + silicate 5695.82438 5733.77262 1.3247032 2000-2009 –> 2010-2019 2.549438 10879.6869 10955.2531
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 5382.17663 5420.12487 1.2060350 2000-2009 –> 2010-2019 2.371027 10406.6055 10482.1717
obs Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 10157.88372 10195.67951 5.3825734 2000-2009 –> 2010-2019 10.261659 19538.0328 19613.2912
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5522.53589 5560.48413 1.2577655 2000-2009 –> 2010-2019 2.462952 10655.1034 10730.6697
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 5527.81431 5565.76255 1.2597536 2000-2009 –> 2010-2019 2.464448 10659.0779 10734.6441
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5516.04255 5553.99080 1.2553241 2000-2009 –> 2010-2019 2.458057 10642.1113 10717.6776
obs Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 10159.45477 10197.25056 5.3851600 2000-2009 –> 2010-2019 10.274591 19546.2078 19621.4662
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate_star 5609.62286 5642.14993 1.2917443 2000-2009 –> 2010-2019 2.549508 10876.3202 10941.0913
obs Atlantic (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 10171.87211 10209.66790 5.4056482 2000-2009 –> 2010-2019 10.264447 19539.0298 19614.2882
obs Atlantic (27.75,27.85] cstar_tref ~ temp + aou + silicate + phosphate 10169.93622 10202.33261 5.4057542 2000-2009 –> 2010-2019 10.284977 19548.1728 19612.6800
obs Atlantic (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 10172.04496 10209.84074 5.4059340 2000-2009 –> 2010-2019 10.273714 19544.9607 19620.2190
obs Atlantic (27.75,27.85] cstar_tref ~ temp + silicate + phosphate + phosphate_star 10170.11273 10202.50912 5.4060460 2000-2009 –> 2010-2019 10.296233 19555.3481 19619.8553
mod Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate 6803.06130 6835.45769 1.9306150 2000-2009 –> 2010-2019 5.004644 15045.8226 15110.5564
obs Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 21442.22503 21485.27830 5.3289671 2000-2009 –> 2010-2019 10.385768 41329.6037 41415.3026
mod Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 6801.11325 6838.90904 1.9282855 2000-2009 –> 2010-2019 4.999330 15042.7294 15118.2521
obs Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 20798.33068 20841.38395 4.8561376 2000-2009 –> 2010-2019 9.631234 40310.9525 40396.6514
mod Atlantic (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate 6552.87937 6585.27576 1.7884161 2000-2009 –> 2010-2019 4.716247 14637.8623 14702.5961
obs Atlantic (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate_star 20802.49333 20839.39613 4.8604579 2000-2009 –> 2010-2019 9.636126 40313.8977 40387.3540
obs Atlantic (27.85,27.95] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 21214.28793 21257.34120 5.1565414 2000-2009 –> 2010-2019 10.101520 40955.4681 41041.1671
obs Atlantic (27.85,27.95] cstar_tref ~ sal + silicate + phosphate + phosphate_star 21254.76410 21291.66691 5.1882445 2000-2009 –> 2010-2019 10.141374 41004.7117 41078.1680
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + silicate 6611.92954 6649.72533 1.8198916 2000-2009 –> 2010-2019 4.873500 14835.1100 14910.6327
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + phosphate 7036.32486 7068.72125 2.0733653 2000-2009 –> 2010-2019 5.008488 15129.3622 15194.0959
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 6192.38262 6230.17841 1.6007546 2000-2009 –> 2010-2019 4.475450 14220.0617 14295.5844
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 20646.00263 20689.05590 4.7505596 2000-2009 –> 2010-2019 9.493384 40114.2888 40199.9878
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 6310.95201 6348.74780 1.6598627 2000-2009 –> 2010-2019 4.571774 14380.2805 14455.8033
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 20624.38409 20667.43737 4.7357631 2000-2009 –> 2010-2019 9.477559 40091.2512 40176.9502
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 20651.09020 20694.14347 4.7540484 2000-2009 –> 2010-2019 9.548870 40190.6644 40276.3633
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 6321.37423 6359.17002 1.6651615 2000-2009 –> 2010-2019 4.598611 14414.5652 14490.0879
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate 20669.67952 20706.58233 4.7681939 2000-2009 –> 2010-2019 9.565542 40210.6968 40284.1531
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 6299.37582 6337.17161 1.6539969 2000-2009 –> 2010-2019 4.559224 14361.2637 14436.7864
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 20643.50134 20686.55461 4.7488453 2000-2009 –> 2010-2019 9.490825 40110.6216 40196.3206
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate_star 6418.81461 6451.21100 1.7165767 2000-2009 –> 2010-2019 4.726845 14593.7084 14658.4421
obs Atlantic (27.95,28.05] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 21887.09216 21930.25170 5.4184030 2000-2009 –> 2010-2019 11.613169 43685.1699 43771.1652
obs Atlantic (27.95,28.05] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 21438.21855 21481.37809 5.0835224 2000-2009 –> 2010-2019 11.114933 43056.7652 43142.7605
obs Atlantic (27.95,28.05] cstar_tref ~ sal + aou + silicate + phosphate_star 21548.70387 21585.69775 5.1654464 2000-2009 –> 2010-2019 11.213498 43183.7603 43257.4705
obs Atlantic (27.95,28.05] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 21579.60858 21622.76812 5.1867103 2000-2009 –> 2010-2019 11.266128 43251.4161 43337.4114
obs Atlantic (27.95,28.05] cstar_tref ~ sal + silicate + phosphate + phosphate_star 21740.15249 21777.14637 5.3079271 2000-2009 –> 2010-2019 11.407545 43432.2445 43505.9548
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate 10350.69305 10384.13732 3.4421389 2000-2009 –> 2010-2019 7.184815 20348.1366 20414.6303
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate_star 10384.00056 10417.44483 3.4717077 2000-2009 –> 2010-2019 7.225275 20392.0388 20458.5325
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 10216.14478 10255.16309 3.3235276 2000-2009 –> 2010-2019 7.024794 20175.0235 20252.5994
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 21595.75859 21638.91813 5.1986292 2000-2009 –> 2010-2019 11.271486 43260.3113 43346.3066
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 10230.23436 10269.25268 3.3355748 2000-2009 –> 2010-2019 7.037946 20190.2013 20267.7773
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 21526.31819 21569.47772 5.1475747 2000-2009 –> 2010-2019 11.195361 43163.0804 43249.0757
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate + phosphate_star 10455.27785 10488.72212 3.5358404 2000-2009 –> 2010-2019 7.336367 20508.6469 20575.1406
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 21616.77545 21659.93498 5.2141810 2000-2009 –> 2010-2019 11.422635 43429.6811 43515.6764
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 10212.31896 10251.33728 3.3202638 2000-2009 –> 2010-2019 7.003669 20153.5605 20231.1365
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + phosphate + phosphate_star 10378.73341 10412.17768 3.4670150 2000-2009 –> 2010-2019 7.219108 20385.3392 20451.8329
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + phosphate_star 10474.20235 10502.07258 3.5548914 2000-2009 –> 2010-2019 7.361306 20531.2156 20586.6270
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate 21745.28658 21782.28047 5.3118017 2000-2009 –> 2010-2019 11.523695 43559.9124 43633.6227
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 10230.20326 10269.22158 3.3355482 2000-2009 –> 2010-2019 7.037896 20190.1472 20267.7232
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 21575.13213 21618.29166 5.1834114 2000-2009 –> 2010-2019 11.246927 43229.3432 43315.3385
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate_star 10228.23485 10261.67912 3.3355752 2000-2009 –> 2010-2019 7.038207 20186.4587 20252.9524
obs Atlantic (28.05,28.1] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 15796.88267 15837.81154 5.2913102 2000-2009 –> 2010-2019 11.063624 30776.4334 30857.7242
obs Atlantic (28.05,28.1] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 15738.85469 15779.78356 5.2316329 2000-2009 –> 2010-2019 11.010497 30723.7556 30805.0464
obs Atlantic (28.05,28.1] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 15742.05204 15782.98090 5.2349035 2000-2009 –> 2010-2019 10.997027 30713.2664 30794.5573
obs Atlantic (28.05,28.1] cstar_tref ~ sal + silicate + phosphate + phosphate_star 15840.24274 15875.32462 5.3384333 2000-2009 –> 2010-2019 11.124496 30829.0175 30898.6953
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 15863.85041 15904.77928 5.3610281 2000-2009 –> 2010-2019 11.113329 30827.0152 30908.3060
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate 3781.06746 3811.18860 1.3036715 2000-2009 –> 2010-2019 2.350725 6725.8341 6785.4077
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate_star 3825.23878 3855.35992 1.3296577 2000-2009 –> 2010-2019 2.399705 6813.4944 6873.0681
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 3780.87083 3816.01216 1.3023926 2000-2009 –> 2010-2019 2.349291 6727.3409 6796.8435
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 3821.89059 3857.03193 1.3264839 2000-2009 –> 2010-2019 2.396461 6812.0157 6881.5184
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 15777.17126 15818.10013 5.2709626 2000-2009 –> 2010-2019 10.977696 30702.8130 30784.1039
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 15856.58649 15897.51535 5.3534216 2000-2009 –> 2010-2019 11.108887 30822.3460 30903.6368
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + phosphate_star 3874.91026 3905.03141 1.3594988 2000-2009 –> 2010-2019 2.461182 6921.4980 6981.0716
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + phosphate_star 15864.93078 15900.01267 5.3642569 2000-2009 –> 2010-2019 11.120156 30829.0454 30898.7232
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 3843.90390 3879.04524 1.3395958 2000-2009 –> 2010-2019 2.431320 6874.3111 6943.8137
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 15740.73727 15781.66614 5.2335584 2000-2009 –> 2010-2019 10.961735 30684.0740 30765.3648
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + phosphate + phosphate_star 3813.53807 3843.65921 1.3227241 2000-2009 –> 2010-2019 2.387091 6791.1382 6850.7119
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 3812.21045 3847.35178 1.3207588 2000-2009 –> 2010-2019 2.385125 6791.8085 6861.3111
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 15823.81429 15864.74315 5.3192382 2000-2009 –> 2010-2019 11.054002 30772.5733 30853.8641
mod Atlantic (28.05,28.1] cstar_tref ~ temp + aou + phosphate 3800.68674 3825.78770 1.3163262 2000-2009 –> 2010-2019 2.430136 6867.1895 6916.8343
mod Atlantic (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate 3802.68644 3832.80758 1.3163260 2000-2009 –> 2010-2019 2.426565 6864.7620 6924.3357
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 17845.40954 17886.94612 5.9100517 2000-2009 –> 2010-2019 13.618707 33712.0827 33793.7734
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + phosphate 4996.04481 5028.69976 1.0417955 2000-2009 –> 2010-2019 2.347915 10551.6094 10616.6937
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + phosphate_star 5049.82506 5082.48002 1.0583367 2000-2009 –> 2010-2019 2.373854 10628.9624 10694.0467
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + silicate + phosphate 4964.13363 5002.23108 1.0314986 2000-2009 –> 2010-2019 2.247694 10287.1563 10363.0880
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5028.09030 5066.18775 1.0510045 2000-2009 –> 2010-2019 2.282709 10392.7753 10468.7070
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 17864.37955 17905.91612 5.9301779 2000-2009 –> 2010-2019 13.661007 33744.2076 33825.8984
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 17827.81097 17869.34755 5.8914415 2000-2009 –> 2010-2019 13.598436 33693.4967 33775.1874
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + phosphate_star 5073.52383 5106.17879 1.0657089 2000-2009 –> 2010-2019 2.377658 10643.7324 10708.8167
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 5074.17080 5112.26825 1.0652866 2000-2009 –> 2010-2019 2.360947 10605.2999 10681.2316
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 17841.91454 17883.45112 5.9063511 2000-2009 –> 2010-2019 13.626772 33715.5728 33797.2635
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + phosphate + phosphate_star 5034.04780 5066.70276 1.0534571 2000-2009 –> 2010-2019 2.363022 10598.2758 10663.3601
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + phosphate + phosphate_star 17867.53795 17903.14073 5.9356626 2000-2009 –> 2010-2019 13.666024 33745.0889 33815.1095
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5006.47406 5044.57151 1.0443709 2000-2009 –> 2010-2019 2.265101 10341.7324 10417.6641
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 17864.93038 17906.46696 5.9307634 2000-2009 –> 2010-2019 13.661123 33744.4807 33826.1714
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + silicate + phosphate_star 17863.05946 17898.66224 5.9309006 2000-2009 –> 2010-2019 13.662634 33741.4236 33811.4442
obs Atlantic (28.1,28.15] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 17847.21050 17888.74708 5.9119595 2000-2009 –> 2010-2019 13.640990 33725.9733 33807.6640
mod Atlantic (28.1,28.15] cstar_tref ~ temp + aou + silicate + phosphate 5105.43581 5138.09077 1.0757172 2000-2009 –> 2010-2019 2.339294 10552.1213 10617.2055
obs Atlantic (28.1,28.15] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 17861.24465 17896.84743 5.9289719 2000-2009 –> 2010-2019 13.663597 33741.3210 33811.3416
obs Atlantic (28.1,28.15] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 17836.18184 17877.71842 5.9002862 2000-2009 –> 2010-2019 13.628642 33714.5446 33796.2354
mod Atlantic (28.1,28.15] cstar_tref ~ temp + silicate + phosphate + phosphate_star 5160.40303 5193.05799 1.0931770 2000-2009 –> 2010-2019 2.363081 10623.5116 10688.5959
mod Atlantic (28.15,28.2] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 7641.76060 7683.96435 0.8384273 2000-2009 –> 2010-2019 1.846449 15657.9753 15741.7507
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate 7678.41684 7714.59148 0.8437242 2000-2009 –> 2010-2019 1.839748 15625.4823 15697.2898
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 6504.60340 6539.27919 5.3692114 2000-2009 –> 2010-2019 11.964491 12990.3665 13059.2481
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate + silicate 7545.28844 7587.49219 0.8253525 2000-2009 –> 2010-2019 1.816283 15465.6030 15549.3784
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + phosphate_star 6505.13165 6534.85376 5.3756980 2000-2009 –> 2010-2019 11.980877 12991.8316 13050.8729
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate 6443.72937 6478.40516 5.2153717 2000-2009 –> 2010-2019 11.814302 12930.5761 12999.4577
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 6386.82322 6421.49901 5.0755482 2000-2009 –> 2010-2019 11.600338 12851.5467 12920.4282
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 6506.27929 6540.95508 5.3735103 2000-2009 –> 2010-2019 11.975151 12993.9301 13062.8116
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + nitrate + phosphate_star 8000.59490 8036.76954 0.8891934 2000-2009 –> 2010-2019 1.913003 16101.9615 16173.7690
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 7959.06630 8001.27005 0.8829099 2000-2009 –> 2010-2019 1.899820 16024.5108 16108.2862
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 6407.61976 6442.29555 5.1262071 2000-2009 –> 2010-2019 11.720165 12892.9904 12961.8719
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 6409.46023 6444.13602 5.1307146 2000-2009 –> 2010-2019 11.668505 12878.0813 12946.9628
mod Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + phosphate 7653.34421 7689.51885 0.8402848 2000-2009 –> 2010-2019 1.835490 15595.7999 15667.6074
mod Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + phosphate_star 7689.10328 7725.27792 0.8451944 2000-2009 –> 2010-2019 1.845209 15658.5932 15730.4007
mod Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate 7567.84101 7604.01565 0.8286606 2000-2009 –> 2010-2019 1.827690 15531.8026 15603.6101
obs Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 6423.03434 6457.71013 5.1640820 2000-2009 –> 2010-2019 11.885506 12945.8936 13014.7751
obs Atlantic (28.15,28.2] cstar_tref ~ temp + aou + silicate + phosphate_star 6435.46334 6465.18544 5.1997886 2000-2009 –> 2010-2019 11.933230 12959.8200 13018.8613
mod Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 7744.82210 7780.99674 0.8529018 2000-2009 –> 2010-2019 1.858576 15745.9586 15817.7660
obs Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 6436.97023 6471.64602 5.1985643 2000-2009 –> 2010-2019 11.949167 12968.3109 13037.1924
mod Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + silicate + phosphate_star 7989.30776 8025.48241 0.8875598 2000-2009 –> 2010-2019 1.904995 16055.6501 16127.4576
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + phosphate 11253.94830 11293.55199 0.6806407 2000-2009 –> 2010-2019 1.355373 20372.3799 20450.3741
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 11168.84756 11215.05186 0.6752086 2000-2009 –> 2010-2019 1.349911 20288.8853 20379.8785
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 12721.07328 12761.62839 3.3236932 2000-2009 –> 2010-2019 6.462349 24892.6791 24973.6347
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + phosphate_star 11174.84495 11214.44864 0.6757055 2000-2009 –> 2010-2019 1.350737 20297.2173 20375.2114
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate 11172.46679 11212.07048 0.6755577 2000-2009 –> 2010-2019 1.336390 20106.0572 20184.0513
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 12715.02857 12755.58367 3.3195534 2000-2009 –> 2010-2019 6.458418 24886.9496 24967.9051
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 10576.65576 10622.86006 0.6394077 2000-2009 –> 2010-2019 1.282698 19273.3346 19364.3278
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + phosphate + phosphate_star 12724.93987 12759.70139 3.3277160 2000-2009 –> 2010-2019 6.466983 24895.4696 24964.8600
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + silicate + phosphate 12714.80715 12749.56867 3.3207709 2000-2009 –> 2010-2019 6.459646 24884.7442 24954.1346
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 12708.67382 12749.22893 3.3152067 2000-2009 –> 2010-2019 6.451440 24876.6169 24957.5724
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate 11227.71776 11267.32145 0.6790003 2000-2009 –> 2010-2019 1.354269 20353.2109 20431.2051
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate 11176.84356 11223.04787 0.6757055 2000-2009 –> 2010-2019 1.350435 20297.2413 20388.2345
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate 12714.85839 12755.41350 3.3194369 2000-2009 –> 2010-2019 6.461299 24891.3081 24972.2636
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 11171.68722 11217.89152 0.6753850 2000-2009 –> 2010-2019 1.350287 20294.3557 20385.3489
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + silicate 11095.67485 11141.87915 0.6706786 2000-2009 –> 2010-2019 1.327839 19981.7861 20072.7793
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + silicate + phosphate 12699.16306 12739.71817 3.3087120 2000-2009 –> 2010-2019 6.447402 24870.8210 24951.7766
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 12721.66882 12762.22393 3.3241014 2000-2009 –> 2010-2019 6.443396 24863.9195 24944.8751
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 11223.12567 11269.32997 0.6785886 2000-2009 –> 2010-2019 1.356786 20389.0458 20480.0390
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 12736.76032 12777.31542 3.3344609 2000-2009 –> 2010-2019 6.454145 24879.6046 24960.5602
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 12699.94943 12740.50454 3.3092486 2000-2009 –> 2010-2019 6.429014 24842.9174 24923.8729
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 25735.97929 25780.23264 5.5177690 2000-2009 –> 2010-2019 11.410848 43482.3306 43568.0879
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 34006.85478 34052.98137 5.7155326 2000-2009 –> 2010-2019 11.915892 57646.8426 57736.3712
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + phosphate + phosphate_star 26081.11664 26119.04809 5.7556019 2000-2009 –> 2010-2019 11.817393 43982.2384 44055.7446
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 25881.65839 25925.91175 5.6163568 2000-2009 –> 2010-2019 11.511739 43630.1795 43715.9368
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 34370.22654 34416.35314 5.9120315 2000-2009 –> 2010-2019 12.370456 58307.2419 58396.7705
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate 34716.83875 34756.37583 6.1068956 2000-2009 –> 2010-2019 12.482821 58558.2101 58634.9489
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 25796.68370 25840.93706 5.5586384 2000-2009 –> 2010-2019 11.350430 43446.7456 43532.5028
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 34350.91918 34397.04578 5.9014229 2000-2009 –> 2010-2019 12.132849 58027.3118 58116.8405
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 25896.69738 25940.95073 5.6266342 2000-2009 –> 2010-2019 11.486814 43611.9549 43697.7122
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 34396.61072 34442.73731 5.9265595 2000-2009 –> 2010-2019 12.168076 58084.7890 58174.3176
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 26426.48591 26470.73927 6.0009367 2000-2009 –> 2010-2019 12.203192 44456.8370 44542.5942
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 34668.15115 34714.27775 6.0781686 2000-2009 –> 2010-2019 12.451974 58509.0998 58598.6284
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 26331.88251 26376.13586 5.9323181 2000-2009 –> 2010-2019 11.983054 44224.8666 44310.6238
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 34747.20076 34793.32736 6.1230289 2000-2009 –> 2010-2019 12.687927 58803.3214 58892.8500
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 26538.22680 26582.48015 6.0830088 2000-2009 –> 2010-2019 12.297091 44579.1590 44664.9162
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 34843.22819 34889.35479 6.1779696 2000-2009 –> 2010-2019 12.780206 58940.6591 59030.1877
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 25998.65254 26042.90589 5.6968061 2000-2009 –> 2010-2019 11.628788 43781.5476 43867.3048
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 34504.34727 34550.47386 5.9862541 2000-2009 –> 2010-2019 12.273130 58245.2690 58334.7976
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 26596.14599 26640.39934 6.1259904 2000-2009 –> 2010-2019 12.377110 44670.0830 44755.8402
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 34879.66824 34925.79484 6.1989470 2000-2009 –> 2010-2019 12.816548 58994.0316 59083.5603
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 19876.39094 19919.46640 4.2128411 2000-2009 –> 2010-2019 8.407165 35368.9819 35453.3981
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 20086.71250 20129.78796 4.3422413 2000-2009 –> 2010-2019 8.614533 35679.2411 35763.6572
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 21480.67654 21523.97984 4.8067285 2000-2009 –> 2010-2019 9.502931 37930.5461 38015.3381
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 20259.00587 20302.08133 4.4512009 2000-2009 –> 2010-2019 8.824036 35977.7498 36062.1660
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 21486.87708 21530.18038 4.8108802 2000-2009 –> 2010-2019 9.519278 37951.1194 38035.9114
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + phosphate 21853.98656 21891.10367 5.0645928 2000-2009 –> 2010-2019 10.076047 38661.9297 38734.6086
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 19664.57354 19707.64900 4.0864177 2000-2009 –> 2010-2019 8.050366 34850.6417 34935.0579
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 21478.88285 21522.18615 4.8055282 2000-2009 –> 2010-2019 9.629820 38077.8864 38162.6785
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 19867.85574 19910.93120 4.2076720 2000-2009 –> 2010-2019 8.323030 35257.3190 35341.7351
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 21530.30953 21573.61284 4.8400616 2000-2009 –> 2010-2019 9.701876 38172.2509 38257.0429
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 20308.58711 20351.66257 4.4830601 2000-2009 –> 2010-2019 8.891200 36070.9624 36155.3785
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 21488.78980 21532.09310 4.8121616 2000-2009 –> 2010-2019 9.525317 37958.6289 38043.4210
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate_star 21566.43694 21603.55406 4.8658246 2000-2009 –> 2010-2019 9.774790 38259.8663 38332.5452
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 21394.95483 21438.25813 4.7496982 2000-2009 –> 2010-2019 9.613451 38039.1047 38123.8967
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 19906.50501 19949.58047 4.2311295 2000-2009 –> 2010-2019 8.382212 35342.8670 35427.2831
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 21529.08798 21572.39128 4.8392385 2000-2009 –> 2010-2019 9.703005 38173.2537 38258.0457
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate_star 21534.12823 21571.24535 4.8439845 2000-2009 –> 2010-2019 9.731131 38202.8706 38275.5495
mod Indo-Pacific (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 20242.22820 20285.30367 4.4404715 2000-2009 –> 2010-2019 8.638956 35740.1986 35824.6147
mod Indo-Pacific (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 20327.00766 20370.08312 4.4949545 2000-2009 –> 2010-2019 8.710544 35847.0392 35931.4553
mod Indo-Pacific (26,26.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 20375.98576 20419.06122 4.5267341 2000-2009 –> 2010-2019 8.769414 35930.7751 36015.1912
obs Indo-Pacific (26.5,26.75] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 21591.70886 21628.45198 5.9239468 2000-2009 –> 2010-2019 11.463218 39674.6954 39747.2443
mod Indo-Pacific (26.5,26.75] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 18778.98142 18821.51485 4.4709268 2000-2009 –> 2010-2019 8.723316 32867.1009 32950.2755
mod Indo-Pacific (26.5,26.75] cstar_tref ~ aou + silicate + phosphate + phosphate_star 18913.80351 18950.26074 4.5670213 2000-2009 –> 2010-2019 9.051760 33261.1315 33332.4240
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 19006.98920 19049.52263 4.6322081 2000-2009 –> 2010-2019 8.981066 33205.2515 33288.4261
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 21227.89045 21270.75743 5.6113511 2000-2009 –> 2010-2019 10.810363 38946.9643 39031.6047
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 19247.12298 19289.65641 4.8083611 2000-2009 –> 2010-2019 9.279870 33581.9460 33665.1206
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 21261.98448 21304.85146 5.6397740 2000-2009 –> 2010-2019 10.892520 39040.4083 39125.0487
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + phosphate 21574.49798 21611.24111 5.9088569 2000-2009 –> 2010-2019 11.476281 39686.7453 39759.2942
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 18572.98511 18615.51854 4.3300491 2000-2009 –> 2010-2019 8.545721 32618.5258 32701.7004
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 21576.45819 21619.32517 5.9088221 2000-2009 –> 2010-2019 11.471309 39685.5850 39770.2254
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 18775.76074 18818.29417 4.4686893 2000-2009 –> 2010-2019 8.806055 32961.0310 33044.2056
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 21612.11375 21654.98073 5.9401261 2000-2009 –> 2010-2019 11.532171 39751.8306 39836.4710
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 21266.06419 21308.93117 5.6431847 2000-2009 –> 2010-2019 10.905378 39054.8596 39139.5000
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate_star 21300.39151 21337.13463 5.6736463 2000-2009 –> 2010-2019 11.066085 39228.3111 39300.8600
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 18923.37923 18965.91266 4.5724019 2000-2009 –> 2010-2019 9.045337 33259.7678 33342.9424
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 21291.29464 21334.16162 5.6643238 2000-2009 –> 2010-2019 10.994090 39153.7380 39238.3784
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 18788.69357 18831.22701 4.4776808 2000-2009 –> 2010-2019 8.825375 32985.6420 33068.8166
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 21604.69955 21647.56653 5.9336032 2000-2009 –> 2010-2019 11.522436 39741.0998 39825.7401
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate_star 18921.51253 18957.96976 4.5724967 2000-2009 –> 2010-2019 9.059044 33270.8201 33342.1126
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 19011.19242 19053.72586 4.6352353 2000-2009 –> 2010-2019 8.991437 33217.7386 33300.9132
mod Indo-Pacific (26.75,27] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 24186.72607 24224.87244 4.1225581 2000-2009 –> 2010-2019 8.029058 44348.4940 44423.8091
obs Indo-Pacific (26.75,27] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 31055.84395 31094.92787 5.4488893 2000-2009 –> 2010-2019 10.716158 58101.8779 58179.2817
mod Indo-Pacific (26.75,27] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 23834.22090 23878.72500 3.9546597 2000-2009 –> 2010-2019 7.836629 43952.3570 44040.2245
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 31434.19299 31479.79091 5.6585490 2000-2009 –> 2010-2019 11.178048 58892.7249 58983.0293
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + phosphate 23655.16211 23699.66620 3.8724720 2000-2009 –> 2010-2019 7.536425 43354.5966 43442.4642
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + phosphate 31304.30877 31349.90669 5.5852958 2000-2009 –> 2010-2019 10.943140 58501.9708 58592.2751
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 23892.53087 23937.03497 3.9817986 2000-2009 –> 2010-2019 7.725117 43747.2022 43835.0697
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate 31532.66606 31578.26398 5.7147264 2000-2009 –> 2010-2019 11.152326 58860.0021 58950.3064
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 23913.98356 23958.48766 3.9918300 2000-2009 –> 2010-2019 7.740100 43778.2318 43866.0994
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + phosphate_star 24201.41782 24239.56419 4.1296681 2000-2009 –> 2010-2019 7.918327 44141.3045 44216.6195
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 24065.89580 24110.39990 4.0635920 2000-2009 –> 2010-2019 7.843056 43990.1804 44078.0479
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 31335.40719 31381.00511 5.6027481 2000-2009 –> 2010-2019 10.935309 58491.5580 58581.8623
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + phosphate 31383.81124 31422.89517 5.6311506 2000-2009 –> 2010-2019 11.062016 58698.2722 58775.6759
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 23959.93443 24004.43853 4.0134021 2000-2009 –> 2010-2019 7.902300 44090.9878 44178.8554
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 31236.99377 31282.59168 5.5477047 2000-2009 –> 2010-2019 10.874994 58384.4644 58474.7688
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 24200.70692 24245.21102 4.1283552 2000-2009 –> 2010-2019 8.118970 44518.7976 44606.6651
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 31390.21496 31435.81288 5.6336390 2000-2009 –> 2010-2019 11.040265 58667.4182 58757.7225
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 24225.19615 24269.70025 4.1402301 2000-2009 –> 2010-2019 8.138158 44556.5491 44644.4167
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 31386.22095 31431.81886 5.6313821 2000-2009 –> 2010-2019 11.039481 58665.8158 58756.1201
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate_star 31388.40121 31427.48514 5.6337442 2000-2009 –> 2010-2019 11.041940 58666.1519 58743.5556
mod Indo-Pacific (27,27.25] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 26051.49202 26090.93120 2.8382622 2000-2009 –> 2010-2019 5.927219 49972.5886 50050.7520
mod Indo-Pacific (27,27.25] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 25895.57074 25941.58310 2.7961960 2000-2009 –> 2010-2019 5.737245 49358.0229 49449.2135
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 37724.90174 37772.07064 4.9715772 2000-2009 –> 2010-2019 9.882117 72411.7963 72505.5748
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 37907.81958 37954.98848 5.0450054 2000-2009 –> 2010-2019 10.071545 72863.6351 72957.4136
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate 37733.58556 37774.01605 4.9758365 2000-2009 –> 2010-2019 9.886379 72418.4850 72498.8666
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate_star 37908.47307 37948.90356 5.0460785 2000-2009 –> 2010-2019 10.081812 72883.3362 72963.7178
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 37842.69786 37889.86677 5.0187403 2000-2009 –> 2010-2019 10.005537 72707.0828 72800.8613
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + silicate + phosphate + phosphate_star 37843.47943 37883.90992 5.0198593 2000-2009 –> 2010-2019 10.019166 72734.7214 72815.1030
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate 26195.12783 26241.14019 2.8765284 2000-2009 –> 2010-2019 5.939360 50038.4849 50129.6756
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 26653.12175 26699.13412 3.0038331 2000-2009 –> 2010-2019 6.161950 50784.0905 50875.2811
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate 37707.07400 37754.24290 4.9644781 2000-2009 –> 2010-2019 9.871768 72386.3412 72480.1197
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 37898.46737 37945.63628 5.0412250 2000-2009 –> 2010-2019 10.006312 72712.5991 72806.3777
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 26687.75340 26733.76577 3.0136854 2000-2009 –> 2010-2019 6.179967 50842.9607 50934.1514
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 38033.70138 38080.87028 5.0961667 2000-2009 –> 2010-2019 10.078768 72888.3920 72982.1705
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 37845.31930 37892.48821 5.0197949 2000-2009 –> 2010-2019 9.974766 72635.9603 72729.7388
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + phosphate 26435.13214 26474.57131 2.9431091 2000-2009 –> 2010-2019 6.076441 50490.1342 50568.2976
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + phosphate_star 26846.26399 26885.70316 3.0597726 2000-2009 –> 2010-2019 6.274321 51141.5036 51219.6670
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 26387.75140 26433.76377 2.9293994 2000-2009 –> 2010-2019 5.973137 50172.4015 50263.5922
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 26835.68779 26881.70016 3.0561363 2000-2009 –> 2010-2019 6.212825 50962.4116 51053.6023
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 26861.01059 26907.02296 3.0634625 2000-2009 –> 2010-2019 6.224708 51001.2735 51092.4642
obs Indo-Pacific (27.25,27.5] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 40201.32110 40242.38106 4.3998083 2000-2009 –> 2010-2019 8.908091 76136.3871 76217.7835
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 40180.31206 40228.21535 4.3925081 2000-2009 –> 2010-2019 8.850825 75980.4977 76075.4602
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 35251.29780 35299.02511 3.2835120 2000-2009 –> 2010-2019 6.337675 65320.4510 65414.9892
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 39613.94929 39661.85258 4.2165846 2000-2009 –> 2010-2019 8.617389 75254.6650 75349.6275
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + silicate + phosphate 39620.71482 39661.77477 4.2192530 2000-2009 –> 2010-2019 8.627237 75279.4507 75360.8472
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 39848.59819 39896.50147 4.2885998 2000-2009 –> 2010-2019 8.784399 75751.6094 75846.5719
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + silicate + phosphate + phosphate_star 39848.53574 39889.59569 4.2891996 2000-2009 –> 2010-2019 8.819518 75843.4901 75924.8865
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 35225.90583 35273.63313 3.2773474 2000-2009 –> 2010-2019 6.447557 65737.1199 65831.6582
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 35524.53607 35572.26337 3.3505868 2000-2009 –> 2010-2019 6.588021 66284.4892 66379.0274
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 39615.35919 39663.26248 4.2170137 2000-2009 –> 2010-2019 8.571082 75124.9431 75219.9056
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 39836.91650 39884.81979 4.2849857 2000-2009 –> 2010-2019 8.688074 75484.0041 75578.9665
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 35548.83556 35596.56286 3.3566177 2000-2009 –> 2010-2019 6.600301 66331.6462 66426.1844
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 34966.63729 35014.36460 3.2150610 2000-2009 –> 2010-2019 6.243620 64935.9955 65030.5337
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 39817.32346 39865.22675 4.2789308 2000-2009 –> 2010-2019 8.679901 75458.4999 75553.4624
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 39776.23486 39824.13814 4.2662609 2000-2009 –> 2010-2019 8.660508 75398.6349 75493.5974
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + phosphate 35266.41000 35307.31912 3.2876730 2000-2009 –> 2010-2019 6.516324 65992.1606 66073.1934
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 34831.63394 34879.36124 3.1830981 2000-2009 –> 2010-2019 6.290710 65106.4423 65200.9806
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 35159.00271 35206.73002 3.2611601 2000-2009 –> 2010-2019 6.454479 65756.3148 65850.8530
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 35159.00352 35206.73083 3.2611603 2000-2009 –> 2010-2019 6.455945 65761.7552 65856.2934
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + nitrate + silicate + phosphate_star 35157.19467 35198.10379 3.2612064 2000-2009 –> 2010-2019 6.457957 65765.2401 65846.2728
obs Indo-Pacific (27.5,27.75] cstar_tref ~ aou + silicate + phosphate + phosphate_star 44847.11838 44889.01752 4.0329986 2000-2009 –> 2010-2019 7.884651 82828.4063 82911.3062
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 39749.60117 39798.13682 3.3237923 2000-2009 –> 2010-2019 6.279729 74777.1295 74873.6359
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 44563.86770 44612.75002 3.9614511 2000-2009 –> 2010-2019 7.926673 82945.8521 83042.5687
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 39968.67029 40017.20594 3.3721521 2000-2009 –> 2010-2019 6.331075 75010.3263 75106.8326
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate 44577.67769 44619.57682 3.9653832 2000-2009 –> 2010-2019 7.954929 83041.5692 83124.4691
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate_star 40377.80651 40419.40850 3.4648184 2000-2009 –> 2010-2019 6.553921 76019.8065 76102.5262
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 39913.33350 39961.86915 3.3598704 2000-2009 –> 2010-2019 6.314333 74933.8858 75030.3921
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + silicate + phosphate + phosphate_star 40373.24019 40414.84218 3.4637753 2000-2009 –> 2010-2019 6.563482 76063.1817 76145.9014
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 40295.62545 40344.16110 3.4456397 2000-2009 –> 2010-2019 6.499755 75780.2710 75876.7773
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 44314.63245 44363.51478 3.8999769 2000-2009 –> 2010-2019 7.643695 81907.9561 82004.6727
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 40369.58230 40418.11795 3.4624833 2000-2009 –> 2010-2019 6.530313 75916.9062 76013.4125
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 44425.31192 44474.19424 3.9271575 2000-2009 –> 2010-2019 7.699951 82124.7800 82221.4965
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 40398.27588 40446.81153 3.4690404 2000-2009 –> 2010-2019 6.517912 75858.8778 75955.3841
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 39893.91465 39942.45030 3.3555711 2000-2009 –> 2010-2019 6.357487 75137.3826 75233.8889
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 44411.13043 44460.01276 3.9236643 2000-2009 –> 2010-2019 7.686626 82074.7986 82171.5151
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 40338.29042 40386.82607 3.4553466 2000-2009 –> 2010-2019 6.515681 75851.3956 75947.9019
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 44409.97227 44458.85459 3.9233791 2000-2009 –> 2010-2019 7.694295 82102.6122 82199.3287
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate_star 44713.87961 44755.77874 3.9994198 2000-2009 –> 2010-2019 7.791430 82481.0533 82563.9532
obs Indo-Pacific (27.5,27.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 44951.61869 45000.50101 4.0590224 2000-2009 –> 2010-2019 7.946967 83063.5763 83160.2929
obs Indo-Pacific (27.5,27.75] cstar_tref ~ temp + aou + silicate + phosphate 44949.97549 44991.87463 4.0591133 2000-2009 –> 2010-2019 7.959028 83102.1094 83185.0093
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate 14331.62218 14367.63232 2.6614068 2000-2009 –> 2010-2019 5.033467 26820.8470 26892.3360
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + phosphate 14039.70947 14081.72130 2.5335963 2000-2009 –> 2010-2019 4.805930 26296.1604 26379.5643
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 14076.38356 14118.39539 2.5492030 2000-2009 –> 2010-2019 4.830282 26353.8311 26437.2351
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + silicate 14305.56820 14347.58003 2.6489340 2000-2009 –> 2010-2019 5.014297 26781.3378 26864.7417
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + phosphate 14378.78412 14414.79426 2.6825076 2000-2009 –> 2010-2019 5.019255 26786.0242 26857.5133
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 14176.80761 14218.81944 2.5924323 2000-2009 –> 2010-2019 4.926685 26580.2078 26663.6117
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 24597.84655 24642.61647 3.8845814 2000-2009 –> 2010-2019 7.457663 43819.1916 43907.2179
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 14215.10815 14257.11997 2.6091120 2000-2009 –> 2010-2019 4.951713 26638.0238 26721.4277
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 24576.23057 24621.00049 3.8751114 2000-2009 –> 2010-2019 7.449419 43800.0232 43888.0495
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 14104.50034 14146.51217 2.5612332 2000-2009 –> 2010-2019 4.848847 26397.5849 26480.9888
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 24569.66468 24614.43460 3.8722394 2000-2009 –> 2010-2019 7.460670 43821.5911 43909.6174
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 14231.51769 14273.52952 2.6162910 2000-2009 –> 2010-2019 4.961633 26660.8240 26744.2279
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 24575.46492 24620.23485 3.8747764 2000-2009 –> 2010-2019 7.449431 43799.9514 43887.9777
mod Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 14193.38680 14235.39863 2.5996393 2000-2009 –> 2010-2019 5.041096 26842.2277 26925.6316
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 24590.26224 24635.03216 3.8812561 2000-2009 –> 2010-2019 7.454084 43811.1023 43899.1286
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 24573.26622 24618.03614 3.8738145 2000-2009 –> 2010-2019 7.447521 43795.8589 43883.8852
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + silicate + phosphate 24606.24846 24644.62268 3.8891468 2000-2009 –> 2010-2019 7.462295 43825.7273 43901.1784
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + silicate + phosphate_star 24583.85950 24622.23372 3.8793270 2000-2009 –> 2010-2019 7.453709 43805.8012 43881.2523
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 24573.57537 24618.34530 3.8739497 2000-2009 –> 2010-2019 7.448046 43796.9476 43884.9739
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + silicate + phosphate + phosphate_star 24582.82085 24621.19507 3.8788720 2000-2009 –> 2010-2019 7.453610 43805.4729 43880.9240
obs Indo-Pacific (27.85,27.95] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 31344.06574 31390.96817 3.2861820 2000-2009 –> 2010-2019 6.225536 55315.5166 55407.7492
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 16283.03577 16326.57483 2.1625595 2000-2009 –> 2010-2019 4.203579 30119.1754 30205.2999
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate 15756.34125 15793.66044 2.0150723 2000-2009 –> 2010-2019 3.968390 29305.7906 29379.6116
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + phosphate 15747.70783 15791.24688 2.0121897 2000-2009 –> 2010-2019 3.961097 29284.5036 29370.6282
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 15740.86477 15784.40382 2.0103368 2000-2009 –> 2010-2019 3.960618 29282.2281 29368.3526
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + silicate 15744.82976 15788.36881 2.0114102 2000-2009 –> 2010-2019 3.942257 29221.2781 29307.4027
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + phosphate 16432.75474 16470.07393 2.2071845 2000-2009 –> 2010-2019 4.281663 30372.2944 30446.1154
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 16295.10722 16338.64627 2.1660767 2000-2009 –> 2010-2019 4.238144 30229.1078 30315.2323
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 31315.29162 31362.19405 3.2783182 2000-2009 –> 2010-2019 6.221034 55297.7083 55389.9410
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 16312.65363 16356.19268 2.1711995 2000-2009 –> 2010-2019 4.249488 30266.0871 30352.2117
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 31305.08392 31351.98635 3.2755330 2000-2009 –> 2010-2019 6.228452 55320.7063 55412.9390
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 15794.22260 15837.76166 2.0248297 2000-2009 –> 2010-2019 3.985273 29369.2729 29455.3974
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 16331.44851 16374.98756 2.1767002 2000-2009 –> 2010-2019 4.257673 30293.2490 30379.3735
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 31305.03196 31351.93439 3.2755189 2000-2009 –> 2010-2019 6.229905 55325.4215 55417.6542
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 31299.32306 31346.22549 3.2739622 2000-2009 –> 2010-2019 6.216813 55282.1784 55374.4111
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 31296.67308 31343.57551 3.2732399 2000-2009 –> 2010-2019 6.226438 55313.2019 55405.4346
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + silicate + phosphate 31318.31001 31358.51209 3.2796883 2000-2009 –> 2010-2019 6.222793 55299.9928 55379.0494
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + silicate + phosphate_star 31310.30447 31350.50655 3.2775029 2000-2009 –> 2010-2019 6.230740 55324.9617 55404.0182
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 31296.87702 31343.77946 3.2732955 2000-2009 –> 2010-2019 6.227783 55317.5941 55409.8267
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + silicate + phosphate + phosphate_star 31310.49522 31350.69731 3.2775550 2000-2009 –> 2010-2019 6.232172 55329.6347 55408.6913
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 52928.89065 52972.18322 3.3651773 2000-2009 –> 2010-2019 6.587433 94473.9511 94559.1826
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 52930.64464 52981.15263 3.3651362 2000-2009 –> 2010-2019 6.587019 94475.8487 94575.2854
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + phosphate + phosphate_star 53067.81444 53103.89157 3.3888516 2000-2009 –> 2010-2019 6.619207 94651.1514 94722.1776
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + silicate + phosphate + phosphate_star 53068.58039 53111.87295 3.3886435 2000-2009 –> 2010-2019 6.618903 94653.4382 94738.6696
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 52930.17379 52980.68178 3.3650573 2000-2009 –> 2010-2019 6.587303 94477.1855 94576.6222
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + aou + phosphate + phosphate_star 53061.08780 53104.38037 3.3873807 2000-2009 –> 2010-2019 6.617245 94643.9837 94729.2151
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate 20495.10865 20534.35059 1.7905768 2000-2009 –> 2010-2019 3.493732 34327.0718 34403.3467
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate + phosphate 53157.05491 53207.56290 3.4032522 2000-2009 –> 2010-2019 6.661144 94880.5582 94979.9949
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate + silicate 20364.34815 20410.13041 1.7674985 2000-2009 –> 2010-2019 3.429208 34023.8434 34112.8308
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate 20875.33255 20914.57449 1.8583533 2000-2009 –> 2010-2019 3.682311 35192.6001 35268.8750
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate_star 20879.03920 20918.28114 1.8590265 2000-2009 –> 2010-2019 3.687693 35214.5654 35290.8404
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 20873.33128 20919.11355 1.8576269 2000-2009 –> 2010-2019 3.657680 35099.1696 35188.1571
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 53198.62015 53249.12814 3.4102964 2000-2009 –> 2010-2019 6.673713 94949.3051 95048.7418
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 20877.91134 20923.69360 1.8584584 2000-2009 –> 2010-2019 3.664790 35128.4077 35217.3951
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 21157.78524 21203.56751 1.9099839 2000-2009 –> 2010-2019 3.761540 35583.4052 35672.3926
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + phosphate + phosphate_star 20888.00852 20927.25046 1.8606565 2000-2009 –> 2010-2019 3.691525 35232.0562 35308.3311
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 20886.61081 20932.39307 1.8600389 2000-2009 –> 2010-2019 3.668099 35143.8802 35232.8676
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate_star 21244.83032 21284.07226 1.9266750 2000-2009 –> 2010-2019 3.797643 35742.3140 35818.5889
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + nitrate + phosphate 53166.54469 53209.83725 3.4051979 2000-2009 –> 2010-2019 6.666337 94904.0296 94989.2611
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 53076.56149 53127.06948 3.3896519 2000-2009 –> 2010-2019 6.634048 94733.4734 94832.9101
obs Indo-Pacific (28.05,28.1] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 38106.86005 38148.48809 2.9506841 2000-2009 –> 2010-2019 5.699848 62882.9593 62963.8035
obs Indo-Pacific (28.05,28.1] cstar_tref ~ aou + phosphate + phosphate_star 38169.62813 38204.31817 2.9632575 2000-2009 –> 2010-2019 5.713450 62947.5390 63014.9091
obs Indo-Pacific (28.05,28.1] cstar_tref ~ aou + silicate + phosphate + phosphate_star 38159.00032 38200.62836 2.9608018 2000-2009 –> 2010-2019 5.709678 62934.0326 63014.8767
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 38132.40203 38180.96807 2.9552481 2000-2009 –> 2010-2019 5.718269 62961.7345 63056.0527
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + aou + silicate + phosphate 38193.12253 38234.75057 2.9674420 2000-2009 –> 2010-2019 5.736010 63040.8913 63121.7354
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + silicate 10688.10269 10729.93886 1.5126152 2000-2009 –> 2010-2019 3.167793 18797.2877 18878.6884
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 10894.42709 10936.26325 1.5671625 2000-2009 –> 2010-2019 3.147854 18809.7600 18891.1607
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 38112.49250 38161.05855 2.9513879 2000-2009 –> 2010-2019 5.680258 62815.0897 62909.4079
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 10900.67112 10942.50728 1.5688436 2000-2009 –> 2010-2019 3.154378 18828.8809 18910.2815
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 38255.02076 38303.58681 2.9791341 2000-2009 –> 2010-2019 5.730020 63039.5001 63133.8182
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 10852.39720 10894.23337 1.5558935 2000-2009 –> 2010-2019 3.141614 18781.1006 18862.5012
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 38241.43419 38290.00024 2.9764780 2000-2009 –> 2010-2019 5.722445 63007.6756 63101.9938
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 10900.11248 10941.94865 1.5686931 2000-2009 –> 2010-2019 3.154304 18828.5276 18909.9283
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 38254.07631 38302.64235 2.9789494 2000-2009 –> 2010-2019 5.729351 63036.7612 63131.0794
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate_star 10903.08968 10938.94925 1.5700343 2000-2009 –> 2010-2019 3.155811 18829.9448 18899.7168
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate_star 38253.02104 38294.64908 2.9791341 2000-2009 –> 2010-2019 5.730192 63036.1380 63116.9821
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate 11241.50838 11277.36795 1.6639679 2000-2009 –> 2010-2019 3.266946 19213.7841 19283.5561
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate_star 11250.74562 11286.60519 1.6666092 2000-2009 –> 2010-2019 3.275848 19239.4308 19309.2028
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + nitrate + silicate + phosphate_star 11235.22616 11271.08573 1.6621740 2000-2009 –> 2010-2019 3.270998 19222.8279 19292.5999
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + silicate + phosphate + phosphate_star 11249.01434 11284.87391 1.6661138 2000-2009 –> 2010-2019 3.275428 19237.8976 19307.6696
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 42216.46305 42265.96060 2.7359819 2000-2009 –> 2010-2019 5.273673 72054.4153 72151.2008
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 22277.77985 22325.08763 1.3917103 2000-2009 –> 2010-2019 2.864287 40940.0462 41033.1989
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 42128.70309 42178.20064 2.7222172 2000-2009 –> 2010-2019 5.257932 71956.7700 72053.5555
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + phosphate + phosphate_star 42251.00094 42293.42741 2.7417331 2000-2009 –> 2010-2019 5.281512 72097.3892 72180.3482
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + silicate + phosphate 42159.43880 42201.86527 2.7273435 2000-2009 –> 2010-2019 5.264518 71992.8045 72075.7635
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 22362.24345 22409.55123 1.4009779 2000-2009 –> 2010-2019 2.889470 41135.5104 41228.6631
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 42248.93214 42298.42968 2.7410921 2000-2009 –> 2010-2019 5.312123 72252.5184 72349.3039
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + silicate + phosphate + phosphate_star 42319.74397 42362.17044 2.7525865 2000-2009 –> 2010-2019 5.327023 72338.1290 72421.0880
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + nitrate + silicate 20833.12389 20880.43168 1.2423607 2000-2009 –> 2010-2019 2.770685 39879.0788 39972.2315
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate 20513.03955 20560.34733 1.2115026 2000-2009 –> 2010-2019 2.603937 38597.4618 38690.6145
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate 42143.75243 42193.24998 2.7245727 2000-2009 –> 2010-2019 5.241325 71876.5651 71973.3506
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 20535.69345 20583.00124 1.2136611 2000-2009 –> 2010-2019 2.608222 38635.8703 38729.0230
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 42240.32863 42289.82618 2.7397371 2000-2009 –> 2010-2019 5.276080 72071.5345 72168.3200
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 21547.89606 21595.20384 1.3141363 2000-2009 –> 2010-2019 2.785400 40200.9487 40294.1014
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 20468.40521 20515.71300 1.2072609 2000-2009 –> 2010-2019 2.610849 38635.2114 38728.3641
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 42344.71294 42394.21049 2.7562224 2000-2009 –> 2010-2019 5.308796 72256.8664 72353.6519
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate 21582.73700 21623.28654 1.3179462 2000-2009 –> 2010-2019 2.888162 40905.9197 40985.7649
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 20531.96830 20579.27609 1.2133059 2000-2009 –> 2010-2019 2.607572 38629.9639 38723.1166
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 42220.85972 42270.35727 2.7366733 2000-2009 –> 2010-2019 5.269110 72032.5086 72129.2941
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate_star 20684.59142 20725.14095 1.2281377 2000-2009 –> 2010-2019 2.635415 38876.5020 38956.3472
obs Indo-Pacific (28.2,28.3] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 13092.61063 13128.63682 2.1501383 2000-2009 –> 2010-2019 4.760201 21343.0864 21411.8547
obs Indo-Pacific (28.2,28.3] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 12996.08445 13038.11501 2.1150496 2000-2009 –> 2010-2019 4.723130 21245.9283 21326.1581
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 13007.41222 13049.44278 2.1190545 2000-2009 –> 2010-2019 4.726617 21256.5680 21336.7978
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 36175.64309 36229.58560 0.7279004 2000-2009 –> 2010-2019 1.399852 60490.0582 60595.6885
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 39835.37087 39889.31338 0.8137277 2000-2009 –> 2010-2019 1.594678 67726.3411 67831.9714
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + phosphate + phosphate_star 13133.39411 13169.42031 2.1648326 2000-2009 –> 2010-2019 4.774031 21382.7222 21451.4906
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 13121.79882 13163.82937 2.1599231 2000-2009 –> 2010-2019 4.768986 21372.9474 21453.1772
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + nitrate + silicate + phosphate 36823.20429 36869.44073 0.7424440 2000-2009 –> 2010-2019 1.498060 63927.5389 64018.0790
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 36527.57910 36581.52161 0.7357445 2000-2009 –> 2010-2019 1.412853 61023.9139 61129.5441
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + nitrate + phosphate 12990.72994 13032.76050 2.1131591 2000-2009 –> 2010-2019 4.771332 21306.4746 21386.7044
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 13110.92787 13152.95843 2.1560054 2000-2009 –> 2010-2019 4.764238 21360.9742 21441.2039
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + silicate + phosphate 39801.84607 39855.78858 0.8128972 2000-2009 –> 2010-2019 1.654263 69465.6564 69571.2866
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 39821.68727 39875.62978 0.8133886 2000-2009 –> 2010-2019 1.655105 69495.4338 69601.0640
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 36696.14256 36750.08507 0.7395314 2000-2009 –> 2010-2019 1.488049 63577.9255 63683.5557
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 39818.33482 39872.27733 0.8133056 2000-2009 –> 2010-2019 1.654978 69490.8274 69596.4576
obs Indo-Pacific (28.2,28.3] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 13043.29191 13085.32247 2.1317898 2000-2009 –> 2010-2019 4.741777 21295.6677 21375.8974
mod Indo-Pacific (28.2,28.3] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 40217.03372 40270.97623 0.8232416 2000-2009 –> 2010-2019 1.662926 69833.2493 69938.8795
obs Indo-Pacific (28.2,28.3] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 13015.18977 13057.22033 2.1218086 2000-2009 –> 2010-2019 4.769773 21317.6055 21397.8353
obs Indo-Pacific (28.2,28.3] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 13140.49796 13182.52852 2.1666786 2000-2009 –> 2010-2019 4.776228 21392.2924 21472.5222
mod Indo-Pacific (28.2,28.3] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 40231.88082 40285.82333 0.8236140 2000-2009 –> 2010-2019 1.664487 69881.7779 69987.4081
mod Indo-Pacific (28.3,28.4] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 23714.05693 23755.30614 1.2694046 2000-2009 –> 2010-2019 2.196126 36751.5724 36831.7420
mod Indo-Pacific (28.3,28.4] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 23711.67563 23759.79971 1.2690157 2000-2009 –> 2010-2019 2.193381 36726.4901 36820.0213
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 23527.49310 23575.61717 1.2527757 2000-2009 –> 2010-2019 2.179041 36562.2342 36655.7654
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 4735.89637 4770.84747 2.1150067 2000-2009 –> 2010-2019 4.763250 8013.7232 8080.3495
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 4732.66017 4767.61128 2.1118664 2000-2009 –> 2010-2019 4.746371 8003.3925 8070.0188
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + silicate + phosphate_star 4735.79524 4765.75333 2.1168514 2000-2009 –> 2010-2019 4.767409 8012.8139 8069.9222
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 4752.69112 4787.64222 2.1313787 2000-2009 –> 2010-2019 4.764816 8022.8709 8089.4972
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + silicate + phosphate + phosphate_star 4752.34091 4782.29900 2.1329938 2000-2009 –> 2010-2019 4.767903 8021.2828 8078.3910
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + nitrate + phosphate 23499.64979 23547.77386 1.2503389 2000-2009 –> 2010-2019 2.173868 36505.6878 36599.2190
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 23577.81423 23625.93830 1.2571920 2000-2009 –> 2010-2019 2.185691 36635.9118 36729.4430
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + silicate + phosphate 4713.07591 4748.02702 2.0929619 2000-2009 –> 2010-2019 4.726490 7983.3024 8049.9287
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4707.02275 4741.97385 2.0871532 2000-2009 –> 2010-2019 4.736865 7985.6061 8052.2324
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 23500.66407 23548.78814 1.2504275 2000-2009 –> 2010-2019 2.175040 36518.0784 36611.6096
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 4731.71919 4766.67030 2.1109542 2000-2009 –> 2010-2019 4.760955 8010.4511 8077.0774
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + silicate + phosphate 4729.83332 4759.79141 2.1110648 2000-2009 –> 2010-2019 4.761954 8007.0223 8064.1305
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4709.72119 4744.67229 2.0897407 2000-2009 –> 2010-2019 4.724104 7980.3801 8047.0064
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + aou + nitrate + phosphate 23689.57671 23730.82592 1.2672334 2000-2009 –> 2010-2019 2.191216 36698.3819 36778.5515
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 23688.54515 23736.66923 1.2669647 2000-2009 –> 2010-2019 2.188261 36671.1003 36764.6315
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 23693.90163 23735.15083 1.2676167 2000-2009 –> 2010-2019 2.192737 36714.6372 36794.8069
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 23692.72080 23740.84487 1.2673347 2000-2009 –> 2010-2019 2.189827 36687.8593 36781.3905
obs Indo-Pacific (28.4, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate 59.49593 63.96933 1.2286501 2000-2009 –> 2010-2019 2.873999 173.1743 186.4543
obs Indo-Pacific (28.4, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 59.30114 63.77455 1.2201325 2000-2009 –> 2010-2019 2.872167 173.1903 186.4704
obs Indo-Pacific (28.4, Inf] cstar_tref ~ sal + temp + aou + silicate + phosphate 57.48219 61.95559 1.1433890 2000-2009 –> 2010-2019 2.550503 163.0271 176.3072
obs Indo-Pacific (28.4, Inf] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 57.51401 61.98742 1.1446895 2000-2009 –> 2010-2019 2.531283 162.2950 175.5751
obs Indo-Pacific (28.4, Inf] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 57.43569 61.90909 1.1414919 2000-2009 –> 2010-2019 2.596529 164.7221 178.0022
obs Indo-Pacific (28.4, Inf] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 57.31795 61.79135 1.1367019 2000-2009 –> 2010-2019 2.652283 166.7243 180.0044
obs Indo-Pacific (28.4, Inf] cstar_tref ~ sal + temp + silicate + phosphate 55.61097 59.44531 1.1486600 2000-2009 –> 2010-2019 2.668207 163.1532 174.5361
obs Indo-Pacific (28.4, Inf] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 57.43984 61.91324 1.1416612 2000-2009 –> 2010-2019 2.538619 162.6081 175.8881
obs Indo-Pacific (28.4, Inf] cstar_tref ~ sal + temp + silicate + phosphate_star 55.58749 59.42183 1.1476971 2000-2009 –> 2010-2019 2.686342 163.7792 175.1621
obs Indo-Pacific (28.4, Inf] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 63.85038 68.32378 1.4353835 2000-2009 –> 2010-2019 2.878122 170.6954 183.9755

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

[[2]]

Version Author Date
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
d3c3fde jens-daniel-mueller 2021-03-26
594d503 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26

8.1.2 Best models

left_join(lm_best_target %>% select(basin, gamma_slab, era, model),
          GLODAP_glanced) %>%
  group_split(era) %>%
  imap(
  ~ ggplot(data = .x,
         aes("",
             !!sym(params_local$MLR_criterion),
             col = data_source)) +
  geom_hline(yintercept = c(0, 10)) +
  geom_boxplot() +
  facet_grid(gamma_slab ~ basin) +
  scale_color_brewer(palette = "Set1") +
  ylim(c(0, NA)) +
  labs(x = "Number of predictors pooled",
           title = paste("Era:", unique(.x$era)))
  )
[[1]]

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

[[2]]

Version Author Date
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
d3c3fde jens-daniel-mueller 2021-03-26
594d503 jens-daniel-mueller 2021-03-26

8.2 RMSE correlation between eras

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

8.2.1 All models

# find max rmse to scale axis
max_rmse <-
  max(c(lm_all_eras$rmse,
        lm_all_eras$rmse_sum - lm_all_eras$rmse))

lm_all_eras %>%
  ggplot(aes(rmse, rmse_sum - rmse, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  coord_equal(xlim = c(0, max_rmse),
              ylim = c(0, max_rmse)) +
  geom_abline(slope = 1,
              col = 'red') +
  facet_grid(data_source ~ basin)

Version Author Date
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
ac1a836 jens-daniel-mueller 2021-02-24
b03fbd3 jens-daniel-mueller 2021-02-24
3d3b4cc jens-daniel-mueller 2021-02-23
7b672f7 jens-daniel-mueller 2021-01-11
33ba23c jens-daniel-mueller 2021-01-07
318609d jens-daniel-mueller 2020-12-23
9d0b2d0 jens-daniel-mueller 2020-12-23
0aa2b50 jens-daniel-mueller 2020-12-23
2886da0 jens-daniel-mueller 2020-12-19
02f0ee9 jens-daniel-mueller 2020-12-18
7bcb4eb jens-daniel-mueller 2020-12-18
7131186 jens-daniel-mueller 2020-12-17
5b48ef5 jens-daniel-mueller 2020-12-17
e4ca289 jens-daniel-mueller 2020-12-16
158fe26 jens-daniel-mueller 2020-12-15
7a9a4cb jens-daniel-mueller 2020-12-15
61b263c jens-daniel-mueller 2020-12-15
984697e jens-daniel-mueller 2020-12-12
3ebff89 jens-daniel-mueller 2020-12-12
ba112d3 jens-daniel-mueller 2020-12-11
24a632f jens-daniel-mueller 2020-12-07
6a8004b jens-daniel-mueller 2020-12-07
70bf1a5 jens-daniel-mueller 2020-12-07
7555355 jens-daniel-mueller 2020-12-07
143d6fa jens-daniel-mueller 2020-12-07
090e4d5 jens-daniel-mueller 2020-12-02
7c25f7a jens-daniel-mueller 2020-12-02
b03ddb8 jens-daniel-mueller 2020-12-02
91435ae jens-daniel-mueller 2020-12-01
196be51 jens-daniel-mueller 2020-11-30 
rm(max_rmse)

8.2.2 Best models

# find max rmse to scale axis
max_rmse <-
  max(c(lm_best_eras$rmse,
        lm_best_eras$rmse_sum - lm_best_eras$rmse))

lm_best_eras %>%
  ggplot(aes(rmse, rmse_sum - rmse, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  coord_equal(xlim = c(0, max_rmse),
              ylim = c(0, max_rmse)) +
  geom_abline(slope = 1,
              col = 'red') +
  facet_grid(data_source ~ basin)

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

8.3 Predictor counts

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

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

# print table
lm_all_stats %>%
  gt(rowname_col = "gamma_slab",
     groupname_col = c("basin", "data_source")) %>% 
  summary_rows(
    groups = TRUE,
    fns = list(total = "sum")
  )
aou nitrate phosphate phosphate_star sal silicate temp
Atlantic - mod
(-Inf,26] 5 7 6 8 10 8 5
(26,26.5] 7 6 6 7 10 6 8
(26.5,26.75] 8 5 7 6 4 10 8
(26.75,27] 8 5 8 6 4 10 7
(27,27.25] 9 7 6 4 6 10 6
(27.25,27.5] 6 8 4 7 9 8 7
(27.5,27.75] 7 5 8 5 10 7 6
(27.75,27.85] 5 4 4 6 10 10 6
(27.85,27.95] 7 4 4 5 10 9 7
(27.95,28.05] 4 2 4 8 10 5 10
(28.05,28.1] 6 2 6 6 8 5 10
(28.1,28.15] 5 2 6 7 8 6 10
(28.15,28.2] 6 10 2 6 5 5 9
(28.2, Inf] 9 10 4 5 10 3 5
total 92.00 77.00 75.00 86.00 114.00 102.00 104.00
Indo-Pacific - mod
(-Inf,26] 8 5 7 7 10 5 7
(26,26.5] 7 6 7 7 7 7 9
(26.5,26.75] 7 5 6 7 7 8 8
(26.75,27] 6 10 6 8 5 5 8
(27,27.25] 8 10 7 7 3 4 8
(27.25,27.5] 6 10 6 6 5 6 9
(27.5,27.75] 5 5 6 7 10 10 5
(27.75,27.85] 8 6 6 4 9 5 10
(27.85,27.95] 8 6 5 5 10 5 9
(27.95,28.05] 6 3 4 6 10 6 10
(28.05,28.1] 5 3 4 7 6 10 10
(28.1,28.2] 4 5 6 5 10 10 8
(28.2,28.3] 5 7 8 5 8 10 6
(28.3,28.4] 7 10 9 7 4 3 7
total 90.00 91.00 87.00 88.00 104.00 94.00 114.00
Atlantic - obs
(-Inf,26] 8 8 6 6 8 5 8
(26,26.5] 8 10 8 7 4 4 7
(26.5,26.75] 7 9 6 7 4 4 10
(26.75,27] 5 5 5 8 4 10 10
(27,27.25] 5 8 5 9 5 6 10
(27.25,27.5] 6 5 7 8 9 4 9
(27.5,27.75] 7 3 9 7 4 10 7
(27.75,27.85] 7 4 10 5 4 10 6
(27.85,27.95] 5 4 7 6 10 10 5
(27.95,28.05] 5 4 7 6 10 10 5
(28.05,28.1] 4 7 5 9 10 7 6
(28.1,28.15] 3 6 5 10 7 6 10
(28.15,28.2] 6 5 4 9 7 7 10
(28.2, Inf] 8 4 9 6 10 6 5
total 84.00 82.00 93.00 103.00 96.00 99.00 108.00
Indo-Pacific - obs
(-Inf,26] 8 6 6 6 10 5 8
(26,26.5] 5 5 5 7 10 5 10
(26.5,26.75] 6 6 6 7 9 4 9
(26.75,27] 7 9 7 5 4 7 8
(27,27.25] 6 4 6 7 10 10 4
(27.25,27.5] 6 5 8 7 9 8 4
(27.5,27.75] 7 3 7 5 7 10 7
(27.75,27.85] 6 4 6 7 4 10 10
(27.85,27.95] 7 4 7 7 3 10 9
(27.95,28.05] 10 6 10 6 4 4 4
(28.05,28.1] 7 3 7 7 7 8 5
(28.1,28.2] 6 4 8 7 10 8 4
(28.2,28.3] 10 6 10 8 5 4 5
(28.3,28.4] 5 3 7 7 10 10 5
(28.4, Inf] 5 5 5 6 9 8 10
total 101.00 73.00 105.00 99.00 111.00 111.00 102.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
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
ac1a836 jens-daniel-mueller 2021-02-24
b03fbd3 jens-daniel-mueller 2021-02-24
3d3b4cc jens-daniel-mueller 2021-02-23
7b672f7 jens-daniel-mueller 2021-01-11
33ba23c jens-daniel-mueller 2021-01-07
318609d jens-daniel-mueller 2020-12-23
9d0b2d0 jens-daniel-mueller 2020-12-23
0aa2b50 jens-daniel-mueller 2020-12-23
2886da0 jens-daniel-mueller 2020-12-19
02f0ee9 jens-daniel-mueller 2020-12-18
7bcb4eb jens-daniel-mueller 2020-12-18
7131186 jens-daniel-mueller 2020-12-17
5b48ef5 jens-daniel-mueller 2020-12-17 
lm_best_eras %>% 
  ggplot(aes(rmse, AIC, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  facet_grid(data_source~basin)

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

8.4.2 AIC vs BIC

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

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

Version Author Date
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20

8.5 RMSE vs VIF

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

Version Author Date
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20

8.6 Residual patterns

8.6.1 Fitted vs actual models

Plotted are fitted vs actual target variable values, here: rparams_local$MLR_target`

GLODAP_augmented_best <- left_join(
  lm_best_target %>% select(data_source, basin, gamma_slab, era, model),
  GLODAP_augmented
)
# calculate equal axis limits and binwidth
axis_lims <- GLODAP_augmented %>%
  summarise(
    max_value = max(
      c(max(.fitted, max(!!sym(params_local$MLR_target))))
      ),
    min_value = min(
      c(min(.fitted, min(!!sym(params_local$MLR_target))))
    )
  )

i_binwidth <- 1
# binwidth_value <- (axis_lims$max_value - axis_lims$min_value) / 40
axis_lims <- c(axis_lims$min_value, axis_lims$max_value)

GLODAP_augmented %>%
  ggplot(aes(!!sym(params_local$MLR_target), .fitted)) +
  geom_bin2d(binwidth = i_binwidth) +
  scale_fill_viridis_c() +
  geom_abline(slope = 1,
              col = 'red') +
  coord_equal(xlim = axis_lims,
              ylim = axis_lims) +
  labs(title = "All models") +
  facet_wrap(~data_source)

Version Author Date
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
27ae473 jens-daniel-mueller 2021-02-24
4bc00ea jens-daniel-mueller 2021-02-24 
GLODAP_augmented_best %>%
  ggplot(aes(!!sym(params_local$MLR_target), .fitted)) +
  geom_bin2d(binwidth = i_binwidth) +
  scale_fill_viridis_c() +
  geom_abline(slope = 1,
              col = 'red') +
  coord_equal(xlim = axis_lims,
              ylim = axis_lims) +
  labs(title = "Selected models") +
  facet_wrap(~data_source)

Version Author Date
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
27ae473 jens-daniel-mueller 2021-02-24
4bc00ea jens-daniel-mueller 2021-02-24 
rm(binwidth_value, axis_lims)

8.6.2 Pooled

In the following, we present residual patterns vs latitude across all domains.

i_ylim <- c(-30,30)

GLODAP_augmented_best %>%
  ggplot(aes(lat, .resid)) +
  geom_bin2d(binwidth = i_binwidth) +
  geom_hline(yintercept = 0, col = "white") +
  scale_fill_viridis_c() +
  labs(
    title = paste(
      "Target variable:",
      params_local$MLR_target,
      "| Selected models",
      "| All domains"
    )
  ) +
  facet_wrap(~data_source)

Version Author Date
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
27ae473 jens-daniel-mueller 2021-02-24
4bc00ea jens-daniel-mueller 2021-02-24

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

GLODAP_augmented_best %>%
  ggplot(aes(lat, .resid)) +
  geom_bin2d(binwidth = i_binwidth) +
  geom_hline(yintercept = 0, col = "white") +
  scale_fill_viridis_c() +
  coord_cartesian(ylim = i_ylim) +
  labs(
    title = paste(
      "Target variable:",
      params_local$MLR_target,
      "| Selected models",
      "| All domains"
    )
  ) +
  facet_wrap(~data_source)

Version Author Date
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
d3c3fde jens-daniel-mueller 2021-03-26
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
27ae473 jens-daniel-mueller 2021-02-24
4bc00ea jens-daniel-mueller 2021-02-24

8.6.3 By model domain

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

GLODAP_augmented_best %>% 
  group_split(basin, era) %>% 
  map(~ ggplot(data = .x, aes(lat, .resid)) +
      geom_bin2d(binwidth = i_binwidth) +
      geom_hline(yintercept = 0, col = "black") +
      scale_fill_viridis_c() +
      facet_grid(gamma_slab ~ data_source) +
      coord_cartesian(ylim = i_ylim) +
      labs(
        title = paste(
          "Target:",
          params_local$MLR_target,
          "| best models | basin:",
          unique(.x$basin),
          "| era:",
          unique(.x$era)
        )
      ))
[[1]]

Version Author Date
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
d3c3fde jens-daniel-mueller 2021-03-26
594d503 jens-daniel-mueller 2021-03-26
06bfaf7 jens-daniel-mueller 2021-03-26
2eb881a jens-daniel-mueller 2021-03-26
03b6009 jens-daniel-mueller 2021-03-23
de6256b jens-daniel-mueller 2021-03-23
f155edd jens-daniel-mueller 2021-03-23
380d215 jens-daniel-mueller 2021-03-21
33b385b jens-daniel-mueller 2021-03-20
108d6b2 jens-daniel-mueller 2021-03-20
83a13de jens-daniel-mueller 2021-03-20
cf98c6d jens-daniel-mueller 2021-03-16
a1d52ff jens-daniel-mueller 2021-03-15
0bade3b jens-daniel-mueller 2021-03-15
27c1f4b jens-daniel-mueller 2021-03-14
af75ebf jens-daniel-mueller 2021-03-14
5017709 jens-daniel-mueller 2021-03-11
585b07f jens-daniel-mueller 2021-03-11
85a5ed2 jens-daniel-mueller 2021-03-10
6c0bec6 jens-daniel-mueller 2021-03-05
3c2ec33 jens-daniel-mueller 2021-03-05
af70b94 jens-daniel-mueller 2021-03-04
27ae473 jens-daniel-mueller 2021-02-24
a1ba577 jens-daniel-mueller 2021-02-24
071743d jens-daniel-mueller 2021-02-24 

[[2]]

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

[[3]]

Version Author Date
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
d3c3fde jens-daniel-mueller 2021-03-26
594d503 jens-daniel-mueller 2021-03-26 

[[4]]

Version Author Date
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
d3c3fde jens-daniel-mueller 2021-03-26
594d503 jens-daniel-mueller 2021-03-26

8.6.4 Latitudinal mean

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

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

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

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

[[2]]

Version Author Date
4b7a5ee jens-daniel-mueller 2021-05-28
25bd183 jens-daniel-mueller 2021-05-26
62bd574 jens-daniel-mueller 2021-05-20
7c56c39 jens-daniel-mueller 2021-05-19
52e7583 jens-daniel-mueller 2021-05-12
969e631 jens-daniel-mueller 2021-05-12
d2a83bc jens-daniel-mueller 2021-04-16
c0a47df jens-daniel-mueller 2021-04-16
50290e8 jens-daniel-mueller 2021-04-16
a00ec94 jens-daniel-mueller 2021-04-16
b6fe355 jens-daniel-mueller 2021-04-16
81b7c6d jens-daniel-mueller 2021-04-16
ddec5b7 jens-daniel-mueller 2021-04-15
29edae5 jens-daniel-mueller 2021-04-14
099d566 jens-daniel-mueller 2021-04-14
bb44686 jens-daniel-mueller 2021-04-14
bf40480 jens-daniel-mueller 2021-04-13
9f31fe3 jens-daniel-mueller 2021-04-13
338dd3c jens-daniel-mueller 2021-04-09
a79ca2c jens-daniel-mueller 2021-04-09
594d503 jens-daniel-mueller 2021-03-26

8.6.5 Latitudinal offset

# calculate residual offset for adjacent eras
lat_residual_offset <- lat_residual  %>%
  select(data_source, basin, gamma_slab, era, lat_grid, .resid_mean) %>% 
  arrange(era) %>% 
  group_by(basin, gamma_slab, lat_grid) %>% 
  mutate(eras = paste(lag(era), era, sep = " --> "),
         .resid_mean_offset = .resid_mean - lag(.resid_mean)
         ) %>% 
  ungroup() %>% 
  select(-era) %>% 
  drop_na() %>% 
  filter(eras != paste(unique(lat_residual$era)[1],
                       unique(lat_residual$era)[3],
                       sep = " --> "))

lat_residual_offset %>%
  ggplot(aes(lat_grid, .resid_mean_offset, col=eras)) +
  geom_line() +
  geom_point() +
  geom_hline(yintercept = 0, col = "black") +
  facet_grid(gamma_slab ~ basin)

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

9 VIF

GLODAP_glanced %>% 
  ggplot(aes(as.factor(n),
             log10(vif_max),
             col=data_source)) + 
  geom_hline(yintercept = log10(params_local$vif_max)) +
  scale_color_brewer(palette = "Set1") +
  geom_violin() +
  facet_grid(gamma_slab~basin)

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