eMLR - model fitting

workflowr

• Summary
• Checks
• Past versions

Last updated: 2021-06-11

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

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

Note that any generated files, e.g. HTML, png, CSS, etc., are not included in this status report because it is ok for generated content to have uncommitted changes.

---

These are the previous versions of the repository in which changes were made to the R Markdown (analysis/eMLR_model_fitting.Rmd) and HTML (docs/eMLR_model_fitting.html) files. If you’ve configured a remote Git repository (see ?wflow_git_remote), click on the hyperlinks in the table below to view the files as they were in that past version.

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

---

1 Version ID

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

2 Required data

Required are:

• cleaned and prepared GLODAPv2_2020 file

GLODAP <-
  read_csv(paste(path_version_data,
                 "GLODAPv2.2020_MLR_fitting_ready.csv",
                 sep = ""))

3 Predictor combinations

Find all possible combinations of following considered predictor variables:

• sal, temp, aou, nitrate, silicate, phosphate, phosphate_star

# the following code is a workaround to find all predictor combinations
# using the olsrr package and fit all models for one era, slab, and basin

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

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

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

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

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

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

# convert to tibble
lm_all <- as_tibble(lm_all)

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

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

4 Apply predictor threshold

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

• Minimum: 2
  
• Maximum: 5

lm_all <- lm_all %>% 
  filter(n >= params_local$MLR_predictors_min,
         n <= params_local$MLR_predictors_max)

This results in a total number of MLR models of:

• 112

5 Fit all models

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

• cstar_tref

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

library(furrr)
library(purrr)

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



library(tictoc)
plan(sequential)

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

plan(multisession, workers = 3)

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


future_map_dbl(1:4, function(x){
  Sys.sleep(1)
  x^2
}, .progress = TRUE)

# prepare nested data frame

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

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

# expand with model definitions

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

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


# test parallel processing

# expand with model definitions


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


### test furrr

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

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

6 Tidy models

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

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

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

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

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

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

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

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


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

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

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

7 Prepare coeffcients

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

7.1 VIF threshold

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

• 1000

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

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

GLODAP_glanced_clean %>% 
  group_by(data_source, era, basin, gamma_slab) %>% 
  count() %>% 
  filter(n < params_local$MLR_number)

# A tibble: 0 x 5
# Groups:   data_source, era, basin, gamma_slab [0]
# … with 5 variables: data_source <chr>, era <chr>, basin <chr>,
#   gamma_slab <chr>, n <int>

7.2 Predictor selection

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

• 10

The criterion used to select the best models was:

• rmse

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

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

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

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

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

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


# print table
lm_best_eras %>% 
  kable() %>%
  add_header_above() %>%
  kable_styling() %>%
  scroll_box(width = "100%", height = "400px")

data_source	basin	gamma_slab	model	AIC	BIC	rmse	eras	rmse_sum	aic_sum	bic_sum
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	111.51613	125.03891	0.6294568	2000-2009 –> 2010-2019	4.679849	767.9637	800.6398
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	850.04343	869.89741	6.6771693	2000-2009 –> 2010-2019	14.634901	1276.2272	1310.6239
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	61.25581	74.77859	0.3845635	2000-2009 –> 2010-2019	4.768801	735.7614	768.4375
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	83.26497	96.78775	0.4771758	2000-2009 –> 2010-2019	4.462325	736.0100	768.6862
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	842.62280	862.47677	6.4834137	2000-2009 –> 2010-2019	15.395343	1282.1697	1316.5664
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	59.54431	73.06709	0.3781646	2000-2009 –> 2010-2019	4.250879	705.7642	738.4404
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate	852.26730	869.28499	6.7900310	2000-2009 –> 2010-2019	14.843352	1277.8600	1307.3429
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	850.65190	870.50588	6.6933112	2000-2009 –> 2010-2019	14.682682	1277.3039	1311.7006
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	850.69933	870.55331	6.6945712	2000-2009 –> 2010-2019	14.680763	1277.3043	1311.7011
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + silicate	56.00567	69.52845	0.3652701	2000-2009 –> 2010-2019	4.657750	725.6888	758.3650
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + silicate	844.99527	864.84925	6.5447406	2000-2009 –> 2010-2019	14.530779	1271.5980	1305.9948
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + silicate + phosphate	81.55505	95.07783	0.4692432	2000-2009 –> 2010-2019	4.346449	728.0392	760.7153
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	82.03144	95.55422	0.4714399	2000-2009 –> 2010-2019	4.758121	751.4063	784.0825
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	110.93667	124.45945	0.6258910	2000-2009 –> 2010-2019	4.662031	766.5805	799.2567
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	850.91332	870.76730	6.7002583	2000-2009 –> 2010-2019	14.710031	1277.8662	1312.2630
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	55.80179	69.32457	0.3645407	2000-2009 –> 2010-2019	4.672665	726.3144	758.9906
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	81.36882	94.89160	0.4683873	2000-2009 –> 2010-2019	4.655780	745.4007	778.0768
obs	Atlantic	(-Inf,26]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	857.70216	877.55613	6.8832155	2000-2009 –> 2010-2019	15.374571	1291.5446	1325.9414
obs	Atlantic	(-Inf,26]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	856.78193	876.63590	6.8581257	2000-2009 –> 2010-2019	15.341317	1290.5109	1324.9076
obs	Atlantic	(-Inf,26]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	853.95890	873.81288	6.7817265	2000-2009 –> 2010-2019	15.228807	1287.1845	1321.5813
obs	Atlantic	(26,26.5]	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	3999.42220	4025.97141	6.1248842	2000-2009 –> 2010-2019	11.594247	7254.2468	7306.3190
mod	Atlantic	(26,26.5]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	2097.89682	2125.38530	3.8945984	2000-2009 –> 2010-2019	8.364974	3459.2787	3510.8641
obs	Atlantic	(26,26.5]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	4001.41992	4032.39400	6.1248729	2000-2009 –> 2010-2019	11.571505	7253.9133	7314.6642
mod	Atlantic	(26,26.5]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	2052.13712	2079.62561	3.6640813	2000-2009 –> 2010-2019	7.670457	3362.8906	3414.4760
obs	Atlantic	(26,26.5]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	3867.46896	3898.44304	5.4948307	2000-2009 –> 2010-2019	10.952592	7122.0852	7182.8361
mod	Atlantic	(26,26.5]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2023.98809	2051.47657	3.5291095	2000-2009 –> 2010-2019	7.920282	3377.1112	3428.6966
mod	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	2055.84501	2083.33349	3.6822408	2000-2009 –> 2010-2019	7.969869	3397.9436	3449.5290
obs	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	3884.71744	3915.69152	5.5721751	2000-2009 –> 2010-2019	11.098855	7152.3843	7213.1352
mod	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	2056.11954	2083.60802	3.6835889	2000-2009 –> 2010-2019	8.040996	3405.6764	3457.2618
obs	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	3885.65990	3916.63398	5.5764325	2000-2009 –> 2010-2019	11.176958	7167.1308	7227.8817
mod	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate	2028.09670	2055.58519	3.5484956	2000-2009 –> 2010-2019	8.286920	3416.3819	3467.9673
mod	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	2028.78633	2056.27481	3.5517600	2000-2009 –> 2010-2019	8.407466	3428.3673	3479.9527
mod	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	2057.12475	2084.61323	3.6885292	2000-2009 –> 2010-2019	8.118506	3414.3125	3465.8979
obs	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	3895.90297	3926.87705	5.6229134	2000-2009 –> 2010-2019	11.284425	7188.6377	7249.3886
obs	Atlantic	(26,26.5]	cstar_tref ~ temp + aou + nitrate + phosphate	3997.95985	4024.50907	6.1176302	2000-2009 –> 2010-2019	11.678795	7270.0958	7322.1680
mod	Atlantic	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	2095.43506	2122.92355	3.8818359	2000-2009 –> 2010-2019	8.392330	3460.9445	3512.5299
obs	Atlantic	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	3999.95985	4030.93393	6.1176302	2000-2009 –> 2010-2019	11.600203	7259.2935	7320.0444
mod	Atlantic	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	2091.75381	2119.24230	3.8628293	2000-2009 –> 2010-2019	8.409738	3460.9782	3512.5636
obs	Atlantic	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	3999.88143	4030.85552	6.1172415	2000-2009 –> 2010-2019	11.642650	7267.3089	7328.0598
obs	Atlantic	(26,26.5]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	4000.40583	4031.37991	6.1198416	2000-2009 –> 2010-2019	11.700267	7278.1375	7338.8884
mod	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	2753.04564	2782.30575	4.1227171	2000-2009 –> 2010-2019	7.290891	4654.9626	4711.5603
obs	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	4797.35857	4829.95566	5.2341233	2000-2009 –> 2010-2019	10.905470	8066.6634	8128.9832
mod	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	2756.11136	2785.37148	4.1358219	2000-2009 –> 2010-2019	7.194530	4632.2189	4688.8166
obs	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	4796.47431	4829.07139	5.2311497	2000-2009 –> 2010-2019	10.931675	8071.0751	8133.3950
mod	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	2802.69363	2831.95375	4.3401455	2000-2009 –> 2010-2019	7.192771	4627.6026	4684.2002
obs	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	4890.99412	4923.59121	5.5587678	2000-2009 –> 2010-2019	11.256637	8165.1139	8227.4338
mod	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	2808.88149	2838.14160	4.3680362	2000-2009 –> 2010-2019	7.166669	4619.7646	4676.3622
mod	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	2762.71088	2791.97100	4.1641737	2000-2009 –> 2010-2019	7.089474	4606.0852	4662.6828
obs	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	4793.58245	4826.17953	5.2214365	2000-2009 –> 2010-2019	10.920607	8067.9380	8130.2578
mod	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	2809.86182	2839.12194	4.3724713	2000-2009 –> 2010-2019	7.167166	4619.7114	4676.3090
obs	Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + phosphate	4858.47142	4886.41178	5.4507888	2000-2009 –> 2010-2019	11.126290	8126.5320	8179.9490
obs	Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + phosphate_star	4855.52835	4883.46870	5.4404887	2000-2009 –> 2010-2019	11.141396	8128.1982	8181.6152
mod	Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	2775.01041	2804.27052	4.2175327	2000-2009 –> 2010-2019	7.285669	4653.3769	4709.9746
obs	Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	4858.64797	4891.24506	5.4444048	2000-2009 –> 2010-2019	11.119184	8128.5772	8190.8970
mod	Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	2778.87604	2808.13616	4.2344438	2000-2009 –> 2010-2019	7.202484	4632.8968	4689.4945
obs	Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	4855.85653	4888.45361	5.4346464	2000-2009 –> 2010-2019	11.127736	8129.1103	8191.4301
mod	Atlantic	(26.5,26.75]	cstar_tref ~ temp + nitrate + phosphate + phosphate_star	2788.12272	2813.20282	4.2840315	2000-2009 –> 2010-2019	7.348945	4663.7178	4712.2301
obs	Atlantic	(26.5,26.75]	cstar_tref ~ temp + nitrate + phosphate + phosphate_star	4853.10056	4881.04092	5.4320067	2000-2009 –> 2010-2019	11.132230	8125.6467	8179.0637
mod	Atlantic	(26.5,26.75]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	2782.64341	2811.90353	4.2509903	2000-2009 –> 2010-2019	7.170927	4624.6708	4681.2685
obs	Atlantic	(26.5,26.75]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	4853.40118	4885.99826	5.4260773	2000-2009 –> 2010-2019	11.119079	8126.6391	8188.9589
mod	Atlantic	(26.75,27]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	4280.42042	4314.14523	2.4967490	2000-2009 –> 2010-2019	4.824464	7128.7303	7193.5306
mod	Atlantic	(26.75,27]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	4650.31885	4679.22584	3.0600615	2000-2009 –> 2010-2019	5.422729	7515.2886	7570.8317
mod	Atlantic	(26.75,27]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	4429.65429	4463.37910	2.7091292	2000-2009 –> 2010-2019	5.238579	7382.0234	7446.8236
mod	Atlantic	(26.75,27]	cstar_tref ~ sal + aou + silicate + phosphate	4578.18207	4607.08905	2.9416566	2000-2009 –> 2010-2019	5.477546	7531.7349	7587.2780
mod	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + aou + silicate + phosphate	4558.99204	4592.71685	2.9077540	2000-2009 –> 2010-2019	5.183143	7378.8362	7443.6365
obs	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + aou + silicate + phosphate	5689.34448	5723.52155	4.4437319	2000-2009 –> 2010-2019	10.822181	9714.6491	9779.7548
mod	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	4579.06299	4612.78780	2.9398563	2000-2009 –> 2010-2019	5.289313	7439.0126	7503.8129
obs	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	5694.19239	5728.36945	4.4547932	2000-2009 –> 2010-2019	10.777112	9708.6604	9773.7661
obs	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	5704.32493	5738.50199	4.4780013	2000-2009 –> 2010-2019	10.786611	9716.1316	9781.2373
mod	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	4586.33923	4620.06404	2.9515815	2000-2009 –> 2010-2019	5.334424	7463.9544	7528.7546
obs	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	5695.23802	5729.41508	4.4571826	2000-2009 –> 2010-2019	10.776375	9709.0997	9774.2053
obs	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + silicate + phosphate_star	5713.88207	5743.17670	4.5046201	2000-2009 –> 2010-2019	10.910392	9742.4273	9798.2321
mod	Atlantic	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	4299.90550	4333.63031	2.5235048	2000-2009 –> 2010-2019	4.799211	7119.9243	7184.7245
obs	Atlantic	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	5694.92386	5729.10092	4.4564646	2000-2009 –> 2010-2019	10.925441	9737.5064	9802.6121
mod	Atlantic	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	4309.88827	4343.61308	2.5373234	2000-2009 –> 2010-2019	4.829968	7139.1914	7203.9917
obs	Atlantic	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	5700.70402	5734.88108	4.4696939	2000-2009 –> 2010-2019	10.851194	9726.5949	9791.7005
mod	Atlantic	(26.75,27]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	4332.34561	4366.07042	2.5686872	2000-2009 –> 2010-2019	4.894538	7179.6520	7244.4523
obs	Atlantic	(26.75,27]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	5702.85912	5737.03619	4.4746365	2000-2009 –> 2010-2019	10.867689	9730.9675	9796.0731
obs	Atlantic	(26.75,27]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	5707.73899	5737.03362	4.4904515	2000-2009 –> 2010-2019	10.914365	9739.7514	9795.5562
obs	Atlantic	(26.75,27]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	5703.08028	5732.37491	4.4797363	2000-2009 –> 2010-2019	10.931155	9740.3308	9796.1356
mod	Atlantic	(27,27.25]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	5425.89355	5461.87672	2.0652160	2000-2009 –> 2010-2019	4.383652	9412.6773	9482.1120
mod	Atlantic	(27,27.25]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	5573.45586	5609.43903	2.1895554	2000-2009 –> 2010-2019	4.668471	9677.8996	9747.3343
mod	Atlantic	(27,27.25]	cstar_tref ~ sal + aou + silicate + phosphate	5720.41995	5751.26267	2.3226704	2000-2009 –> 2010-2019	4.786122	9811.8622	9871.3776
mod	Atlantic	(27,27.25]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	5721.53189	5757.51507	2.3218533	2000-2009 –> 2010-2019	4.666726	9728.2477	9797.6824
obs	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	11086.94927	11125.05901	6.1642466	2000-2009 –> 2010-2019	12.061644	19059.0260	19133.0296
obs	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	11061.13140	11099.24114	6.1178874	2000-2009 –> 2010-2019	11.901259	18984.5532	19058.5568
mod	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + aou + nitrate + silicate	5524.22224	5560.20541	2.1472593	2000-2009 –> 2010-2019	4.412224	9469.9856	9539.4202
mod	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + aou + silicate + phosphate	5720.46731	5756.45048	2.3208742	2000-2009 –> 2010-2019	4.700448	9753.0087	9822.4434
obs	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	11105.93921	11144.04895	6.1985695	2000-2009 –> 2010-2019	11.937655	19010.2055	19084.2091
obs	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	11068.02731	11106.13705	6.1302357	2000-2009 –> 2010-2019	11.946684	19005.6610	19079.6646
obs	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	11124.54003	11162.64977	6.2323744	2000-2009 –> 2010-2019	12.067045	19069.9689	19143.9725
obs	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	11107.60171	11145.71145	6.2015835	2000-2009 –> 2010-2019	11.953959	19017.6317	19091.6353
obs	Atlantic	(27,27.25]	cstar_tref ~ temp + aou + nitrate + phosphate_star	11081.77723	11114.44272	6.1585320	2000-2009 –> 2010-2019	12.069662	19057.6498	19121.0814
mod	Atlantic	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate	5689.28557	5720.12829	2.2941954	2000-2009 –> 2010-2019	4.597162	9662.2995	9721.8149
mod	Atlantic	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	5690.73613	5726.71930	2.2936961	2000-2009 –> 2010-2019	4.554682	9633.4084	9702.8431
obs	Atlantic	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	11047.86453	11085.97427	6.0942009	2000-2009 –> 2010-2019	12.069122	19052.4859	19126.4895
mod	Atlantic	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	5691.06141	5727.04459	2.2939917	2000-2009 –> 2010-2019	4.553802	9632.8189	9702.2535
obs	Atlantic	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	11007.04811	11045.15785	6.0219012	2000-2009 –> 2010-2019	11.864577	18955.8935	19029.8971
mod	Atlantic	(27,27.25]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	5606.79449	5642.77766	2.2186683	2000-2009 –> 2010-2019	4.480607	9550.2075	9619.6422
obs	Atlantic	(27,27.25]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	11018.23444	11056.34418	6.0416302	2000-2009 –> 2010-2019	11.918514	18981.6276	19055.6312
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	13195.89388	13235.34100	5.8421203	2000-2009 –> 2010-2019	10.522926	21889.7516	21966.2260
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + aou + nitrate + silicate	11037.45394	11071.20458	3.5662299	2000-2009 –> 2010-2019	6.894367	20706.5123	20773.3746
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	10896.49756	10935.87331	3.4439689	2000-2009 –> 2010-2019	6.734159	20525.3094	20603.3154
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	10909.32399	10948.69974	3.4547652	2000-2009 –> 2010-2019	6.766793	20562.5068	20640.5128
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	13180.81354	13220.26067	5.8208786	2000-2009 –> 2010-2019	10.423087	21825.0550	21901.5294
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	11157.90501	11197.28076	3.6708146	2000-2009 –> 2010-2019	6.936162	20758.7952	20836.8012
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	13179.38847	13218.83560	5.8188753	2000-2009 –> 2010-2019	10.396076	21807.6651	21884.1395
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + nitrate + silicate	11009.17829	11048.55404	3.5399800	2000-2009 –> 2010-2019	6.788182	20590.6746	20668.6806
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + phosphate	13246.91757	13280.72940	5.9174254	2000-2009 –> 2010-2019	10.534285	21898.4718	21964.0213
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + phosphate_star	13250.32640	13284.13822	5.9222998	2000-2009 –> 2010-2019	10.516521	21887.4783	21953.0277
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate	13243.82005	13283.26718	5.9101439	2000-2009 –> 2010-2019	10.525052	21896.1356	21972.6100
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	10888.63893	10928.01468	3.4373708	2000-2009 –> 2010-2019	7.102061	20914.5809	20992.5869
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	13247.02073	13286.46786	5.9147148	2000-2009 –> 2010-2019	10.508828	21886.1034	21962.5778
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	11153.46905	11192.84480	3.6668432	2000-2009 –> 2010-2019	6.799450	20601.4709	20679.4769
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	13176.34751	13215.79464	5.8146027	2000-2009 –> 2010-2019	10.396936	21807.9078	21884.3822
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	10866.78992	10906.16567	3.4190928	2000-2009 –> 2010-2019	6.956689	20762.6995	20840.7055
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + phosphate + phosphate_star	13247.08801	13280.89983	5.9176690	2000-2009 –> 2010-2019	10.515706	21886.6724	21952.2219
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	10889.91629	10929.29204	3.4384424	2000-2009 –> 2010-2019	7.105563	20918.3009	20996.3069
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	13243.97440	13283.42153	5.9103642	2000-2009 –> 2010-2019	10.507924	21885.2547	21961.7291
mod	Atlantic	(27.25,27.5]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	11164.22607	11203.60182	3.6764811	2000-2009 –> 2010-2019	6.802414	20604.3715	20682.3775
obs	Atlantic	(27.5,27.75]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	14154.22725	14194.47251	5.0963482	2000-2009 –> 2010-2019	10.367132	25278.7373	25357.4631
obs	Atlantic	(27.5,27.75]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	14164.14047	14198.63641	5.1094499	2000-2009 –> 2010-2019	10.381377	25287.4330	25354.9122
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + aou + nitrate + phosphate	18254.11365	18292.09883	2.1790881	2000-2009 –> 2010-2019	4.614463	33879.8363	33954.5770
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	17816.88643	17861.20248	2.0667711	2000-2009 –> 2010-2019	4.495410	33425.8812	33513.0787
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + aou + silicate + phosphate	18566.42047	18604.40566	2.2626435	2000-2009 –> 2010-2019	4.704592	34210.3719	34285.1127
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	18533.60484	18577.92089	2.2531724	2000-2009 –> 2010-2019	4.687930	34159.6128	34246.8103
obs	Atlantic	(27.5,27.75]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	14165.80459	14206.04985	5.1090800	2000-2009 –> 2010-2019	10.305216	25238.8795	25317.6052
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	18142.76979	18187.08583	2.1495330	2000-2009 –> 2010-2019	4.569233	33726.8288	33814.0264
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	18564.92605	18609.24209	2.2616911	2000-2009 –> 2010-2019	4.683981	34156.2205	34243.4180
obs	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	14171.47060	14211.71586	5.1153226	2000-2009 –> 2010-2019	10.318819	25249.6496	25328.3753
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	18373.36336	18417.67940	2.2100894	2000-2009 –> 2010-2019	4.590018	33845.3549	33932.5525
obs	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	14200.28685	14240.53211	5.1471897	2000-2009 –> 2010-2019	10.378305	25297.5553	25376.2810
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	18456.79044	18501.10648	2.2324160	2000-2009 –> 2010-2019	4.621898	33955.8727	34043.0702
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	18483.61399	18527.93004	2.2396423	2000-2009 –> 2010-2019	4.627205	33977.2621	34064.4597
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	18378.76421	18423.08025	2.2115280	2000-2009 –> 2010-2019	4.593222	33855.7700	33942.9675
obs	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	14180.41961	14220.66487	5.1251979	2000-2009 –> 2010-2019	10.341281	25267.3110	25346.0367
obs	Atlantic	(27.5,27.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	14160.90925	14201.15450	5.1036926	2000-2009 –> 2010-2019	10.391462	25297.0210	25375.7467
obs	Atlantic	(27.5,27.75]	cstar_tref ~ temp + aou + silicate + phosphate	14169.72961	14204.22554	5.1156082	2000-2009 –> 2010-2019	10.403983	25304.2543	25371.7335
obs	Atlantic	(27.5,27.75]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	14170.72842	14210.97368	5.1145045	2000-2009 –> 2010-2019	10.426267	25323.1660	25401.8917
obs	Atlantic	(27.5,27.75]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	14179.08775	14213.58369	5.1259360	2000-2009 –> 2010-2019	10.437868	25329.6403	25397.1194
obs	Atlantic	(27.75,27.85]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	9247.26510	9284.40612	5.3738975	2000-2009 –> 2010-2019	10.551043	16974.3833	17047.4908
obs	Atlantic	(27.75,27.85]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	9245.26632	9277.10148	5.3738997	2000-2009 –> 2010-2019	10.575865	16982.4275	17045.0911
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + nitrate + silicate	5218.61114	5250.58176	1.3369472	2000-2009 –> 2010-2019	2.574440	9490.1978	9553.2122
obs	Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	9234.02262	9271.16364	5.3500541	2000-2009 –> 2010-2019	10.543296	16968.9576	17042.0651
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + silicate + phosphate	5005.30788	5037.27850	1.2465273	2000-2009 –> 2010-2019	2.432958	9166.9143	9229.9287
obs	Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + silicate + phosphate	9232.08400	9263.91916	5.3501644	2000-2009 –> 2010-2019	10.571233	16978.4753	17041.1389
obs	Atlantic	(27.75,27.85]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	9235.41382	9272.55484	5.3525540	2000-2009 –> 2010-2019	10.577176	16985.5181	17058.6256
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + nitrate + silicate + phosphate_star	5613.88816	5645.85878	1.5222021	2000-2009 –> 2010-2019	2.887847	10142.6640	10205.6784
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	5605.56069	5637.53132	1.5180462	2000-2009 –> 2010-2019	2.889565	10145.5390	10208.5533
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + aou + nitrate + silicate	5176.14954	5213.44860	1.3175739	2000-2009 –> 2010-2019	2.551624	9442.4665	9515.9832
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + aou + silicate + phosphate	4889.18592	4926.48498	1.1991130	2000-2009 –> 2010-2019	2.380826	9042.3944	9115.9112
obs	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + aou + silicate + phosphate	9232.98948	9270.13050	5.3481983	2000-2009 –> 2010-2019	10.558418	16976.1430	17049.2505
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	5017.51753	5054.81659	1.2507124	2000-2009 –> 2010-2019	2.471550	9255.7384	9329.2551
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	5022.44929	5059.74835	1.2527391	2000-2009 –> 2010-2019	2.472194	9257.7134	9331.2302
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	5011.48395	5048.78301	1.2482374	2000-2009 –> 2010-2019	2.465100	9241.1935	9314.7103
obs	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	9234.63831	9271.77933	5.3511603	2000-2009 –> 2010-2019	10.570896	16982.3864	17055.4939
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + silicate + phosphate_star	5098.72166	5130.69229	1.2853476	2000-2009 –> 2010-2019	2.589623	9507.4923	9570.5067
obs	Atlantic	(27.75,27.85]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	9247.31549	9284.45651	5.3739885	2000-2009 –> 2010-2019	10.558070	16977.8053	17050.9128
obs	Atlantic	(27.75,27.85]	cstar_tref ~ temp + aou + silicate + phosphate	9245.31731	9277.15247	5.3739917	2000-2009 –> 2010-2019	10.584227	16986.4782	17049.1418
obs	Atlantic	(27.75,27.85]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	9247.66438	9284.80540	5.3746181	2000-2009 –> 2010-2019	10.566543	16981.9607	17055.0682
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + aou + nitrate + silicate	6319.86372	6351.74305	1.9811603	2000-2009 –> 2010-2019	4.900152	12826.2451	12889.1636
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	19492.92007	19535.31062	5.3174384	2000-2009 –> 2010-2019	10.588219	36163.4386	36247.1439
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	6320.36916	6357.56170	1.9801736	2000-2009 –> 2010-2019	4.821533	12758.4502	12831.8550
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	18891.13878	18933.52933	4.8333096	2000-2009 –> 2010-2019	9.870505	35316.6088	35400.3141
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + aou + silicate + phosphate	6096.37290	6128.25222	1.8389335	2000-2009 –> 2010-2019	4.573538	12432.5783	12495.4968
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + aou + silicate + phosphate_star	18895.62034	18931.95509	4.8382816	2000-2009 –> 2010-2019	9.875752	35319.3850	35391.1325
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	19291.28314	19333.67369	5.1500485	2000-2009 –> 2010-2019	10.343805	35882.2191	35965.9245
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	19333.12611	19369.46087	5.1859908	2000-2009 –> 2010-2019	10.383272	35925.7293	35997.4767
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + aou + nitrate + silicate	6161.54623	6198.73878	1.8780679	2000-2009 –> 2010-2019	4.784306	12658.5070	12731.9119
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + aou + phosphate	6554.82631	6586.70563	2.1425647	2000-2009 –> 2010-2019	4.894394	12907.4078	12970.3263
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate + phosphate	5798.93091	5836.12345	1.6642453	2000-2009 –> 2010-2019	4.398612	12136.9093	12210.3142
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate + phosphate	18729.46997	18771.86052	4.7109332	2000-2009 –> 2010-2019	9.717362	35121.8196	35205.5250
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	5902.60629	5939.79883	1.7227644	2000-2009 –> 2010-2019	4.481952	12264.1517	12337.5565
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	18710.72029	18753.11083	4.6969425	2000-2009 –> 2010-2019	9.698610	35097.9259	35181.6313
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate	18710.06518	18752.45572	4.6964544	2000-2009 –> 2010-2019	9.741798	35144.2724	35227.9778
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	5910.84027	5948.03281	1.7274993	2000-2009 –> 2010-2019	4.530939	12313.8812	12387.2861
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + silicate + phosphate	18734.22519	18770.55994	4.7159840	2000-2009 –> 2010-2019	9.762702	35167.9049	35239.6523
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	5892.33074	5929.52328	1.7168737	2000-2009 –> 2010-2019	4.467528	12245.7978	12319.2027
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	18728.19883	18770.58938	4.7099834	2000-2009 –> 2010-2019	9.715257	35119.3006	35203.0060
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + silicate + phosphate_star	5992.75693	6024.63625	1.7765035	2000-2009 –> 2010-2019	4.809054	12598.6747	12661.5931
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	19948.18248	19990.70066	5.3982472	2000-2009 –> 2010-2019	12.105512	37715.6436	37799.3959
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	19533.40136	19575.91954	5.0605069	2000-2009 –> 2010-2019	11.622625	37183.9482	37267.7004
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + aou + silicate + phosphate_star	19619.58372	19656.02788	5.1304955	2000-2009 –> 2010-2019	11.706939	37279.7841	37351.5717
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	19669.48093	19711.99911	5.1689153	2000-2009 –> 2010-2019	11.775919	37356.4564	37440.2086
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	19798.27113	19834.71528	5.2752980	2000-2009 –> 2010-2019	11.896184	37494.4636	37566.2513
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + aou + phosphate	9622.93704	9655.85672	3.5776628	2000-2009 –> 2010-2019	7.378045	17467.4274	17531.9060
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + aou + phosphate_star	9652.77286	9685.69254	3.6077049	2000-2009 –> 2010-2019	7.406908	17496.3807	17560.8593
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate + phosphate	9486.46434	9524.87064	3.4414749	2000-2009 –> 2010-2019	7.240996	17332.3105	17407.5355
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate + phosphate	19651.79877	19694.31695	5.1546985	2000-2009 –> 2010-2019	11.753860	37332.4278	37416.1800
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	9498.23017	9536.63647	3.4528422	2000-2009 –> 2010-2019	7.249634	17342.0325	17417.2575
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	19597.22855	19639.74674	5.1110691	2000-2009 –> 2010-2019	11.686971	37258.9887	37342.7409
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + nitrate + phosphate_star	9717.42429	9750.34397	3.6736716	2000-2009 –> 2010-2019	7.459780	17551.2129	17615.6915
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate	19681.05404	19723.57222	5.1782415	2000-2009 –> 2010-2019	11.893560	37454.9286	37538.6808
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	9472.47336	9510.87965	3.4280065	2000-2009 –> 2010-2019	7.178110	17281.0866	17356.3116
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + phosphate + phosphate_star	9647.96734	9680.88702	3.6028491	2000-2009 –> 2010-2019	7.402051	17491.5741	17556.0527
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + phosphate_star	9730.62522	9758.05829	3.6893561	2000-2009 –> 2010-2019	7.491447	17574.3936	17628.1257
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + silicate + phosphate	19784.16854	19820.61270	5.2637226	2000-2009 –> 2010-2019	11.983532	37559.6159	37631.4036
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	9498.47694	9536.88323	3.4530810	2000-2009 –> 2010-2019	7.249878	17342.2832	17417.5082
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	19636.04712	19678.56530	5.1420668	2000-2009 –> 2010-2019	11.732551	37309.6439	37393.3962
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + silicate + phosphate_star	9497.08253	9530.00221	3.4536672	2000-2009 –> 2010-2019	7.250472	17338.8946	17403.3732
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	14408.21610	14448.50946	5.2699204	2000-2009 –> 2010-2019	11.607321	26001.3499	26080.0063
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	14352.59855	14392.89192	5.2075569	2000-2009 –> 2010-2019	11.544754	25945.6190	26024.2754
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	14355.72070	14396.01406	5.2110381	2000-2009 –> 2010-2019	11.535490	25941.6021	26020.2584
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	14441.59079	14476.12796	5.3099738	2000-2009 –> 2010-2019	11.657934	26038.6281	26106.0478
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	14474.73618	14515.02955	5.3454904	2000-2009 –> 2010-2019	11.650688	26049.8055	26128.4619
mod	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + aou + phosphate	3563.50196	3593.17804	1.3366383	2000-2009 –> 2010-2019	2.330493	5731.3960	5788.8957
mod	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + aou + phosphate_star	3602.88882	3632.56490	1.3622148	2000-2009 –> 2010-2019	2.374451	5798.7474	5856.2470
mod	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + aou + silicate + phosphate	3564.13164	3598.75374	1.3357571	2000-2009 –> 2010-2019	2.311004	5705.1838	5772.2667
mod	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	3601.13412	3635.75622	1.3597557	2000-2009 –> 2010-2019	2.358091	5777.8922	5844.9751
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	14389.08209	14429.37546	5.2483817	2000-2009 –> 2010-2019	11.519803	25945.1049	26023.7613
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	14467.67248	14507.96585	5.3374145	2000-2009 –> 2010-2019	11.645893	26044.5867	26123.2431
mod	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + nitrate + phosphate_star	3649.35750	3679.03358	1.3930201	2000-2009 –> 2010-2019	2.416257	5861.7122	5919.2118
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + nitrate + phosphate_star	14476.26769	14510.80486	5.3495325	2000-2009 –> 2010-2019	11.658439	26051.4223	26118.8420
mod	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	3622.54142	3657.16352	1.3738361	2000-2009 –> 2010-2019	2.396251	5835.6684	5902.7513
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	14356.25100	14396.54436	5.2116296	2000-2009 –> 2010-2019	11.502710	25923.3721	26002.0285
mod	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + phosphate + phosphate_star	3591.65992	3621.33600	1.3548737	2000-2009 –> 2010-2019	2.362706	5780.8663	5838.3659
mod	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	3591.54991	3626.17201	1.3534986	2000-2009 –> 2010-2019	2.345352	5758.3667	5825.4496
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	14430.40991	14470.70328	5.2950141	2000-2009 –> 2010-2019	11.592434	26001.1026	26079.7589
mod	Atlantic	(28.05,28.1]	cstar_tref ~ temp + aou + phosphate	3580.45963	3605.18970	1.3488882	2000-2009 –> 2010-2019	2.426716	5870.1295	5918.0459
mod	Atlantic	(28.05,28.1]	cstar_tref ~ temp + aou + silicate + phosphate	3582.44244	3612.11853	1.3488771	2000-2009 –> 2010-2019	2.402198	5839.0147	5896.5143
obs	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	16565.64633	16606.65140	5.9375848	2000-2009 –> 2010-2019	15.026695	27624.4896	27702.7937
mod	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + aou + phosphate	4574.70929	4606.71135	1.0737627	2000-2009 –> 2010-2019	2.400491	8946.3480	9009.2824
mod	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + aou + phosphate_star	4622.60589	4654.60795	1.0906908	2000-2009 –> 2010-2019	2.425717	9010.2181	9073.1526
mod	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + aou + silicate + phosphate	4544.78780	4582.12354	1.0626268	2000-2009 –> 2010-2019	2.270803	8678.6141	8752.0376
mod	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	4601.60905	4638.94479	1.0825300	2000-2009 –> 2010-2019	2.304659	8764.8526	8838.2762
obs	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	16600.12599	16641.13107	5.9773005	2000-2009 –> 2010-2019	15.082661	27664.4105	27742.7146
obs	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	16548.83387	16589.83895	5.9183150	2000-2009 –> 2010-2019	15.001276	27605.6157	27683.9199
mod	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + nitrate + phosphate_star	4642.38079	4674.38285	1.0977575	2000-2009 –> 2010-2019	2.407295	8980.6058	9043.5402
mod	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	4643.22697	4680.56271	1.0973439	2000-2009 –> 2010-2019	2.376674	8923.6614	8997.0849
obs	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	16574.26893	16615.27401	5.9474920	2000-2009 –> 2010-2019	15.056600	27639.8068	27718.1110
mod	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + phosphate + phosphate_star	4607.47040	4639.47246	1.0853128	2000-2009 –> 2010-2019	2.414680	8984.2002	9047.1346
obs	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + phosphate + phosphate_star	16601.32742	16636.47462	5.9810015	2000-2009 –> 2010-2019	15.084799	27663.0888	27730.2066
mod	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	4580.89804	4618.23378	1.0752326	2000-2009 –> 2010-2019	2.286627	8721.5390	8794.9625
obs	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	16600.07060	16641.07567	5.9772365	2000-2009 –> 2010-2019	15.079252	27663.2358	27741.5400
obs	Atlantic	(28.1,28.15]	cstar_tref ~ temp + aou + nitrate + phosphate_star	16592.67751	16627.82472	5.9710070	2000-2009 –> 2010-2019	15.074475	27654.3287	27721.4465
obs	Atlantic	(28.1,28.15]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	16568.32225	16609.32732	5.9406576	2000-2009 –> 2010-2019	15.043379	27631.7238	27710.0279
mod	Atlantic	(28.1,28.15]	cstar_tref ~ temp + aou + silicate + phosphate	4673.93368	4705.93574	1.1091280	2000-2009 –> 2010-2019	2.395484	8966.4007	9029.3351
obs	Atlantic	(28.1,28.15]	cstar_tref ~ temp + nitrate + phosphate + phosphate_star	16579.61253	16614.75974	5.9559427	2000-2009 –> 2010-2019	15.054172	27639.5103	27706.6281
obs	Atlantic	(28.1,28.15]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	16556.94417	16597.94924	5.9276029	2000-2009 –> 2010-2019	15.025234	27618.6420	27696.9461
mod	Atlantic	(28.1,28.15]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	4721.04461	4753.04666	1.1263246	2000-2009 –> 2010-2019	2.418188	9024.4560	9087.3905
mod	Atlantic	(28.15,28.2]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	7262.38485	7304.17519	0.8468827	2000-2009 –> 2010-2019	1.897848	13656.1891	13737.7633
mod	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + aou + nitrate	7322.04748	7357.86778	0.8559563	2000-2009 –> 2010-2019	1.898161	13677.4919	13747.4126
mod	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + aou + nitrate + silicate	7224.84215	7266.63249	0.8414054	2000-2009 –> 2010-2019	1.879774	13566.2670	13647.8411
obs	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + aou + silicate + phosphate	5567.35865	5601.06045	5.0988144	2000-2009 –> 2010-2019	12.844662	10191.2656	10256.4660
obs	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	5527.35203	5561.05383	4.9880772	2000-2009 –> 2010-2019	12.651832	10137.0880	10202.2883
mod	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + nitrate + phosphate_star	7635.69073	7671.51102	0.9036360	2000-2009 –> 2010-2019	1.964038	14066.3263	14136.2470
mod	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	7598.82006	7640.61041	0.8975857	2000-2009 –> 2010-2019	1.951570	14005.0843	14086.6585
obs	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	5512.92096	5546.62276	4.9487254	2000-2009 –> 2010-2019	12.661607	10131.1554	10196.3557
obs	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	5543.74231	5577.44411	5.0331512	2000-2009 –> 2010-2019	12.708591	10155.5046	10220.7049
mod	Atlantic	(28.15,28.2]	cstar_tref ~ temp + aou + nitrate + phosphate	7328.29583	7364.11613	0.8568812	2000-2009 –> 2010-2019	1.898208	13680.0793	13750.0001
mod	Atlantic	(28.15,28.2]	cstar_tref ~ temp + aou + nitrate + phosphate_star	7361.65325	7397.47355	0.8618355	2000-2009 –> 2010-2019	1.908715	13736.5373	13806.4580
mod	Atlantic	(28.15,28.2]	cstar_tref ~ temp + aou + nitrate + silicate	7252.62381	7288.44411	0.8457474	2000-2009 –> 2010-2019	1.896370	13643.0126	13712.9333
obs	Atlantic	(28.15,28.2]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	5528.10408	5561.80588	4.9901365	2000-2009 –> 2010-2019	12.994807	10195.7259	10260.9262
obs	Atlantic	(28.15,28.2]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	5499.72187	5533.42367	4.9130050	2000-2009 –> 2010-2019	12.818986	10150.8441	10216.0444
obs	Atlantic	(28.15,28.2]	cstar_tref ~ temp + aou + silicate + phosphate_star	5527.88762	5556.77488	4.9950237	2000-2009 –> 2010-2019	12.922110	10180.5556	10236.4416
mod	Atlantic	(28.15,28.2]	cstar_tref ~ temp + nitrate + phosphate + phosphate_star	7412.54073	7448.36103	0.8694487	2000-2009 –> 2010-2019	1.920307	13803.9024	13873.8231
obs	Atlantic	(28.15,28.2]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	5506.16921	5539.87101	4.9304210	2000-2009 –> 2010-2019	12.865804	10162.2285	10227.4289
mod	Atlantic	(28.15,28.2]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	7615.87209	7651.69238	0.9005461	2000-2009 –> 2010-2019	1.954677	14020.7409	14090.6616
obs	Atlantic	(28.15,28.2]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	5513.88129	5542.76854	4.9567725	2000-2009 –> 2010-2019	12.982950	10183.0716	10238.9576
obs	Atlantic	(28.15,28.2]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	5543.63821	5572.52547	5.0383912	2000-2009 –> 2010-2019	12.986718	10199.8651	10255.7511
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate	10053.86511	10086.25737	0.6873733	2000-2009 –> 2010-2019	1.438728	15899.1305	15960.7908
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate + phosphate	9991.39085	10030.26156	0.6827819	2000-2009 –> 2010-2019	1.432030	15824.1947	15898.1870
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	10756.56707	10796.06138	3.1850016	2000-2009 –> 2010-2019	6.499524	18334.4252	18410.8506
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate + phosphate_star	9906.75750	9945.62822	0.6768014	2000-2009 –> 2010-2019	1.426066	15739.6801	15813.6725
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate + silicate	9944.15703	9983.02775	0.6794377	2000-2009 –> 2010-2019	1.414546	15678.8438	15752.8361
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	10755.15331	10794.64762	3.1839215	2000-2009 –> 2010-2019	6.500712	18334.9891	18411.4145
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	9302.56565	9347.91481	0.6354695	2000-2009 –> 2010-2019	1.358632	14954.8746	15041.1990
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + phosphate + phosphate_star	10757.88507	10791.73734	3.1875381	2000-2009 –> 2010-2019	6.503475	18334.9763	18400.4837
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + silicate + phosphate	10753.74245	10787.59471	3.1843715	2000-2009 –> 2010-2019	6.501474	18331.8493	18397.3567
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	10735.70753	10775.20185	3.1691015	2000-2009 –> 2010-2019	6.482906	18312.9403	18389.3656
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + aou + nitrate	9972.94617	10011.81689	0.6814740	2000-2009 –> 2010-2019	1.431354	15810.0898	15884.0822
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	9908.76132	9954.11049	0.6768016	2000-2009 –> 2010-2019	1.425827	15742.0384	15828.3628
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	10754.78925	10794.28356	3.1836434	2000-2009 –> 2010-2019	6.502118	18336.0913	18412.5167
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	9903.57275	9948.92192	0.6764367	2000-2009 –> 2010-2019	1.425676	15738.3187	15824.6431
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + aou + nitrate + silicate	9865.00528	9910.35444	0.6737302	2000-2009 –> 2010-2019	1.406501	15585.2894	15671.6138
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + aou + silicate + phosphate	10726.14697	10765.64128	3.1618405	2000-2009 –> 2010-2019	6.478927	18306.2406	18382.6659
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	10769.26832	10808.76263	3.1947222	2000-2009 –> 2010-2019	6.469184	18311.9847	18388.4101
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	9948.83647	9994.18563	0.6796269	2000-2009 –> 2010-2019	1.431765	15803.4660	15889.7904
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	10796.84785	10836.34216	3.2159316	2000-2009 –> 2010-2019	6.497709	18346.0130	18422.4384
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	10742.11184	10781.60615	3.1739747	2000-2009 –> 2010-2019	6.455773	18291.2956	18367.7210
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	22564.13056	22607.49789	5.4318652	2000-2009 –> 2010-2019	10.970950	38075.4317	38159.4970
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	29326.91216	29372.06654	5.5517976	2000-2009 –> 2010-2019	11.850122	50317.2665	50404.9564
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + aou + phosphate + phosphate_star	22895.11923	22932.29123	5.6872380	2000-2009 –> 2010-2019	11.582173	38712.6280	38784.6840
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	22674.33566	22717.70300	5.5150872	2000-2009 –> 2010-2019	11.169340	38287.5009	38371.5662
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	29644.81478	29689.96916	5.7436807	2000-2009 –> 2010-2019	12.265789	50859.8758	50947.5658
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate	30107.80682	30146.51057	6.0363536	2000-2009 –> 2010-2019	12.519698	51282.5014	51357.6642
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	22672.59203	22715.95936	5.5137606	2000-2009 –> 2010-2019	11.330137	38425.6914	38509.7567
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	29786.02010	29831.17448	5.8310246	2000-2009 –> 2010-2019	12.151685	50799.1580	50886.8480
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	22775.10629	22818.47362	5.5923003	2000-2009 –> 2010-2019	11.511459	38614.9141	38698.9794
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	29839.32647	29884.48085	5.8643420	2000-2009 –> 2010-2019	12.195663	50863.3101	50951.0001
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + silicate	23245.90883	23289.27616	5.9676123	2000-2009 –> 2010-2019	12.158391	39307.8062	39391.8715
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + silicate	30026.89762	30072.05200	5.9830981	2000-2009 –> 2010-2019	12.463199	51200.3719	51288.0619
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + silicate + phosphate	23124.75584	23168.12318	5.8686905	2000-2009 –> 2010-2019	11.965542	39110.9753	39195.0406
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + silicate + phosphate	30107.65161	30152.80599	6.0349632	2000-2009 –> 2010-2019	12.651568	51415.2924	51502.9824
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	23322.73428	23366.10161	6.0312027	2000-2009 –> 2010-2019	12.304625	39450.2739	39534.3392
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	30219.37011	30264.52449	6.1074577	2000-2009 –> 2010-2019	12.760547	51562.4018	51650.0917
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	22875.52962	22918.89696	5.6703226	2000-2009 –> 2010-2019	11.681479	38791.6801	38875.7454
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	29967.89496	30013.04934	5.9454850	2000-2009 –> 2010-2019	12.321085	51036.7327	51124.4227
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	23391.71885	23435.08618	6.0888801	2000-2009 –> 2010-2019	12.398712	39547.9042	39631.9695
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	30268.44220	30313.59658	6.1395754	2000-2009 –> 2010-2019	12.808017	51626.3079	51713.9979
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	16275.54665	16317.29331	4.0925021	2000-2009 –> 2010-2019	8.541673	30504.2523	30586.6001
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	16566.57385	16608.32051	4.3049691	2000-2009 –> 2010-2019	8.574255	30593.7943	30676.1421
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	15843.18311	15884.92976	3.7960569	2000-2009 –> 2010-2019	7.863053	29633.4221	29715.7699
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	16255.27061	16297.01727	4.0780962	2000-2009 –> 2010-2019	8.358038	30294.6610	30377.0088
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	17496.78499	17538.72851	4.6516906	2000-2009 –> 2010-2019	9.225392	32183.5771	32266.2777
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	17516.24529	17558.18880	4.6670224	2000-2009 –> 2010-2019	9.257892	32221.7412	32304.4418
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + phosphate	17781.44393	17817.39551	4.8827177	2000-2009 –> 2010-2019	9.746559	32773.2734	32844.1597
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate	15464.97372	15506.72038	3.5544038	2000-2009 –> 2010-2019	7.633367	29269.5563	29351.9041
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate	17359.88071	17401.82422	4.5452446	2000-2009 –> 2010-2019	9.222466	32158.4011	32241.1017
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	15611.91413	15653.66079	3.6464066	2000-2009 –> 2010-2019	7.894356	29614.6745	29697.0223
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	17431.14760	17473.09112	4.6003486	2000-2009 –> 2010-2019	9.313864	32268.2557	32350.9563
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	17519.50586	17561.44937	4.6695962	2000-2009 –> 2010-2019	9.265902	32230.9092	32313.6098
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + nitrate + phosphate_star	17573.99084	17609.94242	4.7144095	2000-2009 –> 2010-2019	9.482402	32466.4637	32537.3500
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	17341.73058	17383.67409	4.5313166	2000-2009 –> 2010-2019	9.245338	32179.3737	32262.0743
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	15654.34909	15696.09575	3.6734166	2000-2009 –> 2010-2019	7.953660	29694.0828	29776.4306
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	17430.56586	17472.50938	4.5998961	2000-2009 –> 2010-2019	9.315174	32269.5398	32352.2405
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + silicate + phosphate_star	17429.63261	17465.58419	4.6007259	2000-2009 –> 2010-2019	9.335597	32287.3071	32358.1933
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	16024.28456	16066.03122	3.9175199	2000-2009 –> 2010-2019	8.195807	30061.7871	30144.1349
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	16088.78451	16130.53117	3.9617116	2000-2009 –> 2010-2019	8.250027	30137.7168	30220.0646
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	16135.50083	16177.24749	3.9940300	2000-2009 –> 2010-2019	8.299901	30204.3787	30286.7265
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	15823.41399	15858.37265	5.6733113	2000-2009 –> 2010-2019	11.054951	32412.1362	32482.4384
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	14207.78878	14248.94382	3.5510526	2000-2009 –> 2010-2019	7.888186	27036.3049	27117.3968
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	14249.46654	14284.74228	3.5805273	2000-2009 –> 2010-2019	8.197080	27353.1916	27422.6989
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	15669.63426	15710.41936	5.4996893	2000-2009 –> 2010-2019	10.527509	31896.9296	31978.9488
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	15694.12302	15734.90812	5.5266267	2000-2009 –> 2010-2019	10.612570	31982.8429	32064.8621
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + aou + phosphate	15875.18189	15910.14055	5.7322135	2000-2009 –> 2010-2019	11.108167	32458.2549	32528.5570
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	14030.45872	14071.61376	3.4338574	2000-2009 –> 2010-2019	7.800796	26889.3825	26970.4744
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	15822.12127	15862.90637	5.6695854	2000-2009 –> 2010-2019	11.036492	32398.1937	32480.2128
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	14200.61385	14241.76889	3.5462340	2000-2009 –> 2010-2019	8.026476	27173.2676	27254.3595
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	15696.61460	15737.39970	5.5293748	2000-2009 –> 2010-2019	10.624593	31995.0707	32077.0899
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + nitrate + phosphate_star	15717.26180	15752.22046	5.5544164	2000-2009 –> 2010-2019	10.762136	32130.4288	32200.7309
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	14714.78113	14755.93617	3.9086523	2000-2009 –> 2010-2019	8.522297	27817.7090	27898.8009
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	15702.91811	15743.70321	5.5363334	2000-2009 –> 2010-2019	10.675833	32047.6170	32129.6362
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	14225.81757	14266.97261	3.5631893	2000-2009 –> 2010-2019	8.051978	27206.9327	27288.0245
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	15858.18344	15898.96854	5.7105261	2000-2009 –> 2010-2019	11.107569	32464.1798	32546.1990
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	14377.68576	14418.84080	3.6670853	2000-2009 –> 2010-2019	8.098241	27301.4248	27382.5167
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	15805.02723	15845.81233	5.6502815	2000-2009 –> 2010-2019	11.069715	32433.1478	32515.1670
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	14594.20185	14635.35689	3.8204681	2000-2009 –> 2010-2019	8.347949	27613.4245	27694.5164
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + aou + silicate + phosphate	14423.08233	14458.35808	3.7001264	2000-2009 –> 2010-2019	8.451175	27654.3112	27723.8186
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	14636.49385	14677.64889	3.8511690	2000-2009 –> 2010-2019	8.395411	27672.1234	27753.2153
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	18480.66702	18517.37090	3.8035477	2000-2009 –> 2010-2019	7.730935	36841.7016	36915.0010
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	24818.73249	24856.52373	5.3016164	2000-2009 –> 2010-2019	10.342916	49408.9928	49484.6184
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	17494.59426	17537.41545	3.2823141	2000-2009 –> 2010-2019	7.202466	35845.4885	35931.0045
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ aou + phosphate + phosphate_star	25127.70744	25159.20014	5.5107982	2000-2009 –> 2010-2019	10.730133	49996.8778	50059.8991
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	17817.76114	17860.58233	3.4444139	2000-2009 –> 2010-2019	7.145450	35789.9598	35875.4758
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	24995.83429	25039.92407	5.4184046	2000-2009 –> 2010-2019	10.556942	49742.7260	49830.9559
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	17992.96232	18035.78351	3.5356162	2000-2009 –> 2010-2019	7.316246	36105.2529	36190.7688
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + aou + silicate + phosphate	25243.04827	25287.13805	5.5876831	2000-2009 –> 2010-2019	10.752477	50031.1939	50119.4238
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	18014.99716	18057.81835	3.5472562	2000-2009 –> 2010-2019	7.331631	36133.8072	36219.3232
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + nitrate + phosphate_star	18501.45669	18538.16057	3.8153611	2000-2009 –> 2010-2019	7.622660	36658.0286	36731.3279
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	17932.06381	17974.88499	3.5036444	2000-2009 –> 2010-2019	7.306846	36083.5465	36169.0625
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	25127.90716	25171.99694	5.5081927	2000-2009 –> 2010-2019	10.595523	49793.7353	49881.9652
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + nitrate + phosphate	25062.14613	25099.93736	5.4646613	2000-2009 –> 2010-2019	10.665323	49904.3075	49979.9331
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	17477.89930	17520.72048	3.2741503	2000-2009 –> 2010-2019	7.203513	35844.2452	35929.7611
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	25033.34083	25077.43061	5.4437531	2000-2009 –> 2010-2019	10.529435	49696.5461	49784.7760
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	17661.72048	17704.54166	3.3651684	2000-2009 –> 2010-2019	7.393713	36192.1909	36277.7068
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	25175.62909	25219.71886	5.5410006	2000-2009 –> 2010-2019	10.705684	49963.6017	50051.8316
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	17689.79391	17732.61510	3.3792898	2000-2009 –> 2010-2019	7.413006	36228.7118	36314.2278
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	25176.16376	25220.25354	5.5413692	2000-2009 –> 2010-2019	10.704699	49962.0148	50050.2447
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	25174.49248	25212.28372	5.5415959	2000-2009 –> 2010-2019	10.706371	49960.6078	50036.2335
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	21219.99254	21258.36947	2.6504782	2000-2009 –> 2010-2019	5.832352	43167.5378	43244.0526
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	21171.92464	21216.69772	2.6355426	2000-2009 –> 2010-2019	5.646992	42652.7850	42742.0524
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	31763.37686	31809.36669	4.9173257	2000-2009 –> 2010-2019	9.670583	62924.8369	63016.7619
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	31969.13693	32015.12676	5.0142455	2000-2009 –> 2010-2019	9.875603	63365.8182	63457.7432
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + aou + silicate + phosphate	31764.42613	31803.84598	4.9187483	2000-2009 –> 2010-2019	9.673068	62926.2231	63005.0159
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + aou + silicate + phosphate_star	31973.67169	32013.09154	5.0173547	2000-2009 –> 2010-2019	9.897905	63409.5674	63488.3603
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	31897.63079	31943.62062	4.9803490	2000-2009 –> 2010-2019	9.800357	63204.9604	63296.8854
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	31903.68159	31943.10144	4.9841539	2000-2009 –> 2010-2019	9.826022	63256.3434	63335.1362
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	21163.78236	21208.55544	2.6331217	2000-2009 –> 2010-2019	5.798112	43068.0291	43157.2964
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	21562.42513	21607.19822	2.7543008	2000-2009 –> 2010-2019	6.015747	43722.2700	43811.5373
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + aou + silicate + phosphate	31763.58471	31809.57454	4.9174227	2000-2009 –> 2010-2019	9.669390	62922.2061	63014.1311
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	31940.30886	31986.29869	5.0005523	2000-2009 –> 2010-2019	9.811504	63227.9673	63319.8923
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	21590.07666	21634.84974	2.7629103	2000-2009 –> 2010-2019	6.032196	43770.3607	43859.6281
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	22054.74603	22099.51911	2.9116802	2000-2009 –> 2010-2019	6.153266	44162.5874	44251.8547
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	32031.09280	32077.08263	5.0438012	2000-2009 –> 2010-2019	9.896568	63409.2719	63501.1969
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	31895.66953	31941.65936	4.9794226	2000-2009 –> 2010-2019	9.779178	63158.9554	63250.8804
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ temp + aou + nitrate + phosphate	21550.00328	21588.38021	2.7510629	2000-2009 –> 2010-2019	5.978938	43619.7552	43696.2700
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	21550.63127	21595.40436	2.7506369	2000-2009 –> 2010-2019	5.870347	43332.1941	43421.4614
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	21910.31252	21955.08560	2.8645996	2000-2009 –> 2010-2019	6.102957	44009.6704	44098.9377
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	21930.96116	21975.73424	2.8712834	2000-2009 –> 2010-2019	6.113399	44040.1938	44129.4611
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	33845.35583	33885.44088	4.2787063	2000-2009 –> 2010-2019	8.643229	66594.0761	66674.0072
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	33842.50450	33889.27039	4.2769442	2000-2009 –> 2010-2019	8.584594	66444.7732	66538.0261
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	29322.74883	29369.25808	3.1976634	2000-2009 –> 2010-2019	6.264292	57017.4867	57110.2196
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	33535.38943	33582.15532	4.1668630	2000-2009 –> 2010-2019	8.405941	65956.0422	66049.2951
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + aou + silicate + phosphate	33540.32546	33580.41051	4.1693176	2000-2009 –> 2010-2019	8.419418	65988.3687	66068.2998
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	33782.62223	33829.38812	4.2552544	2000-2009 –> 2010-2019	8.589281	66453.9825	66547.2354
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	33783.20650	33823.29154	4.2561881	2000-2009 –> 2010-2019	8.633429	66564.8579	66644.7890
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	28896.44535	28942.95460	3.0798282	2000-2009 –> 2010-2019	6.265283	57005.5617	57098.2945
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	29157.32312	29203.83237	3.1514118	2000-2009 –> 2010-2019	6.405939	57500.2644	57592.9973
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate	33517.11849	33563.88438	4.1604040	2000-2009 –> 2010-2019	8.355279	65819.1315	65912.3844
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	33730.73221	33777.49810	4.2365483	2000-2009 –> 2010-2019	8.483595	66172.6409	66265.8938
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	29176.68908	29223.19833	3.1567916	2000-2009 –> 2010-2019	6.417365	57539.8644	57632.5973
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	28784.37818	28830.88743	3.0495789	2000-2009 –> 2010-2019	6.100708	56423.8849	56516.6178
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	33673.49644	33720.26233	4.2160106	2000-2009 –> 2010-2019	8.467790	66128.0094	66221.2624
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	33677.51770	33724.28359	4.2174503	2000-2009 –> 2010-2019	8.455152	66094.4942	66187.7471
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + aou + nitrate + phosphate	28898.22796	28938.09303	3.0808544	2000-2009 –> 2010-2019	6.331044	57224.6336	57304.1189
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	28575.04104	28621.55029	2.9938681	2000-2009 –> 2010-2019	6.126370	56501.4382	56594.1710
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	28856.50848	28903.01772	3.0690141	2000-2009 –> 2010-2019	6.289652	57085.3544	57178.0872
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	28858.00240	28904.51165	3.0694180	2000-2009 –> 2010-2019	6.291171	57090.6221	57183.3549
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	28856.46839	28896.33346	3.0695440	2000-2009 –> 2010-2019	6.292250	57090.3107	57169.7960
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	38445.20447	38486.19649	4.0004420	2000-2009 –> 2010-2019	7.670224	72668.5519	72750.0252
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	38973.97974	39021.80377	4.1572575	2000-2009 –> 2010-2019	7.798112	73099.7566	73194.8088
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	32817.26501	32864.56950	3.1897244	2000-2009 –> 2010-2019	6.102882	64978.4996	65073.2178
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	38211.46218	38259.28621	3.9321934	2000-2009 –> 2010-2019	7.770687	73002.3438	73097.3960
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	32990.70544	33038.00992	3.2335149	2000-2009 –> 2010-2019	6.147023	65153.4945	65248.2127
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + silicate + phosphate_star	33256.28292	33296.82963	3.3022555	2000-2009 –> 2010-2019	6.342383	65966.7092	66047.8963
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	32943.99984	32991.30433	3.2216638	2000-2009 –> 2010-2019	6.130840	65087.5653	65182.2834
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	33183.38762	33230.69211	3.2828690	2000-2009 –> 2010-2019	6.286897	65741.4784	65836.1966
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	37999.44519	38047.26921	3.8718085	2000-2009 –> 2010-2019	7.455853	71927.3518	72022.4039
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	33243.02044	33290.32493	3.2982956	2000-2009 –> 2010-2019	6.316140	65860.4002	65955.1184
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	38097.24697	38145.07099	3.8995475	2000-2009 –> 2010-2019	7.516467	72140.0376	72235.0897
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate	33245.75140	33293.05589	3.2990038	2000-2009 –> 2010-2019	6.304641	65810.7630	65905.4811
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	32932.59040	32979.89488	3.2187754	2000-2009 –> 2010-2019	6.181824	65313.2227	65407.9409
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	38062.59124	38110.41527	3.8896956	2000-2009 –> 2010-2019	7.501032	72085.9457	72180.9979
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + silicate + phosphate	33282.57372	33323.12042	3.3090880	2000-2009 –> 2010-2019	6.345521	65977.2908	66058.4778
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	33213.98625	33261.29074	3.2907756	2000-2009 –> 2010-2019	6.301130	65799.2586	65893.9768
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	38083.03050	38130.85453	3.8955030	2000-2009 –> 2010-2019	7.511014	72120.9202	72215.9723
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + silicate + phosphate_star	38374.35443	38415.34645	3.9798070	2000-2009 –> 2010-2019	7.604961	72443.7580	72525.2313
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	38518.98995	38566.81397	4.0214586	2000-2009 –> 2010-2019	7.737200	72900.9336	72995.9858
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ temp + aou + silicate + phosphate	38519.83536	38560.82739	4.0222939	2000-2009 –> 2010-2019	7.745320	72924.4213	73005.8946
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	11508.36973	11549.09593	2.4444299	2000-2009 –> 2010-2019	4.677545	22634.0509	22715.5454
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	11549.13308	11589.85927	2.4645612	2000-2009 –> 2010-2019	4.705603	22692.5287	22774.0233
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + aou + phosphate	11837.08033	11871.98850	2.6126193	2000-2009 –> 2010-2019	4.894405	23068.5641	23138.4165
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + aou + silicate + phosphate	11592.45381	11633.18001	2.4861374	2000-2009 –> 2010-2019	4.767515	22825.0449	22906.5394
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + aou + silicate + phosphate	21251.47382	21295.20436	3.9079788	2000-2009 –> 2010-2019	6.998440	37649.4438	37735.7055
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	11632.77237	11673.49857	2.5063879	2000-2009 –> 2010-2019	4.795171	22881.5669	22963.0614
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	21230.31766	21274.04820	3.8971636	2000-2009 –> 2010-2019	6.989128	37631.4151	37717.6769
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	11574.44124	11615.16744	2.4771433	2000-2009 –> 2010-2019	4.723791	22730.3295	22811.8240
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	21228.12501	21271.85555	3.8960444	2000-2009 –> 2010-2019	6.996773	37647.4280	37733.6897
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	11649.68979	11690.41598	2.5149339	2000-2009 –> 2010-2019	4.805945	22903.3491	22984.8436
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	21229.75761	21273.48816	3.8968777	2000-2009 –> 2010-2019	6.989247	37631.6959	37717.9577
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	11501.35831	11542.08451	2.4409838	2000-2009 –> 2010-2019	4.838171	22981.5317	23063.0263
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	21256.17809	21299.90863	3.9103877	2000-2009 –> 2010-2019	7.000657	37653.7476	37740.0093
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	11552.02205	11592.74825	2.4659943	2000-2009 –> 2010-2019	4.871868	23050.2815	23131.7761
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	21237.67177	21281.40231	3.9009197	2000-2009 –> 2010-2019	6.992876	37638.7525	37725.0143
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + aou + silicate + phosphate_star	21243.02935	21280.51267	3.9046812	2000-2009 –> 2010-2019	6.996664	37642.1644	37716.1031
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	11586.28395	11627.01014	2.4830529	2000-2009 –> 2010-2019	4.894867	23096.8731	23178.3676
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	21237.88989	21281.62043	3.9010311	2000-2009 –> 2010-2019	6.993376	37639.7789	37726.0407
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	21236.09232	21273.57564	3.9011346	2000-2009 –> 2010-2019	7.001867	37653.4018	37727.3405
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	21242.05814	21279.54146	3.9041844	2000-2009 –> 2010-2019	6.996568	37642.0259	37715.9646
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	27294.76927	27340.68701	3.3057541	2000-2009 –> 2010-2019	6.049702	48697.5535	48788.2032
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	13697.92152	13740.35855	2.0904264	2000-2009 –> 2010-2019	4.107884	26389.9260	26474.3819
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	27321.20879	27367.12653	3.3141414	2000-2009 –> 2010-2019	6.059210	48727.5899	48818.2395
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou + nitrate	13091.54367	13127.91827	1.9005255	2000-2009 –> 2010-2019	3.852592	25584.5748	25656.9655
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	13093.29600	13135.73303	1.9004513	2000-2009 –> 2010-2019	3.847632	25573.3456	25657.8015
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	13093.45109	13135.88812	1.9004978	2000-2009 –> 2010-2019	3.849053	25577.7195	25662.1754
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou + nitrate + silicate	13076.22168	13118.65871	1.8953449	2000-2009 –> 2010-2019	3.815854	25473.8213	25558.2772
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou + phosphate	13721.05078	13757.42538	2.0987206	2000-2009 –> 2010-2019	4.152036	26516.3726	26588.7634
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate + phosphate	13572.57758	13615.01461	2.0495421	2000-2009 –> 2010-2019	4.096647	26351.7903	26436.2462
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate + phosphate	27275.23964	27321.15739	3.2995724	2000-2009 –> 2010-2019	6.048050	48692.5531	48783.2027
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	13590.36587	13632.80290	2.0552952	2000-2009 –> 2010-2019	4.107931	26385.7071	26470.1630
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	27266.15943	27312.07717	3.2967022	2000-2009 –> 2010-2019	6.056281	48718.9782	48809.6278
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	13137.57753	13180.01456	1.9137588	2000-2009 –> 2010-2019	3.871402	25649.6611	25734.1170
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	13606.42644	13648.86347	2.0605033	2000-2009 –> 2010-2019	4.115354	26408.2161	26492.6720
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	27266.04559	27311.96333	3.2966662	2000-2009 –> 2010-2019	6.057919	48724.2035	48814.8532
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	27257.88895	27303.80669	3.2940901	2000-2009 –> 2010-2019	6.042244	48674.1672	48764.8168
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	27254.68271	27300.60045	3.2930780	2000-2009 –> 2010-2019	6.052784	48707.9053	48798.5550
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + aou + silicate + phosphate	27280.33019	27319.68826	3.3018154	2000-2009 –> 2010-2019	6.050962	48697.7871	48775.4867
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + aou + silicate + phosphate_star	27273.31753	27312.67560	3.2995970	2000-2009 –> 2010-2019	6.059753	48725.9773	48803.6769
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	27254.69106	27300.60880	3.2930807	2000-2009 –> 2010-2019	6.054057	48711.9700	48802.6196
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	46540.46390	46582.93845	3.4341433	2000-2009 –> 2010-2019	6.285449	83247.3241	83331.2848
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	46540.79655	46590.35020	3.4338169	2000-2009 –> 2010-2019	6.272729	83184.8573	83282.8115
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ aou + phosphate + phosphate_star	46697.28183	46732.67729	3.4653794	2000-2009 –> 2010-2019	6.320612	83422.6140	83492.5812
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	46696.67724	46739.15179	3.4648648	2000-2009 –> 2010-2019	6.312848	83386.1933	83470.1540
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	46542.44437	46591.99802	3.4341395	2000-2009 –> 2010-2019	6.285418	83251.1608	83349.1150
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + aou + phosphate + phosphate_star	46684.77655	46727.25110	3.4625148	2000-2009 –> 2010-2019	6.317594	83411.3090	83495.2696
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + aou + nitrate	17705.95762	17744.41794	1.7342206	2000-2009 –> 2010-2019	3.460006	30495.8511	30570.8335
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + aou + nitrate + silicate	17579.27697	17624.14734	1.7095579	2000-2009 –> 2010-2019	3.388501	30192.1964	30279.6760
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + aou + phosphate	18006.65898	18045.11930	1.7932486	2000-2009 –> 2010-2019	3.620038	31166.4843	31241.4668
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + aou + phosphate_star	18009.22715	18047.68747	1.7937613	2000-2009 –> 2010-2019	3.624702	31183.8178	31258.8002
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate + phosphate	18002.04340	18046.91378	1.7919286	2000-2009 –> 2010-2019	3.588770	31056.3590	31143.8385
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate + phosphate	46780.81669	46830.37033	3.4811288	2000-2009 –> 2010-2019	6.338918	83523.4648	83621.4189
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	18005.55862	18050.42900	1.7926299	2000-2009 –> 2010-2019	3.595302	31080.9467	31168.4262
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	18226.99414	18271.86451	1.8373632	2000-2009 –> 2010-2019	3.678817	31440.8225	31528.3021
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + phosphate + phosphate_star	18015.82616	18054.28648	1.7950794	2000-2009 –> 2010-2019	3.628039	31197.5821	31272.5646
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	18011.93874	18056.80911	1.7939034	2000-2009 –> 2010-2019	3.598108	31092.8543	31180.3339
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + silicate + phosphate_star	18283.25052	18321.71084	1.8493163	2000-2009 –> 2010-2019	3.707881	31555.2356	31630.2181
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	46629.10884	46678.66248	3.4511495	2000-2009 –> 2010-2019	6.303230	83342.0103	83439.9645
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	46743.38793	46792.94157	3.4737083	2000-2009 –> 2010-2019	6.353840	83601.8867	83699.8409
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	46727.98276	46777.53640	3.4706587	2000-2009 –> 2010-2019	6.348936	83576.8999	83674.8541
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	35851.85667	35893.11175	2.9591257	2000-2009 –> 2010-2019	5.687428	58713.9915	58793.9989
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ aou + phosphate + phosphate_star	35925.96892	35960.34815	2.9749022	2000-2009 –> 2010-2019	5.704864	58791.8378	58858.5106
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	35902.61283	35943.86790	2.9696371	2000-2009 –> 2010-2019	5.698621	58767.1020	58847.1094
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	35868.27730	35916.40822	2.9621083	2000-2009 –> 2010-2019	5.702661	58774.6667	58868.0086
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	35866.89631	35915.02723	2.9618226	2000-2009 –> 2010-2019	5.716464	58821.6488	58914.9908
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + aou + nitrate + silicate	9290.31028	9331.31264	1.4554553	2000-2009 –> 2010-2019	3.126998	16737.4183	16817.3522
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + aou + silicate + phosphate	9477.66197	9518.66433	1.5091657	2000-2009 –> 2010-2019	3.087607	16704.3592	16784.2930
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + aou + silicate + phosphate	35866.07395	35914.20488	2.9616524	2000-2009 –> 2010-2019	5.674405	58676.2969	58769.6388
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	9482.91700	9523.91936	1.5107005	2000-2009 –> 2010-2019	3.093563	16720.3710	16800.3049
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	35986.15693	36034.28786	2.9866028	2000-2009 –> 2010-2019	5.722168	58875.3643	58968.7063
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	9432.87632	9473.87868	1.4961489	2000-2009 –> 2010-2019	3.078990	16670.2792	16750.2131
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	35955.43081	36003.56173	2.9801987	2000-2009 –> 2010-2019	5.709905	58824.4164	58917.7583
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	9482.28253	9523.28489	1.5105151	2000-2009 –> 2010-2019	3.093377	16719.7362	16799.6700
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	35987.63666	36035.76758	2.9869116	2000-2009 –> 2010-2019	5.722113	58875.5928	58968.9347
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + silicate + phosphate_star	9487.80959	9522.95447	1.5127159	2000-2009 –> 2010-2019	3.095579	16723.2671	16791.7819
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	36048.29575	36096.42667	2.9995962	2000-2009 –> 2010-2019	5.722129	58892.4612	58985.8032
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + aou + silicate + phosphate	9886.01203	9921.15692	1.6338324	2000-2009 –> 2010-2019	3.237043	17170.5947	17239.1094
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + aou + silicate + phosphate_star	9894.26552	9929.41040	1.6364427	2000-2009 –> 2010-2019	3.245576	17193.0292	17261.5439
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	9884.18414	9919.32902	1.6332548	2000-2009 –> 2010-2019	3.241924	17181.8389	17250.3536
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	9892.19372	9927.33860	1.6357871	2000-2009 –> 2010-2019	3.244897	17190.9014	17259.4162
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	18973.22438	19019.59163	1.3299060	2000-2009 –> 2010-2019	2.771644	36449.5900	36541.4248
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	37477.96473	37526.62753	2.7369813	2000-2009 –> 2010-2019	5.271126	67006.9208	67102.8028
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + aou + phosphate + phosphate_star	37542.81709	37584.52806	2.7488546	2000-2009 –> 2010-2019	5.288130	67095.1855	67177.3701
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + aou + silicate + phosphate	37512.45862	37554.16960	2.7434564	2000-2009 –> 2010-2019	5.279460	67048.6285	67130.8130
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + nitrate + silicate + phosphate	19063.02637	19102.76973	1.3409246	2000-2009 –> 2010-2019	2.803834	36680.0179	36758.7334
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	19064.85062	19111.21787	1.3409034	2000-2009 –> 2010-2019	2.797721	36643.0166	36734.8513
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	37592.71074	37641.37354	2.7573923	2000-2009 –> 2010-2019	5.326216	67292.4625	67388.3445
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	37636.16615	37684.82895	2.7651618	2000-2009 –> 2010-2019	5.339630	67363.4978	67459.3798
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + aou + nitrate + silicate	18062.38697	18108.75422	1.2253697	2000-2009 –> 2010-2019	2.747791	36071.5195	36163.3542
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + aou + silicate + phosphate	17650.93462	17697.30187	1.1808818	2000-2009 –> 2010-2019	2.556017	34664.5456	34756.3803
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + aou + silicate + phosphate	37497.17657	37545.83937	2.7403882	2000-2009 –> 2010-2019	5.257680	66942.2845	67038.1665
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	17669.99984	17716.36709	1.1829070	2000-2009 –> 2010-2019	2.559938	34697.0893	34788.9240
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	37599.04335	37647.70615	2.7585232	2000-2009 –> 2010-2019	5.295811	67143.5747	67239.4568
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate	18515.19677	18561.56402	1.2762689	2000-2009 –> 2010-2019	2.715247	35972.8164	36064.6512
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	17628.64675	17675.01400	1.1785186	2000-2009 –> 2010-2019	2.565869	34728.7782	34820.6130
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	37675.58961	37724.25241	2.7722294	2000-2009 –> 2010-2019	5.326795	67305.4008	67401.2828
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	17666.44978	17712.81703	1.1825296	2000-2009 –> 2010-2019	2.559293	34691.6358	34783.4705
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	37582.76522	37631.42802	2.7556172	2000-2009 –> 2010-2019	5.288876	67107.3267	67203.2087
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + silicate + phosphate_star	17843.69286	17883.43622	1.2017348	2000-2009 –> 2010-2019	2.592203	34963.7917	35042.5072
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + silicate + phosphate_star	37705.82597	37747.53694	2.7780219	2000-2009 –> 2010-2019	5.340581	67372.8927	67455.0773
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	11304.11049	11339.41569	2.0291885	2000-2009 –> 2010-2019	4.639251	19554.5862	19622.6336
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	11223.56757	11264.75696	1.9978890	2000-2009 –> 2010-2019	4.605970	19473.4114	19552.8000
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	11246.09273	11287.28212	2.0063821	2000-2009 –> 2010-2019	4.613945	19495.2485	19574.6371
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	31795.32980	31848.65628	0.6963125	2000-2009 –> 2010-2019	1.350588	54011.2932	54115.8603
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	34922.76295	34976.08944	0.7726380	2000-2009 –> 2010-2019	1.541949	60753.8398	60858.4069
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + aou + phosphate + phosphate_star	11377.89427	11413.19947	2.0575814	2000-2009 –> 2010-2019	4.666780	19627.2224	19695.2698
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	11368.01035	11409.19975	2.0529816	2000-2009 –> 2010-2019	4.662045	19619.1589	19698.5476
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + nitrate + silicate + phosphate	32249.18898	32294.89739	0.7069497	2000-2009 –> 2010-2019	1.453830	57417.8068	57507.4358
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	32071.05612	32124.38261	0.7027271	2000-2009 –> 2010-2019	1.362640	54478.5056	54583.0727
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	11228.53031	11269.71971	1.9997571	2000-2009 –> 2010-2019	4.657930	19544.2750	19623.6636
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	11360.97678	11402.16618	2.0502640	2000-2009 –> 2010-2019	4.658497	19611.0231	19690.4117
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + temp + aou + silicate + phosphate	34854.17621	34907.50270	0.7708775	2000-2009 –> 2010-2019	1.601780	62404.2610	62508.8281
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	34866.36203	34919.68852	0.7711900	2000-2009 –> 2010-2019	1.602395	62424.5726	62529.1397
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate	32060.89791	32114.22440	0.7024898	2000-2009 –> 2010-2019	1.441528	56995.9115	57100.4786
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	34863.98913	34917.31562	0.7711292	2000-2009 –> 2010-2019	1.602296	62421.1771	62525.7443
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star	11260.30421	11301.49361	2.0117591	2000-2009 –> 2010-2019	4.621747	19512.6800	19592.0686
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	36017.98625	36071.31274	0.8013001	2000-2009 –> 2010-2019	1.638793	63744.4142	63848.9813
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	11231.61396	11272.80336	2.0009188	2000-2009 –> 2010-2019	4.648883	19534.0297	19613.4183
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	11374.26922	11415.45862	2.0554029	2000-2009 –> 2010-2019	4.664953	19626.0636	19705.4523
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	36017.47367	36070.80016	0.8012864	2000-2009 –> 2010-2019	1.640030	63777.1906	63881.7577
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	20732.88002	20780.08362	1.2630738	2000-2009 –> 2010-2019	2.190018	33637.7296	33730.2663
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	20404.07985	20451.28346	1.2303812	2000-2009 –> 2010-2019	2.159145	33327.7509	33420.2876
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	4498.01762	4532.70677	2.0510328	2000-2009 –> 2010-2019	4.699276	7775.8445	7842.2088
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	4481.76681	4516.45596	2.0352071	2000-2009 –> 2010-2019	4.669712	7752.4991	7818.8635
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + aou + silicate + phosphate_star	4501.98944	4531.72299	2.0568792	2000-2009 –> 2010-2019	4.707437	7779.0081	7835.8918
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	4529.52786	4564.21701	2.0820701	2000-2009 –> 2010-2019	4.715508	7799.7076	7866.0720
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	4530.78745	4560.52100	2.0853075	2000-2009 –> 2010-2019	4.720217	7799.7293	7856.6130
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	20381.35774	20428.56134	1.2281534	2000-2009 –> 2010-2019	2.154205	33276.9575	33369.4942
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	20463.78985	20510.99346	1.2362546	2000-2009 –> 2010-2019	2.167270	33410.7022	33503.2389
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + aou + silicate + phosphate	4506.70614	4541.39529	2.0595444	2000-2009 –> 2010-2019	4.693072	7776.9326	7843.2970
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	4480.75815	4515.44730	2.0342289	2000-2009 –> 2010-2019	4.683941	7759.3415	7825.7058
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	20379.26391	20426.46751	1.2279483	2000-2009 –> 2010-2019	2.155069	33285.9390	33378.4758
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate	20354.16131	20401.36491	1.2254923	2000-2009 –> 2010-2019	2.191741	33657.5948	33750.1316
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate	4506.09510	4540.78424	2.0589446	2000-2009 –> 2010-2019	4.708945	7784.8270	7851.1914
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	4500.47447	4535.16362	2.0534360	2000-2009 –> 2010-2019	4.687799	7771.1334	7837.4978
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ temp + aou + nitrate + phosphate	20713.00379	20753.46402	1.2612743	2000-2009 –> 2010-2019	2.187867	33612.2085	33691.5257
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	20714.64697	20761.85057	1.2612384	2000-2009 –> 2010-2019	2.185142	33587.9564	33680.4931
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	4478.37010	4513.05925	2.0319148	2000-2009 –> 2010-2019	4.717310	7775.1996	7841.5640
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ temp + nitrate + phosphate + phosphate_star	20715.76399	20756.22422	1.2615520	2000-2009 –> 2010-2019	2.189271	33626.6296	33705.9469
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	20717.36442	20764.56802	1.2615118	2000-2009 –> 2010-2019	2.186601	33602.9846	33695.5214
obs	Indo-Pacific	(28.4, Inf]	cstar_tref ~ sal + temp + aou + phosphate	-Inf	-Inf	0.0000000	2000-2009 –> 2010-2019	1.656889	-Inf	-Inf
obs	Indo-Pacific	(28.4, Inf]	cstar_tref ~ sal + temp + nitrate + phosphate_star	-Inf	-Inf	0.0000000	2000-2009 –> 2010-2019	1.652663	-Inf	-Inf
obs	Indo-Pacific	(28.4, Inf]	cstar_tref ~ sal + temp + silicate + phosphate	-Inf	-Inf	0.0000000	2000-2009 –> 2010-2019	1.519547	-Inf	-Inf
obs	Indo-Pacific	(28.4, Inf]	cstar_tref ~ sal + temp + silicate + phosphate_star	-Inf	-Inf	0.0000000	2000-2009 –> 2010-2019	1.538645	-Inf	-Inf
obs	Indo-Pacific	(28.4, Inf]	cstar_tref ~ temp + aou + silicate + phosphate	-Inf	-Inf	0.0000000	2000-2009 –> 2010-2019	1.581647	-Inf	-Inf
obs	Indo-Pacific	(28.4, Inf]	cstar_tref ~ temp + aou + silicate + phosphate_star	-Inf	-Inf	0.0000000	2000-2009 –> 2010-2019	1.545706	-Inf	-Inf
obs	Indo-Pacific	(28.4, Inf]	cstar_tref ~ temp + nitrate + phosphate + phosphate_star	-Inf	-Inf	0.0000000	2000-2009 –> 2010-2019	1.853069	-Inf	-Inf
obs	Indo-Pacific	(28.4, Inf]	cstar_tref ~ temp + nitrate + silicate + phosphate	-Inf	-Inf	0.0000000	2000-2009 –> 2010-2019	1.567579	-Inf	-Inf
obs	Indo-Pacific	(28.4, Inf]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	-Inf	-Inf	0.0000000	2000-2009 –> 2010-2019	1.674364	-Inf	-Inf
obs	Indo-Pacific	(28.4, Inf]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	-Inf	-Inf	0.0000000	2000-2009 –> 2010-2019	1.565887	-Inf	-Inf

7.3 Target variable coefficients

A data frame to map the target variable is prepared.

# create table with two era belonging to one eras
eras_forward <- GLODAP_glanced %>%
  arrange(era) %>% 
  group_by(basin, gamma_slab, model) %>% 
  mutate(eras = paste(era, lead(era), sep = " --> ")) %>% 
  ungroup() %>% 
  select(era, eras) %>% 
  unique()

eras_backward <- GLODAP_glanced %>%
  arrange(era) %>% 
  group_by(basin, gamma_slab, model) %>% 
  mutate(eras = paste(lag(era), era, sep = " --> ")) %>% 
  ungroup() %>% 
  select(era, eras) %>% 
  unique()

eras_era <- full_join(eras_backward, eras_forward) %>% 
  filter(str_detect(eras, "NA") == FALSE)

# extend best model selection from eras to era
lm_best <- full_join(
  lm_best_eras %>% select(data_source, basin, gamma_slab, model, eras),
  eras_era) %>% 
  drop_na() %>% 
  arrange(era, basin, gamma_slab, data_source)

lm_best <- left_join(
  lm_best,
  GLODAP_tidy %>% select(data_source, basin, gamma_slab, era, model, term, estimate))

rm(eras_era, eras_forward, eras_backward)

7.4 Cant coeffcients

A data frame of coefficient offsets is prepared to facilitate the direct mapping of Cant.

# subtract coefficients of adjacent era  
lm_best_cant <- lm_best %>%
  arrange(era) %>%
  group_by(data_source, basin, gamma_slab, eras, model, term) %>%
  mutate(delta_coeff = estimate - lag(estimate)) %>%
  ungroup() %>%
  arrange(basin, gamma_slab, model, term, eras) %>%
  drop_na() %>%
  select(-c(era,estimate))

# pivot to wide format
lm_best_cant <- lm_best_cant %>%
  pivot_wider(values_from = delta_coeff,
              names_from = term,
              names_prefix = "delta_coeff_",
              values_fill = 0)

7.5 Write files

# create table of target varaible coefficients in wide format
lm_best_target <- lm_best %>% 
  pivot_wider(names_from = "term",
              names_prefix = "coeff_",
              values_from = "estimate",
              values_fill = 0
              )


lm_best_target %>%
  write_csv(paste(path_version_data,
                  "lm_best_target.csv",
                  sep = ""))

lm_best_cant %>%
  write_csv(paste(path_version_data,
                  "lm_best_cant.csv",
                  sep = ""))

8 Model diagnotics

8.1 Selection criterion vs predictors

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

8.1.1 All models

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

[[1]]

Past versions of eMLR_predictor_model_quality_all-1.png

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

[[2]]

Past versions of eMLR_predictor_model_quality_all-2.png

Version	Author	Date
33ffcab	jens-daniel-mueller	2021-06-10
7e1f407	jens-daniel-mueller	2021-06-10
2cbe18c	jens-daniel-mueller	2021-06-10
69c79d0	jens-daniel-mueller	2021-06-08
1772903	jens-daniel-mueller	2021-06-07
594ed9a	jens-daniel-mueller	2021-06-04
db7df0e	jens-daniel-mueller	2021-06-04
2edc791	jens-daniel-mueller	2021-06-04
207339d	jens-daniel-mueller	2021-06-03
315710b	jens-daniel-mueller	2021-06-03
be90356	jens-daniel-mueller	2021-06-02
d37a85d	jens-daniel-mueller	2021-05-31
4b7a5ee	jens-daniel-mueller	2021-05-28
25bd183	jens-daniel-mueller	2021-05-26
62bd574	jens-daniel-mueller	2021-05-20
7c56c39	jens-daniel-mueller	2021-05-19
52e7583	jens-daniel-mueller	2021-05-12
969e631	jens-daniel-mueller	2021-05-12
d2a83bc	jens-daniel-mueller	2021-04-16
c0a47df	jens-daniel-mueller	2021-04-16
50290e8	jens-daniel-mueller	2021-04-16
a00ec94	jens-daniel-mueller	2021-04-16
b6fe355	jens-daniel-mueller	2021-04-16
81b7c6d	jens-daniel-mueller	2021-04-16
ddec5b7	jens-daniel-mueller	2021-04-15
29edae5	jens-daniel-mueller	2021-04-14
099d566	jens-daniel-mueller	2021-04-14
bb44686	jens-daniel-mueller	2021-04-14
bf40480	jens-daniel-mueller	2021-04-13
9f31fe3	jens-daniel-mueller	2021-04-13
338dd3c	jens-daniel-mueller	2021-04-09
a79ca2c	jens-daniel-mueller	2021-04-09
d3c3fde	jens-daniel-mueller	2021-03-26
594d503	jens-daniel-mueller	2021-03-26
2eb881a	jens-daniel-mueller	2021-03-26

8.1.2 Best models

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

[[1]]

Past versions of eMLR_predictor_model_quality_best-1.png

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

[[2]]

Past versions of eMLR_predictor_model_quality_best-2.png

Version	Author	Date
33ffcab	jens-daniel-mueller	2021-06-10
7e1f407	jens-daniel-mueller	2021-06-10
2cbe18c	jens-daniel-mueller	2021-06-10
69c79d0	jens-daniel-mueller	2021-06-08
1772903	jens-daniel-mueller	2021-06-07
594ed9a	jens-daniel-mueller	2021-06-04
db7df0e	jens-daniel-mueller	2021-06-04
2edc791	jens-daniel-mueller	2021-06-04
207339d	jens-daniel-mueller	2021-06-03
315710b	jens-daniel-mueller	2021-06-03
be90356	jens-daniel-mueller	2021-06-02
d37a85d	jens-daniel-mueller	2021-05-31
4b7a5ee	jens-daniel-mueller	2021-05-28
25bd183	jens-daniel-mueller	2021-05-26
62bd574	jens-daniel-mueller	2021-05-20
7c56c39	jens-daniel-mueller	2021-05-19
52e7583	jens-daniel-mueller	2021-05-12
969e631	jens-daniel-mueller	2021-05-12
d2a83bc	jens-daniel-mueller	2021-04-16
c0a47df	jens-daniel-mueller	2021-04-16
50290e8	jens-daniel-mueller	2021-04-16
a00ec94	jens-daniel-mueller	2021-04-16
b6fe355	jens-daniel-mueller	2021-04-16
81b7c6d	jens-daniel-mueller	2021-04-16
ddec5b7	jens-daniel-mueller	2021-04-15
29edae5	jens-daniel-mueller	2021-04-14
099d566	jens-daniel-mueller	2021-04-14
bb44686	jens-daniel-mueller	2021-04-14
bf40480	jens-daniel-mueller	2021-04-13
9f31fe3	jens-daniel-mueller	2021-04-13
338dd3c	jens-daniel-mueller	2021-04-09
a79ca2c	jens-daniel-mueller	2021-04-09
d3c3fde	jens-daniel-mueller	2021-03-26
594d503	jens-daniel-mueller	2021-03-26

8.2 RMSE correlation between eras

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

8.2.1 All models

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

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

Past versions of rmse_correlation_across_eras_all_models-1.png

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

rm(max_rmse)

8.2.2 Best models

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

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

Past versions of rmse_correlation_across_eras_best_models-1.png

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

rm(max_rmse)

8.3 Predictor counts

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

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

# print table
lm_all_stats %>%
  gt(rowname_col = "gamma_slab",
     groupname_col = c("basin", "data_source")) %>% 
  summary_rows(
    groups = TRUE,
    fns = list(total = "sum")
  )

	aou	nitrate	phosphate	phosphate_star	sal	silicate	temp
Atlantic - mod
(-Inf,26]	6	6	6	8	10	8	6
(26,26.5]	9	7	7	7	7	6	7
(26.5,26.75]	6	7	7	7	6	6	10
(26.75,27]	8	5	8	6	5	10	6
(27,27.25]	9	7	6	4	6	10	6
(27.25,27.5]	6	8	4	7	9	8	7
(27.5,27.75]	7	5	8	5	10	7	6
(27.75,27.85]	5	4	4	6	10	10	6
(27.85,27.95]	7	4	4	5	10	9	7
(27.95,28.05]	4	2	4	8	10	5	10
(28.05,28.1]	6	2	6	6	8	5	10
(28.1,28.15]	5	2	6	7	8	6	10
(28.15,28.2]	6	10	2	6	5	5	9
(28.2, Inf]	9	10	3	4	10	3	5
total	93.00	79.00	75.00	86.00	114.00	98.00	105.00
Indo-Pacific - mod
(-Inf,26]	8	5	7	7	10	5	7
(26,26.5]	8	6	8	7	6	8	7
(26.5,26.75]	7	5	7	7	4	10	8
(26.75,27]	6	10	6	8	5	5	8
(27,27.25]	7	10	7	7	4	5	8
(27.25,27.5]	6	10	6	6	5	6	9
(27.5,27.75]	5	5	6	6	10	10	6
(27.75,27.85]	7	6	7	6	7	6	10
(27.85,27.95]	8	6	5	5	10	5	9
(27.95,28.05]	6	3	4	6	10	6	10
(28.05,28.1]	5	3	4	7	6	10	10
(28.1,28.2]	4	6	6	5	10	10	7
(28.2,28.3]	5	7	8	5	8	10	6
(28.3,28.4]	6	10	9	6	5	4	8
total	88.00	92.00	90.00	88.00	100.00	100.00	113.00
Atlantic - obs
(-Inf,26]	8	9	6	6	7	5	8
(26,26.5]	8	10	8	7	4	4	7
(26.5,26.75]	7	9	7	6	4	4	10
(26.75,27]	4	5	5	8	5	10	10
(27,27.25]	6	8	5	8	6	6	10
(27.25,27.5]	7	4	7	7	10	3	9
(27.5,27.75]	7	3	9	7	4	10	7
(27.75,27.85]	7	5	10	5	5	10	5
(27.85,27.95]	5	4	7	6	10	10	5
(27.95,28.05]	5	4	7	6	10	10	5
(28.05,28.1]	4	7	5	9	10	7	6
(28.1,28.15]	4	7	5	10	6	5	10
(28.15,28.2]	5	5	5	8	4	10	10
(28.2, Inf]	8	4	9	6	10	6	5
total	85.00	84.00	95.00	99.00	95.00	100.00	107.00
Indo-Pacific - obs
(-Inf,26]	8	6	6	6	10	5	8
(26,26.5]	5	5	5	7	10	5	10
(26.5,26.75]	6	7	7	6	8	4	9
(26.75,27]	7	8	7	6	3	6	8
(27,27.25]	6	4	6	7	10	10	4
(27.25,27.5]	6	5	8	7	9	8	4
(27.5,27.75]	7	4	7	6	7	9	7
(27.75,27.85]	5	5	5	8	4	10	10
(27.85,27.95]	8	4	7	7	4	10	8
(27.95,28.05]	9	6	9	8	3	6	4
(28.05,28.1]	8	5	8	6	6	7	6
(28.1,28.2]	6	4	7	6	10	8	6
(28.2,28.3]	10	6	10	8	5	4	5
(28.3,28.4]	6	3	7	8	9	10	5
(28.4, Inf]	3	4	6	6	4	7	10
total	100.00	76.00	105.00	102.00	102.00	109.00	104.00

8.4 RMSE alternatives

8.4.1 AIC

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

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

Past versions of aic_vs_rmse_all_models-1.png

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

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

Past versions of aic_vs_rmse_selected_models-1.png

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

8.4.2 AIC vs BIC

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

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

Past versions of BIC_vs_AIC_all_models-1.png

Version	Author	Date
33ffcab	jens-daniel-mueller	2021-06-10
7e1f407	jens-daniel-mueller	2021-06-10
2cbe18c	jens-daniel-mueller	2021-06-10
69c79d0	jens-daniel-mueller	2021-06-08
1772903	jens-daniel-mueller	2021-06-07
594ed9a	jens-daniel-mueller	2021-06-04
db7df0e	jens-daniel-mueller	2021-06-04
2edc791	jens-daniel-mueller	2021-06-04
207339d	jens-daniel-mueller	2021-06-03
315710b	jens-daniel-mueller	2021-06-03
be90356	jens-daniel-mueller	2021-06-02
d37a85d	jens-daniel-mueller	2021-05-31
4b7a5ee	jens-daniel-mueller	2021-05-28
25bd183	jens-daniel-mueller	2021-05-26
62bd574	jens-daniel-mueller	2021-05-20
7c56c39	jens-daniel-mueller	2021-05-19
52e7583	jens-daniel-mueller	2021-05-12
969e631	jens-daniel-mueller	2021-05-12
d2a83bc	jens-daniel-mueller	2021-04-16
c0a47df	jens-daniel-mueller	2021-04-16
50290e8	jens-daniel-mueller	2021-04-16
a00ec94	jens-daniel-mueller	2021-04-16
b6fe355	jens-daniel-mueller	2021-04-16
81b7c6d	jens-daniel-mueller	2021-04-16
ddec5b7	jens-daniel-mueller	2021-04-15
29edae5	jens-daniel-mueller	2021-04-14
099d566	jens-daniel-mueller	2021-04-14
bb44686	jens-daniel-mueller	2021-04-14
bf40480	jens-daniel-mueller	2021-04-13
9f31fe3	jens-daniel-mueller	2021-04-13
338dd3c	jens-daniel-mueller	2021-04-09
a79ca2c	jens-daniel-mueller	2021-04-09
594d503	jens-daniel-mueller	2021-03-26
06bfaf7	jens-daniel-mueller	2021-03-26
2eb881a	jens-daniel-mueller	2021-03-26
03b6009	jens-daniel-mueller	2021-03-23
de6256b	jens-daniel-mueller	2021-03-23
f155edd	jens-daniel-mueller	2021-03-23
380d215	jens-daniel-mueller	2021-03-21
33b385b	jens-daniel-mueller	2021-03-20
108d6b2	jens-daniel-mueller	2021-03-20

8.5 RMSE vs VIF

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

Past versions of rmse_vs_VIF_all_models-1.png

Version	Author	Date
33ffcab	jens-daniel-mueller	2021-06-10
7e1f407	jens-daniel-mueller	2021-06-10
2cbe18c	jens-daniel-mueller	2021-06-10
69c79d0	jens-daniel-mueller	2021-06-08
1772903	jens-daniel-mueller	2021-06-07
594ed9a	jens-daniel-mueller	2021-06-04
db7df0e	jens-daniel-mueller	2021-06-04
2edc791	jens-daniel-mueller	2021-06-04
207339d	jens-daniel-mueller	2021-06-03
315710b	jens-daniel-mueller	2021-06-03
be90356	jens-daniel-mueller	2021-06-02
d37a85d	jens-daniel-mueller	2021-05-31
4b7a5ee	jens-daniel-mueller	2021-05-28
25bd183	jens-daniel-mueller	2021-05-26
62bd574	jens-daniel-mueller	2021-05-20
7c56c39	jens-daniel-mueller	2021-05-19
52e7583	jens-daniel-mueller	2021-05-12
969e631	jens-daniel-mueller	2021-05-12
d2a83bc	jens-daniel-mueller	2021-04-16
c0a47df	jens-daniel-mueller	2021-04-16
50290e8	jens-daniel-mueller	2021-04-16
a00ec94	jens-daniel-mueller	2021-04-16
b6fe355	jens-daniel-mueller	2021-04-16
81b7c6d	jens-daniel-mueller	2021-04-16
ddec5b7	jens-daniel-mueller	2021-04-15
29edae5	jens-daniel-mueller	2021-04-14
099d566	jens-daniel-mueller	2021-04-14
bb44686	jens-daniel-mueller	2021-04-14
bf40480	jens-daniel-mueller	2021-04-13
9f31fe3	jens-daniel-mueller	2021-04-13
338dd3c	jens-daniel-mueller	2021-04-09
a79ca2c	jens-daniel-mueller	2021-04-09
594d503	jens-daniel-mueller	2021-03-26
2eb881a	jens-daniel-mueller	2021-03-26
03b6009	jens-daniel-mueller	2021-03-23
de6256b	jens-daniel-mueller	2021-03-23
33b385b	jens-daniel-mueller	2021-03-20
108d6b2	jens-daniel-mueller	2021-03-20

8.6 Residual patterns

8.6.1 Fitted vs actual models

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

GLODAP_augmented_best <- left_join(
  lm_best_target %>% select(data_source, basin, gamma_slab, era, model),
  GLODAP_augmented
)

# calculate equal axis limits and binwidth
axis_lims <- GLODAP_augmented %>%
  summarise(
    max_value = max(
      c(max(.fitted, max(!!sym(params_local$MLR_target))))
      ),
    min_value = min(
      c(min(.fitted, min(!!sym(params_local$MLR_target))))
    )
  )

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

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

Past versions of fitted_vs_actual_all-1.png

Version	Author	Date
33ffcab	jens-daniel-mueller	2021-06-10
7e1f407	jens-daniel-mueller	2021-06-10
2cbe18c	jens-daniel-mueller	2021-06-10
69c79d0	jens-daniel-mueller	2021-06-08
1772903	jens-daniel-mueller	2021-06-07
594ed9a	jens-daniel-mueller	2021-06-04
db7df0e	jens-daniel-mueller	2021-06-04
2edc791	jens-daniel-mueller	2021-06-04
207339d	jens-daniel-mueller	2021-06-03
315710b	jens-daniel-mueller	2021-06-03
be90356	jens-daniel-mueller	2021-06-02
d37a85d	jens-daniel-mueller	2021-05-31
4b7a5ee	jens-daniel-mueller	2021-05-28
25bd183	jens-daniel-mueller	2021-05-26
62bd574	jens-daniel-mueller	2021-05-20
7c56c39	jens-daniel-mueller	2021-05-19
52e7583	jens-daniel-mueller	2021-05-12
969e631	jens-daniel-mueller	2021-05-12
d2a83bc	jens-daniel-mueller	2021-04-16
c0a47df	jens-daniel-mueller	2021-04-16
50290e8	jens-daniel-mueller	2021-04-16
a00ec94	jens-daniel-mueller	2021-04-16
b6fe355	jens-daniel-mueller	2021-04-16
81b7c6d	jens-daniel-mueller	2021-04-16
ddec5b7	jens-daniel-mueller	2021-04-15
29edae5	jens-daniel-mueller	2021-04-14
099d566	jens-daniel-mueller	2021-04-14
bb44686	jens-daniel-mueller	2021-04-14
bf40480	jens-daniel-mueller	2021-04-13
9f31fe3	jens-daniel-mueller	2021-04-13
338dd3c	jens-daniel-mueller	2021-04-09
a79ca2c	jens-daniel-mueller	2021-04-09
594d503	jens-daniel-mueller	2021-03-26
2eb881a	jens-daniel-mueller	2021-03-26
03b6009	jens-daniel-mueller	2021-03-23
de6256b	jens-daniel-mueller	2021-03-23
33b385b	jens-daniel-mueller	2021-03-20
108d6b2	jens-daniel-mueller	2021-03-20
83a13de	jens-daniel-mueller	2021-03-20
cf98c6d	jens-daniel-mueller	2021-03-16
a1d52ff	jens-daniel-mueller	2021-03-15
0bade3b	jens-daniel-mueller	2021-03-15
27c1f4b	jens-daniel-mueller	2021-03-14
af75ebf	jens-daniel-mueller	2021-03-14
5017709	jens-daniel-mueller	2021-03-11
585b07f	jens-daniel-mueller	2021-03-11
85a5ed2	jens-daniel-mueller	2021-03-10
6c0bec6	jens-daniel-mueller	2021-03-05
3c2ec33	jens-daniel-mueller	2021-03-05
af70b94	jens-daniel-mueller	2021-03-04
27ae473	jens-daniel-mueller	2021-02-24
4bc00ea	jens-daniel-mueller	2021-02-24

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

Past versions of fitted_vs_actual_best-1.png

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

rm(binwidth_value, axis_lims)

8.6.2 Pooled

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

i_ylim <- c(-30,30)

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

Past versions of residual_patterns_full_range-1.png

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

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

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

Past versions of residual_patterns_y_range_limit-1.png

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

8.6.3 By model domain

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

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

[[1]]

Past versions of residual_patterns_by_model_domain-1.png

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

[[2]]

Past versions of residual_patterns_by_model_domain-2.png

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

[[3]]

Past versions of residual_patterns_by_model_domain-3.png

Version	Author	Date
7e1f407	jens-daniel-mueller	2021-06-10
2cbe18c	jens-daniel-mueller	2021-06-10
1772903	jens-daniel-mueller	2021-06-07
594ed9a	jens-daniel-mueller	2021-06-04
db7df0e	jens-daniel-mueller	2021-06-04
207339d	jens-daniel-mueller	2021-06-03
315710b	jens-daniel-mueller	2021-06-03
be90356	jens-daniel-mueller	2021-06-02
d37a85d	jens-daniel-mueller	2021-05-31
4b7a5ee	jens-daniel-mueller	2021-05-28
25bd183	jens-daniel-mueller	2021-05-26
62bd574	jens-daniel-mueller	2021-05-20
7c56c39	jens-daniel-mueller	2021-05-19
52e7583	jens-daniel-mueller	2021-05-12
969e631	jens-daniel-mueller	2021-05-12
d2a83bc	jens-daniel-mueller	2021-04-16
c0a47df	jens-daniel-mueller	2021-04-16
50290e8	jens-daniel-mueller	2021-04-16
b6fe355	jens-daniel-mueller	2021-04-16
81b7c6d	jens-daniel-mueller	2021-04-16
ddec5b7	jens-daniel-mueller	2021-04-15
29edae5	jens-daniel-mueller	2021-04-14
099d566	jens-daniel-mueller	2021-04-14
bb44686	jens-daniel-mueller	2021-04-14
bf40480	jens-daniel-mueller	2021-04-13
9f31fe3	jens-daniel-mueller	2021-04-13
338dd3c	jens-daniel-mueller	2021-04-09
a79ca2c	jens-daniel-mueller	2021-04-09
d3c3fde	jens-daniel-mueller	2021-03-26
594d503	jens-daniel-mueller	2021-03-26

[[4]]

Past versions of residual_patterns_by_model_domain-4.png

Version	Author	Date
33ffcab	jens-daniel-mueller	2021-06-10
7e1f407	jens-daniel-mueller	2021-06-10
2cbe18c	jens-daniel-mueller	2021-06-10
69c79d0	jens-daniel-mueller	2021-06-08
1772903	jens-daniel-mueller	2021-06-07
594ed9a	jens-daniel-mueller	2021-06-04
db7df0e	jens-daniel-mueller	2021-06-04
2edc791	jens-daniel-mueller	2021-06-04
207339d	jens-daniel-mueller	2021-06-03
315710b	jens-daniel-mueller	2021-06-03
be90356	jens-daniel-mueller	2021-06-02
d37a85d	jens-daniel-mueller	2021-05-31
4b7a5ee	jens-daniel-mueller	2021-05-28
25bd183	jens-daniel-mueller	2021-05-26
62bd574	jens-daniel-mueller	2021-05-20
7c56c39	jens-daniel-mueller	2021-05-19
52e7583	jens-daniel-mueller	2021-05-12
969e631	jens-daniel-mueller	2021-05-12
d2a83bc	jens-daniel-mueller	2021-04-16
c0a47df	jens-daniel-mueller	2021-04-16
50290e8	jens-daniel-mueller	2021-04-16
a00ec94	jens-daniel-mueller	2021-04-16
b6fe355	jens-daniel-mueller	2021-04-16
81b7c6d	jens-daniel-mueller	2021-04-16
ddec5b7	jens-daniel-mueller	2021-04-15
29edae5	jens-daniel-mueller	2021-04-14
099d566	jens-daniel-mueller	2021-04-14
bb44686	jens-daniel-mueller	2021-04-14
bf40480	jens-daniel-mueller	2021-04-13
9f31fe3	jens-daniel-mueller	2021-04-13
338dd3c	jens-daniel-mueller	2021-04-09
a79ca2c	jens-daniel-mueller	2021-04-09
d3c3fde	jens-daniel-mueller	2021-03-26
594d503	jens-daniel-mueller	2021-03-26

8.6.4 Latitudinal mean

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

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

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

[[1]]

Past versions of latitudinal_residual-1.png

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

[[2]]

Past versions of latitudinal_residual-2.png

Version	Author	Date
33ffcab	jens-daniel-mueller	2021-06-10
7e1f407	jens-daniel-mueller	2021-06-10
2cbe18c	jens-daniel-mueller	2021-06-10
69c79d0	jens-daniel-mueller	2021-06-08
1772903	jens-daniel-mueller	2021-06-07
594ed9a	jens-daniel-mueller	2021-06-04
db7df0e	jens-daniel-mueller	2021-06-04
2edc791	jens-daniel-mueller	2021-06-04
207339d	jens-daniel-mueller	2021-06-03
315710b	jens-daniel-mueller	2021-06-03
be90356	jens-daniel-mueller	2021-06-02
d37a85d	jens-daniel-mueller	2021-05-31
4b7a5ee	jens-daniel-mueller	2021-05-28
25bd183	jens-daniel-mueller	2021-05-26
62bd574	jens-daniel-mueller	2021-05-20
7c56c39	jens-daniel-mueller	2021-05-19
52e7583	jens-daniel-mueller	2021-05-12
969e631	jens-daniel-mueller	2021-05-12
d2a83bc	jens-daniel-mueller	2021-04-16
c0a47df	jens-daniel-mueller	2021-04-16
50290e8	jens-daniel-mueller	2021-04-16
a00ec94	jens-daniel-mueller	2021-04-16
b6fe355	jens-daniel-mueller	2021-04-16
81b7c6d	jens-daniel-mueller	2021-04-16
ddec5b7	jens-daniel-mueller	2021-04-15
29edae5	jens-daniel-mueller	2021-04-14
099d566	jens-daniel-mueller	2021-04-14
bb44686	jens-daniel-mueller	2021-04-14
bf40480	jens-daniel-mueller	2021-04-13
9f31fe3	jens-daniel-mueller	2021-04-13
338dd3c	jens-daniel-mueller	2021-04-09
a79ca2c	jens-daniel-mueller	2021-04-09
594d503	jens-daniel-mueller	2021-03-26

8.6.5 Latitudinal offset

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

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

Past versions of latitudinal_residual_offset-1.png

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

9 VIF

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

Past versions of vif_distribution_n_predictors-1.png

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

GLODAP_glanced_vif %>% 
  ggplot(aes(Variables,
             log10(VIF),
             col = data_source)) + 
  geom_hline(yintercept = log10(params_local$vif_max)) +
  geom_violin() +
  scale_color_brewer(palette = "Set1") +
  facet_grid(gamma_slab~basin) +
  coord_flip()

Past versions of vif_distribution_type_predictor-1.png

Version	Author	Date
7e1f407	jens-daniel-mueller	2021-06-10
2cbe18c	jens-daniel-mueller	2021-06-10
594ed9a	jens-daniel-mueller	2021-06-04
db7df0e	jens-daniel-mueller	2021-06-04
207339d	jens-daniel-mueller	2021-06-03
315710b	jens-daniel-mueller	2021-06-03
d37a85d	jens-daniel-mueller	2021-05-31
25bd183	jens-daniel-mueller	2021-05-26
62bd574	jens-daniel-mueller	2021-05-20
7c56c39	jens-daniel-mueller	2021-05-19
52e7583	jens-daniel-mueller	2021-05-12
969e631	jens-daniel-mueller	2021-05-12
d2a83bc	jens-daniel-mueller	2021-04-16
c0a47df	jens-daniel-mueller	2021-04-16
50290e8	jens-daniel-mueller	2021-04-16
b6fe355	jens-daniel-mueller	2021-04-16
81b7c6d	jens-daniel-mueller	2021-04-16
ddec5b7	jens-daniel-mueller	2021-04-15
29edae5	jens-daniel-mueller	2021-04-14
099d566	jens-daniel-mueller	2021-04-14
bb44686	jens-daniel-mueller	2021-04-14
bf40480	jens-daniel-mueller	2021-04-13
9f31fe3	jens-daniel-mueller	2021-04-13
338dd3c	jens-daniel-mueller	2021-04-09
a79ca2c	jens-daniel-mueller	2021-04-09
d3c3fde	jens-daniel-mueller	2021-03-26
594d503	jens-daniel-mueller	2021-03-26
2eb881a	jens-daniel-mueller	2021-03-26
03b6009	jens-daniel-mueller	2021-03-23
de6256b	jens-daniel-mueller	2021-03-23
380d215	jens-daniel-mueller	2021-03-21
33b385b	jens-daniel-mueller	2021-03-20
108d6b2	jens-daniel-mueller	2021-03-20
330dcd0	jens-daniel-mueller	2021-03-20

Session information

sessionInfo()

R version 4.0.3 (2020-10-10)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: openSUSE Leap 15.2

Matrix products: default
BLAS:   /usr/local/R-4.0.3/lib64/R/lib/libRblas.so
LAPACK: /usr/local/R-4.0.3/lib64/R/lib/libRlapack.so

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
 [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8    
 [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8   
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C            
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
 [1] tictoc_1.0.1     gt_0.2.2         corrr_0.4.3      broom_0.7.5     
 [5] kableExtra_1.3.1 knitr_1.30       olsrr_0.5.3      GGally_2.0.0    
 [9] lubridate_1.7.9  metR_0.9.0       scico_1.2.0      patchwork_1.1.1 
[13] collapse_1.5.0   forcats_0.5.0    stringr_1.4.0    dplyr_1.0.5     
[17] purrr_0.3.4      readr_1.4.0      tidyr_1.1.2      tibble_3.0.4    
[21] ggplot2_3.3.3    tidyverse_1.3.0  workflowr_1.6.2 

loaded via a namespace (and not attached):
 [1] fs_1.5.0                 webshot_0.5.2            RColorBrewer_1.1-2      
 [4] httr_1.4.2               rprojroot_2.0.2          tools_4.0.3             
 [7] backports_1.1.10         utf8_1.1.4               R6_2.5.0                
[10] nortest_1.0-4            DBI_1.1.0                colorspace_1.4-1        
[13] withr_2.3.0              gridExtra_2.3            tidyselect_1.1.0        
[16] curl_4.3                 compiler_4.0.3           git2r_0.27.1            
[19] cli_2.1.0                rvest_0.3.6              xml2_1.3.2              
[22] sass_0.2.0               labeling_0.4.2           scales_1.1.1            
[25] checkmate_2.0.0          goftest_1.2-2            digest_0.6.27           
[28] foreign_0.8-80           rmarkdown_2.5            rio_0.5.16              
[31] pkgconfig_2.0.3          htmltools_0.5.0          highr_0.8               
[34] dbplyr_1.4.4             rlang_0.4.10             readxl_1.3.1            
[37] rstudioapi_0.11          farver_2.0.3             generics_0.0.2          
[40] jsonlite_1.7.1           zip_2.1.1                car_3.0-10              
[43] magrittr_1.5             Matrix_1.2-18            Rcpp_1.0.5              
[46] munsell_0.5.0            fansi_0.4.1              abind_1.4-5             
[49] lifecycle_1.0.0          stringi_1.5.3            whisker_0.4             
[52] yaml_2.2.1               carData_3.0-4            plyr_1.8.6              
[55] grid_4.0.3               blob_1.2.1               parallel_4.0.3          
[58] promises_1.1.1           crayon_1.3.4             lattice_0.20-41         
[61] haven_2.3.1              hms_0.5.3                pillar_1.4.7            
[64] reprex_0.3.0             glue_1.4.2               evaluate_0.14           
[67] RcppArmadillo_0.10.1.2.0 data.table_1.13.2        modelr_0.1.8            
[70] vctrs_0.3.5              httpuv_1.5.4             cellranger_1.1.0        
[73] gtable_0.3.0             reshape_0.8.8            assertthat_0.2.1        
[76] xfun_0.18                openxlsx_4.2.3           RcppEigen_0.3.3.7.0     
[79] later_1.1.0.1            viridisLite_0.3.0        ellipsis_0.3.1          
[82] here_0.1