eMLR - model fitting

workflowr

• Summary
• Checks
• Past versions

Last updated: 2021-03-16

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: 323f626

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 323f626. 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	aecbf75	Donghe-Zhu	2021-03-14	Build site.
html	1b2a0c2	Donghe-Zhu	2021-03-14	Build site.
html	6733e48	Donghe-Zhu	2021-03-12	Build site.
html	ba71e6a	Donghe-Zhu	2021-03-12	Build site.
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 + silicate + phosphate	66	0.6591713	146.28600	2.0481768	1982-1997 –> 1998-2012	2.3666301	308.78580
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	66	0.6600730	146.46645	2.0451681	1982-1997 –> 1998-2012	2.3506754	308.21223
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	66	0.7572964	164.60381	2.3032765	1982-1997 –> 1998-2012	2.4358319	325.80518
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	66	0.7572964	164.60381	2.3032765	1982-1997 –> 1998-2012	2.4358319	325.80518
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	66	0.6574801	145.94691	2.0439713	1982-1997 –> 1998-2012	2.3420642	307.42166
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	66	0.6574801	145.94691	2.0439713	1982-1997 –> 1998-2012	2.3420642	307.42166
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	66	0.6574801	145.94691	2.0439713	1982-1997 –> 1998-2012	2.3420642	307.42166
Atlantic	(-Inf,26]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	66	0.6441997	143.25335	1.9226005	1982-1997 –> 1998-2012	2.0560882	291.30722
Atlantic	(-Inf,26]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	66	0.6441997	143.25335	1.9226005	1982-1997 –> 1998-2012	2.0560882	291.30722
Atlantic	(-Inf,26]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	66	0.6441997	143.25335	1.9226005	1982-1997 –> 1998-2012	2.0560882	291.30722
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	25	0.0440426	-71.18296	0.0902802	1998-2012 –> 2013-2019	0.7135750	77.16173
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	25	0.0439239	-71.31794	0.0921668	1998-2012 –> 2013-2019	0.7030952	74.96806
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	25	0.0438963	-71.34931	0.0930618	1998-2012 –> 2013-2019	0.7039693	75.11715
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + aou + silicate	25	0.0481951	-66.67802	0.1236976	1998-2012 –> 2013-2019	0.6913400	76.35903
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	25	0.0439039	-71.34072	0.0953370	1998-2012 –> 2013-2019	0.7013840	74.60619
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	25	0.0439039	-71.34072	0.0953370	1998-2012 –> 2013-2019	0.7013840	74.60619
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	25	0.0439039	-71.34072	0.0953370	1998-2012 –> 2013-2019	0.7013840	74.60619
Atlantic	(-Inf,26]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	25	0.0438948	-71.35101	0.0947611	1998-2012 –> 2013-2019	0.6880945	71.90234
Atlantic	(-Inf,26]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	25	0.0438948	-71.35101	0.0947611	1998-2012 –> 2013-2019	0.6880945	71.90234
Atlantic	(-Inf,26]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	25	0.0438948	-71.35101	0.0947611	1998-2012 –> 2013-2019	0.6880945	71.90234
Atlantic	(26,26.5]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	425	3.8481871	2365.55958	12.0230141	1982-1997 –> 1998-2012	7.0158532	3773.54759
Atlantic	(26,26.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	425	3.8481871	2365.55958	12.0230141	1982-1997 –> 1998-2012	7.0158532	3773.54759
Atlantic	(26,26.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	425	3.8481871	2365.55958	12.0230141	1982-1997 –> 1998-2012	7.0158532	3773.54759
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	425	3.9226946	2381.85975	13.1801335	1982-1997 –> 1998-2012	7.1790485	3804.81405
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	425	3.9375650	2385.07589	13.4229978	1982-1997 –> 1998-2012	7.2176660	3811.96841
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	425	3.9451972	2386.72184	13.4793446	1982-1997 –> 1998-2012	7.2340750	3815.06269
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	425	3.9451972	2386.72184	13.4793446	1982-1997 –> 1998-2012	7.2340750	3815.06269
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	425	3.9451972	2386.72184	13.4793446	1982-1997 –> 1998-2012	7.2340750	3815.06269
Atlantic	(26,26.5]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	425	4.3078640	2461.47361	14.5624304	1982-1997 –> 1998-2012	7.8173433	3925.00178
Atlantic	(26,26.5]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	425	4.3078640	2461.47361	14.5624304	1982-1997 –> 1998-2012	7.8173433	3925.00178
Atlantic	(26,26.5]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	199	2.8025712	988.89137	7.6989339	1998-2012 –> 2013-2019	6.6507582	3354.45095
Atlantic	(26,26.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	199	2.8025712	988.89137	7.6989339	1998-2012 –> 2013-2019	6.6507582	3354.45095
Atlantic	(26,26.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	199	2.8025712	988.89137	7.6989339	1998-2012 –> 2013-2019	6.6507582	3354.45095
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	199	2.7313224	978.64232	7.5548120	1998-2012 –> 2013-2019	6.6540170	3360.50207
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	199	2.7280026	978.15828	7.5408639	1998-2012 –> 2013-2019	6.6655676	3363.23417
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	199	2.7288867	978.28723	7.5449068	1998-2012 –> 2013-2019	6.6740838	3365.00907
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	199	2.7288867	978.28723	7.5449068	1998-2012 –> 2013-2019	6.6740838	3365.00907
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	199	2.7288867	978.28723	7.5449068	1998-2012 –> 2013-2019	6.6740838	3365.00907
Atlantic	(26,26.5]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	199	2.3946444	926.28496	7.2713354	1998-2012 –> 2013-2019	6.7025083	3387.75857
Atlantic	(26,26.5]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	199	2.3946444	926.28496	7.2713354	1998-2012 –> 2013-2019	6.7025083	3387.75857
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	521	3.1403571	2684.93259	16.8261155	1982-1997 –> 1998-2012	5.6573978	4310.24404
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	521	3.1403571	2684.93259	16.8261155	1982-1997 –> 1998-2012	5.6573978	4310.24404
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	521	3.1403571	2684.93259	16.8261155	1982-1997 –> 1998-2012	5.6573978	4310.24404
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	521	3.1508236	2688.39970	16.9873029	1982-1997 –> 1998-2012	5.7213427	4328.17563
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	521	3.1677881	2693.99495	17.2348725	1982-1997 –> 1998-2012	5.7727496	4342.92816
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	521	3.1741349	2696.08057	17.3161280	1982-1997 –> 1998-2012	5.7903227	4347.97239
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	521	3.1741349	2696.08057	17.3161280	1982-1997 –> 1998-2012	5.7903227	4347.97239
Atlantic	(26.5,26.75]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	521	3.1199039	2678.12384	14.1886801	1982-1997 –> 1998-2012	5.5658726	4283.72912
Atlantic	(26.5,26.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	521	3.1199039	2678.12384	14.1886801	1982-1997 –> 1998-2012	5.5658726	4283.72912
Atlantic	(26.5,26.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	521	3.1199039	2678.12384	14.1886801	1982-1997 –> 1998-2012	5.5658726	4283.72912
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	211	2.4752393	995.26630	13.5241049	1998-2012 –> 2013-2019	5.6155963	3680.19889
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	211	2.4752393	995.26630	13.5241049	1998-2012 –> 2013-2019	5.6155963	3680.19889
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	211	2.4752393	995.26630	13.5241049	1998-2012 –> 2013-2019	5.6155963	3680.19889
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	211	2.4729511	994.87602	13.5832429	1998-2012 –> 2013-2019	5.6237747	3683.27572
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	211	2.4726620	994.82668	13.6432683	1998-2012 –> 2013-2019	5.6404501	3688.82162
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	211	2.4726502	994.82466	13.6374758	1998-2012 –> 2013-2019	5.6467851	3690.90523
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	211	2.4726502	994.82466	13.6374758	1998-2012 –> 2013-2019	5.6467851	3690.90523
Atlantic	(26.5,26.75]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	211	2.4733229	994.93946	13.4966740	1998-2012 –> 2013-2019	5.5932268	3673.06331
Atlantic	(26.5,26.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	211	2.4733229	994.93946	13.4966740	1998-2012 –> 2013-2019	5.5932268	3673.06331
Atlantic	(26.5,26.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	211	2.4733229	994.93946	13.4966740	1998-2012 –> 2013-2019	5.5932268	3673.06331
Atlantic	(26.75,27]	4	cstar_tref ~ aou + silicate + phosphate + phosphate_star	855	2.5728463	4054.35676	11.9116385	1982-1997 –> 1998-2012	4.1520158	6805.07968
Atlantic	(26.75,27]	4	cstar_tref ~ oxygen + aou + silicate + phosphate	855	2.5728463	4054.35676	11.9116385	1982-1997 –> 1998-2012	4.1520158	6805.07968
Atlantic	(26.75,27]	5	cstar_tref ~ oxygen + aou + silicate + phosphate + phosphate_star	855	2.5728463	4054.35676	11.9116385	1982-1997 –> 1998-2012	4.1520158	6805.07968
Atlantic	(26.75,27]	4	cstar_tref ~ oxygen + aou + silicate + phosphate_star	855	2.5728463	4054.35676	11.9116385	1982-1997 –> 1998-2012	4.1520158	6805.07968
Atlantic	(26.75,27]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	855	2.5223861	4022.48601	11.5331125	1982-1997 –> 1998-2012	4.0903615	6764.82269
Atlantic	(26.75,27]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	855	2.5223861	4022.48601	11.5331125	1982-1997 –> 1998-2012	4.0903615	6764.82269
Atlantic	(26.75,27]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	855	2.5223861	4022.48601	11.5331125	1982-1997 –> 1998-2012	4.0903615	6764.82269
Atlantic	(26.75,27]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	855	2.4117043	3945.75545	11.0299838	1982-1997 –> 1998-2012	3.8323448	6544.02315
Atlantic	(26.75,27]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	855	2.4117043	3945.75545	11.0299838	1982-1997 –> 1998-2012	3.8323448	6544.02315
Atlantic	(26.75,27]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	855	2.4117043	3945.75545	11.0299838	1982-1997 –> 1998-2012	3.8323448	6544.02315
Atlantic	(26.75,27]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	425	2.7252227	2072.26539	9.0683468	1998-2012 –> 2013-2019	5.2476089	6094.75139
Atlantic	(26.75,27]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	425	2.7252227	2072.26539	9.0683468	1998-2012 –> 2013-2019	5.2476089	6094.75139
Atlantic	(26.75,27]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	425	2.7252227	2072.26539	9.0683468	1998-2012 –> 2013-2019	5.2476089	6094.75139
Atlantic	(26.75,27]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	425	2.7220465	2071.27414	9.2147292	1998-2012 –> 2013-2019	5.3432476	6159.47089
Atlantic	(26.75,27]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	425	2.7247254	2072.11027	9.3935742	1998-2012 –> 2013-2019	5.4059066	6198.99536
Atlantic	(26.75,27]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	425	2.7272732	2072.90471	9.4459729	1998-2012 –> 2013-2019	5.4336237	6215.76737
Atlantic	(26.75,27]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	425	2.7272732	2072.90471	9.4459729	1998-2012 –> 2013-2019	5.4336237	6215.76737
Atlantic	(26.75,27]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	425	2.9272208	2133.04316	9.5822294	1998-2012 –> 2013-2019	5.3389250	6078.79861
Atlantic	(26.75,27]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	425	2.9272208	2133.04316	9.5822294	1998-2012 –> 2013-2019	5.3389250	6078.79861
Atlantic	(26.75,27]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	425	2.9272208	2133.04316	9.5822294	1998-2012 –> 2013-2019	5.3389250	6078.79861
Atlantic	(27,27.25]	4	cstar_tref ~ sal + aou + silicate + phosphate	1022	2.3323993	4643.36879	10.4766556	1982-1997 –> 1998-2012	3.8388673	8093.33683
Atlantic	(27,27.25]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1022	2.3230449	4637.15457	10.8441024	1982-1997 –> 1998-2012	3.6799366	7892.52858
Atlantic	(27,27.25]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1022	2.3230449	4637.15457	10.8441024	1982-1997 –> 1998-2012	3.6799366	7892.52858
Atlantic	(27,27.25]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1022	2.3230449	4637.15457	10.8441024	1982-1997 –> 1998-2012	3.6799366	7892.52858
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1022	2.3311133	4644.24151	10.6142895	1982-1997 –> 1998-2012	3.7252040	7950.46161
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1022	2.3681648	4676.47409	10.6323907	1982-1997 –> 1998-2012	3.8111195	8047.46109
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	1022	2.3768676	4683.97185	10.6325853	1982-1997 –> 1998-2012	3.8334163	8072.58730
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	1022	2.3768676	4683.97185	10.6325853	1982-1997 –> 1998-2012	3.8334163	8072.58730
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1022	2.3768676	4683.97185	10.6325853	1982-1997 –> 1998-2012	3.8334163	8072.58730
Atlantic	(27,27.25]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	1022	2.5735376	4846.46561	12.8845435	1982-1997 –> 1998-2012	4.1068132	8331.59410
Atlantic	(27,27.25]	4	cstar_tref ~ sal + aou + silicate + phosphate	634	1.9025406	2626.77919	6.7820397	1998-2012 –> 2013-2019	4.2349399	7270.14798
Atlantic	(27,27.25]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	634	1.8765819	2611.35920	6.9669260	1998-2012 –> 2013-2019	4.1996268	7248.51377
Atlantic	(27,27.25]	4	cstar_tref ~ sal + aou + silicate + phosphate_star	634	1.8956251	2622.16178	7.3094081	1998-2012 –> 2013-2019	4.3682851	7384.89408
Atlantic	(27,27.25]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	634	1.8765819	2611.35920	6.9669260	1998-2012 –> 2013-2019	4.1996268	7248.51377
Atlantic	(27,27.25]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	634	1.8765819	2611.35920	6.9669260	1998-2012 –> 2013-2019	4.1996268	7248.51377
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	634	1.8681960	2605.68014	6.9934579	1998-2012 –> 2013-2019	4.1993093	7249.92165
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	634	1.8829582	2615.66029	7.0527384	1998-2012 –> 2013-2019	4.2511230	7292.13438
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	634	1.8854575	2617.34222	7.0582726	1998-2012 –> 2013-2019	4.2623251	7301.31407
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	634	1.8854575	2617.34222	7.0582726	1998-2012 –> 2013-2019	4.2623251	7301.31407
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	634	1.8854575	2617.34222	7.0582726	1998-2012 –> 2013-2019	4.2623251	7301.31407
Atlantic	(27.25,27.5]	4	cstar_tref ~ sal + aou + silicate + phosphate	2082	3.4731827	11104.93731	19.5274910	1982-1997 –> 1998-2012	4.9671179	15762.47618
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2082	3.4730787	11106.81262	19.5923495	1982-1997 –> 1998-2012	4.7730048	15411.35229
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	2082	3.4730787	11106.81262	19.5923495	1982-1997 –> 1998-2012	4.7730048	15411.35229
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	2082	3.4730787	11106.81262	19.5923495	1982-1997 –> 1998-2012	4.7730048	15411.35229
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2082	3.5071477	11147.46006	18.9318630	1982-1997 –> 1998-2012	5.2438538	16191.27042
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	2082	3.4603008	11091.46449	18.9670851	1982-1997 –> 1998-2012	4.6912436	15256.85110
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	2082	3.3757537	10988.45998	20.2192665	1982-1997 –> 1998-2012	4.6620669	15266.13396
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	2082	3.3762152	10989.02917	20.1706040	1982-1997 –> 1998-2012	4.6792476	15299.65985
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	2082	3.3762152	10989.02917	20.1706040	1982-1997 –> 1998-2012	4.6792476	15299.65985
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	2082	3.3762152	10989.02917	20.1706040	1982-1997 –> 1998-2012	4.6792476	15299.65985
Atlantic	(27.25,27.5]	4	cstar_tref ~ sal + aou + silicate + phosphate	852	3.4508931	4540.50201	9.5952246	1998-2012 –> 2013-2019	6.9240758	15645.43932
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	852	3.3766522	4505.44288	9.8885682	1998-2012 –> 2013-2019	6.8497310	15612.25550
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	852	3.3766522	4505.44288	9.8885682	1998-2012 –> 2013-2019	6.8497310	15612.25550
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	852	3.3766522	4505.44288	9.8885682	1998-2012 –> 2013-2019	6.8497310	15612.25550
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	852	3.4454043	4539.78957	9.7298098	1998-2012 –> 2013-2019	6.9057051	15631.25405
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	852	3.3079282	4470.40406	9.2707326	1998-2012 –> 2013-2019	6.6836819	15458.86405
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	852	3.3080498	4470.46673	9.2842314	1998-2012 –> 2013-2019	6.6842650	15459.49591
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	852	3.3080498	4470.46673	9.2842314	1998-2012 –> 2013-2019	6.6842650	15459.49591
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	852	3.3080498	4470.46673	9.2842314	1998-2012 –> 2013-2019	6.6842650	15459.49591
Atlantic	(27.25,27.5]	4	cstar_tref ~ sal + temp + silicate + phosphate_star	852	3.3080558	4468.46982	9.2819015	1998-2012 –> 2013-2019	6.7278557	15508.90939
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	3707	2.4580826	17202.02548	19.3309233	1982-1997 –> 1998-2012	4.1013383	27230.87652
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	3707	2.4580826	17202.02548	19.3309233	1982-1997 –> 1998-2012	4.1013383	27230.87652
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	3707	2.4580826	17202.02548	19.3309233	1982-1997 –> 1998-2012	4.1013383	27230.87652
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	3707	2.5042001	17339.83500	20.3137604	1982-1997 –> 1998-2012	4.1663862	27428.56862
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	3707	2.4472624	17169.31788	19.5434003	1982-1997 –> 1998-2012	3.9740426	26813.81477
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	3707	2.3940793	17006.42276	21.7852238	1982-1997 –> 1998-2012	3.8185595	26288.33655
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	3707	2.3962181	17013.04328	21.6346402	1982-1997 –> 1998-2012	3.8199746	26292.30000
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	3707	2.3962181	17013.04328	21.6346402	1982-1997 –> 1998-2012	3.8199746	26292.30000
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	3707	2.3962181	17013.04328	21.6346402	1982-1997 –> 1998-2012	3.8199746	26292.30000
Atlantic	(27.5,27.75]	4	cstar_tref ~ sal + temp + silicate + phosphate_star	3707	2.5209257	17387.18868	24.5685526	1982-1997 –> 1998-2012	3.9527386	26693.94504
Atlantic	(27.5,27.75]	4	cstar_tref ~ sal + aou + silicate + phosphate	1724	2.2252473	7662.44510	8.8019209	1998-2012 –> 2013-2019	4.6864231	24871.79444
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1724	2.1714617	7580.08113	8.8496375	1998-2012 –> 2013-2019	4.6295443	24782.10661
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1724	2.1714617	7580.08113	8.8496375	1998-2012 –> 2013-2019	4.6295443	24782.10661
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1724	2.1714617	7580.08113	8.8496375	1998-2012 –> 2013-2019	4.6295443	24782.10661
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1724	2.1837456	7599.53137	9.2746814	1998-2012 –> 2013-2019	4.6310080	24768.84925
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1724	2.1696869	7577.26187	8.5365707	1998-2012 –> 2013-2019	4.5637662	24583.68463
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	1724	2.1686654	7575.63808	8.4699729	1998-2012 –> 2013-2019	4.5648835	24588.68137
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	1724	2.1686654	7575.63808	8.4699729	1998-2012 –> 2013-2019	4.5648835	24588.68137
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1724	2.1686654	7575.63808	8.4699729	1998-2012 –> 2013-2019	4.5648835	24588.68137
Atlantic	(27.5,27.75]	4	cstar_tref ~ sal + temp + silicate + phosphate_star	1724	2.2734135	7736.28194	9.9683769	1998-2012 –> 2013-2019	4.7943392	25123.47062
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1493	1.2038819	4805.00659	6.7593452	1982-1997 –> 1998-2012	1.8399229	6647.52284
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1493	1.2038819	4805.00659	6.7593452	1982-1997 –> 1998-2012	1.8399229	6647.52284
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1493	1.2038819	4805.00659	6.7593452	1982-1997 –> 1998-2012	1.8399229	6647.52284
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	1493	1.1524640	4674.67061	6.1587580	1982-1997 –> 1998-2012	1.8716249	6749.56626
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1493	1.2387270	4890.20595	6.8720197	1982-1997 –> 1998-2012	1.8981353	6800.98561
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	1493	1.1524640	4674.67061	6.1587580	1982-1997 –> 1998-2012	1.8716249	6749.56626
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	1493	1.1524640	4674.67061	6.1587580	1982-1997 –> 1998-2012	1.8716249	6749.56626
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	1493	1.2990503	5032.18804	7.1731452	1982-1997 –> 1998-2012	1.9956401	7046.75077
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	1493	1.2990503	5032.18804	7.1731452	1982-1997 –> 1998-2012	1.9956401	7046.75077
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1493	1.2990503	5032.18804	7.1731452	1982-1997 –> 1998-2012	1.9956401	7046.75077
Atlantic	(27.75,27.85]	4	cstar_tref ~ sal + aou + silicate + phosphate	602	1.1923439	1932.21089	5.1084443	1998-2012 –> 2013-2019	2.4743820	6923.03623
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	602	1.1324026	1872.10925	3.7968077	1998-2012 –> 2013-2019	2.3362845	6677.11584
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	602	1.1324026	1872.10925	3.7968077	1998-2012 –> 2013-2019	2.3362845	6677.11584
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	602	1.1324026	1872.10925	3.7968077	1998-2012 –> 2013-2019	2.3362845	6677.11584
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	602	1.0172544	1742.99909	5.1649960	1998-2012 –> 2013-2019	2.1697184	6417.66970
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	602	1.1542981	1895.16688	4.2322700	1998-2012 –> 2013-2019	2.3930251	6785.37283
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	602	1.0172544	1742.99909	5.1649960	1998-2012 –> 2013-2019	2.1697184	6417.66970
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	602	1.0172544	1742.99909	5.1649960	1998-2012 –> 2013-2019	2.1697184	6417.66970
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	602	1.1892935	1931.12678	4.5202548	1998-2012 –> 2013-2019	2.4883439	6963.31482
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	602	1.1892935	1931.12678	4.5202548	1998-2012 –> 2013-2019	2.4883439	6963.31482
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	1432	1.2769885	4778.10113	8.7799676	1982-1997 –> 1998-2012	2.2286642	7519.96094
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	1432	1.2769885	4778.10113	8.7799676	1982-1997 –> 1998-2012	2.2286642	7519.96094
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	1432	1.2769885	4778.10113	8.7799676	1982-1997 –> 1998-2012	2.2286642	7519.96094
Atlantic	(27.85,27.95]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	1432	1.5733116	5373.75519	8.5379549	1982-1997 –> 1998-2012	2.6351889	8331.87718
Atlantic	(27.85,27.95]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1432	1.5715734	5372.58921	8.5422972	1982-1997 –> 1998-2012	2.6327735	8331.44050
Atlantic	(27.85,27.95]	4	cstar_tref ~ temp + oxygen + aou + phosphate	1432	1.5733116	5373.75519	8.5379549	1982-1997 –> 1998-2012	2.6351889	8331.87718
Atlantic	(27.85,27.95]	5	cstar_tref ~ temp + oxygen + aou + phosphate + phosphate_star	1432	1.5733116	5373.75519	8.5379549	1982-1997 –> 1998-2012	2.6351889	8331.87718
Atlantic	(27.85,27.95]	4	cstar_tref ~ temp + oxygen + aou + phosphate_star	1432	1.5733116	5373.75519	8.5379549	1982-1997 –> 1998-2012	2.6351889	8331.87718
Atlantic	(27.85,27.95]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	1432	1.5715734	5372.58921	8.5422972	1982-1997 –> 1998-2012	2.6327735	8331.44050
Atlantic	(27.85,27.95]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	1432	1.5715734	5372.58921	8.5422972	1982-1997 –> 1998-2012	2.6327735	8331.44050
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	526	1.1217760	1627.61197	8.7596639	1998-2012 –> 2013-2019	2.6670066	6951.78786
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	526	1.1217760	1627.61197	8.7596639	1998-2012 –> 2013-2019	2.6670066	6951.78786
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	526	1.1217760	1627.61197	8.7596639	1998-2012 –> 2013-2019	2.6670066	6951.78786
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	526	0.8686718	1358.61237	3.9573036	1998-2012 –> 2013-2019	2.1456603	6136.71350
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	526	1.1840833	1684.47865	9.5292198	1998-2012 –> 2013-2019	2.7914549	7121.57379
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	526	0.8686718	1358.61237	3.9573036	1998-2012 –> 2013-2019	2.1456603	6136.71350
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	526	0.8686718	1358.61237	3.9573036	1998-2012 –> 2013-2019	2.1456603	6136.71350
Atlantic	(27.85,27.95]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	526	1.2607414	1750.47170	7.4103180	1998-2012 –> 2013-2019	2.8323148	7123.06091
Atlantic	(27.85,27.95]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	526	1.2607414	1750.47170	7.4103180	1998-2012 –> 2013-2019	2.8323148	7123.06091
Atlantic	(27.85,27.95]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	526	1.2607414	1750.47170	7.4103180	1998-2012 –> 2013-2019	2.8323148	7123.06091
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1742	4.5193193	10212.71290	18.1638465	1982-1997 –> 1998-2012	8.1656819	16959.55652
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1742	4.5193193	10212.71290	18.1638465	1982-1997 –> 1998-2012	8.1656819	16959.55652
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1742	4.5193193	10212.71290	18.1638465	1982-1997 –> 1998-2012	8.1656819	16959.55652
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1742	4.5254696	10217.45101	18.1089389	1982-1997 –> 1998-2012	8.1614693	16957.23082
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1742	4.5312636	10221.90875	18.0789190	1982-1997 –> 1998-2012	8.1651862	16960.27024
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	1742	4.5367376	10226.11511	18.3917284	1982-1997 –> 1998-2012	8.3231307	17066.48978
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	1742	4.5367376	10226.11511	18.3917284	1982-1997 –> 1998-2012	8.3231307	17066.48978
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	1742	4.5314418	10222.04577	18.0837165	1982-1997 –> 1998-2012	8.1648671	16960.06761
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	1742	4.5314418	10222.04577	18.0837165	1982-1997 –> 1998-2012	8.1648671	16960.06761
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1742	4.5314418	10222.04577	18.0837165	1982-1997 –> 1998-2012	8.1648671	16960.06761
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	652	0.8773869	1693.72328	4.9823040	1998-2012 –> 2013-2019	5.3967062	11906.43618
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	652	0.8773869	1693.72328	4.9823040	1998-2012 –> 2013-2019	5.3967062	11906.43618
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	652	0.8773869	1693.72328	4.9823040	1998-2012 –> 2013-2019	5.3967062	11906.43618
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	652	0.6079228	1215.28542	3.0285389	1998-2012 –> 2013-2019	5.1446604	11441.40054
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	652	0.9365914	1778.87315	5.3839435	1998-2012 –> 2013-2019	5.4620609	11996.32416
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	652	0.6079228	1215.28542	3.0285389	1998-2012 –> 2013-2019	5.1446604	11441.40054
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	652	0.6079228	1215.28542	3.0285389	1998-2012 –> 2013-2019	5.1446604	11441.40054
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	652	0.9819451	1840.53710	5.5789771	1998-2012 –> 2013-2019	5.5133869	12062.58287
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	652	0.9819451	1840.53710	5.5789771	1998-2012 –> 2013-2019	5.5133869	12062.58287
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	652	0.9819451	1840.53710	5.5789771	1998-2012 –> 2013-2019	5.5133869	12062.58287
Atlantic	(28.05,28.1]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	1022	1.1584146	3214.88539	8.4362920	1982-1997 –> 1998-2012	1.8713099	4756.73761
Atlantic	(28.05,28.1]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	1022	1.1584146	3214.88539	8.4362920	1982-1997 –> 1998-2012	1.8713099	4756.73761
Atlantic	(28.05,28.1]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	1022	1.1584146	3214.88539	8.4362920	1982-1997 –> 1998-2012	1.8713099	4756.73761
Atlantic	(28.05,28.1]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	1022	1.1596697	3215.09869	8.4449733	1982-1997 –> 1998-2012	1.8856301	4780.63038
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1022	1.0483946	3010.90983	7.8529237	1982-1997 –> 1998-2012	1.7164427	4460.88821
Atlantic	(28.05,28.1]	4	cstar_tref ~ temp + oxygen + aou + phosphate	1022	1.1596697	3215.09869	8.4449733	1982-1997 –> 1998-2012	1.8856301	4780.63038
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + phosphate + phosphate_star	1022	1.1596697	3215.09869	8.4449733	1982-1997 –> 1998-2012	1.8856301	4780.63038
Atlantic	(28.05,28.1]	4	cstar_tref ~ temp + oxygen + aou + phosphate_star	1022	1.1596697	3215.09869	8.4449733	1982-1997 –> 1998-2012	1.8856301	4780.63038
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	1022	1.0483946	3010.90983	7.8529237	1982-1997 –> 1998-2012	1.7164427	4460.88821
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	1022	1.0483946	3010.90983	7.8529237	1982-1997 –> 1998-2012	1.7164427	4460.88821
Atlantic	(28.05,28.1]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	384	0.5158719	595.40804	4.4019146	1998-2012 –> 2013-2019	1.6742865	3810.29343
Atlantic	(28.05,28.1]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	384	0.5158719	595.40804	4.4019146	1998-2012 –> 2013-2019	1.6742865	3810.29343
Atlantic	(28.05,28.1]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	384	0.5158719	595.40804	4.4019146	1998-2012 –> 2013-2019	1.6742865	3810.29343
Atlantic	(28.05,28.1]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	384	0.5706321	670.88866	3.5681959	1998-2012 –> 2013-2019	1.7303018	3885.98735
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	384	0.4418284	476.41658	3.6655050	1998-2012 –> 2013-2019	1.4902230	3487.32641
Atlantic	(28.05,28.1]	4	cstar_tref ~ temp + oxygen + aou + phosphate	384	0.5706321	670.88866	3.5681959	1998-2012 –> 2013-2019	1.7303018	3885.98735
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + phosphate + phosphate_star	384	0.5706321	670.88866	3.5681959	1998-2012 –> 2013-2019	1.7303018	3885.98735
Atlantic	(28.05,28.1]	4	cstar_tref ~ temp + oxygen + aou + phosphate_star	384	0.5706321	670.88866	3.5681959	1998-2012 –> 2013-2019	1.7303018	3885.98735
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	384	0.4418284	476.41658	3.6655050	1998-2012 –> 2013-2019	1.4902230	3487.32641
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	384	0.4418284	476.41658	3.6655050	1998-2012 –> 2013-2019	1.4902230	3487.32641
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1466	1.0636370	4355.21507	6.1657158	1982-1997 –> 1998-2012	1.6894643	5931.43863
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1466	1.0636370	4355.21507	6.1657158	1982-1997 –> 1998-2012	1.6894643	5931.43863
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1466	1.0636370	4355.21507	6.1657158	1982-1997 –> 1998-2012	1.6894643	5931.43863
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	1466	1.0199237	4232.16982	6.3788982	1982-1997 –> 1998-2012	1.6919482	5925.47973
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	1466	1.0199237	4232.16982	6.3788982	1982-1997 –> 1998-2012	1.6919482	5925.47973
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	1466	1.0199237	4232.16982	6.3788982	1982-1997 –> 1998-2012	1.6919482	5925.47973
Atlantic	(28.1,28.15]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	1466	1.0290196	4256.20200	6.4204600	1982-1997 –> 1998-2012	1.7079187	5964.24424
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1466	0.8822868	3807.12962	7.7969805	1982-1997 –> 1998-2012	1.4012070	5075.39221
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	1466	0.8822868	3807.12962	7.7969805	1982-1997 –> 1998-2012	1.4012070	5075.39221
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	1466	0.8822868	3807.12962	7.7969805	1982-1997 –> 1998-2012	1.4012070	5075.39221
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	582	0.5162736	896.10271	4.4951585	1998-2012 –> 2013-2019	1.5361974	5128.27253
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	582	0.5162736	896.10271	4.4951585	1998-2012 –> 2013-2019	1.5361974	5128.27253
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	582	0.5162736	896.10271	4.4951585	1998-2012 –> 2013-2019	1.5361974	5128.27253
Atlantic	(28.1,28.15]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	582	0.5232831	909.79993	4.3024426	1998-2012 –> 2013-2019	1.5523026	5166.00193
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	582	0.3897878	568.97846	4.2817939	1998-2012 –> 2013-2019	1.2720745	4376.10808
Atlantic	(28.1,28.15]	4	cstar_tref ~ temp + oxygen + aou + phosphate	582	0.5232831	909.79993	4.3024426	1998-2012 –> 2013-2019	1.5523026	5166.00193
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + oxygen + aou + phosphate + phosphate_star	582	0.5232831	909.79993	4.3024426	1998-2012 –> 2013-2019	1.5523026	5166.00193
Atlantic	(28.1,28.15]	4	cstar_tref ~ temp + oxygen + aou + phosphate_star	582	0.5232831	909.79993	4.3024426	1998-2012 –> 2013-2019	1.5523026	5166.00193
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	582	0.3897878	568.97846	4.2817939	1998-2012 –> 2013-2019	1.2720745	4376.10808
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	582	0.3897878	568.97846	4.2817939	1998-2012 –> 2013-2019	1.2720745	4376.10808
Atlantic	(28.15,28.2]	4	cstar_tref ~ aou + silicate + phosphate + phosphate_star	2690	0.9002057	7080.27927	9.9903677	1982-1997 –> 1998-2012	1.7745126	10761.13891
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2690	0.8864214	6999.26154	11.1469166	1982-1997 –> 1998-2012	1.7448531	10629.78680
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	2690	0.8864214	6999.26154	11.1469166	1982-1997 –> 1998-2012	1.7448531	10629.78680
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	2690	0.8864214	6999.26154	11.1469166	1982-1997 –> 1998-2012	1.7448531	10629.78680
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2690	0.7128535	5826.87019	21.9287876	1982-1997 –> 1998-2012	1.5534334	9397.37642
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	2690	0.7128535	5826.87019	21.9287876	1982-1997 –> 1998-2012	1.5534334	9397.37642
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	2690	0.7128535	5826.87019	21.9287876	1982-1997 –> 1998-2012	1.5534334	9397.37642
Atlantic	(28.15,28.2]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	2690	0.7562818	6145.03358	12.7582986	1982-1997 –> 1998-2012	1.6025400	9734.76772
Atlantic	(28.15,28.2]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	2690	0.7562818	6145.03358	12.7582986	1982-1997 –> 1998-2012	1.6025400	9734.76772
Atlantic	(28.15,28.2]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	2690	0.7562818	6145.03358	12.7582986	1982-1997 –> 1998-2012	1.6025400	9734.76772
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1293	0.6856262	2707.35981	2.6406083	1998-2012 –> 2013-2019	1.5720476	9706.62135
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1293	0.6856262	2707.35981	2.6406083	1998-2012 –> 2013-2019	1.5720476	9706.62135
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1293	0.6856262	2707.35981	2.6406083	1998-2012 –> 2013-2019	1.5720476	9706.62135
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	1293	0.5462649	2119.74708	1.7217459	1998-2012 –> 2013-2019	1.2591184	7946.61728
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1293	0.7095935	2796.21414	2.7904553	1998-2012 –> 2013-2019	1.6178907	9926.63561
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	1293	0.5462649	2119.74708	1.7217459	1998-2012 –> 2013-2019	1.2591184	7946.61728
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	1293	0.5462649	2119.74708	1.7217459	1998-2012 –> 2013-2019	1.2591184	7946.61728
Atlantic	(28.15,28.2]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1293	0.6469827	2557.33823	2.3542884	1998-2012 –> 2013-2019	1.4032645	8702.37180
Atlantic	(28.15,28.2]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	1293	0.6469827	2557.33823	2.3542884	1998-2012 –> 2013-2019	1.4032645	8702.37180
Atlantic	(28.15,28.2]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	1293	0.6469827	2557.33823	2.3542884	1998-2012 –> 2013-2019	1.4032645	8702.37180
Atlantic	(28.2, Inf]	4	cstar_tref ~ sal + aou + silicate + phosphate	5021	0.8620888	12770.77857	3.4164614	1982-1997 –> 1998-2012	1.7133239	18918.67618
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	5021	0.8453747	12576.17212	4.0966513	1982-1997 –> 1998-2012	1.6963174	18724.39445
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	5021	0.8453747	12576.17212	4.0966513	1982-1997 –> 1998-2012	1.6963174	18724.39445
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	5021	0.8453747	12576.17212	4.0966513	1982-1997 –> 1998-2012	1.6963174	18724.39445
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	5021	0.8461863	12585.80883	4.0848831	1982-1997 –> 1998-2012	1.6971298	18734.03507
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	5021	0.8496374	12626.68064	4.1936700	1982-1997 –> 1998-2012	1.7008248	18776.30513
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	5021	0.8480923	12608.40230	4.1405739	1982-1997 –> 1998-2012	1.6991109	18757.05915
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	5021	0.8480923	12608.40230	4.1405739	1982-1997 –> 1998-2012	1.6991109	18757.05915
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	5021	0.8480923	12608.40230	4.1405739	1982-1997 –> 1998-2012	1.6991109	18757.05915
Atlantic	(28.2, Inf]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	5021	0.8398056	12509.79926	5.5423887	1982-1997 –> 1998-2012	1.7187880	18816.16649
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	1688	0.5346185	2690.27904	2.7359339	1998-2012 –> 2013-2019	1.4035026	15541.90173
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	1688	0.5346185	2690.27904	2.7359339	1998-2012 –> 2013-2019	1.4035026	15541.90173
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	1688	0.5346185	2690.27904	2.7359339	1998-2012 –> 2013-2019	1.4035026	15541.90173
Atlantic	(28.2, Inf]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	1688	0.6551946	3374.88610	3.2951802	1998-2012 –> 2013-2019	1.5311631	16306.05366
Atlantic	(28.2, Inf]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1688	0.6183071	3181.25678	2.4866152	1998-2012 –> 2013-2019	1.4581126	15691.05604
Atlantic	(28.2, Inf]	4	cstar_tref ~ temp + oxygen + aou + phosphate	1688	0.6551946	3374.88610	3.2951802	1998-2012 –> 2013-2019	1.5311631	16306.05366
Atlantic	(28.2, Inf]	5	cstar_tref ~ temp + oxygen + aou + phosphate + phosphate_star	1688	0.6551946	3374.88610	3.2951802	1998-2012 –> 2013-2019	1.5311631	16306.05366
Atlantic	(28.2, Inf]	4	cstar_tref ~ temp + oxygen + aou + phosphate_star	1688	0.6551946	3374.88610	3.2951802	1998-2012 –> 2013-2019	1.5311631	16306.05366
Atlantic	(28.2, Inf]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	1688	0.6183071	3181.25678	2.4866152	1998-2012 –> 2013-2019	1.4581126	15691.05604
Atlantic	(28.2, Inf]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	1688	0.6183071	3181.25678	2.4866152	1998-2012 –> 2013-2019	1.4581126	15691.05604
Indo-Pacific	(-Inf,26]	4	cstar_tref ~ sal + aou + phosphate + phosphate_star	2630	5.9123026	16822.82270	24.6255914	1982-1997 –> 1998-2012	11.2200388	23857.17110
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2630	5.6660966	16601.08906	22.3298123	1982-1997 –> 1998-2012	10.8607596	23588.46990
Indo-Pacific	(-Inf,26]	4	cstar_tref ~ sal + oxygen + aou + phosphate	2630	5.9123026	16822.82270	24.6255914	1982-1997 –> 1998-2012	11.2200388	23857.17110
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + oxygen + aou + phosphate + phosphate_star	2630	5.9123026	16822.82270	24.6255914	1982-1997 –> 1998-2012	11.2200388	23857.17110
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	2630	5.6660966	16601.08906	22.3298123	1982-1997 –> 1998-2012	10.8607596	23588.46990
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	2630	5.6660966	16601.08906	22.3298123	1982-1997 –> 1998-2012	10.8607596	23588.46990
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2630	5.8434473	16763.20470	26.1150528	1982-1997 –> 1998-2012	10.8516180	23667.44571
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	2630	6.0450601	16941.62654	20.7139099	1982-1997 –> 1998-2012	10.9788124	23811.82376
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	2630	5.8434473	16763.20470	26.1150528	1982-1997 –> 1998-2012	10.8516180	23667.44571
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	2630	5.8434473	16763.20470	26.1150528	1982-1997 –> 1998-2012	10.8516180	23667.44571
Indo-Pacific	(-Inf,26]	4	cstar_tref ~ sal + aou + phosphate + phosphate_star	2149	5.6899852	13583.56331	23.7304247	1998-2012 –> 2013-2019	11.6022878	30406.38601
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2149	5.3293238	13304.11614	23.1544130	1998-2012 –> 2013-2019	10.9954204	29905.20521
Indo-Pacific	(-Inf,26]	4	cstar_tref ~ sal + oxygen + aou + phosphate	2149	5.6899852	13583.56331	23.7304247	1998-2012 –> 2013-2019	11.6022878	30406.38601
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + oxygen + aou + phosphate + phosphate_star	2149	5.6899852	13583.56331	23.7304247	1998-2012 –> 2013-2019	11.6022878	30406.38601
Indo-Pacific	(-Inf,26]	4	cstar_tref ~ sal + oxygen + aou + phosphate_star	2149	5.6899852	13583.56331	23.7304247	1998-2012 –> 2013-2019	11.6022878	30406.38601
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	2149	5.3293238	13304.11614	23.1544130	1998-2012 –> 2013-2019	10.9954204	29905.20521
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	2149	5.3293238	13304.11614	23.1544130	1998-2012 –> 2013-2019	10.9954204	29905.20521
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2149	5.5249053	13459.02356	24.6947128	1998-2012 –> 2013-2019	11.3683526	30222.22825
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	2149	5.5249053	13459.02356	24.6947128	1998-2012 –> 2013-2019	11.3683526	30222.22825
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	2149	5.5249053	13459.02356	24.6947128	1998-2012 –> 2013-2019	11.3683526	30222.22825
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2648	4.4975277	15491.38191	40.6448953	1982-1997 –> 1998-2012	8.2407309	24269.79877
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	2648	4.4975277	15491.38191	40.6448953	1982-1997 –> 1998-2012	8.2407309	24269.79877
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	2648	4.4975277	15491.38191	40.6448953	1982-1997 –> 1998-2012	8.2407309	24269.79877
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	2648	4.6093689	15621.46809	44.1528270	1982-1997 –> 1998-2012	8.4217951	24458.52193
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	2648	4.7414430	15771.08321	44.0432510	1982-1997 –> 1998-2012	8.6380295	24678.00939
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	2648	5.1430871	16201.71144	37.8878230	1982-1997 –> 1998-2012	9.2026886	25239.78618
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	2648	5.1430871	16201.71144	37.8878230	1982-1997 –> 1998-2012	9.2026886	25239.78618
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	2648	4.7680588	15800.72885	44.0132410	1982-1997 –> 1998-2012	8.6766104	24717.46614
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	2648	4.7680588	15800.72885	44.0132410	1982-1997 –> 1998-2012	8.6766104	24717.46614
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	2648	4.7680588	15800.72885	44.0132410	1982-1997 –> 1998-2012	8.6766104	24717.46614
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1471	3.6943928	8033.17047	20.9215804	1998-2012 –> 2013-2019	8.1919205	23524.55238
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1471	3.6943928	8033.17047	20.9215804	1998-2012 –> 2013-2019	8.1919205	23524.55238
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1471	3.6943928	8033.17047	20.9215804	1998-2012 –> 2013-2019	8.1919205	23524.55238
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1471	3.7518193	8078.54979	20.8916072	1998-2012 –> 2013-2019	8.3611881	23700.01788
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1471	3.7775516	8098.65902	20.7211356	1998-2012 –> 2013-2019	8.5189947	23869.74223
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	1471	3.7899589	8108.30605	20.7125254	1998-2012 –> 2013-2019	8.5580177	23909.03490
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	1471	3.7899589	8108.30605	20.7125254	1998-2012 –> 2013-2019	8.5580177	23909.03490
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1471	3.7899589	8108.30605	20.7125254	1998-2012 –> 2013-2019	8.5580177	23909.03490
Indo-Pacific	(26,26.5]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	1471	4.0665720	8315.55589	22.3051985	1998-2012 –> 2013-2019	9.1013723	24404.57036
Indo-Pacific	(26,26.5]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	1471	4.0665720	8315.55589	22.3051985	1998-2012 –> 2013-2019	9.1013723	24404.57036
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2385	4.5775343	14038.27234	18.8115832	1982-1997 –> 1998-2012	7.7168084	20849.16305
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	2385	4.5775343	14038.27234	18.8115832	1982-1997 –> 1998-2012	7.7168084	20849.16305
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	2385	4.5775343	14038.27234	18.8115832	1982-1997 –> 1998-2012	7.7168084	20849.16305
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2385	4.6354693	14098.26442	15.2743281	1982-1997 –> 1998-2012	7.9967846	21090.39478
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	2385	4.7914538	14256.13438	19.2727838	1982-1997 –> 1998-2012	8.1470105	21243.71751
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	2385	4.6354693	14098.26442	15.2743281	1982-1997 –> 1998-2012	7.9967846	21090.39478
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	2385	4.6354693	14098.26442	15.2743281	1982-1997 –> 1998-2012	7.9967846	21090.39478
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	2385	4.7271062	14191.64089	15.7889852	1982-1997 –> 1998-2012	8.1018364	21194.33424
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	2385	4.7271062	14191.64089	15.7889852	1982-1997 –> 1998-2012	8.1018364	21194.33424
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	2385	4.7271062	14191.64089	15.7889852	1982-1997 –> 1998-2012	8.1018364	21194.33424
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1329	3.2330242	6904.48361	16.5464518	1998-2012 –> 2013-2019	7.8105585	20942.75595
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1329	3.2330242	6904.48361	16.5464518	1998-2012 –> 2013-2019	7.8105585	20942.75595
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1329	3.2330242	6904.48361	16.5464518	1998-2012 –> 2013-2019	7.8105585	20942.75595
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	1329	3.4239406	7056.98403	14.2361648	1998-2012 –> 2013-2019	8.0594099	21155.24845
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1329	3.4057129	7042.79608	18.0529528	1998-2012 –> 2013-2019	8.1971667	21298.93046
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	1329	3.4239406	7056.98403	14.2361648	1998-2012 –> 2013-2019	8.0594099	21155.24845
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	1329	3.4239406	7056.98403	14.2361648	1998-2012 –> 2013-2019	8.0594099	21155.24845
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1329	3.0927721	6786.60086	18.0189689	1998-2012 –> 2013-2019	7.8198783	20978.24175
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	1329	3.0927721	6786.60086	18.0189689	1998-2012 –> 2013-2019	7.8198783	20978.24175
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	1329	3.0927721	6786.60086	18.0189689	1998-2012 –> 2013-2019	7.8198783	20978.24175
Indo-Pacific	(26.75,27]	4	cstar_tref ~ sal + aou + phosphate + phosphate_star	3460	4.3846898	20059.63730	16.0687755	1982-1997 –> 1998-2012	7.6038289	29787.18638
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	3460	4.3823944	20058.01370	16.3220378	1982-1997 –> 1998-2012	7.5541159	29731.85541
Indo-Pacific	(26.75,27]	4	cstar_tref ~ sal + oxygen + aou + phosphate	3460	4.3846898	20059.63730	16.0687755	1982-1997 –> 1998-2012	7.6038289	29787.18638
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + oxygen + aou + phosphate + phosphate_star	3460	4.3846898	20059.63730	16.0687755	1982-1997 –> 1998-2012	7.6038289	29787.18638
Indo-Pacific	(26.75,27]	4	cstar_tref ~ sal + oxygen + aou + phosphate_star	3460	4.3846898	20059.63730	16.0687755	1982-1997 –> 1998-2012	7.6038289	29787.18638
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	3460	4.3823944	20058.01370	16.3220378	1982-1997 –> 1998-2012	7.5541159	29731.85541
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	3460	4.3823944	20058.01370	16.3220378	1982-1997 –> 1998-2012	7.5541159	29731.85541
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	3460	4.2922314	19914.15728	21.8298908	1982-1997 –> 1998-2012	7.4333954	29551.65631
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	3460	4.2922314	19914.15728	21.8298908	1982-1997 –> 1998-2012	7.4333954	29551.65631
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	3460	4.2922314	19914.15728	21.8298908	1982-1997 –> 1998-2012	7.4333954	29551.65631
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1765	3.8037880	9738.92388	12.2952264	1998-2012 –> 2013-2019	8.1861824	29796.93758
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1765	3.8037880	9738.92388	12.2952264	1998-2012 –> 2013-2019	8.1861824	29796.93758
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1765	3.8037880	9738.92388	12.2952264	1998-2012 –> 2013-2019	8.1861824	29796.93758
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	1765	3.8936510	9821.34889	15.8177101	1998-2012 –> 2013-2019	8.1858824	29735.50618
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	1765	3.8936510	9821.34889	15.8177101	1998-2012 –> 2013-2019	8.1858824	29735.50618
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	1765	3.8936510	9821.34889	15.8177101	1998-2012 –> 2013-2019	8.1858824	29735.50618
Indo-Pacific	(26.75,27]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1765	3.7693127	9706.78414	14.7360756	1998-2012 –> 2013-2019	8.1868024	29819.99430
Indo-Pacific	(26.75,27]	4	cstar_tref ~ temp + oxygen + aou + phosphate_star	1765	3.9598723	9878.88056	16.8972842	1998-2012 –> 2013-2019	8.3793405	29993.18936
Indo-Pacific	(26.75,27]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	1765	3.7693127	9706.78414	14.7360756	1998-2012 –> 2013-2019	8.1868024	29819.99430
Indo-Pacific	(26.75,27]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	1765	3.7693127	9706.78414	14.7360756	1998-2012 –> 2013-2019	8.1868024	29819.99430
Indo-Pacific	(27,27.25]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	4602	4.2227644	26332.18001	26.0074533	1982-1997 –> 1998-2012	7.0994868	36852.58460
Indo-Pacific	(27,27.25]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	4602	4.2227644	26332.18001	26.0074533	1982-1997 –> 1998-2012	7.0994868	36852.58460
Indo-Pacific	(27,27.25]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	4602	4.2227644	26332.18001	26.0074533	1982-1997 –> 1998-2012	7.0994868	36852.58460
Indo-Pacific	(27,27.25]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	4602	4.2647374	26421.21330	24.5515466	1982-1997 –> 1998-2012	7.1781771	36993.44379
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	4602	3.6512504	24993.73164	16.5259190	1982-1997 –> 1998-2012	6.5408888	35533.14836
Indo-Pacific	(27,27.25]	4	cstar_tref ~ temp + oxygen + aou + phosphate	4602	4.2647374	26421.21330	24.5515466	1982-1997 –> 1998-2012	7.1781771	36993.44379
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + oxygen + aou + phosphate + phosphate_star	4602	4.2647374	26421.21330	24.5515466	1982-1997 –> 1998-2012	7.1781771	36993.44379
Indo-Pacific	(27,27.25]	4	cstar_tref ~ temp + oxygen + aou + phosphate_star	4602	4.2647374	26421.21330	24.5515466	1982-1997 –> 1998-2012	7.1781771	36993.44379
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	4602	3.6512504	24993.73164	16.5259190	1982-1997 –> 1998-2012	6.5408888	35533.14836
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	4602	3.6512504	24993.73164	16.5259190	1982-1997 –> 1998-2012	6.5408888	35533.14836
Indo-Pacific	(27,27.25]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2578	4.2265824	14761.87306	19.4557796	1998-2012 –> 2013-2019	8.4493467	41094.05307
Indo-Pacific	(27,27.25]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	2578	4.2265824	14761.87306	19.4557796	1998-2012 –> 2013-2019	8.4493467	41094.05307
Indo-Pacific	(27,27.25]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	2578	4.2265824	14761.87306	19.4557796	1998-2012 –> 2013-2019	8.4493467	41094.05307
Indo-Pacific	(27,27.25]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	2578	4.3231966	14876.40573	18.1291512	1998-2012 –> 2013-2019	8.5879340	41297.61902
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	2578	3.9209008	14374.80029	15.9055025	1998-2012 –> 2013-2019	7.5721512	39368.53194
Indo-Pacific	(27,27.25]	4	cstar_tref ~ temp + oxygen + aou + phosphate	2578	4.3231966	14876.40573	18.1291512	1998-2012 –> 2013-2019	8.5879340	41297.61902
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + oxygen + aou + phosphate + phosphate_star	2578	4.3231966	14876.40573	18.1291512	1998-2012 –> 2013-2019	8.5879340	41297.61902
Indo-Pacific	(27,27.25]	4	cstar_tref ~ temp + oxygen + aou + phosphate_star	2578	4.3231966	14876.40573	18.1291512	1998-2012 –> 2013-2019	8.5879340	41297.61902
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	2578	3.9209008	14374.80029	15.9055025	1998-2012 –> 2013-2019	7.5721512	39368.53194
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	2578	3.9209008	14374.80029	15.9055025	1998-2012 –> 2013-2019	7.5721512	39368.53194
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	5832	3.6008356	31508.01850	12.0046897	1982-1997 –> 1998-2012	6.1632969	45138.68057
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	5832	3.5430879	31319.44308	11.5165495	1982-1997 –> 1998-2012	6.0894346	44913.70472
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + aou + silicate	5832	3.0686742	29642.70795	14.0261468	1982-1997 –> 1998-2012	5.9485752	43947.23489
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	5832	3.5621002	31381.86527	11.7330693	1982-1997 –> 1998-2012	6.1386018	45044.05615
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	5832	3.5621002	31381.86527	11.7330693	1982-1997 –> 1998-2012	6.1386018	45044.05615
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	5832	3.5621002	31381.86527	11.7330693	1982-1997 –> 1998-2012	6.1386018	45044.05615
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	5832	3.1116328	29804.86077	18.6001083	1982-1997 –> 1998-2012	5.3490296	42652.79231
Indo-Pacific	(27.25,27.5]	4	cstar_tref ~ temp + oxygen + aou + silicate	5832	3.1889293	30089.06734	17.8089669	1982-1997 –> 1998-2012	6.0700064	44393.95017
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	5832	3.1116328	29804.86077	18.6001083	1982-1997 –> 1998-2012	5.3490296	42652.79231
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	5832	3.1116328	29804.86077	18.6001083	1982-1997 –> 1998-2012	5.3490296	42652.79231
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	3422	3.4387736	18178.34160	13.7873631	1998-2012 –> 2013-2019	6.9818614	49497.78468
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	3422	3.4626480	18225.69346	14.1337349	1998-2012 –> 2013-2019	7.0482836	49684.37110
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + aou + silicate	3422	2.9674206	17169.38673	11.9858738	1998-2012 –> 2013-2019	6.0360948	46812.09468
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	3422	3.4424625	18185.67954	13.9244017	1998-2012 –> 2013-2019	7.0045627	49567.54480
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	3422	3.4424625	18185.67954	13.9244017	1998-2012 –> 2013-2019	7.0045627	49567.54480
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	3422	3.4424625	18185.67954	13.9244017	1998-2012 –> 2013-2019	7.0045627	49567.54480
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	3422	3.3521746	18003.78126	14.9487549	1998-2012 –> 2013-2019	6.4638075	47808.64203
Indo-Pacific	(27.25,27.5]	4	cstar_tref ~ temp + oxygen + aou + silicate	3422	3.3739326	18046.06006	15.7011643	1998-2012 –> 2013-2019	6.5628619	48135.12741
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	3422	3.3521746	18003.78126	14.9487549	1998-2012 –> 2013-2019	6.4638075	47808.64203
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	3422	3.3521746	18003.78126	14.9487549	1998-2012 –> 2013-2019	6.4638075	47808.64203
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	6762	3.0099602	34106.18362	10.7452293	1982-1997 –> 1998-2012	5.3300247	48732.27214
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	6762	2.9740300	33943.77474	10.2070917	1982-1997 –> 1998-2012	5.2852457	48545.16245
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	6762	2.9827371	33983.31141	10.3326712	1982-1997 –> 1998-2012	5.3212343	48660.55302
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + aou + silicate	6762	2.8404458	33322.25358	16.9101113	1982-1997 –> 1998-2012	5.1388209	47887.62840
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	6762	2.9745224	33946.01404	10.2300353	1982-1997 –> 1998-2012	5.3087642	48611.48222
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	6762	2.9745224	33946.01404	10.2300353	1982-1997 –> 1998-2012	5.3087642	48611.48222
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	6762	2.9745224	33946.01404	10.2300353	1982-1997 –> 1998-2012	5.3087642	48611.48222
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	6762	3.0424685	34251.46277	20.7180320	1982-1997 –> 1998-2012	5.1315768	48199.75160
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	6762	3.0424685	34251.46277	20.7180320	1982-1997 –> 1998-2012	5.1315768	48199.75160
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	6762	3.0424685	34251.46277	20.7180320	1982-1997 –> 1998-2012	5.1315768	48199.75160
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	3904	2.7574902	19012.88985	12.6480309	1998-2012 –> 2013-2019	5.7674504	53119.07347
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	3904	2.7574902	19012.88985	12.6480309	1998-2012 –> 2013-2019	5.7674504	53119.07347
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	3904	2.7574902	19012.88985	12.6480309	1998-2012 –> 2013-2019	5.7674504	53119.07347
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	3904	2.7447079	18976.61186	12.1505620	1998-2012 –> 2013-2019	5.7187379	52920.38660
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	3904	2.7455471	18978.99873	12.2093988	1998-2012 –> 2013-2019	5.7282842	52962.31014
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + aou + silicate	3904	2.7499572	18991.53044	11.4294859	1998-2012 –> 2013-2019	5.5904030	52313.78403
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	3904	2.7432404	18972.43604	12.0931425	1998-2012 –> 2013-2019	5.7177629	52918.45008
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	3904	2.7432404	18972.43604	12.0931425	1998-2012 –> 2013-2019	5.7177629	52918.45008
Indo-Pacific	(27.5,27.75]	4	cstar_tref ~ sal + temp + silicate + phosphate	3904	2.9944894	19654.68141	8.4412807	1998-2012 –> 2013-2019	5.9701436	53603.83986
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	3904	2.7432404	18972.43604	12.0931425	1998-2012 –> 2013-2019	5.7177629	52918.45008
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2653	3.1677939	13660.89373	9.1876445	1982-1997 –> 1998-2012	4.7779117	18678.34293
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	2653	3.1677939	13660.89373	9.1876445	1982-1997 –> 1998-2012	4.7779117	18678.34293
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	2653	3.1677939	13660.89373	9.1876445	1982-1997 –> 1998-2012	4.7779117	18678.34293
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + aou + silicate	2653	3.0964084	13539.95648	9.9082757	1982-1997 –> 1998-2012	4.7755004	18668.14254
Indo-Pacific	(27.75,27.85]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	2653	3.2256771	13754.97200	8.5626959	1982-1997 –> 1998-2012	4.8367490	18771.98496
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	2653	2.9812792	13338.91009	7.7029615	1982-1997 –> 1998-2012	4.5264040	18247.58416
Indo-Pacific	(27.75,27.85]	4	cstar_tref ~ temp + oxygen + aou + phosphate	2653	3.2256771	13754.97200	8.5626959	1982-1997 –> 1998-2012	4.8367490	18771.98496
Indo-Pacific	(27.75,27.85]	4	cstar_tref ~ temp + oxygen + aou + silicate	2653	3.1216319	13581.00425	9.8951151	1982-1997 –> 1998-2012	4.8064931	18716.24560
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	2653	2.9812792	13338.91009	7.7029615	1982-1997 –> 1998-2012	4.5264040	18247.58416
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	2653	2.9812792	13338.91009	7.7029615	1982-1997 –> 1998-2012	4.5264040	18247.58416
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1520	2.4782579	7086.54298	9.5043297	1998-2012 –> 2013-2019	5.3977167	20314.26988
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1520	2.4782579	7086.54298	9.5043297	1998-2012 –> 2013-2019	5.3977167	20314.26988
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1520	2.4782579	7086.54298	9.5043297	1998-2012 –> 2013-2019	5.3977167	20314.26988
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1520	2.4923267	7103.75193	8.9701919	1998-2012 –> 2013-2019	5.4103885	20328.93927
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1520	2.4967804	7109.17935	8.8090180	1998-2012 –> 2013-2019	5.4146842	20334.07960
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	1520	2.4996524	7112.67420	8.7188128	1998-2012 –> 2013-2019	5.4175605	20337.58222
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	1520	2.4996524	7112.67420	8.7188128	1998-2012 –> 2013-2019	5.4175605	20337.58222
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1520	2.4996524	7112.67420	8.7188128	1998-2012 –> 2013-2019	5.4175605	20337.58222
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	1520	2.8722373	7535.05058	8.9679710	1998-2012 –> 2013-2019	5.8535165	20873.96068
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	1520	2.8722373	7535.05058	8.9679710	1998-2012 –> 2013-2019	5.8535165	20873.96068
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	3224	3.2866854	16835.65911	18.0593035	1982-1997 –> 1998-2012	4.7711922	21986.96703
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	3224	3.2866854	16835.65911	18.0593035	1982-1997 –> 1998-2012	4.7711922	21986.96703
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	3224	3.2866854	16835.65911	18.0593035	1982-1997 –> 1998-2012	4.7711922	21986.96703
Indo-Pacific	(27.85,27.95]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	3224	3.3004999	16860.70441	17.8787169	1982-1997 –> 1998-2012	4.7871682	22014.13280
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	3224	3.1399295	16541.11951	17.2962244	1982-1997 –> 1998-2012	4.5928957	21631.60872
Indo-Pacific	(27.85,27.95]	4	cstar_tref ~ temp + oxygen + aou + phosphate	3224	3.3004999	16860.70441	17.8787169	1982-1997 –> 1998-2012	4.7871682	22014.13280
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ temp + oxygen + aou + phosphate + phosphate_star	3224	3.3004999	16860.70441	17.8787169	1982-1997 –> 1998-2012	4.7871682	22014.13280
Indo-Pacific	(27.85,27.95]	4	cstar_tref ~ temp + oxygen + aou + phosphate_star	3224	3.3004999	16860.70441	17.8787169	1982-1997 –> 1998-2012	4.7871682	22014.13280
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	3224	3.1399295	16541.11951	17.2962244	1982-1997 –> 1998-2012	4.5928957	21631.60872
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	3224	3.1399295	16541.11951	17.2962244	1982-1997 –> 1998-2012	4.5928957	21631.60872
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1823	3.3105153	9552.09060	9.3111647	1998-2012 –> 2013-2019	6.4466295	26085.37037
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1823	3.3105153	9552.09060	9.3111647	1998-2012 –> 2013-2019	6.4466295	26085.37037
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1823	3.3105153	9552.09060	9.3111647	1998-2012 –> 2013-2019	6.4466295	26085.37037
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1823	3.3444727	9589.29868	9.4041251	1998-2012 –> 2013-2019	6.4832831	26128.11960
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1823	3.3596636	9605.82162	9.4523390	1998-2012 –> 2013-2019	6.5012991	26150.44347
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	1823	3.3664542	9613.18354	9.4679657	1998-2012 –> 2013-2019	6.5093236	26160.33752
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	1823	3.3664542	9613.18354	9.4679657	1998-2012 –> 2013-2019	6.5093236	26160.33752
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1823	3.3664542	9613.18354	9.4679657	1998-2012 –> 2013-2019	6.5093236	26160.33752
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1823	3.4615597	9714.75830	9.8097748	1998-2012 –> 2013-2019	6.6014892	26255.87781
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	1823	3.4615597	9714.75830	9.8097748	1998-2012 –> 2013-2019	6.6014892	26255.87781
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	3905	1.6902874	15195.36784	8.6171219	1982-1997 –> 1998-2012	3.3469292	21969.35876
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	3905	1.6902874	15195.36784	8.6171219	1982-1997 –> 1998-2012	3.3469292	21969.35876
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	3905	1.6902874	15195.36784	8.6171219	1982-1997 –> 1998-2012	3.3469292	21969.35876
Indo-Pacific	(27.95,28.05]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	3905	1.7423179	15430.14983	10.5368511	1982-1997 –> 1998-2012	3.4014979	22207.52049
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	3905	1.6583558	15046.41573	9.0958786	1982-1997 –> 1998-2012	3.1955837	21557.51817
Indo-Pacific	(27.95,28.05]	4	cstar_tref ~ temp + oxygen + aou + phosphate	3905	1.7423179	15430.14983	10.5368511	1982-1997 –> 1998-2012	3.4014979	22207.52049
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ temp + oxygen + aou + phosphate + phosphate_star	3905	1.7423179	15430.14983	10.5368511	1982-1997 –> 1998-2012	3.4014979	22207.52049
Indo-Pacific	(27.95,28.05]	4	cstar_tref ~ temp + oxygen + aou + phosphate_star	3905	1.7423179	15430.14983	10.5368511	1982-1997 –> 1998-2012	3.4014979	22207.52049
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	3905	1.6583558	15046.41573	9.0958786	1982-1997 –> 1998-2012	3.1955837	21557.51817
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	3905	1.6583558	15046.41573	9.0958786	1982-1997 –> 1998-2012	3.1955837	21557.51817
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2207	1.5684161	8263.78728	9.1231543	1998-2012 –> 2013-2019	3.2519158	23427.72904
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	2207	1.5684161	8263.78728	9.1231543	1998-2012 –> 2013-2019	3.2519158	23427.72904
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	2207	1.5684161	8263.78728	9.1231543	1998-2012 –> 2013-2019	3.2519158	23427.72904
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	2207	1.5711791	8271.55624	9.0746907	1998-2012 –> 2013-2019	3.2682088	23498.01482
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	2207	1.5719412	8273.69666	9.0469574	1998-2012 –> 2013-2019	3.2735410	23521.15955
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	2207	1.5719951	8273.84817	9.0399522	1998-2012 –> 2013-2019	3.2749070	23527.33053
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	2207	1.5719951	8273.84817	9.0399522	1998-2012 –> 2013-2019	3.2749070	23527.33053
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	2207	1.5719951	8273.84817	9.0399522	1998-2012 –> 2013-2019	3.2749070	23527.33053
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	2207	1.6205132	8408.02173	9.7347045	1998-2012 –> 2013-2019	3.2788689	23454.43746
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	2207	1.6205132	8408.02173	9.7347045	1998-2012 –> 2013-2019	3.2788689	23454.43746
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2237	1.5455270	8310.15381	7.2890494	1982-1997 –> 1998-2012	3.1117339	12365.62994
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	2237	1.5455270	8310.15381	7.2890494	1982-1997 –> 1998-2012	3.1117339	12365.62994
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	2237	1.5455270	8310.15381	7.2890494	1982-1997 –> 1998-2012	3.1117339	12365.62994
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	2237	1.5516233	8327.76665	7.3523826	1982-1997 –> 1998-2012	3.1283321	12397.70458
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	2237	1.5533951	8332.87273	7.3766845	1982-1997 –> 1998-2012	3.1337101	12407.75443
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	2237	1.5533951	8332.87273	7.3766845	1982-1997 –> 1998-2012	3.1337101	12407.75443
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	2237	1.5533951	8332.87273	7.3766845	1982-1997 –> 1998-2012	3.1337101	12407.75443
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	2237	1.5214362	8239.86632	6.1811399	1982-1997 –> 1998-2012	2.9802631	12141.64629
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	2237	1.5214362	8239.86632	6.1811399	1982-1997 –> 1998-2012	2.9802631	12141.64629
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	2237	1.5214362	8239.86632	6.1811399	1982-1997 –> 1998-2012	2.9802631	12141.64629
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1247	1.2902622	4188.41727	5.8262405	1998-2012 –> 2013-2019	2.8357892	12498.57108
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1247	1.2902622	4188.41727	5.8262405	1998-2012 –> 2013-2019	2.8357892	12498.57108
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1247	1.2902622	4188.41727	5.8262405	1998-2012 –> 2013-2019	2.8357892	12498.57108
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1247	1.2934477	4194.56700	5.7870805	1998-2012 –> 2013-2019	2.8450710	12522.33365
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	1247	1.2941037	4195.83169	5.7686036	1998-2012 –> 2013-2019	2.8474989	12528.70442
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	1247	1.2941037	4195.83169	5.7686036	1998-2012 –> 2013-2019	2.8474989	12528.70442
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1247	1.2941037	4195.83169	5.7686036	1998-2012 –> 2013-2019	2.8474989	12528.70442
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1247	1.2940968	4195.81825	6.4737894	1998-2012 –> 2013-2019	2.8155329	12435.68457
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	1247	1.2940968	4195.81825	6.4737894	1998-2012 –> 2013-2019	2.8155329	12435.68457
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	1247	1.2940968	4195.81825	6.4737894	1998-2012 –> 2013-2019	2.8155329	12435.68457
Indo-Pacific	(28.1, Inf]	4	cstar_tref ~ sal + aou + silicate + phosphate	24288	1.5795391	91144.05477	10.1415860	1982-1997 –> 1998-2012	3.0059144	135330.52210
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	24288	1.5388978	89879.84132	9.6263790	1982-1997 –> 1998-2012	2.8571433	132105.32317
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	24288	1.5388978	89879.84132	9.6263790	1982-1997 –> 1998-2012	2.8571433	132105.32317
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	24288	1.5388978	89879.84132	9.6263790	1982-1997 –> 1998-2012	2.8571433	132105.32317
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	24288	1.5390476	89884.57128	9.6276937	1982-1997 –> 1998-2012	2.8564425	132093.97968
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	24288	1.5392039	89889.50373	9.6288941	1982-1997 –> 1998-2012	2.8555507	132079.09450
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	24288	1.5391920	89889.12692	9.6291400	1982-1997 –> 1998-2012	2.8555803	132079.50355
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	24288	1.5391920	89889.12692	9.6291400	1982-1997 –> 1998-2012	2.8555803	132079.50355
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	24288	1.5391920	89889.12692	9.6291400	1982-1997 –> 1998-2012	2.8555803	132079.50355
Indo-Pacific	(28.1, Inf]	4	cstar_tref ~ sal + temp + silicate + phosphate_star	24288	1.5464842	90116.72102	9.2837867	1982-1997 –> 1998-2012	2.8747338	132528.45475
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	14824	1.2904044	49641.61538	9.2901765	1998-2012 –> 2013-2019	2.8293022	139521.45670
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	14824	1.2904044	49641.61538	9.2901765	1998-2012 –> 2013-2019	2.8293022	139521.45670
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	14824	1.2904044	49641.61538	9.2901765	1998-2012 –> 2013-2019	2.8293022	139521.45670
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	14824	1.2907334	49649.17197	9.2905994	1998-2012 –> 2013-2019	2.8297810	139533.74325
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	14824	1.2911776	49659.37507	9.2937847	1998-2012 –> 2013-2019	2.8303815	139548.87880
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	14824	1.2911324	49658.33768	9.2927061	1998-2012 –> 2013-2019	2.8303244	139547.46460
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	14824	1.2911324	49658.33768	9.2927061	1998-2012 –> 2013-2019	2.8303244	139547.46460
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	14824	1.2911324	49658.33768	9.2927061	1998-2012 –> 2013-2019	2.8303244	139547.46460
Indo-Pacific	(28.1, Inf]	4	cstar_tref ~ sal + temp + silicate + phosphate_star	14824	1.3186734	50282.10445	10.1511497	1998-2012 –> 2013-2019	2.8651576	140398.82547
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	14824	1.3443327	50855.46644	9.7413986	1998-2012 –> 2013-2019	2.9195345	141867.95131

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-1997 –> 1998-2012	2.2783426	306.47950
Atlantic	(-Inf,26]	1998-2012 –> 2013-2019	0.6980415	74.31315
Atlantic	(26,26.5]	1982-1997 –> 1998-2012	7.2781186	3823.26169
Atlantic	(26,26.5]	1998-2012 –> 2013-2019	6.6699127	3365.76334
Atlantic	(26.5,26.75]	1982-1997 –> 1998-2012	5.6744549	4314.89681
Atlantic	(26.5,26.75]	1998-2012 –> 2013-2019	5.6184265	3681.36944
Atlantic	(26.75,27]	1982-1997 –> 1998-2012	4.0376182	6714.68562
Atlantic	(26.75,27]	1998-2012 –> 2013-2019	5.3376003	6131.06510
Atlantic	(27,27.25]	1982-1997 –> 1998-2012	3.8022063	8031.82013
Atlantic	(27,27.25]	1998-2012 –> 2013-2019	4.2439513	7284.65816
Atlantic	(27.25,27.5]	1982-1997 –> 1998-2012	4.7921039	15460.97681
Atlantic	(27.25,27.5]	1998-2012 –> 2013-2019	6.7843306	15545.97210
Atlantic	(27.5,27.75]	1982-1997 –> 1998-2012	3.9675666	26779.41945
Atlantic	(27.5,27.75]	1998-2012 –> 2013-2019	4.6258820	24746.01629
Atlantic	(27.75,27.85]	1982-1997 –> 1998-2012	1.9019699	6813.25053
Atlantic	(27.75,27.85]	1998-2012 –> 2013-2019	2.3362104	6691.93953
Atlantic	(27.85,27.95]	1982-1997 –> 1998-2012	2.5125069	8088.17131
Atlantic	(27.85,27.95]	1998-2012 –> 2013-2019	2.5726400	6775.62606
Atlantic	(27.95,28.05]	1982-1997 –> 1998-2012	8.1964564	16980.93530
Atlantic	(27.95,28.05]	1998-2012 –> 2013-2019	5.3626321	11822.75829
Atlantic	(28.05,28.1]	1982-1997 –> 1998-2012	1.8305778	4677.53990
Atlantic	(28.05,28.1]	1998-2012 –> 2013-2019	1.6414736	3743.68089
Atlantic	(28.1,28.15]	1982-1997 –> 1998-2012	1.6055777	5676.11760
Atlantic	(28.1,28.15]	1998-2012 –> 2013-2019	1.4634026	4917.71495
Atlantic	(28.15,28.2]	1982-1997 –> 1998-2012	1.6476992	10004.69317
Atlantic	(28.15,28.2]	1998-2012 –> 2013-2019	1.4321182	8899.34669
Atlantic	(28.2, Inf]	1982-1997 –> 1998-2012	1.7016351	18768.95437
Atlantic	(28.2, Inf]	1998-2012 –> 2013-2019	1.4709498	15892.30879
Indo-Pacific	(-Inf,26]	1982-1997 –> 1998-2012	10.9776062	23715.10839
Indo-Pacific	(-Inf,26]	1998-2012 –> 2013-2019	11.3500470	30200.78444
Indo-Pacific	(26,26.5]	1982-1997 –> 1998-2012	8.6217226	24657.78984
Indo-Pacific	(26,26.5]	1998-2012 –> 2013-2019	8.5332742	23867.96627
Indo-Pacific	(26.5,26.75]	1982-1997 –> 1998-2012	7.9593299	21064.53938
Indo-Pacific	(26.5,26.75]	1998-2012 –> 2013-2019	7.9266707	21052.76689
Indo-Pacific	(26.75,27]	1982-1997 –> 1998-2012	7.5377850	29699.92807
Indo-Pacific	(26.75,27]	1998-2012 –> 2013-2019	8.2055942	29805.05035
Indo-Pacific	(27,27.25]	1982-1997 –> 1998-2012	6.9633836	36513.09740
Indo-Pacific	(27,27.25]	1998-2012 –> 2013-2019	8.2416230	40657.82311
Indo-Pacific	(27.25,27.5]	1982-1997 –> 1998-2012	5.8734207	44148.41157
Indo-Pacific	(27.25,27.5]	1998-2012 –> 2013-2019	6.7034212	48625.79384
Indo-Pacific	(27.5,27.75]	1982-1997 –> 1998-2012	5.2396349	48425.93175
Indo-Pacific	(27.5,27.75]	1998-2012 –> 2013-2019	5.7463209	52991.28913
Indo-Pacific	(27.75,27.85]	1982-1997 –> 1998-2012	4.7168439	18570.61393
Indo-Pacific	(27.75,27.85]	1998-2012 –> 2013-2019	5.4977937	20436.64965
Indo-Pacific	(27.85,27.95]	1982-1997 –> 1998-2012	4.7240936	21891.22585
Indo-Pacific	(27.85,27.95]	1998-2012 –> 2013-2019	6.5055420	26152.74423
Indo-Pacific	(27.95,28.05]	1982-1997 –> 1998-2012	3.3233530	21941.07127
Indo-Pacific	(27.95,28.05]	1998-2012 –> 2013-2019	3.2679956	23479.32280
Indo-Pacific	(28.05,28.1]	1982-1997 –> 1998-2012	3.0805453	12314.27965
Indo-Pacific	(28.05,28.1]	1998-2012 –> 2013-2019	2.8341534	12491.12139
Indo-Pacific	(28.1, Inf]	1982-1997 –> 1998-2012	2.8730812	132458.65312
Indo-Pacific	(28.1, Inf]	1998-2012 –> 2013-2019	2.8423734	139855.61627

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
aecbf75	Donghe-Zhu	2021-03-14
1b2a0c2	Donghe-Zhu	2021-03-14
6733e48	Donghe-Zhu	2021-03-12
ba71e6a	Donghe-Zhu	2021-03-12
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
aecbf75	Donghe-Zhu	2021-03-14
1b2a0c2	Donghe-Zhu	2021-03-14
6733e48	Donghe-Zhu	2021-03-12
ba71e6a	Donghe-Zhu	2021-03-12
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
aecbf75	Donghe-Zhu	2021-03-14
1b2a0c2	Donghe-Zhu	2021-03-14
6733e48	Donghe-Zhu	2021-03-12
ba71e6a	Donghe-Zhu	2021-03-12
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
aecbf75	Donghe-Zhu	2021-03-14
1b2a0c2	Donghe-Zhu	2021-03-14
6733e48	Donghe-Zhu	2021-03-12
ba71e6a	Donghe-Zhu	2021-03-12
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-1997 --> 1998-2012
(-Inf,26]	7	6	6	6	7	8	10
(26,26.5]	7	5	7	6	8	10	7
(26.5,26.75]	8	5	6	7	7	10	7
(26.75,27]	10	7	7	6	3	10	3
(27,27.25]	7	5	6	6	9	10	6
(27.25,27.5]	7	4	7	6	10	9	6
(27.5,27.75]	6	4	6	7	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	6	6	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	7	6	7	6
(28.2, Inf]	7	4	7	6	9	10	6
total	113.00	78.00	93.00	88.00	101.00	110.00	99.00
Atlantic - 1998-2012 --> 2013-2019
(-Inf,26]	7	6	5	6	7	10	9
(26,26.5]	7	5	6	7	8	10	7
(26.5,26.75]	8	5	6	7	7	10	7
(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	6	7	6	10	7	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]	10	7	7	6	3	3	10
total	114.00	78.00	90.00	86.00	102.00	107.00	104.00
Indo-Pacific - 1982-1997 --> 1998-2012
(-Inf,26]	10	6	8	5	10	4	4
(26,26.5]	7	6	6	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]	7	7	5	5	6	10	9
(27.5,27.75]	7	5	6	6	7	10	9
(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	104.00	75.00	79.00	69.00	80.00	76.00	92.00
Indo-Pacific - 1998-2012 --> 2013-2019
(-Inf,26]	10	7	7	6	10	3	3
(26,26.5]	7	6	6	6	8	10	7
(26.5,26.75]	10	6	7	6	7	7	7
(26.75,27]	10	7	6	7	6	6	7
(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	6	8	10	7
(27.85,27.95]	7	5	7	6	8	10	7
(27.95,28.05]	7	5	7	6	8	10	7
(28.05,28.1]	7	6	7	6	7	10	7
(28.1, Inf]	6	4	6	7	9	10	7
total	94.00	70.00	77.00	72.00	90.00	99.00	86.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
aecbf75	Donghe-Zhu	2021-03-14
1b2a0c2	Donghe-Zhu	2021-03-14
6733e48	Donghe-Zhu	2021-03-12
ba71e6a	Donghe-Zhu	2021-03-12
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
aecbf75	Donghe-Zhu	2021-03-14
1b2a0c2	Donghe-Zhu	2021-03-14
6733e48	Donghe-Zhu	2021-03-12
ba71e6a	Donghe-Zhu	2021-03-12
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