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: 49d7033

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 49d7033. 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	f745381	Donghe-Zhu	2021-03-16	Build site.
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 + phosphate + phosphate_star	26	0.6987083	69.14167	2.2173182	1982-1999 –> 2000-2009	2.2138763	273.9260
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	26	0.6987083	69.14167	2.2173182	1982-1999 –> 2000-2009	2.2138763	273.9260
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	26	0.6987083	69.14167	2.2173182	1982-1999 –> 2000-2009	2.2138763	273.9260
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + aou + silicate	26	0.4865676	50.32507	1.1237168	1982-1999 –> 2000-2009	2.1237705	263.1656
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	26	0.4584940	47.23478	1.1190183	1982-1999 –> 2000-2009	2.2860715	271.5155
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	26	0.4584940	47.23478	1.1190183	1982-1999 –> 2000-2009	2.2860715	271.5155
Atlantic	(-Inf,26]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	26	0.5133095	53.10723	0.9152479	1982-1999 –> 2000-2009	1.8797173	247.1441
Atlantic	(-Inf,26]	4	cstar_tref ~ temp + oxygen + aou + silicate	26	0.5135527	51.13187	0.9115320	1982-1999 –> 2000-2009	2.2338454	267.1209
Atlantic	(-Inf,26]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	26	0.5133095	53.10723	0.9152479	1982-1999 –> 2000-2009	1.8797173	247.1441
Atlantic	(-Inf,26]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	26	0.5133095	53.10723	0.9152479	1982-1999 –> 2000-2009	1.8797173	247.1441
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	51	0.4665504	80.96802	1.1999270	2000-2009 –> 2010-2019	0.9278469	128.5197
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	51	0.4665504	80.96802	1.1999270	2000-2009 –> 2010-2019	0.9278469	128.5197
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	51	0.4585669	79.20752	1.1603831	2000-2009 –> 2010-2019	0.9174727	126.4890
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	51	0.4609553	79.73739	1.1668857	2000-2009 –> 2010-2019	0.9194428	126.9714
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	51	0.4577744	79.03110	1.1530624	2000-2009 –> 2010-2019	0.9162684	126.2659
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	51	0.4577744	79.03110	1.1530624	2000-2009 –> 2010-2019	0.9162684	126.2659
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	51	0.4577744	79.03110	1.1530624	2000-2009 –> 2010-2019	0.9162684	126.2659
Atlantic	(-Inf,26]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	51	0.4032687	66.10020	0.9135637	2000-2009 –> 2010-2019	0.9165782	119.2074
Atlantic	(-Inf,26]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	51	0.4032687	66.10020	0.9135637	2000-2009 –> 2010-2019	0.9165782	119.2074
Atlantic	(-Inf,26]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	51	0.4032687	66.10020	0.9135637	2000-2009 –> 2010-2019	0.9165782	119.2074
Atlantic	(26,26.5]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	195	4.2121666	1128.19712	10.4631952	1982-1999 –> 2000-2009	7.2952003	2877.7866
Atlantic	(26,26.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	195	4.2121666	1128.19712	10.4631952	1982-1999 –> 2000-2009	7.2952003	2877.7866
Atlantic	(26,26.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	195	4.2121666	1128.19712	10.4631952	1982-1999 –> 2000-2009	7.2952003	2877.7866
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	195	4.3050687	1136.70533	10.9027558	1982-1999 –> 2000-2009	7.4906122	2908.6022
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	195	4.3272566	1138.71019	11.0127213	1982-1999 –> 2000-2009	7.5320449	2914.7149
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	195	4.3360442	1139.50139	10.9822045	1982-1999 –> 2000-2009	7.5483565	2917.1053
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	195	4.3360442	1139.50139	10.9822045	1982-1999 –> 2000-2009	7.5483565	2917.1053
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	195	4.3360442	1139.50139	10.9822045	1982-1999 –> 2000-2009	7.5483565	2917.1053
Atlantic	(26,26.5]	4	cstar_tref ~ sal + temp + silicate + phosphate_star	195	4.5800495	1158.85285	10.2362288	1982-1999 –> 2000-2009	8.2098184	3017.7822
Atlantic	(26,26.5]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	195	4.5249170	1156.12973	12.0222034	1982-1999 –> 2000-2009	8.2642811	3037.3462
Atlantic	(26,26.5]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	359	3.5711147	1946.72411	12.5903544	2000-2009 –> 2010-2019	7.7832813	3074.9212
Atlantic	(26,26.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	359	3.5711147	1946.72411	12.5903544	2000-2009 –> 2010-2019	7.7832813	3074.9212
Atlantic	(26,26.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	359	3.5711147	1946.72411	12.5903544	2000-2009 –> 2010-2019	7.7832813	3074.9212
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	359	3.5848982	1949.49007	12.5733578	2000-2009 –> 2010-2019	7.8899669	3086.1954
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	359	3.5871576	1949.94245	12.5685911	2000-2009 –> 2010-2019	7.9144142	3088.6526
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + aou + silicate	359	3.7496452	1981.75063	12.8302731	2000-2009 –> 2010-2019	8.3164141	3141.4710
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	359	3.5886836	1950.24782	12.5558472	2000-2009 –> 2010-2019	7.9247278	3089.7492
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	359	3.5886836	1950.24782	12.5558472	2000-2009 –> 2010-2019	7.9247278	3089.7492
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	359	3.5886836	1950.24782	12.5558472	2000-2009 –> 2010-2019	7.9247278	3089.7492
Atlantic	(26,26.5]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	359	3.8229271	1995.64763	13.7338573	2000-2009 –> 2010-2019	8.3478441	3151.7774
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	266	2.9772616	1349.28940	7.0206396	1982-1999 –> 2000-2009	5.6952281	3187.0816
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	266	2.9772616	1349.28940	7.0206396	1982-1999 –> 2000-2009	5.6952281	3187.0816
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	266	2.9772616	1349.28940	7.0206396	1982-1999 –> 2000-2009	5.6952281	3187.0816
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	266	2.9759262	1349.05073	7.0393155	1982-1999 –> 2000-2009	5.7498404	3202.2060
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	266	2.9763191	1349.12095	7.1250744	1982-1999 –> 2000-2009	5.7840058	3211.4008
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	266	2.9768698	1349.21938	7.1462711	1982-1999 –> 2000-2009	5.7965255	3214.7066
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	266	2.9768698	1349.21938	7.1462711	1982-1999 –> 2000-2009	5.7965255	3214.7066
Atlantic	(26.5,26.75]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	266	3.1255830	1375.15358	10.1985469	1982-1999 –> 2000-2009	5.7614442	3189.8176
Atlantic	(26.5,26.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	266	3.1255830	1375.15358	10.1985469	1982-1999 –> 2000-2009	5.7614442	3189.8176
Atlantic	(26.5,26.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	266	3.1255830	1375.15358	10.1985469	1982-1999 –> 2000-2009	5.7614442	3189.8176
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	433	4.3897736	2523.85521	11.0515632	2000-2009 –> 2010-2019	7.3670352	3873.1446
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	433	4.3897736	2523.85521	11.0515632	2000-2009 –> 2010-2019	7.3670352	3873.1446
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	433	4.3897736	2523.85521	11.0515632	2000-2009 –> 2010-2019	7.3670352	3873.1446
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	433	4.5002079	2545.37180	11.0513937	2000-2009 –> 2010-2019	7.4761341	3894.4225
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	433	4.5250345	2550.13620	11.0968834	2000-2009 –> 2010-2019	7.5013536	3899.2572
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	433	4.5291329	2550.92019	11.1095627	2000-2009 –> 2010-2019	7.5060027	3900.1396
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	433	4.5291329	2550.92019	11.1095627	2000-2009 –> 2010-2019	7.5060027	3900.1396
Atlantic	(26.5,26.75]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	433	4.3096222	2507.89707	11.0498463	2000-2009 –> 2010-2019	7.4352052	3883.0507
Atlantic	(26.5,26.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	433	4.3096222	2507.89707	11.0498463	2000-2009 –> 2010-2019	7.4352052	3883.0507
Atlantic	(26.5,26.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	433	4.3096222	2507.89707	11.0498463	2000-2009 –> 2010-2019	7.4352052	3883.0507
Atlantic	(26.75,27]	4	cstar_tref ~ aou + silicate + phosphate + phosphate_star	437	2.3258848	1989.89614	10.1912322	1982-1999 –> 2000-2009	4.0464360	5050.1965
Atlantic	(26.75,27]	4	cstar_tref ~ oxygen + aou + silicate + phosphate	437	2.3258848	1989.89614	10.1912322	1982-1999 –> 2000-2009	4.0464360	5050.1965
Atlantic	(26.75,27]	5	cstar_tref ~ oxygen + aou + silicate + phosphate + phosphate_star	437	2.3258848	1989.89614	10.1912322	1982-1999 –> 2000-2009	4.0464360	5050.1965
Atlantic	(26.75,27]	4	cstar_tref ~ oxygen + aou + silicate + phosphate_star	437	2.3258848	1989.89614	10.1912322	1982-1999 –> 2000-2009	4.0464360	5050.1965
Atlantic	(26.75,27]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	437	2.2874229	1977.32245	9.5012148	1982-1999 –> 2000-2009	4.0072407	5038.9602
Atlantic	(26.75,27]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	437	2.2874229	1977.32245	9.5012148	1982-1999 –> 2000-2009	4.0072407	5038.9602
Atlantic	(26.75,27]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	437	2.2874229	1977.32245	9.5012148	1982-1999 –> 2000-2009	4.0072407	5038.9602
Atlantic	(26.75,27]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	437	2.2192662	1950.88462	9.5895072	1982-1999 –> 2000-2009	3.8043227	4885.7187
Atlantic	(26.75,27]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	437	2.2192662	1950.88462	9.5895072	1982-1999 –> 2000-2009	3.8043227	4885.7187
Atlantic	(26.75,27]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	437	2.2192662	1950.88462	9.5895072	1982-1999 –> 2000-2009	3.8043227	4885.7187
Atlantic	(26.75,27]	4	cstar_tref ~ aou + silicate + phosphate + phosphate_star	796	3.0110279	4025.78231	14.3038724	2000-2009 –> 2010-2019	5.3369127	6015.6785
Atlantic	(26.75,27]	4	cstar_tref ~ oxygen + aou + silicate + phosphate	796	3.0110279	4025.78231	14.3038724	2000-2009 –> 2010-2019	5.3369127	6015.6785
Atlantic	(26.75,27]	5	cstar_tref ~ oxygen + aou + silicate + phosphate + phosphate_star	796	3.0110279	4025.78231	14.3038724	2000-2009 –> 2010-2019	5.3369127	6015.6785
Atlantic	(26.75,27]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	796	2.8636822	3947.90654	13.8821051	2000-2009 –> 2010-2019	5.1511051	5925.2290
Atlantic	(26.75,27]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	796	2.8636822	3947.90654	13.8821051	2000-2009 –> 2010-2019	5.1511051	5925.2290
Atlantic	(26.75,27]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	796	2.8636822	3947.90654	13.8821051	2000-2009 –> 2010-2019	5.1511051	5925.2290
Atlantic	(26.75,27]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	796	2.8777666	3955.71724	13.9254924	2000-2009 –> 2010-2019	5.2914306	5979.9911
Atlantic	(26.75,27]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	796	2.9964500	4020.05592	14.0893004	2000-2009 –> 2010-2019	5.2157162	5970.9405
Atlantic	(26.75,27]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	796	2.9964500	4020.05592	14.0893004	2000-2009 –> 2010-2019	5.2157162	5970.9405
Atlantic	(26.75,27]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	796	2.9964500	4020.05592	14.0893004	2000-2009 –> 2010-2019	5.2157162	5970.9405
Atlantic	(27,27.25]	4	cstar_tref ~ sal + aou + silicate + phosphate	534	2.1918137	2365.51732	9.9205104	1982-1999 –> 2000-2009	3.7044341	5962.4390
Atlantic	(27,27.25]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	534	2.1795915	2361.54519	10.2201822	1982-1999 –> 2000-2009	3.5571409	5777.5077
Atlantic	(27,27.25]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	534	2.1795915	2361.54519	10.2201822	1982-1999 –> 2000-2009	3.5571409	5777.5077
Atlantic	(27,27.25]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	534	2.1795915	2361.54519	10.2201822	1982-1999 –> 2000-2009	3.5571409	5777.5077
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	534	2.1878645	2365.59127	10.0924119	1982-1999 –> 2000-2009	3.5972742	5826.2775
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	534	2.2230194	2382.61561	10.1018365	1982-1999 –> 2000-2009	3.6760617	5902.9374
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	534	2.2315267	2386.69498	10.1078228	1982-1999 –> 2000-2009	3.6973902	5924.2001
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	534	2.2315267	2386.69498	10.1078228	1982-1999 –> 2000-2009	3.6973902	5924.2001
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	534	2.2315267	2386.69498	10.1078228	1982-1999 –> 2000-2009	3.6973902	5924.2001
Atlantic	(27,27.25]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	534	2.6509858	2570.65327	12.5638517	1982-1999 –> 2000-2009	4.2369298	6262.1648
Atlantic	(27,27.25]	4	cstar_tref ~ sal + aou + silicate + phosphate	1084	2.2404859	4837.16868	11.8788182	2000-2009 –> 2010-2019	4.4322996	7202.6860
Atlantic	(27,27.25]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1084	2.2339645	4832.84908	11.4935634	2000-2009 –> 2010-2019	4.4135561	7194.3943
Atlantic	(27,27.25]	4	cstar_tref ~ sal + aou + silicate + phosphate_star	1084	2.2654702	4861.21082	10.6934635	2000-2009 –> 2010-2019	4.6404155	7312.4299
Atlantic	(27,27.25]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1084	2.2339645	4832.84908	11.4935634	2000-2009 –> 2010-2019	4.4135561	7194.3943
Atlantic	(27,27.25]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1084	2.2339645	4832.84908	11.4935634	2000-2009 –> 2010-2019	4.4135561	7194.3943
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1084	2.2239435	4823.10204	11.2760258	2000-2009 –> 2010-2019	4.4118080	7188.6933
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1084	2.2487572	4847.15759	11.2679049	2000-2009 –> 2010-2019	4.4717765	7229.7732
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	1084	2.2531889	4851.42600	11.2583050	2000-2009 –> 2010-2019	4.4847157	7238.1210
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	1084	2.2531889	4851.42600	11.2583050	2000-2009 –> 2010-2019	4.4847157	7238.1210
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1084	2.2531889	4851.42600	11.2583050	2000-2009 –> 2010-2019	4.4847157	7238.1210
Atlantic	(27.25,27.5]	4	cstar_tref ~ sal + aou + silicate + phosphate	1226	3.6686880	6678.43054	15.0922088	1982-1999 –> 2000-2009	5.1560949	11553.6057
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1226	3.6634817	6676.94836	14.6638844	1982-1999 –> 2000-2009	4.9798573	11227.0002
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1226	3.6634817	6676.94836	14.6638844	1982-1999 –> 2000-2009	4.9798573	11227.0002
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1226	3.6634817	6676.94836	14.6638844	1982-1999 –> 2000-2009	4.9798573	11227.0002
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1226	3.6135880	6643.32461	14.0194988	1982-1999 –> 2000-2009	4.8586781	11044.2807
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1226	3.4687265	6543.00427	15.3660669	1982-1999 –> 2000-2009	4.7579524	11037.2461
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	1226	3.5273438	6584.09391	14.7080632	1982-1999 –> 2000-2009	5.2510796	11856.1660
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	1226	3.4735395	6546.40415	15.2642225	1982-1999 –> 2000-2009	4.7786430	11073.4257
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	1226	3.4735395	6546.40415	15.2642225	1982-1999 –> 2000-2009	4.7786430	11073.4257
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1226	3.4735395	6546.40415	15.2642225	1982-1999 –> 2000-2009	4.7786430	11073.4257
Atlantic	(27.25,27.5]	4	cstar_tref ~ sal + aou + silicate + phosphate	1645	3.2253845	8533.06937	9.8912852	2000-2009 –> 2010-2019	6.8940725	15211.4999
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1645	3.1350336	8441.59301	8.9810334	2000-2009 –> 2010-2019	6.7985153	15118.5414
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1645	3.1350336	8441.59301	8.9810334	2000-2009 –> 2010-2019	6.7985153	15118.5414
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1645	3.1350336	8441.59301	8.9810334	2000-2009 –> 2010-2019	6.7985153	15118.5414
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1645	3.1971310	8506.12288	9.6884441	2000-2009 –> 2010-2019	6.8107190	15149.4475
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1645	3.0963777	8400.77419	8.8960601	2000-2009 –> 2010-2019	6.5651042	14943.7785
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	1645	3.0979778	8402.47392	8.9128976	2000-2009 –> 2010-2019	6.5715173	14948.8781
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	1645	3.0979778	8402.47392	8.9128976	2000-2009 –> 2010-2019	6.5715173	14948.8781
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1645	3.0979778	8402.47392	8.9128976	2000-2009 –> 2010-2019	6.5715173	14948.8781
Atlantic	(27.25,27.5]	4	cstar_tref ~ sal + temp + silicate + phosphate_star	1645	3.1097254	8412.92608	9.0183388	2000-2009 –> 2010-2019	6.7366857	15063.3077
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2434	2.4481074	11279.78726	18.7421600	1982-1999 –> 2000-2009	4.0929523	21823.5032
Atlantic	(27.5,27.75]	4	cstar_tref ~ sal + temp + aou + phosphate_star	2434	2.4698976	11320.92492	16.7871668	1982-1999 –> 2000-2009	4.1194297	21878.2743
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	2434	2.5483683	11475.17936	15.2959747	1982-1999 –> 2000-2009	4.0723901	21599.7398
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	2434	2.4492495	11282.05792	17.8121100	1982-1999 –> 2000-2009	3.8712952	21026.1241
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	2434	2.4481074	11279.78726	18.7421600	1982-1999 –> 2000-2009	4.0929523	21823.5032
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	2434	2.4481074	11279.78726	18.7421600	1982-1999 –> 2000-2009	4.0929523	21823.5032
Atlantic	(27.5,27.75]	4	cstar_tref ~ sal + temp + oxygen + phosphate	2434	2.4802681	11341.32168	16.4005360	1982-1999 –> 2000-2009	4.1302755	21900.2539
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	2434	2.4585077	11300.42425	17.5298502	1982-1999 –> 2000-2009	3.8797455	21041.3684
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	2434	2.4585077	11300.42425	17.5298502	1982-1999 –> 2000-2009	3.8797455	21041.3684
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	2434	2.4585077	11300.42425	17.5298502	1982-1999 –> 2000-2009	3.8797455	21041.3684
Atlantic	(27.5,27.75]	4	cstar_tref ~ sal + aou + silicate + phosphate	2864	2.0768817	12326.08939	9.7450563	2000-2009 –> 2010-2019	4.6638558	23872.4622
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2864	2.0588081	12278.02456	9.7244235	2000-2009 –> 2010-2019	4.6363854	23808.6829
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	2864	2.0588081	12278.02456	9.7244235	2000-2009 –> 2010-2019	4.6363854	23808.6829
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	2864	2.0588081	12278.02456	9.7244235	2000-2009 –> 2010-2019	4.6363854	23808.6829
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2864	2.3764718	13099.93385	12.4810064	2000-2009 –> 2010-2019	4.8245792	24379.7211
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	2864	2.0621401	12287.28732	10.0344514	2000-2009 –> 2010-2019	4.6105084	23762.4667
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	2864	2.0682044	12304.10759	9.2632678	2000-2009 –> 2010-2019	4.5174540	23586.1655
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	2864	2.0663935	12299.08990	9.2140915	2000-2009 –> 2010-2019	4.5249012	23599.5142
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	2864	2.0663935	12299.08990	9.2140915	2000-2009 –> 2010-2019	4.5249012	23599.5142
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	2864	2.0663935	12299.08990	9.2140915	2000-2009 –> 2010-2019	4.5249012	23599.5142
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1015	1.0732592	3037.96620	5.4035129	1982-1999 –> 2000-2009	1.7330935	4988.0882
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1015	1.0732592	3037.96620	5.4035129	1982-1999 –> 2000-2009	1.7330935	4988.0882
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1015	1.0732592	3037.96620	5.4035129	1982-1999 –> 2000-2009	1.7330935	4988.0882
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	1015	0.9198056	2724.75148	3.9962002	1982-1999 –> 2000-2009	1.6418127	4848.6627
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1015	1.0868867	3063.57940	5.4279354	1982-1999 –> 2000-2009	1.7706597	5082.4815
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	1015	0.9198056	2724.75148	3.9962002	1982-1999 –> 2000-2009	1.6418127	4848.6627
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	1015	0.9198056	2724.75148	3.9962002	1982-1999 –> 2000-2009	1.6418127	4848.6627
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	1015	1.1326704	3147.33869	5.6481941	1982-1999 –> 2000-2009	1.8550917	5272.3568
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	1015	1.1326704	3147.33869	5.6481941	1982-1999 –> 2000-2009	1.8550917	5272.3568
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1015	1.1326704	3147.33869	5.6481941	1982-1999 –> 2000-2009	1.8550917	5272.3568
Atlantic	(27.75,27.85]	4	cstar_tref ~ sal + aou + silicate + phosphate	1061	1.3115731	3598.53180	4.8846368	2000-2009 –> 2010-2019	2.4085958	6678.9548
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1061	1.1979664	3408.27472	5.3951070	2000-2009 –> 2010-2019	2.2712256	6446.2409
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1061	1.1979664	3408.27472	5.3951070	2000-2009 –> 2010-2019	2.2712256	6446.2409
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1061	1.1979664	3408.27472	5.3951070	2000-2009 –> 2010-2019	2.2712256	6446.2409
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	1061	1.3717855	3695.77976	8.1534242	2000-2009 –> 2010-2019	2.2915911	6420.5312
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1061	1.2387874	3479.37776	5.3974073	2000-2009 –> 2010-2019	2.3256740	6542.9572
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	1061	1.3717855	3695.77976	8.1534242	2000-2009 –> 2010-2019	2.2915911	6420.5312
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	1061	1.3717855	3695.77976	8.1534242	2000-2009 –> 2010-2019	2.2915911	6420.5312
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	1061	1.2929761	3570.22829	5.5437439	2000-2009 –> 2010-2019	2.4256465	6717.5670
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1061	1.2929761	3570.22829	5.5437439	2000-2009 –> 2010-2019	2.4256465	6717.5670
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	883	1.0908703	2673.44482	12.0799454	1982-1999 –> 2000-2009	2.2098926	5845.1790
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	883	1.0908703	2673.44482	12.0799454	1982-1999 –> 2000-2009	2.2098926	5845.1790
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	883	1.0908703	2673.44482	12.0799454	1982-1999 –> 2000-2009	2.2098926	5845.1790
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	883	0.8474205	2227.46961	7.9978561	1982-1999 –> 2000-2009	1.7873954	5039.6787
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	883	0.8474205	2227.46961	7.9978561	1982-1999 –> 2000-2009	1.7873954	5039.6787
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	883	0.8474205	2227.46961	7.9978561	1982-1999 –> 2000-2009	1.7873954	5039.6787
Atlantic	(27.85,27.95]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	883	1.1646554	2789.02837	7.1272880	1982-1999 –> 2000-2009	2.2104382	5821.1866
Atlantic	(27.85,27.95]	4	cstar_tref ~ temp + oxygen + aou + phosphate_star	883	1.1774324	2806.29710	7.6770524	1982-1999 –> 2000-2009	2.2238178	5837.6429
Atlantic	(27.85,27.95]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	883	1.1646554	2789.02837	7.1272880	1982-1999 –> 2000-2009	2.2104382	5821.1866
Atlantic	(27.85,27.95]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	883	1.1646554	2789.02837	7.1272880	1982-1999 –> 2000-2009	2.2104382	5821.1866
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1040	1.6145253	3961.79737	10.0335883	2000-2009 –> 2010-2019	2.7053956	6635.2422
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1040	1.6145253	3961.79737	10.0335883	2000-2009 –> 2010-2019	2.7053956	6635.2422
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1040	1.6145253	3961.79737	10.0335883	2000-2009 –> 2010-2019	2.7053956	6635.2422
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	1040	1.4948809	3801.64901	9.1499979	2000-2009 –> 2010-2019	2.3423015	6029.1186
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1040	1.6906888	4057.67502	10.1431480	2000-2009 –> 2010-2019	2.8064457	6770.9559
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	1040	1.4948809	3801.64901	9.1499979	2000-2009 –> 2010-2019	2.3423015	6029.1186
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	1040	1.4948809	3801.64901	9.1499979	2000-2009 –> 2010-2019	2.3423015	6029.1186
Atlantic	(27.85,27.95]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1040	1.7401327	4117.63182	9.4726647	2000-2009 –> 2010-2019	2.9047881	6906.6602
Atlantic	(27.85,27.95]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	1040	1.7401327	4117.63182	9.4726647	2000-2009 –> 2010-2019	2.9047881	6906.6602
Atlantic	(27.85,27.95]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	1040	1.7401327	4117.63182	9.4726647	2000-2009 –> 2010-2019	2.9047881	6906.6602
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1041	4.1999932	5956.07263	21.0748063	1982-1999 –> 2000-2009	7.8264517	12916.7352
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1041	4.1999932	5956.07263	21.0748063	1982-1999 –> 2000-2009	7.8264517	12916.7352
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1041	4.1999932	5956.07263	21.0748063	1982-1999 –> 2000-2009	7.8264517	12916.7352
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	1041	4.1999578	5956.05510	21.0710192	1982-1999 –> 2000-2009	7.9525309	13004.4369
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1041	4.1980652	5955.11667	21.0993510	1982-1999 –> 2000-2009	7.8181564	12911.2700
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1041	4.1986396	5955.40153	21.0886701	1982-1999 –> 2000-2009	7.8181061	12911.1120
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	1041	4.1999578	5956.05510	21.0710192	1982-1999 –> 2000-2009	7.9525309	13004.4369
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	1041	4.1986339	5955.39871	21.0804178	1982-1999 –> 2000-2009	7.8181374	12911.1353
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	1041	4.1986339	5955.39871	21.0804178	1982-1999 –> 2000-2009	7.8181374	12911.1353
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1041	4.1986339	5955.39871	21.0804178	1982-1999 –> 2000-2009	7.8181374	12911.1353
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1311	3.8464312	7266.67302	19.3635120	2000-2009 –> 2010-2019	8.0464244	13222.7456
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1311	3.8464312	7266.67302	19.3635120	2000-2009 –> 2010-2019	8.0464244	13222.7456
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1311	3.8464312	7266.67302	19.3635120	2000-2009 –> 2010-2019	8.0464244	13222.7456
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	1311	3.8792380	7288.94158	20.0144303	2000-2009 –> 2010-2019	8.0791958	13244.9967
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1311	3.8623389	7277.49447	19.3127339	2000-2009 –> 2010-2019	8.0604041	13232.6111
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1311	3.8731778	7284.84227	19.2843736	2000-2009 –> 2010-2019	8.0718174	13240.2438
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	1311	3.8737398	7285.22274	19.3002607	2000-2009 –> 2010-2019	8.0723737	13240.6214
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	1311	3.8737398	7285.22274	19.3002607	2000-2009 –> 2010-2019	8.0723737	13240.6214
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1311	3.8737398	7285.22274	19.3002607	2000-2009 –> 2010-2019	8.0723737	13240.6214
Atlantic	(27.95,28.05]	4	cstar_tref ~ sal + temp + silicate + phosphate_star	1311	3.8778049	7285.97280	19.5651239	2000-2009 –> 2010-2019	8.0764679	13239.3860
Atlantic	(28.05,28.1]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	595	0.6407021	1172.75998	6.5485308	1982-1999 –> 2000-2009	1.3554845	2763.8429
Atlantic	(28.05,28.1]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	595	0.6407021	1172.75998	6.5485308	1982-1999 –> 2000-2009	1.3554845	2763.8429
Atlantic	(28.05,28.1]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	595	0.6407021	1172.75998	6.5485308	1982-1999 –> 2000-2009	1.3554845	2763.8429
Atlantic	(28.05,28.1]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	595	0.6676843	1219.84843	6.2086559	1982-1999 –> 2000-2009	1.3936148	2831.4647
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	595	0.5327088	953.09832	5.9087175	1982-1999 –> 2000-2009	1.1960354	2435.4026
Atlantic	(28.05,28.1]	4	cstar_tref ~ temp + oxygen + aou + phosphate	595	0.6676843	1219.84843	6.2086559	1982-1999 –> 2000-2009	1.3936148	2831.4647
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + phosphate + phosphate_star	595	0.6676843	1219.84843	6.2086559	1982-1999 –> 2000-2009	1.3936148	2831.4647
Atlantic	(28.05,28.1]	4	cstar_tref ~ temp + oxygen + aou + phosphate_star	595	0.6676843	1219.84843	6.2086559	1982-1999 –> 2000-2009	1.3936148	2831.4647
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	595	0.5327088	953.09832	5.9087175	1982-1999 –> 2000-2009	1.1960354	2435.4026
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	595	0.5327088	953.09832	5.9087175	1982-1999 –> 2000-2009	1.1960354	2435.4026
Atlantic	(28.05,28.1]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	790	1.1664733	2499.21910	7.0545137	2000-2009 –> 2010-2019	1.8071754	3671.9791
Atlantic	(28.05,28.1]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	790	1.1664733	2499.21910	7.0545137	2000-2009 –> 2010-2019	1.8071754	3671.9791
Atlantic	(28.05,28.1]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	790	1.1664733	2499.21910	7.0545137	2000-2009 –> 2010-2019	1.8071754	3671.9791
Atlantic	(28.05,28.1]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	790	1.1699917	2501.97756	7.1949879	2000-2009 –> 2010-2019	1.8376760	3721.8260
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	790	1.0316516	2305.15734	6.1391174	2000-2009 –> 2010-2019	1.5643605	3258.2557
Atlantic	(28.05,28.1]	4	cstar_tref ~ temp + oxygen + aou + phosphate	790	1.1699917	2501.97756	7.1949879	2000-2009 –> 2010-2019	1.8376760	3721.8260
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + phosphate + phosphate_star	790	1.1699917	2501.97756	7.1949879	2000-2009 –> 2010-2019	1.8376760	3721.8260
Atlantic	(28.05,28.1]	4	cstar_tref ~ temp + oxygen + aou + phosphate_star	790	1.1699917	2501.97756	7.1949879	2000-2009 –> 2010-2019	1.8376760	3721.8260
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	790	1.0316516	2305.15734	6.1391174	2000-2009 –> 2010-2019	1.5643605	3258.2557
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	790	1.0316516	2305.15734	6.1391174	2000-2009 –> 2010-2019	1.5643605	3258.2557
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	959	0.8327978	2384.59845	7.5542203	1982-1999 –> 2000-2009	1.4714411	4048.4983
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	959	0.8327978	2384.59845	7.5542203	1982-1999 –> 2000-2009	1.4714411	4048.4983
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	959	0.8327978	2384.59845	7.5542203	1982-1999 –> 2000-2009	1.4714411	4048.4983
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	959	0.8147213	2342.50834	7.6467043	1982-1999 –> 2000-2009	1.4979936	4121.2387
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	959	0.8552215	2435.55886	7.5579020	1982-1999 –> 2000-2009	1.5031691	4124.0471
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	959	0.8147213	2342.50834	7.6467043	1982-1999 –> 2000-2009	1.4979936	4121.2387
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	959	0.8147213	2342.50834	7.6467043	1982-1999 –> 2000-2009	1.4979936	4121.2387
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	959	0.6987635	2048.03054	6.8521937	1982-1999 –> 2000-2009	1.2189835	3363.2702
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	959	0.6987635	2048.03054	6.8521937	1982-1999 –> 2000-2009	1.2189835	3363.2702
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	959	0.6987635	2048.03054	6.8521937	1982-1999 –> 2000-2009	1.2189835	3363.2702
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	1061	0.9531402	2923.14589	4.9392201	2000-2009 –> 2010-2019	1.7678615	5265.6542
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	1061	0.9531402	2923.14589	4.9392201	2000-2009 –> 2010-2019	1.7678615	5265.6542
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	1061	0.9531402	2923.14589	4.9392201	2000-2009 –> 2010-2019	1.7678615	5265.6542
Atlantic	(28.1,28.15]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	1061	0.9906195	3002.98816	5.5794292	2000-2009 –> 2010-2019	1.8141263	5364.0684
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1061	0.8048415	2564.28035	3.9867737	2000-2009 –> 2010-2019	1.5036050	4612.3109
Atlantic	(28.1,28.15]	4	cstar_tref ~ temp + oxygen + aou + phosphate	1061	0.9906195	3002.98816	5.5794292	2000-2009 –> 2010-2019	1.8141263	5364.0684
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + oxygen + aou + phosphate + phosphate_star	1061	0.9906195	3002.98816	5.5794292	2000-2009 –> 2010-2019	1.8141263	5364.0684
Atlantic	(28.1,28.15]	4	cstar_tref ~ temp + oxygen + aou + phosphate_star	1061	0.9906195	3002.98816	5.5794292	2000-2009 –> 2010-2019	1.8141263	5364.0684
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	1061	0.8048415	2564.28035	3.9867737	2000-2009 –> 2010-2019	1.5036050	4612.3109
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	1061	0.8048415	2564.28035	3.9867737	2000-2009 –> 2010-2019	1.5036050	4612.3109
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1667	0.8734430	4293.60966	9.2162370	1982-1999 –> 2000-2009	1.7251067	8097.8350
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1667	0.8734430	4293.60966	9.2162370	1982-1999 –> 2000-2009	1.7251067	8097.8350
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1667	0.8734430	4293.60966	9.2162370	1982-1999 –> 2000-2009	1.7251067	8097.8350
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	1667	0.7221335	3659.37329	18.2088194	1982-1999 –> 2000-2009	1.5482754	7371.9580
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1667	0.8959831	4378.55550	9.0810710	1982-1999 –> 2000-2009	1.7606399	8228.3853
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	1667	0.7221335	3659.37329	18.2088194	1982-1999 –> 2000-2009	1.5482754	7371.9580
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	1667	0.7221335	3659.37329	18.2088194	1982-1999 –> 2000-2009	1.5482754	7371.9580
Atlantic	(28.15,28.2]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1667	0.7651903	3852.46004	11.6924731	1982-1999 –> 2000-2009	1.6044413	7612.4639
Atlantic	(28.15,28.2]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	1667	0.7651903	3852.46004	11.6924731	1982-1999 –> 2000-2009	1.6044413	7612.4639
Atlantic	(28.15,28.2]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	1667	0.7651903	3852.46004	11.6924731	1982-1999 –> 2000-2009	1.6044413	7612.4639
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2238	0.8687676	5735.48687	5.2450626	2000-2009 –> 2010-2019	1.7422106	10029.0965
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	2238	0.8687676	5735.48687	5.2450626	2000-2009 –> 2010-2019	1.7422106	10029.0965
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	2238	0.8687676	5735.48687	5.2450626	2000-2009 –> 2010-2019	1.7422106	10029.0965
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2238	0.6019360	4093.13293	2.8170056	2000-2009 –> 2010-2019	1.3240696	7752.5062
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	2238	0.8911532	5849.35942	5.1983598	2000-2009 –> 2010-2019	1.7871363	10227.9149
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	2238	0.6019360	4093.13293	2.8170056	2000-2009 –> 2010-2019	1.3240696	7752.5062
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	2238	0.6019360	4093.13293	2.8170056	2000-2009 –> 2010-2019	1.3240696	7752.5062
Atlantic	(28.15,28.2]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	2238	0.7454151	5050.05712	3.4737189	2000-2009 –> 2010-2019	1.5106053	8902.5172
Atlantic	(28.15,28.2]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	2238	0.7454151	5050.05712	3.4737189	2000-2009 –> 2010-2019	1.5106053	8902.5172
Atlantic	(28.15,28.2]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	2238	0.7454151	5050.05712	3.4737189	2000-2009 –> 2010-2019	1.5106053	8902.5172
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2311	0.8851602	6008.51191	4.6384145	1982-1999 –> 2000-2009	1.7423575	12273.4070
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	2311	0.8851602	6008.51191	4.6384145	1982-1999 –> 2000-2009	1.7423575	12273.4070
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	2311	0.8851602	6008.51191	4.6384145	1982-1999 –> 2000-2009	1.7423575	12273.4070
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	2311	0.8862703	6014.30480	4.6236874	1982-1999 –> 2000-2009	1.7434841	12279.2950
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	2311	0.8892097	6029.60871	4.5903959	1982-1999 –> 2000-2009	1.7467003	12296.1944
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	2311	0.8892097	6029.60871	4.5903959	1982-1999 –> 2000-2009	1.7467003	12296.1944
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	2311	0.8892097	6029.60871	4.5903959	1982-1999 –> 2000-2009	1.7467003	12296.1944
Atlantic	(28.2, Inf]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	2311	0.8650146	5902.10277	5.9710497	1982-1999 –> 2000-2009	1.7452142	12297.8653
Atlantic	(28.2, Inf]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	2311	0.8650146	5902.10277	5.9710497	1982-1999 –> 2000-2009	1.7452142	12297.8653
Atlantic	(28.2, Inf]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	2311	0.8650146	5902.10277	5.9710497	1982-1999 –> 2000-2009	1.7452142	12297.8653
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	4366	0.6795930	9031.33911	4.7237651	2000-2009 –> 2010-2019	1.6238272	15338.4578
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	4366	0.6795930	9031.33911	4.7237651	2000-2009 –> 2010-2019	1.6238272	15338.4578
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	4366	0.6795930	9031.33911	4.7237651	2000-2009 –> 2010-2019	1.6238272	15338.4578
Atlantic	(28.2, Inf]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	4366	0.7098172	9409.29796	4.2065339	2000-2009 –> 2010-2019	1.6597042	15742.0043
Atlantic	(28.2, Inf]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	4366	0.7095127	9407.55125	3.9551698	2000-2009 –> 2010-2019	1.5745273	15309.6540
Atlantic	(28.2, Inf]	4	cstar_tref ~ temp + oxygen + aou + phosphate	4366	0.7098172	9409.29796	4.2065339	2000-2009 –> 2010-2019	1.6597042	15742.0043
Atlantic	(28.2, Inf]	5	cstar_tref ~ temp + oxygen + aou + phosphate + phosphate_star	4366	0.7098172	9409.29796	4.2065339	2000-2009 –> 2010-2019	1.6597042	15742.0043
Atlantic	(28.2, Inf]	4	cstar_tref ~ temp + oxygen + aou + phosphate_star	4366	0.7098172	9409.29796	4.2065339	2000-2009 –> 2010-2019	1.6597042	15742.0043
Atlantic	(28.2, Inf]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	4366	0.7095127	9407.55125	3.9551698	2000-2009 –> 2010-2019	1.5745273	15309.6540
Atlantic	(28.2, Inf]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	4366	0.7095127	9407.55125	3.9551698	2000-2009 –> 2010-2019	1.5745273	15309.6540
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1872	5.6822839	11831.15640	22.8801608	1982-1999 –> 2000-2009	10.8872183	18908.9484
Indo-Pacific	(-Inf,26]	4	cstar_tref ~ sal + oxygen + aou + phosphate	1872	5.9373577	11993.55905	24.3470956	1982-1999 –> 2000-2009	11.2671157	19123.9056
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + oxygen + aou + phosphate + phosphate_star	1872	5.9373577	11993.55905	24.3470956	1982-1999 –> 2000-2009	11.2671157	19123.9056
Indo-Pacific	(-Inf,26]	4	cstar_tref ~ sal + oxygen + aou + phosphate_star	1872	5.9373577	11993.55905	24.3470956	1982-1999 –> 2000-2009	11.2671157	19123.9056
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1872	5.6822839	11831.15640	22.8801608	1982-1999 –> 2000-2009	10.8872183	18908.9484
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1872	5.6822839	11831.15640	22.8801608	1982-1999 –> 2000-2009	10.8872183	18908.9484
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	1872	5.9122176	11979.67244	25.2385828	1982-1999 –> 2000-2009	10.9729889	18992.8055
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1872	6.0106677	12041.50403	22.8742380	1982-1999 –> 2000-2009	10.9917060	19018.0801
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	1872	5.9122176	11979.67244	25.2385828	1982-1999 –> 2000-2009	10.9729889	18992.8055
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	1872	5.9122176	11979.67244	25.2385828	1982-1999 –> 2000-2009	10.9729889	18992.8055
Indo-Pacific	(-Inf,26]	4	cstar_tref ~ sal + aou + phosphate + phosphate_star	2893	5.4886471	18073.69522	23.3309144	2000-2009 –> 2010-2019	11.4260048	30067.2543
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2893	5.1249717	17679.02558	23.0380319	2000-2009 –> 2010-2019	10.8072555	29510.1820
Indo-Pacific	(-Inf,26]	4	cstar_tref ~ sal + oxygen + aou + phosphate	2893	5.4886471	18073.69522	23.3309144	2000-2009 –> 2010-2019	11.4260048	30067.2543
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + oxygen + aou + phosphate + phosphate_star	2893	5.4886471	18073.69522	23.3309144	2000-2009 –> 2010-2019	11.4260048	30067.2543
Indo-Pacific	(-Inf,26]	4	cstar_tref ~ sal + oxygen + aou + phosphate_star	2893	5.4886471	18073.69522	23.3309144	2000-2009 –> 2010-2019	11.4260048	30067.2543
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	2893	5.1249717	17679.02558	23.0380319	2000-2009 –> 2010-2019	10.8072555	29510.1820
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	2893	5.1249717	17679.02558	23.0380319	2000-2009 –> 2010-2019	10.8072555	29510.1820
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2893	5.3226609	17898.01610	24.4672557	2000-2009 –> 2010-2019	11.2348784	29877.6885
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	2893	5.3226609	17898.01610	24.4672557	2000-2009 –> 2010-2019	11.2348784	29877.6885
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	2893	5.3226609	17898.01610	24.4672557	2000-2009 –> 2010-2019	11.2348784	29877.6885
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2139	4.1702624	12193.11301	15.3930835	1982-1999 –> 2000-2009	7.9514575	21185.2773
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	2139	4.1702624	12193.11301	15.3930835	1982-1999 –> 2000-2009	7.9514575	21185.2773
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	2139	4.1702624	12193.11301	15.3930835	1982-1999 –> 2000-2009	7.9514575	21185.2773
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2139	4.7771447	12774.33952	18.6855989	1982-1999 –> 2000-2009	8.9352076	22076.8038
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	2139	4.1791157	12202.18542	17.9668003	1982-1999 –> 2000-2009	8.0345264	21257.8323
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	2139	4.3520488	12375.64574	18.1417394	1982-1999 –> 2000-2009	8.2915588	21501.7690
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	2139	4.7771447	12774.33952	18.6855989	1982-1999 –> 2000-2009	8.9352076	22076.8038
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	2139	4.3843700	12407.29961	18.1148196	1982-1999 –> 2000-2009	8.3367729	21544.0941
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	2139	4.3843700	12407.29961	18.1148196	1982-1999 –> 2000-2009	8.3367729	21544.0941
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	2139	4.3843700	12407.29961	18.1148196	1982-1999 –> 2000-2009	8.3367729	21544.0941
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1947	3.6532273	10584.45579	20.9867736	2000-2009 –> 2010-2019	7.8234897	22777.5688
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1947	3.6532273	10584.45579	20.9867736	2000-2009 –> 2010-2019	7.8234897	22777.5688
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1947	3.6532273	10584.45579	20.9867736	2000-2009 –> 2010-2019	7.8234897	22777.5688
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1947	3.6623177	10594.13321	21.0767765	2000-2009 –> 2010-2019	7.8414334	22796.3186
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1947	3.6740110	10606.54643	20.9930871	2000-2009 –> 2010-2019	8.0260598	22982.1922
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	1947	3.6807990	10613.73430	20.9801522	2000-2009 –> 2010-2019	8.0651690	23021.0339
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	1947	3.6807990	10613.73430	20.9801522	2000-2009 –> 2010-2019	8.0651690	23021.0339
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1947	3.6807990	10613.73430	20.9801522	2000-2009 –> 2010-2019	8.0651690	23021.0339
Indo-Pacific	(26,26.5]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	1947	4.0070421	10944.42633	21.8072971	2000-2009 –> 2010-2019	8.7014939	23644.0647
Indo-Pacific	(26,26.5]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	1947	4.0070421	10944.42633	21.8072971	2000-2009 –> 2010-2019	8.7014939	23644.0647
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1901	4.2351957	10896.72338	17.3937161	1982-1999 –> 2000-2009	7.4098743	17881.5079
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1901	4.2351957	10896.72338	17.3937161	1982-1999 –> 2000-2009	7.4098743	17881.5079
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1901	4.2351957	10896.72338	17.3937161	1982-1999 –> 2000-2009	7.4098743	17881.5079
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	1901	4.5576256	11175.68470	14.7702482	1982-1999 –> 2000-2009	7.9690576	18355.2441
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1901	4.4528143	11087.22944	17.2381741	1982-1999 –> 2000-2009	7.8462385	18252.4563
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	1901	4.5576256	11175.68470	14.7702482	1982-1999 –> 2000-2009	7.9690576	18355.2441
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	1901	4.5576256	11175.68470	14.7702482	1982-1999 –> 2000-2009	7.9690576	18355.2441
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1901	4.5666130	11183.17465	18.0583552	1982-1999 –> 2000-2009	8.0054122	18384.3715
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	1901	4.5666130	11183.17465	18.0583552	1982-1999 –> 2000-2009	8.0054122	18384.3715
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	1901	4.5666130	11183.17465	18.0583552	1982-1999 –> 2000-2009	8.0054122	18384.3715
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1785	3.0303394	9037.57903	15.9688736	2000-2009 –> 2010-2019	7.2655351	19934.3024
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1785	3.0303394	9037.57903	15.9688736	2000-2009 –> 2010-2019	7.2655351	19934.3024
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1785	3.0303394	9037.57903	15.9688736	2000-2009 –> 2010-2019	7.2655351	19934.3024
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1785	3.1784717	9207.96026	17.5350490	2000-2009 –> 2010-2019	7.6312860	20295.1897
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1785	3.2232603	9257.91483	18.0029415	2000-2009 –> 2010-2019	7.7835268	20435.8019
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	1785	3.2280591	9263.22595	17.9863256	2000-2009 –> 2010-2019	7.7953185	20446.9388
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1785	3.2280591	9263.22595	17.9863256	2000-2009 –> 2010-2019	7.7953185	20446.9388
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1785	2.9887080	8988.19359	17.6209674	2000-2009 –> 2010-2019	7.5553210	20171.3682
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	1785	2.9887080	8988.19359	17.6209674	2000-2009 –> 2010-2019	7.5553210	20171.3682
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	1785	2.9887080	8988.19359	17.6209674	2000-2009 –> 2010-2019	7.5553210	20171.3682
Indo-Pacific	(26.75,27]	4	cstar_tref ~ sal + aou + phosphate + phosphate_star	2866	4.2437161	16430.61384	13.8117555	1982-1999 –> 2000-2009	7.4474419	26574.1321
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2866	4.2203297	16400.93839	14.7620525	1982-1999 –> 2000-2009	7.3794116	26491.4193
Indo-Pacific	(26.75,27]	4	cstar_tref ~ sal + oxygen + aou + phosphate	2866	4.2437161	16430.61384	13.8117555	1982-1999 –> 2000-2009	7.4474419	26574.1321
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + oxygen + aou + phosphate + phosphate_star	2866	4.2437161	16430.61384	13.8117555	1982-1999 –> 2000-2009	7.4474419	26574.1321
Indo-Pacific	(26.75,27]	4	cstar_tref ~ sal + oxygen + aou + phosphate_star	2866	4.2437161	16430.61384	13.8117555	1982-1999 –> 2000-2009	7.4474419	26574.1321
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	2866	4.2203297	16400.93839	14.7620525	1982-1999 –> 2000-2009	7.3794116	26491.4193
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	2866	4.2203297	16400.93839	14.7620525	1982-1999 –> 2000-2009	7.3794116	26491.4193
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2866	4.2380743	16424.98847	15.5224483	1982-1999 –> 2000-2009	7.3779392	26491.5384
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	2866	4.2380743	16424.98847	15.5224483	1982-1999 –> 2000-2009	7.3779392	26491.5384
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	2866	4.2380743	16424.98847	15.5224483	1982-1999 –> 2000-2009	7.3779392	26491.5384
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2275	3.5556056	12241.96088	12.2803939	2000-2009 –> 2010-2019	7.7759353	28642.8993
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	2275	3.5556056	12241.96088	12.2803939	2000-2009 –> 2010-2019	7.7759353	28642.8993
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	2275	3.5556056	12241.96088	12.2803939	2000-2009 –> 2010-2019	7.7759353	28642.8993
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2275	3.6557549	12368.34726	15.0871224	2000-2009 –> 2010-2019	7.8938293	28793.3357
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	2275	3.6534193	12365.43935	13.6360052	2000-2009 –> 2010-2019	8.0072254	28944.8561
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	2275	3.6557549	12368.34726	15.0871224	2000-2009 –> 2010-2019	7.8938293	28793.3357
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	2275	3.6557549	12368.34726	15.0871224	2000-2009 –> 2010-2019	7.8938293	28793.3357
Indo-Pacific	(26.75,27]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	2275	3.5292919	12208.16289	14.2256050	2000-2009 –> 2010-2019	7.9128368	28826.5990
Indo-Pacific	(26.75,27]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	2275	3.5292919	12208.16289	14.2256050	2000-2009 –> 2010-2019	7.9128368	28826.5990
Indo-Pacific	(26.75,27]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	2275	3.5292919	12208.16289	14.2256050	2000-2009 –> 2010-2019	7.9128368	28826.5990
Indo-Pacific	(27,27.25]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	3705	3.6826305	20188.21280	11.9756687	1982-1999 –> 2000-2009	6.6055492	31653.2443
Indo-Pacific	(27,27.25]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	3705	3.6826305	20188.21280	11.9756687	1982-1999 –> 2000-2009	6.6055492	31653.2443
Indo-Pacific	(27,27.25]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	3705	3.6826305	20188.21280	11.9756687	1982-1999 –> 2000-2009	6.6055492	31653.2443
Indo-Pacific	(27,27.25]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	3705	3.7519112	20324.32091	12.4156537	1982-1999 –> 2000-2009	6.7286779	31871.2525
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	3705	3.5173624	19847.97597	12.5447011	1982-1999 –> 2000-2009	6.4689804	31357.9141
Indo-Pacific	(27,27.25]	4	cstar_tref ~ temp + oxygen + aou + phosphate	3705	3.7519112	20324.32091	12.4156537	1982-1999 –> 2000-2009	6.7286779	31871.2525
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + oxygen + aou + phosphate + phosphate_star	3705	3.7519112	20324.32091	12.4156537	1982-1999 –> 2000-2009	6.7286779	31871.2525
Indo-Pacific	(27,27.25]	4	cstar_tref ~ temp + oxygen + aou + phosphate_star	3705	3.7519112	20324.32091	12.4156537	1982-1999 –> 2000-2009	6.7286779	31871.2525
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	3705	3.5173624	19847.97597	12.5447011	1982-1999 –> 2000-2009	6.4689804	31357.9141
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	3705	3.5173624	19847.97597	12.5447011	1982-1999 –> 2000-2009	6.4689804	31357.9141
Indo-Pacific	(27,27.25]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	3299	4.0577959	18617.57877	19.8533642	2000-2009 –> 2010-2019	7.7404263	38805.7916
Indo-Pacific	(27,27.25]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	3299	4.0577959	18617.57877	19.8533642	2000-2009 –> 2010-2019	7.7404263	38805.7916
Indo-Pacific	(27,27.25]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	3299	4.0577959	18617.57877	19.8533642	2000-2009 –> 2010-2019	7.7404263	38805.7916
Indo-Pacific	(27,27.25]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	3299	4.1465717	18758.37243	18.4409321	2000-2009 –> 2010-2019	7.8984829	39082.6933
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	3299	3.7181208	18040.77134	16.0517451	2000-2009 –> 2010-2019	7.2354832	37888.7473
Indo-Pacific	(27,27.25]	4	cstar_tref ~ temp + oxygen + aou + phosphate	3299	4.1465717	18758.37243	18.4409321	2000-2009 –> 2010-2019	7.8984829	39082.6933
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + oxygen + aou + phosphate + phosphate_star	3299	4.1465717	18758.37243	18.4409321	2000-2009 –> 2010-2019	7.8984829	39082.6933
Indo-Pacific	(27,27.25]	4	cstar_tref ~ temp + oxygen + aou + phosphate_star	3299	4.1465717	18758.37243	18.4409321	2000-2009 –> 2010-2019	7.8984829	39082.6933
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	3299	3.7181208	18040.77134	16.0517451	2000-2009 –> 2010-2019	7.2354832	37888.7473
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	3299	3.7181208	18040.77134	16.0517451	2000-2009 –> 2010-2019	7.2354832	37888.7473
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	4717	2.9771414	23692.41654	11.0261154	1982-1999 –> 2000-2009	5.3593669	37620.3642
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	4717	2.9771414	23692.41654	11.0261154	1982-1999 –> 2000-2009	5.3593669	37620.3642
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	4717	2.9771414	23692.41654	11.0261154	1982-1999 –> 2000-2009	5.3593669	37620.3642
Indo-Pacific	(27.25,27.5]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	4717	2.9797311	23698.61932	10.9030999	1982-1999 –> 2000-2009	5.3706774	37646.7984
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	4717	2.7781045	23039.63341	9.6159708	1982-1999 –> 2000-2009	5.0244379	36610.2324
Indo-Pacific	(27.25,27.5]	4	cstar_tref ~ temp + oxygen + aou + phosphate	4717	2.9797311	23698.61932	10.9030999	1982-1999 –> 2000-2009	5.3706774	37646.7984
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ temp + oxygen + aou + phosphate + phosphate_star	4717	2.9797311	23698.61932	10.9030999	1982-1999 –> 2000-2009	5.3706774	37646.7984
Indo-Pacific	(27.25,27.5]	4	cstar_tref ~ temp + oxygen + aou + phosphate_star	4717	2.9797311	23698.61932	10.9030999	1982-1999 –> 2000-2009	5.3706774	37646.7984
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	4717	2.7781045	23039.63341	9.6159708	1982-1999 –> 2000-2009	5.0244379	36610.2324
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	4717	2.7781045	23039.63341	9.6159708	1982-1999 –> 2000-2009	5.0244379	36610.2324
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	4380	3.3158089	22944.52774	15.1365225	2000-2009 –> 2010-2019	6.8118097	48152.5739
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	4380	3.2683330	22818.19500	14.7103876	2000-2009 –> 2010-2019	6.7502185	47988.0737
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + aou + silicate	4380	2.8740352	21691.98274	11.6788375	2000-2009 –> 2010-2019	5.9726227	45761.5969
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	4380	3.2720282	22828.09338	14.8560785	2000-2009 –> 2010-2019	6.7870846	48087.4219
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	4380	3.2720282	22828.09338	14.8560785	2000-2009 –> 2010-2019	6.7870846	48087.4219
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	4380	3.2720282	22828.09338	14.8560785	2000-2009 –> 2010-2019	6.7870846	48087.4219
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	4380	3.2608540	22798.12623	14.3919287	2000-2009 –> 2010-2019	6.0389584	45837.7596
Indo-Pacific	(27.25,27.5]	4	cstar_tref ~ temp + oxygen + aou + silicate	4380	3.2969240	22892.49320	15.3937428	2000-2009 –> 2010-2019	6.4210326	47037.4908
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	4380	3.2608540	22798.12623	14.3919287	2000-2009 –> 2010-2019	6.0389584	45837.7596
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	4380	3.2608540	22798.12623	14.3919287	2000-2009 –> 2010-2019	6.0389584	45837.7596
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	5654	3.0263866	28581.48987	19.0026905	1982-1999 –> 2000-2009	5.2316982	43524.9483
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	5654	3.0059510	28504.87389	9.5282956	1982-1999 –> 2000-2009	5.3075226	43736.9714
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	5654	3.0263866	28581.48987	19.0026905	1982-1999 –> 2000-2009	5.2316982	43524.9483
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	5654	3.0263866	28581.48987	19.0026905	1982-1999 –> 2000-2009	5.2316982	43524.9483
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + aou + silicate	5654	2.8521390	27910.92404	19.1656654	1982-1999 –> 2000-2009	5.1576728	43154.6420
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	5654	3.0135913	28533.57915	9.6219707	1982-1999 –> 2000-2009	5.3375570	43831.0943
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	5654	3.0135913	28533.57915	9.6219707	1982-1999 –> 2000-2009	5.3375570	43831.0943
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	5654	2.8773091	28010.27904	24.8611792	1982-1999 –> 2000-2009	4.9751708	42616.2738
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	5654	2.8773091	28010.27904	24.8611792	1982-1999 –> 2000-2009	4.9751708	42616.2738
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	5654	2.8773091	28010.27904	24.8611792	1982-1999 –> 2000-2009	4.9751708	42616.2738
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	4866	2.6254661	23217.00496	13.0172487	2000-2009 –> 2010-2019	5.6587725	51824.3209
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	4866	2.6254661	23217.00496	13.0172487	2000-2009 –> 2010-2019	5.6587725	51824.3209
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	4866	2.6254661	23217.00496	13.0172487	2000-2009 –> 2010-2019	5.6587725	51824.3209
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	4866	2.6169790	23185.49419	12.5242942	2000-2009 –> 2010-2019	5.6229300	51690.3681
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	4866	2.6175753	23187.71135	12.5830541	2000-2009 –> 2010-2019	5.6390648	51750.8887
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + aou + silicate	4866	2.6406992	23273.30723	12.4994113	2000-2009 –> 2010-2019	5.4928382	51184.2313
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	4866	2.6161546	23182.42798	12.4685846	2000-2009 –> 2010-2019	5.6297459	51716.0071
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	4866	2.6161546	23182.42798	12.4685846	2000-2009 –> 2010-2019	5.6297459	51716.0071
Indo-Pacific	(27.5,27.75]	4	cstar_tref ~ sal + temp + silicate + phosphate	4866	2.8069056	23865.33815	8.5937438	2000-2009 –> 2010-2019	5.8541589	52522.5283
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	4866	2.6161546	23182.42798	12.4685846	2000-2009 –> 2010-2019	5.6297459	51716.0071
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2196	2.2008118	9710.50312	10.4253247	1982-1999 –> 2000-2009	3.8172501	14920.6149
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	2196	2.2008118	9710.50312	10.4253247	1982-1999 –> 2000-2009	3.8172501	14920.6149
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	2196	2.2008118	9710.50312	10.4253247	1982-1999 –> 2000-2009	3.8172501	14920.6149
Indo-Pacific	(27.75,27.85]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	2196	2.3244657	9948.58704	11.8291448	1982-1999 –> 2000-2009	3.9422031	15158.8968
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	2196	2.3202416	9942.59848	11.6925947	1982-1999 –> 2000-2009	3.8562712	15013.1080
Indo-Pacific	(27.75,27.85]	4	cstar_tref ~ temp + oxygen + aou + phosphate	2196	2.3244657	9948.58704	11.8291448	1982-1999 –> 2000-2009	3.9422031	15158.8968
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ temp + oxygen + aou + phosphate + phosphate_star	2196	2.3244657	9948.58704	11.8291448	1982-1999 –> 2000-2009	3.9422031	15158.8968
Indo-Pacific	(27.75,27.85]	4	cstar_tref ~ temp + oxygen + aou + phosphate_star	2196	2.3244657	9948.58704	11.8291448	1982-1999 –> 2000-2009	3.9422031	15158.8968
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	2196	2.3202416	9942.59848	11.6925947	1982-1999 –> 2000-2009	3.8562712	15013.1080
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	2196	2.3202416	9942.59848	11.6925947	1982-1999 –> 2000-2009	3.8562712	15013.1080
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1929	2.4934912	9013.25710	10.5369790	2000-2009 –> 2010-2019	4.8283160	18983.3737
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1929	2.4934912	9013.25710	10.5369790	2000-2009 –> 2010-2019	4.8283160	18983.3737
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1929	2.4934912	9013.25710	10.5369790	2000-2009 –> 2010-2019	4.8283160	18983.3737
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1929	2.5073628	9034.66012	9.9861229	2000-2009 –> 2010-2019	4.8617211	19041.3680
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1929	2.5117660	9041.42921	9.8162742	2000-2009 –> 2010-2019	4.8733691	19061.6314
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	1929	2.8110038	9475.66847	17.5332143	2000-2009 –> 2010-2019	5.0118156	19186.1716
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	1929	2.8110038	9475.66847	17.5332143	2000-2009 –> 2010-2019	5.0118156	19186.1716
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	1929	2.5144202	9045.50382	9.7291295	2000-2009 –> 2010-2019	4.8784939	19070.2983
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	1929	2.5144202	9045.50382	9.7291295	2000-2009 –> 2010-2019	4.8784939	19070.2983
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1929	2.5144202	9045.50382	9.7291295	2000-2009 –> 2010-2019	4.8784939	19070.2983
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2633	1.9048340	10879.51369	8.6514753	1982-1999 –> 2000-2009	3.3754351	16227.9402
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	2633	1.9048340	10879.51369	8.6514753	1982-1999 –> 2000-2009	3.3754351	16227.9402
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	2633	1.9048340	10879.51369	8.6514753	1982-1999 –> 2000-2009	3.3754351	16227.9402
Indo-Pacific	(27.85,27.95]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	2633	1.9436110	10983.63790	9.5632442	1982-1999 –> 2000-2009	3.4163184	16334.2951
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	2633	1.9132780	10902.80597	9.1535103	1982-1999 –> 2000-2009	3.3474582	16177.1023
Indo-Pacific	(27.85,27.95]	4	cstar_tref ~ temp + oxygen + aou + phosphate	2633	1.9436110	10983.63790	9.5632442	1982-1999 –> 2000-2009	3.4163184	16334.2951
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ temp + oxygen + aou + phosphate + phosphate_star	2633	1.9436110	10983.63790	9.5632442	1982-1999 –> 2000-2009	3.4163184	16334.2951
Indo-Pacific	(27.85,27.95]	4	cstar_tref ~ temp + oxygen + aou + phosphate_star	2633	1.9436110	10983.63790	9.5632442	1982-1999 –> 2000-2009	3.4163184	16334.2951
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	2633	1.9132780	10902.80597	9.1535103	1982-1999 –> 2000-2009	3.3474582	16177.1023
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	2633	1.9132780	10902.80597	9.1535103	1982-1999 –> 2000-2009	3.3474582	16177.1023
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2352	3.3320991	12350.43280	12.8043059	2000-2009 –> 2010-2019	5.4091529	23685.7483
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	2352	3.3320991	12350.43280	12.8043059	2000-2009 –> 2010-2019	5.4091529	23685.7483
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	2352	3.3320991	12350.43280	12.8043059	2000-2009 –> 2010-2019	5.4091529	23685.7483
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2352	3.4602416	12527.94281	12.2401246	2000-2009 –> 2010-2019	5.3650756	23407.4565
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	2352	3.3530878	12379.97018	12.6031341	2000-2009 –> 2010-2019	5.4537076	23774.6966
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	2352	3.4602416	12527.94281	12.2401246	2000-2009 –> 2010-2019	5.3650756	23407.4565
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	2352	3.4602416	12527.94281	12.2401246	2000-2009 –> 2010-2019	5.3650756	23407.4565
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	2352	3.4418767	12502.91022	11.7228160	2000-2009 –> 2010-2019	5.3551547	23405.7162
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	2352	3.4418767	12502.91022	11.7228160	2000-2009 –> 2010-2019	5.3551547	23405.7162
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	2352	3.4418767	12502.91022	11.7228160	2000-2009 –> 2010-2019	5.3551547	23405.7162
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2926	1.7066174	11445.60009	10.2618901	1982-1999 –> 2000-2009	3.3595364	18369.2470
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	2926	1.7066174	11445.60009	10.2618901	1982-1999 –> 2000-2009	3.3595364	18369.2470
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	2926	1.7066174	11445.60009	10.2618901	1982-1999 –> 2000-2009	3.3595364	18369.2470
Indo-Pacific	(27.95,28.05]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	2926	1.7270870	11513.37277	11.5148887	1982-1999 –> 2000-2009	3.3826378	18440.7408
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	2926	1.6379271	11205.18929	9.6731446	1982-1999 –> 2000-2009	3.1654735	17845.1837
Indo-Pacific	(27.95,28.05]	4	cstar_tref ~ temp + oxygen + aou + phosphate	2926	1.7270870	11513.37277	11.5148887	1982-1999 –> 2000-2009	3.3826378	18440.7408
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ temp + oxygen + aou + phosphate + phosphate_star	2926	1.7270870	11513.37277	11.5148887	1982-1999 –> 2000-2009	3.3826378	18440.7408
Indo-Pacific	(27.95,28.05]	4	cstar_tref ~ temp + oxygen + aou + phosphate_star	2926	1.7270870	11513.37277	11.5148887	1982-1999 –> 2000-2009	3.3826378	18440.7408
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	2926	1.6379271	11205.18929	9.6731446	1982-1999 –> 2000-2009	3.1654735	17845.1837
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	2926	1.6379271	11205.18929	9.6731446	1982-1999 –> 2000-2009	3.1654735	17845.1837
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	3145	1.4736678	11378.09846	9.5879259	2000-2009 –> 2010-2019	3.2269126	22981.4389
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	3145	1.4736678	11378.09846	9.5879259	2000-2009 –> 2010-2019	3.2269126	22981.4389
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	3145	1.4736678	11378.09846	9.5879259	2000-2009 –> 2010-2019	3.2269126	22981.4389
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	3145	1.5425852	11665.58479	10.1169456	2000-2009 –> 2010-2019	3.2492026	23111.1849
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	3145	1.4739302	11379.21853	9.5485982	2000-2009 –> 2010-2019	3.2475391	23050.1386
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate	3145	1.5425852	11665.58479	10.1169456	2000-2009 –> 2010-2019	3.2492026	23111.1849
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + oxygen + aou + phosphate_star	3145	1.5425852	11665.58479	10.1169456	2000-2009 –> 2010-2019	3.2492026	23111.1849
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	3145	1.5961996	11880.48824	10.3375834	2000-2009 –> 2010-2019	3.2341267	23085.6775
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	3145	1.5961996	11880.48824	10.3375834	2000-2009 –> 2010-2019	3.2341267	23085.6775
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	3145	1.5961996	11880.48824	10.3375834	2000-2009 –> 2010-2019	3.2341267	23085.6775
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1646	1.5899131	6211.57818	7.0924822	1982-1999 –> 2000-2009	3.1343160	10366.4667
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1646	1.5899131	6211.57818	7.0924822	1982-1999 –> 2000-2009	3.1343160	10366.4667
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1646	1.5899131	6211.57818	7.0924822	1982-1999 –> 2000-2009	3.1343160	10366.4667
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1646	1.6034508	6239.48999	7.1712774	1982-1999 –> 2000-2009	3.1581753	10409.2592
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	1646	1.6070635	6246.89870	7.2017473	1982-1999 –> 2000-2009	3.1653670	10421.8047
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	1646	1.6070635	6246.89870	7.2017473	1982-1999 –> 2000-2009	3.1653670	10421.8047
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1646	1.6070635	6246.89870	7.2017473	1982-1999 –> 2000-2009	3.1653670	10421.8047
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1646	1.4954317	6009.89555	7.3348351	1982-1999 –> 2000-2009	2.9364594	10010.0107
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	1646	1.4954317	6009.89555	7.3348351	1982-1999 –> 2000-2009	2.9364594	10010.0107
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	1646	1.4954317	6009.89555	7.3348351	1982-1999 –> 2000-2009	2.9364594	10010.0107
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1803	1.3088524	6101.24979	6.5941017	2000-2009 –> 2010-2019	2.8987655	12312.8280
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	1803	1.3088524	6101.24979	6.5941017	2000-2009 –> 2010-2019	2.8987655	12312.8280
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	1803	1.3088524	6101.24979	6.5941017	2000-2009 –> 2010-2019	2.8987655	12312.8280
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1803	1.3088434	6101.22511	6.5828659	2000-2009 –> 2010-2019	2.9122942	12340.7151
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	1803	1.3088698	6101.29776	6.5795089	2000-2009 –> 2010-2019	2.9159333	12348.1965
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	1803	1.3088698	6101.29776	6.5795089	2000-2009 –> 2010-2019	2.9159333	12348.1965
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1803	1.3088698	6101.29776	6.5795089	2000-2009 –> 2010-2019	2.9159333	12348.1965
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1803	1.3859730	6307.69950	6.7452924	2000-2009 –> 2010-2019	2.8814047	12317.5950
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate	1803	1.3859730	6307.69950	6.7452924	2000-2009 –> 2010-2019	2.8814047	12317.5950
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star	1803	1.3859730	6307.69950	6.7452924	2000-2009 –> 2010-2019	2.8814047	12317.5950
Indo-Pacific	(28.1, Inf]	4	cstar_tref ~ sal + aou + silicate + phosphate	19259	1.5800945	72288.06748	9.4630717	1982-1999 –> 2000-2009	2.9959453	117275.9310
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	19259	1.5326582	71115.99851	8.7999921	1982-1999 –> 2000-2009	2.8402576	114080.9914
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	19259	1.5326582	71115.99851	8.7999921	1982-1999 –> 2000-2009	2.8402576	114080.9914
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	19259	1.5326582	71115.99851	8.7999921	1982-1999 –> 2000-2009	2.8402576	114080.9914
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	19259	1.5328856	71121.71315	8.8027424	1982-1999 –> 2000-2009	2.8396687	114070.8089
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	19259	1.5331400	71128.10654	8.8052272	1982-1999 –> 2000-2009	2.8389172	114057.5974
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	19259	1.5331227	71127.67220	8.8055135	1982-1999 –> 2000-2009	2.8389409	114057.9615
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	19259	1.5331227	71127.67220	8.8055135	1982-1999 –> 2000-2009	2.8389409	114057.9615
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	19259	1.5331227	71127.67220	8.8055135	1982-1999 –> 2000-2009	2.8389409	114057.9615
Indo-Pacific	(28.1, Inf]	4	cstar_tref ~ sal + temp + silicate + phosphate_star	19259	1.5412842	71330.17739	8.4479351	1982-1999 –> 2000-2009	2.8601533	114511.6427
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	19574	1.3883679	68408.19401	8.7391845	2000-2009 –> 2010-2019	2.9210260	139524.1925
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate	19574	1.3883679	68408.19401	8.7391845	2000-2009 –> 2010-2019	2.9210260	139524.1925
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star	19574	1.3883679	68408.19401	8.7391845	2000-2009 –> 2010-2019	2.9210260	139524.1925
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	19574	1.3885810	68414.20407	8.7393905	2000-2009 –> 2010-2019	2.9214666	139535.9172
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	19574	1.3888244	68421.06581	8.7404371	2000-2009 –> 2010-2019	2.9219644	139549.1723
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate	19574	1.3887997	68420.36924	8.7407568	2000-2009 –> 2010-2019	2.9219224	139548.0414
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star	19574	1.3887997	68420.36924	8.7407568	2000-2009 –> 2010-2019	2.9219224	139548.0414
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	19574	1.3887997	68420.36924	8.7407568	2000-2009 –> 2010-2019	2.9219224	139548.0414
Indo-Pacific	(28.1, Inf]	4	cstar_tref ~ sal + temp + silicate + phosphate_star	19574	1.4025767	68804.80513	8.7327719	2000-2009 –> 2010-2019	2.9438609	140134.9825
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	19574	1.4243319	69409.36568	8.6185167	2000-2009 –> 2010-2019	2.9951675	141473.0664

5.3.2 per fitting unit

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

basin	gamma_slab	eras	rmse_sum_mean	ais_sum_mean
Atlantic	(-Inf,26]	1982-1999 –> 2000-2009	2.1210540	263.6528
Atlantic	(-Inf,26]	2000-2009 –> 2010-2019	0.9191149	124.6920
Atlantic	(26,26.5]	1982-1999 –> 2000-2009	7.6027427	2926.3121
Atlantic	(26,26.5]	2000-2009 –> 2010-2019	7.9592667	3096.2108
Atlantic	(26.5,26.75]	1982-1999 –> 2000-2009	5.7496914	3197.3717
Atlantic	(26.5,26.75]	2000-2009 –> 2010-2019	7.4396214	3886.2545
Atlantic	(26.75,27]	1982-1999 –> 2000-2009	3.9620434	4997.4823
Atlantic	(26.75,27]	2000-2009 –> 2010-2019	5.2402633	5971.5535
Atlantic	(27,27.25]	1982-1999 –> 2000-2009	3.6978293	5905.8942
Atlantic	(27,27.25]	2000-2009 –> 2010-2019	4.4651115	7223.1128
Atlantic	(27.25,27.5]	1982-1999 –> 2000-2009	4.9299306	11239.2576
Atlantic	(27.25,27.5]	2000-2009 –> 2010-2019	6.7116679	15057.0292
Atlantic	(27.5,27.75]	1982-1999 –> 2000-2009	4.0111484	21499.9007
Atlantic	(27.5,27.75]	2000-2009 –> 2010-2019	4.6100257	23782.5407
Atlantic	(27.75,27.85]	1982-1999 –> 2000-2009	1.7460654	5040.9805
Atlantic	(27.75,27.85]	2000-2009 –> 2010-2019	2.3274013	6525.7362
Atlantic	(27.85,27.95]	1982-1999 –> 2000-2009	2.0846996	5595.5776
Atlantic	(27.85,27.95]	2000-2009 –> 2010-2019	2.6663901	6548.4019
Atlantic	(27.95,28.05]	1982-1999 –> 2000-2009	7.8475091	12931.4867
Atlantic	(27.95,28.05]	2000-2009 –> 2010-2019	8.0644280	13234.7339
Atlantic	(28.05,28.1]	1982-1999 –> 2000-2009	1.3229019	2692.3595
Atlantic	(28.05,28.1]	2000-2009 –> 2010-2019	1.7465312	3567.8008
Atlantic	(28.1,28.15]	1982-1999 –> 2000-2009	1.4068424	3872.3069
Atlantic	(28.1,28.15]	2000-2009 –> 2010-2019	1.7070905	5109.0169
Atlantic	(28.15,28.2]	1982-1999 –> 2000-2009	1.6394110	7747.5156
Atlantic	(28.15,28.2]	2000-2009 –> 2010-2019	1.5517793	9028.0275
Atlantic	(28.2, Inf]	1982-1999 –> 2000-2009	1.7446300	12288.1695
Atlantic	(28.2, Inf]	2000-2009 –> 2010-2019	1.6233881	15491.2353
Indo-Pacific	(-Inf,26]	1982-1999 –> 2000-2009	11.0373675	19009.5059
Indo-Pacific	(-Inf,26]	2000-2009 –> 2010-2019	11.1830421	29843.2629
Indo-Pacific	(26,26.5]	1982-1999 –> 2000-2009	8.3061191	21510.1323
Indo-Pacific	(26,26.5]	2000-2009 –> 2010-2019	8.0936457	23046.2448
Indo-Pacific	(26.5,26.75]	1982-1999 –> 2000-2009	7.7999271	18211.5827
Indo-Pacific	(26.5,26.75]	2000-2009 –> 2010-2019	7.5468018	20194.1881
Indo-Pacific	(26.75,27]	1982-1999 –> 2000-2009	7.4061820	26524.5401
Indo-Pacific	(26.75,27]	2000-2009 –> 2010-2019	7.8755029	28773.3358
Indo-Pacific	(27,27.25]	1982-1999 –> 2000-2009	6.6138300	31651.8485
Indo-Pacific	(27,27.25]	2000-2009 –> 2010-2019	7.6521661	38641.4390
Indo-Pacific	(27.25,27.5]	1982-1999 –> 2000-2009	5.2634124	37327.8984
Indo-Pacific	(27.25,27.5]	2000-2009 –> 2010-2019	6.4433813	47071.5280
Indo-Pacific	(27.5,27.75]	1982-1999 –> 2000-2009	5.1760916	43297.7468
Indo-Pacific	(27.5,27.75]	2000-2009 –> 2010-2019	5.6474547	51776.9000
Indo-Pacific	(27.75,27.85]	1982-1999 –> 2000-2009	3.8789376	15043.6756
Indo-Pacific	(27.75,27.85]	2000-2009 –> 2010-2019	4.8879151	19063.6359
Indo-Pacific	(27.85,27.95]	1982-1999 –> 2000-2009	3.3833954	16255.2308
Indo-Pacific	(27.85,27.95]	2000-2009 –> 2010-2019	5.3841857	23527.1460
Indo-Pacific	(27.95,28.05]	1982-1999 –> 2000-2009	3.3105581	18240.6255
Indo-Pacific	(27.95,28.05]	2000-2009 –> 2010-2019	3.2378265	23058.5042
Indo-Pacific	(28.05,28.1]	1982-1999 –> 2000-2009	3.0866602	10280.4106
Indo-Pacific	(28.05,28.1]	2000-2009 –> 2010-2019	2.9000605	12327.6574
Indo-Pacific	(28.1, Inf]	1982-1999 –> 2000-2009	2.8572280	114433.2839
Indo-Pacific	(28.1, Inf]	2000-2009 –> 2010-2019	2.9311305	139790.9840

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
f745381	Donghe-Zhu	2021-03-16
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
f745381	Donghe-Zhu	2021-03-16
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
f745381	Donghe-Zhu	2021-03-16
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
f745381	Donghe-Zhu	2021-03-16
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-1999 --> 2000-2009
(-Inf,26]	8	8	5	5	6	7	10
(26,26.5]	6	4	6	7	9	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	4	7	6	9	10	6
(27.25,27.5]	7	5	6	6	10	9	6
(27.5,27.75]	6	5	6	6	10	5	10
(27.75,27.85]	7	6	7	6	10	7	7
(27.85,27.95]	10	7	6	7	6	6	7
(27.95,28.05]	7	5	7	6	10	8	7
(28.05,28.1]	10	7	7	6	3	3	10
(28.1,28.15]	10	6	7	6	7	7	7
(28.15,28.2]	10	6	7	6	7	7	7
(28.2, Inf]	7	6	7	6	7	10	7
total	113.00	81.00	91.00	86.00	104.00	109.00	101.00
Atlantic - 2000-2009 --> 2010-2019
(-Inf,26]	7	5	7	6	7	10	8
(26,26.5]	7	6	6	5	9	10	7
(26.5,26.75]	8	5	7	6	7	10	7
(26.75,27]	10	6	8	5	4	10	4
(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]	7	4	7	6	10	9	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]	6	4	6	7	10	9	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	116.00	76.00	95.00	83.00	100.00	108.00	100.00
Indo-Pacific - 1982-1999 --> 2000-2009
(-Inf,26]	10	7	7	5	10	4	4
(26,26.5]	7	5	7	6	10	8	7
(26.5,26.75]	10	6	7	6	7	7	7
(26.75,27]	10	7	7	6	10	3	3
(27,27.25]	10	7	7	6	3	3	10
(27.25,27.5]	10	7	7	6	3	3	10
(27.5,27.75]	8	6	6	6	7	7	10
(27.75,27.85]	10	7	7	6	3	3	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	108.00	76.00	82.00	71.00	76.00	64.00	94.00
Indo-Pacific - 2000-2009 --> 2010-2019
(-Inf,26]	10	7	7	6	10	3	3
(26,26.5]	7	6	6	6	8	10	7
(26.5,26.75]	8	5	6	7	7	10	7
(26.75,27]	10	6	7	6	7	7	7
(27,27.25]	10	7	7	6	3	3	10
(27.25,27.5]	7	6	6	5	6	10	9
(27.5,27.75]	6	5	6	5	10	10	7
(27.75,27.85]	7	6	6	6	10	8	7
(27.85,27.95]	10	6	7	6	7	7	7
(27.95,28.05]	10	6	7	6	7	7	7
(28.05,28.1]	7	6	7	6	7	10	7
(28.1, Inf]	6	4	6	7	9	10	7
total	98.00	70.00	78.00	72.00	91.00	95.00	85.00

6.4 RMSE alternatives

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

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

Past versions of aic_vs_rmse_all_models-1.png

Version	Author	Date
f745381	Donghe-Zhu	2021-03-16
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
f745381	Donghe-Zhu	2021-03-16
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