eMLR - model fitting

workflowr

• Summary
• Checks
• Past versions

Last updated: 2021-03-12

Checks: 7 0

Knit directory: emlr_mod_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: 3ed96fb

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 3ed96fb. 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:
    Deleted:    a.png
    Modified:   data/auxillary/params_local.rds

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

---

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

File	Version	Author	Date	Message
html	9dbf5bd	Donghe-Zhu	2021-03-11	Build site.
html	a49df30	Donghe-Zhu	2021-03-11	Build site.
html	b3348a5	Donghe-Zhu	2021-03-11	Build site.
Rmd	026c061	Donghe-Zhu	2021-03-11	local rebuild
html	1c24ff7	Donghe-Zhu	2021-03-10	Build site.
Rmd	dbd8ff8	Donghe-Zhu	2021-03-10	local rebuild
html	94ce5a8	Donghe-Zhu	2021-03-10	Build site.
html	db33928	Donghe-Zhu	2021-03-10	Build site.
html	3d1f470	Donghe-Zhu	2021-03-10	Build site.
html	9b7bc66	Donghe-Zhu	2021-03-10	Build site.
Rmd	498be1d	Donghe-Zhu	2021-03-10	local rebuild
html	4168b43	Donghe-Zhu	2021-03-10	Build site.
html	5365f80	Donghe-Zhu	2021-03-10	Build site.
html	2b6c392	Donghe-Zhu	2021-03-10	Build site.
html	9f58753	Donghe-Zhu	2021-03-10	Build site.
Rmd	e255177	Donghe-Zhu	2021-03-10	local rebuild
html	f2a7146	Donghe-Zhu	2021-03-10	Build site.
html	fd528ed	Donghe-Zhu	2021-03-10	Build site.
html	9d5a62c	Donghe-Zhu	2021-03-10	Build site.
Rmd	09b27ef	Donghe-Zhu	2021-03-10	local rebuild
html	5d1e70b	Donghe-Zhu	2021-03-10	Build site.
Rmd	6125b2d	Donghe-Zhu	2021-03-10	local rebuild
html	2093979	Donghe-Zhu	2021-03-10	Build site.
Rmd	558d285	Donghe-Zhu	2021-03-10	local rebuild
html	b865899	Donghe-Zhu	2021-03-10	Build site.
html	cc2a956	Donghe-Zhu	2021-03-10	Build site.
Rmd	2a633bd	Donghe-Zhu	2021-03-10	local rebuild
html	60689fb	Donghe-Zhu	2021-03-10	Build site.
html	9230b52	Donghe-Zhu	2021-03-10	Build site.
html	dba33c8	Donghe-Zhu	2021-03-09	Build site.
Rmd	6a60e02	Donghe-Zhu	2021-03-09	local rebuild
html	17f1c4a	Donghe-Zhu	2021-03-09	Build site.
Rmd	4a4a400	Donghe-Zhu	2021-03-09	local rebuild
html	c024d1a	Donghe-Zhu	2021-03-09	Build site.
Rmd	34744ec	Donghe-Zhu	2021-03-09	local rebuild
html	02f7242	Donghe-Zhu	2021-03-09	Build site.
html	6f50bc6	Donghe-Zhu	2021-03-09	Build site.
Rmd	88f0729	Donghe-Zhu	2021-03-09	local rebuild
html	1691156	Donghe-Zhu	2021-03-08	Build site.
html	c0ceaf8	Donghe-Zhu	2021-03-08	Build site.
html	058e0a1	Donghe-Zhu	2021-03-08	Build site.
html	112dea0	Donghe-Zhu	2021-03-08	Build site.
html	1843412	Donghe-Zhu	2021-03-08	Build site.
Rmd	d709f6c	Donghe-Zhu	2021-03-08	local rebuild
html	65b0cef	Donghe-Zhu	2021-03-07	Build site.
html	4083a6c	Donghe-Zhu	2021-03-07	Build site.
html	3fbbfa4	Donghe-Zhu	2021-03-07	Build site.
html	627c8fb	Donghe-Zhu	2021-03-07	Build site.
html	3607f4d	Donghe-Zhu	2021-03-07	Build site.
html	9ef3222	Donghe-Zhu	2021-03-05	Build site.
Rmd	2be54f9	Donghe-Zhu	2021-03-05	local rebuild
html	8c1e978	Donghe-Zhu	2021-03-05	Build site.
html	865f68c	Donghe-Zhu	2021-03-05	Build site.
html	ee69bc1	Donghe-Zhu	2021-03-05	Build site.
html	a79291f	Donghe-Zhu	2021-03-05	Build site.
html	e8c6f30	Donghe-Zhu	2021-03-04	Build site.
html	59288fe	Donghe-Zhu	2021-03-04	Build site.
html	731abc8	Donghe-Zhu	2021-03-04	Build site.
html	e2a5a33	Donghe-Zhu	2021-03-04	Build site.
html	c7892c1	Donghe-Zhu	2021-03-04	Build site.
html	924430b	Donghe-Zhu	2021-03-03	Build site.
html	0d0bca1	Donghe-Zhu	2021-03-03	Build site.
html	cb63c16	Donghe-Zhu	2021-03-03	Build site.
html	ffda45a	Donghe-Zhu	2021-03-03	Build site.
html	691ba81	Donghe-Zhu	2021-03-03	Build site.
html	c5e45a2	Donghe-Zhu	2021-03-03	Build site.
html	89c3e58	Donghe-Zhu	2021-03-03	Build site.
html	c407a50	Donghe-Zhu	2021-03-03	Build site.
html	c911669	Donghe-Zhu	2021-03-03	Build site.
Rmd	e1f9d21	Donghe-Zhu	2021-03-03	local rebuild
html	b71c719	Donghe-Zhu	2021-03-01	Build site.
html	13666ca	Donghe-Zhu	2021-03-01	Build site.
Rmd	c3195f4	Donghe-Zhu	2021-03-01	local rebuild
html	c6e60fe	Donghe-Zhu	2021-03-01	Build site.
Rmd	8ef147d	Donghe-Zhu	2021-03-01	local rebuild
html	7a388f7	Donghe-Zhu	2021-03-01	Build site.
html	799e913	Donghe-Zhu	2021-03-01	Build site.
Rmd	aa780a7	Donghe-Zhu	2021-03-01	local rebuild
html	66ff99f	Donghe-Zhu	2021-03-01	Build site.
html	ac9bb7a	Donghe-Zhu	2021-02-28	Build site.
html	efdc047	Donghe-Zhu	2021-02-28	Build site.
html	e9a7418	Donghe-Zhu	2021-02-28	Build site.
Rmd	eac1b0e	Donghe-Zhu	2021-02-28	local rebuild
html	e152917	Donghe-Zhu	2021-02-28	Build site.
html	feb991c	Donghe-Zhu	2021-02-27	Build site.
Rmd	4563314	Donghe-Zhu	2021-02-27	local rebuild
html	287123c	Donghe-Zhu	2021-02-27	Build site.
Rmd	764a9c2	Donghe-Zhu	2021-02-27	local rebuild
html	54d5b5b	Donghe-Zhu	2021-02-27	Build site.
Rmd	2e836bf	Donghe-Zhu	2021-02-27	local rebuild
html	330f064	Donghe-Zhu	2021-02-27	Build site.
Rmd	f360282	Donghe-Zhu	2021-02-27	local rebuild
html	adbc9bc	Donghe-Zhu	2021-02-27	Build site.
Rmd	a8a8f9c	Donghe-Zhu	2021-02-27	local rebuild
html	5937141	Donghe-Zhu	2021-02-27	Build site.
Rmd	5f0bfca	Donghe-Zhu	2021-02-27	local rebuild
html	4414bbf	Donghe-Zhu	2021-02-27	Build site.
html	a265efb	Donghe-Zhu	2021-02-27	Build site.
html	19edd1e	Donghe-Zhu	2021-02-27	Build site.
Rmd	1ae5bb4	Donghe-Zhu	2021-02-27	local rebuild
html	f20483f	Donghe-Zhu	2021-02-26	Build site.
html	6a2c7b3	Donghe-Zhu	2021-02-25	Build site.
html	02b976d	Donghe-Zhu	2021-02-24	Build site.
html	354c224	Donghe-Zhu	2021-02-24	Build site.
Rmd	d910d39	Donghe-Zhu	2021-02-24	local rebuild
html	1a0a88a	Donghe-Zhu	2021-02-24	Build site.
Rmd	663e891	Donghe-Zhu	2021-02-24	local rebuild
html	57f701e	Donghe-Zhu	2021-02-24	Build site.
html	06f3149	Donghe-Zhu	2021-02-16	Build site.
html	401eab3	Donghe-Zhu	2021-02-15	Build site.
html	e3bba84	Donghe-Zhu	2021-02-15	Build site.
html	5dce4b1	Donghe-Zhu	2021-02-15	Build site.
Rmd	35b0f2e	Donghe-Zhu	2021-02-15	local rebuild
html	4469a0c	Donghe-Zhu	2021-02-13	Build site.
Rmd	8186d57	Donghe-Zhu	2021-02-10	apply nr_obs
html	5ae6a69	Donghe-Zhu	2021-02-10	Build site.
Rmd	eeb6557	Donghe-Zhu	2021-02-10	local rebuild
html	05385dc	Donghe-Zhu	2021-02-10	Build site.
html	f791ae4	Donghe-Zhu	2021-02-09	Build site.
html	f71ae34	Donghe-Zhu	2021-02-09	Build site.
html	c011832	Donghe-Zhu	2021-02-09	Build site.
html	a145fa7	Donghe-Zhu	2021-02-09	Build site.
Rmd	37a41c8	Donghe-Zhu	2021-02-09	local rebuild
html	c344e42	Donghe-Zhu	2021-02-08	Build site.
Rmd	8f5fa79	Donghe-Zhu	2021-02-08	local rebuild
html	2f095d7	Donghe-Zhu	2021-02-07	Build site.
html	2305044	Donghe-Zhu	2021-02-07	Build site.
Rmd	c3638c1	Donghe-Zhu	2021-02-07	local rebuild
html	1fad5f1	Donghe-Zhu	2021-02-07	Build site.
html	ca03c39	Donghe-Zhu	2021-02-07	Build site.
html	e2ffc14	Donghe-Zhu	2021-02-05	Build site.
html	cd7c52c	Donghe-Zhu	2021-02-04	Build site.
Rmd	bcf91a8	jens-daniel-mueller	2021-02-04	separate MLR tables, added plot residuals vs location
html	bcf84f4	Donghe-Zhu	2021-02-02	Build site.
html	a518739	Donghe-Zhu	2021-02-01	Build site.
html	61666de	Donghe-Zhu	2021-01-31	Build site.
html	865b582	Donghe-Zhu	2021-01-31	Build site.
html	3e68089	Donghe-Zhu	2021-01-31	Build site.
html	ecf335c	Donghe-Zhu	2021-01-31	Build site.
html	a618965	Donghe-Zhu	2021-01-31	Build site.
html	59e006e	Donghe-Zhu	2021-01-31	Build site.
html	a1c8f87	Donghe-Zhu	2021-01-31	Build site.
html	ae5c18f	Donghe-Zhu	2021-01-31	Build site.
html	b50fe52	Donghe-Zhu	2021-01-31	Build site.
Rmd	ac99ae5	jens-daniel-mueller	2021-01-29	code review
html	ac99ae5	jens-daniel-mueller	2021-01-29	code review
html	b5bdcaf	Donghe-Zhu	2021-01-29	Build site.
Rmd	b234505	Donghe-Zhu	2021-01-29	MLR approach across all depth
html	442010d	Donghe-Zhu	2021-01-29	Build site.
Rmd	e67e7dd	Donghe-Zhu	2021-01-29	surface equilibrium approach across all latitudes irrespective of gamma
html	372adf5	Donghe-Zhu	2021-01-29	Build site.
html	af8788e	Donghe-Zhu	2021-01-29	Build site.
html	21c91c9	Donghe-Zhu	2021-01-29	Build site.
html	eded038	Donghe-Zhu	2021-01-29	Build site.
html	541d4dd	Donghe-Zhu	2021-01-29	Build site.
html	6a75576	Donghe-Zhu	2021-01-28	Build site.
html	16fba40	Donghe-Zhu	2021-01-28	Build site.
Rmd	aecf0c6	Donghe-Zhu	2021-01-28	diagnostic ploting
Rmd	a02684e	Donghe-Zhu	2021-01-28	error
Rmd	4584be9	Donghe-Zhu	2021-01-28	latest
html	12bc567	Donghe-Zhu	2021-01-27	Build site.
html	ceed31b	Donghe-Zhu	2021-01-27	Build site.
html	342402d	Donghe-Zhu	2021-01-27	Build site.
html	5bad5c2	Donghe-Zhu	2021-01-27	Build site.
Rmd	c2c9529	Donghe-Zhu	2021-01-27	random subsetting based on lat
html	61efb56	Donghe-Zhu	2021-01-25	Build site.
html	48f638e	Donghe-Zhu	2021-01-25	Build site.
html	c1cec47	Donghe-Zhu	2021-01-25	Build site.
html	05ffb0c	Donghe-Zhu	2021-01-25	Build site.
html	8b97165	Donghe-Zhu	2021-01-25	Build site.
html	c569946	Donghe-Zhu	2021-01-24	Build site.
html	a2f0d56	Donghe-Zhu	2021-01-23	Build site.
html	28509fc	Donghe-Zhu	2021-01-23	Build site.
html	4c28e4a	Donghe-Zhu	2021-01-22	Build site.
html	24cc264	jens-daniel-mueller	2021-01-22	cleaned /docs before creating copies
html	88eb28f	Donghe-Zhu	2021-01-21	Build site.
html	2679490	Donghe-Zhu	2021-01-21	Build site.
html	7891955	Donghe-Zhu	2021-01-21	Build site.
html	d4cf1cb	Donghe-Zhu	2021-01-21	Build site.
Rmd	167eeec	Donghe-Zhu	2021-01-21	surface DIC calculation with atmospheric equilibrium option
html	1f3e5b6	jens-daniel-mueller	2021-01-20	Build site.
html	0e7bdf1	jens-daniel-mueller	2021-01-15	cleaning template repository
html	73cbef3	jens-daniel-mueller	2021-01-15	Build site.
html	4571843	jens-daniel-mueller	2021-01-14	revision and html deleted for template copying
html	23151cd	jens-daniel-mueller	2021-01-14	Build site.
html	b3564aa	jens-daniel-mueller	2021-01-14	Build site.
html	8d032c3	jens-daniel-mueller	2021-01-14	Build site.
html	022871c	Donghe-Zhu	2021-01-13	Build site.
Rmd	d44f36f	Donghe-Zhu	2021-01-13	reorder analysis final
html	17dee1d	jens-daniel-mueller	2021-01-13	Build site.
Rmd	9e04fd7	jens-daniel-mueller	2021-01-13	local rebuild after revision
html	a076226	Donghe-Zhu	2021-01-11	Build site.
Rmd	52eff18	Donghe-Zhu	2021-01-09	Implemet model_run and subsetting
html	7cdea0c	jens-daniel-mueller	2021-01-06	Build site.
Rmd	b5934dd	jens-daniel-mueller	2021-01-06	local rebuild after revision
html	fa85b93	jens-daniel-mueller	2021-01-06	Build site.
html	e5cb81a	Donghe-Zhu	2021-01-05	Build site.
Rmd	608cc45	Donghe-Zhu	2021-01-05	modification of analysis
html	a499f10	Donghe-Zhu	2021-01-05	Build site.
Rmd	715bdb4	Donghe-Zhu	2021-01-02	model modification
html	fb8a752	Donghe-Zhu	2020-12-23	Build site.
Rmd	82e3c9c	Donghe-Zhu	2020-12-23	first build after creating model template
html	8fae0b2	Donghe-Zhu	2020-12-21	Build site.
Rmd	00a1322	Donghe-Zhu	2020-12-21	first build after creating model template
Rmd	d73ae35	Donghe-Zhu	2020-12-21	first version with lm error
html	c8b76b3	jens-daniel-mueller	2020-12-19	Build site.
Rmd	b5fedce	jens-daniel-mueller	2020-12-19	first build after creating model template
Rmd	8e8abf5	Jens Müller	2020-12-18	Initial commit

---

1 Required data

Required are:

• cleaned and prepared GLODAP-based synthetic model subsetting file

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

2 Predictor combinations

Find all possible combinations of following considered predictor variables:

• sal, temp, oxygen, aou, 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)

i_gamma_slab <- unique(GLODAP_basin_era$gamma_slab)[1]
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)

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

# remove certain predictor combinations
# lm_rm_ph <- lm_all %>%
#   filter(str_detect(lm_coeff, "phosphate_star")) %>%
#   mutate(lm_coeff_filter = str_remove(lm_coeff, "phosphate_star")) %>%
#   filter(
#     str_detect(lm_coeff_filter, "oxygen") &
#       str_detect(lm_coeff_filter, "phosphate")
#   )

# lm_rm_si <- lm_all %>%
#   filter(str_detect(lm_coeff, "silicate_star")) %>%
#   mutate(lm_coeff_filter = str_remove(lm_coeff, "silicate_star")) %>%
#   filter(str_detect(lm_coeff_filter, "silicate"))

# lm_rm_o2 <- lm_all %>%
#   filter(str_detect(lm_coeff, "phosphate_star")) %>%
#   mutate(lm_coeff_filter = str_remove(lm_coeff, "phosphate_star")) %>%
#   filter(
#     str_detect(lm_coeff_filter, "phosphate") &
#       str_detect(lm_coeff_filter, "oxygen")
#   )

# lm_rm <- bind_rows(lm_rm_ph, lm_rm_o2) %>%
#   select(-lm_coeff_filter) %>%
#   unique()

# lm_rm <- lm_rm_ph %>%
#   select(-lm_coeff_filter) %>%
#   unique()
# 
# lm_all <- anti_join(lm_all, lm_rm)

# remove temp sal predictor combination
# lm_all <- lm_all %>%
# # filter(!(
# #   str_detect(lm_coeff, "temp") &
# #     str_detect(lm_coeff, "phosphate_star")
# # )) %>%
# mutate(lm_coeff_filter = str_remove(lm_coeff, "phosphate_star")) %>%
# # filter(!(str_detect(lm_coeff_filter, "nitrate") &
# #            str_detect(lm_coeff_filter, "phosphate")
# # )) %>%
# filter(!(
#   str_detect(lm_coeff_filter, "temp") &
#     str_detect(lm_coeff_filter, "sal")
# )) %>%
# filter(!(
#   str_detect(lm_coeff_filter, "oxygen") &
#     str_detect(lm_coeff_filter, "aou")
# )) %>%
# select(-lm_coeff_filter)

# remove helper objects
rm(
  i_gamma_slab,
  i_era,
  i_basin,
  GLODAP_basin_era,
  GLODAP_basin_era_slab,
  lm_full,
  lm_rm_ph,
  lm_rm_si,
  lm_rm_o2,
  lm_rm
)

# lm_full <- lm(paste(
#   params_local$MLR_target,
#   paste(params_local$MLR_predictors, collapse = " + "),
#   sep = " ~ "
# ),
# data = GLODAP_basin_era_slab)
# 
# lm_full <-
#   loess(
#     cstar_tref ~ sal + temp + aou + phosphate,
#     span = 0.1,
#     degree = 1,
#     data = GLODAP_basin_era_slab
#   )

3 Apply predictor threshold

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

• Minimum: 2
  
• Maximum: 5

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

This results in a total number of MLR models of:

• 112

4 Fit all models

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

• cstar_tref

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

# loop across all basins, era, gamma slabs, and MLRs
# fit all MLR models
for (i_basin in unique(GLODAP$basin)) {
  for (i_era in unique(GLODAP$era)) {
    #i_basin <- unique(GLODAP$basin)[1]
    #i_era   <- unique(GLODAP$era)[1]
    print(i_basin)
    print(i_era)
    
    GLODAP_basin_era <- GLODAP %>%
      filter(basin == i_basin, era == i_era)
    
    for (i_gamma_slab in unique(GLODAP_basin_era$gamma_slab)) {
      #i_gamma_slab <- unique(GLODAP_basin_era$gamma_slab)[1]
      print(i_gamma_slab)
      
      GLODAP_basin_era_slab <- GLODAP_basin_era %>%
        filter(gamma_slab == i_gamma_slab)
      
      # number of observations used for each fitting model
      i_nr_obs = nrow(GLODAP_basin_era_slab)
      
      for (i_predictors in unique(lm_all$predictors)) {
        #i_predictors <- unique(lm_all$predictors)[1]
        
        # extract one model definition
        i_lm <- lm_all %>%
          filter(predictors == i_predictors) %>%
          select(lm_coeff) %>%
          pull()
        
        # extract number of predictors
        i_n_predictors <- lm_all %>%
          filter(predictors == i_predictors) %>%
          select(n) %>%
          pull()
        
        if (i_nr_obs > i_n_predictors) {
          # fit model
          if (params_local$MLR_type == "rlm") {
            i_lm_fit <- MASS::rlm(as.formula(i_lm),
                                  data = GLODAP_basin_era_slab)
          }
          
          if (params_local$MLR_type == "lm") {
            i_lm_fit <- lm(as.formula(i_lm),
                           data = GLODAP_basin_era_slab)
          }
          
          # find max predictor correlation
          i_cor_max <- GLODAP_basin_era_slab %>%
            select(!!!syms(str_split(i_predictors, " ",
                                     simplify = TRUE))) %>%
            correlate(quiet = TRUE) %>%
            select(-term) %>%
            abs() %>%
            max(na.rm = TRUE)
          
          # calculate root mean squared error
          i_rmse <- sqrt(c(crossprod(i_lm_fit$residuals)) /
                           length(i_lm_fit$residuals))
          
          # calculate maximum residual
          i_resid_max <- max(abs(i_lm_fit$residuals))
          
          # calculate Akaike information criterion aic
          i_aic <- AIC(i_lm_fit)
          
          # collect model coefficients and diagnostics
          coefficients <- tidy(i_lm_fit)
          
          coefficients <- coefficients %>%
            mutate(
              basin = i_basin,
              era = i_era,
              gamma_slab = i_gamma_slab,
              model = i_lm,
              nr_obs = i_nr_obs,
              rmse = i_rmse,
              aic = i_aic,
              resid_max = i_resid_max,
              n_predictors = i_n_predictors,
              na_predictor = anyNA(coefficients$estimate),
              cor_max = i_cor_max
            )
          
          if (exists("lm_all_fitted")) {
            lm_all_fitted <- bind_rows(lm_all_fitted, coefficients)
          }
          
          if (!exists("lm_all_fitted")) {
            lm_all_fitted <- coefficients
          }
        }
      }
    }
  }
}

rm(
  i_lm_fit,
  coefficients,
  i_rmse,
  GLODAP_basin_era,
  GLODAP_basin_era_slab,
  i_lm,
  i_basin,
  i_era,
  i_gamma_slab,
  i_nr_obs,
  i_predictors,
  i_aic,
  i_n_predictors,
  i_resid_max
)

5 Prepare coeffcients

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

5.1 Formatting

# select relevant columns
lm_all_fitted <- lm_all_fitted %>%
  select(
    basin,
    gamma_slab,
    era,
    model,
    nr_obs,
    n_predictors,
    term,
    estimate,
    rmse,
    aic,
    resid_max,
    na_predictor,
    cor_max
  )

# set coefficient to zero if not fitted (=NA)
lm_all_fitted <- lm_all_fitted %>%
  mutate(estimate = if_else(is.na(estimate), 0, estimate))

# Prepare model coefficients for mapping of target variable
lm_all_fitted_wide <- lm_all_fitted %>%
  pivot_wider(
    values_from = estimate,
    names_from = term,
    names_prefix = "coeff_",
    values_fill = 0
  )

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

Please note, that currently the lm() function produces NAs for some predictors. It is not yet entirely clear when this happens, but presumably it is caused by some form of collinearity between predictors, such that including another predictor does not help to explain the target variable any better. The issues also expresses as exactly identical rmse values of different models. As an interim solution, models with fitted NA predictors were not included.

# remove models with predictors fitted as NA
lm_all_fitted_wide <- lm_all_fitted_wide %>%
  filter(na_predictor == FALSE)

# calculate RMSE sum for adjacent eras
lm_all_fitted_wide_eras <- lm_all_fitted_wide %>%
  select(basin, gamma_slab, n_predictors, model, era, nr_obs, rmse, aic, resid_max) %>%
  arrange(era) %>%
  group_by(basin, gamma_slab, model) %>%
  mutate(
    eras = paste(lag(era), era, sep = " --> "),
    rmse_sum = rmse + lag(rmse),
    aic_sum = aic + lag(aic)
  ) %>%
  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 <- lm_all_fitted_wide_eras %>%
    group_by(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)
} else {
  lm_best <- lm_all_fitted_wide_eras %>%
    group_by(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)
}

5.3 RMSE tables

5.3.1 per model

lm_best %>%
  kable() %>%
  add_header_above() %>%
  kable_styling() %>%
  scroll_box(width = "100%", height = "400px")

basin	gamma_slab	n_predictors	model	nr_obs	rmse	aic	resid_max	eras	rmse_sum	aic_sum
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	37	0.6421586	86.22538	2.3091622	1982-1999 –> 2000-2010	2.1573266	291.0097
Atlantic	(-Inf,26]	4	cstar_tref ~ sal + temp + oxygen + aou	37	0.6935708	89.92471	2.1803209	1982-1999 –> 2000-2010	2.3308825	300.7721
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	37	0.6421586	86.22538	2.3091622	1982-1999 –> 2000-2010	2.1573266	291.0097
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	37	0.6421586	86.22538	2.3091622	1982-1999 –> 2000-2010	2.1573266	291.0097
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + aou + silicate	37	0.6455841	86.61907	1.4926622	1982-1999 –> 2000-2010	2.2827871	299.4596
Atlantic	(-Inf,26]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	37	0.6855656	89.06563	2.0085839	1982-1999 –> 2000-2010	2.3503343	301.6426
Atlantic	(-Inf,26]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	37	0.5798011	78.66626	1.3702897	1982-1999 –> 2000-2010	1.9462089	272.7031
Atlantic	(-Inf,26]	4	cstar_tref ~ temp + oxygen + aou + phosphate	37	0.6855656	89.06563	2.0085839	1982-1999 –> 2000-2010	2.3503343	301.6426
Atlantic	(-Inf,26]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	37	0.5798011	78.66626	1.3702897	1982-1999 –> 2000-2010	1.9462089	272.7031
Atlantic	(-Inf,26]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	37	0.5798011	78.66626	1.3702897	1982-1999 –> 2000-2010	1.9462089	272.7031
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	40	0.4807074	68.91537	1.1464555	2000-2010 –> 2011-2019	1.0603777	147.5649
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	40	0.4857089	69.74342	1.1641047	2000-2010 –> 2011-2019	1.0670403	148.6047
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + aou + silicate	40	0.4110883	56.39931	0.7634964	2000-2010 –> 2011-2019	1.0566725	143.0184
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	40	0.4806721	68.90948	1.1457847	2000-2010 –> 2011-2019	1.0599812	147.5129
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	40	0.4806721	68.90948	1.1457847	2000-2010 –> 2011-2019	1.0599812	147.5129
Atlantic	(-Inf,26]	4	cstar_tref ~ sal + temp + silicate + phosphate	40	0.4807134	66.91635	1.1463628	2000-2010 –> 2011-2019	1.0648497	144.1339
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	40	0.4806721	68.90948	1.1457847	2000-2010 –> 2011-2019	1.0599812	147.5129
Atlantic	(-Inf,26]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	40	0.3533587	44.29335	0.8537963	2000-2010 –> 2011-2019	0.9331598	122.9596
Atlantic	(-Inf,26]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	40	0.3533587	44.29335	0.8537963	2000-2010 –> 2011-2019	0.9331598	122.9596
Atlantic	(-Inf,26]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	40	0.3533587	44.29335	0.8537963	2000-2010 –> 2011-2019	0.9331598	122.9596
Atlantic	(26,26.5]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	277	3.8308935	1544.16828	10.6008456	1982-1999 –> 2000-2010	6.9139272	3293.7578
Atlantic	(26,26.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	277	3.8308935	1544.16828	10.6008456	1982-1999 –> 2000-2010	6.9139272	3293.7578
Atlantic	(26,26.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	277	3.8308935	1544.16828	10.6008456	1982-1999 –> 2000-2010	6.9139272	3293.7578
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	277	3.9272741	1557.93379	11.4491971	1982-1999 –> 2000-2010	7.1128176	3329.8307
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	277	3.9410826	1559.87828	11.6034836	1982-1999 –> 2000-2010	7.1458710	3335.8830
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + aou + silicate	277	4.2143882	1597.02340	13.2878840	1982-1999 –> 2000-2010	8.0775467	3500.4523
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	277	3.9459887	1560.56749	11.6156157	1982-1999 –> 2000-2010	7.1583010	3338.1714
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	277	3.9459887	1560.56749	11.6156157	1982-1999 –> 2000-2010	7.1583010	3338.1714
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	277	3.9459887	1560.56749	11.6156157	1982-1999 –> 2000-2010	7.1583010	3338.1714
Atlantic	(26,26.5]	4	cstar_tref ~ sal + temp + silicate + phosphate_star	277	4.0858296	1577.86069	10.6904513	1982-1999 –> 2000-2010	7.7155986	3436.7900
Atlantic	(26,26.5]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	277	3.5625544	1503.93666	10.5402083	2000-2010 –> 2011-2019	7.3934479	3048.1049
Atlantic	(26,26.5]	4	cstar_tref ~ sal + oxygen + aou + silicate	277	3.6836996	1520.46229	11.9948483	2000-2010 –> 2011-2019	7.9343341	3120.2301
Atlantic	(26,26.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	277	3.5625544	1503.93666	10.5402083	2000-2010 –> 2011-2019	7.3934479	3048.1049
Atlantic	(26,26.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	277	3.5625544	1503.93666	10.5402083	2000-2010 –> 2011-2019	7.3934479	3048.1049
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	277	3.5626119	1503.94559	10.4908398	2000-2010 –> 2011-2019	7.4898860	3061.8794
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	277	3.5627927	1503.97371	10.4837990	2000-2010 –> 2011-2019	7.5038753	3063.8520
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + aou + silicate	277	3.5981804	1509.44921	10.9506822	2000-2010 –> 2011-2019	7.8125686	3106.4726
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	277	3.5630586	1504.01505	10.4751441	2000-2010 –> 2011-2019	7.5090473	3064.5825
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	277	3.5630586	1504.01505	10.4751441	2000-2010 –> 2011-2019	7.5090473	3064.5825
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	277	3.5630586	1504.01505	10.4751441	2000-2010 –> 2011-2019	7.5090473	3064.5825
Atlantic	(26.5,26.75]	4	cstar_tref ~ aou + silicate + phosphate + phosphate_star	432	4.2864097	2495.47124	11.7859659	1982-1999 –> 2000-2010	7.0044886	4331.2946
Atlantic	(26.5,26.75]	4	cstar_tref ~ oxygen + aou + silicate + phosphate	432	4.2864097	2495.47124	11.7859659	1982-1999 –> 2000-2010	7.0044886	4331.2946
Atlantic	(26.5,26.75]	5	cstar_tref ~ oxygen + aou + silicate + phosphate + phosphate_star	432	4.2864097	2495.47124	11.7859659	1982-1999 –> 2000-2010	7.0044886	4331.2946
Atlantic	(26.5,26.75]	4	cstar_tref ~ oxygen + aou + silicate + phosphate_star	432	4.2864097	2495.47124	11.7859659	1982-1999 –> 2000-2010	7.0044886	4331.2946
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	432	4.2782338	2495.82168	12.0215878	1982-1999 –> 2000-2010	6.9962003	4333.6139
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	432	4.2782338	2495.82168	12.0215878	1982-1999 –> 2000-2010	6.9962003	4333.6139
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	432	4.2782338	2495.82168	12.0215878	1982-1999 –> 2000-2010	6.9962003	4333.6139
Atlantic	(26.5,26.75]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	432	4.1715575	2474.00499	14.1348234	1982-1999 –> 2000-2010	6.8074187	4288.6690
Atlantic	(26.5,26.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	432	4.1715575	2474.00499	14.1348234	1982-1999 –> 2000-2010	6.8074187	4288.6690
Atlantic	(26.5,26.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	432	4.1715575	2474.00499	14.1348234	1982-1999 –> 2000-2010	6.8074187	4288.6690
Atlantic	(26.5,26.75]	4	cstar_tref ~ aou + silicate + phosphate + phosphate_star	267	4.0727826	1519.62350	13.9353007	2000-2010 –> 2011-2019	8.3591922	4015.0947
Atlantic	(26.5,26.75]	4	cstar_tref ~ oxygen + aou + silicate + phosphate	267	4.0727826	1519.62350	13.9353007	2000-2010 –> 2011-2019	8.3591922	4015.0947
Atlantic	(26.5,26.75]	5	cstar_tref ~ oxygen + aou + silicate + phosphate + phosphate_star	267	4.0727826	1519.62350	13.9353007	2000-2010 –> 2011-2019	8.3591922	4015.0947
Atlantic	(26.5,26.75]	4	cstar_tref ~ oxygen + aou + silicate + phosphate_star	267	4.0727826	1519.62350	13.9353007	2000-2010 –> 2011-2019	8.3591922	4015.0947
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	267	3.9257615	1501.99040	12.7011089	2000-2010 –> 2011-2019	8.2039953	3997.8121
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	267	3.9257615	1501.99040	12.7011089	2000-2010 –> 2011-2019	8.2039953	3997.8121
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	267	3.9257615	1501.99040	12.7011089	2000-2010 –> 2011-2019	8.2039953	3997.8121
Atlantic	(26.5,26.75]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	267	3.7649862	1479.66057	12.8051319	2000-2010 –> 2011-2019	7.9365436	3953.6656
Atlantic	(26.5,26.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	267	3.7649862	1479.66057	12.8051319	2000-2010 –> 2011-2019	7.9365436	3953.6656
Atlantic	(26.5,26.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	267	3.7649862	1479.66057	12.8051319	2000-2010 –> 2011-2019	7.9365436	3953.6656
Atlantic	(26.75,27]	4	cstar_tref ~ aou + silicate + phosphate + phosphate_star	780	2.3706716	3572.09444	12.3175612	1982-1999 –> 2000-2010	4.0912228	6632.3948
Atlantic	(26.75,27]	4	cstar_tref ~ oxygen + aou + silicate + phosphate	780	2.3706716	3572.09444	12.3175612	1982-1999 –> 2000-2010	4.0912228	6632.3948
Atlantic	(26.75,27]	5	cstar_tref ~ oxygen + aou + silicate + phosphate + phosphate_star	780	2.3706716	3572.09444	12.3175612	1982-1999 –> 2000-2010	4.0912228	6632.3948
Atlantic	(26.75,27]	4	cstar_tref ~ oxygen + aou + silicate + phosphate_star	780	2.3706716	3572.09444	12.3175612	1982-1999 –> 2000-2010	4.0912228	6632.3948
Atlantic	(26.75,27]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	780	2.3149833	3537.01197	11.5369885	1982-1999 –> 2000-2010	4.0348011	6598.6498
Atlantic	(26.75,27]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	780	2.3149833	3537.01197	11.5369885	1982-1999 –> 2000-2010	4.0348011	6598.6498
Atlantic	(26.75,27]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	780	2.3149833	3537.01197	11.5369885	1982-1999 –> 2000-2010	4.0348011	6598.6498
Atlantic	(26.75,27]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	780	2.3369656	3551.75533	11.8773130	1982-1999 –> 2000-2010	3.9220221	6486.5894
Atlantic	(26.75,27]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	780	2.3369656	3551.75533	11.8773130	1982-1999 –> 2000-2010	3.9220221	6486.5894
Atlantic	(26.75,27]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	780	2.3369656	3551.75533	11.8773130	1982-1999 –> 2000-2010	3.9220221	6486.5894
Atlantic	(26.75,27]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	453	2.8285090	2241.57156	11.2544240	2000-2010 –> 2011-2019	5.1434923	5778.5835
Atlantic	(26.75,27]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	453	2.8285090	2241.57156	11.2544240	2000-2010 –> 2011-2019	5.1434923	5778.5835
Atlantic	(26.75,27]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	453	2.8285090	2241.57156	11.2544240	2000-2010 –> 2011-2019	5.1434923	5778.5835
Atlantic	(26.75,27]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	453	2.8420420	2245.89597	10.9693234	2000-2010 –> 2011-2019	5.2238581	5827.3068
Atlantic	(26.75,27]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	453	2.8676549	2254.02440	10.9236197	2000-2010 –> 2011-2019	5.2996123	5867.9350
Atlantic	(26.75,27]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	453	2.8755839	2256.52601	10.8919500	2000-2010 –> 2011-2019	5.3273137	5883.0685
Atlantic	(26.75,27]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	453	2.8755839	2256.52601	10.8919500	2000-2010 –> 2011-2019	5.3273137	5883.0685
Atlantic	(26.75,27]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	453	2.9474945	2278.90397	13.8655751	2000-2010 –> 2011-2019	5.2844601	5830.6593
Atlantic	(26.75,27]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	453	2.9474945	2278.90397	13.8655751	2000-2010 –> 2011-2019	5.2844601	5830.6593
Atlantic	(26.75,27]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	453	2.9474945	2278.90397	13.8655751	2000-2010 –> 2011-2019	5.2844601	5830.6593
Atlantic	(27,27.25]	4	cstar_tref ~ sal + aou + silicate + phosphate	973	1.9233240	4046.04534	10.3027234	1982-1999 –> 2000-2010	3.4359444	7642.9670
Atlantic	(27,27.25]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	973	1.9202035	4044.88542	10.5016518	1982-1999 –> 2000-2010	3.2977529	7460.8480
Atlantic	(27,27.25]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	973	1.9202035	4044.88542	10.5016518	1982-1999 –> 2000-2010	3.2977529	7460.8480
Atlantic	(27,27.25]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	973	1.9202035	4044.88542	10.5016518	1982-1999 –> 2000-2010	3.2977529	7460.8480
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	973	1.9232342	4047.95447	10.3367133	1982-1999 –> 2000-2010	3.3326440	7508.6407
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	973	1.9554415	4080.27314	10.3569411	1982-1999 –> 2000-2010	3.4084839	7600.5950
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	973	1.9636666	4088.44136	10.3528847	1982-1999 –> 2000-2010	3.4295301	7625.9465
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	973	1.9636666	4088.44136	10.3528847	1982-1999 –> 2000-2010	3.4295301	7625.9465
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	973	1.9636666	4088.44136	10.3528847	1982-1999 –> 2000-2010	3.4295301	7625.9465
Atlantic	(27,27.25]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	973	2.4014607	4480.10054	13.2489503	1982-1999 –> 2000-2010	3.9874047	8171.6120
Atlantic	(27,27.25]	4	cstar_tref ~ sal + aou + silicate + phosphate	645	2.4758396	3011.91836	11.8623904	2000-2010 –> 2011-2019	4.3991636	7057.9637
Atlantic	(27,27.25]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	645	2.4686766	3010.18079	11.4746669	2000-2010 –> 2011-2019	4.3888801	7055.0662
Atlantic	(27,27.25]	4	cstar_tref ~ sal + aou + silicate + phosphate_star	645	2.4864169	3017.41777	10.8792317	2000-2010 –> 2011-2019	4.5525323	7202.8264
Atlantic	(27,27.25]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	645	2.4686766	3010.18079	11.4746669	2000-2010 –> 2011-2019	4.3888801	7055.0662
Atlantic	(27,27.25]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	645	2.4686766	3010.18079	11.4746669	2000-2010 –> 2011-2019	4.3888801	7055.0662
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	645	2.4633411	3007.38970	11.3683085	2000-2010 –> 2011-2019	4.3865753	7055.3442
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	645	2.4753411	3013.65859	11.3336158	2000-2010 –> 2011-2019	4.4307826	7093.9317
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	645	2.4767418	3014.38832	11.3279466	2000-2010 –> 2011-2019	4.4404084	7102.8297
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	645	2.4767418	3014.38832	11.3279466	2000-2010 –> 2011-2019	4.4404084	7102.8297
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	645	2.4767418	3014.38832	11.3279466	2000-2010 –> 2011-2019	4.4404084	7102.8297
Atlantic	(27.25,27.5]	4	cstar_tref ~ sal + aou + silicate + phosphate	1778	3.4618040	9473.55006	18.2896694	1982-1999 –> 2000-2010	4.9492109	14348.7252
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1778	3.4578124	9471.44751	17.8687120	1982-1999 –> 2000-2010	4.7741879	14021.4994
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1778	3.4578124	9471.44751	17.8687120	1982-1999 –> 2000-2010	4.7741879	14021.4994
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1778	3.4578124	9471.44751	17.8687120	1982-1999 –> 2000-2010	4.7741879	14021.4994
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1778	3.4264261	9439.02263	17.3757054	1982-1999 –> 2000-2010	4.6715162	13839.9787
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1778	3.3443670	9352.82402	18.9445315	1982-1999 –> 2000-2010	4.6335930	13847.0659
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	1778	3.5287573	9543.66856	17.5801085	1982-1999 –> 2000-2010	5.2524931	14815.7406
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	1778	3.3444455	9352.90743	18.9111532	1982-1999 –> 2000-2010	4.6495490	13879.9290
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	1778	3.3444455	9352.90743	18.9111532	1982-1999 –> 2000-2010	4.6495490	13879.9290
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1778	3.3444455	9352.90743	18.9111532	1982-1999 –> 2000-2010	4.6495490	13879.9290
Atlantic	(27.25,27.5]	4	cstar_tref ~ sal + aou + silicate + phosphate	1093	3.6258704	5929.57392	9.3970334	2000-2010 –> 2011-2019	7.0876743	15403.1240
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1093	3.4931988	5850.08749	9.9713630	2000-2010 –> 2011-2019	6.9510111	15321.5350
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1093	3.4931988	5850.08749	9.9713630	2000-2010 –> 2011-2019	6.9510111	15321.5350
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1093	3.4931988	5850.08749	9.9713630	2000-2010 –> 2011-2019	6.9510111	15321.5350
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1093	3.5765021	5901.60580	9.8741916	2000-2010 –> 2011-2019	7.0029282	15340.6284
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1093	3.4778623	5840.46898	9.6020522	2000-2010 –> 2011-2019	6.8222293	15193.2930
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	1093	3.4824405	5843.34473	9.6087676	2000-2010 –> 2011-2019	6.8268860	15196.2522
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	1093	3.4824405	5843.34473	9.6087676	2000-2010 –> 2011-2019	6.8268860	15196.2522
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1093	3.4824405	5843.34473	9.6087676	2000-2010 –> 2011-2019	6.8268860	15196.2522
Atlantic	(27.25,27.5]	4	cstar_tref ~ sal + temp + silicate + phosphate_star	1093	3.5384197	5876.20456	10.3237086	2000-2010 –> 2011-2019	6.8973118	15242.4394
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	3382	2.4664801	15718.18571	18.3162078	1982-1999 –> 2000-2010	4.1066384	26246.2254
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	3382	2.4664801	15718.18571	18.3162078	1982-1999 –> 2000-2010	4.1066384	26246.2254
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	3382	2.4664801	15718.18571	18.3162078	1982-1999 –> 2000-2010	4.1066384	26246.2254
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	3382	2.4584414	15696.10462	18.5445827	1982-1999 –> 2000-2010	3.9824632	25820.6651
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	3382	2.4201064	15589.80119	20.9365516	1982-1999 –> 2000-2010	3.8421521	25333.8674
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	3382	2.6182981	16122.21610	18.4947739	1982-1999 –> 2000-2010	4.2631430	26665.9321
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	3382	2.4210265	15592.37232	20.7844419	1982-1999 –> 2000-2010	3.8422643	25333.3164
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	3382	2.4210265	15592.37232	20.7844419	1982-1999 –> 2000-2010	3.8422643	25333.3164
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	3382	2.4210265	15592.37232	20.7844419	1982-1999 –> 2000-2010	3.8422643	25333.3164
Atlantic	(27.5,27.75]	4	cstar_tref ~ sal + temp + silicate + phosphate_star	3382	2.5438022	15924.97565	23.5564243	1982-1999 –> 2000-2010	3.9752396	25703.2072
Atlantic	(27.5,27.75]	4	cstar_tref ~ sal + aou + silicate + phosphate	1916	2.1467794	8376.90071	9.8632018	2000-2010 –> 2011-2019	4.6163147	24101.4600
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1916	2.1104367	8313.47361	9.8641429	2000-2010 –> 2011-2019	4.5769167	24031.6593
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1916	2.1104367	8313.47361	9.8641429	2000-2010 –> 2011-2019	4.5769167	24031.6593
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1916	2.1104367	8313.47361	9.8641429	2000-2010 –> 2011-2019	4.5769167	24031.6593
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1916	2.1242812	8338.52954	10.2208593	2000-2010 –> 2011-2019	4.5827225	24034.6342
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1916	2.0993600	8293.30830	10.0544725	2000-2010 –> 2011-2019	4.5194663	23883.1095
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	1916	2.0989582	8292.57491	10.0419290	2000-2010 –> 2011-2019	4.5199847	23884.9472
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	1916	2.0989582	8292.57491	10.0419290	2000-2010 –> 2011-2019	4.5199847	23884.9472
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1916	2.0989582	8292.57491	10.0419290	2000-2010 –> 2011-2019	4.5199847	23884.9472
Atlantic	(27.5,27.75]	4	cstar_tref ~ sal + temp + silicate + phosphate_star	1916	2.1761319	8428.94002	9.9926146	2000-2010 –> 2011-2019	4.7199341	24353.9157
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1397	1.1124513	4276.25944	5.8958942	1982-1999 –> 2000-2010	1.7722856	6226.3815
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1397	1.1124513	4276.25944	5.8958942	1982-1999 –> 2000-2010	1.7722856	6226.3815
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1397	1.1124513	4276.25944	5.8958942	1982-1999 –> 2000-2010	1.7722856	6226.3815
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	1397	1.0838308	4203.43640	6.2969151	1982-1999 –> 2000-2010	1.8058380	6327.3477
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1397	1.1326939	4326.64326	5.9509154	1982-1999 –> 2000-2010	1.8164670	6345.5454
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	1397	1.0838308	4203.43640	6.2969151	1982-1999 –> 2000-2010	1.8058380	6327.3477
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	1397	1.0838308	4203.43640	6.2969151	1982-1999 –> 2000-2010	1.8058380	6327.3477
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	1397	1.1830722	4448.22645	6.1874703	1982-1999 –> 2000-2010	1.9054935	6573.2445
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	1397	1.1830722	4448.22645	6.1874703	1982-1999 –> 2000-2010	1.9054935	6573.2445
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1397	1.1830722	4448.22645	6.1874703	1982-1999 –> 2000-2010	1.9054935	6573.2445
Atlantic	(27.75,27.85]	4	cstar_tref ~ sal + aou + silicate + phosphate	679	1.2482470	2240.04170	4.7629173	2000-2010 –> 2011-2019	2.4017958	6615.6594
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	679	1.1584359	2140.64061	5.2234979	2000-2010 –> 2011-2019	2.2708872	6416.9000
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	679	1.1584359	2140.64061	5.2234979	2000-2010 –> 2011-2019	2.2708872	6416.9000
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	679	1.1584359	2140.64061	5.2234979	2000-2010 –> 2011-2019	2.2708872	6416.9000
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	679	1.2920643	2288.89405	5.3305555	2000-2010 –> 2011-2019	2.3758952	6492.3305
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	679	1.1934737	2181.10546	5.2208713	2000-2010 –> 2011-2019	2.3261676	6507.7487
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	679	1.2920643	2288.89405	5.3305555	2000-2010 –> 2011-2019	2.3758952	6492.3305
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	679	1.2920643	2288.89405	5.3305555	2000-2010 –> 2011-2019	2.3758952	6492.3305
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	679	1.2421421	2235.38373	5.3625128	2000-2010 –> 2011-2019	2.4252143	6683.6102
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	679	1.2421421	2235.38373	5.3625128	2000-2010 –> 2011-2019	2.4252143	6683.6102
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	1286	1.2237783	4182.90744	8.9948698	1982-1999 –> 2000-2010	2.1637532	6995.1165
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	1286	1.2237783	4182.90744	8.9948698	1982-1999 –> 2000-2010	2.1637532	6995.1165
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	1286	1.2237783	4182.90744	8.9948698	1982-1999 –> 2000-2010	2.1637532	6995.1165
Atlantic	(27.85,27.95]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	1286	1.4992620	4703.10037	8.4817728	1982-1999 –> 2000-2010	2.5456473	7734.4462
Atlantic	(27.85,27.95]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1286	1.4935303	4695.24883	8.5353431	1982-1999 –> 2000-2010	2.5393132	7727.4070
Atlantic	(27.85,27.95]	4	cstar_tref ~ temp + oxygen + aou + phosphate	1286	1.4992620	4703.10037	8.4817728	1982-1999 –> 2000-2010	2.5456473	7734.4462
Atlantic	(27.85,27.95]	5	cstar_tref ~ temp + oxygen + aou + phosphate + phosphate_star	1286	1.4992620	4703.10037	8.4817728	1982-1999 –> 2000-2010	2.5456473	7734.4462
Atlantic	(27.85,27.95]	4	cstar_tref ~ temp + oxygen + aou + phosphate_star	1286	1.4992620	4703.10037	8.4817728	1982-1999 –> 2000-2010	2.5456473	7734.4462
Atlantic	(27.85,27.95]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	1286	1.4935303	4695.24883	8.5353431	1982-1999 –> 2000-2010	2.5393132	7727.4070
Atlantic	(27.85,27.95]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	1286	1.4935303	4695.24883	8.5353431	1982-1999 –> 2000-2010	2.5393132	7727.4070
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	637	1.3463948	2200.65420	8.9888677	2000-2010 –> 2011-2019	2.8498131	6912.8748
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	637	1.3463948	2200.65420	8.9888677	2000-2010 –> 2011-2019	2.8498131	6912.8748
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	637	1.3463948	2200.65420	8.9888677	2000-2010 –> 2011-2019	2.8498131	6912.8748
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	637	1.1507074	2000.56787	7.1023411	2000-2010 –> 2011-2019	2.3744857	6183.4753
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	637	1.3905923	2241.80341	8.9790188	2000-2010 –> 2011-2019	2.9642498	7071.4646
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	637	1.1507074	2000.56787	7.1023411	2000-2010 –> 2011-2019	2.3744857	6183.4753
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	637	1.1507074	2000.56787	7.1023411	2000-2010 –> 2011-2019	2.3744857	6183.4753
Atlantic	(27.85,27.95]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	637	1.5045308	2342.13259	8.4438768	2000-2010 –> 2011-2019	2.9980611	7037.3814
Atlantic	(27.85,27.95]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	637	1.5045308	2342.13259	8.4438768	2000-2010 –> 2011-2019	2.9980611	7037.3814
Atlantic	(27.85,27.95]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	637	1.5045308	2342.13259	8.4438768	2000-2010 –> 2011-2019	2.9980611	7037.3814
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1575	3.8701042	8746.49289	20.1232859	1982-1999 –> 2000-2010	7.4965627	15707.1555
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1575	3.8701042	8746.49289	20.1232859	1982-1999 –> 2000-2010	7.4965627	15707.1555
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1575	3.8701042	8746.49289	20.1232859	1982-1999 –> 2000-2010	7.4965627	15707.1555
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	1575	3.8929764	8765.05452	20.4674052	1982-1999 –> 2000-2010	7.6455495	15813.4363
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1575	3.8768233	8751.95700	20.0601522	1982-1999 –> 2000-2010	7.4969145	15708.1103
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1575	3.8831047	8757.05670	20.0203713	1982-1999 –> 2000-2010	7.5025713	15712.7672
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	1575	3.8929764	8765.05452	20.4674052	1982-1999 –> 2000-2010	7.6455495	15813.4363
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	1575	3.8831198	8757.06889	20.0271436	1982-1999 –> 2000-2010	7.5026232	15712.8055
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	1575	3.8831198	8757.06889	20.0271436	1982-1999 –> 2000-2010	7.5026232	15712.8055
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1575	3.8831198	8757.06889	20.0271436	1982-1999 –> 2000-2010	7.5026232	15712.8055
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	777	4.1715844	4438.60232	17.6910597	2000-2010 –> 2011-2019	8.0416886	13185.0952
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	777	4.1715844	4438.60232	17.6910597	2000-2010 –> 2011-2019	8.0416886	13185.0952
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	777	4.1715844	4438.60232	17.6910597	2000-2010 –> 2011-2019	8.0416886	13185.0952
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	777	4.2228128	4457.56973	16.9979082	2000-2010 –> 2011-2019	8.1157892	13222.6243
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	777	4.1965212	4447.86415	17.7042002	2000-2010 –> 2011-2019	8.0733445	13199.8211
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	777	4.2091020	4452.51592	17.7214344	2000-2010 –> 2011-2019	8.0922067	13209.5726
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	777	4.2228128	4457.56973	16.9979082	2000-2010 –> 2011-2019	8.1157892	13222.6243
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	777	4.2124016	4453.73366	17.6662298	2000-2010 –> 2011-2019	8.0955213	13210.8025
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	777	4.2124016	4453.73366	17.6662298	2000-2010 –> 2011-2019	8.0955213	13210.8025
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	777	4.2124016	4453.73366	17.6662298	2000-2010 –> 2011-2019	8.0955213	13210.8025
Atlantic	(28.05,28.1]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	895	0.8088375	2174.13850	5.7350315	1982-1999 –> 2000-2010	1.5236199	3765.2215
Atlantic	(28.05,28.1]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	895	0.8088375	2174.13850	5.7350315	1982-1999 –> 2000-2010	1.5236199	3765.2215
Atlantic	(28.05,28.1]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	895	0.8088375	2174.13850	5.7350315	1982-1999 –> 2000-2010	1.5236199	3765.2215
Atlantic	(28.05,28.1]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	895	0.8320381	2222.76004	5.4054206	1982-1999 –> 2000-2010	1.5579686	3834.3763
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	895	0.7500543	2039.07861	5.0602354	1982-1999 –> 2000-2010	1.4133809	3521.3829
Atlantic	(28.05,28.1]	4	cstar_tref ~ temp + oxygen + aou + phosphate	895	0.8320381	2222.76004	5.4054206	1982-1999 –> 2000-2010	1.5579686	3834.3763
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + phosphate + phosphate_star	895	0.8320381	2222.76004	5.4054206	1982-1999 –> 2000-2010	1.5579686	3834.3763
Atlantic	(28.05,28.1]	4	cstar_tref ~ temp + oxygen + aou + phosphate_star	895	0.8320381	2222.76004	5.4054206	1982-1999 –> 2000-2010	1.5579686	3834.3763
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	895	0.7500543	2039.07861	5.0602354	1982-1999 –> 2000-2010	1.4133809	3521.3829
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	895	0.7500543	2039.07861	5.0602354	1982-1999 –> 2000-2010	1.4133809	3521.3829
Atlantic	(28.05,28.1]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	490	1.2253550	1603.72576	6.4569819	2000-2010 –> 2011-2019	2.0341925	3777.8643
Atlantic	(28.05,28.1]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	490	1.2253550	1603.72576	6.4569819	2000-2010 –> 2011-2019	2.0341925	3777.8643
Atlantic	(28.05,28.1]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	490	1.2253550	1603.72576	6.4569819	2000-2010 –> 2011-2019	2.0341925	3777.8643
Atlantic	(28.05,28.1]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	490	1.2399787	1613.35212	6.6132976	2000-2010 –> 2011-2019	2.0720168	3836.1122
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	490	1.0343603	1437.66731	5.3808096	2000-2010 –> 2011-2019	1.7844146	3476.7459
Atlantic	(28.05,28.1]	4	cstar_tref ~ temp + oxygen + aou + phosphate	490	1.2399787	1613.35212	6.6132976	2000-2010 –> 2011-2019	2.0720168	3836.1122
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + phosphate + phosphate_star	490	1.2399787	1613.35212	6.6132976	2000-2010 –> 2011-2019	2.0720168	3836.1122
Atlantic	(28.05,28.1]	4	cstar_tref ~ temp + oxygen + aou + phosphate_star	490	1.2399787	1613.35212	6.6132976	2000-2010 –> 2011-2019	2.0720168	3836.1122
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	490	1.0343603	1437.66731	5.3808096	2000-2010 –> 2011-2019	1.7844146	3476.7459
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	490	1.0343603	1437.66731	5.3808096	2000-2010 –> 2011-2019	1.7844146	3476.7459
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1399	0.8434408	3507.78698	7.8049894	1982-1999 –> 2000-2010	1.4820841	5171.6869
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1399	0.8434408	3507.78698	7.8049894	1982-1999 –> 2000-2010	1.4820841	5171.6869
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1399	0.8434408	3507.78698	7.8049894	1982-1999 –> 2000-2010	1.4820841	5171.6869
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	1399	0.8031821	3370.94179	7.7722686	1982-1999 –> 2000-2010	1.4864544	5149.6722
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	1399	0.8031821	3370.94179	7.7722686	1982-1999 –> 2000-2010	1.4864544	5149.6722
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	1399	0.8031821	3370.94179	7.7722686	1982-1999 –> 2000-2010	1.4864544	5149.6722
Atlantic	(28.1,28.15]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	1399	0.8073850	3383.54475	7.8844218	1982-1999 –> 2000-2010	1.4942244	5169.1272
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1399	0.7100896	3026.25528	7.0219700	1982-1999 –> 2000-2010	1.2303096	4341.4949
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	1399	0.7100896	3026.25528	7.0219700	1982-1999 –> 2000-2010	1.2303096	4341.4949
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	1399	0.7100896	3026.25528	7.0219700	1982-1999 –> 2000-2010	1.2303096	4341.4949
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	621	0.9007216	1646.45926	3.0458261	2000-2010 –> 2011-2019	1.7039037	5017.4011
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	621	0.9007216	1646.45926	3.0458261	2000-2010 –> 2011-2019	1.7039037	5017.4011
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	621	0.9007216	1646.45926	3.0458261	2000-2010 –> 2011-2019	1.7039037	5017.4011
Atlantic	(28.1,28.15]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	621	0.9453157	1704.47612	3.3408409	2000-2010 –> 2011-2019	1.7527006	5088.0209
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	621	0.6260613	1194.68440	2.8562581	2000-2010 –> 2011-2019	1.3361509	4220.9397
Atlantic	(28.1,28.15]	4	cstar_tref ~ temp + oxygen + aou + phosphate	621	0.9453157	1704.47612	3.3408409	2000-2010 –> 2011-2019	1.7527006	5088.0209
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + oxygen + aou + phosphate + phosphate_star	621	0.9453157	1704.47612	3.3408409	2000-2010 –> 2011-2019	1.7527006	5088.0209
Atlantic	(28.1,28.15]	4	cstar_tref ~ temp + oxygen + aou + phosphate_star	621	0.9453157	1704.47612	3.3408409	2000-2010 –> 2011-2019	1.7527006	5088.0209
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	621	0.6260613	1194.68440	2.8562581	2000-2010 –> 2011-2019	1.3361509	4220.9397
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	621	0.6260613	1194.68440	2.8562581	2000-2010 –> 2011-2019	1.3361509	4220.9397
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2536	0.8423490	6340.69977	9.9812550	1982-1999 –> 2000-2010	1.6940127	10144.9251
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	2536	0.8423490	6340.69977	9.9812550	1982-1999 –> 2000-2010	1.6940127	10144.9251
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	2536	0.8423490	6340.69977	9.9812550	1982-1999 –> 2000-2010	1.6940127	10144.9251
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2536	0.7216666	5556.41046	20.2286210	1982-1999 –> 2000-2010	1.5478085	9268.9952
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	2536	0.8665054	6484.10506	9.8044278	1982-1999 –> 2000-2010	1.7311622	10333.9349
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	2536	0.7216666	5556.41046	20.2286210	1982-1999 –> 2000-2010	1.5478085	9268.9952
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	2536	0.7216666	5556.41046	20.2286210	1982-1999 –> 2000-2010	1.5478085	9268.9952
Atlantic	(28.15,28.2]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	2536	0.7833288	5972.26004	10.2652894	1982-1999 –> 2000-2010	1.6225799	9732.2639
Atlantic	(28.15,28.2]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	2536	0.7833288	5972.26004	10.2652894	1982-1999 –> 2000-2010	1.6225799	9732.2639
Atlantic	(28.15,28.2]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	2536	0.7833288	5972.26004	10.2652894	1982-1999 –> 2000-2010	1.6225799	9732.2639
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1369	0.8708265	3520.35402	4.1276825	2000-2010 –> 2011-2019	1.7131755	9861.0538
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1369	0.8708265	3520.35402	4.1276825	2000-2010 –> 2011-2019	1.7131755	9861.0538
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1369	0.8708265	3520.35402	4.1276825	2000-2010 –> 2011-2019	1.7131755	9861.0538
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	1369	0.5743379	2380.73050	2.1828004	2000-2010 –> 2011-2019	1.2960045	7937.1410
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1369	0.8877361	3573.01053	4.0437786	2000-2010 –> 2011-2019	1.7542415	10057.1156
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	1369	0.5743379	2380.73050	2.1828004	2000-2010 –> 2011-2019	1.2960045	7937.1410
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	1369	0.5743379	2380.73050	2.1828004	2000-2010 –> 2011-2019	1.2960045	7937.1410
Atlantic	(28.15,28.2]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1369	0.6460458	2702.86275	2.4391714	2000-2010 –> 2011-2019	1.4293746	8675.1228
Atlantic	(28.15,28.2]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	1369	0.6460458	2702.86275	2.4391714	2000-2010 –> 2011-2019	1.4293746	8675.1228
Atlantic	(28.15,28.2]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	1369	0.6460458	2702.86275	2.4391714	2000-2010 –> 2011-2019	1.4293746	8675.1228
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	4896	0.8314257	12100.51237	5.4365382	1982-1999 –> 2000-2010	1.6886230	18365.4075
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	4896	0.8314257	12100.51237	5.4365382	1982-1999 –> 2000-2010	1.6886230	18365.4075
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	4896	0.8314257	12100.51237	5.4365382	1982-1999 –> 2000-2010	1.6886230	18365.4075
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	4896	0.8324248	12112.27198	5.4172634	1982-1999 –> 2000-2010	1.6896385	18377.2622
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	4896	0.8356635	12150.29538	5.3797997	1982-1999 –> 2000-2010	1.6931540	18416.8811
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	4896	0.8356635	12150.29538	5.3797997	1982-1999 –> 2000-2010	1.6931540	18416.8811
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	4896	0.8356635	12150.29538	5.3797997	1982-1999 –> 2000-2010	1.6931540	18416.8811
Atlantic	(28.2, Inf]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	4896	0.8092275	11835.52225	6.5515722	1982-1999 –> 2000-2010	1.6894271	18231.2848
Atlantic	(28.2, Inf]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	4896	0.8092275	11835.52225	6.5515722	1982-1999 –> 2000-2010	1.6894271	18231.2848
Atlantic	(28.2, Inf]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	4896	0.8092275	11835.52225	6.5515722	1982-1999 –> 2000-2010	1.6894271	18231.2848
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1781	0.6973573	3784.31007	3.4220217	2000-2010 –> 2011-2019	1.5287830	15884.8224
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1781	0.6973573	3784.31007	3.4220217	2000-2010 –> 2011-2019	1.5287830	15884.8224
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1781	0.6973573	3784.31007	3.4220217	2000-2010 –> 2011-2019	1.5287830	15884.8224
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1781	0.6965792	3780.33326	3.4137373	2000-2010 –> 2011-2019	1.5290040	15892.6052
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	1781	0.6959392	3777.05900	3.4205564	2000-2010 –> 2011-2019	1.5316027	15927.3544
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	1781	0.6959392	3777.05900	3.4205564	2000-2010 –> 2011-2019	1.5316027	15927.3544
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1781	0.6959392	3777.05900	3.4205564	2000-2010 –> 2011-2019	1.5316027	15927.3544
Atlantic	(28.2, Inf]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1781	0.6992300	3793.86275	2.3254934	2000-2010 –> 2011-2019	1.5084575	15629.3850
Atlantic	(28.2, Inf]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	1781	0.6992300	3793.86275	2.3254934	2000-2010 –> 2011-2019	1.5084575	15629.3850
Atlantic	(28.2, Inf]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	1781	0.6992300	3793.86275	2.3254934	2000-2010 –> 2011-2019	1.5084575	15629.3850
Indo-Pacific	(-Inf,26]	4	cstar_tref ~ sal + aou + phosphate + phosphate_star	2050	5.9412796	13135.53858	24.3257524	1982-1999 –> 2000-2010	11.2710377	20265.8851
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2050	5.6930889	12962.58515	22.5350844	1982-1999 –> 2000-2010	10.8980233	20040.3772
Indo-Pacific	(-Inf,26]	4	cstar_tref ~ sal + oxygen + aou + phosphate	2050	5.9412796	13135.53858	24.3257524	1982-1999 –> 2000-2010	11.2710377	20265.8851
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + oxygen + aou + phosphate + phosphate_star	2050	5.9412796	13135.53858	24.3257524	1982-1999 –> 2000-2010	11.2710377	20265.8851
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	2050	5.6930889	12962.58515	22.5350844	1982-1999 –> 2000-2010	10.8980233	20040.3772
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	2050	5.6930889	12962.58515	22.5350844	1982-1999 –> 2000-2010	10.8980233	20040.3772
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2050	5.9108189	13116.46399	25.3089384	1982-1999 –> 2000-2010	10.9715903	20129.5970
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	2050	6.0292385	13197.79291	22.1680830	1982-1999 –> 2000-2010	11.0102769	20174.3690
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	2050	5.9108189	13116.46399	25.3089384	1982-1999 –> 2000-2010	10.9715903	20129.5970
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	2050	5.9108189	13116.46399	25.3089384	1982-1999 –> 2000-2010	10.9715903	20129.5970
Indo-Pacific	(-Inf,26]	4	cstar_tref ~ sal + aou + phosphate + phosphate_star	2715	5.5074880	16981.00605	23.3761365	2000-2010 –> 2011-2019	11.4487676	30116.5446
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2715	5.1397077	16607.72567	22.9684927	2000-2010 –> 2011-2019	10.8327966	29570.3108
Indo-Pacific	(-Inf,26]	4	cstar_tref ~ sal + oxygen + aou + phosphate	2715	5.5074880	16981.00605	23.3761365	2000-2010 –> 2011-2019	11.4487676	30116.5446
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + oxygen + aou + phosphate + phosphate_star	2715	5.5074880	16981.00605	23.3761365	2000-2010 –> 2011-2019	11.4487676	30116.5446
Indo-Pacific	(-Inf,26]	4	cstar_tref ~ sal + oxygen + aou + phosphate_star	2715	5.5074880	16981.00605	23.3761365	2000-2010 –> 2011-2019	11.4487676	30116.5446
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	2715	5.1397077	16607.72567	22.9684927	2000-2010 –> 2011-2019	10.8327966	29570.3108
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	2715	5.1397077	16607.72567	22.9684927	2000-2010 –> 2011-2019	10.8327966	29570.3108
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2715	5.3282374	16803.33754	24.5273202	2000-2010 –> 2011-2019	11.2390564	29919.8015
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	2715	5.3282374	16803.33754	24.5273202	2000-2010 –> 2011-2019	11.2390564	29919.8015
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	2715	5.3282374	16803.33754	24.5273202	2000-2010 –> 2011-2019	11.2390564	29919.8015
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2283	4.4485564	13307.99201	44.1202214	1982-1999 –> 2000-2010	8.2297515	22300.1563
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	2283	4.4485564	13307.99201	44.1202214	1982-1999 –> 2000-2010	8.2297515	22300.1563
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	2283	4.4485564	13307.99201	44.1202214	1982-1999 –> 2000-2010	8.2297515	22300.1563
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2283	5.0339336	13872.45037	41.8370566	1982-1999 –> 2000-2010	9.1919965	23174.9147
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	2283	4.5185718	13379.29612	48.6125461	1982-1999 –> 2000-2010	8.3739825	22434.9430
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	2283	4.6682063	13528.05161	48.6049889	1982-1999 –> 2000-2010	8.6077164	22654.1749
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	2283	5.0339336	13872.45037	41.8370566	1982-1999 –> 2000-2010	9.1919965	23174.9147
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	2283	4.6948097	13553.99864	48.5144649	1982-1999 –> 2000-2010	8.6472125	22690.7931
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	2283	4.6948097	13553.99864	48.5144649	1982-1999 –> 2000-2010	8.6472125	22690.7931
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	2283	4.6948097	13553.99864	48.5144649	1982-1999 –> 2000-2010	8.6472125	22690.7931
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1803	3.5591758	9708.61392	21.3851029	2000-2010 –> 2011-2019	8.0077322	23016.6059
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1803	3.5591758	9708.61392	21.3851029	2000-2010 –> 2011-2019	8.0077322	23016.6059
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1803	3.5591758	9708.61392	21.3851029	2000-2010 –> 2011-2019	8.0077322	23016.6059
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1803	3.6035978	9753.34179	21.4028549	2000-2010 –> 2011-2019	8.1221696	23132.6379
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1803	3.6244598	9774.15753	21.2670831	2000-2010 –> 2011-2019	8.2926662	23302.2091
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	1803	3.6353226	9784.94882	21.2582420	2000-2010 –> 2011-2019	8.3301323	23338.9475
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	1803	3.6353226	9784.94882	21.2582420	2000-2010 –> 2011-2019	8.3301323	23338.9475
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1803	3.6353226	9784.94882	21.2582420	2000-2010 –> 2011-2019	8.3301323	23338.9475
Indo-Pacific	(26,26.5]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1803	3.9251310	10061.53578	23.2882408	2000-2010 –> 2011-2019	8.8686765	23851.2551
Indo-Pacific	(26,26.5]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	1803	3.9251310	10061.53578	23.2882408	2000-2010 –> 2011-2019	8.8686765	23851.2551
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2021	4.5190123	11845.87167	31.8167001	1982-1999 –> 2000-2010	7.6936909	18830.6562
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	2021	4.5190123	11845.87167	31.8167001	1982-1999 –> 2000-2010	7.6936909	18830.6562
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	2021	4.5190123	11845.87167	31.8167001	1982-1999 –> 2000-2010	7.6936909	18830.6562
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2021	4.6895028	11995.55957	21.1750599	1982-1999 –> 2000-2010	8.1009348	19175.1189
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	2021	4.7452341	12043.31254	32.7360586	1982-1999 –> 2000-2010	8.1386583	19208.5394
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	2021	4.6895028	11995.55957	21.1750599	1982-1999 –> 2000-2010	8.1009348	19175.1189
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	2021	4.6895028	11995.55957	21.1750599	1982-1999 –> 2000-2010	8.1009348	19175.1189
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	2021	4.7322021	12032.19657	26.1627880	1982-1999 –> 2000-2010	8.1710013	19233.3935
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	2021	4.7322021	12032.19657	26.1627880	1982-1999 –> 2000-2010	8.1710013	19233.3935
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	2021	4.7322021	12032.19657	26.1627880	1982-1999 –> 2000-2010	8.1710013	19233.3935
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1665	3.0628528	8466.49008	17.2092816	2000-2010 –> 2011-2019	7.5818651	20312.3618
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1665	3.0628528	8466.49008	17.2092816	2000-2010 –> 2011-2019	7.5818651	20312.3618
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1665	3.0628528	8466.49008	17.2092816	2000-2010 –> 2011-2019	7.5818651	20312.3618
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	1665	3.2280370	8641.40642	14.0907509	2000-2010 –> 2011-2019	7.9175398	20636.9660
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1665	3.2269744	8640.31007	18.7487239	2000-2010 –> 2011-2019	7.9722085	20683.6226
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	1665	3.2280370	8641.40642	14.0907509	2000-2010 –> 2011-2019	7.9175398	20636.9660
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	1665	3.2280370	8641.40642	14.0907509	2000-2010 –> 2011-2019	7.9175398	20636.9660
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1665	2.9313951	8320.40861	18.8470683	2000-2010 –> 2011-2019	7.6635972	20352.6052
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	1665	2.9313951	8320.40861	18.8470683	2000-2010 –> 2011-2019	7.6635972	20352.6052
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	1665	2.9313951	8320.40861	18.8470683	2000-2010 –> 2011-2019	7.6635972	20352.6052
Indo-Pacific	(26.75,27]	4	cstar_tref ~ sal + aou + phosphate + phosphate_star	3039	4.3771932	17609.91435	21.3309545	1982-1999 –> 2000-2010	7.5809190	27753.4326
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	3039	4.3722797	17605.08785	22.4382098	1982-1999 –> 2000-2010	7.5313617	27695.5687
Indo-Pacific	(26.75,27]	4	cstar_tref ~ sal + oxygen + aou + phosphate	3039	4.3771932	17609.91435	21.3309545	1982-1999 –> 2000-2010	7.5809190	27753.4326
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + oxygen + aou + phosphate + phosphate_star	3039	4.3771932	17609.91435	21.3309545	1982-1999 –> 2000-2010	7.5809190	27753.4326
Indo-Pacific	(26.75,27]	4	cstar_tref ~ sal + oxygen + aou + phosphate_star	3039	4.3771932	17609.91435	21.3309545	1982-1999 –> 2000-2010	7.5809190	27753.4326
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	3039	4.3722797	17605.08785	22.4382098	1982-1999 –> 2000-2010	7.5313617	27695.5687
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	3039	4.3722797	17605.08785	22.4382098	1982-1999 –> 2000-2010	7.5313617	27695.5687
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	3039	4.3403071	17560.47879	18.5913861	1982-1999 –> 2000-2010	7.4801719	27627.0287
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	3039	4.3403071	17560.47879	18.5913861	1982-1999 –> 2000-2010	7.4801719	27627.0287
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	3039	4.3403071	17560.47879	18.5913861	1982-1999 –> 2000-2010	7.4801719	27627.0287
Indo-Pacific	(26.75,27]	4	cstar_tref ~ oxygen + aou + silicate + phosphate_star	2102	3.6869201	11462.56090	13.3263233	2000-2010 –> 2011-2019	8.2852018	29371.9687
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2102	3.6401250	11410.86134	12.4834697	2000-2010 –> 2011-2019	8.0124047	29015.9492
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	2102	3.6401250	11410.86134	12.4834697	2000-2010 –> 2011-2019	8.0124047	29015.9492
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	2102	3.6401250	11410.86134	12.4834697	2000-2010 –> 2011-2019	8.0124047	29015.9492
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2102	3.7488465	11534.58576	13.8872680	2000-2010 –> 2011-2019	8.0891536	29095.0645
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	2102	3.7488465	11534.58576	13.8872680	2000-2010 –> 2011-2019	8.0891536	29095.0645
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	2102	3.7488465	11534.58576	13.8872680	2000-2010 –> 2011-2019	8.0891536	29095.0645
Indo-Pacific	(26.75,27]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	2102	3.6267593	11395.39682	12.9256801	2000-2010 –> 2011-2019	8.0924123	29128.9184
Indo-Pacific	(26.75,27]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	2102	3.6267593	11395.39682	12.9256801	2000-2010 –> 2011-2019	8.0924123	29128.9184
Indo-Pacific	(26.75,27]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	2102	3.6267593	11395.39682	12.9256801	2000-2010 –> 2011-2019	8.0924123	29128.9184
Indo-Pacific	(27,27.25]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	3979	3.9502666	22238.47846	38.9195516	1982-1999 –> 2000-2010	6.8731853	33703.5099
Indo-Pacific	(27,27.25]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	3979	3.9502666	22238.47846	38.9195516	1982-1999 –> 2000-2010	6.8731853	33703.5099
Indo-Pacific	(27,27.25]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	3979	3.9502666	22238.47846	38.9195516	1982-1999 –> 2000-2010	6.8731853	33703.5099
Indo-Pacific	(27,27.25]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	3979	4.0065357	22349.03545	37.2253026	1982-1999 –> 2000-2010	6.9833024	33895.9671
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	3979	3.5759607	21446.26606	21.6002826	1982-1999 –> 2000-2010	6.5275787	32956.2042
Indo-Pacific	(27,27.25]	4	cstar_tref ~ temp + oxygen + aou + phosphate	3979	4.0065357	22349.03545	37.2253026	1982-1999 –> 2000-2010	6.9833024	33895.9671
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + oxygen + aou + phosphate + phosphate_star	3979	4.0065357	22349.03545	37.2253026	1982-1999 –> 2000-2010	6.9833024	33895.9671
Indo-Pacific	(27,27.25]	4	cstar_tref ~ temp + oxygen + aou + phosphate_star	3979	4.0065357	22349.03545	37.2253026	1982-1999 –> 2000-2010	6.9833024	33895.9671
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	3979	3.5759607	21446.26606	21.6002826	1982-1999 –> 2000-2010	6.5275787	32956.2042
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	3979	3.5759607	21446.26606	21.6002826	1982-1999 –> 2000-2010	6.5275787	32956.2042
Indo-Pacific	(27,27.25]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	3025	4.1547530	17215.30870	19.9509812	2000-2010 –> 2011-2019	8.1050196	39453.7872
Indo-Pacific	(27,27.25]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	3025	4.1547530	17215.30870	19.9509812	2000-2010 –> 2011-2019	8.1050196	39453.7872
Indo-Pacific	(27,27.25]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	3025	4.1547530	17215.30870	19.9509812	2000-2010 –> 2011-2019	8.1050196	39453.7872
Indo-Pacific	(27,27.25]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	3025	4.2648053	17371.47709	18.4788061	2000-2010 –> 2011-2019	8.2713409	39720.5125
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	3025	3.8673130	16781.56590	16.2728164	2000-2010 –> 2011-2019	7.4432738	38227.8320
Indo-Pacific	(27,27.25]	4	cstar_tref ~ temp + oxygen + aou + phosphate	3025	4.2648053	17371.47709	18.4788061	2000-2010 –> 2011-2019	8.2713409	39720.5125
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + oxygen + aou + phosphate + phosphate_star	3025	4.2648053	17371.47709	18.4788061	2000-2010 –> 2011-2019	8.2713409	39720.5125
Indo-Pacific	(27,27.25]	4	cstar_tref ~ temp + oxygen + aou + phosphate_star	3025	4.2648053	17371.47709	18.4788061	2000-2010 –> 2011-2019	8.2713409	39720.5125
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	3025	3.8673130	16781.56590	16.2728164	2000-2010 –> 2011-2019	7.4432738	38227.8320
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	3025	3.8673130	16781.56590	16.2728164	2000-2010 –> 2011-2019	7.4432738	38227.8320
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	5002	3.6304492	27107.78232	11.3237050	1982-1999 –> 2000-2010	6.2427401	41596.6290
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	5002	3.6304492	27107.78232	11.3237050	1982-1999 –> 2000-2010	6.2427401	41596.6290
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	5002	3.6304492	27107.78232	11.3237050	1982-1999 –> 2000-2010	6.2427401	41596.6290
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	5002	3.5766279	26958.36317	10.8565914	1982-1999 –> 2000-2010	6.1889070	41447.1822
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + aou + silicate	5002	3.1524271	25695.38036	16.3912125	1982-1999 –> 2000-2010	6.1565545	41034.5242
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	5002	3.5972791	27015.95942	11.0622260	1982-1999 –> 2000-2010	6.2501335	41598.5521
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	5002	3.0933933	25506.26447	23.9809308	1982-1999 –> 2000-2010	5.3397267	39076.8635
Indo-Pacific	(27.25,27.5]	4	cstar_tref ~ temp + oxygen + aou + silicate	5002	3.2183035	25900.28015	20.8693095	1982-1999 –> 2000-2010	6.2250445	41242.7148
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	5002	3.0933933	25506.26447	23.9809308	1982-1999 –> 2000-2010	5.3397267	39076.8635
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	5002	3.0933933	25506.26447	23.9809308	1982-1999 –> 2000-2010	5.3397267	39076.8635
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	4095	3.3211328	21465.61209	14.7808467	2000-2010 –> 2011-2019	6.8977607	48423.9753
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	4095	3.3429015	21519.11914	15.1346525	2000-2010 –> 2011-2019	6.9631698	48598.8078
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + aou + silicate	4095	2.8828039	20306.37872	11.6411635	2000-2010 –> 2011-2019	6.0352310	46001.7591
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	4095	3.3235812	21471.64772	14.9133477	2000-2010 –> 2011-2019	6.9208603	48487.6071
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	4095	3.3235812	21471.64772	14.9133477	2000-2010 –> 2011-2019	6.9208603	48487.6071
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	4095	3.3235812	21471.64772	14.9133477	2000-2010 –> 2011-2019	6.9208603	48487.6071
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	4095	3.3096339	21437.20634	14.6982873	2000-2010 –> 2011-2019	6.4030272	46943.4708
Indo-Pacific	(27.25,27.5]	4	cstar_tref ~ temp + oxygen + aou + silicate	4095	3.3326268	21491.90777	15.4739517	2000-2010 –> 2011-2019	6.5509303	47392.1879
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	4095	3.3096339	21437.20634	14.6982873	2000-2010 –> 2011-2019	6.4030272	46943.4708
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	4095	3.3096339	21437.20634	14.6982873	2000-2010 –> 2011-2019	6.4030272	46943.4708
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	5970	3.1096083	30502.01736	10.4062902	1982-1999 –> 2000-2010	5.4111798	45734.1148
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	5970	3.1198812	30541.39758	10.0999699	1982-1999 –> 2000-2010	5.4484959	45852.4143
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + aou + silicate	5970	2.9188701	29746.21332	18.7090746	1982-1999 –> 2000-2010	5.2244038	44989.9313
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	5970	3.1112033	30508.14028	10.1635385	1982-1999 –> 2000-2010	5.4351690	45805.6555
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	5970	3.1112033	30508.14028	10.1635385	1982-1999 –> 2000-2010	5.4351690	45805.6555
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	5970	3.1112033	30508.14028	10.1635385	1982-1999 –> 2000-2010	5.4351690	45805.6555
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	5970	3.0564183	30296.01629	22.5240891	1982-1999 –> 2000-2010	5.1542800	44902.0111
Indo-Pacific	(27.5,27.75]	4	cstar_tref ~ temp + oxygen + aou + silicate	5970	3.0945669	30442.12259	25.4156429	1982-1999 –> 2000-2010	5.4357560	45787.5233
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	5970	3.0564183	30296.01629	22.5240891	1982-1999 –> 2000-2010	5.1542800	44902.0111
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	5970	3.0564183	30296.01629	22.5240891	1982-1999 –> 2000-2010	5.1542800	44902.0111
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	4550	2.6059992	21642.46777	12.7070322	2000-2010 –> 2011-2019	5.7533577	52288.5626
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	4550	2.6059992	21642.46777	12.7070322	2000-2010 –> 2011-2019	5.7533577	52288.5626
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	4550	2.6059992	21642.46777	12.7070322	2000-2010 –> 2011-2019	5.7533577	52288.5626
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	4550	2.5973042	21612.05448	12.2057260	2000-2010 –> 2011-2019	5.7069124	52114.0718
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	4550	2.5978634	21614.01348	12.2609763	2000-2010 –> 2011-2019	5.7177446	52155.4111
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + aou + silicate	4550	2.6187797	21686.98756	12.1127746	2000-2010 –> 2011-2019	5.5376498	51433.2009
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	4550	2.5963985	21608.88080	12.1423519	2000-2010 –> 2011-2019	5.7076018	52117.0211
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	4550	2.5963985	21608.88080	12.1423519	2000-2010 –> 2011-2019	5.7076018	52117.0211
Indo-Pacific	(27.5,27.75]	4	cstar_tref ~ sal + temp + silicate + phosphate	4550	2.8405890	22424.84454	8.4616010	2000-2010 –> 2011-2019	5.9584560	52956.5310
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	4550	2.5963985	21608.88080	12.1423519	2000-2010 –> 2011-2019	5.7076018	52117.0211
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2343	3.0753182	11927.43757	17.8219200	1982-1999 –> 2000-2010	4.6917565	17137.5494
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	2343	3.0753182	11927.43757	17.8219200	1982-1999 –> 2000-2010	4.6917565	17137.5494
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	2343	3.0753182	11927.43757	17.8219200	1982-1999 –> 2000-2010	4.6917565	17137.5494
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + aou + silicate	2343	2.9855813	11788.66684	19.3320940	1982-1999 –> 2000-2010	4.6699548	17111.4159
Indo-Pacific	(27.75,27.85]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	2343	3.1310570	12009.60885	17.8868431	1982-1999 –> 2000-2010	4.7487944	17219.9186
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	2343	2.9770264	11775.22031	14.7864653	1982-1999 –> 2000-2010	4.5130561	16845.7299
Indo-Pacific	(27.75,27.85]	4	cstar_tref ~ temp + oxygen + aou + phosphate	2343	3.1310570	12009.60885	17.8868431	1982-1999 –> 2000-2010	4.7487944	17219.9186
Indo-Pacific	(27.75,27.85]	4	cstar_tref ~ temp + oxygen + aou + silicate	2343	3.0233756	11845.61433	18.7438046	1982-1999 –> 2000-2010	4.7144255	17177.1866
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	2343	2.9770264	11775.22031	14.7864653	1982-1999 –> 2000-2010	4.5130561	16845.7299
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	2343	2.9770264	11775.22031	14.7864653	1982-1999 –> 2000-2010	4.5130561	16845.7299
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1782	2.3553713	8124.36996	9.7028538	2000-2010 –> 2011-2019	5.2857023	19825.5067
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1782	2.3553713	8124.36996	9.7028538	2000-2010 –> 2011-2019	5.2857023	19825.5067
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1782	2.3553713	8124.36996	9.7028538	2000-2010 –> 2011-2019	5.2857023	19825.5067
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1782	2.3721654	8149.69155	9.1446736	2000-2010 –> 2011-2019	5.2995538	19846.1202
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1782	2.3771856	8157.22612	8.9795391	2000-2010 –> 2011-2019	5.3033999	19851.7750
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + aou + silicate	1782	2.7388443	8661.95535	12.8709329	2000-2010 –> 2011-2019	5.7244256	20450.6222
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	1782	2.3805408	8162.25276	8.8868318	2000-2010 –> 2011-2019	5.3061879	19855.8933
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	1782	2.3805408	8162.25276	8.8868318	2000-2010 –> 2011-2019	5.3061879	19855.8933
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1782	2.3805408	8162.25276	8.8868318	2000-2010 –> 2011-2019	5.3061879	19855.8933
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	1782	2.7841169	8720.38609	9.1110338	2000-2010 –> 2011-2019	5.7611433	20495.6064
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2770	2.8175305	13613.58821	31.5875052	1982-1999 –> 2000-2010	4.2881316	18962.0147
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	2770	2.8175305	13613.58821	31.5875052	1982-1999 –> 2000-2010	4.2881316	18962.0147
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	2770	2.8175305	13613.58821	31.5875052	1982-1999 –> 2000-2010	4.2881316	18962.0147
Indo-Pacific	(27.85,27.95]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	2770	2.8491816	13673.47564	31.8970705	1982-1999 –> 2000-2010	4.3218890	19024.1328
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	2770	2.7403816	13459.77780	29.3662247	1982-1999 –> 2000-2010	4.1745617	18734.0741
Indo-Pacific	(27.85,27.95]	4	cstar_tref ~ temp + oxygen + aou + phosphate	2770	2.8491816	13673.47564	31.8970705	1982-1999 –> 2000-2010	4.3218890	19024.1328
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ temp + oxygen + aou + phosphate + phosphate_star	2770	2.8491816	13673.47564	31.8970705	1982-1999 –> 2000-2010	4.3218890	19024.1328
Indo-Pacific	(27.85,27.95]	4	cstar_tref ~ temp + oxygen + aou + phosphate_star	2770	2.8491816	13673.47564	31.8970705	1982-1999 –> 2000-2010	4.3218890	19024.1328
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	2770	2.7403816	13459.77780	29.3662247	1982-1999 –> 2000-2010	4.1745617	18734.0741
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	2770	2.7403816	13459.77780	29.3662247	1982-1999 –> 2000-2010	4.1745617	18734.0741
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2215	3.2760268	11556.67453	12.9319322	2000-2010 –> 2011-2019	5.9927242	24968.3639
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	2215	3.2760268	11556.67453	12.9319322	2000-2010 –> 2011-2019	5.9927242	24968.3639
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	2215	3.2760268	11556.67453	12.9319322	2000-2010 –> 2011-2019	5.9927242	24968.3639
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	2215	3.2996846	11588.55085	12.7253385	2000-2010 –> 2011-2019	6.0166615	25000.8100
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	2215	3.3118980	11604.91774	12.6978259	2000-2010 –> 2011-2019	6.0282931	25015.9904
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	2215	3.3172216	11612.03283	12.6560850	2000-2010 –> 2011-2019	6.0333948	25022.6531
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	2215	3.3172216	11612.03283	12.6560850	2000-2010 –> 2011-2019	6.0333948	25022.6531
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	2215	3.3172216	11612.03283	12.6560850	2000-2010 –> 2011-2019	6.0333948	25022.6531
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	2215	3.3987146	11719.54775	11.7616526	2000-2010 –> 2011-2019	6.1390962	25179.3255
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	2215	3.3987146	11719.54775	11.7616526	2000-2010 –> 2011-2019	6.1390962	25179.3255
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	3088	1.7138892	12104.77821	10.5968274	1982-1999 –> 2000-2010	3.3668082	19028.4251
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	3088	1.7138892	12104.77821	10.5968274	1982-1999 –> 2000-2010	3.3668082	19028.4251
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	3088	1.7138892	12104.77821	10.5968274	1982-1999 –> 2000-2010	3.3668082	19028.4251
Indo-Pacific	(27.95,28.05]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	3088	1.7296121	12159.17738	11.7000007	1982-1999 –> 2000-2010	3.3851629	19086.5454
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	3088	1.6471169	11859.35156	9.8766937	1982-1999 –> 2000-2010	3.1746633	18499.3460
Indo-Pacific	(27.95,28.05]	4	cstar_tref ~ temp + oxygen + aou + phosphate	3088	1.7296121	12159.17738	11.7000007	1982-1999 –> 2000-2010	3.3851629	19086.5454
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ temp + oxygen + aou + phosphate + phosphate_star	3088	1.7296121	12159.17738	11.7000007	1982-1999 –> 2000-2010	3.3851629	19086.5454
Indo-Pacific	(27.95,28.05]	4	cstar_tref ~ temp + oxygen + aou + phosphate_star	3088	1.7296121	12159.17738	11.7000007	1982-1999 –> 2000-2010	3.3851629	19086.5454
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	3088	1.6471169	11859.35156	9.8766937	1982-1999 –> 2000-2010	3.1746633	18499.3460
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	3088	1.6471169	11859.35156	9.8766937	1982-1999 –> 2000-2010	3.1746633	18499.3460
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2983	1.4176586	10561.56080	9.0040887	2000-2010 –> 2011-2019	3.1929797	22883.8338
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	2983	1.4176586	10561.56080	9.0040887	2000-2010 –> 2011-2019	3.1929797	22883.8338
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	2983	1.4176586	10561.56080	9.0040887	2000-2010 –> 2011-2019	3.1929797	22883.8338
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2983	1.5021010	10906.74260	9.5241591	2000-2010 –> 2011-2019	3.2159902	23011.5208
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	2983	1.4174570	10560.71241	8.9760277	2000-2010 –> 2011-2019	3.2142133	22957.1077
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	2983	1.5021010	10906.74260	9.5241591	2000-2010 –> 2011-2019	3.2159902	23011.5208
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	2983	1.5021010	10906.74260	9.5241591	2000-2010 –> 2011-2019	3.2159902	23011.5208
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	2983	1.5713386	11175.58889	9.7654130	2000-2010 –> 2011-2019	3.2184555	23034.9404
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	2983	1.5713386	11175.58889	9.7654130	2000-2010 –> 2011-2019	3.2184555	23034.9404
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	2983	1.5713386	11175.58889	9.7654130	2000-2010 –> 2011-2019	3.2184555	23034.9404
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1757	1.6129961	6680.17025	7.2479230	1982-1999 –> 2000-2010	3.1573990	10835.0588
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1757	1.6129961	6680.17025	7.2479230	1982-1999 –> 2000-2010	3.1573990	10835.0588
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1757	1.6129961	6680.17025	7.2479230	1982-1999 –> 2000-2010	3.1573990	10835.0588
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1757	1.6278501	6712.38221	7.3300843	1982-1999 –> 2000-2010	3.1825746	10882.1514
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	1757	1.6317590	6720.81019	7.3626201	1982-1999 –> 2000-2010	3.1900625	10895.7162
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	1757	1.6317590	6720.81019	7.3626201	1982-1999 –> 2000-2010	3.1900625	10895.7162
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1757	1.6317590	6720.81019	7.3626201	1982-1999 –> 2000-2010	3.1900625	10895.7162
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1757	1.5042995	6435.01218	7.7793210	1982-1999 –> 2000-2010	2.9453272	10435.1274
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	1757	1.5042995	6435.01218	7.7793210	1982-1999 –> 2000-2010	2.9453272	10435.1274
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	1757	1.5042995	6435.01218	7.7793210	1982-1999 –> 2000-2010	2.9453272	10435.1274
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1692	1.2298740	5515.87718	6.2087546	2000-2010 –> 2011-2019	2.8428701	12196.0474
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1692	1.2298740	5515.87718	6.2087546	2000-2010 –> 2011-2019	2.8428701	12196.0474
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1692	1.2298740	5515.87718	6.2087546	2000-2010 –> 2011-2019	2.8428701	12196.0474
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1692	1.2294618	5514.74292	6.2020103	2000-2010 –> 2011-2019	2.8573119	12227.1251
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1692	1.2293592	5514.46051	6.1998240	2000-2010 –> 2011-2019	2.8610950	12235.2208
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	1692	1.2293967	5514.56371	6.2015877	2000-2010 –> 2011-2019	2.8611557	12235.3739
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1692	1.2293967	5514.56371	6.2015877	2000-2010 –> 2011-2019	2.8611557	12235.3739
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1692	1.3298241	5780.28602	6.3333437	2000-2010 –> 2011-2019	2.8341236	12215.2982
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	1692	1.3298241	5780.28602	6.3333437	2000-2010 –> 2011-2019	2.8341236	12215.2982
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	1692	1.3298241	5780.28602	6.3333437	2000-2010 –> 2011-2019	2.8341236	12215.2982
Indo-Pacific	(28.1, Inf]	4	cstar_tref ~ sal + aou + silicate + phosphate	19958	1.5925342	75224.32566	9.6327805	1982-1999 –> 2000-2010	3.0083850	120212.1892
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	19958	1.5448820	74013.71690	8.9542419	1982-1999 –> 2000-2010	2.8524814	116978.7098
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	19958	1.5448820	74013.71690	8.9542419	1982-1999 –> 2000-2010	2.8524814	116978.7098
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	19958	1.5448820	74013.71690	8.9542419	1982-1999 –> 2000-2010	2.8524814	116978.7098
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	19958	1.5451511	74020.66934	8.9573978	1982-1999 –> 2000-2010	2.8519343	116969.7651
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	19958	1.5454426	74028.19785	8.9601785	1982-1999 –> 2000-2010	2.8512198	116957.6887
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	19958	1.5454214	74027.64935	8.9605169	1982-1999 –> 2000-2010	2.8512395	116957.9386
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	19958	1.5454214	74027.64935	8.9605169	1982-1999 –> 2000-2010	2.8512395	116957.9386
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	19958	1.5454214	74027.64935	8.9605169	1982-1999 –> 2000-2010	2.8512395	116957.9386
Indo-Pacific	(28.1, Inf]	4	cstar_tref ~ sal + temp + silicate + phosphate_star	19958	1.5546763	74263.97734	8.5745731	1982-1999 –> 2000-2010	2.8735453	117445.4427
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	18875	1.4120953	66605.49622	8.7905474	2000-2010 –> 2011-2019	2.9569773	140619.2131
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	18875	1.4120953	66605.49622	8.7905474	2000-2010 –> 2011-2019	2.9569773	140619.2131
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	18875	1.4120953	66605.49622	8.7905474	2000-2010 –> 2011-2019	2.9569773	140619.2131
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	18875	1.4122891	66610.67876	8.7906147	2000-2010 –> 2011-2019	2.9574403	140631.3481
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	18875	1.4125181	66616.79683	8.7914516	2000-2010 –> 2011-2019	2.9579606	140644.9947
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	18875	1.4124941	66616.15740	8.7917726	2000-2010 –> 2011-2019	2.9579155	140643.8067
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	18875	1.4124941	66616.15740	8.7917726	2000-2010 –> 2011-2019	2.9579155	140643.8067
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	18875	1.4124941	66616.15740	8.7917726	2000-2010 –> 2011-2019	2.9579155	140643.8067
Indo-Pacific	(28.1, Inf]	4	cstar_tref ~ sal + temp + silicate + phosphate_star	18875	1.4266329	66990.14771	8.7576692	2000-2010 –> 2011-2019	2.9813091	141254.1251
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	18875	1.4525020	67670.53618	8.6682222	2000-2010 –> 2011-2019	3.0326233	142584.5202

5.3.2 per fitting unit

lm_best %>%
  select(-n_predictors) %>%
  group_by(basin, gamma_slab, eras) %>% 
  summarise(rmse_sum_mean = mean(rmse_sum),
            ais_sum_mean = mean(aic_sum)) %>% 
  ungroup() %>% 
  kable() %>%
  add_header_above() %>%
  kable_styling() %>%
  scroll_box(width = "100%", height = "400px")

basin	gamma_slab	eras	rmse_sum_mean	ais_sum_mean
Atlantic	(-Inf,26]	1982-1999 –> 2000-2010	2.162495	289.4655
Atlantic	(-Inf,26]	2000-2010 –> 2011-2019	1.022836	139.4740
Atlantic	(26,26.5]	1982-1999 –> 2000-2010	7.226852	3349.8744
Atlantic	(26,26.5]	2000-2010 –> 2011-2019	7.544815	3069.0497
Atlantic	(26.5,26.75]	1982-1999 –> 2000-2010	6.942881	4319.2027
Atlantic	(26.5,26.75]	2000-2010 –> 2011-2019	8.185839	3991.4812
Atlantic	(26.75,27]	1982-1999 –> 2000-2010	4.023536	6578.5296
Atlantic	(26.75,27]	2000-2010 –> 2011-2019	5.246195	5828.9107
Atlantic	(27,27.25]	1982-1999 –> 2000-2010	3.434633	7618.4198
Atlantic	(27,27.25]	2000-2010 –> 2011-2019	4.425692	7088.3754
Atlantic	(27.25,27.5]	1982-1999 –> 2000-2010	4.777802	14055.5795
Atlantic	(27.25,27.5]	2000-2010 –> 2011-2019	6.914383	15273.2846
Atlantic	(27.5,27.75]	1982-1999 –> 2000-2010	3.990971	25826.2297
Atlantic	(27.5,27.75]	2000-2010 –> 2011-2019	4.572914	24012.2939
Atlantic	(27.75,27.85]	1982-1999 –> 2000-2010	1.826732	6372.6466
Atlantic	(27.75,27.85]	2000-2010 –> 2011-2019	2.351874	6521.8320
Atlantic	(27.85,27.95]	1982-1999 –> 2000-2010	2.429179	7510.5355
Atlantic	(27.85,27.95]	2000-2010 –> 2011-2019	2.763133	6747.2659
Atlantic	(27.95,28.05]	1982-1999 –> 2000-2010	7.528814	15730.7633
Atlantic	(27.95,28.05]	2000-2010 –> 2011-2019	8.080876	13204.2336
Atlantic	(28.05,28.1]	1982-1999 –> 2000-2010	1.504288	3719.7318
Atlantic	(28.05,28.1]	2000-2010 –> 2011-2019	1.974389	3710.8279
Atlantic	(28.1,28.15]	1982-1999 –> 2000-2010	1.409077	4915.7689
Atlantic	(28.1,28.15]	2000-2010 –> 2011-2019	1.613097	4806.7106
Atlantic	(28.15,28.2]	1982-1999 –> 2000-2010	1.632437	9777.2487
Atlantic	(28.15,28.2]	2000-2010 –> 2011-2019	1.506991	8947.7068
Atlantic	(28.2, Inf]	1982-1999 –> 2000-2010	1.690325	18341.7982
Atlantic	(28.2, Inf]	2000-2010 –> 2011-2019	1.523553	15821.7291
Indo-Pacific	(-Inf,26]	1982-1999 –> 2000-2010	11.043223	20148.1947
Indo-Pacific	(-Inf,26]	2000-2010 –> 2011-2019	11.201063	29893.6516
Indo-Pacific	(26,26.5]	1982-1999 –> 2000-2010	8.599658	22641.1796
Indo-Pacific	(26,26.5]	2000-2010 –> 2011-2019	8.316578	23320.4017
Indo-Pacific	(26.5,26.75]	1982-1999 –> 2000-2010	8.003554	19092.6045
Indo-Pacific	(26.5,26.75]	2000-2010 –> 2011-2019	7.746122	20458.9421
Indo-Pacific	(26.75,27]	1982-1999 –> 2000-2010	7.535828	27698.1523
Indo-Pacific	(26.75,27]	2000-2010 –> 2011-2019	8.086711	29109.1765
Indo-Pacific	(27,27.25]	1982-1999 –> 2000-2010	6.813550	33556.3011
Indo-Pacific	(27,27.25]	2000-2010 –> 2011-2019	7.973024	39192.6908
Indo-Pacific	(27.25,27.5]	1982-1999 –> 2000-2010	5.956804	40734.3451
Indo-Pacific	(27.25,27.5]	2000-2010 –> 2011-2019	6.641875	47670.9964
Indo-Pacific	(27.5,27.75]	1982-1999 –> 2000-2010	5.328818	45448.6983
Indo-Pacific	(27.5,27.75]	2000-2010 –> 2011-2019	5.730364	52187.5966
Indo-Pacific	(27.75,27.85]	1982-1999 –> 2000-2010	4.649641	17067.8277
Indo-Pacific	(27.75,27.85]	2000-2010 –> 2011-2019	5.386419	19968.8324
Indo-Pacific	(27.85,27.95]	1982-1999 –> 2000-2010	4.267564	18918.4798
Indo-Pacific	(27.85,27.95]	2000-2010 –> 2011-2019	6.040150	25034.8503
Indo-Pacific	(27.95,28.05]	1982-1999 –> 2000-2010	3.316507	18892.9495
Indo-Pacific	(27.95,28.05]	2000-2010 –> 2011-2019	3.209649	22974.7993
Indo-Pacific	(28.05,28.1]	1982-1999 –> 2000-2010	3.106094	10737.9859
Indo-Pacific	(28.05,28.1]	2000-2010 –> 2011-2019	2.847170	12216.7131
Indo-Pacific	(28.1, Inf]	1982-1999 –> 2000-2010	2.869625	117339.5031
Indo-Pacific	(28.1, Inf]	2000-2010 –> 2011-2019	2.967401	140890.4048

5.4 Target variable coefficients

A data frame to map the target variable is prepared.

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

eras_backward <- lm_all_fitted_wide %>%
  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_target <- full_join(
  lm_best %>% select(basin, gamma_slab, model, eras),
  eras_era)

lm_best_target <- left_join(lm_best_target, lm_all_fitted_wide)

rm(eras_era, eras_forward, eras_backward,
   lm_all_fitted)

5.5 Plot selected model residuals

# plot model diagnostics, if activated
if (params_local$plot_all_figures == "y") {
  # mutate predictors column
  lm_best_plot <- lm_best_target %>%
    select(basin, gamma_slab, model, eras, era) %>%
    mutate(
      predictors = str_remove(model, paste(params_local$MLR_target, "~ ")),
      predictors = str_replace_all(predictors, "\\+ ", "")
    )
  
  # loop across all basins, era, gamma slabs, and MLRs
  # fit all MLR models
  for (i_basin in unique(GLODAP$basin)) {
    for (i_era in unique(GLODAP$era)) {
      #i_basin <- unique(GLODAP$basin)[1]
      #i_era   <- unique(GLODAP$era)[2]
      print(i_basin)
      print(i_era)
      
      GLODAP_basin_era <- GLODAP %>%
        filter(basin == i_basin, era == i_era)
      
      for (i_gamma_slab in unique(GLODAP_basin_era$gamma_slab)) {
        #i_gamma_slab <- unique(GLODAP_basin_era$gamma_slab)[1]
        print(i_gamma_slab)
        
        GLODAP_basin_era_slab <- GLODAP_basin_era %>%
          filter(gamma_slab == i_gamma_slab)
        
        lm_best_basin_era_slab <- lm_best_plot %>%
          filter(basin == i_basin, era == i_era, gamma_slab == i_gamma_slab)
        
        for (i_eras in unique(lm_best_basin_era_slab$eras)) {
          #i_eras <- unique(lm_best_basin_era_slab$eras)[1]
          print(i_eras)
          
          lm_best_basin_era_slab_eras <- lm_best_basin_era_slab %>%
            filter(eras == i_eras)
          
          for (i_predictors in unique(lm_best_basin_era_slab_eras$predictors)) {
            #i_predictors <- unique(lm_best_basin_era_slab$predictors)[1]
            print(i_predictors)
            # extract one model definition
            i_lm <- lm_all %>%
              filter(predictors == i_predictors) %>%
              select(lm_coeff) %>%
              pull()
            
             # fit model
          if (params_local$MLR_type == "rlm") {
            i_lm_fit <- MASS::rlm(as.formula(i_lm),
                                  data = GLODAP_basin_era_slab)
          }
          
          if (params_local$MLR_type == "lm") {
            i_lm_fit <- lm(as.formula(i_lm),
                           data = GLODAP_basin_era_slab)
          }
            
            # plot model diagnostics vs predictors
            
            p_model_predictors <- ggnostic(
              i_lm_fit,
              columnsY = c(params_local$MLR_target, ".fitted", ".resid"),
              title = paste(
                "era:",
                i_era,
                "| eras:",
                i_eras,
                "| basin:",
                i_basin,
                "| gamma slab:",
                i_gamma_slab,
                "| predictors:",
                i_predictors
              )
            )
            
            ggsave(
              plot = p_model_predictors,
              path = paste(path_version_figures, "eMLR_diagnostics/", sep = ""),
              filename = paste(
                "MLR_residuals",
                i_era,
                i_eras,
                i_basin,
                i_gamma_slab,
                i_predictors,
                "predictors.png",
                sep = "_"
              ),
              width = 14,
              height = 8
            )

            rm(p_model_predictors)
            
            
            # plot model diagnostics vs location
            
            GLODAP_basin_era_slab <- GLODAP_basin_era_slab %>%
              mutate(fitted = i_lm_fit$fitted.values,
                     residuals = i_lm_fit$residuals)
            
            GLODAP_basin_era_slab_long <- GLODAP_basin_era_slab %>% 
              pivot_longer(cols = c(params_local$MLR_target , fitted, residuals),
                           names_to = "estimate",
                           values_to = "value"
              ) %>% 
              pivot_longer(cols = c(lat, lon, depth),
                           names_to = "coordinate_type",
                           values_to = "coordinate_value"
              )

             p_model_coordinate <- GLODAP_basin_era_slab_long %>%
              ggplot(aes(coordinate_value, value)) +
              geom_bin2d() +
              scale_fill_viridis_c() +
              labs(
                title = paste(
                  "era:",
                  i_era,
                  "| eras:",
                  i_eras,
                  "| basin:",
                  i_basin,
                  "| gamma slab:",
                  i_gamma_slab,
                  "| predictors:",
                  i_predictors
                )
              ) +
              facet_grid(estimate~coordinate_type,
                         scales = "free")

            ggsave(
              plot = p_model_coordinate,
              path = paste(path_version_figures, "eMLR_diagnostics", sep = ""),
              filename = paste(
                "Location_MLR_residuals",
                i_era,
                i_eras,
                i_basin,
                i_gamma_slab,
                i_predictors,
                "predictors.png",
                sep = "_"
              ),
              width = 14,
              height = 8
            )

            rm(p_model_coordinate)
            
          }
          
        }
        
      }
    }
  }
  
  rm(
    lm_best_plot,
    lm_best_basin_era_slab,
    i_rmse,
    GLODAP_basin_era,
    GLODAP_basin_era_slab,
    i_lm,
    lm_all_fitted,
    i_basin,
    i_era,
    i_gamma_slab,
    i_predictors,
    lm_all,
    i_aic,
    i_n_predictors,
    i_resid_max
  )
  
}

Individual residual plots of the MLR models for each basin, era, eras and neutral density (gamma) slab are available at:

/nfs/kryo/work/jenmueller/emlr_cant/model/v_XXX/figures/eMLR_diagnostics/

5.6 Cant coeffcients

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

# pivot long format
lm_best_long <- lm_best_target %>%
  pivot_longer(cols = starts_with("coeff_"),
               names_to = "term",
               values_to = "estimate",
               names_prefix = "coeff_")

# subtract coefficients of adjacent era
lm_best_long <- lm_best_long %>%
  arrange(era) %>%
  group_by(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 back to wide format
lm_best_cant <- lm_best_long %>%
  pivot_wider(values_from = delta_coeff,
              names_from = term,
              names_prefix = "delta_coeff_",
              values_fill = 0)

5.7 Write files

lm_best_target %>%
  select(
    basin,
    gamma_slab,
    model,
    eras,
    era,
    starts_with("coeff_")
  ) %>%
  write_csv(paste(path_version_data,
                  "lm_best_target.csv",
                  sep = ""))

lm_best_cant %>%
  select(
    basin,
    gamma_slab,
    model,
    eras,
    starts_with("delta_coeff_")
  ) %>%
  write_csv(paste(path_version_data,
                  "lm_best_cant.csv",
                  sep = ""))

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

6 Model diagnotics

6.1 Selection criterion vs predictors

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

6.1.1 All models

lm_all_fitted_wide %>%
  ggplot(aes(as.factor(n_predictors),
             !!sym(params_local$MLR_criterion),
             col = basin)) +
  geom_hline(yintercept = 10) +
  geom_boxplot() +
  facet_grid(gamma_slab~era) +
  scale_color_brewer(palette = "Set1") +
  labs(x="Number of predictors")

Past versions of eMLR_predictor_model_quality_all-1.png

Version	Author	Date
9dbf5bd	Donghe-Zhu	2021-03-11
a49df30	Donghe-Zhu	2021-03-11
b3348a5	Donghe-Zhu	2021-03-11
1c24ff7	Donghe-Zhu	2021-03-10
94ce5a8	Donghe-Zhu	2021-03-10
db33928	Donghe-Zhu	2021-03-10
3d1f470	Donghe-Zhu	2021-03-10
9b7bc66	Donghe-Zhu	2021-03-10
4168b43	Donghe-Zhu	2021-03-10
5365f80	Donghe-Zhu	2021-03-10
2b6c392	Donghe-Zhu	2021-03-10
9f58753	Donghe-Zhu	2021-03-10
9d5a62c	Donghe-Zhu	2021-03-10
5d1e70b	Donghe-Zhu	2021-03-10
2093979	Donghe-Zhu	2021-03-10
b865899	Donghe-Zhu	2021-03-10
cc2a956	Donghe-Zhu	2021-03-10
60689fb	Donghe-Zhu	2021-03-10
dba33c8	Donghe-Zhu	2021-03-09
17f1c4a	Donghe-Zhu	2021-03-09
c024d1a	Donghe-Zhu	2021-03-09
02f7242	Donghe-Zhu	2021-03-09
6f50bc6	Donghe-Zhu	2021-03-09
1691156	Donghe-Zhu	2021-03-08
c0ceaf8	Donghe-Zhu	2021-03-08
058e0a1	Donghe-Zhu	2021-03-08
112dea0	Donghe-Zhu	2021-03-08
1843412	Donghe-Zhu	2021-03-08
3fbbfa4	Donghe-Zhu	2021-03-07
627c8fb	Donghe-Zhu	2021-03-07
9ef3222	Donghe-Zhu	2021-03-05
8c1e978	Donghe-Zhu	2021-03-05
865f68c	Donghe-Zhu	2021-03-05
ee69bc1	Donghe-Zhu	2021-03-05
a79291f	Donghe-Zhu	2021-03-05
e8c6f30	Donghe-Zhu	2021-03-04
59288fe	Donghe-Zhu	2021-03-04
731abc8	Donghe-Zhu	2021-03-04
e2a5a33	Donghe-Zhu	2021-03-04
c7892c1	Donghe-Zhu	2021-03-04
924430b	Donghe-Zhu	2021-03-03
0d0bca1	Donghe-Zhu	2021-03-03
cb63c16	Donghe-Zhu	2021-03-03
ffda45a	Donghe-Zhu	2021-03-03
691ba81	Donghe-Zhu	2021-03-03
c5e45a2	Donghe-Zhu	2021-03-03
89c3e58	Donghe-Zhu	2021-03-03
c911669	Donghe-Zhu	2021-03-03
b71c719	Donghe-Zhu	2021-03-01
13666ca	Donghe-Zhu	2021-03-01
c6e60fe	Donghe-Zhu	2021-03-01
7a388f7	Donghe-Zhu	2021-03-01
799e913	Donghe-Zhu	2021-03-01
66ff99f	Donghe-Zhu	2021-03-01
ac9bb7a	Donghe-Zhu	2021-02-28
efdc047	Donghe-Zhu	2021-02-28
e9a7418	Donghe-Zhu	2021-02-28
e152917	Donghe-Zhu	2021-02-28
287123c	Donghe-Zhu	2021-02-27
54d5b5b	Donghe-Zhu	2021-02-27
330f064	Donghe-Zhu	2021-02-27
adbc9bc	Donghe-Zhu	2021-02-27
5937141	Donghe-Zhu	2021-02-27
4414bbf	Donghe-Zhu	2021-02-27
a265efb	Donghe-Zhu	2021-02-27
19edd1e	Donghe-Zhu	2021-02-27
f20483f	Donghe-Zhu	2021-02-26
6a2c7b3	Donghe-Zhu	2021-02-25
02b976d	Donghe-Zhu	2021-02-24
354c224	Donghe-Zhu	2021-02-24
1a0a88a	Donghe-Zhu	2021-02-24
57f701e	Donghe-Zhu	2021-02-24
06f3149	Donghe-Zhu	2021-02-16
401eab3	Donghe-Zhu	2021-02-15
e3bba84	Donghe-Zhu	2021-02-15
5dce4b1	Donghe-Zhu	2021-02-15
4469a0c	Donghe-Zhu	2021-02-13
5ae6a69	Donghe-Zhu	2021-02-10
05385dc	Donghe-Zhu	2021-02-10
f791ae4	Donghe-Zhu	2021-02-09
f71ae34	Donghe-Zhu	2021-02-09
c011832	Donghe-Zhu	2021-02-09
a145fa7	Donghe-Zhu	2021-02-09
c344e42	Donghe-Zhu	2021-02-08
2f095d7	Donghe-Zhu	2021-02-07
1fad5f1	Donghe-Zhu	2021-02-07
ca03c39	Donghe-Zhu	2021-02-07
cd7c52c	Donghe-Zhu	2021-02-04
bcf84f4	Donghe-Zhu	2021-02-02
a518739	Donghe-Zhu	2021-02-01
61666de	Donghe-Zhu	2021-01-31
865b582	Donghe-Zhu	2021-01-31
3e68089	Donghe-Zhu	2021-01-31
ecf335c	Donghe-Zhu	2021-01-31
a618965	Donghe-Zhu	2021-01-31
59e006e	Donghe-Zhu	2021-01-31
a1c8f87	Donghe-Zhu	2021-01-31
ae5c18f	Donghe-Zhu	2021-01-31
b50fe52	Donghe-Zhu	2021-01-31
ac99ae5	jens-daniel-mueller	2021-01-29
b5bdcaf	Donghe-Zhu	2021-01-29
442010d	Donghe-Zhu	2021-01-29
372adf5	Donghe-Zhu	2021-01-29
af8788e	Donghe-Zhu	2021-01-29
21c91c9	Donghe-Zhu	2021-01-29
eded038	Donghe-Zhu	2021-01-29
541d4dd	Donghe-Zhu	2021-01-29
6a75576	Donghe-Zhu	2021-01-28
16fba40	Donghe-Zhu	2021-01-28
12bc567	Donghe-Zhu	2021-01-27
ceed31b	Donghe-Zhu	2021-01-27
342402d	Donghe-Zhu	2021-01-27
5bad5c2	Donghe-Zhu	2021-01-27
61efb56	Donghe-Zhu	2021-01-25
48f638e	Donghe-Zhu	2021-01-25
c1cec47	Donghe-Zhu	2021-01-25
05ffb0c	Donghe-Zhu	2021-01-25
8b97165	Donghe-Zhu	2021-01-25
c569946	Donghe-Zhu	2021-01-24
a2f0d56	Donghe-Zhu	2021-01-23
28509fc	Donghe-Zhu	2021-01-23
4c28e4a	Donghe-Zhu	2021-01-22
24cc264	jens-daniel-mueller	2021-01-22
7891955	Donghe-Zhu	2021-01-21
d4cf1cb	Donghe-Zhu	2021-01-21
1f3e5b6	jens-daniel-mueller	2021-01-20
0e7bdf1	jens-daniel-mueller	2021-01-15
4571843	jens-daniel-mueller	2021-01-14
b3564aa	jens-daniel-mueller	2021-01-14
8d032c3	jens-daniel-mueller	2021-01-14
17dee1d	jens-daniel-mueller	2021-01-13
7cdea0c	jens-daniel-mueller	2021-01-06
fa85b93	jens-daniel-mueller	2021-01-06
e5cb81a	Donghe-Zhu	2021-01-05
a499f10	Donghe-Zhu	2021-01-05
fb8a752	Donghe-Zhu	2020-12-23
8fae0b2	Donghe-Zhu	2020-12-21
c8b76b3	jens-daniel-mueller	2020-12-19

6.1.2 Best models

lm_best_target %>%
  ggplot(aes("",
             !!sym(params_local$MLR_criterion),
             col = basin)) +
  geom_hline(yintercept = 10) +
  geom_boxplot() +
  facet_grid(gamma_slab~era) +
  scale_color_brewer(palette = "Set1") +
  labs(x="Number of predictors pooled")

Past versions of eMLR_predictor_model_quality_best-1.png

Version	Author	Date
9dbf5bd	Donghe-Zhu	2021-03-11
a49df30	Donghe-Zhu	2021-03-11
b3348a5	Donghe-Zhu	2021-03-11
1c24ff7	Donghe-Zhu	2021-03-10
94ce5a8	Donghe-Zhu	2021-03-10
db33928	Donghe-Zhu	2021-03-10
3d1f470	Donghe-Zhu	2021-03-10
9b7bc66	Donghe-Zhu	2021-03-10
4168b43	Donghe-Zhu	2021-03-10
5365f80	Donghe-Zhu	2021-03-10
2b6c392	Donghe-Zhu	2021-03-10
9f58753	Donghe-Zhu	2021-03-10
fd528ed	Donghe-Zhu	2021-03-10
9d5a62c	Donghe-Zhu	2021-03-10
5d1e70b	Donghe-Zhu	2021-03-10
2093979	Donghe-Zhu	2021-03-10
b865899	Donghe-Zhu	2021-03-10
cc2a956	Donghe-Zhu	2021-03-10
60689fb	Donghe-Zhu	2021-03-10
dba33c8	Donghe-Zhu	2021-03-09
17f1c4a	Donghe-Zhu	2021-03-09
c024d1a	Donghe-Zhu	2021-03-09
02f7242	Donghe-Zhu	2021-03-09
6f50bc6	Donghe-Zhu	2021-03-09
1691156	Donghe-Zhu	2021-03-08
c0ceaf8	Donghe-Zhu	2021-03-08
058e0a1	Donghe-Zhu	2021-03-08
112dea0	Donghe-Zhu	2021-03-08
1843412	Donghe-Zhu	2021-03-08
3fbbfa4	Donghe-Zhu	2021-03-07
627c8fb	Donghe-Zhu	2021-03-07
9ef3222	Donghe-Zhu	2021-03-05
8c1e978	Donghe-Zhu	2021-03-05
865f68c	Donghe-Zhu	2021-03-05
ee69bc1	Donghe-Zhu	2021-03-05
a79291f	Donghe-Zhu	2021-03-05
e8c6f30	Donghe-Zhu	2021-03-04
59288fe	Donghe-Zhu	2021-03-04
731abc8	Donghe-Zhu	2021-03-04
e2a5a33	Donghe-Zhu	2021-03-04
c7892c1	Donghe-Zhu	2021-03-04
924430b	Donghe-Zhu	2021-03-03
0d0bca1	Donghe-Zhu	2021-03-03
cb63c16	Donghe-Zhu	2021-03-03
ffda45a	Donghe-Zhu	2021-03-03
691ba81	Donghe-Zhu	2021-03-03
c5e45a2	Donghe-Zhu	2021-03-03
89c3e58	Donghe-Zhu	2021-03-03
c407a50	Donghe-Zhu	2021-03-03
c911669	Donghe-Zhu	2021-03-03
b71c719	Donghe-Zhu	2021-03-01
13666ca	Donghe-Zhu	2021-03-01
c6e60fe	Donghe-Zhu	2021-03-01
7a388f7	Donghe-Zhu	2021-03-01
799e913	Donghe-Zhu	2021-03-01
66ff99f	Donghe-Zhu	2021-03-01
ac9bb7a	Donghe-Zhu	2021-02-28
efdc047	Donghe-Zhu	2021-02-28
e9a7418	Donghe-Zhu	2021-02-28
e152917	Donghe-Zhu	2021-02-28
287123c	Donghe-Zhu	2021-02-27
54d5b5b	Donghe-Zhu	2021-02-27
330f064	Donghe-Zhu	2021-02-27
adbc9bc	Donghe-Zhu	2021-02-27
5937141	Donghe-Zhu	2021-02-27
4414bbf	Donghe-Zhu	2021-02-27
a265efb	Donghe-Zhu	2021-02-27
19edd1e	Donghe-Zhu	2021-02-27
f20483f	Donghe-Zhu	2021-02-26
6a2c7b3	Donghe-Zhu	2021-02-25
354c224	Donghe-Zhu	2021-02-24
1a0a88a	Donghe-Zhu	2021-02-24
57f701e	Donghe-Zhu	2021-02-24
06f3149	Donghe-Zhu	2021-02-16
401eab3	Donghe-Zhu	2021-02-15
e3bba84	Donghe-Zhu	2021-02-15
5dce4b1	Donghe-Zhu	2021-02-15
4469a0c	Donghe-Zhu	2021-02-13
5ae6a69	Donghe-Zhu	2021-02-10
05385dc	Donghe-Zhu	2021-02-10
f791ae4	Donghe-Zhu	2021-02-09
f71ae34	Donghe-Zhu	2021-02-09
c011832	Donghe-Zhu	2021-02-09
a145fa7	Donghe-Zhu	2021-02-09
c344e42	Donghe-Zhu	2021-02-08
2f095d7	Donghe-Zhu	2021-02-07
1fad5f1	Donghe-Zhu	2021-02-07
ca03c39	Donghe-Zhu	2021-02-07
e2ffc14	Donghe-Zhu	2021-02-05
cd7c52c	Donghe-Zhu	2021-02-04
bcf84f4	Donghe-Zhu	2021-02-02
a518739	Donghe-Zhu	2021-02-01
61666de	Donghe-Zhu	2021-01-31
865b582	Donghe-Zhu	2021-01-31
3e68089	Donghe-Zhu	2021-01-31
ecf335c	Donghe-Zhu	2021-01-31
a618965	Donghe-Zhu	2021-01-31
59e006e	Donghe-Zhu	2021-01-31
a1c8f87	Donghe-Zhu	2021-01-31
ae5c18f	Donghe-Zhu	2021-01-31
b50fe52	Donghe-Zhu	2021-01-31
ac99ae5	jens-daniel-mueller	2021-01-29
b5bdcaf	Donghe-Zhu	2021-01-29
442010d	Donghe-Zhu	2021-01-29
372adf5	Donghe-Zhu	2021-01-29
af8788e	Donghe-Zhu	2021-01-29
21c91c9	Donghe-Zhu	2021-01-29
eded038	Donghe-Zhu	2021-01-29
541d4dd	Donghe-Zhu	2021-01-29
6a75576	Donghe-Zhu	2021-01-28
16fba40	Donghe-Zhu	2021-01-28
12bc567	Donghe-Zhu	2021-01-27
ceed31b	Donghe-Zhu	2021-01-27
342402d	Donghe-Zhu	2021-01-27
5bad5c2	Donghe-Zhu	2021-01-27
61efb56	Donghe-Zhu	2021-01-25
48f638e	Donghe-Zhu	2021-01-25
c1cec47	Donghe-Zhu	2021-01-25
05ffb0c	Donghe-Zhu	2021-01-25
8b97165	Donghe-Zhu	2021-01-25
c569946	Donghe-Zhu	2021-01-24
a2f0d56	Donghe-Zhu	2021-01-23
28509fc	Donghe-Zhu	2021-01-23
4c28e4a	Donghe-Zhu	2021-01-22
24cc264	jens-daniel-mueller	2021-01-22
7891955	Donghe-Zhu	2021-01-21
d4cf1cb	Donghe-Zhu	2021-01-21
1f3e5b6	jens-daniel-mueller	2021-01-20
0e7bdf1	jens-daniel-mueller	2021-01-15
4571843	jens-daniel-mueller	2021-01-14
b3564aa	jens-daniel-mueller	2021-01-14
8d032c3	jens-daniel-mueller	2021-01-14
17dee1d	jens-daniel-mueller	2021-01-13
7cdea0c	jens-daniel-mueller	2021-01-06
fa85b93	jens-daniel-mueller	2021-01-06
e5cb81a	Donghe-Zhu	2021-01-05
a499f10	Donghe-Zhu	2021-01-05
fb8a752	Donghe-Zhu	2020-12-23
8fae0b2	Donghe-Zhu	2020-12-21
c8b76b3	jens-daniel-mueller	2020-12-19

6.2 RMSE correlation between eras

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

6.2.1 All models

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

lm_all_fitted_wide_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)) +
  facet_grid(eras ~ basin)

Past versions of rmse_correlation_across_eras_all_models-1.png

Version	Author	Date
9dbf5bd	Donghe-Zhu	2021-03-11
a49df30	Donghe-Zhu	2021-03-11
b3348a5	Donghe-Zhu	2021-03-11
1c24ff7	Donghe-Zhu	2021-03-10
94ce5a8	Donghe-Zhu	2021-03-10
db33928	Donghe-Zhu	2021-03-10
3d1f470	Donghe-Zhu	2021-03-10
9b7bc66	Donghe-Zhu	2021-03-10
5365f80	Donghe-Zhu	2021-03-10
2b6c392	Donghe-Zhu	2021-03-10
9f58753	Donghe-Zhu	2021-03-10
9d5a62c	Donghe-Zhu	2021-03-10
5d1e70b	Donghe-Zhu	2021-03-10
2093979	Donghe-Zhu	2021-03-10
b865899	Donghe-Zhu	2021-03-10
cc2a956	Donghe-Zhu	2021-03-10
dba33c8	Donghe-Zhu	2021-03-09
17f1c4a	Donghe-Zhu	2021-03-09
c024d1a	Donghe-Zhu	2021-03-09
02f7242	Donghe-Zhu	2021-03-09
1691156	Donghe-Zhu	2021-03-08
c0ceaf8	Donghe-Zhu	2021-03-08
058e0a1	Donghe-Zhu	2021-03-08
112dea0	Donghe-Zhu	2021-03-08
1843412	Donghe-Zhu	2021-03-08
3fbbfa4	Donghe-Zhu	2021-03-07
627c8fb	Donghe-Zhu	2021-03-07
9ef3222	Donghe-Zhu	2021-03-05
8c1e978	Donghe-Zhu	2021-03-05
865f68c	Donghe-Zhu	2021-03-05
ee69bc1	Donghe-Zhu	2021-03-05
a79291f	Donghe-Zhu	2021-03-05
e8c6f30	Donghe-Zhu	2021-03-04
59288fe	Donghe-Zhu	2021-03-04
731abc8	Donghe-Zhu	2021-03-04
e2a5a33	Donghe-Zhu	2021-03-04
c7892c1	Donghe-Zhu	2021-03-04
924430b	Donghe-Zhu	2021-03-03
0d0bca1	Donghe-Zhu	2021-03-03
cb63c16	Donghe-Zhu	2021-03-03
ffda45a	Donghe-Zhu	2021-03-03
691ba81	Donghe-Zhu	2021-03-03
c5e45a2	Donghe-Zhu	2021-03-03
89c3e58	Donghe-Zhu	2021-03-03
c911669	Donghe-Zhu	2021-03-03
13666ca	Donghe-Zhu	2021-03-01
c6e60fe	Donghe-Zhu	2021-03-01
7a388f7	Donghe-Zhu	2021-03-01
799e913	Donghe-Zhu	2021-03-01
66ff99f	Donghe-Zhu	2021-03-01
ac9bb7a	Donghe-Zhu	2021-02-28
efdc047	Donghe-Zhu	2021-02-28
e9a7418	Donghe-Zhu	2021-02-28
287123c	Donghe-Zhu	2021-02-27
54d5b5b	Donghe-Zhu	2021-02-27
330f064	Donghe-Zhu	2021-02-27
adbc9bc	Donghe-Zhu	2021-02-27
5937141	Donghe-Zhu	2021-02-27
4414bbf	Donghe-Zhu	2021-02-27
a265efb	Donghe-Zhu	2021-02-27
19edd1e	Donghe-Zhu	2021-02-27
f20483f	Donghe-Zhu	2021-02-26
6a2c7b3	Donghe-Zhu	2021-02-25
02b976d	Donghe-Zhu	2021-02-24
354c224	Donghe-Zhu	2021-02-24
1a0a88a	Donghe-Zhu	2021-02-24
57f701e	Donghe-Zhu	2021-02-24
06f3149	Donghe-Zhu	2021-02-16
e3bba84	Donghe-Zhu	2021-02-15
5dce4b1	Donghe-Zhu	2021-02-15
4469a0c	Donghe-Zhu	2021-02-13
5ae6a69	Donghe-Zhu	2021-02-10
05385dc	Donghe-Zhu	2021-02-10
f791ae4	Donghe-Zhu	2021-02-09
f71ae34	Donghe-Zhu	2021-02-09
c011832	Donghe-Zhu	2021-02-09
a145fa7	Donghe-Zhu	2021-02-09
c344e42	Donghe-Zhu	2021-02-08
2f095d7	Donghe-Zhu	2021-02-07
1fad5f1	Donghe-Zhu	2021-02-07
ca03c39	Donghe-Zhu	2021-02-07
cd7c52c	Donghe-Zhu	2021-02-04
bcf84f4	Donghe-Zhu	2021-02-02
61666de	Donghe-Zhu	2021-01-31
865b582	Donghe-Zhu	2021-01-31
3e68089	Donghe-Zhu	2021-01-31
ecf335c	Donghe-Zhu	2021-01-31
a618965	Donghe-Zhu	2021-01-31
59e006e	Donghe-Zhu	2021-01-31
a1c8f87	Donghe-Zhu	2021-01-31
ae5c18f	Donghe-Zhu	2021-01-31
b50fe52	Donghe-Zhu	2021-01-31
ac99ae5	jens-daniel-mueller	2021-01-29
b5bdcaf	Donghe-Zhu	2021-01-29
442010d	Donghe-Zhu	2021-01-29
372adf5	Donghe-Zhu	2021-01-29
af8788e	Donghe-Zhu	2021-01-29
21c91c9	Donghe-Zhu	2021-01-29
eded038	Donghe-Zhu	2021-01-29
541d4dd	Donghe-Zhu	2021-01-29
6a75576	Donghe-Zhu	2021-01-28
16fba40	Donghe-Zhu	2021-01-28
12bc567	Donghe-Zhu	2021-01-27
ceed31b	Donghe-Zhu	2021-01-27
342402d	Donghe-Zhu	2021-01-27
5bad5c2	Donghe-Zhu	2021-01-27
61efb56	Donghe-Zhu	2021-01-25
48f638e	Donghe-Zhu	2021-01-25
c1cec47	Donghe-Zhu	2021-01-25
05ffb0c	Donghe-Zhu	2021-01-25
8b97165	Donghe-Zhu	2021-01-25
c569946	Donghe-Zhu	2021-01-24
a2f0d56	Donghe-Zhu	2021-01-23
28509fc	Donghe-Zhu	2021-01-23
4c28e4a	Donghe-Zhu	2021-01-22
24cc264	jens-daniel-mueller	2021-01-22
7891955	Donghe-Zhu	2021-01-21
d4cf1cb	Donghe-Zhu	2021-01-21
1f3e5b6	jens-daniel-mueller	2021-01-20
0e7bdf1	jens-daniel-mueller	2021-01-15
4571843	jens-daniel-mueller	2021-01-14
b3564aa	jens-daniel-mueller	2021-01-14
8d032c3	jens-daniel-mueller	2021-01-14
17dee1d	jens-daniel-mueller	2021-01-13
7cdea0c	jens-daniel-mueller	2021-01-06
fa85b93	jens-daniel-mueller	2021-01-06
e5cb81a	Donghe-Zhu	2021-01-05
a499f10	Donghe-Zhu	2021-01-05
fb8a752	Donghe-Zhu	2020-12-23
8fae0b2	Donghe-Zhu	2020-12-21
c8b76b3	jens-daniel-mueller	2020-12-19

rm(max_rmse)

6.2.2 Best models

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

lm_best %>%
  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)) +
  facet_grid(eras ~ basin)

Past versions of rmse_correlation_across_eras_best_models-1.png

Version	Author	Date
9dbf5bd	Donghe-Zhu	2021-03-11
a49df30	Donghe-Zhu	2021-03-11
b3348a5	Donghe-Zhu	2021-03-11
1c24ff7	Donghe-Zhu	2021-03-10
94ce5a8	Donghe-Zhu	2021-03-10
db33928	Donghe-Zhu	2021-03-10
3d1f470	Donghe-Zhu	2021-03-10
9b7bc66	Donghe-Zhu	2021-03-10
4168b43	Donghe-Zhu	2021-03-10
5365f80	Donghe-Zhu	2021-03-10
2b6c392	Donghe-Zhu	2021-03-10
9f58753	Donghe-Zhu	2021-03-10
fd528ed	Donghe-Zhu	2021-03-10
9d5a62c	Donghe-Zhu	2021-03-10
5d1e70b	Donghe-Zhu	2021-03-10
2093979	Donghe-Zhu	2021-03-10
b865899	Donghe-Zhu	2021-03-10
cc2a956	Donghe-Zhu	2021-03-10
60689fb	Donghe-Zhu	2021-03-10
dba33c8	Donghe-Zhu	2021-03-09
17f1c4a	Donghe-Zhu	2021-03-09
c024d1a	Donghe-Zhu	2021-03-09
02f7242	Donghe-Zhu	2021-03-09
6f50bc6	Donghe-Zhu	2021-03-09
1691156	Donghe-Zhu	2021-03-08
c0ceaf8	Donghe-Zhu	2021-03-08
058e0a1	Donghe-Zhu	2021-03-08
112dea0	Donghe-Zhu	2021-03-08
1843412	Donghe-Zhu	2021-03-08
3fbbfa4	Donghe-Zhu	2021-03-07
627c8fb	Donghe-Zhu	2021-03-07
9ef3222	Donghe-Zhu	2021-03-05
8c1e978	Donghe-Zhu	2021-03-05
865f68c	Donghe-Zhu	2021-03-05
ee69bc1	Donghe-Zhu	2021-03-05
a79291f	Donghe-Zhu	2021-03-05
e8c6f30	Donghe-Zhu	2021-03-04
59288fe	Donghe-Zhu	2021-03-04
731abc8	Donghe-Zhu	2021-03-04
e2a5a33	Donghe-Zhu	2021-03-04
c7892c1	Donghe-Zhu	2021-03-04
924430b	Donghe-Zhu	2021-03-03
0d0bca1	Donghe-Zhu	2021-03-03
cb63c16	Donghe-Zhu	2021-03-03
ffda45a	Donghe-Zhu	2021-03-03
691ba81	Donghe-Zhu	2021-03-03
c5e45a2	Donghe-Zhu	2021-03-03
89c3e58	Donghe-Zhu	2021-03-03
c407a50	Donghe-Zhu	2021-03-03
c911669	Donghe-Zhu	2021-03-03
b71c719	Donghe-Zhu	2021-03-01
13666ca	Donghe-Zhu	2021-03-01
c6e60fe	Donghe-Zhu	2021-03-01
7a388f7	Donghe-Zhu	2021-03-01
799e913	Donghe-Zhu	2021-03-01
66ff99f	Donghe-Zhu	2021-03-01
ac9bb7a	Donghe-Zhu	2021-02-28
efdc047	Donghe-Zhu	2021-02-28
e9a7418	Donghe-Zhu	2021-02-28
e152917	Donghe-Zhu	2021-02-28
287123c	Donghe-Zhu	2021-02-27
54d5b5b	Donghe-Zhu	2021-02-27
330f064	Donghe-Zhu	2021-02-27
adbc9bc	Donghe-Zhu	2021-02-27
5937141	Donghe-Zhu	2021-02-27
4414bbf	Donghe-Zhu	2021-02-27
a265efb	Donghe-Zhu	2021-02-27
19edd1e	Donghe-Zhu	2021-02-27
f20483f	Donghe-Zhu	2021-02-26
6a2c7b3	Donghe-Zhu	2021-02-25
354c224	Donghe-Zhu	2021-02-24
1a0a88a	Donghe-Zhu	2021-02-24
57f701e	Donghe-Zhu	2021-02-24
06f3149	Donghe-Zhu	2021-02-16
401eab3	Donghe-Zhu	2021-02-15
e3bba84	Donghe-Zhu	2021-02-15
5dce4b1	Donghe-Zhu	2021-02-15
4469a0c	Donghe-Zhu	2021-02-13
5ae6a69	Donghe-Zhu	2021-02-10
05385dc	Donghe-Zhu	2021-02-10
f791ae4	Donghe-Zhu	2021-02-09
f71ae34	Donghe-Zhu	2021-02-09
c011832	Donghe-Zhu	2021-02-09
a145fa7	Donghe-Zhu	2021-02-09
c344e42	Donghe-Zhu	2021-02-08
2f095d7	Donghe-Zhu	2021-02-07
1fad5f1	Donghe-Zhu	2021-02-07
ca03c39	Donghe-Zhu	2021-02-07
e2ffc14	Donghe-Zhu	2021-02-05
cd7c52c	Donghe-Zhu	2021-02-04
bcf84f4	Donghe-Zhu	2021-02-02
a518739	Donghe-Zhu	2021-02-01
61666de	Donghe-Zhu	2021-01-31
865b582	Donghe-Zhu	2021-01-31
3e68089	Donghe-Zhu	2021-01-31
ecf335c	Donghe-Zhu	2021-01-31
a618965	Donghe-Zhu	2021-01-31
59e006e	Donghe-Zhu	2021-01-31
a1c8f87	Donghe-Zhu	2021-01-31
ae5c18f	Donghe-Zhu	2021-01-31
b50fe52	Donghe-Zhu	2021-01-31
ac99ae5	jens-daniel-mueller	2021-01-29
b5bdcaf	Donghe-Zhu	2021-01-29
442010d	Donghe-Zhu	2021-01-29
372adf5	Donghe-Zhu	2021-01-29
af8788e	Donghe-Zhu	2021-01-29
21c91c9	Donghe-Zhu	2021-01-29
eded038	Donghe-Zhu	2021-01-29
541d4dd	Donghe-Zhu	2021-01-29
6a75576	Donghe-Zhu	2021-01-28
16fba40	Donghe-Zhu	2021-01-28
12bc567	Donghe-Zhu	2021-01-27
ceed31b	Donghe-Zhu	2021-01-27
342402d	Donghe-Zhu	2021-01-27
5bad5c2	Donghe-Zhu	2021-01-27
61efb56	Donghe-Zhu	2021-01-25
48f638e	Donghe-Zhu	2021-01-25
c1cec47	Donghe-Zhu	2021-01-25
05ffb0c	Donghe-Zhu	2021-01-25
8b97165	Donghe-Zhu	2021-01-25
c569946	Donghe-Zhu	2021-01-24
a2f0d56	Donghe-Zhu	2021-01-23
28509fc	Donghe-Zhu	2021-01-23
4c28e4a	Donghe-Zhu	2021-01-22
24cc264	jens-daniel-mueller	2021-01-22
7891955	Donghe-Zhu	2021-01-21
d4cf1cb	Donghe-Zhu	2021-01-21
1f3e5b6	jens-daniel-mueller	2021-01-20
0e7bdf1	jens-daniel-mueller	2021-01-15
4571843	jens-daniel-mueller	2021-01-14
b3564aa	jens-daniel-mueller	2021-01-14
8d032c3	jens-daniel-mueller	2021-01-14
17dee1d	jens-daniel-mueller	2021-01-13
7cdea0c	jens-daniel-mueller	2021-01-06
fa85b93	jens-daniel-mueller	2021-01-06
e5cb81a	Donghe-Zhu	2021-01-05
a499f10	Donghe-Zhu	2021-01-05
fb8a752	Donghe-Zhu	2020-12-23
8fae0b2	Donghe-Zhu	2020-12-21
c8b76b3	jens-daniel-mueller	2020-12-19

rm(max_rmse)

6.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_long %>% 
  filter(term != "(Intercept)",
         delta_coeff != 0) %>% 
  group_by(basin, eras, gamma_slab) %>% 
  count(term) %>% 
  ungroup() %>% 
  pivot_wider(values_from = n,
              names_from = term,
              values_fill = 0)

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

	aou	oxygen	phosphate	phosphate_star	sal	silicate	temp
Atlantic - 1982-1999 --> 2000-2010
(-Inf,26]	10	7	6	5	5	4	10
(26,26.5]	6	5	5	6	10	10	7
(26.5,26.75]	10	7	7	6	3	10	3
(26.75,27]	10	7	7	6	3	10	3
(27,27.25]	7	4	7	6	9	10	6
(27.25,27.5]	7	5	6	6	10	9	6
(27.5,27.75]	6	5	6	6	10	9	7
(27.75,27.85]	7	6	7	6	10	7	7
(27.85,27.95]	10	7	7	6	3	3	10
(27.95,28.05]	7	5	7	6	10	8	7
(28.05,28.1]	10	7	7	6	3	3	10
(28.1,28.15]	10	6	7	7	6	6	7
(28.15,28.2]	10	6	7	6	7	7	7
(28.2, Inf]	7	6	7	6	7	10	7
total	117.00	83.00	93.00	84.00	96.00	106.00	97.00
Atlantic - 2000-2010 --> 2011-2019
(-Inf,26]	6	5	6	5	7	10	10
(26,26.5]	7	6	5	5	10	10	6
(26.5,26.75]	10	7	7	6	3	10	3
(26.75,27]	8	6	6	6	7	10	7
(27,27.25]	7	4	6	6	10	10	5
(27.25,27.5]	6	4	6	6	10	10	6
(27.5,27.75]	6	4	6	6	10	10	6
(27.75,27.85]	8	5	7	6	10	7	6
(27.85,27.95]	10	6	7	6	7	7	7
(27.95,28.05]	7	5	7	6	10	8	7
(28.05,28.1]	10	7	7	6	3	3	10
(28.1,28.15]	10	7	7	6	3	3	10
(28.15,28.2]	10	6	7	6	7	7	7
(28.2, Inf]	7	6	7	6	7	10	7
total	112.00	78.00	91.00	82.00	104.00	115.00	97.00
Indo-Pacific - 1982-1999 --> 2000-2010
(-Inf,26]	10	6	8	5	10	4	4
(26,26.5]	7	5	7	6	10	8	7
(26.5,26.75]	10	6	7	6	7	7	7
(26.75,27]	10	7	7	6	10	3	3
(27,27.25]	10	7	7	6	3	3	10
(27.25,27.5]	9	7	5	5	6	10	7
(27.5,27.75]	7	6	5	5	6	10	10
(27.75,27.85]	10	7	6	5	4	5	10
(27.85,27.95]	10	7	7	6	3	3	10
(27.95,28.05]	10	7	7	6	3	3	10
(28.05,28.1]	7	6	7	6	7	10	7
(28.1, Inf]	6	4	6	6	10	10	6
total	106.00	75.00	79.00	68.00	79.00	76.00	91.00
Indo-Pacific - 2000-2010 --> 2011-2019
(-Inf,26]	10	7	7	6	10	3	3
(26,26.5]	7	5	7	6	8	10	7
(26.5,26.75]	10	6	7	6	7	7	7
(26.75,27]	10	7	6	7	6	7	6
(27,27.25]	10	7	7	6	3	3	10
(27.25,27.5]	7	6	5	5	6	10	10
(27.5,27.75]	6	5	6	5	10	10	7
(27.75,27.85]	7	6	6	5	9	10	7
(27.85,27.95]	7	6	6	6	8	10	7
(27.95,28.05]	10	6	7	6	7	7	7
(28.05,28.1]	8	5	7	6	7	10	7
(28.1, Inf]	6	4	6	7	9	10	7
total	98.00	70.00	77.00	71.00	90.00	97.00	85.00

6.4 RMSE alternatives

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

lm_all_fitted_wide_eras %>% 
  ggplot(aes(rmse, aic, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  facet_grid(eras~basin)

Past versions of aic_vs_rmse_all_models-1.png

Version	Author	Date
9dbf5bd	Donghe-Zhu	2021-03-11
a49df30	Donghe-Zhu	2021-03-11
b3348a5	Donghe-Zhu	2021-03-11
1c24ff7	Donghe-Zhu	2021-03-10
94ce5a8	Donghe-Zhu	2021-03-10
db33928	Donghe-Zhu	2021-03-10
3d1f470	Donghe-Zhu	2021-03-10
9b7bc66	Donghe-Zhu	2021-03-10
5365f80	Donghe-Zhu	2021-03-10
2b6c392	Donghe-Zhu	2021-03-10
9f58753	Donghe-Zhu	2021-03-10
9d5a62c	Donghe-Zhu	2021-03-10
5d1e70b	Donghe-Zhu	2021-03-10
2093979	Donghe-Zhu	2021-03-10
b865899	Donghe-Zhu	2021-03-10
cc2a956	Donghe-Zhu	2021-03-10
dba33c8	Donghe-Zhu	2021-03-09
17f1c4a	Donghe-Zhu	2021-03-09
c024d1a	Donghe-Zhu	2021-03-09
02f7242	Donghe-Zhu	2021-03-09
1691156	Donghe-Zhu	2021-03-08
c0ceaf8	Donghe-Zhu	2021-03-08
058e0a1	Donghe-Zhu	2021-03-08
112dea0	Donghe-Zhu	2021-03-08
1843412	Donghe-Zhu	2021-03-08
3fbbfa4	Donghe-Zhu	2021-03-07
627c8fb	Donghe-Zhu	2021-03-07
9ef3222	Donghe-Zhu	2021-03-05
8c1e978	Donghe-Zhu	2021-03-05
865f68c	Donghe-Zhu	2021-03-05
ee69bc1	Donghe-Zhu	2021-03-05
a79291f	Donghe-Zhu	2021-03-05
e8c6f30	Donghe-Zhu	2021-03-04
59288fe	Donghe-Zhu	2021-03-04
731abc8	Donghe-Zhu	2021-03-04
e2a5a33	Donghe-Zhu	2021-03-04
c7892c1	Donghe-Zhu	2021-03-04
924430b	Donghe-Zhu	2021-03-03
0d0bca1	Donghe-Zhu	2021-03-03
cb63c16	Donghe-Zhu	2021-03-03
ffda45a	Donghe-Zhu	2021-03-03
691ba81	Donghe-Zhu	2021-03-03
c5e45a2	Donghe-Zhu	2021-03-03
89c3e58	Donghe-Zhu	2021-03-03
c911669	Donghe-Zhu	2021-03-03
13666ca	Donghe-Zhu	2021-03-01
c6e60fe	Donghe-Zhu	2021-03-01
7a388f7	Donghe-Zhu	2021-03-01
799e913	Donghe-Zhu	2021-03-01
66ff99f	Donghe-Zhu	2021-03-01
ac9bb7a	Donghe-Zhu	2021-02-28
efdc047	Donghe-Zhu	2021-02-28
e9a7418	Donghe-Zhu	2021-02-28
287123c	Donghe-Zhu	2021-02-27
54d5b5b	Donghe-Zhu	2021-02-27
330f064	Donghe-Zhu	2021-02-27
adbc9bc	Donghe-Zhu	2021-02-27
5937141	Donghe-Zhu	2021-02-27
4414bbf	Donghe-Zhu	2021-02-27
a265efb	Donghe-Zhu	2021-02-27
19edd1e	Donghe-Zhu	2021-02-27
f20483f	Donghe-Zhu	2021-02-26
6a2c7b3	Donghe-Zhu	2021-02-25
02b976d	Donghe-Zhu	2021-02-24
354c224	Donghe-Zhu	2021-02-24
1a0a88a	Donghe-Zhu	2021-02-24
57f701e	Donghe-Zhu	2021-02-24
06f3149	Donghe-Zhu	2021-02-16
e3bba84	Donghe-Zhu	2021-02-15
5dce4b1	Donghe-Zhu	2021-02-15
4469a0c	Donghe-Zhu	2021-02-13
5ae6a69	Donghe-Zhu	2021-02-10
05385dc	Donghe-Zhu	2021-02-10
f791ae4	Donghe-Zhu	2021-02-09
f71ae34	Donghe-Zhu	2021-02-09
c011832	Donghe-Zhu	2021-02-09
a145fa7	Donghe-Zhu	2021-02-09
c344e42	Donghe-Zhu	2021-02-08
2f095d7	Donghe-Zhu	2021-02-07
1fad5f1	Donghe-Zhu	2021-02-07
ca03c39	Donghe-Zhu	2021-02-07
cd7c52c	Donghe-Zhu	2021-02-04
bcf84f4	Donghe-Zhu	2021-02-02
61666de	Donghe-Zhu	2021-01-31
865b582	Donghe-Zhu	2021-01-31
3e68089	Donghe-Zhu	2021-01-31
ecf335c	Donghe-Zhu	2021-01-31
a618965	Donghe-Zhu	2021-01-31
59e006e	Donghe-Zhu	2021-01-31
a1c8f87	Donghe-Zhu	2021-01-31
ae5c18f	Donghe-Zhu	2021-01-31
b50fe52	Donghe-Zhu	2021-01-31
ac99ae5	jens-daniel-mueller	2021-01-29
b5bdcaf	Donghe-Zhu	2021-01-29
442010d	Donghe-Zhu	2021-01-29
372adf5	Donghe-Zhu	2021-01-29
af8788e	Donghe-Zhu	2021-01-29
21c91c9	Donghe-Zhu	2021-01-29
eded038	Donghe-Zhu	2021-01-29
541d4dd	Donghe-Zhu	2021-01-29
6a75576	Donghe-Zhu	2021-01-28
16fba40	Donghe-Zhu	2021-01-28
12bc567	Donghe-Zhu	2021-01-27
ceed31b	Donghe-Zhu	2021-01-27
342402d	Donghe-Zhu	2021-01-27
5bad5c2	Donghe-Zhu	2021-01-27
61efb56	Donghe-Zhu	2021-01-25
48f638e	Donghe-Zhu	2021-01-25
c1cec47	Donghe-Zhu	2021-01-25
05ffb0c	Donghe-Zhu	2021-01-25
8b97165	Donghe-Zhu	2021-01-25
c569946	Donghe-Zhu	2021-01-24
a2f0d56	Donghe-Zhu	2021-01-23
28509fc	Donghe-Zhu	2021-01-23
4c28e4a	Donghe-Zhu	2021-01-22
24cc264	jens-daniel-mueller	2021-01-22
7891955	Donghe-Zhu	2021-01-21
d4cf1cb	Donghe-Zhu	2021-01-21
1f3e5b6	jens-daniel-mueller	2021-01-20
0e7bdf1	jens-daniel-mueller	2021-01-15
4571843	jens-daniel-mueller	2021-01-14
b3564aa	jens-daniel-mueller	2021-01-14
8d032c3	jens-daniel-mueller	2021-01-14
17dee1d	jens-daniel-mueller	2021-01-13
7cdea0c	jens-daniel-mueller	2021-01-06
fa85b93	jens-daniel-mueller	2021-01-06
e5cb81a	Donghe-Zhu	2021-01-05
a499f10	Donghe-Zhu	2021-01-05
fb8a752	Donghe-Zhu	2020-12-23
8fae0b2	Donghe-Zhu	2020-12-21
c8b76b3	jens-daniel-mueller	2020-12-19

lm_best %>% 
  ggplot(aes(rmse, aic, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  facet_grid(eras~basin)

Past versions of aic_vs_rmse_selected_models-1.png

Version	Author	Date
9dbf5bd	Donghe-Zhu	2021-03-11
a49df30	Donghe-Zhu	2021-03-11
b3348a5	Donghe-Zhu	2021-03-11
1c24ff7	Donghe-Zhu	2021-03-10
94ce5a8	Donghe-Zhu	2021-03-10
db33928	Donghe-Zhu	2021-03-10
3d1f470	Donghe-Zhu	2021-03-10
9b7bc66	Donghe-Zhu	2021-03-10
4168b43	Donghe-Zhu	2021-03-10
5365f80	Donghe-Zhu	2021-03-10
2b6c392	Donghe-Zhu	2021-03-10
9f58753	Donghe-Zhu	2021-03-10
fd528ed	Donghe-Zhu	2021-03-10
9d5a62c	Donghe-Zhu	2021-03-10
5d1e70b	Donghe-Zhu	2021-03-10
2093979	Donghe-Zhu	2021-03-10
b865899	Donghe-Zhu	2021-03-10
cc2a956	Donghe-Zhu	2021-03-10
60689fb	Donghe-Zhu	2021-03-10
dba33c8	Donghe-Zhu	2021-03-09
17f1c4a	Donghe-Zhu	2021-03-09
c024d1a	Donghe-Zhu	2021-03-09
02f7242	Donghe-Zhu	2021-03-09
6f50bc6	Donghe-Zhu	2021-03-09
1691156	Donghe-Zhu	2021-03-08
c0ceaf8	Donghe-Zhu	2021-03-08
058e0a1	Donghe-Zhu	2021-03-08
112dea0	Donghe-Zhu	2021-03-08
1843412	Donghe-Zhu	2021-03-08
3fbbfa4	Donghe-Zhu	2021-03-07
627c8fb	Donghe-Zhu	2021-03-07
9ef3222	Donghe-Zhu	2021-03-05
8c1e978	Donghe-Zhu	2021-03-05
865f68c	Donghe-Zhu	2021-03-05
ee69bc1	Donghe-Zhu	2021-03-05
a79291f	Donghe-Zhu	2021-03-05
e8c6f30	Donghe-Zhu	2021-03-04
59288fe	Donghe-Zhu	2021-03-04
731abc8	Donghe-Zhu	2021-03-04
e2a5a33	Donghe-Zhu	2021-03-04
c7892c1	Donghe-Zhu	2021-03-04
924430b	Donghe-Zhu	2021-03-03
0d0bca1	Donghe-Zhu	2021-03-03
cb63c16	Donghe-Zhu	2021-03-03
ffda45a	Donghe-Zhu	2021-03-03
691ba81	Donghe-Zhu	2021-03-03
c5e45a2	Donghe-Zhu	2021-03-03
89c3e58	Donghe-Zhu	2021-03-03
c407a50	Donghe-Zhu	2021-03-03
c911669	Donghe-Zhu	2021-03-03
b71c719	Donghe-Zhu	2021-03-01
13666ca	Donghe-Zhu	2021-03-01
c6e60fe	Donghe-Zhu	2021-03-01
7a388f7	Donghe-Zhu	2021-03-01
799e913	Donghe-Zhu	2021-03-01
66ff99f	Donghe-Zhu	2021-03-01
ac9bb7a	Donghe-Zhu	2021-02-28
efdc047	Donghe-Zhu	2021-02-28
e9a7418	Donghe-Zhu	2021-02-28
e152917	Donghe-Zhu	2021-02-28
287123c	Donghe-Zhu	2021-02-27
54d5b5b	Donghe-Zhu	2021-02-27
330f064	Donghe-Zhu	2021-02-27
adbc9bc	Donghe-Zhu	2021-02-27
5937141	Donghe-Zhu	2021-02-27
4414bbf	Donghe-Zhu	2021-02-27
a265efb	Donghe-Zhu	2021-02-27
19edd1e	Donghe-Zhu	2021-02-27
f20483f	Donghe-Zhu	2021-02-26
6a2c7b3	Donghe-Zhu	2021-02-25
354c224	Donghe-Zhu	2021-02-24
1a0a88a	Donghe-Zhu	2021-02-24
57f701e	Donghe-Zhu	2021-02-24
06f3149	Donghe-Zhu	2021-02-16
401eab3	Donghe-Zhu	2021-02-15
e3bba84	Donghe-Zhu	2021-02-15
5dce4b1	Donghe-Zhu	2021-02-15
4469a0c	Donghe-Zhu	2021-02-13
5ae6a69	Donghe-Zhu	2021-02-10
05385dc	Donghe-Zhu	2021-02-10
f791ae4	Donghe-Zhu	2021-02-09
f71ae34	Donghe-Zhu	2021-02-09
c011832	Donghe-Zhu	2021-02-09
a145fa7	Donghe-Zhu	2021-02-09
c344e42	Donghe-Zhu	2021-02-08
2f095d7	Donghe-Zhu	2021-02-07
1fad5f1	Donghe-Zhu	2021-02-07
ca03c39	Donghe-Zhu	2021-02-07
e2ffc14	Donghe-Zhu	2021-02-05
cd7c52c	Donghe-Zhu	2021-02-04
bcf84f4	Donghe-Zhu	2021-02-02
a518739	Donghe-Zhu	2021-02-01
61666de	Donghe-Zhu	2021-01-31
865b582	Donghe-Zhu	2021-01-31
3e68089	Donghe-Zhu	2021-01-31
ecf335c	Donghe-Zhu	2021-01-31
a618965	Donghe-Zhu	2021-01-31
59e006e	Donghe-Zhu	2021-01-31
a1c8f87	Donghe-Zhu	2021-01-31
ae5c18f	Donghe-Zhu	2021-01-31
b50fe52	Donghe-Zhu	2021-01-31
ac99ae5	jens-daniel-mueller	2021-01-29
b5bdcaf	Donghe-Zhu	2021-01-29
442010d	Donghe-Zhu	2021-01-29
372adf5	Donghe-Zhu	2021-01-29
af8788e	Donghe-Zhu	2021-01-29
21c91c9	Donghe-Zhu	2021-01-29
eded038	Donghe-Zhu	2021-01-29
541d4dd	Donghe-Zhu	2021-01-29
6a75576	Donghe-Zhu	2021-01-28
16fba40	Donghe-Zhu	2021-01-28
12bc567	Donghe-Zhu	2021-01-27
ceed31b	Donghe-Zhu	2021-01-27
342402d	Donghe-Zhu	2021-01-27
5bad5c2	Donghe-Zhu	2021-01-27
61efb56	Donghe-Zhu	2021-01-25
48f638e	Donghe-Zhu	2021-01-25
c1cec47	Donghe-Zhu	2021-01-25
05ffb0c	Donghe-Zhu	2021-01-25
8b97165	Donghe-Zhu	2021-01-25
c569946	Donghe-Zhu	2021-01-24
a2f0d56	Donghe-Zhu	2021-01-23
28509fc	Donghe-Zhu	2021-01-23
4c28e4a	Donghe-Zhu	2021-01-22
24cc264	jens-daniel-mueller	2021-01-22
7891955	Donghe-Zhu	2021-01-21
d4cf1cb	Donghe-Zhu	2021-01-21
1f3e5b6	jens-daniel-mueller	2021-01-20
0e7bdf1	jens-daniel-mueller	2021-01-15
4571843	jens-daniel-mueller	2021-01-14
b3564aa	jens-daniel-mueller	2021-01-14
8d032c3	jens-daniel-mueller	2021-01-14
17dee1d	jens-daniel-mueller	2021-01-13
7cdea0c	jens-daniel-mueller	2021-01-06
fa85b93	jens-daniel-mueller	2021-01-06
e5cb81a	Donghe-Zhu	2021-01-05
a499f10	Donghe-Zhu	2021-01-05
fb8a752	Donghe-Zhu	2020-12-23
8fae0b2	Donghe-Zhu	2020-12-21
c8b76b3	jens-daniel-mueller	2020-12-19

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] gt_0.2.2         corrr_0.4.3      broom_0.7.5      kableExtra_1.3.1
 [5] knitr_1.30       olsrr_0.5.3      GGally_2.0.0     lubridate_1.7.9 
 [9] metR_0.9.0       scico_1.2.0      patchwork_1.1.1  collapse_1.5.0  
[13] forcats_0.5.0    stringr_1.4.0    dplyr_1.0.2      purrr_0.3.4     
[17] readr_1.4.0      tidyr_1.1.2      tibble_3.0.4     ggplot2_3.3.3   
[21] 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         R6_2.5.0                 nortest_1.0-4           
[10] DBI_1.1.0                colorspace_2.0-0         withr_2.3.0             
[13] gridExtra_2.3            tidyselect_1.1.0         curl_4.3                
[16] compiler_4.0.3           git2r_0.27.1             cli_2.2.0               
[19] rvest_0.3.6              xml2_1.3.2               sass_0.2.0              
[22] labeling_0.4.2           scales_1.1.1             checkmate_2.0.0         
[25] goftest_1.2-2            digest_0.6.27            foreign_0.8-80          
[28] rmarkdown_2.5            rio_0.5.16               pkgconfig_2.0.3         
[31] htmltools_0.5.0          highr_0.8                dbplyr_1.4.4            
[34] rlang_0.4.10             readxl_1.3.1             rstudioapi_0.13         
[37] farver_2.0.3             generics_0.1.0           jsonlite_1.7.2          
[40] zip_2.1.1                car_3.0-10               magrittr_2.0.1          
[43] Matrix_1.2-18            Rcpp_1.0.5               munsell_0.5.0           
[46] fansi_0.4.1              abind_1.4-5              lifecycle_0.2.0         
[49] stringi_1.5.3            whisker_0.4              yaml_2.2.1              
[52] carData_3.0-4            plyr_1.8.6               grid_4.0.3              
[55] blob_1.2.1               parallel_4.0.3           promises_1.1.1          
[58] crayon_1.3.4             lattice_0.20-41          haven_2.3.1             
[61] hms_0.5.3                pillar_1.4.7             reprex_0.3.0            
[64] glue_1.4.2               evaluate_0.14            RcppArmadillo_0.10.1.2.2
[67] data.table_1.13.6        modelr_0.1.8             vctrs_0.3.6             
[70] httpuv_1.5.4             cellranger_1.1.0         gtable_0.3.0            
[73] reshape_0.8.8            assertthat_0.2.1         xfun_0.20               
[76] openxlsx_4.2.3           RcppEigen_0.3.3.9.1      later_1.1.0.1           
[79] viridisLite_0.3.0        ellipsis_0.3.1           here_1.0.1