eMLR - model fitting

workflowr

• Summary
• Checks
• Past versions

Last updated: 2021-03-27

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: 88206bb

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 88206bb. 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/

Untracked files:
    Untracked:  2_Pac_zon_em.png
    Untracked:  so_Pac_zon_off.png

Unstaged changes:
    Deleted:    README.md
    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	b883157	Donghe-Zhu	2021-03-27	Build site.
html	0c20513	Donghe-Zhu	2021-03-26	Build site.
html	d19654d	Donghe-Zhu	2021-03-26	Build site.
Rmd	921c9db	Donghe-Zhu	2021-03-26	local rebuild
html	6c53dbf	Donghe-Zhu	2021-03-25	Build site.
Rmd	195db95	Donghe-Zhu	2021-03-24	local rebuild
html	1914a11	Donghe-Zhu	2021-03-24	Build site.
Rmd	e32af14	Donghe-Zhu	2021-03-24	local rebuild
html	8be810e	Donghe-Zhu	2021-03-23	Build site.
Rmd	85274db	Donghe-Zhu	2021-03-23	local rebuild
html	bf19764	Donghe-Zhu	2021-03-22	Build site.
html	3ec9d3d	Donghe-Zhu	2021-03-22	Build site.
Rmd	073fa6a	Donghe-Zhu	2021-03-22	local rebuild
html	134ace1	Donghe-Zhu	2021-03-22	Build site.
html	f6d70a4	Donghe-Zhu	2021-03-22	Build site.
html	183443b	Donghe-Zhu	2021-03-21	Build site.
Rmd	adca794	Donghe-Zhu	2021-03-21	local rebuild
html	2e6976b	Donghe-Zhu	2021-03-21	Build site.
Rmd	4e7de35	Donghe-Zhu	2021-03-21	local rebuild
html	51a42bd	Donghe-Zhu	2021-03-16	Build site.
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, aou, nitrate, silicate, phosphate, phosphate_star

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

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

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

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)[100]
        
        # 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)
          
          # calculate max VIF
          if (is.na(max(i_lm_fit$coefficients))){
            i_vif_max <- as.double("NA")
          } else{
          i_vif_max <- max(vif(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,
              vif_max = i_vif_max,
              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_vif_max,
  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,
    vif_max,
    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)

5.3 VIF threshold

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

• 10

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

# remove models with predictors fitted as NA
lm_all_fitted_wide <- lm_all_fitted_wide %>%
  filter(vif_max <= params_local$vif_max)

lm_all_fitted_wide_check <- lm_all_fitted_wide %>%
  group_by(era, basin, gamma_slab) %>% 
  count() %>% 
  filter(n < params_local$MLR_number)

# 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() %>%
  filter(eras != "1982-1999 --> 2010-2019")

# 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.4 RMSE tables

5.4.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 + aou + silicate + phosphate + phosphate_star	26	0.4612966	47.551668	1.1237758	1982-1999 –> 2000-2009	2.4644429	281.37206
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	26	0.6987083	69.141666	2.2173182	1982-1999 –> 2000-2009	2.2138763	273.92601
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	26	0.4589059	47.281475	1.1203109	1982-1999 –> 2000-2009	2.3155101	273.20099
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	26	0.4584876	47.234056	1.1189189	1982-1999 –> 2000-2009	2.2947015	272.00507
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + nitrate + silicate + phosphate	26	0.4678743	48.287914	1.5811594	1982-1999 –> 2000-2009	2.4153823	279.17878
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	26	0.4450286	45.684733	1.0682985	1982-1999 –> 2000-2009	2.2784176	270.29562
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	26	0.4584940	47.234783	1.1190183	1982-1999 –> 2000-2009	2.2860715	271.51549
Atlantic	(-Inf,26]	5	cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star	26	0.6974382	69.047053	2.0723505	1982-1999 –> 2000-2009	2.1344991	268.32706
Atlantic	(-Inf,26]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	26	0.7086122	67.873568	2.0369830	1982-1999 –> 2000-2009	2.3733810	280.45054
Atlantic	(-Inf,26]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	26	0.5133095	53.107234	0.9152479	1982-1999 –> 2000-2009	1.8797173	247.14411
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	26	0.2946619	24.244614	0.7713187	2000-2009 –> 2010-2019	0.7559584	71.79628
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + nitrate + silicate	26	0.2506018	15.822521	0.4743830	2000-2009 –> 2010-2019	0.7519133	67.69989
Atlantic	(-Inf,26]	4	cstar_tref ~ sal + temp + aou + silicate	26	0.2833559	20.210130	0.7252885	2000-2009 –> 2010-2019	0.8169883	73.33643
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	26	0.2790915	21.421584	0.7403142	2000-2009 –> 2010-2019	0.7379973	68.70306
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	26	0.2782978	21.273503	0.7394145	2000-2009 –> 2010-2019	0.7367854	68.50756
Atlantic	(-Inf,26]	4	cstar_tref ~ sal + temp + nitrate + silicate	26	0.2743139	18.523733	0.7065044	2000-2009 –> 2010-2019	0.7756628	68.40498
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + nitrate + silicate + phosphate	26	0.2727615	20.228613	0.6440915	2000-2009 –> 2010-2019	0.7406358	68.51653
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	26	0.2609157	17.919799	0.6909568	2000-2009 –> 2010-2019	0.7059443	63.60453
Atlantic	(-Inf,26]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	26	0.2749403	20.642332	0.7285522	2000-2009 –> 2010-2019	0.7334343	67.87711
Atlantic	(-Inf,26]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	26	0.2080469	6.145234	0.6119006	2000-2009 –> 2010-2019	0.7213564	59.25247
Atlantic	(26,26.5]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	195	4.2358016	1130.379338	9.5121560	1982-1999 –> 2000-2009	7.3628248	2879.39492
Atlantic	(26,26.5]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	195	4.2121666	1128.197119	10.4631952	1982-1999 –> 2000-2009	7.2964370	2867.87909
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate	195	4.3374769	1139.630229	9.9260162	1982-1999 –> 2000-2009	7.5720102	2911.56434
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	195	4.3511486	1140.857569	10.0334319	1982-1999 –> 2000-2009	7.5962934	2915.01237
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + aou + nitrate + silicate	195	4.4744559	1151.756086	10.7896397	1982-1999 –> 2000-2009	7.9164457	2965.83812
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	195	4.3050687	1136.705331	10.9027558	1982-1999 –> 2000-2009	7.4929724	2898.79415
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	195	4.3272566	1138.710193	11.0127213	1982-1999 –> 2000-2009	7.5347245	2904.94719
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	195	4.3784146	1143.293838	10.0170066	1982-1999 –> 2000-2009	7.6468844	2922.30444
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	195	4.3443157	1140.244645	11.3072321	1982-1999 –> 2000-2009	7.5851708	2913.50260
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	195	4.3360442	1139.501385	10.9822045	1982-1999 –> 2000-2009	7.5511513	2907.35125
Atlantic	(26,26.5]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	280	3.6791944	1538.114111	12.9487717	2000-2009 –> 2010-2019	7.9149960	2668.49345
Atlantic	(26,26.5]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	280	3.5179481	1513.017200	14.3863133	2000-2009 –> 2010-2019	7.7301147	2641.21432
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate	280	3.6885346	1539.533960	12.9504155	2000-2009 –> 2010-2019	8.0260115	2679.16419
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	280	3.6904295	1539.821576	12.9472472	2000-2009 –> 2010-2019	8.0415781	2680.67915
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	280	3.5280245	1514.618911	14.3935647	2000-2009 –> 2010-2019	7.8330932	2651.32424
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	280	3.5316531	1515.194567	14.3791024	2000-2009 –> 2010-2019	7.8589096	2653.90476
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	280	3.6950679	1540.524975	12.9230879	2000-2009 –> 2010-2019	8.0734825	2683.81881
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	280	3.4879104	1508.215153	15.1596208	2000-2009 –> 2010-2019	7.8322261	2648.45980
Atlantic	(26,26.5]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	280	3.5337546	1515.527700	14.3625309	2000-2009 –> 2010-2019	7.8697988	2655.02908
Atlantic	(26,26.5]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	280	3.5480863	1517.794279	17.0800491	2000-2009 –> 2010-2019	8.0752334	2674.11617
Atlantic	(26.5,26.75]	5	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	266	3.0832651	1367.901514	9.7437457	1982-1999 –> 2000-2009	5.7507884	3181.96826
Atlantic	(26.5,26.75]	4	cstar_tref ~ aou + silicate + phosphate + phosphate_star	266	3.1720764	1381.008861	7.8194486	1982-1999 –> 2000-2009	5.8827065	3205.09858
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	266	2.9959557	1352.619372	8.6471002	1982-1999 –> 2000-2009	5.9234581	3236.43548
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	266	2.9772616	1349.289400	7.0206396	1982-1999 –> 2000-2009	5.6878307	3175.36228
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	266	2.9759262	1349.050729	7.0393155	1982-1999 –> 2000-2009	5.7419127	3190.30265
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	266	2.9763191	1349.120955	7.1250744	1982-1999 –> 2000-2009	5.7760658	3199.47158
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	266	2.9768698	1349.219384	7.1462711	1982-1999 –> 2000-2009	5.7885455	3202.75877
Atlantic	(26.5,26.75]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	266	3.1148181	1373.318142	7.7324317	1982-1999 –> 2000-2009	5.8443928	3204.63140
Atlantic	(26.5,26.75]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	266	3.1295698	1375.831733	7.7973121	1982-1999 –> 2000-2009	5.8943301	3216.75108
Atlantic	(26.5,26.75]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	266	3.1255830	1375.153581	10.1985469	1982-1999 –> 2000-2009	5.7547879	3178.36858
Atlantic	(26.5,26.75]	5	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	399	4.0481592	2262.126237	11.6682058	2000-2009 –> 2010-2019	7.1314243	3630.02775
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	399	3.9454953	2241.627400	10.3542399	2000-2009 –> 2010-2019	7.1174889	3624.62238
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	399	4.1747037	2286.689569	12.6310803	2000-2009 –> 2010-2019	7.1519653	3635.97897
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate	399	4.0497790	2262.445478	11.1047560	2000-2009 –> 2010-2019	7.1990919	3641.62282
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	399	4.0610781	2264.668832	11.1620397	2000-2009 –> 2010-2019	7.2024003	3642.49464
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	399	4.2468642	2300.365289	13.1771635	2000-2009 –> 2010-2019	7.2227904	3649.41602
Atlantic	(26.5,26.75]	5	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	399	4.0887332	2270.084630	11.3809634	2000-2009 –> 2010-2019	7.2207378	3646.33011
Atlantic	(26.5,26.75]	5	cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star	399	4.0014728	2252.869624	9.4335770	2000-2009 –> 2010-2019	7.2156011	3642.88480
Atlantic	(26.5,26.75]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	399	4.0879171	2269.925340	11.7699890	2000-2009 –> 2010-2019	7.2027352	3643.24348
Atlantic	(26.5,26.75]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	399	4.0446741	2261.438923	11.1809128	2000-2009 –> 2010-2019	7.1702571	3636.59250
Atlantic	(26.75,27]	5	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	437	2.2563607	1965.372563	6.7544677	1982-1999 –> 2000-2009	3.9467597	5000.19804
Atlantic	(26.75,27]	4	cstar_tref ~ aou + silicate + phosphate + phosphate_star	437	2.3258848	1989.896145	10.1912322	1982-1999 –> 2000-2009	4.0464360	5050.19648
Atlantic	(26.75,27]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	437	2.2874229	1977.322446	9.5012148	1982-1999 –> 2000-2009	4.0072407	5038.96023
Atlantic	(26.75,27]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	437	2.4136640	2024.273815	9.4271827	1982-1999 –> 2000-2009	4.2924413	5223.28960
Atlantic	(26.75,27]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	437	2.4901171	2051.528462	9.5746019	1982-1999 –> 2000-2009	4.4450789	5312.31493
Atlantic	(26.75,27]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	437	2.3904831	2015.839335	6.1759019	1982-1999 –> 2000-2009	4.2399778	5190.44376
Atlantic	(26.75,27]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	437	2.4414961	2034.294310	5.9087131	1982-1999 –> 2000-2009	4.3637303	5268.84532
Atlantic	(26.75,27]	4	cstar_tref ~ temp + aou + silicate + phosphate	437	2.5022807	2053.787344	10.6705256	1982-1999 –> 2000-2009	4.3866544	5255.42523
Atlantic	(26.75,27]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	437	2.2192662	1950.884622	9.5895072	1982-1999 –> 2000-2009	3.8043227	4885.71866
Atlantic	(26.75,27]	5	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	437	2.4936948	2052.783295	6.1766706	1982-1999 –> 2000-2009	4.4554475	5318.95848
Atlantic	(26.75,27]	5	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	718	2.4769462	3354.085665	8.0320187	2000-2009 –> 2010-2019	4.7333068	5319.45823
Atlantic	(26.75,27]	4	cstar_tref ~ aou + silicate + phosphate + phosphate_star	718	2.9696399	3612.596599	13.9488001	2000-2009 –> 2010-2019	5.2955248	5602.49274
Atlantic	(26.75,27]	5	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	718	2.6719290	3462.897516	8.7435497	2000-2009 –> 2010-2019	5.3010271	5561.89389
Atlantic	(26.75,27]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	718	2.8360788	3548.514300	13.5935133	2000-2009 –> 2010-2019	5.1235017	5525.83675
Atlantic	(26.75,27]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	718	2.8589493	3560.047931	13.6254278	2000-2009 –> 2010-2019	5.2726133	5584.32175
Atlantic	(26.75,27]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	718	2.8870239	3574.080516	13.7078252	2000-2009 –> 2010-2019	5.3771410	5625.60898
Atlantic	(26.75,27]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	718	2.5100961	3373.176717	8.2219439	2000-2009 –> 2010-2019	4.9005792	5389.01605
Atlantic	(26.75,27]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	718	2.5228723	3380.467341	8.2625661	2000-2009 –> 2010-2019	4.9643684	5414.76165
Atlantic	(26.75,27]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	718	2.9258385	3593.258214	13.5518755	2000-2009 –> 2010-2019	5.1451047	5544.14284
Atlantic	(26.75,27]	5	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	718	2.5532499	3397.654789	8.4070927	2000-2009 –> 2010-2019	5.0469447	5450.43808
Atlantic	(27,27.25]	5	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	532	1.8720745	2190.928805	6.0350504	1982-1999 –> 2000-2009	3.2982588	5657.01820
Atlantic	(27,27.25]	5	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	532	2.0076108	2265.300465	7.1133256	1982-1999 –> 2000-2009	3.5036649	5824.46402
Atlantic	(27,27.25]	4	cstar_tref ~ sal + aou + silicate + phosphate	532	2.1943693	2357.942538	9.9168850	1982-1999 –> 2000-2009	3.7017445	5929.77663
Atlantic	(27,27.25]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	532	2.1821648	2354.008305	10.2133301	1982-1999 –> 2000-2009	3.5584434	5750.78241
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	532	2.1904752	2358.052687	10.0860683	1982-1999 –> 2000-2009	3.5981713	5798.75032
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	532	2.2258811	2375.113201	10.0970598	1982-1999 –> 2000-2009	3.6771068	5875.07460
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	532	2.2344064	2379.180642	10.1032253	1982-1999 –> 2000-2009	3.6983967	5896.18371
Atlantic	(27,27.25]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	532	1.9807149	2250.949758	6.7823323	1982-1999 –> 2000-2009	3.5121789	5855.63632
Atlantic	(27,27.25]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	532	1.9864705	2254.037066	6.7937677	1982-1999 –> 2000-2009	3.5533153	5903.16980
Atlantic	(27,27.25]	5	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	532	2.0677086	2296.683886	7.2196765	1982-1999 –> 2000-2009	3.6774913	5998.42730
Atlantic	(27,27.25]	5	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	993	1.9788983	4187.536960	9.0909913	2000-2009 –> 2010-2019	3.8509729	6378.46576
Atlantic	(27,27.25]	5	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	993	2.0230370	4231.347209	9.0217246	2000-2009 –> 2010-2019	4.0306478	6496.64767
Atlantic	(27,27.25]	5	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	993	2.1242832	4328.332679	9.5947536	2000-2009 –> 2010-2019	4.3456159	6701.26948
Atlantic	(27,27.25]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	993	2.2411881	4434.726072	11.6181637	2000-2009 –> 2010-2019	4.4233528	6788.73438
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + aou + nitrate + silicate	993	2.0125378	4221.013403	9.7400104	2000-2009 –> 2010-2019	3.9982529	6474.64579
Atlantic	(27,27.25]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	993	2.2307291	4425.436317	11.3758512	2000-2009 –> 2010-2019	4.4212043	6783.48900
Atlantic	(27,27.25]	4	cstar_tref ~ temp + aou + nitrate + silicate	993	2.1023896	4305.758085	9.6134739	2000-2009 –> 2010-2019	4.1043825	6566.07696
Atlantic	(27,27.25]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	993	2.0934467	4299.292271	9.7133649	2000-2009 –> 2010-2019	4.0741616	6550.24203
Atlantic	(27,27.25]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	993	2.0888173	4294.895559	9.7051217	2000-2009 –> 2010-2019	4.0752878	6548.93262
Atlantic	(27,27.25]	5	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	993	2.1740649	4374.336811	10.0889611	2000-2009 –> 2010-2019	4.2417735	6671.02070
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	1146	3.4217842	6085.738650	10.2052540	1982-1999 –> 2000-2009	4.8918674	10920.71486
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1146	3.7606715	6302.184677	15.0121679	1982-1999 –> 2000-2009	5.0768650	10841.70406
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	1146	3.3336920	6025.959379	9.4485997	1982-1999 –> 2000-2009	5.1325438	11400.22410
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1146	3.7223836	6278.729879	14.2697860	1982-1999 –> 2000-2009	4.9676432	10670.22095
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1146	3.5777878	6187.921929	15.5407138	1982-1999 –> 2000-2009	4.8677210	10673.59125
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	1146	3.1771583	5915.729771	9.9721428	1982-1999 –> 2000-2009	5.0203569	11355.06793
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	1146	3.4629436	6113.143799	13.8886795	1982-1999 –> 2000-2009	4.8280311	10750.12352
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1146	3.5824413	6190.901109	15.4457626	1982-1999 –> 2000-2009	4.8883080	10709.37338
Atlantic	(27.25,27.5]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	1146	3.5123400	6145.606545	8.3861841	1982-1999 –> 2000-2009	5.0947381	11177.30205
Atlantic	(27.25,27.5]	5	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	1146	3.0609354	5830.314603	7.5696476	1982-1999 –> 2000-2009	4.5828490	10757.87428
Atlantic	(27.25,27.5]	4	cstar_tref ~ sal + aou + nitrate + silicate	1481	3.3112919	7761.412285	9.4536252	2000-2009 –> 2010-2019	6.7875362	13881.34251
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	1481	3.1212077	7588.303651	9.2308580	2000-2009 –> 2010-2019	6.5429919	13674.04230
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	1481	3.0852387	7553.971095	9.4787069	2000-2009 –> 2010-2019	6.5352142	13658.51572
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + aou + nitrate + silicate	1481	3.2982541	7751.726777	8.9832350	2000-2009 –> 2010-2019	6.6723393	13805.29066
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1481	3.0105839	7481.416943	8.6069108	2000-2009 –> 2010-2019	6.5883717	13669.33887
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	1481	3.4931167	7921.748892	11.9393875	2000-2009 –> 2010-2019	6.6702751	13837.47866
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	1481	3.0103932	7481.229324	8.5526966	2000-2009 –> 2010-2019	6.4733367	13594.37312
Atlantic	(27.25,27.5]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1481	3.0117753	7482.588911	8.6149081	2000-2009 –> 2010-2019	6.5942166	13673.49002
Atlantic	(27.25,27.5]	4	cstar_tref ~ sal + temp + silicate + phosphate_star	1481	3.0202896	7488.950709	8.9268723	2000-2009 –> 2010-2019	6.7623769	13777.78086
Atlantic	(27.25,27.5]	5	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	1481	3.5290741	7952.083217	9.6500374	2000-2009 –> 2010-2019	6.5900095	13782.39782
Atlantic	(27.5,27.75]	4	cstar_tref ~ sal + aou + nitrate + phosphate	2232	2.5957161	10604.183493	15.8720202	1982-1999 –> 2000-2009	4.1571778	20848.89393
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	2232	2.5847821	10587.339973	15.3700496	1982-1999 –> 2000-2009	3.9662612	20161.94796
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	2232	2.5753150	10570.959912	16.3852428	1982-1999 –> 2000-2009	4.1253119	20777.22688
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate	2232	2.5697378	10561.282028	15.4758679	1982-1999 –> 2000-2009	3.9036484	19943.59140
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	2232	2.5179329	10470.370144	17.5512775	1982-1999 –> 2000-2009	3.8649132	19906.18956
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	2232	2.5503283	10527.437055	15.8715208	1982-1999 –> 2000-2009	4.0741275	20640.15398
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	2232	2.4941516	10428.008227	18.2048947	1982-1999 –> 2000-2009	3.9159150	20160.35876
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	2232	2.5227542	10478.909507	17.0120152	1982-1999 –> 2000-2009	3.8226608	19719.50487
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	2232	2.4966766	10432.525252	17.4199840	1982-1999 –> 2000-2009	3.9130570	20144.05787
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	2232	2.4981354	10435.132827	17.9164450	1982-1999 –> 2000-2009	3.9190869	20164.34856
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	2445	2.2522681	10922.985093	10.5282563	2000-2009 –> 2010-2019	4.8370502	21510.32507
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	2445	2.0380307	10434.211067	9.2449209	2000-2009 –> 2010-2019	4.6133456	21005.17098
Atlantic	(27.5,27.75]	4	cstar_tref ~ sal + aou + silicate + phosphate	2445	2.1082755	10597.915276	9.1080725	2000-2009 –> 2010-2019	4.6863009	21171.57087
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2445	2.0911881	10560.120822	9.0148908	2000-2009 –> 2010-2019	4.6584814	21117.15444
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate	2445	2.2258838	10865.362628	10.5219647	2000-2009 –> 2010-2019	4.7956215	21426.64466
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	2445	2.0929814	10564.312368	9.3559664	2000-2009 –> 2010-2019	4.6433097	21091.74942
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	2445	2.1072455	10597.525692	9.2920764	2000-2009 –> 2010-2019	4.6013971	21025.53392
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	2445	2.3393694	11108.529575	11.4438471	2000-2009 –> 2010-2019	4.8621236	21587.43908
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	2445	2.1183958	10623.332544	9.3062709	2000-2009 –> 2010-2019	4.6150725	21055.85780
Atlantic	(27.5,27.75]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	2445	2.1049524	10592.201378	9.2385623	2000-2009 –> 2010-2019	4.6030878	21027.33421
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	943	1.2008724	3035.347091	6.6927010	1982-1999 –> 2000-2009	1.8069102	4817.65706
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	943	1.1308599	2922.055262	5.2848712	1982-1999 –> 2000-2009	1.9439582	5270.43497
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	943	1.1108068	2888.311457	5.4067728	1982-1999 –> 2000-2009	1.7690986	4829.91288
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate	943	1.2270444	3076.009445	6.8197469	1982-1999 –> 2000-2009	1.8530080	4920.62529
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	943	0.9505889	2594.547618	4.0066623	1982-1999 –> 2000-2009	1.6724408	4713.67189
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	943	1.1248948	2912.080473	5.4321788	1982-1999 –> 2000-2009	1.8070002	4922.12417
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	943	1.1767772	2997.120258	5.6967216	1982-1999 –> 2000-2009	1.8977700	5113.95094
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	943	1.2683518	3138.454705	7.1217548	1982-1999 –> 2000-2009	1.9172956	5052.51037
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	943	1.1727790	2990.701513	5.6610359	1982-1999 –> 2000-2009	1.8943764	5109.14657
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	943	1.1722742	2989.889508	5.6523079	1982-1999 –> 2000-2009	1.8929717	5105.93118
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	949	1.1285278	2936.639978	4.8014315	2000-2009 –> 2010-2019	2.3294001	5971.98707
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	949	1.2286944	3098.042522	4.9562206	2000-2009 –> 2010-2019	2.3595544	6020.09778
Atlantic	(27.75,27.85]	4	cstar_tref ~ sal + aou + silicate + phosphate	949	1.3353100	3253.977583	4.7733283	2000-2009 –> 2010-2019	2.4695161	6179.60523
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	949	1.2134912	3074.411203	5.3458160	2000-2009 –> 2010-2019	2.3242981	5962.72266
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate	949	1.1654208	2997.695327	4.9226766	2000-2009 –> 2010-2019	2.3924652	6073.70477
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	949	1.4216877	3374.946579	8.1359293	2000-2009 –> 2010-2019	2.3722766	5969.49420
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	949	1.2572889	3141.707054	5.3430480	2000-2009 –> 2010-2019	2.3821836	6053.78753
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	949	1.1939155	3043.543488	4.9754674	2000-2009 –> 2010-2019	2.4622673	6181.99819
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	949	1.3208536	3235.317303	5.5135714	2000-2009 –> 2010-2019	2.4936326	6226.01882
Atlantic	(27.75,27.85]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	949	1.3146656	3226.404572	5.4986088	2000-2009 –> 2010-2019	2.4869398	6216.29408
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	836	0.9149792	2237.901342	10.3082182	1982-1999 –> 2000-2009	1.8187176	4958.50496
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	836	1.1188958	2574.301628	12.0400737	1982-1999 –> 2000-2009	2.1624689	5590.40546
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate	836	0.9455565	2292.864005	10.6986475	1982-1999 –> 2000-2009	1.8715518	5063.43952
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	836	1.0499028	2467.887544	11.8633623	1982-1999 –> 2000-2009	2.0416164	5379.29625
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	836	0.8651218	2144.217577	7.9471856	1982-1999 –> 2000-2009	1.7859937	4903.39716
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	836	1.1444639	2612.078713	12.2962059	1982-1999 –> 2000-2009	2.2214064	5692.83339
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	836	0.9729963	2340.694221	10.9020013	1982-1999 –> 2000-2009	1.9208367	5159.16298
Atlantic	(27.85,27.95]	5	cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star	836	1.1889522	2675.842271	7.7193161	1982-1999 –> 2000-2009	2.2314409	5689.81073
Atlantic	(27.85,27.95]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	836	1.1954694	2682.982349	7.6077967	1982-1999 –> 2000-2009	2.2409866	5700.90896
Atlantic	(27.85,27.95]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	836	1.1819417	2665.954331	7.0484875	1982-1999 –> 2000-2009	2.2267576	5684.50299
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	941	1.9787747	3968.865535	16.1519445	2000-2009 –> 2010-2019	2.8937538	6206.76688
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	941	1.6597049	3637.938467	9.9823244	2000-2009 –> 2010-2019	2.7786008	6212.24009
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate	941	2.0815787	4064.186586	17.5322112	2000-2009 –> 2010-2019	3.0271353	6357.05059
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	941	1.5550503	3515.360193	9.1470242	2000-2009 –> 2010-2019	2.4201721	5659.57777
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	941	1.7396021	3726.423680	10.0969727	2000-2009 –> 2010-2019	2.8840660	6338.50239
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	941	1.8046688	3795.532038	10.1543304	2000-2009 –> 2010-2019	3.0018193	6482.86381
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	941	1.8100710	3801.157217	10.1705079	2000-2009 –> 2010-2019	3.0245175	6512.47281
Atlantic	(27.85,27.95]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	941	1.7969381	3787.452717	10.0977272	2000-2009 –> 2010-2019	2.9865419	6464.21110
Atlantic	(27.85,27.95]	5	cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star	941	1.8161834	3807.501869	8.8789568	2000-2009 –> 2010-2019	3.0051356	6483.34414
Atlantic	(27.85,27.95]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	941	1.7962850	3786.768593	9.6316612	2000-2009 –> 2010-2019	2.9782267	6452.72292
Atlantic	(27.95,28.05]	5	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	995	2.1335555	4345.689433	15.5189572	1982-1999 –> 2000-2009	3.8596003	9385.57383
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	995	2.5489717	4699.710857	11.4984169	1982-1999 –> 2000-2009	4.8255993	10447.80995
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	995	2.6477576	4775.376741	9.4725070	1982-1999 –> 2000-2009	5.0596858	10671.15255
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	995	2.8023488	4888.298948	18.2726781	1982-1999 –> 2000-2009	4.9450405	10481.30042
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	995	2.9651998	5000.707070	18.9729074	1982-1999 –> 2000-2009	5.2041151	10706.07815
Atlantic	(27.95,28.05]	5	cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star	995	2.7050518	4817.978610	18.5390472	1982-1999 –> 2000-2009	4.9404853	10519.36854
Atlantic	(27.95,28.05]	4	cstar_tref ~ temp + aou + nitrate + phosphate_star	995	2.8747252	4937.042306	19.1756649	1982-1999 –> 2000-2009	5.1155386	10642.58100
Atlantic	(27.95,28.05]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	995	2.0564931	4272.481971	15.6866947	1982-1999 –> 2000-2009	3.7084393	9200.12533
Atlantic	(27.95,28.05]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	995	1.8564843	4068.869896	15.4184543	1982-1999 –> 2000-2009	3.3794431	8788.55019
Atlantic	(27.95,28.05]	5	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	995	2.0307903	4247.453480	15.8306654	1982-1999 –> 2000-2009	3.6569656	9134.87672
Atlantic	(27.95,28.05]	5	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	1159	2.4687603	5397.913511	10.1800328	2000-2009 –> 2010-2019	4.6023158	9743.60294
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	1159	2.7421718	5641.382592	9.9061783	2000-2009 –> 2010-2019	5.2911436	10341.09345
Atlantic	(27.95,28.05]	4	cstar_tref ~ sal + nitrate + silicate + phosphate	1159	2.7851801	5675.456020	10.9222188	2000-2009 –> 2010-2019	5.5663214	10546.63784
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	1159	2.7746533	5668.678298	10.5530196	2000-2009 –> 2010-2019	5.4224109	10444.05504
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	1159	2.8287698	5713.453071	11.8030835	2000-2009 –> 2010-2019	5.6311186	10601.75202
Atlantic	(27.95,28.05]	5	cstar_tref ~ sal + temp + nitrate + silicate + phosphate	1159	2.6988957	5604.508905	10.5210115	2000-2009 –> 2010-2019	5.4524757	10457.87133
Atlantic	(27.95,28.05]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	1159	2.4431995	5373.788483	10.2445534	2000-2009 –> 2010-2019	4.4996926	9646.27045
Atlantic	(27.95,28.05]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	1159	2.4283795	5359.685133	10.1657196	2000-2009 –> 2010-2019	4.2848638	9428.55503
Atlantic	(27.95,28.05]	4	cstar_tref ~ temp + nitrate + silicate + phosphate	1159	2.7132164	5614.775963	10.6009794	2000-2009 –> 2010-2019	5.5070650	10495.02957
Atlantic	(27.95,28.05]	5	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	1159	2.4383406	5369.174023	10.2715256	2000-2009 –> 2010-2019	4.4691309	9616.62750
Atlantic	(28.05,28.1]	5	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	565	0.9530368	1563.045607	8.4238686	1982-1999 –> 2000-2009	1.6865569	3185.15099
Atlantic	(28.05,28.1]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	565	1.0988011	1723.868825	11.0329995	1982-1999 –> 2000-2009	1.8124660	3306.18412
Atlantic	(28.05,28.1]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	565	1.0696196	1693.453007	10.3483683	1982-1999 –> 2000-2009	1.7837730	3276.76044
Atlantic	(28.05,28.1]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	565	1.0302574	1651.084320	11.3456194	1982-1999 –> 2000-2009	1.7081982	3158.93705
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star	565	0.9596324	1570.838916	11.2298811	1982-1999 –> 2000-2009	1.5986805	2993.02521
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	565	0.9667475	1579.186225	8.2426612	1982-1999 –> 2000-2009	1.7192799	3238.39598
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	565	1.0089752	1627.497311	8.0254745	1982-1999 –> 2000-2009	1.8094728	3376.30159
Atlantic	(28.05,28.1]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	565	1.0495706	1670.071178	11.3060108	1982-1999 –> 2000-2009	1.7462070	3215.36911
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	565	0.9041957	1503.598893	9.9780290	1982-1999 –> 2000-2009	1.5342295	2905.18597
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	565	0.9727422	1586.171613	8.1484037	1982-1999 –> 2000-2009	1.7354398	3264.83680
Atlantic	(28.05,28.1]	5	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	709	1.3459604	2447.353664	7.2268997	2000-2009 –> 2010-2019	2.2989972	4010.39927
Atlantic	(28.05,28.1]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	709	1.3604520	2462.539361	7.2737180	2000-2009 –> 2010-2019	2.4592532	4186.40819
Atlantic	(28.05,28.1]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	709	1.3665441	2468.874983	7.2076672	2000-2009 –> 2010-2019	2.4361637	4162.32799
Atlantic	(28.05,28.1]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	709	1.2024266	2287.451398	6.8358347	2000-2009 –> 2010-2019	2.2326840	3938.53572
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star	709	1.0632146	2112.973983	6.1715212	2000-2009 –> 2010-2019	2.0228470	3683.81290
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	709	1.3751209	2477.746909	7.2525742	2000-2009 –> 2010-2019	2.3418684	4056.93313
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	709	1.4142077	2517.490280	7.1695910	2000-2009 –> 2010-2019	2.4231829	4144.98759
Atlantic	(28.05,28.1]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	709	1.2089146	2293.081900	7.0646333	2000-2009 –> 2010-2019	2.2584851	3963.15308
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	709	1.0823156	2138.222627	6.1089412	2000-2009 –> 2010-2019	1.9865113	3641.82152
Atlantic	(28.05,28.1]	5	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	709	1.3840146	2486.888411	7.2353901	2000-2009 –> 2010-2019	2.3567568	4073.06002
Atlantic	(28.1,28.15]	5	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	908	1.1067176	2774.932137	9.0447692	1982-1999 –> 2000-2009	1.7262965	4379.79001
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	908	1.2080347	2934.006965	8.7563730	1982-1999 –> 2000-2009	1.8735137	4659.78699
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	908	1.2881407	3050.603389	8.5301866	1982-1999 –> 2000-2009	1.9011658	4637.46797
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	908	1.2931994	3057.721025	8.4302400	1982-1999 –> 2000-2009	1.9156353	4670.36306
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	908	1.2958815	3061.483572	8.4132194	1982-1999 –> 2000-2009	1.9270548	4697.71178
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star	908	1.3558501	3143.634772	12.7091836	1982-1999 –> 2000-2009	1.8730694	4442.96294
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	908	1.1097470	2779.896202	8.9216934	1982-1999 –> 2000-2009	1.7463943	4430.73554
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	908	1.1263375	2806.844064	8.4411038	1982-1999 –> 2000-2009	1.7963633	4544.14533
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	908	1.2036210	2927.359882	12.8240665	1982-1999 –> 2000-2009	1.7089992	4187.50153
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	908	1.1101453	2780.547793	8.9121384	1982-1999 –> 2000-2009	1.7553281	4453.92109
Atlantic	(28.1,28.15]	5	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	935	1.0702994	2794.459956	4.9598195	2000-2009 –> 2010-2019	2.1770171	5569.39209
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	935	1.1033410	2851.316199	5.3405839	2000-2009 –> 2010-2019	2.3914818	5901.91959
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	935	0.9649802	2600.753814	4.6024595	2000-2009 –> 2010-2019	2.3587381	5794.46479
Atlantic	(28.1,28.15]	5	cstar_tref ~ sal + temp + nitrate + silicate + phosphate	935	1.1493847	2927.769060	3.1696545	2000-2009 –> 2010-2019	2.2949818	5765.40304
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star	935	0.8646157	2395.385782	4.4165643	2000-2009 –> 2010-2019	2.2204657	5539.02055
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	935	1.0941208	2835.623656	4.9066811	2000-2009 –> 2010-2019	2.2038678	5615.51986
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	935	1.1286149	2893.668538	5.0116212	2000-2009 –> 2010-2019	2.2549524	5700.51260
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	935	0.8367170	2334.051331	3.9023945	2000-2009 –> 2010-2019	2.0403380	5261.41121
Atlantic	(28.1,28.15]	4	cstar_tref ~ temp + nitrate + silicate + phosphate	935	1.1505673	2927.692132	3.1816529	2000-2009 –> 2010-2019	2.2972426	5765.03447
Atlantic	(28.1,28.15]	5	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	935	1.0996996	2845.134285	4.8277561	2000-2009 –> 2010-2019	2.2098448	5625.68208
Atlantic	(28.15,28.2]	5	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	1598	0.9948721	4532.496509	11.7853333	1982-1999 –> 2000-2009	1.6885761	7693.54369
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	1598	1.0000672	4549.142268	11.3734931	1982-1999 –> 2000-2009	1.6471742	7502.42275
Atlantic	(28.15,28.2]	4	cstar_tref ~ sal + temp + aou + nitrate	1598	1.0207059	4612.427873	10.3934017	1982-1999 –> 2000-2009	1.6748479	7596.01682
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate	1598	1.0145632	4595.135952	11.2048624	1982-1999 –> 2000-2009	1.6618833	7549.40005
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	1598	1.0170476	4602.952383	10.9937643	1982-1999 –> 2000-2009	1.6626483	7549.26956
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + aou + nitrate + silicate	1598	1.0204592	4613.655376	10.5439663	1982-1999 –> 2000-2009	1.6729567	7591.72296
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	1598	0.8241320	3930.742454	22.0243055	1982-1999 –> 2000-2009	1.6078329	7456.27226
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	1598	1.0219820	4618.421102	11.2930549	1982-1999 –> 2000-2009	1.6753683	7600.55588
Atlantic	(28.15,28.2]	5	cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star	1598	0.7410099	3590.954221	14.1622973	1982-1999 –> 2000-2009	1.3094493	6156.93943
Atlantic	(28.15,28.2]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1598	0.8299492	3953.222509	13.4004578	1982-1999 –> 2000-2009	1.6328807	7551.18698
Atlantic	(28.15,28.2]	4	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	1799	0.8421752	4499.322320	5.0205298	2000-2009 –> 2010-2019	1.8442221	9052.78520
Atlantic	(28.15,28.2]	5	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	1799	0.8356059	4473.146562	5.0080576	2000-2009 –> 2010-2019	1.8304780	9005.64307
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	1799	0.8412541	4497.384982	5.0120045	2000-2009 –> 2010-2019	1.8413213	9046.52725
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate	1799	0.8620040	4585.054546	4.9749893	2000-2009 –> 2010-2019	1.8765672	9180.19050
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + aou + nitrate + silicate	1799	0.8535794	4549.717377	4.9644937	2000-2009 –> 2010-2019	1.8740387	9163.37275
Atlantic	(28.15,28.2]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	1799	0.6180280	3387.905622	2.7434494	2000-2009 –> 2010-2019	1.4421599	7318.64808
Atlantic	(28.15,28.2]	5	cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star	1799	0.5590174	3026.834931	3.5666234	2000-2009 –> 2010-2019	1.3000272	6617.78915
Atlantic	(28.15,28.2]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	1799	0.8536781	4550.133405	4.9663029	2000-2009 –> 2010-2019	1.8673856	9142.57256
Atlantic	(28.15,28.2]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	1799	0.8528647	4546.703651	4.9278645	2000-2009 –> 2010-2019	1.8754000	9166.85441
Atlantic	(28.15,28.2]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1799	0.7989017	4311.527439	3.2000413	2000-2009 –> 2010-2019	1.6288509	8264.74995
Atlantic	(28.2, Inf]	4	cstar_tref ~ sal + aou + nitrate + phosphate	2503	0.7433965	5630.798439	3.8567442	1982-1999 –> 2000-2009	1.2539903	9285.48189
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	2503	0.7426981	5628.093397	3.7381585	1982-1999 –> 2000-2009	1.2519256	9271.70654
Atlantic	(28.2, Inf]	4	cstar_tref ~ sal + aou + nitrate + phosphate_star	2503	0.7427187	5626.231880	3.7225076	1982-1999 –> 2000-2009	1.2612431	9356.09953
Atlantic	(28.2, Inf]	4	cstar_tref ~ sal + aou + nitrate + silicate	2503	0.7280124	5526.115632	3.7500934	1982-1999 –> 2000-2009	1.2626263	9405.04325
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	2503	0.7078047	5387.197168	3.8952298	1982-1999 –> 2000-2009	1.1769241	8630.63190
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	2503	0.7188280	5464.559417	3.7938297	1982-1999 –> 2000-2009	1.2261529	9089.91294
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate	2503	0.7424889	5626.682955	3.7234726	1982-1999 –> 2000-2009	1.2521063	9274.02968
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	2503	0.7427126	5628.191003	3.7143011	1982-1999 –> 2000-2009	1.2533461	9285.25363
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + aou + nitrate + silicate	2503	0.7265701	5518.188502	3.7520993	1982-1999 –> 2000-2009	1.2610625	9398.00673
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	2503	0.7456953	5648.254909	3.7066617	1982-1999 –> 2000-2009	1.2517208	9261.09802
Atlantic	(28.2, Inf]	4	cstar_tref ~ sal + aou + nitrate + phosphate	4184	0.6260611	7966.882435	3.5083580	2000-2009 –> 2010-2019	1.3694576	13597.68087
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	4184	0.6210576	7901.736290	3.7727480	2000-2009 –> 2010-2019	1.3637557	13529.82969
Atlantic	(28.2, Inf]	4	cstar_tref ~ sal + aou + nitrate + phosphate_star	4184	0.6210622	7899.797458	3.7812284	2000-2009 –> 2010-2019	1.3637808	13526.02934
Atlantic	(28.2, Inf]	4	cstar_tref ~ sal + aou + nitrate + silicate	4184	0.6072181	7711.156949	3.5915793	2000-2009 –> 2010-2019	1.3352305	13237.27258
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	4184	0.5268904	6525.775488	2.7635469	2000-2009 –> 2010-2019	1.2346951	11912.97266
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	4184	0.5658093	7122.116242	3.8433141	2000-2009 –> 2010-2019	1.2846372	12586.67566
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate	4184	0.6214339	7906.805352	3.7601977	2000-2009 –> 2010-2019	1.3639228	13533.48831
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	4184	0.6210308	7901.375522	3.7592164	2000-2009 –> 2010-2019	1.3637434	13529.56652
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + aou + nitrate + silicate	4184	0.6055580	7690.247086	3.7505322	2000-2009 –> 2010-2019	1.3321281	13208.43559
Atlantic	(28.2, Inf]	5	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	4184	0.6236731	7936.903077	3.7401585	2000-2009 –> 2010-2019	1.3693685	13585.15799
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	1870	5.5723358	11745.455701	27.9741821	1982-1999 –> 2000-2009	10.7930679	18559.92534
Indo-Pacific	(-Inf,26]	4	cstar_tref ~ sal + aou + phosphate + phosphate_star	1870	5.9369868	11980.524616	24.3591049	1982-1999 –> 2000-2009	11.2229843	18820.50034
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1870	5.6773350	11815.272547	22.8831084	1982-1999 –> 2000-2009	10.8455704	18607.36680
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate	1870	5.8719382	11941.321155	32.0912255	1982-1999 –> 2000-2009	10.7938984	18625.31805
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	1870	5.9779183	12008.220906	35.1896101	1982-1999 –> 2000-2009	11.0594274	18762.84753
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	1870	5.9112576	11966.281277	25.2616093	1982-1999 –> 2000-2009	10.9471159	18700.92812
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1870	6.0078479	12026.899164	22.9150820	1982-1999 –> 2000-2009	10.9871046	18736.52039
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1870	6.1502900	12114.537248	23.4211186	1982-1999 –> 2000-2009	11.3358115	18914.02425
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	1870	6.0414750	12047.774327	33.1165944	1982-1999 –> 2000-2009	11.1846578	18829.11028
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1870	6.1766158	12130.511836	23.4218468	1982-1999 –> 2000-2009	11.4008529	18946.46738
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	2886	4.1302229	16390.721920	15.7392806	2000-2009 –> 2010-2019	9.7025587	28136.17762
Indo-Pacific	(-Inf,26]	4	cstar_tref ~ sal + aou + phosphate + phosphate_star	2886	5.4725451	18013.034340	23.3019252	2000-2009 –> 2010-2019	11.4095319	29993.55896
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2886	5.1162063	17626.402210	23.0261491	2000-2009 –> 2010-2019	10.7935413	29441.67476
Indo-Pacific	(-Inf,26]	4	cstar_tref ~ sal + temp + aou + nitrate	2886	4.2991951	16620.158632	17.2633536	2000-2009 –> 2010-2019	11.2916590	29212.66410
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate	2886	4.2974429	16619.805608	17.5543020	2000-2009 –> 2010-2019	10.1693811	28561.12676
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	2886	4.2991757	16622.132523	17.2990517	2000-2009 –> 2010-2019	10.2770940	28630.35343
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + nitrate + silicate	2886	3.9649307	16154.976195	17.2402482	2000-2009 –> 2010-2019	10.4675221	28477.83717
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2886	5.3090882	17840.006038	24.4320335	2000-2009 –> 2010-2019	11.2203458	29806.28732
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	2886	5.4794782	18022.342174	30.5719314	2000-2009 –> 2010-2019	11.4873260	30049.24134
Indo-Pacific	(-Inf,26]	5	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	2886	4.4339077	16800.244808	18.6649548	2000-2009 –> 2010-2019	10.4753827	28848.01913
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	2140	4.3337619	12363.402871	14.6466597	1982-1999 –> 2000-2009	7.7925132	20489.93362
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2140	4.1648510	12193.249497	15.3956493	1982-1999 –> 2000-2009	8.0042505	20638.22573
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate	2140	4.3118926	12341.750126	15.9832553	1982-1999 –> 2000-2009	7.7358537	20437.44688
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	2140	4.4676017	12493.582365	16.7086034	1982-1999 –> 2000-2009	7.9870443	20673.16785
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	2140	4.1715467	12200.124791	17.9633425	1982-1999 –> 2000-2009	8.0675851	20689.76602
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	2140	4.3435829	12373.091000	18.1460689	1982-1999 –> 2000-2009	8.3218557	20926.43924
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	2140	4.5070941	12531.250121	16.8344788	1982-1999 –> 2000-2009	8.0461941	20727.82363
Indo-Pacific	(26,26.5]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	2140	4.3195906	12349.384446	13.5717701	1982-1999 –> 2000-2009	8.1850964	20815.02921
Indo-Pacific	(26,26.5]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	2140	4.3221480	12351.917578	13.3774609	1982-1999 –> 2000-2009	8.2039083	20830.36079
Indo-Pacific	(26,26.5]	5	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	2140	4.3381491	12367.733396	13.4826494	1982-1999 –> 2000-2009	8.2309728	20854.85705
Indo-Pacific	(26,26.5]	5	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	1932	4.0143112	10867.219888	21.9696794	2000-2009 –> 2010-2019	8.5256766	23402.52413
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1932	3.6420910	10491.222444	20.8529455	2000-2009 –> 2010-2019	7.8069420	22684.47194
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1932	3.6529167	10502.690757	20.9430510	2000-2009 –> 2010-2019	7.8244634	22702.81555
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1932	3.6649500	10515.398499	20.8589812	2000-2009 –> 2010-2019	8.0085329	22888.48950
Indo-Pacific	(26,26.5]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1932	3.6718783	10522.696187	20.8468210	2000-2009 –> 2010-2019	8.0477148	22927.45125
Indo-Pacific	(26,26.5]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	1932	3.9948587	10848.450179	21.6703492	2000-2009 –> 2010-2019	8.3144493	23197.83463
Indo-Pacific	(26,26.5]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	1932	3.9935119	10847.147322	21.5798292	2000-2009 –> 2010-2019	8.3156599	23199.06490
Indo-Pacific	(26,26.5]	4	cstar_tref ~ temp + aou + silicate + phosphate	1932	3.9958611	10847.419627	21.7933287	2000-2009 –> 2010-2019	8.6940367	23554.38245
Indo-Pacific	(26,26.5]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1932	3.9949527	10848.541120	21.7413033	2000-2009 –> 2010-2019	8.6856701	23550.70418
Indo-Pacific	(26,26.5]	5	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	1932	3.9939916	10847.611452	21.5268630	2000-2009 –> 2010-2019	8.3321407	23215.34485
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	1907	4.3964132	11073.560872	22.1807811	1982-1999 –> 2000-2009	7.6526864	17404.39912
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	1907	4.3134175	11000.871731	19.7327836	1982-1999 –> 2000-2009	7.4081318	17208.05264
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1907	4.2262421	10923.000148	17.2669617	1982-1999 –> 2000-2009	7.4411332	17222.75895
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate	1907	4.4329880	11105.159243	21.2070279	1982-1999 –> 2000-2009	7.6190576	17383.03488
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	1907	4.5113568	11171.996158	21.9247837	1982-1999 –> 2000-2009	7.8170664	17539.44909
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	1907	4.5506168	11205.043797	14.6476541	1982-1999 –> 2000-2009	7.9809067	17662.39123
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1907	4.4411063	11112.137530	17.1466081	1982-1999 –> 2000-2009	7.8841920	17578.53259
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	1907	4.5194424	11178.825773	21.9357588	1982-1999 –> 2000-2009	7.8400210	17557.18426
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	1907	4.4700813	11136.940303	23.4566610	1982-1999 –> 2000-2009	7.8243682	17539.84223
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	1907	4.5398664	11196.022976	24.1822569	1982-1999 –> 2000-2009	8.0116413	17682.58191
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	1788	3.1983994	9245.762358	19.6031693	2000-2009 –> 2010-2019	7.5948125	20319.32323
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1788	3.0234331	9044.585514	15.9755692	2000-2009 –> 2010-2019	7.2496753	19967.58566
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1788	3.1717048	9215.790973	17.5413553	2000-2009 –> 2010-2019	7.6128111	20327.92850
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1788	3.2164329	9265.868320	18.0082162	2000-2009 –> 2010-2019	7.7639673	20468.32775
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	1788	3.1513395	9192.755729	18.5764950	2000-2009 –> 2010-2019	7.7751464	20458.65411
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1788	3.2212411	9271.209993	17.9917661	2000-2009 –> 2010-2019	7.7756969	20479.47005
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	1788	3.2850676	9341.372894	20.2319469	2000-2009 –> 2010-2019	7.7551488	20478.31320
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	1788	3.3286702	9388.524810	20.6847218	2000-2009 –> 2010-2019	7.8685366	20584.54779
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1788	2.9841422	8997.809050	17.6374803	2000-2009 –> 2010-2019	7.5437976	20210.42092
Indo-Pacific	(26.5,26.75]	5	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	1788	3.3322430	9392.361052	20.6885733	2000-2009 –> 2010-2019	7.8837100	20598.11735
Indo-Pacific	(26.75,27]	4	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	2861	3.9172204	15943.883876	14.4947878	1982-1999 –> 2000-2009	6.9448118	25077.29226
Indo-Pacific	(26.75,27]	5	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	2861	3.9128245	15939.459133	14.4028520	1982-1999 –> 2000-2009	6.8819828	25004.51258
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	2861	3.7621437	15714.752932	16.2549700	1982-1999 –> 2000-2009	6.6261609	24649.65417
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate	2861	3.7117385	15637.571342	13.5830674	1982-1999 –> 2000-2009	6.6551111	24671.13693
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	2861	3.8026290	15775.999679	14.7400015	1982-1999 –> 2000-2009	6.9137102	25009.61073
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	2861	3.8057482	15780.691478	14.7076553	1982-1999 –> 2000-2009	6.9304968	25030.12684
Indo-Pacific	(26.75,27]	4	cstar_tref ~ sal + temp + nitrate + phosphate_star	2861	3.8189064	15798.440789	14.1423806	1982-1999 –> 2000-2009	7.0308304	25145.21013
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	2861	3.8161980	15796.381381	13.9784932	1982-1999 –> 2000-2009	7.0264397	25143.25930
Indo-Pacific	(26.75,27]	5	cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star	2861	3.9147755	15942.311507	14.9696270	1982-1999 –> 2000-2009	6.9332076	25066.78210
Indo-Pacific	(26.75,27]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	2861	3.9444009	15985.450216	16.7734682	1982-1999 –> 2000-2009	7.0377716	25198.45166
Indo-Pacific	(26.75,27]	4	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	2267	3.8095407	12504.218555	13.5902815	2000-2009 –> 2010-2019	7.7267611	28448.10243
Indo-Pacific	(26.75,27]	5	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	2267	3.5821384	12227.279888	12.7836849	2000-2009 –> 2010-2019	7.4949629	28166.73902
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	2267	3.6640147	12329.700887	14.2498065	2000-2009 –> 2010-2019	7.4261584	28044.45382
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2267	3.5568366	12195.155322	12.2576424	2000-2009 –> 2010-2019	7.7488796	28529.02451
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate	2267	3.8135377	12516.486054	14.2848980	2000-2009 –> 2010-2019	7.5252762	28154.05740
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	2267	3.8818301	12591.411424	14.3709517	2000-2009 –> 2010-2019	7.6844590	28367.41110
Indo-Pacific	(26.75,27]	5	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	2267	3.8614215	12567.521737	14.3444599	2000-2009 –> 2010-2019	7.6671697	28348.21321
Indo-Pacific	(26.75,27]	5	cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star	2267	3.6325097	12290.564002	14.0693671	2000-2009 –> 2010-2019	7.5472852	28232.87551
Indo-Pacific	(26.75,27]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	2267	3.6605854	12330.890327	13.7104753	2000-2009 –> 2010-2019	7.6049863	28316.34054
Indo-Pacific	(26.75,27]	5	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	2267	3.6937394	12366.318891	14.0296546	2000-2009 –> 2010-2019	7.7491331	28510.55854
Indo-Pacific	(27,27.25]	4	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	3710	3.2264502	19232.182168	10.1381957	1982-1999 –> 2000-2009	5.7083438	29231.13210
Indo-Pacific	(27,27.25]	5	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	3710	3.0493469	18815.285407	9.6613802	1982-1999 –> 2000-2009	5.4937635	28750.96155
Indo-Pacific	(27,27.25]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	3710	3.1455747	19045.818758	9.8881879	1982-1999 –> 2000-2009	5.5815741	28966.69702
Indo-Pacific	(27,27.25]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate	3710	3.2378490	19260.350422	10.4943300	1982-1999 –> 2000-2009	5.6438987	29128.15780
Indo-Pacific	(27,27.25]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	3710	3.3461856	19504.556245	11.4221962	1982-1999 –> 2000-2009	5.8304734	29509.64256
Indo-Pacific	(27,27.25]	5	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	3710	3.3545843	19523.156524	11.4541169	1982-1999 –> 2000-2009	5.8553349	29556.57794
Indo-Pacific	(27,27.25]	4	cstar_tref ~ temp + aou + nitrate + phosphate	3710	3.2930111	19383.697618	10.3992334	1982-1999 –> 2000-2009	5.7554645	29348.91243
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star	3710	3.2246429	19230.024805	10.4742323	1982-1999 –> 2000-2009	5.6766463	29178.99521
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	3710	3.1795160	19125.452987	10.3415987	1982-1999 –> 2000-2009	5.6274001	29067.21032
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	3710	3.3103902	19424.754097	10.9031059	1982-1999 –> 2000-2009	5.8374016	29502.99066
Indo-Pacific	(27,27.25]	4	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	3289	4.6262266	19421.573569	20.1115374	2000-2009 –> 2010-2019	7.8526768	38653.75574
Indo-Pacific	(27,27.25]	5	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	3289	4.4199106	19123.471494	17.2929892	2000-2009 –> 2010-2019	7.4692574	37938.75690
Indo-Pacific	(27,27.25]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	3289	4.4801041	19212.450872	21.8200567	2000-2009 –> 2010-2019	7.6256788	38258.26963
Indo-Pacific	(27,27.25]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate	3289	4.5165758	19265.784350	22.2689788	2000-2009 –> 2010-2019	7.7544248	38526.13477
Indo-Pacific	(27,27.25]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	3289	4.0620674	18568.108077	19.8504699	2000-2009 –> 2010-2019	7.7306061	38755.09906
Indo-Pacific	(27,27.25]	5	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	3289	4.5251549	19278.267325	22.3261346	2000-2009 –> 2010-2019	7.8797392	38801.42385
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star	3289	4.1402456	18693.504714	18.5080805	2000-2009 –> 2010-2019	7.3648885	37923.52952
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	3289	4.5772822	19353.609090	17.5135872	2000-2009 –> 2010-2019	7.7567982	38479.06208
Indo-Pacific	(27,27.25]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	3289	4.1511820	18708.857619	18.4368540	2000-2009 –> 2010-2019	7.8909729	39036.58189
Indo-Pacific	(27,27.25]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	3289	3.7204199	17990.194401	16.0235554	2000-2009 –> 2010-2019	7.2242503	37836.33549
Indo-Pacific	(27.25,27.5]	4	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	4723	3.1472975	24245.549506	12.7971055	1982-1999 –> 2000-2009	5.2338284	35926.00389
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	4723	3.0403422	23920.962762	10.9240181	1982-1999 –> 2000-2009	5.0961258	35523.01259
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	4723	3.1176462	24158.135017	12.2210141	1982-1999 –> 2000-2009	5.2035572	35838.97997
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	4723	3.0694002	24010.813825	9.5273998	1982-1999 –> 2000-2009	5.4022359	36297.59492
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	4723	2.9688129	23696.073359	11.0074317	1982-1999 –> 2000-2009	5.2076293	35760.05511
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	4723	3.0922603	24080.904496	10.2493784	1982-1999 –> 2000-2009	5.3547979	36201.96944
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star	4723	2.7287433	22899.577214	9.6595329	1982-1999 –> 2000-2009	4.7058288	34290.22279
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	4723	3.1867515	24365.226871	11.3039234	1982-1999 –> 2000-2009	5.3647053	36279.93560
Indo-Pacific	(27.25,27.5]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	4723	2.9712777	23701.912485	10.8876563	1982-1999 –> 2000-2009	5.2339631	35821.33109
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	4723	2.7707049	23043.728643	9.6043714	1982-1999 –> 2000-2009	4.9423566	34942.74297
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	4368	3.7845495	24036.824144	17.3333745	2000-2009 –> 2010-2019	6.8248918	47957.78691
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	4368	3.5590327	23500.101225	13.2428414	2000-2009 –> 2010-2019	6.6284329	47510.91505
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	4368	3.2335121	22662.144832	11.7854301	2000-2009 –> 2010-2019	6.4965054	47250.70260
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	4368	3.3167994	22884.313559	15.1119374	2000-2009 –> 2010-2019	6.8052638	48104.02160
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	4368	3.4009681	23103.236198	13.2049563	2000-2009 –> 2010-2019	6.6046059	47518.38443
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	4368	3.2689247	22757.299252	14.6826416	2000-2009 –> 2010-2019	6.7439132	47940.44665
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	4368	3.2907765	22815.502366	15.0448476	2000-2009 –> 2010-2019	6.8201812	48145.42242
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + nitrate + silicate + phosphate	4368	3.1983178	22566.538815	15.7597511	2000-2009 –> 2010-2019	6.9017129	48351.00779
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	4368	3.2726452	22767.236369	14.8286143	2000-2009 –> 2010-2019	6.7811356	48041.01524
Indo-Pacific	(27.25,27.5]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	4368	3.2622072	22739.328511	14.3843788	2000-2009 –> 2010-2019	6.0329121	45783.05715
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	5652	3.1825663	29140.184926	9.8980753	1982-1999 –> 2000-2009	5.1996815	42140.78176
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	5652	2.9103964	28129.623409	9.2295589	1982-1999 –> 2000-2009	4.8731636	40962.95482
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	5652	2.9165840	28153.630326	9.4572980	1982-1999 –> 2000-2009	4.9471957	41195.06640
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	5652	2.9110502	28132.162520	9.4452271	1982-1999 –> 2000-2009	4.9499242	41198.46542
Indo-Pacific	(27.5,27.75]	4	cstar_tref ~ sal + nitrate + silicate + phosphate_star	5652	3.0338245	28597.132184	13.9982016	1982-1999 –> 2000-2009	5.0867953	41703.63097
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	5652	3.0181919	28540.734643	19.1407962	1982-1999 –> 2000-2009	5.1344754	41835.24124
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	5652	2.9979523	28464.675993	9.5171895	1982-1999 –> 2000-2009	5.2187162	42054.29666
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	5652	2.9622712	28329.331037	10.6633756	1982-1999 –> 2000-2009	5.0129371	41430.95045
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star	5652	3.1528062	29033.984421	26.5711843	1982-1999 –> 2000-2009	5.2203833	42185.89923
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	5652	2.8663266	27957.146644	24.9178880	1982-1999 –> 2000-2009	4.9029830	41016.78536
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	4860	3.2606307	25294.351323	11.5635098	2000-2009 –> 2010-2019	6.1772147	53447.98165
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	4860	2.6257985	23189.625178	13.0042990	2000-2009 –> 2010-2019	5.6512071	51757.35614
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + nitrate + silicate	4860	2.6570208	23304.519827	14.6256682	2000-2009 –> 2010-2019	5.8687844	52547.93663
Indo-Pacific	(27.5,27.75]	4	cstar_tref ~ sal + temp + aou + silicate	4860	2.6732581	23361.738900	15.2584478	2000-2009 –> 2010-2019	5.8960627	52641.94847
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	4860	2.6172499	23157.928707	12.5120436	2000-2009 –> 2010-2019	5.6152021	51622.60470
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	4860	2.6178554	23160.177429	12.5711647	2000-2009 –> 2010-2019	5.6313174	51683.18336
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + nitrate + silicate + phosphate	4860	2.7107047	23498.950540	8.2454565	2000-2009 –> 2010-2019	5.7254255	52026.67731
Indo-Pacific	(27.5,27.75]	4	cstar_tref ~ sal + temp + silicate + phosphate	4860	2.8062588	23833.686009	8.5870530	2000-2009 –> 2010-2019	5.8450927	52449.46766
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	4860	2.6164184	23154.840191	12.4566518	2000-2009 –> 2010-2019	5.6219578	51648.08812
Indo-Pacific	(27.5,27.75]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	4860	3.3662154	25604.112326	14.4834017	2000-2009 –> 2010-2019	6.2325419	53561.25897
Indo-Pacific	(27.75,27.85]	4	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	2200	2.4195854	10143.152896	9.0483206	1982-1999 –> 2000-2009	3.8744003	14711.44848
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	2200	2.4178792	10142.049090	8.8306720	1982-1999 –> 2000-2009	3.8561906	14683.36614
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	2200	2.2805231	9884.710922	8.2886047	1982-1999 –> 2000-2009	3.6657422	14330.49485
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate	2200	2.3061575	9933.893449	8.6243669	1982-1999 –> 2000-2009	3.7410840	14469.22595
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	2200	2.3138702	9948.584258	8.6965226	1982-1999 –> 2000-2009	3.7740742	14528.27159
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2200	2.1949560	9716.442286	10.4214802	1982-1999 –> 2000-2009	3.4822528	13976.00934
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	2200	2.3194630	9959.206584	8.7663173	1982-1999 –> 2000-2009	3.7938657	14563.47306
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star	2200	2.2620542	9848.932269	12.0172199	1982-1999 –> 2000-2009	3.4012273	13798.00480
Indo-Pacific	(27.75,27.85]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	2200	2.3189494	9956.232205	11.8352213	1982-1999 –> 2000-2009	3.6065197	14214.33882
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	2200	2.3143046	9949.410198	11.6909072	1982-1999 –> 2000-2009	3.5975529	14200.97650
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	1919	2.8157007	9433.026407	10.3322542	2000-2009 –> 2010-2019	5.0962238	19317.73733
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1919	2.5018255	8979.411438	10.5185792	2000-2009 –> 2010-2019	4.8302028	18955.49591
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate	1919	2.8365637	9461.359341	9.9177022	2000-2009 –> 2010-2019	5.1427211	19395.25279
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	1919	2.8464872	9474.762921	9.7134205	2000-2009 –> 2010-2019	5.1603574	19423.34718
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + nitrate + silicate	1919	2.7906176	9398.683269	11.2303409	2000-2009 –> 2010-2019	5.1719930	19473.79693
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	1919	2.8183936	9436.695308	17.5273535	2000-2009 –> 2010-2019	5.0133496	19153.13759
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1919	2.5158120	9000.808062	9.9660519	2000-2009 –> 2010-2019	4.8635856	19013.39443
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1919	2.5202431	9007.561974	9.7958752	2000-2009 –> 2010-2019	4.8752656	19033.71280
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	1919	2.8518373	9481.969854	9.6262169	2000-2009 –> 2010-2019	5.1713003	19441.17644
Indo-Pacific	(27.75,27.85]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1919	2.5229163	9011.630795	9.7084621	2000-2009 –> 2010-2019	4.8803662	19042.31448
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	2628	1.9679106	11030.107586	6.9896634	1982-1999 –> 2000-2009	3.5036057	16133.28385
Indo-Pacific	(27.85,27.95]	4	cstar_tref ~ sal + temp + aou + nitrate	2628	1.9983204	11108.706679	7.1967445	1982-1999 –> 2000-2009	3.6464026	16404.39566
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate	2628	1.9958604	11104.232269	7.3529297	1982-1999 –> 2000-2009	3.5746210	16283.57595
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	2628	1.9968602	11106.864617	7.3287790	1982-1999 –> 2000-2009	3.5855210	16303.42430
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + nitrate + silicate	2628	1.9693862	11034.047413	8.0404805	1982-1999 –> 2000-2009	3.5960146	16295.65098
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2628	1.9015317	10849.760492	8.6525819	1982-1999 –> 2000-2009	3.2732335	15641.92485
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	2628	2.0066584	11132.591647	7.4654486	1982-1999 –> 2000-2009	3.6115182	16357.09074
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star	2628	1.8797016	10789.070946	9.8070025	1982-1999 –> 2000-2009	3.1777181	15429.17372
Indo-Pacific	(27.85,27.95]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	2628	1.9399668	10952.938997	9.5635078	1982-1999 –> 2000-2009	3.3160397	15751.86534
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	2628	1.9089746	10870.293216	9.1395530	1982-1999 –> 2000-2009	3.2587329	15618.04473
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	2345	2.7122119	11348.337208	14.5138688	2000-2009 –> 2010-2019	4.6801224	22378.44479
Indo-Pacific	(27.85,27.95]	4	cstar_tref ~ sal + aou + nitrate + phosphate_star	2345	2.7134036	11348.397443	14.8775353	2000-2009 –> 2010-2019	4.8658692	22847.66147
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	2345	2.7217829	11364.858366	15.2270936	2000-2009 –> 2010-2019	4.8158194	22721.47493
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	2345	2.6883428	11306.879668	15.2138380	2000-2009 –> 2010-2019	4.7488413	22578.63514
Indo-Pacific	(27.85,27.95]	4	cstar_tref ~ sal + temp + aou + nitrate	2345	2.7428291	11398.984319	11.9615508	2000-2009 –> 2010-2019	4.7411495	22507.69100
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate	2345	2.7149271	11353.030020	14.6215644	2000-2009 –> 2010-2019	4.7107874	22457.26229
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	2345	2.7129795	11349.664384	14.6721974	2000-2009 –> 2010-2019	4.7098397	22456.52900
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + aou + nitrate + silicate	2345	2.6850320	11301.100218	11.8834488	2000-2009 –> 2010-2019	4.6544182	22335.14763
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	2345	2.7193969	11360.745197	14.7544507	2000-2009 –> 2010-2019	4.7260553	22493.33684
Indo-Pacific	(27.85,27.95]	5	cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star	2345	2.8146176	11522.157640	14.2602429	2000-2009 –> 2010-2019	4.6943192	22311.22859
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	2918	1.6912785	11361.654264	8.9619075	1982-1999 –> 2000-2009	3.4260542	18073.03096
Indo-Pacific	(27.95,28.05]	4	cstar_tref ~ sal + temp + aou + nitrate	2918	1.6993817	11387.548690	8.7993916	1982-1999 –> 2000-2009	3.4717165	18169.75214
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate	2918	1.6980639	11385.021406	9.0021398	1982-1999 –> 2000-2009	3.4521589	18134.05278
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	2918	1.6983084	11385.861849	8.9835131	1982-1999 –> 2000-2009	3.4548598	18139.65108
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + nitrate + silicate	2918	1.6631032	11263.612580	8.8182407	1982-1999 –> 2000-2009	3.3840560	17947.78900
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	2918	1.7070955	11415.979337	10.2194785	1982-1999 –> 2000-2009	3.3855056	18015.05011
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	2918	1.7014919	11396.791138	9.0006823	1982-1999 –> 2000-2009	3.4659648	18165.87900
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star	2918	1.7260762	11480.510098	11.4581605	1982-1999 –> 2000-2009	3.4005818	18071.66227
Indo-Pacific	(27.95,28.05]	4	cstar_tref ~ temp + aou + phosphate + phosphate_star	2918	1.7277226	11484.074069	11.4795795	1982-1999 –> 2000-2009	3.4084709	18085.87820
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	2918	1.6381933	11175.539950	9.6305648	1982-1999 –> 2000-2009	3.1813329	17488.91571
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	3148	1.5260597	11608.887494	10.5059765	2000-2009 –> 2010-2019	3.2173382	22970.54176
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	3148	1.4714317	11379.378179	9.5780200	2000-2009 –> 2010-2019	3.2248879	22951.73575
Indo-Pacific	(27.95,28.05]	4	cstar_tref ~ sal + temp + aou + nitrate	3148	1.5262886	11607.831747	10.9819324	2000-2009 –> 2010-2019	3.2256703	22995.38044
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate	3148	1.5170142	11571.457600	10.4502664	2000-2009 –> 2010-2019	3.2150781	22956.47901
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	3148	1.5183530	11577.011664	10.5053206	2000-2009 –> 2010-2019	3.2166615	22962.87351
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + nitrate + silicate	3148	1.4377505	11233.586883	9.8466889	2000-2009 –> 2010-2019	3.1008537	22497.19946
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	3148	1.5403869	11667.720702	10.1002219	2000-2009 –> 2010-2019	3.2474824	23083.70004
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	3148	1.4716705	11380.399586	9.5390125	2000-2009 –> 2010-2019	3.2454564	23020.03114
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	3148	1.5147700	11562.136889	10.4644530	2000-2009 –> 2010-2019	3.2162620	22958.92803
Indo-Pacific	(27.95,28.05]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	3148	1.5941545	11883.735586	10.3246913	2000-2009 –> 2010-2019	3.2323478	23059.27554
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	1642	1.5719412	6159.184327	7.1081093	1982-1999 –> 2000-2009	3.1353019	10154.75112
Indo-Pacific	(28.05,28.1]	4	cstar_tref ~ aou + silicate + phosphate + phosphate_star	1642	1.5777348	6169.265879	7.0240945	1982-1999 –> 2000-2009	3.1423492	10164.54320
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1642	1.5581124	6130.166312	7.0145523	1982-1999 –> 2000-2009	3.1011379	10097.79321
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1642	1.5728107	6161.000385	7.0282181	1982-1999 –> 2000-2009	3.1266067	10143.47107
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1642	1.5769953	6169.726267	7.0660555	1982-1999 –> 2000-2009	3.1352743	10158.34505
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	1642	1.5771329	6170.012730	7.0503737	1982-1999 –> 2000-2009	3.1240409	10143.00228
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1642	1.5771001	6169.944339	7.0640396	1982-1999 –> 2000-2009	3.1344740	10157.32338
Indo-Pacific	(28.05,28.1]	4	cstar_tref ~ sal + temp + silicate + phosphate_star	1642	1.5774940	6168.764586	7.0426149	1982-1999 –> 2000-2009	3.1589266	10186.85809
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	1642	1.5955807	6208.202922	7.1183867	1982-1999 –> 2000-2009	3.1749353	10225.49055
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1642	1.4497738	5893.496442	7.1264250	1982-1999 –> 2000-2009	2.8737959	9689.84897
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	1799	1.4954226	6567.207918	7.8875747	2000-2009 –> 2010-2019	3.0673638	12726.39224
Indo-Pacific	(28.05,28.1]	4	cstar_tref ~ aou + silicate + phosphate + phosphate_star	1799	1.4956731	6565.810455	7.8899099	2000-2009 –> 2010-2019	3.0734079	12735.07633
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1799	1.3068998	6082.373571	6.6993870	2000-2009 –> 2010-2019	2.8650122	12212.53988
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + aou + nitrate + silicate	1799	1.3177498	6112.120962	6.4893078	2000-2009 –> 2010-2019	2.9667025	12428.37606
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	1799	1.3069266	6082.447358	6.6873236	2000-2009 –> 2010-2019	2.8797373	12243.44774
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1799	1.3069792	6082.592171	6.6826474	2000-2009 –> 2010-2019	2.8839746	12252.31844
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	1799	1.3071988	6083.196580	6.7069823	2000-2009 –> 2010-2019	2.8843317	12253.20931
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1799	1.3069632	6082.547930	6.6835097	2000-2009 –> 2010-2019	2.8840632	12252.49227
Indo-Pacific	(28.05,28.1]	4	cstar_tref ~ sal + temp + silicate + phosphate_star	1799	1.3237771	6126.540712	7.3000093	2000-2009 –> 2010-2019	2.9012712	12295.30530
Indo-Pacific	(28.05,28.1]	5	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1799	1.3849403	6291.054830	6.8391738	2000-2009 –> 2010-2019	2.8347141	12184.55127
Indo-Pacific	(28.1, Inf]	4	cstar_tref ~ aou + nitrate + silicate + phosphate	19077	1.4742839	68960.509250	7.7290890	1982-1999 –> 2000-2009	2.8336159	110498.43221
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	19077	1.4673168	68781.777146	7.7847613	1982-1999 –> 2000-2009	2.7659481	109222.57658
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	19077	1.3607566	65905.171996	6.6602134	1982-1999 –> 2000-2009	2.5971623	105164.57013
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	19077	1.3799626	66439.921165	6.9379940	1982-1999 –> 2000-2009	2.7239495	107706.69283
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate	19077	1.5229608	70201.902111	8.6576732	1982-1999 –> 2000-2009	2.8372375	110930.83615
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	19077	1.5232047	70208.013249	8.6600432	1982-1999 –> 2000-2009	2.8364841	110918.68165
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	19077	1.5231891	70207.622864	8.6603162	1982-1999 –> 2000-2009	2.8365137	110919.11899
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	19077	1.4696196	68841.606815	7.7899186	1982-1999 –> 2000-2009	2.7735575	109380.52450
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	19077	1.4846262	69229.231531	7.8885521	1982-1999 –> 2000-2009	2.8015459	110006.50124
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	19077	1.4713268	68885.903253	7.7801027	1982-1999 –> 2000-2009	2.7804013	109519.41296
Indo-Pacific	(28.1, Inf]	4	cstar_tref ~ aou + nitrate + silicate + phosphate	19505	1.4075670	68700.854742	9.8671496	2000-2009 –> 2010-2019	2.8818509	137661.36399
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	19505	1.4020505	68549.668314	10.0753048	2000-2009 –> 2010-2019	2.8693674	137331.44546
Indo-Pacific	(28.1, Inf]	4	cstar_tref ~ nitrate + silicate + phosphate + phosphate_star	19505	1.4111982	68801.363733	9.6610597	2000-2009 –> 2010-2019	2.8917523	137923.79835
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	19505	1.2746677	64833.952896	7.9553093	2000-2009 –> 2010-2019	2.6354243	130739.12489
Indo-Pacific	(28.1, Inf]	4	cstar_tref ~ sal + nitrate + silicate + phosphate	19505	1.3507758	67094.285031	7.3489055	2000-2009 –> 2010-2019	2.8345035	136298.41577
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	19505	1.2876306	65228.667385	7.6306703	2000-2009 –> 2010-2019	2.6675932	131668.58855
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ sal + temp + nitrate + silicate + phosphate	19505	1.3425830	66858.956750	7.0512062	2000-2009 –> 2010-2019	2.8071297	135568.63524
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	19505	1.4040420	68605.040376	10.0458278	2000-2009 –> 2010-2019	2.8736616	137446.64719
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	19505	1.4202989	69054.127097	10.1604997	2000-2009 –> 2010-2019	2.9049251	138283.35863
Indo-Pacific	(28.1, Inf]	5	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	19505	1.4052286	68637.993821	10.0153330	2000-2009 –> 2010-2019	2.8765554	137523.89707

5.4.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.2656000	271.74157
Atlantic	(-Inf,26]	2000-2009 –> 2010-2019	0.7476676	67.76988
Atlantic	(26,26.5]	1982-1999 –> 2000-2009	7.5554914	2908.65884
Atlantic	(26,26.5]	2000-2009 –> 2010-2019	7.9255444	2663.62040
Atlantic	(26.5,26.75]	1982-1999 –> 2000-2009	5.8044818	3199.11487
Atlantic	(26.5,26.75]	2000-2009 –> 2010-2019	7.1834492	3639.32135
Atlantic	(26.75,27]	1982-1999 –> 2000-2009	4.1988089	5154.43507
Atlantic	(26.75,27]	2000-2009 –> 2010-2019	5.1160112	5501.79710
Atlantic	(27,27.25]	1982-1999 –> 2000-2009	3.5778772	5848.92833
Atlantic	(27,27.25]	2000-2009 –> 2010-2019	4.1565652	6595.95244
Atlantic	(27.25,27.5]	1982-1999 –> 2000-2009	4.9350924	10925.61964
Atlantic	(27.25,27.5]	2000-2009 –> 2010-2019	6.6216668	13735.40505
Atlantic	(27.5,27.75]	1982-1999 –> 2000-2009	3.9662160	20246.62738
Atlantic	(27.5,27.75]	2000-2009 –> 2010-2019	4.6915790	21201.87804
Atlantic	(27.75,27.85]	1982-1999 –> 2000-2009	1.8454830	4985.59653
Atlantic	(27.75,27.85]	2000-2009 –> 2010-2019	2.4072534	6085.57103
Atlantic	(27.85,27.95]	1982-1999 –> 2000-2009	2.0521777	5382.22624
Atlantic	(27.85,27.95]	2000-2009 –> 2010-2019	2.8999969	6316.97525
Atlantic	(27.95,28.05]	1982-1999 –> 2000-2009	4.4694913	9997.74167
Atlantic	(27.95,28.05]	2000-2009 –> 2010-2019	5.0726538	10132.14952
Atlantic	(28.05,28.1]	1982-1999 –> 2000-2009	1.7134304	3192.01473
Atlantic	(28.05,28.1]	2000-2009 –> 2010-2019	2.2816750	3986.14394
Atlantic	(28.1,28.15]	1982-1999 –> 2000-2009	1.8223821	4510.43862
Atlantic	(28.1,28.15]	2000-2009 –> 2010-2019	2.2448930	5653.83603
Atlantic	(28.15,28.2]	1982-1999 –> 2000-2009	1.6233618	7424.73304
Atlantic	(28.15,28.2]	2000-2009 –> 2010-2019	1.7380451	8595.91329
Atlantic	(28.2, Inf]	1982-1999 –> 2000-2009	1.2451098	9225.72641
Atlantic	(28.2, Inf]	2000-2009 –> 2010-2019	1.3380720	13224.71092
Indo-Pacific	(-Inf,26]	1982-1999 –> 2000-2009	11.0570491	18750.30085
Indo-Pacific	(-Inf,26]	2000-2009 –> 2010-2019	10.7294343	29115.69406
Indo-Pacific	(26,26.5]	1982-1999 –> 2000-2009	8.0575274	20708.30500
Indo-Pacific	(26,26.5]	2000-2009 –> 2010-2019	8.2555286	23132.30834
Indo-Pacific	(26.5,26.75]	1982-1999 –> 2000-2009	7.7479205	17477.82269
Indo-Pacific	(26.5,26.75]	2000-2009 –> 2010-2019	7.6823302	20389.26886
Indo-Pacific	(26.75,27]	1982-1999 –> 2000-2009	6.8980523	24999.60367
Indo-Pacific	(26.75,27]	2000-2009 –> 2010-2019	7.6175071	28311.77761
Indo-Pacific	(27,27.25]	1982-1999 –> 2000-2009	5.7010301	29224.12776
Indo-Pacific	(27,27.25]	2000-2009 –> 2010-2019	7.6549293	38420.89489
Indo-Pacific	(27.25,27.5]	1982-1999 –> 2000-2009	5.1745028	35688.18484
Indo-Pacific	(27.25,27.5]	2000-2009 –> 2010-2019	6.6639555	47660.27598
Indo-Pacific	(27.5,27.75]	1982-1999 –> 2000-2009	5.0546255	41572.40723
Indo-Pacific	(27.5,27.75]	2000-2009 –> 2010-2019	5.8264806	52338.65030
Indo-Pacific	(27.75,27.85]	1982-1999 –> 2000-2009	3.6792910	14347.56095
Indo-Pacific	(27.75,27.85]	2000-2009 –> 2010-2019	5.0205365	19224.93659
Indo-Pacific	(27.85,27.95]	1982-1999 –> 2000-2009	3.4543407	16021.84301
Indo-Pacific	(27.85,27.95]	2000-2009 –> 2010-2019	4.7347222	22508.74117
Indo-Pacific	(27.95,28.05]	1982-1999 –> 2000-2009	3.4030701	18029.16612
Indo-Pacific	(27.95,28.05]	2000-2009 –> 2010-2019	3.2142038	22945.61447
Indo-Pacific	(28.05,28.1]	1982-1999 –> 2000-2009	3.1106843	10112.14269
Indo-Pacific	(28.05,28.1]	2000-2009 –> 2010-2019	2.9240578	12358.37088
Indo-Pacific	(28.1, Inf]	1982-1999 –> 2000-2009	2.7786416	109426.73472
Indo-Pacific	(28.1, Inf]	2000-2009 –> 2010-2019	2.8242763	136044.52751

5.5 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.6 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.7 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.8 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
b883157	Donghe-Zhu	2021-03-27
d19654d	Donghe-Zhu	2021-03-26
6c53dbf	Donghe-Zhu	2021-03-25
1914a11	Donghe-Zhu	2021-03-24
8be810e	Donghe-Zhu	2021-03-23
bf19764	Donghe-Zhu	2021-03-22
3ec9d3d	Donghe-Zhu	2021-03-22
134ace1	Donghe-Zhu	2021-03-22
f6d70a4	Donghe-Zhu	2021-03-22
183443b	Donghe-Zhu	2021-03-21
2e6976b	Donghe-Zhu	2021-03-21
51a42bd	Donghe-Zhu	2021-03-16
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
b883157	Donghe-Zhu	2021-03-27
d19654d	Donghe-Zhu	2021-03-26
6c53dbf	Donghe-Zhu	2021-03-25
1914a11	Donghe-Zhu	2021-03-24
8be810e	Donghe-Zhu	2021-03-23
bf19764	Donghe-Zhu	2021-03-22
3ec9d3d	Donghe-Zhu	2021-03-22
134ace1	Donghe-Zhu	2021-03-22
f6d70a4	Donghe-Zhu	2021-03-22
183443b	Donghe-Zhu	2021-03-21
2e6976b	Donghe-Zhu	2021-03-21
51a42bd	Donghe-Zhu	2021-03-16
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
b883157	Donghe-Zhu	2021-03-27
d19654d	Donghe-Zhu	2021-03-26
6c53dbf	Donghe-Zhu	2021-03-25
1914a11	Donghe-Zhu	2021-03-24
8be810e	Donghe-Zhu	2021-03-23
bf19764	Donghe-Zhu	2021-03-22
3ec9d3d	Donghe-Zhu	2021-03-22
134ace1	Donghe-Zhu	2021-03-22
f6d70a4	Donghe-Zhu	2021-03-22
183443b	Donghe-Zhu	2021-03-21
2e6976b	Donghe-Zhu	2021-03-21
51a42bd	Donghe-Zhu	2021-03-16
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
b883157	Donghe-Zhu	2021-03-27
d19654d	Donghe-Zhu	2021-03-26
6c53dbf	Donghe-Zhu	2021-03-25
1914a11	Donghe-Zhu	2021-03-24
8be810e	Donghe-Zhu	2021-03-23
bf19764	Donghe-Zhu	2021-03-22
3ec9d3d	Donghe-Zhu	2021-03-22
134ace1	Donghe-Zhu	2021-03-22
f6d70a4	Donghe-Zhu	2021-03-22
183443b	Donghe-Zhu	2021-03-21
2e6976b	Donghe-Zhu	2021-03-21
51a42bd	Donghe-Zhu	2021-03-16
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	nitrate	phosphate	phosphate_star	sal	silicate	temp
Atlantic - 1982-1999 --> 2000-2009
(-Inf,26]	7	3	8	8	7	7	9
(26,26.5]	7	6	6	7	10	6	8
(26.5,26.75]	9	4	8	7	5	10	6
(26.75,27]	9	4	8	7	3	10	7
(27,27.25]	8	5	8	6	6	10	6
(27.25,27.5]	6	6	6	8	8	8	8
(27.5,27.75]	7	7	7	6	10	5	7
(27.75,27.85]	7	5	8	7	10	6	7
(27.85,27.95]	9	5	9	8	7	3	8
(27.95,28.05]	7	10	7	8	4	6	7
(28.05,28.1]	9	6	9	9	3	6	7
(28.1,28.15]	8	6	9	7	4	9	7
(28.15,28.2]	9	8	7	7	7	3	8
(28.2, Inf]	9	10	5	5	10	4	4
total	111.00	85.00	105.00	100.00	94.00	93.00	99.00
Atlantic - 2000-2009 --> 2010-2019
(-Inf,26]	6	4	5	5	9	10	9
(26,26.5]	7	6	6	8	9	6	8
(26.5,26.75]	9	7	9	7	6	5	7
(26.75,27]	9	5	8	7	4	10	6
(27,27.25]	9	8	6	5	5	10	6
(27.25,27.5]	5	7	4	7	9	9	7
(27.5,27.75]	7	5	8	6	10	7	6
(27.75,27.85]	7	5	9	6	10	6	6
(27.85,27.95]	8	4	8	8	8	6	8
(27.95,28.05]	5	10	8	5	5	9	6
(28.05,28.1]	9	6	9	9	3	6	7
(28.1,28.15]	7	7	9	7	3	8	8
(28.15,28.2]	10	8	8	7	4	5	7
(28.2, Inf]	9	10	5	5	10	4	4
total	107.00	92.00	102.00	92.00	95.00	101.00	95.00
Indo-Pacific - 1982-1999 --> 2000-2009
(-Inf,26]	8	4	8	8	10	4	7
(26,26.5]	8	7	7	7	7	6	8
(26.5,26.75]	9	7	8	7	7	5	7
(26.75,27]	7	10	7	8	6	3	7
(27,27.25]	9	10	8	7	4	3	7
(27.25,27.5]	9	7	9	8	4	5	6
(27.5,27.75]	7	7	7	8	8	7	5
(27.75,27.85]	9	7	9	9	5	2	7
(27.85,27.95]	9	7	7	7	7	2	9
(27.95,28.05]	9	7	7	7	7	2	9
(28.05,28.1]	7	3	7	8	6	10	7
(28.1, Inf]	7	7	8	6	5	10	6
total	98.00	83.00	92.00	90.00	76.00	59.00	85.00
Indo-Pacific - 2000-2009 --> 2010-2019
(-Inf,26]	9	6	7	6	10	3	7
(26,26.5]	8	4	8	7	4	10	8
(26.5,26.75]	7	5	7	8	5	10	8
(26.75,27]	8	9	9	8	5	4	6
(27,27.25]	9	7	10	8	4	3	7
(27.25,27.5]	7	5	8	7	8	10	5
(27.5,27.75]	7	3	6	5	9	10	8
(27.75,27.85]	8	5	7	7	10	5	8
(27.85,27.95]	9	10	5	6	9	3	6
(27.95,28.05]	9	6	7	6	9	4	8
(28.05,28.1]	7	3	6	8	7	10	7
(28.1, Inf]	5	10	9	5	4	10	4
total	93.00	73.00	89.00	81.00	84.00	82.00	82.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
b883157	Donghe-Zhu	2021-03-27
d19654d	Donghe-Zhu	2021-03-26
6c53dbf	Donghe-Zhu	2021-03-25
1914a11	Donghe-Zhu	2021-03-24
8be810e	Donghe-Zhu	2021-03-23
bf19764	Donghe-Zhu	2021-03-22
3ec9d3d	Donghe-Zhu	2021-03-22
134ace1	Donghe-Zhu	2021-03-22
f6d70a4	Donghe-Zhu	2021-03-22
183443b	Donghe-Zhu	2021-03-21
2e6976b	Donghe-Zhu	2021-03-21
51a42bd	Donghe-Zhu	2021-03-16
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
b883157	Donghe-Zhu	2021-03-27
d19654d	Donghe-Zhu	2021-03-26
6c53dbf	Donghe-Zhu	2021-03-25
1914a11	Donghe-Zhu	2021-03-24
8be810e	Donghe-Zhu	2021-03-23
bf19764	Donghe-Zhu	2021-03-22
3ec9d3d	Donghe-Zhu	2021-03-22
134ace1	Donghe-Zhu	2021-03-22
f6d70a4	Donghe-Zhu	2021-03-22
183443b	Donghe-Zhu	2021-03-21
2e6976b	Donghe-Zhu	2021-03-21
51a42bd	Donghe-Zhu	2021-03-16
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] car_3.0-10       carData_3.0-4    gt_0.2.2         corrr_0.4.3     
 [5] broom_0.7.5      kableExtra_1.3.1 knitr_1.30       olsrr_0.5.3     
 [9] GGally_2.0.0     lubridate_1.7.9  metR_0.9.0       scico_1.2.0     
[13] patchwork_1.1.1  collapse_1.5.0   forcats_0.5.0    stringr_1.4.0   
[17] dplyr_1.0.5      purrr_0.3.4      readr_1.4.0      tidyr_1.1.2     
[21] tibble_3.0.4     ggplot2_3.3.3    tidyverse_1.3.0  workflowr_1.6.2 

loaded via a namespace (and not attached):
 [1] fs_1.5.0                 webshot_0.5.2            RColorBrewer_1.1-2      
 [4] httr_1.4.2               rprojroot_2.0.2          tools_4.0.3             
 [7] backports_1.1.10         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                magrittr_2.0.1           Matrix_1.2-18           
[43] Rcpp_1.0.5               munsell_0.5.0            fansi_0.4.1             
[46] abind_1.4-5              lifecycle_1.0.0          stringi_1.5.3           
[49] whisker_0.4              yaml_2.2.1               plyr_1.8.6              
[52] grid_4.0.3               blob_1.2.1               parallel_4.0.3          
[55] promises_1.1.1           crayon_1.3.4             lattice_0.20-41         
[58] haven_2.3.1              hms_0.5.3                pillar_1.4.7            
[61] reprex_0.3.0             glue_1.4.2               evaluate_0.14           
[64] RcppArmadillo_0.10.1.2.2 data.table_1.13.6        modelr_0.1.8            
[67] vctrs_0.3.6              httpuv_1.5.4             cellranger_1.1.0        
[70] gtable_0.3.0             reshape_0.8.8            assertthat_0.2.1        
[73] xfun_0.20                openxlsx_4.2.3           RcppEigen_0.3.3.9.1     
[76] later_1.1.0.1            viridisLite_0.3.0        ellipsis_0.3.1          
[79] here_1.0.1