Last updated: 2021-05-27

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

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 dc8e4e1. 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
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 Required data

Required are:

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

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

3 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

4 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")

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

6 Prepare coeffcients

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

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

6.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 855.93651 875.90072 6.4874916 2000-2009 –> 2010-2019 14.061879 1241.8573 1275.7444
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 846.14064 866.10485 6.2439365 2000-2009 –> 2010-2019 14.476974 1241.0667 1274.9538
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 859.06429 876.17647 6.6187488 2000-2009 –> 2010-2019 14.244629 1243.7168 1272.7629
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 856.29891 876.26312 6.4966819 2000-2009 –> 2010-2019 14.120021 1242.9154 1276.8025
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 856.33854 876.30276 6.4976879 2000-2009 –> 2010-2019 14.120171 1242.9429 1276.8300
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 846.05404 866.01825 6.2418246 2000-2009 –> 2010-2019 13.796482 1231.6931 1265.5802
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
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 856.54163 876.50584 6.5028446 2000-2009 –> 2010-2019 14.153839 1243.5492 1277.4363
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 853.00265 872.96686 6.4135668 2000-2009 –> 2010-2019 14.426117 1244.9970 1278.8841
obs Atlantic (-Inf,26] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 852.50209 872.46630 6.4010386 2000-2009 –> 2010-2019 14.406223 1244.3971 1278.2842
obs Atlantic (-Inf,26] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 850.49061 870.45483 6.3509406 2000-2009 –> 2010-2019 14.336086 1242.1149 1276.0020
obs Atlantic (26,26.5] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 3915.67054 3942.18073 5.8423438 2000-2009 –> 2010-2019 11.642731 7333.2679 7385.4828
obs Atlantic (26,26.5] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 3917.62772 3948.55628 5.8421398 2000-2009 –> 2010-2019 11.602194 7329.7449 7390.6624
mod Atlantic (26,26.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 2210.91637 2238.99524 3.5728405 2000-2009 –> 2010-2019 7.471804 3798.1979 3851.7949
obs Atlantic (26,26.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 3838.95896 3869.88752 5.4790412 2000-2009 –> 2010-2019 11.274879 7257.7152 7318.6327
mod Atlantic (26,26.5] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2178.03498 2206.11385 3.4317321 2000-2009 –> 2010-2019 7.734211 3821.0565 3874.6535
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 2214.28557 2242.36444 3.5876230 2000-2009 –> 2010-2019 7.847795 3851.7139 3905.3109
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 3856.26642 3887.19497 5.5569374 2000-2009 –> 2010-2019 11.534414 7308.1030 7369.0205
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 2214.59707 2242.67594 3.5889928 2000-2009 –> 2010-2019 7.904835 3859.3738 3912.9708
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 3858.42581 3889.35437 5.5667336 2000-2009 –> 2010-2019 11.648452 7328.7958 7389.7133
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + silicate 2238.11440 2266.19327 3.6939334 2000-2009 –> 2010-2019 8.325392 3922.8391 3976.4361
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 2181.44058 2209.51945 3.4460845 2000-2009 –> 2010-2019 8.191343 3879.9044 3933.5014
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2182.12221 2210.20108 3.4489643 2000-2009 –> 2010-2019 8.296411 3892.6453 3946.2423
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 2215.62678 2243.70565 3.5935246 2000-2009 –> 2010-2019 7.977876 3869.3176 3922.9146
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 3868.43461 3899.36317 5.6123652 2000-2009 –> 2010-2019 11.817885 7360.4076 7421.3251
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 2174.96895 2203.04782 3.4188619 2000-2009 –> 2010-2019 8.375699 3898.1228 3951.7198
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 2182.48250 2210.56137 3.4504875 2000-2009 –> 2010-2019 8.320517 3895.6363 3949.2333
obs Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + phosphate 3917.59541 3944.10560 5.8515237 2000-2009 –> 2010-2019 11.833571 7368.2513 7420.4663
obs Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 3919.56328 3950.49183 5.8513704 2000-2009 –> 2010-2019 11.725613 7352.7240 7413.6415
obs Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 3921.46009 3952.38864 5.8604303 2000-2009 –> 2010-2019 11.810201 7368.3163 7429.2338
obs Atlantic (26,26.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 3924.27684 3955.20539 5.8739102 2000-2009 –> 2010-2019 11.917860 7387.9688 7448.8863
mod Atlantic (26.5,26.75] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 2916.31233 2946.08913 3.9422294 2000-2009 –> 2010-2019 7.511114 5303.9432 5362.3434
mod Atlantic (26.5,26.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 3001.30660 3031.08340 4.2779403 2000-2009 –> 2010-2019 7.566306 5316.7351 5375.1352
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 5546.20710 5579.80049 5.2845578 2000-2009 –> 2010-2019 10.894338 9992.4717 10057.9725
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 5543.81946 5577.41285 5.2775292 2000-2009 –> 2010-2019 10.890335 9990.8444 10056.3452
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 2969.10837 2998.88517 4.1475252 2000-2009 –> 2010-2019 7.396349 5273.8665 5332.2666
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2976.24524 3006.02205 4.1760850 2000-2009 –> 2010-2019 7.425795 5281.2441 5339.6442
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5641.65156 5675.24495 5.5733203 2000-2009 –> 2010-2019 11.184126 10088.1740 10153.6747
obs Atlantic (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 5539.06154 5572.65493 5.2635511 2000-2009 –> 2010-2019 10.879367 9986.8423 10052.3431
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 2977.42863 3007.20543 4.1808395 2000-2009 –> 2010-2019 7.441673 5285.4413 5343.8414
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + phosphate 5612.64989 5641.44423 5.4900639 2000-2009 –> 2010-2019 11.114312 10060.5463 10116.6898
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + phosphate_star 5607.58871 5636.38305 5.4745973 2000-2009 –> 2010-2019 11.109066 10058.0450 10114.1885
mod Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 2928.82612 2958.60292 3.9899507 2000-2009 –> 2010-2019 7.375531 5269.9512 5328.3513
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 5614.55298 5648.14637 5.4897673 2000-2009 –> 2010-2019 11.112600 10064.0945 10129.5953
mod Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2933.72561 2963.50241 4.0087919 2000-2009 –> 2010-2019 7.360938 5266.0973 5324.4974
obs Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5609.54210 5643.13549 5.4744551 2000-2009 –> 2010-2019 11.103829 10060.7230 10126.2238
mod Atlantic (26.5,26.75] cstar_tref ~ temp + aou + silicate + phosphate 2971.74781 2997.27078 4.1660686 2000-2009 –> 2010-2019 7.569586 5315.5335 5365.5907
mod Atlantic (26.5,26.75] cstar_tref ~ temp + aou + silicate + phosphate_star 2978.32167 3003.84464 4.1924859 2000-2009 –> 2010-2019 7.582531 5318.6093 5368.6665
obs Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 5604.24322 5633.03756 5.4643977 2000-2009 –> 2010-2019 11.100364 10055.0744 10111.2179
mod Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 2938.58362 2968.36042 4.0275615 2000-2009 –> 2010-2019 7.374759 5269.6523 5328.0524
obs Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 5606.19549 5639.78888 5.4642523 2000-2009 –> 2010-2019 11.094561 10057.6104 10123.1112
mod Atlantic (26.75,27] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 4779.27629 4813.85106 2.4346898 2000-2009 –> 2010-2019 4.637183 8326.9223 8394.2894
mod Atlantic (26.75,27] cstar_tref ~ aou + silicate + phosphate + phosphate_star 5182.60518 5212.24071 2.9629922 2000-2009 –> 2010-2019 5.207903 8758.7724 8816.5156
mod Atlantic (26.75,27] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 4934.82351 4969.39829 2.6252638 2000-2009 –> 2010-2019 5.078101 8654.7178 8722.0848
mod Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate 5080.76411 5115.33889 2.8176101 2000-2009 –> 2010-2019 4.978419 8597.8377 8665.2048
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate 6872.73619 6907.97074 4.9788333 2000-2009 –> 2010-2019 10.875739 12353.7792 12422.2796
mod Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5106.87088 5141.44566 2.8534753 2000-2009 –> 2010-2019 5.082230 8673.4810 8740.8480
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 6864.40531 6899.63986 4.9605785 2000-2009 –> 2010-2019 10.791479 12326.1776 12394.6780
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 6865.04852 6900.28307 4.9619855 2000-2009 –> 2010-2019 10.784315 12324.3023 12392.8028
mod Atlantic (26.75,27] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5116.84600 5151.42078 2.8672993 2000-2009 –> 2010-2019 5.129458 8707.2588 8774.6258
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 6864.00626 6899.24081 4.9597058 2000-2009 –> 2010-2019 10.787440 12324.8487 12393.3492
obs Atlantic (26.75,27] cstar_tref ~ sal + temp + silicate + phosphate_star 6880.93910 6911.14014 5.0012817 2000-2009 –> 2010-2019 10.884871 12356.1118 12414.8265
mod Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 4787.35268 4821.92745 2.4442354 2000-2009 –> 2010-2019 4.580554 8286.1889 8353.5560
mod Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 4796.68824 4831.26302 2.4553158 2000-2009 –> 2010-2019 4.605142 8305.6093 8372.9763
obs Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 6865.77586 6901.01041 4.9635771 2000-2009 –> 2010-2019 10.857783 12346.0349 12414.5353
mod Atlantic (26.75,27] cstar_tref ~ temp + aou + silicate + phosphate 5234.07394 5263.70946 3.0378075 2000-2009 –> 2010-2019 5.297121 8820.4733 8878.2165
obs Atlantic (26.75,27] cstar_tref ~ temp + aou + silicate + phosphate_star 6864.50334 6894.70438 4.9651694 2000-2009 –> 2010-2019 10.930932 12363.4213 12422.1359
mod Atlantic (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 4820.01557 4854.59035 2.4832232 2000-2009 –> 2010-2019 4.665328 8352.7812 8420.1483
obs Atlantic (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 6865.66986 6900.90441 4.9633451 2000-2009 –> 2010-2019 10.871685 12350.0295 12418.5299
obs Atlantic (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate_star 6863.85669 6894.05773 4.9637540 2000-2009 –> 2010-2019 10.906955 12356.2877 12415.0024
obs Atlantic (26.75,27] cstar_tref ~ temp + silicate + phosphate + phosphate_star 6864.20180 6894.40284 4.9645093 2000-2009 –> 2010-2019 10.927709 12362.3840 12421.0987
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 5917.48938 5954.08302 2.0640435 2000-2009 –> 2010-2019 4.243082 10872.1935 10943.9596
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 6082.63377 6119.22741 2.1916007 2000-2009 –> 2010-2019 4.571860 11235.8944 11307.6606
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate 6212.77164 6244.13761 2.2993179 2000-2009 –> 2010-2019 4.608250 11295.6382 11357.1521
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 6214.01262 6250.60626 2.2986843 2000-2009 –> 2010-2019 4.587638 11279.3433 11351.1095
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 11940.94980 11979.68191 5.8813819 2000-2009 –> 2010-2019 11.645745 22297.5551 22374.0659
mod Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + silicate 6020.96232 6057.55596 2.1430687 2000-2009 –> 2010-2019 4.357511 11011.8977 11083.6639
mod Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate 6213.20159 6249.79523 2.2980075 2000-2009 –> 2010-2019 4.571698 11263.4925 11335.2587
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate 11964.30624 12003.03835 5.9182460 2000-2009 –> 2010-2019 11.651486 22303.2520 22379.7627
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 11951.56690 11990.29901 5.8981106 2000-2009 –> 2010-2019 11.555488 22247.0614 22323.5722
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 11964.36085 12003.09296 5.9183325 2000-2009 –> 2010-2019 11.612274 22280.8700 22357.3807
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 11951.81887 11990.55098 5.8985082 2000-2009 –> 2010-2019 11.562152 22250.9242 22327.4350
mod Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate 6185.34883 6216.71480 2.2765362 2000-2009 –> 2010-2019 4.556485 11239.8191 11301.3330
mod Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 6186.96492 6223.55856 2.2762189 2000-2009 –> 2010-2019 4.532896 11220.3711 11292.1373
obs Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 11903.93052 11942.66264 5.8234231 2000-2009 –> 2010-2019 11.648579 22294.7577 22371.2685
mod Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 6186.94260 6223.53624 2.2762004 2000-2009 –> 2010-2019 4.512111 11199.5664 11271.3326
obs Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 11877.79262 11916.52473 5.7828450 2000-2009 –> 2010-2019 11.510702 22213.6739 22290.1847
obs Atlantic (27,27.25] cstar_tref ~ temp + aou + silicate + phosphate_star 11953.01306 11986.21201 5.9035508 2000-2009 –> 2010-2019 11.658048 22302.0310 22367.6116
mod Atlantic (27,27.25] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 6095.90116 6132.49480 2.2021842 2000-2009 –> 2010-2019 4.345215 11013.1480 11084.9142
obs Atlantic (27,27.25] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 11883.48846 11922.22057 5.7916634 2000-2009 –> 2010-2019 11.536241 22228.8783 22305.3890
obs Atlantic (27,27.25] cstar_tref ~ temp + silicate + phosphate + phosphate_star 11953.20485 11986.40380 5.9038537 2000-2009 –> 2010-2019 11.666405 22306.7847 22372.3653
mod Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate 12706.56633 12741.15211 3.5833513 2000-2009 –> 2010-2019 6.895846 24320.7484 24389.5577
mod Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 12489.79018 12530.14025 3.4207136 2000-2009 –> 2010-2019 6.702198 24064.2671 24144.5445
mod Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 12513.20634 12553.55642 3.4377623 2000-2009 –> 2010-2019 6.743681 24120.5771 24200.8545
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 14723.81243 14764.01531 5.8739448 2000-2009 –> 2010-2019 10.383847 26157.8483 26237.0874
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 12767.50886 12807.85893 3.6284763 2000-2009 –> 2010-2019 6.923502 24360.2452 24440.5226
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 14734.05498 14774.25787 5.8870044 2000-2009 –> 2010-2019 10.373733 26147.9785 26227.2177
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + silicate 12653.47336 12693.82343 3.5416809 2000-2009 –> 2010-2019 6.751302 24129.7682 24210.0456
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + phosphate 14766.16521 14800.62483 5.9307060 2000-2009 –> 2010-2019 10.451580 26207.6866 26275.6059
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + phosphate_star 14772.54841 14807.00803 5.9389200 2000-2009 –> 2010-2019 10.443323 26199.8208 26267.7401
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 14755.52954 14795.73243 5.9144796 2000-2009 –> 2010-2019 10.430856 26195.1654 26274.4045
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 12390.39414 12430.74421 3.3492820 2000-2009 –> 2010-2019 7.049668 24497.5784 24577.8558
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 14763.78229 14803.98518 5.9250725 2000-2009 –> 2010-2019 10.429135 26192.7591 26271.9982
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 12766.70374 12807.05381 3.6278561 2000-2009 –> 2010-2019 6.802376 24194.2412 24274.5186
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 14724.48023 14764.68312 5.8747954 2000-2009 –> 2010-2019 10.366512 26142.7418 26221.9810
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 12373.31871 12413.66878 3.3371617 2000-2009 –> 2010-2019 6.963442 24390.8036 24471.0810
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + phosphate + phosphate_star 14767.50684 14801.96646 5.9324315 2000-2009 –> 2010-2019 10.438891 26196.5607 26264.4800
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 12393.60818 12433.95825 3.3515683 2000-2009 –> 2010-2019 7.051785 24500.5884 24580.8658
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 14758.37501 14798.57790 5.9181298 2000-2009 –> 2010-2019 10.423475 26188.4637 26267.7029
mod Atlantic (27.25,27.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 12730.53554 12770.88561 3.6001044 2000-2009 –> 2010-2019 6.696329 24047.3426 24127.6201
obs Atlantic (27.25,27.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 14742.93632 14783.13921 5.8983520 2000-2009 –> 2010-2019 10.511039 26264.9486 26344.1877
obs Atlantic (27.5,27.75] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 16012.35526 16053.41427 5.2100254 2000-2009 –> 2010-2019 10.158101 29884.7235 29965.9550
obs Atlantic (27.5,27.75] cstar_tref ~ aou + silicate + phosphate + phosphate_star 16036.86375 16072.05718 5.2365914 2000-2009 –> 2010-2019 10.184758 29907.3164 29976.9434
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + nitrate + phosphate 20350.95851 20389.46142 2.2909280 2000-2009 –> 2010-2019 4.771004 39534.6597 39611.1028
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 19863.29181 19908.21187 2.1702409 2000-2009 –> 2010-2019 4.617648 38939.7571 39028.9407
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate 20455.07998 20493.58289 2.3174435 2000-2009 –> 2010-2019 4.772596 39555.5766 39632.0197
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 20381.40759 20426.32765 2.2981425 2000-2009 –> 2010-2019 4.750106 39473.1955 39562.3791
obs Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 16038.32488 16079.38388 5.2360500 2000-2009 –> 2010-2019 10.131980 29862.0434 29943.2749
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 20182.72773 20227.64780 2.2482290 2000-2009 –> 2010-2019 4.723342 39351.9282 39441.1118
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 20440.56071 20485.48077 2.3132163 2000-2009 –> 2010-2019 4.759705 39513.9330 39603.1166
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 16039.33259 16080.39159 5.2370624 2000-2009 –> 2010-2019 10.118367 29849.3122 29930.5437
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 20181.63957 20226.55963 2.2479586 2000-2009 –> 2010-2019 4.674025 39185.9565 39275.1401
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 16061.78642 16102.84542 5.2596729 2000-2009 –> 2010-2019 10.154873 29884.8195 29966.0510
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 20408.41390 20453.33397 2.3050123 2000-2009 –> 2010-2019 4.773158 39554.3934 39643.5771
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 20285.48666 20330.40672 2.2739078 2000-2009 –> 2010-2019 4.701707 39295.6845 39384.8682
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 20191.22970 20236.14976 2.2503425 2000-2009 –> 2010-2019 4.676079 39194.4279 39283.6115
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 16044.64614 16085.70515 5.2424043 2000-2009 –> 2010-2019 10.126632 29857.3750 29938.6065
obs Atlantic (27.5,27.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 16015.19939 16056.25840 5.2128692 2000-2009 –> 2010-2019 10.158440 29885.2424 29966.4739
obs Atlantic (27.5,27.75] cstar_tref ~ temp + aou + silicate + phosphate 16037.97571 16073.16914 5.2377087 2000-2009 –> 2010-2019 10.183434 29906.1616 29975.7886
obs Atlantic (27.5,27.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 16022.23306 16063.29207 5.2199088 2000-2009 –> 2010-2019 10.171989 29898.3160 29979.5474
obs Atlantic (27.5,27.75] cstar_tref ~ temp + silicate + phosphate + phosphate_star 16043.35349 16078.54692 5.2431158 2000-2009 –> 2010-2019 10.195459 29917.6803 29987.3073
mod Atlantic (27.75,27.85] cstar_tref ~ sal + aou + nitrate + silicate 5627.69238 5660.15083 1.3240911 2000-2009 –> 2010-2019 2.556874 10843.4610 10908.1785
obs Atlantic (27.75,27.85] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 10071.67129 10109.28051 5.6966474 2000-2009 –> 2010-2019 10.524423 19059.6578 19134.4549
mod Atlantic (27.75,27.85] cstar_tref ~ sal + aou + silicate + phosphate 5368.04650 5400.50496 1.2240205 2000-2009 –> 2010-2019 2.417942 10481.4441 10546.1616
obs Atlantic (27.75,27.85] cstar_tref ~ sal + aou + silicate + phosphate 10070.85008 10103.08656 5.6987568 2000-2009 –> 2010-2019 10.547027 19069.5364 19133.6482
obs Atlantic (27.75,27.85] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 10063.24884 10100.85807 5.6815983 2000-2009 –> 2010-2019 10.524826 19060.8154 19135.6125
mod Atlantic (27.75,27.85] cstar_tref ~ sal + nitrate + silicate + phosphate_star 6043.61702 6076.07547 1.5017204 2000-2009 –> 2010-2019 2.913760 11693.2776 11757.9951
mod Atlantic (27.75,27.85] cstar_tref ~ sal + silicate + phosphate + phosphate_star 6032.95563 6065.41408 1.4968825 2000-2009 –> 2010-2019 2.913014 11691.8646 11756.5821
obs Atlantic (27.75,27.85] cstar_tref ~ sal + silicate + phosphate + phosphate_star 10062.17382 10094.41030 5.6832491 2000-2009 –> 2010-2019 10.551807 19073.3791 19137.4909
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + silicate 5585.27302 5623.14121 1.3064089 2000-2009 –> 2010-2019 2.539187 10803.0276 10878.5314
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 5248.42249 5286.29068 1.1797821 2000-2009 –> 2010-2019 2.360039 10327.0304 10402.5341
obs Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 10066.09046 10103.69969 5.6866712 2000-2009 –> 2010-2019 10.534409 19066.4477 19141.2448
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5389.84710 5427.71529 1.2313780 2000-2009 –> 2010-2019 2.453053 10578.6858 10654.1896
obs Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 10054.33727 10091.94649 5.6657185 2000-2009 –> 2010-2019 10.527500 19063.3672 19138.1643
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 5396.59613 5434.46432 1.2338958 2000-2009 –> 2010-2019 2.454554 10582.7745 10658.2783
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5382.89703 5420.76522 1.2287904 2000-2009 –> 2010-2019 2.447954 10565.1607 10640.6645
obs Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 10063.06383 10100.67306 5.6812682 2000-2009 –> 2010-2019 10.533376 19066.1225 19140.9196
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate_star 5490.51826 5522.97671 1.2702435 2000-2009 –> 2010-2019 2.543045 10808.3850 10873.1025
obs Atlantic (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 10086.01226 10123.62148 5.7223634 2000-2009 –> 2010-2019 10.560689 19080.5431 19155.3402
obs Atlantic (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 10074.85101 10112.46024 5.7023392 2000-2009 –> 2010-2019 10.549235 19074.6876 19149.4847
obs Atlantic (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 10083.07909 10120.68831 5.7170943 2000-2009 –> 2010-2019 10.557185 19078.7034 19153.5005
mod Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate 6524.77440 6557.07833 1.8288040 2000-2009 –> 2010-2019 4.893552 14808.5228 14873.2011
obs Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 21105.98312 21148.97349 5.2179947 2000-2009 –> 2010-2019 10.575047 41364.2789 41449.9128
mod Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 6521.97070 6559.65863 1.8260778 2000-2009 –> 2010-2019 4.887852 14804.5570 14880.0149
obs Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 20515.91373 20558.90409 4.7884052 2000-2009 –> 2010-2019 9.864777 40422.4611 40508.0950
mod Atlantic (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate 6281.28696 6313.59090 1.6956141 2000-2009 –> 2010-2019 4.614472 14406.1333 14470.8115
obs Atlantic (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate_star 20525.95915 20562.80804 4.7968107 2000-2009 –> 2010-2019 9.873678 40431.1440 40504.5445
obs Atlantic (27.85,27.95] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 20888.11414 20931.10451 5.0550657 2000-2009 –> 2010-2019 10.296218 41003.4513 41089.0852
obs Atlantic (27.85,27.95] cstar_tref ~ sal + silicate + phosphate + phosphate_star 20936.43278 20973.28166 5.0922378 2000-2009 –> 2010-2019 10.341654 41060.0669 41133.4674
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + silicate 6328.66014 6366.34807 1.7196762 2000-2009 –> 2010-2019 4.764014 14592.6393 14668.0973
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + phosphate 6915.65154 6947.95548 2.0648398 2000-2009 –> 2010-2019 4.991456 15049.1465 15113.8247
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 5908.11183 5945.79975 1.5091265 2000-2009 –> 2010-2019 4.375039 13975.3191 14050.7771
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 20418.01591 20461.00628 4.7206345 2000-2009 –> 2010-2019 9.766219 40284.8025 40370.4364
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 6031.37344 6069.06136 1.5680158 2000-2009 –> 2010-2019 4.471025 14140.4814 14215.9393
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 20381.86901 20424.85938 4.6958547 2000-2009 –> 2010-2019 9.731547 40235.8288 40321.4627
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 20439.17495 20482.16532 4.7352004 2000-2009 –> 2010-2019 9.846822 40390.9504 40476.5843
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 6041.10684 6078.79477 1.5727627 2000-2009 –> 2010-2019 4.497235 14174.2145 14249.6725
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate 20457.10382 20493.95271 4.7489604 2000-2009 –> 2010-2019 9.865317 40412.9302 40486.3307
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 6018.85922 6056.54715 1.5619337 2000-2009 –> 2010-2019 4.458285 14120.4863 14195.9443
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 20413.01305 20456.00342 4.7171971 2000-2009 –> 2010-2019 9.760316 40276.6054 40362.2393
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate_star 6154.12835 6186.43228 1.6299588 2000-2009 –> 2010-2019 4.631261 14369.7230 14434.4013
obs Atlantic (27.95,28.05] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 23079.07636 23122.49374 5.7011104 2000-2009 –> 2010-2019 12.011295 46253.4333 46340.0756
obs Atlantic (27.95,28.05] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 22653.59293 22697.01030 5.3783174 2000-2009 –> 2010-2019 11.529932 45647.2014 45733.8436
obs Atlantic (27.95,28.05] cstar_tref ~ sal + aou + silicate + phosphate_star 22733.44888 22770.66377 5.4389646 2000-2009 –> 2010-2019 11.608712 45745.9612 45820.2260
obs Atlantic (27.95,28.05] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 22790.12443 22833.54180 5.4798543 2000-2009 –> 2010-2019 11.667869 45825.6331 45912.2754
obs Atlantic (27.95,28.05] cstar_tref ~ sal + silicate + phosphate + phosphate_star 22897.00664 22934.22153 5.5622009 2000-2009 –> 2010-2019 11.766695 45949.4035 46023.6683
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate 10204.23590 10237.59950 3.4346254 2000-2009 –> 2010-2019 7.187067 20255.0419 20321.4812
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate_star 10238.03284 10271.39645 3.4649721 2000-2009 –> 2010-2019 7.227611 20298.7768 20365.2161
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 10050.65045 10089.57466 3.2983156 2000-2009 –> 2010-2019 7.018468 20071.8091 20149.3216
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 22800.21797 22843.63535 5.4874364 2000-2009 –> 2010-2019 11.684945 45846.6144 45933.2567
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 10064.39307 10103.31728 3.3101346 2000-2009 –> 2010-2019 7.031867 20087.1068 20164.6194
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 22721.43420 22764.85158 5.4285328 2000-2009 –> 2010-2019 11.596435 45733.8217 45820.4639
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate + phosphate_star 10312.24886 10345.61247 3.5325557 2000-2009 –> 2010-2019 7.348338 20424.3534 20490.7927
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 22870.51744 22913.93482 5.5405360 2000-2009 –> 2010-2019 11.952255 46158.2400 46244.8822
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 10044.58724 10083.51145 3.2931145 2000-2009 –> 2010-2019 7.006972 20059.5444 20137.0569
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 22974.30257 23017.71995 5.6198692 2000-2009 –> 2010-2019 12.006185 46233.8307 46320.4729
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + phosphate + phosphate_star 10232.48977 10265.85338 3.4599766 2000-2009 –> 2010-2019 7.221099 20291.7576 20358.1969
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + phosphate_star 10332.77267 10360.57568 3.5533262 2000-2009 –> 2010-2019 7.379001 20452.3584 20507.7245
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 10064.43823 10103.36244 3.3101736 2000-2009 –> 2010-2019 7.031815 20087.0622 20164.5747
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 22774.17995 22817.59732 5.4678984 2000-2009 –> 2010-2019 11.651942 45805.1295 45891.7717
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate_star 10062.44302 10095.80662 3.3101777 2000-2009 –> 2010-2019 7.034663 20085.8647 20152.3040
obs Atlantic (28.05,28.1] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 18909.54780 18951.62496 5.5605174 2000-2009 –> 2010-2019 12.223292 36575.1341 36658.4191
obs Atlantic (28.05,28.1] cstar_tref ~ sal + aou + silicate + phosphate_star 18973.85815 19009.92429 5.6220229 2000-2009 –> 2010-2019 12.292218 36643.3707 36714.7578
obs Atlantic (28.05,28.1] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 18925.52622 18967.60338 5.5752762 2000-2009 –> 2010-2019 12.244042 36595.8976 36679.1826
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 18962.46362 19004.54078 5.6095444 2000-2009 –> 2010-2019 12.249112 36609.4741 36692.7591
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate 3674.65579 3704.68508 1.2749201 2000-2009 –> 2010-2019 2.318935 6634.0735 6693.5971
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate_star 3719.46084 3749.49013 1.3011032 2000-2009 –> 2010-2019 2.368132 6722.8349 6782.3585
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate_star 18991.59343 19027.65957 5.6385881 2000-2009 –> 2010-2019 12.279542 36637.7153 36709.1025
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 3671.34739 3706.38157 1.2718531 2000-2009 –> 2010-2019 2.315665 6632.3723 6701.8166
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 3711.60230 3746.63647 1.2952963 2000-2009 –> 2010-2019 2.362292 6716.9146 6786.3589
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 18938.04099 18980.11815 5.5868631 2000-2009 –> 2010-2019 12.208871 36570.9520 36654.2370
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 18943.38523 18985.46240 5.5918184 2000-2009 –> 2010-2019 12.236173 36594.2322 36677.5172
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + phosphate_star 3772.54847 3802.57777 1.3328232 2000-2009 –> 2010-2019 2.431368 6834.6050 6894.1286
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + phosphate_star 18966.17329 19002.23943 5.6148602 2000-2009 –> 2010-2019 12.259576 36615.3103 36686.6975
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 3730.22324 3765.25741 1.3062862 2000-2009 –> 2010-2019 2.395543 6777.1628 6846.6070
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 18898.11587 18940.19304 5.5499820 2000-2009 –> 2010-2019 12.188588 36544.3545 36627.6395
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + phosphate + phosphate_star 3706.75408 3736.78337 1.2936235 2000-2009 –> 2010-2019 2.354958 6699.3411 6758.8648
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 3701.78270 3736.81688 1.2895381 2000-2009 –> 2010-2019 2.350863 6696.3507 6765.7950
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 18986.57501 19028.65218 5.6320269 2000-2009 –> 2010-2019 12.277618 36638.4129 36721.6979
mod Atlantic (28.05,28.1] cstar_tref ~ temp + aou + phosphate 3715.13491 3740.15932 1.2997308 2000-2009 –> 2010-2019 2.412821 6801.7085 6851.3115
mod Atlantic (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate 3717.05011 3747.07940 1.2996808 2000-2009 –> 2010-2019 2.409764 6800.1717 6859.6954
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 15187.90736 15228.27228 6.0247281 2000-2009 –> 2010-2019 13.443224 28996.2541 29075.8777
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + phosphate 4861.91033 4894.41941 1.0372954 2000-2009 –> 2010-2019 2.365563 10479.5737 10544.5194
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + phosphate_star 4914.61049 4947.11958 1.0538321 2000-2009 –> 2010-2019 2.390224 10552.3490 10617.2947
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + silicate + phosphate 4826.39235 4864.31961 1.0256811 2000-2009 –> 2010-2019 2.271056 10233.9233 10309.6933
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4889.82177 4927.74903 1.0453935 2000-2009 –> 2010-2019 2.304333 10333.0152 10408.7852
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 15195.87608 15236.24100 6.0349082 2000-2009 –> 2010-2019 13.445576 28999.9683 29079.5919
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 15168.66446 15209.02938 6.0002160 2000-2009 –> 2010-2019 13.414628 28974.7923 29054.4158
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + phosphate_star 4931.70690 4964.21598 1.0592532 2000-2009 –> 2010-2019 2.394022 10565.4431 10630.3889
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 4933.25032 4971.17759 1.0591081 2000-2009 –> 2010-2019 2.379536 10533.4275 10609.1974
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 15186.70585 15227.07076 6.0231946 2000-2009 –> 2010-2019 13.430223 28988.8183 29068.4418
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + phosphate + phosphate_star 4899.18577 4931.69485 1.0489649 2000-2009 –> 2010-2019 2.380412 10524.7211 10589.6668
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + phosphate + phosphate_star 15195.06087 15229.65937 6.0364232 2000-2009 –> 2010-2019 13.455108 29001.5103 29069.7591
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4868.23827 4906.16554 1.0386437 2000-2009 –> 2010-2019 2.288381 10287.2804 10363.0504
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 15194.49960 15234.86451 6.0331485 2000-2009 –> 2010-2019 13.445436 28999.4726 29079.0961
obs Atlantic (28.1,28.15] cstar_tref ~ sal + temp + silicate + phosphate_star 15193.90209 15228.50059 6.0349414 2000-2009 –> 2010-2019 13.447868 28997.2222 29065.4709
obs Atlantic (28.1,28.15] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 15188.41694 15228.78186 6.0253786 2000-2009 –> 2010-2019 13.453928 29002.2216 29081.8451
mod Atlantic (28.1,28.15] cstar_tref ~ temp + aou + silicate + phosphate 4983.46035 5015.96944 1.0758342 2000-2009 –> 2010-2019 2.384904 10553.1966 10618.1423
obs Atlantic (28.1,28.15] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 15181.00676 15215.60526 6.0184761 2000-2009 –> 2010-2019 13.458418 28998.9872 29067.2360
obs Atlantic (28.1,28.15] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 15175.08214 15215.44706 6.0083799 2000-2009 –> 2010-2019 13.436483 28988.6448 29068.2684
mod Atlantic (28.1,28.15] cstar_tref ~ temp + silicate + phosphate + phosphate_star 5038.43157 5070.94066 1.0937305 2000-2009 –> 2010-2019 2.409528 10625.0431 10689.9888
mod Atlantic (28.15,28.2] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 7558.29195 7600.47057 0.8307679 2000-2009 –> 2010-2019 1.838790 15574.5067 15658.2569
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate 7592.72529 7628.87839 0.8357316 2000-2009 –> 2010-2019 1.831756 15539.7908 15611.5767
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate + phosphate 4786.65060 4820.13388 3.6095446 2000-2009 –> 2010-2019 7.236273 9754.9535 9822.1693
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate + silicate 7477.54925 7519.72787 0.8198723 2000-2009 –> 2010-2019 1.810803 15397.8638 15481.6140
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate 4714.83478 4748.31806 3.4657041 2000-2009 –> 2010-2019 7.084558 9679.1605 9746.3762
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4706.93017 4740.41344 3.4502263 2000-2009 –> 2010-2019 7.105705 9689.6833 9756.8990
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + nitrate + phosphate_star 7963.18549 7999.33859 0.8879183 2000-2009 –> 2010-2019 1.911728 16064.5521 16136.3380
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 7912.95293 7955.13155 0.8803675 2000-2009 –> 2010-2019 1.897277 15978.3974 16062.1477
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4707.76988 4741.25315 3.4518672 2000-2009 –> 2010-2019 7.097083 9685.3781 9752.5938
mod Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + phosphate 7593.43483 7629.58792 0.8358285 2000-2009 –> 2010-2019 1.831034 15535.8906 15607.6765
mod Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + phosphate_star 7630.10581 7666.25891 0.8408551 2000-2009 –> 2010-2019 1.840870 15599.5957 15671.3816
mod Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate 7512.85191 7549.00501 0.8248881 2000-2009 –> 2010-2019 1.823918 15476.8135 15548.5995
obs Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 4730.18098 4763.66426 3.4959516 2000-2009 –> 2010-2019 7.119746 9697.0032 9764.2190
obs Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 4731.64572 4765.12900 3.4988524 2000-2009 –> 2010-2019 7.164103 9719.2842 9786.4999
obs Atlantic (28.15,28.2] cstar_tref ~ temp + aou + silicate + phosphate 4732.27134 4760.97129 3.5040582 2000-2009 –> 2010-2019 7.128043 9697.1898 9754.8033
obs Atlantic (28.15,28.2] cstar_tref ~ temp + aou + silicate + phosphate_star 4734.43878 4763.13873 3.5083615 2000-2009 –> 2010-2019 7.173706 9720.1243 9777.7378
mod Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 7687.01248 7723.16558 0.8487154 2000-2009 –> 2010-2019 1.854389 15688.1490 15759.9349
obs Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 4728.99325 4762.47653 3.4936012 2000-2009 –> 2010-2019 7.147143 9710.7764 9777.9922
mod Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + silicate + phosphate_star 7935.07937 7971.23246 0.8838473 2000-2009 –> 2010-2019 1.901283 16001.4217 16073.2076
obs Atlantic (28.15,28.2] cstar_tref ~ temp + silicate + phosphate + phosphate_star 4731.11881 4759.81877 3.5017722 2000-2009 –> 2010-2019 7.155879 9711.1849 9768.7984
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate 10252.48910 10285.37113 0.6353516 2000-2009 –> 2010-2019 1.313525 19434.7046 19499.5899
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + phosphate 10153.24015 10192.69858 0.6293174 2000-2009 –> 2010-2019 1.305186 19307.1586 19385.0210
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 10981.43096 11020.96212 3.3154281 2000-2009 –> 2010-2019 6.527253 22036.7016 22115.8963
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + phosphate_star 10065.35519 10104.81362 0.6241262 2000-2009 –> 2010-2019 1.300159 19221.4420 19299.3043
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate 10098.86519 10138.32362 0.6261005 2000-2009 –> 2010-2019 1.287030 19053.8620 19131.7244
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 9279.40154 9325.43637 0.5794546 2000-2009 –> 2010-2019 1.226098 18041.9035 18132.7430
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + phosphate + phosphate_star 10979.43110 11013.31495 3.3154282 2000-2009 –> 2010-2019 6.535465 22043.6055 22111.4867
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 10941.92201 10981.45317 3.2843127 2000-2009 –> 2010-2019 6.504145 22007.8262 22087.0209
obs Atlantic (28.2, Inf] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 11008.63629 11048.16745 3.3370249 2000-2009 –> 2010-2019 6.549616 22064.9256 22144.1203
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate 10171.05568 10210.51411 0.6303750 2000-2009 –> 2010-2019 1.307142 19336.8007 19414.6631
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate 10072.55493 10118.58977 0.6244321 2000-2009 –> 2010-2019 1.300145 19226.4199 19317.2594
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 10066.55334 10112.58817 0.6240790 2000-2009 –> 2010-2019 1.300040 19223.6897 19314.5291
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 11028.86031 11068.39147 3.3531708 2000-2009 –> 2010-2019 6.548899 22062.6776 22141.8723
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + silicate 10017.08964 10063.12447 0.6211763 2000-2009 –> 2010-2019 1.279243 18935.4453 19026.2847
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + silicate + phosphate 10944.49486 10984.02602 3.2863301 2000-2009 –> 2010-2019 6.538331 22052.8475 22132.0422
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 10941.75482 10981.28598 3.2841817 2000-2009 –> 2010-2019 6.481648 21977.8936 22057.0883
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 10114.85865 10160.89349 0.6269268 2000-2009 –> 2010-2019 1.306182 19315.2553 19406.0947
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 11000.15093 11039.68209 3.3302738 2000-2009 –> 2010-2019 6.536045 22047.3663 22126.5610
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + phosphate + phosphate_star 11000.32749 11034.21134 3.3320042 2000-2009 –> 2010-2019 6.544658 22054.7003 22122.5815
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 10927.58482 10967.11598 3.2730938 2000-2009 –> 2010-2019 6.485700 21983.8943 22063.0891
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 22636.52619 22680.33515 4.5332563 2000-2009 –> 2010-2019 9.981143 39927.9210 40013.2263
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 30092.73735 30138.27036 5.0870655 2000-2009 –> 2010-2019 11.275555 53427.0200 53515.8680
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + phosphate + phosphate_star 24276.58065 24314.13119 5.6077068 2000-2009 –> 2010-2019 11.425391 41930.3377 42003.4565
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 23836.08781 23879.89677 5.2953238 2000-2009 –> 2010-2019 10.885285 41270.3132 41355.6185
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 31247.94025 31293.47326 5.7183017 2000-2009 –> 2010-2019 12.179809 54892.7997 54981.6477
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate 22704.74267 22742.29321 4.5746761 2000-2009 –> 2010-2019 11.016988 40924.3132 40997.4320
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate 30435.53114 30474.55943 5.2678032 2000-2009 –> 2010-2019 11.620417 53956.1190 54032.2745
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 22618.42982 22662.23877 4.5226424 2000-2009 –> 2010-2019 10.254099 40191.3406 40276.6458
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 30305.03183 30350.56484 5.1976008 2000-2009 –> 2010-2019 11.411928 53669.2885 53758.1365
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 22631.21979 22675.02875 4.5301414 2000-2009 –> 2010-2019 10.355592 40294.3779 40379.6831
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 30325.88600 30371.41901 5.2085876 2000-2009 –> 2010-2019 11.432371 53701.0811 53789.9291
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 22296.11384 22339.92279 4.3377060 2000-2009 –> 2010-2019 10.367840 40150.8602 40236.1655
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 30431.69870 30477.23171 5.2646931 2000-2009 –> 2010-2019 11.615225 53951.9284 54040.7764
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 24508.88903 24552.69799 5.7775202 2000-2009 –> 2010-2019 11.826424 42380.8782 42466.1835
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 31505.68022 31551.21323 5.8695001 2000-2009 –> 2010-2019 12.433425 55263.6736 55352.5216
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 24655.24818 24699.05713 5.8880978 2000-2009 –> 2010-2019 12.104983 42679.2077 42764.5130
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 31579.85157 31625.38458 5.9137475 2000-2009 –> 2010-2019 12.514090 55377.6478 55466.4958
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 22761.40692 22805.21588 4.6071837 2000-2009 –> 2010-2019 10.526160 40512.9294 40598.2347
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 30465.14396 30510.67697 5.2825523 2000-2009 –> 2010-2019 11.550873 53891.6362 53980.4842
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 31621.00954 31666.54255 5.9384443 2000-2009 –> 2010-2019 12.554245 55435.6347 55524.4827
mod Indo-Pacific (26,26.5] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 16238.32782 16279.99613 4.1973068 2000-2009 –> 2010-2019 8.537006 31863.8375 31946.8232
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 15994.31847 16035.98678 4.0209935 2000-2009 –> 2010-2019 8.191875 31405.2515 31488.2372
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 15360.08537 15401.75367 3.5965897 2000-2009 –> 2010-2019 7.528212 30451.5383 30534.5240
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 16398.40878 16440.07709 4.3171549 2000-2009 –> 2010-2019 8.506118 31832.7356 31915.7213
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 16465.06596 16506.65752 4.5096315 2000-2009 –> 2010-2019 9.220876 32707.0224 32790.0060
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 16486.49131 16528.08287 4.5268443 2000-2009 –> 2010-2019 9.254318 32747.2443 32830.2280
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + phosphate 16846.51984 16882.16974 4.8278302 2000-2009 –> 2010-2019 9.862062 33448.9203 33520.0491
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 15809.01367 15850.68197 3.8920630 2000-2009 –> 2010-2019 7.838223 30920.4275 31003.4133
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 16671.82029 16713.41185 4.6785035 2000-2009 –> 2010-2019 9.515625 33057.9022 33140.8858
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 16018.41191 16060.08022 4.0380679 2000-2009 –> 2010-2019 8.155206 31359.2072 31442.1930
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 16706.22527 16747.81683 4.7072121 2000-2009 –> 2010-2019 9.582038 33134.7485 33217.7321
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 16491.57315 16533.16471 4.5309366 2000-2009 –> 2010-2019 9.263927 32758.7007 32841.6843
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate_star 16605.80231 16641.45222 4.6255501 2000-2009 –> 2010-2019 9.554787 33092.9976 33164.1265
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 16561.87789 16603.46945 4.5879326 2000-2009 –> 2010-2019 9.467035 32995.1947 33078.1783
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 16061.53107 16103.19937 4.0688066 2000-2009 –> 2010-2019 8.219409 31446.1056 31529.0914
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 16704.94373 16746.53529 4.7061396 2000-2009 –> 2010-2019 9.582985 33135.7306 33218.7143
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate_star 16714.97845 16750.62836 4.7162210 2000-2009 –> 2010-2019 9.617605 33171.2002 33242.3291
mod Indo-Pacific (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 16347.82083 16389.48914 4.2789159 2000-2009 –> 2010-2019 8.467573 31781.7523 31864.7380
mod Indo-Pacific (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 16440.11015 16481.77846 4.3489335 2000-2009 –> 2010-2019 8.559790 31902.6275 31985.6132
mod Indo-Pacific (26,26.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 16488.51568 16530.18399 4.3861145 2000-2009 –> 2010-2019 8.618071 31978.0645 32061.0502
mod Indo-Pacific (26.5,26.75] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 15221.76733 15262.80211 4.5222539 2000-2009 –> 2010-2019 8.759629 29137.9247 29219.5233
mod Indo-Pacific (26.5,26.75] cstar_tref ~ aou + silicate + phosphate + phosphate_star 15308.54203 15343.71471 4.6002115 2000-2009 –> 2010-2019 9.088425 29501.9722 29571.9138
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 15341.70243 15382.73721 4.6278926 2000-2009 –> 2010-2019 8.973335 29380.1510 29461.7496
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 16778.78261 16820.20537 5.1282062 2000-2009 –> 2010-2019 10.228058 34283.0668 34366.2216
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 15497.47446 15538.50924 4.7687888 2000-2009 –> 2010-2019 9.214617 29646.8283 29728.4268
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 16821.03562 16862.45837 5.1678270 2000-2009 –> 2010-2019 10.328424 34393.2626 34476.4174
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + phosphate 17161.33252 17196.83774 5.5002935 2000-2009 –> 2010-2019 10.954720 35049.3023 35120.5778
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 14942.49914 14983.53393 4.2855249 2000-2009 –> 2010-2019 8.505445 28838.6126 28920.2111
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 17146.04192 17187.46467 5.4829977 2000-2009 –> 2010-2019 10.922949 35020.7642 35103.9189
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 15095.08345 15136.11824 4.4132884 2000-2009 –> 2010-2019 8.749293 29122.9740 29204.5725
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 17184.11080 17225.53355 5.5211502 2000-2009 –> 2010-2019 10.997911 35097.5287 35180.6835
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 16826.85844 16868.28119 5.1733110 2000-2009 –> 2010-2019 10.343556 34409.8029 34492.9576
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate_star 16873.96111 16909.46633 5.2197888 2000-2009 –> 2010-2019 10.498220 34573.7336 34645.0091
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 15234.78594 15275.82072 4.5336030 2000-2009 –> 2010-2019 9.006957 29414.1122 29495.7107
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 16818.30583 16859.72859 5.1652580 2000-2009 –> 2010-2019 10.364639 34433.4947 34516.6495
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + phosphate + phosphate_star 17197.34030 17232.84552 5.5364873 2000-2009 –> 2010-2019 11.035638 35132.1683 35203.4438
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 15105.50021 15146.53499 4.4221483 2000-2009 –> 2010-2019 8.767020 29143.3117 29224.9103
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 17174.10612 17215.52887 5.5110979 2000-2009 –> 2010-2019 10.983211 35082.6528 35165.8076
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate_star 15235.26627 15270.43894 4.5357684 2000-2009 –> 2010-2019 9.024942 29429.7341 29499.6757
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 15390.80223 15431.83701 4.6718483 2000-2009 –> 2010-2019 9.012346 29423.7218 29505.3204
obs Indo-Pacific (26.75,27] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 24706.64851 24744.56682 4.9021248 2000-2009 –> 2010-2019 10.014670 51265.2813 51341.4685
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 25072.33193 25116.56995 5.1242140 2000-2009 –> 2010-2019 10.453812 51994.7809 52083.6661
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + phosphate 18975.28803 19018.25594 3.8603199 2000-2009 –> 2010-2019 7.512165 38542.1770 38628.4696
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + phosphate 24897.47385 24941.71187 5.0162055 2000-2009 –> 2010-2019 10.239573 51644.7207 51733.6059
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 19182.26324 19225.23115 3.9788111 2000-2009 –> 2010-2019 7.713575 38910.8955 38997.1882
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 25069.40149 25113.63951 5.1223849 2000-2009 –> 2010-2019 10.462910 52009.6748 52098.5600
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 19196.63115 19239.59906 3.9871703 2000-2009 –> 2010-2019 7.726889 38934.8144 39021.1070
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + phosphate_star 19318.13946 19354.96910 4.0597554 2000-2009 –> 2010-2019 7.837853 39127.8783 39201.8434
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 19227.64573 19270.61364 4.0052745 2000-2009 –> 2010-2019 7.774451 39022.3532 39108.6459
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 24937.66301 24981.90103 5.0408268 2000-2009 –> 2010-2019 10.228168 51624.6841 51713.5692
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + phosphate 19650.06532 19686.89495 4.2614411 2000-2009 –> 2010-2019 8.189644 39740.4825 39814.4476
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + phosphate 24945.85107 24983.76937 5.0470875 2000-2009 –> 2010-2019 10.309155 51755.3341 51831.5214
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 19408.25275 19451.22066 4.1123453 2000-2009 –> 2010-2019 7.990205 39407.7476 39494.0402
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 24851.49504 24895.73306 4.9881848 2000-2009 –> 2010-2019 10.145146 51487.4020 51576.2871
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 19581.70762 19624.67553 4.2178697 2000-2009 –> 2010-2019 8.200445 39773.1551 39859.4477
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 24972.82618 25017.06420 5.0624681 2000-2009 –> 2010-2019 10.291477 51729.4678 51818.3529
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 19591.49237 19634.46028 4.2239024 2000-2009 –> 2010-2019 8.213769 39796.1174 39882.4101
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 24978.94190 25023.17992 5.0662415 2000-2009 –> 2010-2019 10.296885 51738.3027 51827.1878
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate_star 19648.81345 19685.64308 4.2606619 2000-2009 –> 2010-2019 8.294070 39929.6293 40003.5944
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate_star 24979.40126 25017.31956 5.0677597 2000-2009 –> 2010-2019 10.298710 51737.2721 51813.4594
mod Indo-Pacific (27,27.25] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 21880.03532 21918.45812 2.8022935 2000-2009 –> 2010-2019 5.890306 45726.9576 45804.0868
mod Indo-Pacific (27,27.25] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 21598.12522 21642.95182 2.7145826 2000-2009 –> 2010-2019 5.655170 44989.1719 45079.1558
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 31385.23616 31431.34972 4.5030522 2000-2009 –> 2010-2019 9.311013 66725.7391 66818.6423
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 31569.81626 31615.92982 4.5811849 2000-2009 –> 2010-2019 9.493390 67163.8129 67256.7161
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate 31385.28390 31424.80981 4.5039116 2000-2009 –> 2010-2019 9.312330 66724.9122 66804.5435
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate_star 31613.92361 31653.44952 4.6009128 2000-2009 –> 2010-2019 9.519742 67221.8480 67301.4793
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 31500.72478 31546.83834 4.5517810 2000-2009 –> 2010-2019 9.427618 67006.9013 67099.8044
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + silicate + phosphate + phosphate_star 31545.57099 31585.09690 4.5716971 2000-2009 –> 2010-2019 9.456238 67070.8244 67150.4557
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate 21841.01534 21885.84195 2.7894478 2000-2009 –> 2010-2019 5.852412 45613.7067 45703.6907
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 22142.10246 22186.92906 2.8851231 2000-2009 –> 2010-2019 6.043872 46203.0145 46292.9985
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate 31385.48589 31431.59945 4.5031570 2000-2009 –> 2010-2019 9.311237 66726.2820 66819.1852
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 31483.04093 31529.15449 4.5442855 2000-2009 –> 2010-2019 9.406681 66956.5931 67049.4963
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 22161.57969 22206.40629 2.8914241 2000-2009 –> 2010-2019 6.058364 46246.7310 46336.7149
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 22238.60491 22283.43151 2.9164774 2000-2009 –> 2010-2019 6.090364 46344.2675 46434.2515
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 31558.50127 31604.61483 4.5763565 2000-2009 –> 2010-2019 9.461186 67086.4535 67179.3567
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 31467.09232 31513.20588 4.5375360 2000-2009 –> 2010-2019 9.389558 66915.4052 67008.3084
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + phosphate 22054.12220 22092.54500 2.8574716 2000-2009 –> 2010-2019 5.990738 46037.2180 46114.3471
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 21922.28061 21967.10721 2.8149541 2000-2009 –> 2010-2019 5.859162 45637.4793 45727.4633
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 22234.23459 22279.06119 2.9150502 2000-2009 –> 2010-2019 6.072640 46291.7129 46381.6969
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 22242.46835 22287.29495 2.9177398 2000-2009 –> 2010-2019 6.079915 46313.5271 46403.5111
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 29335.26711 29381.93967 3.0161225 2000-2009 –> 2010-2019 6.073220 59344.7866 59438.2535
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 32603.91114 32650.73975 3.7561580 2000-2009 –> 2010-2019 8.115891 67782.5023 67876.3213
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 32908.95874 32955.78735 3.8538221 2000-2009 –> 2010-2019 8.394277 68581.5274 68675.3465
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + silicate + phosphate 32677.34307 32717.48187 3.7800755 2000-2009 –> 2010-2019 8.146178 67871.6895 67952.1058
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + silicate + phosphate_star 33044.09666 33084.23547 3.8985516 2000-2009 –> 2010-2019 8.457266 68763.4746 68843.8909
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 29710.13017 29756.80273 3.1149923 2000-2009 –> 2010-2019 6.320847 60281.5369 60375.0038
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 32736.05170 32782.88031 3.7981566 2000-2009 –> 2010-2019 8.248382 68164.4985 68258.3176
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + silicate + phosphate + phosphate_star 32900.17701 32940.31581 3.8516235 2000-2009 –> 2010-2019 8.335736 68418.8833 68499.2996
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 29689.31870 29735.99125 3.1094193 2000-2009 –> 2010-2019 6.282208 60138.1141 60231.5810
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 29905.57419 29952.24675 3.1678192 2000-2009 –> 2010-2019 6.407843 60602.3577 60695.8246
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 32603.20390 32650.03251 3.7559345 2000-2009 –> 2010-2019 8.078762 67678.4056 67772.2247
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 32638.69331 32685.52192 3.7671677 2000-2009 –> 2010-2019 8.136679 67844.5300 67938.3491
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 29339.96057 29386.63312 3.0173408 2000-2009 –> 2010-2019 6.048578 59249.0158 59342.4827
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 32777.45608 32824.28469 3.8114126 2000-2009 –> 2010-2019 8.182385 67987.3612 68081.1803
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 32631.99606 32678.82467 3.7650453 2000-2009 –> 2010-2019 8.125732 67813.2462 67907.0653
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + phosphate 29848.80560 29888.81065 3.1529260 2000-2009 –> 2010-2019 6.383824 60510.2348 60590.3493
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 29440.04166 29486.71421 3.0434364 2000-2009 –> 2010-2019 6.152806 59650.0609 59743.5278
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 29707.75900 29754.43155 3.1143568 2000-2009 –> 2010-2019 6.309413 60239.2638 60332.7307
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 29705.85621 29752.52876 3.1138470 2000-2009 –> 2010-2019 6.310340 60242.6765 60336.1434
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + nitrate + silicate + phosphate_star 29706.83119 29746.83623 3.1146442 2000-2009 –> 2010-2019 6.313008 60248.5725 60328.6870
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 37678.54025 37726.23015 3.9887688 2000-2009 –> 2010-2019 7.779385 75316.1003 75411.6009
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 31764.44334 31811.80919 2.8730640 2000-2009 –> 2010-2019 5.830636 66739.6333 66834.9578
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 36153.81863 36201.50854 3.5609810 2000-2009 –> 2010-2019 7.478289 74240.9284 74336.4291
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 31658.89709 31706.26294 2.8495293 2000-2009 –> 2010-2019 5.810375 66649.5353 66744.8598
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate 36394.47890 36435.35596 3.6258587 2000-2009 –> 2010-2019 7.569199 74570.1590 74652.0167
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate_star 32364.11279 32404.71209 3.0109850 2000-2009 –> 2010-2019 6.102631 67956.4821 68038.1888
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 31672.77784 31720.14370 2.8526134 2000-2009 –> 2010-2019 5.808974 66642.2467 66737.5712
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 36452.93158 36500.62148 3.6411206 2000-2009 –> 2010-2019 7.764753 75241.9245 75337.4251
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 37553.14222 37600.83212 3.9517258 2000-2009 –> 2010-2019 7.772152 75297.8191 75393.3197
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 31947.51098 31994.87683 2.9143463 2000-2009 –> 2010-2019 5.971106 67383.3907 67478.7151
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 36088.56967 36136.25957 3.5437349 2000-2009 –> 2010-2019 7.250042 73418.5678 73514.0684
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 31966.14235 32013.50820 2.9185809 2000-2009 –> 2010-2019 5.989064 67464.5857 67559.9102
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 36030.79668 36078.48658 3.5285346 2000-2009 –> 2010-2019 7.262841 73463.7052 73559.2058
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 31668.88550 31716.25135 2.8517482 2000-2009 –> 2010-2019 5.903170 67080.3563 67175.6808
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 32253.36720 32300.73305 2.9846455 2000-2009 –> 2010-2019 5.987338 67440.0111 67535.3356
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 36186.77385 36234.46375 3.5697233 2000-2009 –> 2010-2019 7.294461 73584.6019 73680.1025
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate 31972.26567 32012.86497 2.9204291 2000-2009 –> 2010-2019 6.007284 67542.9786 67624.6853
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 31938.53718 31985.90303 2.9123090 2000-2009 –> 2010-2019 5.975301 67402.8655 67498.1900
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 36038.82600 36086.51590 3.5306432 2000-2009 –> 2010-2019 7.263047 73464.7685 73560.2692
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate_star 36404.12969 36445.00675 3.6284632 2000-2009 –> 2010-2019 7.383700 73911.4668 73993.3245
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 10423.95529 10464.63627 1.9921080 2000-2009 –> 2010-2019 4.350527 22860.8864 22942.9466
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + phosphate 10365.68404 10406.36502 1.9687381 2000-2009 –> 2010-2019 4.242829 22603.9544 22686.0147
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 10371.46297 10412.14395 1.9710435 2000-2009 –> 2010-2019 4.253878 22630.6712 22712.7314
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + phosphate 10520.09474 10554.96415 2.0320958 2000-2009 –> 2010-2019 4.370044 22907.4611 22977.7984
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + phosphate_star 10524.95093 10559.82035 2.0340953 2000-2009 –> 2010-2019 4.380311 22931.5775 23001.9149
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 10520.89051 10561.57149 2.0316003 2000-2009 –> 2010-2019 4.366959 22904.2099 22986.2701
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 19270.70331 19314.01246 3.5256634 2000-2009 –> 2010-2019 6.788495 37183.2576 37269.5694
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 10525.72284 10566.40382 2.0335894 2000-2009 –> 2010-2019 4.377300 22928.5207 23010.5810
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 19209.94922 19253.25836 3.4959895 2000-2009 –> 2010-2019 6.758183 37121.1565 37207.4682
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 10375.24704 10415.92802 1.9725545 2000-2009 –> 2010-2019 4.261924 22650.0513 22732.1115
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 19192.26994 19235.57909 3.4874015 2000-2009 –> 2010-2019 6.758787 37122.8366 37209.1483
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + phosphate + phosphate_star 10526.09176 10560.96118 2.0345652 2000-2009 –> 2010-2019 4.382990 22937.8517 23008.1891
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 10526.70612 10567.38710 2.0339944 2000-2009 –> 2010-2019 4.380457 22935.9051 23017.9654
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 19208.74656 19252.05571 3.4954046 2000-2009 –> 2010-2019 6.757674 37120.1137 37206.4254
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate_star 19210.92861 19248.05074 3.4974389 2000-2009 –> 2010-2019 6.782703 37168.6223 37242.6038
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 19263.02805 19306.33719 3.5219007 2000-2009 –> 2010-2019 6.784406 37174.8926 37261.2043
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 19215.05566 19258.36481 3.4984740 2000-2009 –> 2010-2019 6.760459 37125.8236 37212.1354
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 19220.11258 19263.42172 3.5009361 2000-2009 –> 2010-2019 6.763051 37131.1554 37217.4671
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate_star 19238.07245 19275.19458 3.5106711 2000-2009 –> 2010-2019 6.783638 37169.9640 37243.9455
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + silicate + phosphate + phosphate_star 19271.28123 19308.40335 3.5269280 2000-2009 –> 2010-2019 6.789934 37182.2030 37256.1845
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 12167.64139 12209.79708 1.7767201 2000-2009 –> 2010-2019 3.818931 25986.2330 26070.9633
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate 11992.39636 12028.52980 1.7269375 2000-2009 –> 2010-2019 3.681391 25524.7237 25597.3497
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + phosphate 11848.53553 11890.69122 1.6861069 2000-2009 –> 2010-2019 3.636190 25368.3736 25453.1039
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 11851.04399 11893.19967 1.6868009 2000-2009 –> 2010-2019 3.638252 25375.4204 25460.1507
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + silicate 11884.45864 11926.61432 1.6960723 2000-2009 –> 2010-2019 3.628111 25344.1308 25428.8612
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + phosphate 12106.58233 12142.71578 1.7595901 2000-2009 –> 2010-2019 3.834650 26026.4483 26099.0743
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 12090.06406 12132.21975 1.7542530 2000-2009 –> 2010-2019 3.826988 26004.6734 26089.4037
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 26672.37463 26718.18011 3.2456616 2000-2009 –> 2010-2019 6.156964 50322.9704 50414.0387
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 12096.87963 12139.03532 1.7562154 2000-2009 –> 2010-2019 3.835163 26030.8585 26115.5889
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 26634.15347 26679.95895 3.2336026 2000-2009 –> 2010-2019 6.153137 50311.5786 50402.6468
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 11855.67468 11897.83036 1.6880827 2000-2009 –> 2010-2019 3.649702 25413.6862 25498.4165
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 26616.26899 26662.07448 3.2279753 2000-2009 –> 2010-2019 6.157442 50325.9675 50417.0358
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 12098.81164 12140.96733 1.7567721 2000-2009 –> 2010-2019 3.838412 26041.1686 26125.8990
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 26633.09192 26678.89740 3.2332683 2000-2009 –> 2010-2019 6.153923 50314.1638 50405.2321
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 26655.12859 26700.93407 3.2402148 2000-2009 –> 2010-2019 6.151581 50305.9337 50397.0019
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 26631.19676 26677.00224 3.2326716 2000-2009 –> 2010-2019 6.152227 50308.6918 50399.7601
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + silicate + phosphate 26692.40833 26731.67018 3.2526338 2000-2009 –> 2010-2019 6.164057 50341.3987 50419.4572
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + silicate + phosphate_star 26665.38248 26704.64432 3.2440840 2000-2009 –> 2010-2019 6.163640 50340.8777 50418.9362
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 26631.97717 26677.78266 3.2329173 2000-2009 –> 2010-2019 6.153558 50313.0035 50404.0717
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + nitrate + silicate + phosphate_star 26634.99432 26674.25617 3.2344973 2000-2009 –> 2010-2019 6.164041 50342.9428 50421.0013
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 44065.94289 44108.29301 3.1434437 2000-2009 –> 2010-2019 6.007280 84613.8760 84698.2989
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 44027.22206 44076.63054 3.1360018 2000-2009 –> 2010-2019 5.989990 84520.6485 84619.1419
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + phosphate + phosphate_star 44107.58211 44142.87388 3.1514385 2000-2009 –> 2010-2019 6.015949 84657.3806 84727.7331
obs Indo-Pacific (27.95,28.05] cstar_tref ~ aou + silicate + phosphate + phosphate_star 44065.92027 44108.27040 3.1434395 2000-2009 –> 2010-2019 6.000771 84576.5552 84660.9781
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 44067.78477 44117.19325 3.1434147 2000-2009 –> 2010-2019 6.007187 84617.3515 84715.8449
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + aou + phosphate + phosphate_star 44109.11690 44151.46703 3.1513532 2000-2009 –> 2010-2019 6.015638 84659.6233 84744.0463
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate 15868.16424 15906.15087 1.6339207 2000-2009 –> 2010-2019 3.337076 29700.1274 29775.1470
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate + silicate 15550.44086 15594.75859 1.5721921 2000-2009 –> 2010-2019 3.233902 29209.9361 29297.4590
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate 16053.21505 16091.20168 1.6707496 2000-2009 –> 2010-2019 3.494708 30370.4826 30445.5022
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate_star 16053.00334 16090.98996 1.6707070 2000-2009 –> 2010-2019 3.499374 30388.5296 30463.5492
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 15778.22693 15822.54466 1.6159263 2000-2009 –> 2010-2019 3.415980 30004.0653 30091.5882
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 43980.89245 44030.30093 3.1275562 2000-2009 –> 2010-2019 5.989642 84520.7965 84619.2899
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 15781.87555 15826.19328 1.6166367 2000-2009 –> 2010-2019 3.422969 30032.3719 30119.8948
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 15827.22686 15871.54459 1.6254920 2000-2009 –> 2010-2019 3.477048 30252.8468 30340.3697
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 44051.97484 44101.38332 3.1405233 2000-2009 –> 2010-2019 6.022761 84706.9747 84805.4682
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + phosphate + phosphate_star 16052.92195 16090.90857 1.6706906 2000-2009 –> 2010-2019 3.501559 30396.9696 30471.9893
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 15782.60160 15826.91933 1.6167780 2000-2009 –> 2010-2019 3.424838 30039.8710 30127.3939
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate_star 15844.02908 15882.01570 1.6291776 2000-2009 –> 2010-2019 3.500146 30341.5127 30416.5324
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 44057.25030 44106.65877 3.1414878 2000-2009 –> 2010-2019 6.007652 84620.5140 84719.0075
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + silicate + phosphate 44103.85760 44146.20773 3.1503886 2000-2009 –> 2010-2019 6.028609 84733.9786 84818.4015
obs Indo-Pacific (28.05,28.1] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 39619.29455 39661.11791 2.9983066 2000-2009 –> 2010-2019 5.678208 69041.0010 69123.1354
obs Indo-Pacific (28.05,28.1] cstar_tref ~ aou + silicate + phosphate + phosphate_star 39724.57490 39766.39826 3.0184337 2000-2009 –> 2010-2019 5.702290 69164.3170 69246.4514
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 39428.53893 39477.33284 2.9618033 2000-2009 –> 2010-2019 5.632574 68810.5079 68906.3313
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + aou + silicate + phosphate 39477.59188 39519.41523 2.9714280 2000-2009 –> 2010-2019 5.649172 68889.4499 68971.5843
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 39516.99573 39565.78964 2.9784994 2000-2009 –> 2010-2019 5.700094 69129.5083 69225.3317
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + silicate 8387.44603 8427.70636 1.4649228 2000-2009 –> 2010-2019 3.121809 16477.8879 16557.6927
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 8326.92162 8367.18195 1.4459789 2000-2009 –> 2010-2019 3.028021 16223.3164 16303.1213
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 39476.60467 39525.39858 2.9708640 2000-2009 –> 2010-2019 5.615354 68737.6308 68833.4542
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 8327.67381 8367.93414 1.4462128 2000-2009 –> 2010-2019 3.033081 16236.8549 16316.6598
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 39563.85347 39612.64738 2.9873818 2000-2009 –> 2010-2019 5.656105 68936.4440 69032.2674
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 8330.71757 8370.97790 1.4471598 2000-2009 –> 2010-2019 3.034217 16240.3996 16320.2044
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 39722.84573 39771.63965 3.0177185 2000-2009 –> 2010-2019 5.693764 69128.9450 69224.7684
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 8327.51235 8367.77268 1.4461626 2000-2009 –> 2010-2019 3.033111 16236.9060 16316.7109
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 39543.73946 39592.53337 2.9835657 2000-2009 –> 2010-2019 5.649952 68905.6130 69001.4364
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate_star 8366.46916 8400.97801 1.4589566 2000-2009 –> 2010-2019 3.046085 16274.3413 16342.7455
obs Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 39776.65283 39825.44674 3.0280548 2000-2009 –> 2010-2019 5.696742 69149.0772 69244.9006
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate 8816.34278 8850.85163 1.6071598 2000-2009 –> 2010-2019 3.211521 16769.5496 16837.9537
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate_star 8823.04069 8857.54954 1.6094764 2000-2009 –> 2010-2019 3.220081 16792.5525 16860.9566
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + nitrate + silicate + phosphate_star 8821.27698 8855.78583 1.6088661 2000-2009 –> 2010-2019 3.219026 16789.6292 16858.0334
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + silicate + phosphate + phosphate_star 8821.43333 8855.94219 1.6089202 2000-2009 –> 2010-2019 3.219604 16791.1524 16859.5566
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 26722.59172 26768.90338 2.7203385 2000-2009 –> 2010-2019 5.371768 53722.1812 53814.9513
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 26576.50574 26622.81741 2.6845725 2000-2009 –> 2010-2019 5.354116 53652.8405 53745.6107
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + phosphate + phosphate_star 26785.48787 26825.18358 2.7363794 2000-2009 –> 2010-2019 5.391982 53800.8049 53880.3222
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + silicate + phosphate 17536.69366 17575.96951 1.3281813 2000-2009 –> 2010-2019 2.901202 36863.3331 36941.8999
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + aou + silicate + phosphate 26593.90763 26633.60334 2.6892954 2000-2009 –> 2010-2019 5.362703 53684.5458 53764.0631
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + nitrate + silicate + phosphate 17547.01126 17586.28711 1.3295134 2000-2009 –> 2010-2019 2.825233 36353.7219 36432.2887
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 26530.49971 26576.81138 2.6734066 2000-2009 –> 2010-2019 5.389925 53803.4611 53896.2313
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + silicate + phosphate + phosphate_star 26553.48785 26593.18356 2.6794655 2000-2009 –> 2010-2019 5.403197 53854.3376 53933.8549
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + nitrate + silicate 16044.68809 16090.50992 1.1487205 2000-2009 –> 2010-2019 2.677515 35079.0719 35170.7332
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 26733.51411 26779.82577 2.7230317 2000-2009 –> 2010-2019 5.377131 53744.4477 53837.2179
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate 15972.75060 16018.57243 1.1407194 2000-2009 –> 2010-2019 2.533680 34047.0657 34138.7270
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate 26576.58612 26622.89778 2.6845920 2000-2009 –> 2010-2019 5.333669 53566.1718 53658.9419
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 15988.06133 16033.88316 1.1424176 2000-2009 –> 2010-2019 2.537504 34078.1123 34169.7736
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 26559.44237 26605.75403 2.6804257 2000-2009 –> 2010-2019 5.360206 53678.9188 53771.6890
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + nitrate + silicate 16801.44615 16840.72200 1.2366076 2000-2009 –> 2010-2019 2.855845 36426.8684 36505.4352
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 15953.05828 15998.88011 1.1385388 2000-2009 –> 2010-2019 2.542634 34109.5232 34201.1845
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate 16791.28179 16830.55764 1.2353869 2000-2009 –> 2010-2019 2.806176 36103.2757 36181.8426
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 15985.87494 16031.69677 1.1421749 2000-2009 –> 2010-2019 2.536967 34073.7495 34165.4108
obs Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 26555.45453 26601.76619 2.6794575 2000-2009 –> 2010-2019 5.353246 53649.7027 53742.4728
mod Indo-Pacific (28.1,28.2] cstar_tref ~ sal + temp + silicate + phosphate_star 16078.38613 16117.66198 1.1527118 2000-2009 –> 2010-2019 2.560463 34259.6844 34338.2512
obs Indo-Pacific (28.2,28.3] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 8089.10240 8127.99739 1.9970030 2000-2009 –> 2010-2019 4.706424 15311.4719 15387.5220
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 7971.02236 8009.91736 1.9363118 2000-2009 –> 2010-2019 4.638357 15185.2582 15261.3083
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 29586.69393 29639.57343 0.6903484 2000-2009 –> 2010-2019 1.362340 53857.5987 53962.1530
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 33226.89558 33279.77508 0.7854407 2000-2009 –> 2010-2019 1.566787 61078.3397 61182.8940
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + phosphate + phosphate_star 7970.01488 8003.35344 1.9368142 2000-2009 –> 2010-2019 4.638906 15182.3021 15247.4879
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 7971.49141 8010.38641 1.9365492 2000-2009 –> 2010-2019 4.622332 15167.7132 15243.7633
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + nitrate + silicate + phosphate 30522.43865 30567.76393 0.7136840 2000-2009 –> 2010-2019 1.469043 57568.5843 57658.2023
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 29909.15273 29962.03222 0.6982854 2000-2009 –> 2010-2019 1.375467 54362.7695 54467.3237
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + nitrate + phosphate 7996.57305 8035.46805 1.9492861 2000-2009 –> 2010-2019 4.653020 15212.6724 15288.7225
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + phosphate 8001.42220 8034.76076 1.9527788 2000-2009 –> 2010-2019 4.656513 15215.5218 15280.7076
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 7970.13862 8009.03362 1.9358646 2000-2009 –> 2010-2019 4.626647 15171.9097 15247.9598
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + silicate + phosphate 31158.98991 31211.86941 0.7299206 2000-2009 –> 2010-2019 1.571705 60779.8025 60884.3568
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + silicate + phosphate 7994.97517 8033.87016 1.9484722 2000-2009 –> 2010-2019 4.651946 15210.7878 15286.8379
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 31195.23378 31248.11328 0.7308591 2000-2009 –> 2010-2019 1.572993 60825.8841 60930.4383
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 8048.45676 8087.35175 1.9759001 2000-2009 –> 2010-2019 4.698769 15285.6004 15361.6505
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 30517.46878 30570.34828 0.7135077 2000-2009 –> 2010-2019 1.461873 57344.7211 57449.2753
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 31310.30413 31363.18363 0.7338466 2000-2009 –> 2010-2019 1.592989 61415.8244 61520.3786
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 31192.48597 31245.36546 0.7307879 2000-2009 –> 2010-2019 1.572877 60821.8913 60926.4455
obs Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 8036.79925 8075.69425 1.9698888 2000-2009 –> 2010-2019 4.688973 15269.7922 15345.8423
mod Indo-Pacific (28.2,28.3] cstar_tref ~ sal + temp + silicate + phosphate_star 31382.05536 31427.38064 0.7357678 2000-2009 –> 2010-2019 1.600101 61628.6337 61718.2516
mod Indo-Pacific (28.3,28.4] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 16564.43023 16603.70491 1.2088328 2000-2009 –> 2010-2019 2.135555 29601.9457 29680.1408
mod Indo-Pacific (28.3,28.4] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 16565.10853 16610.92899 1.2086775 2000-2009 –> 2010-2019 2.133042 29579.9230 29671.1506
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 16399.58800 16445.40846 1.1893908 2000-2009 –> 2010-2019 2.115657 29434.3291 29525.5567
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 2447.22720 2477.49771 2.1412132 2000-2009 –> 2010-2019 4.708203 5092.1302 5152.6587
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 2450.35180 2480.62231 2.1472166 2000-2009 –> 2010-2019 4.632251 5059.1086 5119.6370
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 2443.68487 2473.95538 2.1344275 2000-2009 –> 2010-2019 4.614406 5050.1726 5110.7011
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2489.29979 2519.57031 2.2234767 2000-2009 –> 2010-2019 4.712661 5099.9154 5160.4439
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 2444.44388 2474.71439 2.1358796 2000-2009 –> 2010-2019 4.616124 5051.0512 5111.5797
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + silicate + phosphate + phosphate_star 2487.31241 2513.25857 2.2235018 2000-2009 –> 2010-2019 4.713825 5096.4377 5148.3192
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + nitrate + phosphate 16382.15518 16427.97565 1.1873775 2000-2009 –> 2010-2019 2.110906 29388.1932 29479.4208
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 16441.76677 16487.58724 1.1942761 2000-2009 –> 2010-2019 2.122775 29499.8643 29591.0919
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 2447.80423 2478.07474 2.1423206 2000-2009 –> 2010-2019 4.706400 5091.4434 5151.9719
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + silicate + phosphate 2477.94089 2508.21140 2.2009604 2000-2009 –> 2010-2019 4.693730 5090.1599 5150.6883
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2485.27495 2515.54546 2.2154722 2000-2009 –> 2010-2019 4.707743 5097.2707 5157.7992
mod Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 16382.25384 16428.07430 1.1873889 2000-2009 –> 2010-2019 2.112002 29399.6681 29490.8957
obs Indo-Pacific (28.3,28.4] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 2480.39633 2510.66684 2.2058083 2000-2009 –> 2010-2019 4.695864 5091.4018 5151.9303
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + aou + nitrate + phosphate 16549.97356 16589.24825 1.2071357 2000-2009 –> 2010-2019 2.131119 29558.7788 29636.9738
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 16550.04573 16595.86619 1.2069095 2000-2009 –> 2010-2019 2.128205 29532.6009 29623.8285
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 16552.71791 16591.99260 1.2074577 2000-2009 –> 2010-2019 2.132578 29573.4535 29651.6486
mod Indo-Pacific (28.3,28.4] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 16552.85741 16598.67787 1.2072394 2000-2009 –> 2010-2019 2.129731 29547.9959 29639.2235

6.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)

6.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)

6.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 = ""))

7 Model diagnotics

7.1 Selection criterion vs predictors

The selection criterion (rmse) was plotted against the number of predictors (limited to 2 - 5).

7.1.1 All models

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

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

[[2]]

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

7.1.2 Best models

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

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

[[2]]

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

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

7.2.1 All models

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

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

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

7.2.2 Best models

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

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

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

7.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 3 4 10 3 5
total 92.00 77.00 74.00 85.00 114.00 102.00 104.00
Indo-Pacific - mod
(-Inf,26] 9 6 6 6 10 4 7
(26,26.5] 8 6 8 7 6 8 7
(26.5,26.75] 7 5 6 7 7 8 8
(26.75,27] 5 10 5 7 5 5 10
(27,27.25] 7 10 7 7 4 5 8
(27.25,27.5] 6 10 6 6 5 7 8
(27.5,27.75] 5 5 6 6 10 10 6
(27.75,27.85] 7 4 7 6 10 4 9
(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 4 5 4 10 10 8
(28.2,28.3] 4 6 6 6 10 10 6
(28.3,28.4] 7 10 9 7 4 3 7
total 88.00 88.00 84.00 87.00 107.00 95.00 113.00
Atlantic - obs
(-Inf,26] 8 9 6 6 7 5 8
(26,26.5] 8 10 8 7 4 4 7
(26.5,26.75] 7 9 6 7 4 4 10
(26.75,27] 4 4 4 9 5 10 10
(27,27.25] 6 5 5 8 5 9 10
(27.25,27.5] 6 4 7 7 9 4 10
(27.5,27.75] 7 3 9 7 4 10 7
(27.75,27.85] 6 5 8 6 7 10 6
(27.85,27.95] 5 4 7 6 10 10 5
(27.95,28.05] 5 5 6 7 10 10 5
(28.05,28.1] 5 6 3 10 10 6 7
(28.1,28.15] 3 6 5 10 7 6 10
(28.15,28.2] 7 4 7 6 4 9 10
(28.2, Inf] 6 4 8 9 10 5 6
total 83.00 78.00 89.00 105.00 96.00 102.00 111.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] 5 5 6 7 10 4 10
(26.75,27] 7 10 6 6 4 6 8
(27,27.25] 6 4 6 7 10 10 4
(27.25,27.5] 6 4 6 7 10 10 4
(27.5,27.75] 6 5 7 6 10 8 6
(27.75,27.85] 4 5 5 8 5 10 10
(27.85,27.95] 6 5 5 7 4 10 10
(27.95,28.05] 9 5 9 7 4 6 4
(28.05,28.1] 7 5 8 6 7 9 5
(28.1,28.2] 7 3 9 7 10 7 4
(28.2,28.3] 8 4 10 7 9 4 6
(28.3,28.4] 7 5 7 8 10 8 4
total 91.00 71.00 95.00 96.00 113.00 102.00 93.00

7.4 RMSE alternatives

7.4.1 AIC

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

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

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

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

7.4.2 AIC vs BIC

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

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

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

7.5 RMSE vs VIF

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

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

7.6 Residual patterns

7.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(cstar_tref, .fitted)) +
  geom_bin2d(binwidth = i_binwidth) +
  scale_fill_viridis_c() +
  geom_abline(slope = 1,
              col = 'red') +
  coord_equal(xlim = axis_lims,
              ylim = axis_lims) +
  labs(title = "All models") +
  facet_wrap(~data_source)

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

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

7.6.2 Pooled

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

i_ylim <- c(-30,30)

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

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

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

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

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

7.6.3 By model domain

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

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

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

[[2]]

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

[[3]]

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

[[4]]

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

7.6.4 Latitudinal mean

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

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

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

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

[[2]]

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

7.6.5 Latitudinal offset

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

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

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

8 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.13          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