eMLR - model fitting

workflowr

• Summary
• Checks
• Past versions

Last updated: 2021-06-07

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

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 7c32ee8. 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	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	108.13805	121.66083	0.6089516	2000-2009 –> 2010-2019	4.659344	764.5856	797.2617
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	850.18709	870.04106	6.6809767	2000-2009 –> 2010-2019	14.638708	1276.3708	1310.7676
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	62.52964	76.05242	0.3893963	2000-2009 –> 2010-2019	4.773633	737.0352	769.7113
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	79.98794	93.51072	0.4620889	2000-2009 –> 2010-2019	4.447238	732.7330	765.4092
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	842.78382	862.63780	6.4875578	2000-2009 –> 2010-2019	15.399487	1282.3307	1316.7274
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	61.72820	75.25098	0.3863487	2000-2009 –> 2010-2019	4.259063	707.9481	740.6243
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate	852.42457	869.44226	6.7942699	2000-2009 –> 2010-2019	14.847591	1278.0173	1307.5002
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	850.81035	870.66433	6.6975211	2000-2009 –> 2010-2019	14.686892	1277.4623	1311.8591
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	850.85870	870.71267	6.6988061	2000-2009 –> 2010-2019	14.684998	1277.4637	1311.8604
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + silicate	57.78180	71.30458	0.3716863	2000-2009 –> 2010-2019	4.664166	727.4649	760.1411
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + silicate	845.11761	864.97158	6.5479185	2000-2009 –> 2010-2019	14.533957	1271.7203	1306.1171
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + silicate + phosphate	78.02865	91.55142	0.4532975	2000-2009 –> 2010-2019	4.330503	724.5128	757.1889
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	78.61728	92.14006	0.4559210	2000-2009 –> 2010-2019	4.742602	747.9922	780.6683
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	107.39110	120.91388	0.6045085	2000-2009 –> 2010-2019	4.640648	763.0349	795.7111
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	851.07825	870.93223	6.7046450	2000-2009 –> 2010-2019	14.714418	1278.0312	1312.4279
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	57.27825	70.80103	0.3698559	2000-2009 –> 2010-2019	4.677980	727.7908	760.4670
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	77.83079	91.35357	0.4524190	2000-2009 –> 2010-2019	4.639812	741.8626	774.5388
obs	Atlantic	(-Inf,26]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	858.15126	878.00523	6.8954933	2000-2009 –> 2010-2019	15.386848	1291.9937	1326.3905
obs	Atlantic	(-Inf,26]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	857.22505	877.07902	6.8701958	2000-2009 –> 2010-2019	15.353387	1290.9540	1325.3508
obs	Atlantic	(-Inf,26]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	854.37113	874.22510	6.7928293	2000-2009 –> 2010-2019	15.239910	1287.5967	1321.9935
obs	Atlantic	(26,26.5]	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	4006.80533	4033.35454	6.1616397	2000-2009 –> 2010-2019	11.631002	7261.6299	7313.7021
mod	Atlantic	(26,26.5]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	2101.55149	2129.03997	3.9136226	2000-2009 –> 2010-2019	8.383998	3462.9334	3514.5188
obs	Atlantic	(26,26.5]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	4008.79247	4039.76655	6.1615755	2000-2009 –> 2010-2019	11.608207	7261.2858	7322.0367
mod	Atlantic	(26,26.5]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	2054.02664	2081.51513	3.6733241	2000-2009 –> 2010-2019	7.679700	3364.7801	3416.3655
obs	Atlantic	(26,26.5]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	3872.28083	3903.25491	5.5162992	2000-2009 –> 2010-2019	10.974060	7126.8970	7187.6479
mod	Atlantic	(26,26.5]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2027.54557	2055.03405	3.5458890	2000-2009 –> 2010-2019	7.937061	3380.6686	3432.2540
mod	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	2057.69518	2085.18366	3.6913357	2000-2009 –> 2010-2019	7.978964	3399.7937	3451.3792
obs	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	3889.57780	3920.55188	5.5941655	2000-2009 –> 2010-2019	11.120846	7157.2447	7217.9956
mod	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	2057.96867	2085.45715	3.6926820	2000-2009 –> 2010-2019	8.050089	3407.5255	3459.1109
obs	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	3890.47223	3921.44631	5.5982218	2000-2009 –> 2010-2019	11.198747	7171.9431	7232.6940
mod	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate	2031.57368	2059.06217	3.5649845	2000-2009 –> 2010-2019	8.303409	3419.8588	3471.4442
mod	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	2032.24756	2059.73604	3.5681891	2000-2009 –> 2010-2019	8.423895	3431.8285	3483.4139
mod	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	2058.95997	2086.44845	3.6975660	2000-2009 –> 2010-2019	8.127543	3416.1477	3467.7331
obs	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	3900.69201	3931.66609	5.6447778	2000-2009 –> 2010-2019	11.306290	7193.4267	7254.1776
obs	Atlantic	(26,26.5]	cstar_tref ~ temp + aou + nitrate + phosphate	4005.28481	4031.83402	6.1540521	2000-2009 –> 2010-2019	11.715217	7277.4208	7329.4930
mod	Atlantic	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	2099.13199	2126.62047	3.9010176	2000-2009 –> 2010-2019	8.411511	3464.6415	3516.2269
obs	Atlantic	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	4007.25821	4038.23229	6.1539195	2000-2009 –> 2010-2019	11.636493	7266.5918	7327.3427
mod	Atlantic	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	2095.54435	2123.03283	3.8824016	2000-2009 –> 2010-2019	8.429310	3464.7688	3516.3542
obs	Atlantic	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	4007.12183	4038.09591	6.1532394	2000-2009 –> 2010-2019	11.678648	7274.5493	7335.3002
obs	Atlantic	(26,26.5]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	4007.57272	4038.54680	6.1554881	2000-2009 –> 2010-2019	11.735914	7285.3044	7346.0553
mod	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	2760.49024	2789.75035	4.1546121	2000-2009 –> 2010-2019	7.322786	4662.4072	4719.0049
obs	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	4802.70508	4835.30216	5.2521390	2000-2009 –> 2010-2019	10.923485	8072.0099	8134.3297
mod	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	2763.54585	2792.80596	4.1677746	2000-2009 –> 2010-2019	7.226483	4639.6534	4696.2511
obs	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	4801.86186	4834.45895	5.2492936	2000-2009 –> 2010-2019	10.949819	8076.4627	8138.7825
mod	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	2806.21723	2835.47734	4.3560056	2000-2009 –> 2010-2019	7.208631	4631.1262	4687.7238
obs	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	4897.57126	4930.16834	5.5823141	2000-2009 –> 2010-2019	11.280183	8171.6911	8234.0109
mod	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	2812.37714	2841.63726	4.3838714	2000-2009 –> 2010-2019	7.182504	4623.2602	4679.8579
mod	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	2769.99929	2799.25941	4.1957110	2000-2009 –> 2010-2019	7.121011	4613.3736	4669.9713
obs	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	4798.96588	4831.56297	5.2395328	2000-2009 –> 2010-2019	10.938704	8073.3214	8135.6413
mod	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	2813.34771	2842.60782	4.3882782	2000-2009 –> 2010-2019	7.182973	4623.1972	4679.7949
obs	Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + phosphate	4864.60195	4892.54231	5.4723068	2000-2009 –> 2010-2019	11.147808	8132.6625	8186.0795
obs	Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + phosphate_star	4861.70188	4889.64224	5.4621171	2000-2009 –> 2010-2019	11.163024	8134.3717	8187.7887
mod	Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	2781.42810	2810.68822	4.2456455	2000-2009 –> 2010-2019	7.313782	4659.7946	4716.3923
obs	Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	4864.84587	4897.44296	5.4661344	2000-2009 –> 2010-2019	11.140914	8134.7751	8197.0949
mod	Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	2785.31491	2814.57502	4.2627627	2000-2009 –> 2010-2019	7.230803	4639.3357	4695.9333
obs	Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	4862.09426	4894.69135	5.4564767	2000-2009 –> 2010-2019	11.149566	8135.3480	8197.6679
mod	Atlantic	(26.5,26.75]	cstar_tref ~ temp + nitrate + phosphate + phosphate_star	2795.02730	2820.10740	4.3147618	2000-2009 –> 2010-2019	7.379675	4670.6224	4719.1347
obs	Atlantic	(26.5,26.75]	cstar_tref ~ temp + nitrate + phosphate + phosphate_star	4859.26611	4887.20647	5.4535733	2000-2009 –> 2010-2019	11.153797	8131.8123	8185.2293
mod	Atlantic	(26.5,26.75]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	2788.95970	2818.21981	4.2788768	2000-2009 –> 2010-2019	7.198814	4630.9871	4687.5847
obs	Atlantic	(26.5,26.75]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	4859.63165	4892.22874	5.4478478	2000-2009 –> 2010-2019	11.140850	8132.8695	8195.1894
mod	Atlantic	(26.75,27]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	4293.88861	4327.61342	2.5152123	2000-2009 –> 2010-2019	4.842928	7142.1985	7206.9988
mod	Atlantic	(26.75,27]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	4655.76363	4684.67062	3.0691896	2000-2009 –> 2010-2019	5.431857	7520.7334	7576.2765
mod	Atlantic	(26.75,27]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	4441.86457	4475.58938	2.7272857	2000-2009 –> 2010-2019	5.256735	7394.2336	7459.0339
mod	Atlantic	(26.75,27]	cstar_tref ~ sal + aou + silicate + phosphate	4586.00211	4614.90909	2.9542677	2000-2009 –> 2010-2019	5.490157	7539.5550	7595.0980
mod	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + aou + silicate + phosphate	4566.51769	4600.24251	2.9197495	2000-2009 –> 2010-2019	5.195139	7386.3619	7451.1621
obs	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + aou + silicate + phosphate	5698.78086	5732.95792	4.4652879	2000-2009 –> 2010-2019	10.843737	9724.0855	9789.1911
mod	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	4586.62456	4620.34937	2.9520423	2000-2009 –> 2010-2019	5.301499	7446.5742	7511.3745
obs	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	5703.92253	5738.09959	4.4770773	2000-2009 –> 2010-2019	10.799396	9718.3906	9783.4962
obs	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	5714.63928	5748.81634	4.5017501	2000-2009 –> 2010-2019	10.810359	9726.4460	9791.5516
mod	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	4593.90296	4627.62777	2.9638196	2000-2009 –> 2010-2019	5.346662	7471.5181	7536.3183
obs	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	5705.00716	5739.18422	4.4795683	2000-2009 –> 2010-2019	10.798760	9718.8688	9783.9744
obs	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + silicate + phosphate_star	5724.40337	5753.69799	4.5289906	2000-2009 –> 2010-2019	10.934763	9752.9486	9808.7534
mod	Atlantic	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	4313.92110	4347.64592	2.5429273	2000-2009 –> 2010-2019	4.818634	7133.9399	7198.7401
obs	Atlantic	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	5704.66179	5738.83885	4.4787750	2000-2009 –> 2010-2019	10.947751	9747.2444	9812.3500
mod	Atlantic	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	4323.99200	4357.71681	2.5569756	2000-2009 –> 2010-2019	4.849621	7153.2951	7218.0954
obs	Atlantic	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	5710.67703	5744.85409	4.4926121	2000-2009 –> 2010-2019	10.874112	9736.5679	9801.6735
mod	Atlantic	(26.75,27]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	4346.36171	4380.08652	2.5884582	2000-2009 –> 2010-2019	4.914309	7193.6681	7258.4684
obs	Atlantic	(26.75,27]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	5712.95889	5747.13595	4.4978724	2000-2009 –> 2010-2019	10.890925	9741.0673	9806.1729
obs	Atlantic	(26.75,27]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	5718.45447	5747.74909	4.5151950	2000-2009 –> 2010-2019	10.939109	9750.4669	9806.2717
obs	Atlantic	(26.75,27]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	5713.49479	5742.78942	4.5037255	2000-2009 –> 2010-2019	10.955144	9750.7453	9806.5501
mod	Atlantic	(27,27.25]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	5434.86920	5470.85238	2.0725732	2000-2009 –> 2010-2019	4.391010	9421.6530	9491.0876
mod	Atlantic	(27,27.25]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	5582.39939	5618.38256	2.1973276	2000-2009 –> 2010-2019	4.676243	9686.8432	9756.2778
mod	Atlantic	(27,27.25]	cstar_tref ~ sal + aou + silicate + phosphate	5729.68935	5760.53207	2.3312161	2000-2009 –> 2010-2019	4.794667	9821.1316	9880.6470
mod	Atlantic	(27,27.25]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	5730.51251	5766.49568	2.3301294	2000-2009 –> 2010-2019	4.675001	9737.2284	9806.6630
obs	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	11097.65645	11135.76619	6.1835756	2000-2009 –> 2010-2019	12.080973	19069.7331	19143.7367
obs	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	11072.00964	11110.11938	6.1373780	2000-2009 –> 2010-2019	11.920749	18995.4314	19069.4350
mod	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + aou + nitrate + silicate	5526.90071	5562.88388	2.1495392	2000-2009 –> 2010-2019	4.414504	9472.6640	9542.0987
mod	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + aou + silicate + phosphate	5729.98485	5765.96803	2.3296423	2000-2009 –> 2010-2019	4.709216	9762.5263	9831.9609
obs	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	11114.16998	11152.27972	6.2135053	2000-2009 –> 2010-2019	11.952591	19018.4363	19092.4399
obs	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	11078.84673	11116.95647	6.1496598	2000-2009 –> 2010-2019	11.966108	19016.4804	19090.4840
obs	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	11132.45128	11170.56102	6.2468080	2000-2009 –> 2010-2019	12.081479	19077.8801	19151.8837
obs	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	11115.81632	11153.92607	6.2164972	2000-2009 –> 2010-2019	11.968872	19025.8463	19099.8499
obs	Atlantic	(27,27.25]	cstar_tref ~ temp + aou + nitrate + phosphate_star	11091.68233	11124.34782	6.1763943	2000-2009 –> 2010-2019	12.087524	19067.5549	19130.9865
mod	Atlantic	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate	5690.88639	5721.72910	2.2956510	2000-2009 –> 2010-2019	4.598617	9663.9003	9723.4158
mod	Atlantic	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	5692.20130	5728.18447	2.2950280	2000-2009 –> 2010-2019	4.556014	9634.8736	9704.3082
obs	Atlantic	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	11057.76681	11095.87655	6.1118716	2000-2009 –> 2010-2019	12.086793	19062.3882	19136.3918
mod	Atlantic	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	5692.75334	5728.73651	2.2955300	2000-2009 –> 2010-2019	4.555340	9634.5108	9703.9455
obs	Atlantic	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	11017.26870	11055.37844	6.0399244	2000-2009 –> 2010-2019	11.882600	18966.1141	19040.1177
mod	Atlantic	(27,27.25]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	5611.04156	5647.02473	2.2224047	2000-2009 –> 2010-2019	4.484344	9554.4546	9623.8893
obs	Atlantic	(27,27.25]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	11028.32390	11066.43364	6.0594801	2000-2009 –> 2010-2019	11.936364	18991.7171	19065.7207
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	13209.26121	13248.70833	5.8610140	2000-2009 –> 2010-2019	10.541820	21903.1189	21979.5933
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + aou + nitrate + silicate	11037.42240	11071.17305	3.5662025	2000-2009 –> 2010-2019	6.894340	20706.4808	20773.3431
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	10897.40600	10936.78175	3.4447325	2000-2009 –> 2010-2019	6.734922	20526.2179	20604.2239
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	10910.83596	10950.21171	3.4560400	2000-2009 –> 2010-2019	6.768068	20564.0188	20642.0248
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	13194.31240	13233.75953	5.8398891	2000-2009 –> 2010-2019	10.442098	21838.5539	21915.0283
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	11158.80841	11198.18416	3.6716239	2000-2009 –> 2010-2019	6.936971	20759.6986	20837.7046
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	13193.03942	13232.48654	5.8380937	2000-2009 –> 2010-2019	10.415294	21821.3160	21897.7904
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + nitrate + silicate	11008.65905	11048.03480	3.5395315	2000-2009 –> 2010-2019	6.787734	20590.1554	20668.1614
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + phosphate	13259.81841	13293.63023	5.9358937	2000-2009 –> 2010-2019	10.552754	21911.3726	21976.9221
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + phosphate_star	13263.30383	13297.11565	5.9408932	2000-2009 –> 2010-2019	10.535115	21900.4557	21966.0052
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate	13256.51464	13295.96177	5.9282941	2000-2009 –> 2010-2019	10.543203	21908.8302	21985.3046
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	10891.15562	10930.53137	3.4394825	2000-2009 –> 2010-2019	7.104173	20917.0976	20995.1036
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	13259.78652	13299.23365	5.9329811	2000-2009 –> 2010-2019	10.527094	21898.8692	21975.3436
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	11154.21614	11193.59189	3.6675117	2000-2009 –> 2010-2019	6.800118	20602.2180	20680.2240
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	13189.94864	13229.39577	5.8337368	2000-2009 –> 2010-2019	10.416070	21821.5089	21897.9833
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	10869.06445	10908.44020	3.4209911	2000-2009 –> 2010-2019	6.958587	20764.9741	20842.9801
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + phosphate + phosphate_star	13260.04993	13293.86175	5.9362257	2000-2009 –> 2010-2019	10.534263	21899.6343	21965.1838
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	10892.41195	10931.78770	3.4405371	2000-2009 –> 2010-2019	7.107658	20920.7965	20998.8025
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	13256.72890	13296.17603	5.9286009	2000-2009 –> 2010-2019	10.526161	21898.0092	21974.4836
mod	Atlantic	(27.25,27.5]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	11161.91335	11201.28910	3.6744068	2000-2009 –> 2010-2019	6.800340	20602.0588	20680.0648
obs	Atlantic	(27.5,27.75]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	14167.22822	14207.47348	5.1106478	2000-2009 –> 2010-2019	10.381432	25291.7383	25370.4640
obs	Atlantic	(27.5,27.75]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	14176.83816	14211.33409	5.1234514	2000-2009 –> 2010-2019	10.395379	25300.1307	25367.6099
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + aou + nitrate + phosphate	18258.38137	18296.36655	2.1802088	2000-2009 –> 2010-2019	4.615583	33884.1040	33958.8448
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	17826.82625	17871.14229	2.0692477	2000-2009 –> 2010-2019	4.497887	33435.8210	33523.0185
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + aou + silicate + phosphate	18578.30431	18616.28949	2.2658854	2000-2009 –> 2010-2019	4.707834	34222.2558	34296.9965
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	18545.41015	18589.72619	2.2563794	2000-2009 –> 2010-2019	4.691136	34171.4181	34258.6156
obs	Atlantic	(27.5,27.75]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	14178.43811	14218.68337	5.1230097	2000-2009 –> 2010-2019	10.319146	25251.5131	25330.2388
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	18147.10046	18191.41651	2.1506548	2000-2009 –> 2010-2019	4.570354	33731.1595	33818.3570
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	18576.85334	18621.16938	2.2649435	2000-2009 –> 2010-2019	4.687234	34168.1478	34255.3453
obs	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	14184.21185	14224.45710	5.1293884	2000-2009 –> 2010-2019	10.332885	25262.3908	25341.1165
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	18383.91905	18428.23510	2.2129019	2000-2009 –> 2010-2019	4.592830	33855.9106	33943.1082
obs	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	14213.13506	14253.38032	5.1614621	2000-2009 –> 2010-2019	10.392577	25310.4035	25389.1292
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	18459.44219	18503.75824	2.2331294	2000-2009 –> 2010-2019	4.622611	33958.5244	34045.7219
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	18494.34935	18538.66540	2.2425410	2000-2009 –> 2010-2019	4.630103	33987.9975	34075.1950
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	18389.43056	18433.74661	2.2143718	2000-2009 –> 2010-2019	4.596065	33866.4363	33953.6338
obs	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	14193.20615	14233.45140	5.1393410	2000-2009 –> 2010-2019	10.355424	25280.0975	25358.8232
obs	Atlantic	(27.5,27.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	14173.99297	14214.23823	5.1181042	2000-2009 –> 2010-2019	10.405873	25310.1047	25388.8304
obs	Atlantic	(27.5,27.75]	cstar_tref ~ temp + aou + silicate + phosphate	14182.51797	14217.01390	5.1297268	2000-2009 –> 2010-2019	10.418101	25317.0427	25384.5219
obs	Atlantic	(27.5,27.75]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	14183.87290	14224.11815	5.1290137	2000-2009 –> 2010-2019	10.440777	25336.3104	25415.0361
obs	Atlantic	(27.5,27.75]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	14191.94840	14226.44433	5.1401632	2000-2009 –> 2010-2019	10.452096	25342.5009	25409.9801
obs	Atlantic	(27.75,27.85]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	9252.90290	9290.04392	5.3840807	2000-2009 –> 2010-2019	10.561226	16980.0211	17053.1286
obs	Atlantic	(27.75,27.85]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	9250.90677	9282.74193	5.3840877	2000-2009 –> 2010-2019	10.586053	16988.0679	17050.7315
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + nitrate + silicate	5231.59887	5263.56950	1.3426600	2000-2009 –> 2010-2019	2.580152	9503.1855	9566.1999
obs	Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	9239.90552	9277.04654	5.3606333	2000-2009 –> 2010-2019	10.553875	16974.8405	17047.9480
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + silicate + phosphate	5020.60430	5052.57493	1.2528028	2000-2009 –> 2010-2019	2.439233	9182.2107	9245.2251
obs	Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + silicate + phosphate	9237.98024	9269.81540	5.3607678	2000-2009 –> 2010-2019	10.581837	16984.3716	17047.0352
obs	Atlantic	(27.75,27.85]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	9241.28362	9278.42464	5.3631146	2000-2009 –> 2010-2019	10.587737	16991.3879	17064.4954
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + nitrate + silicate + phosphate_star	5624.44596	5656.41658	1.5274874	2000-2009 –> 2010-2019	2.893133	10153.2218	10216.2362
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	5616.19752	5648.16814	1.5233566	2000-2009 –> 2010-2019	2.894876	10156.1758	10219.1902
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + aou + nitrate + silicate	5189.45548	5226.75454	1.3233421	2000-2009 –> 2010-2019	2.557392	9455.7724	9529.2892
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + aou + silicate + phosphate	4905.61417	4942.91323	1.2055978	2000-2009 –> 2010-2019	2.387311	9058.8227	9132.3395
obs	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + aou + silicate + phosphate	9238.91361	9276.05463	5.3588481	2000-2009 –> 2010-2019	10.569068	16982.0671	17055.1746
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	5033.13903	5070.43809	1.2571432	2000-2009 –> 2010-2019	2.477980	9271.3599	9344.8766
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	5037.97734	5075.27640	1.2591416	2000-2009 –> 2010-2019	2.478597	9273.2414	9346.7582
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	5027.15135	5064.45041	1.2546744	2000-2009 –> 2010-2019	2.471537	9256.8609	9330.3777
obs	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	9240.53374	9277.67476	5.3617643	2000-2009 –> 2010-2019	10.581500	16988.2818	17061.3894
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + silicate + phosphate_star	5113.24630	5145.21692	1.2914913	2000-2009 –> 2010-2019	2.595767	9522.0170	9585.0313
obs	Atlantic	(27.75,27.85]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	9252.95473	9290.09575	5.3841745	2000-2009 –> 2010-2019	10.568256	16983.4445	17056.5521
obs	Atlantic	(27.75,27.85]	cstar_tref ~ temp + aou + silicate + phosphate	9250.95962	9282.79478	5.3841833	2000-2009 –> 2010-2019	10.594418	16992.1205	17054.7841
obs	Atlantic	(27.75,27.85]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	9253.29963	9290.44065	5.3847981	2000-2009 –> 2010-2019	10.576723	16987.5959	17060.7035
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + aou + nitrate + silicate	6331.04137	6362.92069	1.9885557	2000-2009 –> 2010-2019	4.907548	12837.4227	12900.3412
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	19500.54451	19542.93505	5.3238735	2000-2009 –> 2010-2019	10.594654	36171.0630	36254.7684
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	6331.48753	6368.68007	1.9875260	2000-2009 –> 2010-2019	4.828886	12769.5685	12842.9734
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	18900.48421	18942.87476	4.8404801	2000-2009 –> 2010-2019	9.877676	35325.9542	35409.6596
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + aou + silicate + phosphate	6109.68897	6141.56829	1.8471141	2000-2009 –> 2010-2019	4.581719	12445.8944	12508.8129
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + aou + silicate + phosphate_star	18905.25782	18941.59258	4.8456840	2000-2009 –> 2010-2019	9.883154	35329.0225	35400.7700
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	19299.12897	19341.51952	5.1564621	2000-2009 –> 2010-2019	10.350219	35890.0650	35973.7703
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	19341.63546	19377.97022	5.1929958	2000-2009 –> 2010-2019	10.390277	35934.2386	36005.9861
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + aou + nitrate + silicate	6173.37628	6210.56882	1.8854884	2000-2009 –> 2010-2019	4.791726	12670.3371	12743.7420
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + aou + phosphate	6564.83763	6596.71695	2.1497266	2000-2009 –> 2010-2019	4.901556	12917.4191	12980.3376
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate + phosphate	5814.24796	5851.44050	1.6727642	2000-2009 –> 2010-2019	4.407131	12152.2263	12225.6312
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate + phosphate	18738.86221	18781.25276	4.7179572	2000-2009 –> 2010-2019	9.724386	35131.2118	35214.9172
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	5917.26304	5954.45558	1.7312017	2000-2009 –> 2010-2019	4.490389	12278.8084	12352.2133
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	18720.54199	18762.93254	4.7042661	2000-2009 –> 2010-2019	9.705934	35107.7476	35191.4530
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate	18718.84725	18761.23780	4.7030016	2000-2009 –> 2010-2019	9.748345	35153.0545	35236.7598
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	5925.36453	5962.55708	1.7358832	2000-2009 –> 2010-2019	4.539323	12328.4055	12401.8104
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + silicate + phosphate	18743.57803	18779.91278	4.7229860	2000-2009 –> 2010-2019	9.769704	35177.2577	35249.0052
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	5907.05267	5944.24521	1.7253197	2000-2009 –> 2010-2019	4.475974	12260.5197	12333.9246
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	18737.61461	18780.00516	4.7170236	2000-2009 –> 2010-2019	9.722297	35128.7164	35212.4218
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + silicate + phosphate_star	6006.50708	6038.38641	1.7846646	2000-2009 –> 2010-2019	4.817215	12612.4248	12675.3433
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	19953.95824	19996.47642	5.4031060	2000-2009 –> 2010-2019	12.110370	37721.4194	37805.1716
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	19539.96560	19582.48379	5.0656838	2000-2009 –> 2010-2019	11.627802	37190.5124	37274.2647
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + aou + silicate + phosphate_star	19626.46629	19662.91045	5.1359986	2000-2009 –> 2010-2019	11.712442	37286.6667	37358.4543
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	19675.87661	19718.39479	5.1740672	2000-2009 –> 2010-2019	11.781071	37362.8521	37446.6043
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	19805.11370	19841.55786	5.2809235	2000-2009 –> 2010-2019	11.901810	37501.3062	37573.0939
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + aou + phosphate	9626.98497	9659.90465	3.5817240	2000-2009 –> 2010-2019	7.382106	17471.4754	17535.9540
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + aou + phosphate_star	9656.80741	9689.72709	3.6117866	2000-2009 –> 2010-2019	7.410990	17500.4152	17564.8938
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate + phosphate	9490.42882	9528.83511	3.4453009	2000-2009 –> 2010-2019	7.244822	17336.2750	17411.5000
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate + phosphate	19658.53666	19701.05485	5.1601113	2000-2009 –> 2010-2019	11.759273	37339.1657	37422.9179
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	9502.20210	9540.60840	3.4566881	2000-2009 –> 2010-2019	7.253480	17346.0044	17421.2295
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	19604.04193	19646.56011	5.1164962	2000-2009 –> 2010-2019	11.692398	37265.8020	37349.5543
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + nitrate + phosphate_star	9721.69859	9754.61827	3.6780751	2000-2009 –> 2010-2019	7.464184	17555.4872	17619.9658
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate	19686.11632	19728.63450	5.1823263	2000-2009 –> 2010-2019	11.897645	37459.9909	37543.7431
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	9476.62296	9515.02926	3.4319956	2000-2009 –> 2010-2019	7.182099	17285.2362	17360.4612
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + phosphate + phosphate_star	9651.98429	9684.90397	3.6069076	2000-2009 –> 2010-2019	7.406109	17495.5910	17560.0696
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + phosphate_star	9735.09741	9762.53047	3.6939833	2000-2009 –> 2010-2019	7.496074	17578.8657	17632.5979
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + silicate + phosphate	19789.78958	19826.23374	5.2683333	2000-2009 –> 2010-2019	11.988143	37565.2370	37637.0246
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	9502.44386	9540.85016	3.4569223	2000-2009 –> 2010-2019	7.253720	17346.2501	17421.4751
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	19642.84935	19685.36754	5.1475179	2000-2009 –> 2010-2019	11.738002	37316.4462	37400.1984
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + silicate + phosphate_star	9501.02024	9533.93992	3.4574808	2000-2009 –> 2010-2019	7.254286	17342.8323	17407.3109
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	14409.73358	14450.02695	5.2716323	2000-2009 –> 2010-2019	11.609033	26002.8674	26081.5238
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	14354.39055	14394.68391	5.2095547	2000-2009 –> 2010-2019	11.546752	25947.4110	26026.0674
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	14357.43256	14397.72592	5.2129478	2000-2009 –> 2010-2019	11.537400	25943.3139	26021.9703
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	14443.37146	14477.90863	5.3119980	2000-2009 –> 2010-2019	11.659958	26040.4088	26107.8285
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	14475.81145	14516.10482	5.3467208	2000-2009 –> 2010-2019	11.651918	26050.8808	26129.5372
mod	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + aou + phosphate	3570.08772	3599.76380	1.3408811	2000-2009 –> 2010-2019	2.334736	5737.9818	5795.4814
mod	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + aou + phosphate_star	3609.39251	3639.06860	1.3664849	2000-2009 –> 2010-2019	2.378721	5805.2511	5862.7507
mod	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + aou + silicate + phosphate	3570.61954	3605.24164	1.3399341	2000-2009 –> 2010-2019	2.315181	5711.6717	5778.7546
mod	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	3607.48895	3642.11105	1.3639204	2000-2009 –> 2010-2019	2.362256	5784.2471	5851.3300
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	14390.77390	14431.06727	5.2502826	2000-2009 –> 2010-2019	11.521704	25946.7967	26025.4531
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	14468.71506	14509.00843	5.3386057	2000-2009 –> 2010-2019	11.647085	26045.6293	26124.2856
mod	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + nitrate + phosphate_star	3655.86044	3685.53653	1.3973863	2000-2009 –> 2010-2019	2.420624	5868.2152	5925.7148
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + nitrate + phosphate_star	14477.37822	14511.91539	5.3508042	2000-2009 –> 2010-2019	11.659711	26052.5328	26119.9525
mod	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	3628.76634	3663.38844	1.3779578	2000-2009 –> 2010-2019	2.400372	5841.8934	5908.9763
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	14357.57509	14397.86846	5.2131068	2000-2009 –> 2010-2019	11.504188	25924.6962	26003.3526
mod	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + phosphate + phosphate_star	3598.19640	3627.87248	1.3591422	2000-2009 –> 2010-2019	2.366975	5787.4028	5844.9024
mod	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	3597.97098	3632.59308	1.3576874	2000-2009 –> 2010-2019	2.349540	5764.7878	5831.8707
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	14432.04775	14472.34112	5.2968707	2000-2009 –> 2010-2019	11.594291	26002.7404	26081.3968
mod	Atlantic	(28.05,28.1]	cstar_tref ~ temp + aou + phosphate	3586.41903	3611.14910	1.3527622	2000-2009 –> 2010-2019	2.430590	5876.0889	5924.0053
mod	Atlantic	(28.05,28.1]	cstar_tref ~ temp + aou + silicate + phosphate	3588.41394	3618.09003	1.3527589	2000-2009 –> 2010-2019	2.406080	5844.9862	5902.4858
obs	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	16565.88872	16606.89379	5.9378630	2000-2009 –> 2010-2019	15.026973	27624.7320	27703.0361
mod	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + aou + phosphate	4579.69962	4611.70167	1.0755141	2000-2009 –> 2010-2019	2.402242	8951.3383	9014.2727
mod	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + aou + phosphate_star	4627.62406	4659.62612	1.0924798	2000-2009 –> 2010-2019	2.427506	9015.2363	9078.1707
mod	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + aou + silicate + phosphate	4550.65705	4587.99278	1.0646656	2000-2009 –> 2010-2019	2.272842	8684.4834	8757.9069
mod	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	4607.34783	4644.68356	1.0845607	2000-2009 –> 2010-2019	2.306690	8770.5914	8844.0149
obs	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	16600.57872	16641.58379	5.9778237	2000-2009 –> 2010-2019	15.083184	27664.8632	27743.1673
obs	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	16549.07677	16590.08185	5.9185929	2000-2009 –> 2010-2019	15.001554	27605.8586	27684.1628
mod	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + nitrate + phosphate_star	4648.10126	4680.10331	1.0998102	2000-2009 –> 2010-2019	2.409347	8986.3262	9049.2607
mod	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	4648.79852	4686.13425	1.0993424	2000-2009 –> 2010-2019	2.378672	8929.2329	9002.6564
obs	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	16574.83696	16615.84203	5.9481452	2000-2009 –> 2010-2019	15.057254	27640.3749	27718.6790
mod	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + phosphate + phosphate_star	4612.47147	4644.47353	1.0870869	2000-2009 –> 2010-2019	2.416454	8989.2012	9052.1357
obs	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + phosphate + phosphate_star	16601.64855	16636.79576	5.9813729	2000-2009 –> 2010-2019	15.085170	27663.4099	27730.5278
mod	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	4586.71107	4624.04680	1.0772758	2000-2009 –> 2010-2019	2.288670	8727.3521	8800.7756
obs	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	16600.48218	16641.48725	5.9777122	2000-2009 –> 2010-2019	15.079728	27663.6474	27741.9515
obs	Atlantic	(28.1,28.15]	cstar_tref ~ temp + aou + nitrate + phosphate_star	16592.85309	16628.00029	5.9712097	2000-2009 –> 2010-2019	15.074677	27654.5042	27721.6221
obs	Atlantic	(28.1,28.15]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	16568.76564	16609.77071	5.9411669	2000-2009 –> 2010-2019	15.043889	27632.1672	27710.4713
mod	Atlantic	(28.1,28.15]	cstar_tref ~ temp + aou + silicate + phosphate	4678.91640	4710.91846	1.1109343	2000-2009 –> 2010-2019	2.397291	8971.3834	9034.3179
obs	Atlantic	(28.1,28.15]	cstar_tref ~ temp + nitrate + phosphate + phosphate_star	16579.78663	16614.93384	5.9561432	2000-2009 –> 2010-2019	15.054372	27639.6844	27706.8022
obs	Atlantic	(28.1,28.15]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	16557.36961	16598.37469	5.9280905	2000-2009 –> 2010-2019	15.025722	27619.0674	27697.3715
mod	Atlantic	(28.1,28.15]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	4725.93675	4757.93881	1.1281256	2000-2009 –> 2010-2019	2.419989	9029.3482	9092.2826
mod	Atlantic	(28.15,28.2]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	7269.90172	7311.69206	0.8479836	2000-2009 –> 2010-2019	1.898949	13663.7060	13745.2802
mod	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + aou + nitrate	7329.78543	7365.60573	0.8571018	2000-2009 –> 2010-2019	1.899307	13685.2298	13755.1506
mod	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + aou + nitrate + silicate	7232.18307	7273.97342	0.8424736	2000-2009 –> 2010-2019	1.880842	13573.6079	13655.1821
obs	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + aou + silicate + phosphate	5567.19947	5600.90127	5.0983689	2000-2009 –> 2010-2019	12.844217	10191.1065	10256.3068
obs	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	5527.22379	5560.92559	4.9877261	2000-2009 –> 2010-2019	12.651481	10136.9597	10202.1600
mod	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + nitrate + phosphate_star	7641.57719	7677.39749	0.9045558	2000-2009 –> 2010-2019	1.964958	14072.2127	14142.1334
mod	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	7605.25953	7647.04988	0.8985852	2000-2009 –> 2010-2019	1.952569	14011.5238	14093.0980
obs	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	5512.58820	5546.29000	4.9478217	2000-2009 –> 2010-2019	12.660703	10130.8227	10196.0230
obs	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	5543.57605	5577.27785	5.0326919	2000-2009 –> 2010-2019	12.708132	10155.3384	10220.5387
mod	Atlantic	(28.15,28.2]	cstar_tref ~ temp + aou + nitrate + phosphate	7335.49274	7371.31304	0.8579477	2000-2009 –> 2010-2019	1.899275	13687.2763	13757.1970
mod	Atlantic	(28.15,28.2]	cstar_tref ~ temp + aou + nitrate + phosphate_star	7368.79651	7404.61681	0.8629002	2000-2009 –> 2010-2019	1.909780	13743.6806	13813.6013
mod	Atlantic	(28.15,28.2]	cstar_tref ~ temp + aou + nitrate + silicate	7259.51626	7295.33656	0.8467555	2000-2009 –> 2010-2019	1.897378	13649.9051	13719.8258
obs	Atlantic	(28.15,28.2]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	5527.89950	5561.60130	4.9895762	2000-2009 –> 2010-2019	12.994247	10195.5213	10260.7216
obs	Atlantic	(28.15,28.2]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	5499.56379	5533.26559	4.9125787	2000-2009 –> 2010-2019	12.818560	10150.6860	10215.8864
obs	Atlantic	(28.15,28.2]	cstar_tref ~ temp + aou + silicate + phosphate_star	5527.79304	5556.68029	4.9947644	2000-2009 –> 2010-2019	12.921851	10180.4610	10236.3470
mod	Atlantic	(28.15,28.2]	cstar_tref ~ temp + nitrate + phosphate + phosphate_star	7419.56795	7455.38825	0.8705054	2000-2009 –> 2010-2019	1.921363	13810.9296	13880.8503
obs	Atlantic	(28.15,28.2]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	5505.97307	5539.67487	4.9298902	2000-2009 –> 2010-2019	12.865273	10162.0324	10227.2327
mod	Atlantic	(28.15,28.2]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	7622.50659	7658.32688	0.9015793	2000-2009 –> 2010-2019	1.955710	14027.3754	14097.2961
obs	Atlantic	(28.15,28.2]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	5513.58796	5542.47522	4.9559745	2000-2009 –> 2010-2019	12.982152	10182.7783	10238.6643
obs	Atlantic	(28.15,28.2]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	5543.50167	5572.38893	5.0380137	2000-2009 –> 2010-2019	12.986340	10199.7286	10255.6145
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate	10073.55970	10105.95196	0.6887820	2000-2009 –> 2010-2019	1.440137	15918.8251	15980.4854
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate + phosphate	10010.76617	10049.63688	0.6841584	2000-2009 –> 2010-2019	1.433406	15843.5700	15917.5623
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	10756.31174	10795.80605	3.1848065	2000-2009 –> 2010-2019	6.499328	18334.1699	18410.5952
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate + phosphate_star	9925.87865	9964.74937	0.6781479	2000-2009 –> 2010-2019	1.427413	15758.8013	15832.7936
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate + silicate	9964.42392	10003.29464	0.6808706	2000-2009 –> 2010-2019	1.415979	15699.1107	15773.1030
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	10754.95091	10794.44522	3.1837669	2000-2009 –> 2010-2019	6.500558	18334.7867	18411.2121
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	9322.47112	9367.82029	0.6367858	2000-2009 –> 2010-2019	1.359948	14974.7801	15061.1045
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + phosphate + phosphate_star	10757.52063	10791.37290	3.1872594	2000-2009 –> 2010-2019	6.503196	18334.6118	18400.1193
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + silicate + phosphate	10753.49968	10787.35194	3.1841861	2000-2009 –> 2010-2019	6.501288	18331.6065	18397.1140
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	10735.51552	10775.00983	3.1689555	2000-2009 –> 2010-2019	6.482760	18312.7483	18389.1736
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + aou + nitrate	9992.41449	10031.28521	0.6828545	2000-2009 –> 2010-2019	1.432734	15829.5582	15903.5505
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	9927.89933	9973.24849	0.6781494	2000-2009 –> 2010-2019	1.427175	15761.1764	15847.5008
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	10754.42541	10793.91972	3.1833655	2000-2009 –> 2010-2019	6.501840	18335.7275	18412.1529
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	9922.70198	9968.05114	0.6777831	2000-2009 –> 2010-2019	1.427023	15757.4480	15843.7724
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + aou + nitrate + silicate	9885.04604	9930.39521	0.6751352	2000-2009 –> 2010-2019	1.407906	15605.3302	15691.6546
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + aou + silicate + phosphate	10725.84487	10765.33918	3.1616114	2000-2009 –> 2010-2019	6.478698	18305.9385	18382.3638
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	10768.48288	10807.97719	3.1941202	2000-2009 –> 2010-2019	6.468582	18311.1993	18387.6246
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	9967.86135	10013.21052	0.6809723	2000-2009 –> 2010-2019	1.433110	15822.4909	15908.8153
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	10796.19888	10835.69319	3.2154309	2000-2009 –> 2010-2019	6.497208	18345.3640	18421.7894
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	10741.47978	10780.97409	3.1734934	2000-2009 –> 2010-2019	6.455292	18290.6636	18367.0889
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	22577.80870	22621.17603	5.4421257	2000-2009 –> 2010-2019	10.981211	38089.1098	38173.1751
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	29332.57567	29377.73006	5.5551593	2000-2009 –> 2010-2019	11.853484	50322.9300	50410.6200
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + aou + phosphate + phosphate_star	22904.98847	22942.16047	5.6949873	2000-2009 –> 2010-2019	11.589923	38722.4972	38794.5532
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	22681.80055	22725.16789	5.5207703	2000-2009 –> 2010-2019	11.175023	38294.9658	38379.0311
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	29643.78771	29688.94209	5.7430502	2000-2009 –> 2010-2019	12.265159	50858.8488	50946.5388
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate	30113.58068	30152.28444	6.0400799	2000-2009 –> 2010-2019	12.523425	51288.2753	51363.4381
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	22688.30903	22731.67637	5.5257300	2000-2009 –> 2010-2019	11.342107	38441.4084	38525.4737
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	29791.26556	29836.41994	5.8342947	2000-2009 –> 2010-2019	12.154955	50804.4035	50892.0935
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	22790.61272	22833.98006	5.6042773	2000-2009 –> 2010-2019	11.523436	38630.4205	38714.4858
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	29844.66095	29889.81533	5.8676866	2000-2009 –> 2010-2019	12.199008	50868.6446	50956.3346
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + silicate	23253.75762	23297.12496	5.9740781	2000-2009 –> 2010-2019	12.164857	39315.6550	39399.7203
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + silicate	30031.18860	30076.34298	5.9858428	2000-2009 –> 2010-2019	12.465944	51204.6629	51292.3529
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + silicate + phosphate	23130.77601	23174.14334	5.8735671	2000-2009 –> 2010-2019	11.970418	39116.9955	39201.0608
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + silicate + phosphate	30107.05093	30152.20531	6.0345758	2000-2009 –> 2010-2019	12.651181	51414.6917	51502.3817
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	23327.70471	23371.07205	6.0353401	2000-2009 –> 2010-2019	12.308762	39455.2444	39539.3097
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	30219.00048	30264.15487	6.1072164	2000-2009 –> 2010-2019	12.760306	51562.0321	51649.7221
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	22890.67874	22934.04607	5.6821866	2000-2009 –> 2010-2019	11.693343	38806.8292	38890.8945
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	29973.24482	30018.39920	5.9488857	2000-2009 –> 2010-2019	12.324485	51042.0826	51129.7726
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	23397.05830	23440.42564	6.0933673	2000-2009 –> 2010-2019	12.403199	39553.2437	39637.3090
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	30268.35043	30313.50481	6.1395152	2000-2009 –> 2010-2019	12.807956	51626.2161	51713.9061
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	16270.39696	16312.14361	4.0888385	2000-2009 –> 2010-2019	8.538009	30499.1026	30581.4504
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	16568.79934	16610.54600	4.3066356	2000-2009 –> 2010-2019	8.575922	30596.0198	30678.3676
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	15846.86587	15888.61253	3.7984890	2000-2009 –> 2010-2019	7.865486	29637.1049	29719.4527
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	16255.85194	16297.59859	4.0785085	2000-2009 –> 2010-2019	8.358451	30295.2424	30377.5902
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	17484.17368	17526.11719	4.6417817	2000-2009 –> 2010-2019	9.215483	32170.9658	32253.6664
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	17503.40832	17545.35183	4.6569031	2000-2009 –> 2010-2019	9.247773	32208.9042	32291.6048
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + phosphate	17771.33591	17807.28750	4.8743794	2000-2009 –> 2010-2019	9.738221	32763.1654	32834.0517
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate	15467.08445	15508.83111	3.5557088	2000-2009 –> 2010-2019	7.634672	29271.6670	29354.0148
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate	17348.69776	17390.64127	4.5366580	2000-2009 –> 2010-2019	9.213879	32147.2181	32229.9188
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	15612.75806	15654.50472	3.6469418	2000-2009 –> 2010-2019	7.894891	29615.5184	29697.8662
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	17419.93080	17461.87432	4.5916316	2000-2009 –> 2010-2019	9.305147	32257.0388	32339.7395
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	17506.67624	17548.61975	4.6594771	2000-2009 –> 2010-2019	9.255783	32218.0796	32300.7802
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + nitrate + phosphate_star	17562.55236	17598.50395	4.7053000	2000-2009 –> 2010-2019	9.473293	32455.0252	32525.9115
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	17330.46134	17372.40485	4.5226903	2000-2009 –> 2010-2019	9.236711	32168.1045	32250.8051
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	15654.97186	15696.71852	3.6738145	2000-2009 –> 2010-2019	7.954058	29694.7055	29777.0533
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	17419.38770	17461.33122	4.5912099	2000-2009 –> 2010-2019	9.306488	32258.3617	32341.0623
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + silicate + phosphate_star	17418.34256	17454.29415	4.5919513	2000-2009 –> 2010-2019	9.326823	32276.0170	32346.9033
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	16016.05963	16057.80629	3.9119202	2000-2009 –> 2010-2019	8.190207	30053.5622	30135.9100
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	16079.81026	16121.55691	3.9555333	2000-2009 –> 2010-2019	8.243849	30128.7425	30211.0903
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	16126.32765	16168.07430	3.9876632	2000-2009 –> 2010-2019	8.293534	30195.2055	30277.5533
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	15814.96516	15849.92382	5.6637558	2000-2009 –> 2010-2019	11.045395	32403.6874	32473.9895
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	14197.40173	14238.55677	3.5440789	2000-2009 –> 2010-2019	7.881213	27025.9179	27107.0098
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	14238.86160	14274.13735	3.5733485	2000-2009 –> 2010-2019	8.189901	27342.5866	27412.0940
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	15663.77521	15704.56031	5.4932639	2000-2009 –> 2010-2019	10.521084	31891.0705	31973.0897
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	15687.81113	15728.59623	5.5196711	2000-2009 –> 2010-2019	10.605615	31976.5310	32058.5502
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + aou + phosphate	15869.03284	15903.99150	5.7251851	2000-2009 –> 2010-2019	11.101138	32452.1058	32522.4080
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	14020.91333	14062.06837	3.4276598	2000-2009 –> 2010-2019	7.794599	26879.8371	26960.9290
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	15815.73675	15856.52185	5.6623678	2000-2009 –> 2010-2019	11.029274	32391.8091	32473.8283
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	14190.87893	14232.03397	3.5397066	2000-2009 –> 2010-2019	8.019949	27163.5327	27244.6246
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	15690.27077	15731.05587	5.5223805	2000-2009 –> 2010-2019	10.617599	31988.7268	32070.7460
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + nitrate + phosphate_star	15711.01170	15745.97036	5.5474942	2000-2009 –> 2010-2019	10.755214	32124.1787	32194.4808
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	14706.24595	14747.40099	3.9023438	2000-2009 –> 2010-2019	8.515989	27809.1739	27890.2657
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	15696.56014	15737.34525	5.5293147	2000-2009 –> 2010-2019	10.668814	32041.2591	32123.2783
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	14216.09867	14257.25371	3.5566415	2000-2009 –> 2010-2019	8.045430	27197.2138	27278.3056
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	15851.46154	15892.24664	5.7028725	2000-2009 –> 2010-2019	11.099915	32457.4579	32539.4771
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	14366.46043	14407.61547	3.6593032	2000-2009 –> 2010-2019	8.090459	27290.1995	27371.2913
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	15795.65578	15836.44088	5.6397265	2000-2009 –> 2010-2019	11.059160	32423.7764	32505.7956
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	14582.84818	14624.00322	3.8122679	2000-2009 –> 2010-2019	8.339749	27602.0708	27683.1627
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + aou + silicate + phosphate	14411.73431	14447.01006	3.6921885	2000-2009 –> 2010-2019	8.443237	27642.9632	27712.4705
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	14625.13830	14666.29334	3.8429015	2000-2009 –> 2010-2019	8.387144	27660.7679	27741.8598
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	18476.46288	18513.16675	3.8011632	2000-2009 –> 2010-2019	7.728550	36837.4975	36910.7968
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	24805.80078	24843.59201	5.2930918	2000-2009 –> 2010-2019	10.334391	49396.0610	49471.6867
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	17488.54473	17531.36592	3.2793535	2000-2009 –> 2010-2019	7.199506	35839.4390	35924.9550
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ aou + phosphate + phosphate_star	25113.78552	25145.27821	5.5012593	2000-2009 –> 2010-2019	10.720594	49982.9558	50045.9772
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	17819.09418	17861.91537	3.4450989	2000-2009 –> 2010-2019	7.146135	35791.2929	35876.8088
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	24983.65780	25027.74758	5.4102006	2000-2009 –> 2010-2019	10.548737	49730.5495	49818.7794
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	17993.06258	18035.88377	3.5356690	2000-2009 –> 2010-2019	7.316299	36105.3532	36190.8691
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + aou + silicate + phosphate	25230.24130	25274.33107	5.5787851	2000-2009 –> 2010-2019	10.743579	50018.3869	50106.6168
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	18014.93694	18057.75813	3.5472244	2000-2009 –> 2010-2019	7.331599	36133.7470	36219.2629
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + nitrate + phosphate_star	18498.40483	18535.10870	3.8136246	2000-2009 –> 2010-2019	7.620923	36654.9767	36728.2761
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	17925.89872	17968.71991	3.5004239	2000-2009 –> 2010-2019	7.303625	36077.3815	36162.8974
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	25115.69066	25159.78044	5.4998254	2000-2009 –> 2010-2019	10.587155	49781.5188	49869.7487
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + nitrate + phosphate	25048.56880	25086.36004	5.4554362	2000-2009 –> 2010-2019	10.656098	49890.7301	49966.3557
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	17469.33490	17512.15608	3.2699702	2000-2009 –> 2010-2019	7.199333	35835.6808	35921.1967
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	25020.11040	25064.20018	5.4347979	2000-2009 –> 2010-2019	10.520480	49683.3157	49771.5456
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	17651.61637	17694.43756	3.3601003	2000-2009 –> 2010-2019	7.388645	36182.0868	36267.6027
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	25162.21878	25206.30855	5.5317616	2000-2009 –> 2010-2019	10.696445	49950.1914	50038.4213
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	17679.52477	17722.34596	3.3741174	2000-2009 –> 2010-2019	7.407834	36218.4427	36303.9586
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	25162.83981	25206.92958	5.5321891	2000-2009 –> 2010-2019	10.695519	49948.6909	50036.9208
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	25161.13345	25198.92469	5.5323912	2000-2009 –> 2010-2019	10.697166	49947.2488	50022.8744
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	21209.86378	21248.24071	2.6474499	2000-2009 –> 2010-2019	5.829324	43157.4090	43233.9238
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	21162.70045	21207.47353	2.6328001	2000-2009 –> 2010-2019	5.644250	42643.5608	42732.8282
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	31749.81850	31795.80832	4.9110055	2000-2009 –> 2010-2019	9.664263	62911.2786	63003.2036
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	31955.49319	32001.48302	5.0077601	2000-2009 –> 2010-2019	9.869118	63352.1745	63444.0994
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + aou + silicate + phosphate	31750.83988	31790.25973	4.9124132	2000-2009 –> 2010-2019	9.666733	62912.6368	62991.4297
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + aou + silicate + phosphate_star	31960.11323	31999.53309	5.0109058	2000-2009 –> 2010-2019	9.891456	63396.0090	63474.8018
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	31884.10259	31930.09242	4.9739620	2000-2009 –> 2010-2019	9.793970	63191.4322	63283.3572
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	31890.24974	31929.66959	4.9778075	2000-2009 –> 2010-2019	9.819676	63242.9115	63321.7044
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	21163.38399	21208.15707	2.6330033	2000-2009 –> 2010-2019	5.797993	43067.6307	43156.8981
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	21562.23147	21607.00456	2.7542406	2000-2009 –> 2010-2019	6.015687	43722.0763	43811.3436
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + aou + silicate + phosphate	31749.97594	31795.96577	4.9110788	2000-2009 –> 2010-2019	9.663046	62908.5974	63000.5224
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	31926.66304	31972.65287	4.9940836	2000-2009 –> 2010-2019	9.805036	63214.3214	63306.2464
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	21589.78732	21634.56040	2.7628200	2000-2009 –> 2010-2019	6.032106	43770.0714	43859.3387
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	22051.23409	22096.00718	2.9105263	2000-2009 –> 2010-2019	6.152113	44159.0755	44248.3428
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	32017.38031	32063.37014	5.0372448	2000-2009 –> 2010-2019	9.890011	63395.5594	63487.4844
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	31882.16526	31928.15509	4.9730480	2000-2009 –> 2010-2019	9.772803	63145.4511	63237.3761
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ temp + aou + nitrate + phosphate	21541.61232	21579.98925	2.7484587	2000-2009 –> 2010-2019	5.976333	43611.3642	43687.8791
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	21542.37996	21587.15305	2.7480764	2000-2009 –> 2010-2019	5.867787	43323.9428	43413.2101
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	21902.41525	21947.18834	2.8620474	2000-2009 –> 2010-2019	6.100405	44001.7732	44091.0405
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	21923.08493	21967.85802	2.8687321	2000-2009 –> 2010-2019	6.110848	44032.3175	44121.5849
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	33835.55869	33875.64374	4.2751486	2000-2009 –> 2010-2019	8.639672	66584.2790	66664.2100
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	33832.47225	33879.23814	4.2733028	2000-2009 –> 2010-2019	8.580953	66434.7410	66527.9939
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	29297.14425	29343.65350	3.1904604	2000-2009 –> 2010-2019	6.257089	56991.8821	57084.6150
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	33522.48330	33569.24919	4.1622995	2000-2009 –> 2010-2019	8.401378	65943.1361	66036.3890
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + aou + silicate + phosphate	33527.31949	33567.40454	4.1647161	2000-2009 –> 2010-2019	8.414817	65975.3628	66055.2938
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	33769.91191	33816.67780	4.2506648	2000-2009 –> 2010-2019	8.584691	66441.2722	66534.5251
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	33770.58232	33810.66737	4.2516286	2000-2009 –> 2010-2019	8.628869	66552.2337	66632.1648
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	28866.89293	28913.40217	3.0718224	2000-2009 –> 2010-2019	6.257277	56976.0092	57068.7421
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	29128.41637	29174.92562	3.1433986	2000-2009 –> 2010-2019	6.397926	57471.3577	57564.0905
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate	33504.10700	33550.87289	4.1558105	2000-2009 –> 2010-2019	8.350685	65806.1200	65899.3729
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	33717.51158	33764.27747	4.2317956	2000-2009 –> 2010-2019	8.478843	66159.4203	66252.6732
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	29147.84050	29194.34975	3.1487809	2000-2009 –> 2010-2019	6.409355	57511.0159	57603.7487
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	28756.56923	28803.07848	3.0421189	2000-2009 –> 2010-2019	6.093248	56396.0760	56488.8088
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	33660.69136	33707.45725	4.2114294	2000-2009 –> 2010-2019	8.463209	66115.2044	66208.4573
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	33664.46750	33711.23339	4.2127798	2000-2009 –> 2010-2019	8.450481	66081.4440	66174.6969
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + aou + nitrate + phosphate	28868.97652	28908.84159	3.0729274	2000-2009 –> 2010-2019	6.323117	57195.3822	57274.8675
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	28547.89345	28594.40269	2.9867183	2000-2009 –> 2010-2019	6.119220	56474.2906	56567.0234
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	28830.17317	28876.68242	3.0619039	2000-2009 –> 2010-2019	6.282542	57059.0191	57151.7519
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	28831.80777	28878.31702	3.0623447	2000-2009 –> 2010-2019	6.284098	57064.4274	57157.1603
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	28830.41487	28870.27994	3.0625085	2000-2009 –> 2010-2019	6.285214	57064.2572	57143.7425
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	38435.63019	38476.62222	3.9976473	2000-2009 –> 2010-2019	7.667429	72658.9777	72740.4509
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	38965.40719	39013.23122	4.1546570	2000-2009 –> 2010-2019	7.795512	73091.1841	73186.2362
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	32780.60963	32827.91411	3.1805457	2000-2009 –> 2010-2019	6.093703	64941.8443	65036.5624
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	38202.23957	38250.06360	3.9295472	2000-2009 –> 2010-2019	7.768041	72993.1212	73088.1734
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	32954.72313	33002.02762	3.2243809	2000-2009 –> 2010-2019	6.137889	65117.5122	65212.2304
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + silicate + phosphate_star	33220.13798	33260.68468	3.2928852	2000-2009 –> 2010-2019	6.333013	65930.5643	66011.7513
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	32907.91230	32955.21678	3.2125367	2000-2009 –> 2010-2019	6.121713	65051.4777	65146.1959
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	33147.13682	33194.44130	3.2735264	2000-2009 –> 2010-2019	6.277554	65705.2276	65799.9458
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	37989.68778	38037.51181	3.8690519	2000-2009 –> 2010-2019	7.453097	71917.5944	72012.6465
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	33207.02546	33254.32995	3.2889753	2000-2009 –> 2010-2019	6.306820	65824.4052	65919.1234
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	38087.65327	38135.47730	3.8968177	2000-2009 –> 2010-2019	7.513738	72130.4439	72225.4960
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate	33211.58819	33258.89267	3.2901552	2000-2009 –> 2010-2019	6.295793	65776.5997	65871.3179
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	32896.10128	32943.40576	3.2095551	2000-2009 –> 2010-2019	6.172604	65276.7336	65371.4518
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	38052.98049	38100.80452	3.8869679	2000-2009 –> 2010-2019	7.498305	72076.3350	72171.3871
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + silicate + phosphate	33247.82575	33288.37246	3.3000607	2000-2009 –> 2010-2019	6.336494	65942.5429	66023.7299
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	33177.98745	33225.29193	3.2814756	2000-2009 –> 2010-2019	6.291830	65763.2598	65857.9780
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	38073.43332	38121.25735	3.8927751	2000-2009 –> 2010-2019	7.508286	72111.3230	72206.3752
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + silicate + phosphate_star	38364.34272	38405.33475	3.9768997	2000-2009 –> 2010-2019	7.602054	72433.7463	72515.2196
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	38509.57276	38557.39679	4.0186953	2000-2009 –> 2010-2019	7.734437	72891.5164	72986.5686
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ temp + aou + silicate + phosphate	38510.41237	38551.40439	4.0195283	2000-2009 –> 2010-2019	7.742554	72914.9983	72996.4716
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	11491.27921	11532.00541	2.4360386	2000-2009 –> 2010-2019	4.669154	22616.9604	22698.4549
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	11532.29321	11573.01940	2.4562247	2000-2009 –> 2010-2019	4.697266	22675.6889	22757.1834
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + aou + phosphate	11821.09773	11856.00590	2.6042311	2000-2009 –> 2010-2019	4.886017	23052.5815	23122.4339
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + aou + silicate + phosphate	11575.96656	11616.69275	2.4779036	2000-2009 –> 2010-2019	4.759282	22808.5576	22890.0521
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + aou + silicate + phosphate	21246.71314	21290.44368	3.9055425	2000-2009 –> 2010-2019	6.996004	37644.6831	37730.9449
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	11616.48064	11657.20684	2.4981854	2000-2009 –> 2010-2019	4.786969	22865.2751	22946.7697
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	21225.63055	21269.36109	3.8947716	2000-2009 –> 2010-2019	6.986736	37626.7280	37712.9898
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	11557.67601	11598.40220	2.4688012	2000-2009 –> 2010-2019	4.715449	22713.5642	22795.0588
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	21223.46950	21267.20004	3.8936692	2000-2009 –> 2010-2019	6.994398	37642.7724	37729.0342
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	11633.43967	11674.16587	2.5067244	2000-2009 –> 2010-2019	4.797736	22887.0990	22968.5935
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	21225.07693	21268.80747	3.8944891	2000-2009 –> 2010-2019	6.986858	37627.0152	37713.2770
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	11484.72550	11525.45170	2.4328284	2000-2009 –> 2010-2019	4.830016	22964.8989	23046.3934
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	21251.61315	21295.34369	3.9080501	2000-2009 –> 2010-2019	6.998319	37649.1826	37735.4444
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	11535.64302	11576.36922	2.4578808	2000-2009 –> 2010-2019	4.863755	23033.9025	23115.3970
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	21233.17337	21276.90391	3.8986217	2000-2009 –> 2010-2019	6.990578	37634.2541	37720.5159
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + aou + silicate + phosphate_star	21238.51686	21276.00018	3.9023738	2000-2009 –> 2010-2019	6.994357	37637.6519	37711.5906
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	11569.98681	11610.71300	2.4749241	2000-2009 –> 2010-2019	4.886738	23080.5760	23162.0705
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	21233.39086	21277.12140	3.8987328	2000-2009 –> 2010-2019	6.991078	37635.2799	37721.5416
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	21231.59499	21269.07831	3.8988370	2000-2009 –> 2010-2019	6.999569	37648.9044	37722.8431
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	21237.55148	21275.03480	3.9018803	2000-2009 –> 2010-2019	6.994263	37637.5193	37711.4580
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	27292.02057	27337.93832	3.3048834	2000-2009 –> 2010-2019	6.048831	48694.8049	48785.4545
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	13688.73191	13731.16894	2.0874014	2000-2009 –> 2010-2019	4.104860	26380.7364	26465.1923
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	27318.23923	27364.15697	3.3131983	2000-2009 –> 2010-2019	6.058267	48724.6203	48815.2699
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou + nitrate	13072.05233	13108.42692	1.8946971	2000-2009 –> 2010-2019	3.846764	25565.0834	25637.4742
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	13073.61700	13116.05403	1.8945672	2000-2009 –> 2010-2019	3.841748	25553.6666	25638.1225
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	13074.03071	13116.46774	1.8946907	2000-2009 –> 2010-2019	3.843246	25558.2991	25642.7550
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou + nitrate + silicate	13056.34126	13098.77829	1.8894166	2000-2009 –> 2010-2019	3.809926	25453.9409	25538.3968
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou + phosphate	13704.68510	13741.05970	2.0933151	2000-2009 –> 2010-2019	4.146631	26500.0070	26572.3977
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate + phosphate	13554.22697	13596.66400	2.0436240	2000-2009 –> 2010-2019	4.090729	26333.4397	26417.8956
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate + phosphate	27272.24898	27318.16672	3.2986268	2000-2009 –> 2010-2019	6.047104	48689.5624	48780.2121
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	13572.37960	13614.81663	2.0494781	2000-2009 –> 2010-2019	4.102114	26367.7209	26452.1768
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	27263.13718	27309.05492	3.2957474	2000-2009 –> 2010-2019	6.055326	48715.9559	48806.6055
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	13118.58384	13161.02087	1.9080395	2000-2009 –> 2010-2019	3.865683	25630.6674	25715.1233
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	13588.61105	13631.04808	2.0547269	2000-2009 –> 2010-2019	4.109578	26390.4007	26474.8566
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	27263.02165	27308.93940	3.2957109	2000-2009 –> 2010-2019	6.056963	48721.1796	48811.8292
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	27255.09905	27301.01680	3.2932094	2000-2009 –> 2010-2019	6.041364	48671.3773	48762.0269
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	27251.88143	27297.79918	3.2921940	2000-2009 –> 2010-2019	6.051900	48705.1041	48795.7537
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + aou + silicate + phosphate	27277.49825	27316.85632	3.3009194	2000-2009 –> 2010-2019	6.050066	48694.9551	48772.6548
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + aou + silicate + phosphate_star	27270.47676	27309.83483	3.2986988	2000-2009 –> 2010-2019	6.058855	48723.1365	48800.8362
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	27251.88962	27297.80737	3.2921966	2000-2009 –> 2010-2019	6.053173	48709.1686	48799.8182
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	46537.06810	46579.54266	3.4334785	2000-2009 –> 2010-2019	6.284784	83243.9283	83327.8890
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	46537.47267	46587.02632	3.4331663	2000-2009 –> 2010-2019	6.272079	83181.5335	83279.4876
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ aou + phosphate + phosphate_star	46693.91112	46729.30659	3.4647135	2000-2009 –> 2010-2019	6.319946	83419.2433	83489.2105
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	46693.39571	46735.87026	3.4642167	2000-2009 –> 2010-2019	6.312200	83382.9117	83466.8724
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	46539.04727	46588.60091	3.4334745	2000-2009 –> 2010-2019	6.284752	83247.7637	83345.7179
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + aou + phosphate + phosphate_star	46681.43663	46723.91118	3.4618555	2000-2009 –> 2010-2019	6.316934	83407.9690	83491.9297
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + aou + nitrate	17701.01498	17739.47530	1.7332668	2000-2009 –> 2010-2019	3.459052	30490.9084	30565.8909
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + aou + nitrate + silicate	17576.36399	17621.23437	1.7090036	2000-2009 –> 2010-2019	3.387947	30189.2835	30276.7630
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + aou + phosphate	18003.55024	18042.01056	1.7926282	2000-2009 –> 2010-2019	3.619417	31163.3756	31238.3581
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + aou + phosphate_star	18006.17960	18044.63992	1.7931529	2000-2009 –> 2010-2019	3.624094	31180.7702	31255.7527
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate + phosphate	17999.38164	18044.25201	1.7913978	2000-2009 –> 2010-2019	3.588239	31053.6972	31141.1768
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate + phosphate	46777.24132	46826.79497	3.4804193	2000-2009 –> 2010-2019	6.338208	83519.8894	83617.8436
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	18002.93030	18047.80067	1.7921055	2000-2009 –> 2010-2019	3.594777	31078.3183	31165.7979
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	18226.53825	18271.40862	1.8372700	2000-2009 –> 2010-2019	3.678724	31440.3667	31527.8462
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + phosphate + phosphate_star	18012.84184	18051.30216	1.7944832	2000-2009 –> 2010-2019	3.627442	31194.5978	31269.5803
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	18009.37981	18054.25018	1.7933925	2000-2009 –> 2010-2019	3.597597	31090.2954	31177.7749
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + silicate + phosphate_star	18284.60205	18323.06237	1.8495945	2000-2009 –> 2010-2019	3.708159	31556.5871	31631.5696
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	46625.75681	46675.31046	3.4504900	2000-2009 –> 2010-2019	6.302570	83338.6583	83436.6125
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	46740.08688	46789.64052	3.4730546	2000-2009 –> 2010-2019	6.353186	83598.5857	83696.5398
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	46724.66873	46774.22238	3.4700030	2000-2009 –> 2010-2019	6.348280	83573.5859	83671.5400
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	35850.76198	35892.01706	2.9588994	2000-2009 –> 2010-2019	5.687202	58712.8968	58792.9042
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ aou + phosphate + phosphate_star	35924.96731	35959.34654	2.9746941	2000-2009 –> 2010-2019	5.704656	58790.8362	58857.5090
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	35901.75719	35943.01227	2.9694596	2000-2009 –> 2010-2019	5.698443	58766.2464	58846.2538
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	35866.67300	35914.80392	2.9617764	2000-2009 –> 2010-2019	5.702329	58773.0624	58866.4043
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	35865.83563	35913.96655	2.9616031	2000-2009 –> 2010-2019	5.716244	58820.5881	58913.9301
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + aou + nitrate + silicate	9288.58627	9329.58864	1.4549700	2000-2009 –> 2010-2019	3.126513	16735.6943	16815.6282
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + aou + silicate + phosphate	9477.25346	9518.25582	1.5090465	2000-2009 –> 2010-2019	3.087488	16703.9507	16783.8845
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + aou + silicate + phosphate	35864.90307	35913.03400	2.9614102	2000-2009 –> 2010-2019	5.674163	58675.1260	58768.4680
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	9482.52855	9523.53092	1.5105870	2000-2009 –> 2010-2019	3.093449	16719.9826	16799.9164
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	35985.04204	36033.17296	2.9863702	2000-2009 –> 2010-2019	5.721935	58874.2494	58967.5914
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	9432.24841	9473.25078	1.4959672	2000-2009 –> 2010-2019	3.078808	16669.6513	16749.5852
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	35954.33180	36002.46272	2.9799699	2000-2009 –> 2010-2019	5.709676	58823.3174	58916.6593
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	9481.89149	9522.89386	1.5104009	2000-2009 –> 2010-2019	3.093263	16719.3451	16799.2790
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	35986.51489	36034.64581	2.9866775	2000-2009 –> 2010-2019	5.721879	58874.4710	58967.8129
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + silicate + phosphate_star	9487.36224	9522.50713	1.5125850	2000-2009 –> 2010-2019	3.095448	16722.8198	16791.3345
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	36047.39886	36095.52978	2.9994083	2000-2009 –> 2010-2019	5.721941	58891.5644	58984.9063
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + aou + silicate + phosphate	9885.75293	9920.89781	1.6337505	2000-2009 –> 2010-2019	3.236962	17170.3356	17238.8503
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + aou + silicate + phosphate_star	9893.99673	9929.14162	1.6363577	2000-2009 –> 2010-2019	3.245491	17192.7604	17261.2751
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	9883.71077	9918.85565	1.6331053	2000-2009 –> 2010-2019	3.241775	17181.3655	17249.8803
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	9891.92137	9927.06625	1.6357009	2000-2009 –> 2010-2019	3.244811	17190.6291	17259.1438
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	18970.42873	19016.79598	1.3295718	2000-2009 –> 2010-2019	2.771310	36446.7944	36538.6291
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	37478.34600	37527.00880	2.7370489	2000-2009 –> 2010-2019	5.271194	67007.3021	67103.1841
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + aou + phosphate + phosphate_star	37543.43565	37585.14663	2.7489647	2000-2009 –> 2010-2019	5.288240	67095.8041	67177.9886
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + aou + silicate + phosphate	37512.76129	37554.47226	2.7435102	2000-2009 –> 2010-2019	5.279514	67048.9311	67131.1157
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + nitrate + silicate + phosphate	19061.26194	19101.00530	1.3407119	2000-2009 –> 2010-2019	2.803621	36678.2535	36756.9690
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	19063.09964	19109.46689	1.3406924	2000-2009 –> 2010-2019	2.797510	36641.2656	36733.1004
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	37593.06657	37641.72937	2.7574558	2000-2009 –> 2010-2019	5.326280	67292.8183	67388.7003
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	37637.22128	37685.88408	2.7653508	2000-2009 –> 2010-2019	5.339819	67364.5529	67460.4350
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + aou + nitrate + silicate	18049.03500	18095.40225	1.2239000	2000-2009 –> 2010-2019	2.746322	36058.1675	36150.0022
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + aou + silicate + phosphate	17637.34583	17683.71308	1.1794404	2000-2009 –> 2010-2019	2.554576	34650.9568	34742.7915
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + aou + silicate + phosphate	37497.64075	37546.30355	2.7404705	2000-2009 –> 2010-2019	5.257763	66942.7487	67038.6307
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	17656.52571	17702.89296	1.1814753	2000-2009 –> 2010-2019	2.558507	34683.6152	34775.4499
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	37599.66180	37648.32460	2.7586336	2000-2009 –> 2010-2019	5.295922	67144.1932	67240.0752
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate	18509.74995	18556.11720	1.2756443	2000-2009 –> 2010-2019	2.714623	35967.3696	36059.2044
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	17614.77115	17661.13840	1.1770497	2000-2009 –> 2010-2019	2.564400	34714.9026	34806.7374
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	37676.28833	37724.95113	2.7723548	2000-2009 –> 2010-2019	5.326921	67306.0995	67401.9815
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	17652.98941	17699.35666	1.1810999	2000-2009 –> 2010-2019	2.557863	34678.1754	34770.0102
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	37583.35403	37632.01683	2.7557223	2000-2009 –> 2010-2019	5.288981	67107.9155	67203.7975
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + silicate + phosphate_star	17827.69114	17867.43450	1.2000077	2000-2009 –> 2010-2019	2.590476	34947.7900	35026.5055
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + silicate + phosphate_star	37706.63652	37748.34749	2.7781677	2000-2009 –> 2010-2019	5.340727	67373.7033	67455.8879
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	11305.36811	11340.67330	2.0296691	2000-2009 –> 2010-2019	4.639732	19555.8438	19623.8912
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	11224.59929	11265.78869	1.9982773	2000-2009 –> 2010-2019	4.606358	19474.4431	19553.8318
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	11247.18979	11288.37919	2.0067967	2000-2009 –> 2010-2019	4.614359	19496.3456	19575.7342
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	31794.67076	31847.99725	0.6962973	2000-2009 –> 2010-2019	1.350573	54010.6341	54115.2012
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	34927.28215	34980.60863	0.7727541	2000-2009 –> 2010-2019	1.542066	60758.3590	60862.9261
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + aou + phosphate + phosphate_star	11378.91591	11414.22111	2.0579773	2000-2009 –> 2010-2019	4.667176	19628.2440	19696.2914
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	11368.97598	11410.16538	2.0533550	2000-2009 –> 2010-2019	4.662418	19620.1246	19699.5132
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + nitrate + silicate + phosphate	32251.20979	32296.91821	0.7069972	2000-2009 –> 2010-2019	1.453878	57419.8276	57509.4566
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	32071.80334	32125.12983	0.7027446	2000-2009 –> 2010-2019	1.362657	54479.2528	54583.8199
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	11229.50383	11270.69323	2.0001238	2000-2009 –> 2010-2019	4.658297	19545.2485	19624.6371
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	11361.88753	11403.07692	2.0506157	2000-2009 –> 2010-2019	4.658849	19611.9338	19691.3224
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + temp + aou + silicate + phosphate	34853.64561	34906.97210	0.7708639	2000-2009 –> 2010-2019	1.601766	62403.7304	62508.2975
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	34865.84468	34919.17116	0.7711768	2000-2009 –> 2010-2019	1.602382	62424.0553	62528.6224
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate	32062.42355	32115.75004	0.7025254	2000-2009 –> 2010-2019	1.441564	56997.4371	57102.0042
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	34863.47379	34916.80027	0.7711160	2000-2009 –> 2010-2019	1.602283	62420.6618	62525.2289
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star	11261.34733	11302.53673	2.0121544	2000-2009 –> 2010-2019	4.622142	19513.7231	19593.1117
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	36006.70565	36060.03214	0.8009995	2000-2009 –> 2010-2019	1.638492	63733.1336	63837.7007
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	11232.75012	11273.93952	2.0013470	2000-2009 –> 2010-2019	4.649311	19535.1658	19614.5545
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	11375.19307	11416.38247	2.0557605	2000-2009 –> 2010-2019	4.665310	19626.9875	19706.3761
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	36006.41096	36059.73745	0.8009916	2000-2009 –> 2010-2019	1.639735	63766.1279	63870.6950
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	20739.58571	20786.78931	1.2637496	2000-2009 –> 2010-2019	2.190694	33644.4353	33736.9720
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	20408.05818	20455.26178	1.2307716	2000-2009 –> 2010-2019	2.159535	33331.7292	33424.2659
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	4496.86862	4531.55777	2.0499098	2000-2009 –> 2010-2019	4.698153	7774.6955	7841.0598
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	4480.78683	4515.47598	2.0342567	2000-2009 –> 2010-2019	4.668761	7751.5191	7817.8835
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + aou + silicate + phosphate_star	4500.77426	4530.50781	2.0556882	2000-2009 –> 2010-2019	4.706246	7777.7929	7834.6767
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	4528.29766	4562.98681	2.0808496	2000-2009 –> 2010-2019	4.714287	7798.4774	7864.8418
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	4529.51509	4559.24865	2.0840432	2000-2009 –> 2010-2019	4.718953	7798.4570	7855.3407
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	20385.41080	20432.61441	1.2285505	2000-2009 –> 2010-2019	2.154602	33281.0106	33373.5473
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	20467.80967	20515.01328	1.2366510	2000-2009 –> 2010-2019	2.167667	33414.7220	33507.2587
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + aou + silicate + phosphate	4505.76622	4540.45537	2.0586219	2000-2009 –> 2010-2019	4.692150	7775.9927	7842.3571
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	4479.82603	4514.51518	2.0333253	2000-2009 –> 2010-2019	4.683037	7758.4093	7824.7737
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	20383.28744	20430.49105	1.2283425	2000-2009 –> 2010-2019	2.155463	33289.9626	33382.4993
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate	20358.68974	20405.89334	1.2259350	2000-2009 –> 2010-2019	2.192184	33662.1233	33754.6600
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate	4505.40641	4540.09556	2.0582689	2000-2009 –> 2010-2019	4.708270	7784.1383	7850.5027
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	4499.49454	4534.18369	2.0524771	2000-2009 –> 2010-2019	4.686840	7770.1535	7836.5178
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ temp + aou + nitrate + phosphate	20719.82618	20760.28641	1.2619608	2000-2009 –> 2010-2019	2.188553	33619.0309	33698.3481
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	20721.45892	20768.66252	1.2619238	2000-2009 –> 2010-2019	2.185827	33594.7683	33687.3051
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	4477.47920	4512.16835	2.0310521	2000-2009 –> 2010-2019	4.716448	7774.3087	7840.6731
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ temp + nitrate + phosphate + phosphate_star	20722.57078	20763.03101	1.2622371	2000-2009 –> 2010-2019	2.189956	33633.4364	33712.7536
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	20724.16023	20771.36383	1.2621957	2000-2009 –> 2010-2019	2.187285	33609.7804	33702.3172
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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