eMLR - model fitting

workflowr

• Summary
• Checks
• Past versions

Last updated: 2021-03-10

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: 38fa816

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 38fa816. See the Past versions tab to see a history of the changes made to the R Markdown and HTML files.

Note that you need to be careful to ensure that all relevant files for the analysis have been committed to Git prior to generating the results (you can use wflow_publish or wflow_git_commit). workflowr only checks the R Markdown file, but you know if there are other scripts or data files that it depends on. Below is the status of the Git repository when the results were generated:

Ignored files:
    Ignored:    .Rhistory
    Ignored:    .Rproj.user/

Unstaged changes:
    Deleted:    a.png
    Modified:   data/auxillary/params_local.rds

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

---

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

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

---

1 Required data

Required are:

• cleaned and prepared GLODAP-based synthetic model subsetting file

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

2 Predictor combinations

Find all possible combinations of following considered predictor variables:

• sal, temp, oxygen, aou, nitrate, silicate, silicate_star, 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()


# 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)

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

# 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: 9

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:

• 502

4 Fit all models

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

• cstar_tref

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

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

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

5 Prepare coeffcients

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

5.1 Formatting

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

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

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

5.2 Predictor selection

Within each basin and slab, the following number of best linear regression models was selected:

• 10

The criterion used to select the best models was:

• rmse

The criterion was summed up for two adjacent eras, and the models with lowest summed values were selected.

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

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

# calculate RMSE sum for adjacent eras
lm_all_fitted_wide_eras <- lm_all_fitted_wide %>%
  select(basin, gamma_slab, n_predictors, model, era, nr_obs, rmse, aic, resid_max) %>%
  arrange(era) %>%
  group_by(basin, gamma_slab, model) %>%
  mutate(
    eras = paste(lag(era), era, sep = " --> "),
    rmse_sum = rmse + lag(rmse),
    aic_sum = aic + lag(aic)
  ) %>%
  ungroup() %>%
  select(-c(era)) %>%
  drop_na()

# subset models with lowest summed criterion
# chose which criterion is applied

if (params_local$MLR_criterion == "aic") {
  lm_best <- lm_all_fitted_wide_eras %>%
    group_by(basin, gamma_slab, eras) %>%
    slice_min(order_by = aic_sum,
              with_ties = FALSE,
              n = params_local$MLR_number) %>%
    ungroup() %>%
    arrange(basin, gamma_slab, eras, model)
} else {
  lm_best <- lm_all_fitted_wide_eras %>%
    group_by(basin, gamma_slab, eras) %>%
    slice_min(order_by = rmse_sum,
              with_ties = FALSE,
              n = params_local$MLR_number) %>%
    ungroup() %>%
    arrange(basin, gamma_slab, eras, model)
}

5.3 RMSE tables

5.3.1 per model

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

basin	gamma_slab	n_predictors	model	nr_obs	rmse	aic	resid_max	eras	rmse_sum	aic_sum
Atlantic	(-Inf,26]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	158	1.0093505	469.3256	3.392692	1982-1999 –> 2000-2012	2.3449212	1197.9773
Atlantic	(-Inf,26]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	158	1.0093505	469.3256	3.392692	1982-1999 –> 2000-2012	2.3449212	1197.9773
Atlantic	(-Inf,26]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	158	1.0093505	469.3256	3.392692	1982-1999 –> 2000-2012	2.3449212	1197.9773
Atlantic	(-Inf,26]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	158	1.0093505	469.3256	3.392692	1982-1999 –> 2000-2012	2.3449212	1197.9773
Atlantic	(-Inf,26]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	158	1.0093505	469.3256	3.392692	1982-1999 –> 2000-2012	2.3449212	1197.9773
Atlantic	(-Inf,26]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate_star	158	1.0093505	469.3256	3.392692	1982-1999 –> 2000-2012	2.3449212	1197.9773
Atlantic	(-Inf,26]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	158	1.0093505	469.3256	3.392692	1982-1999 –> 2000-2012	2.3449212	1197.9773
Atlantic	(-Inf,26]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate_star	158	1.0093505	469.3256	3.392692	1982-1999 –> 2000-2012	2.3449212	1197.9773
Atlantic	(-Inf,26]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	158	1.0093505	469.3256	3.392692	1982-1999 –> 2000-2012	2.3449212	1197.9773
Atlantic	(-Inf,26]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	158	1.0093505	469.3256	3.392692	1982-1999 –> 2000-2012	2.3449212	1197.9773
Atlantic	(-Inf,26]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	72	1.0760286	232.8790	2.653279	2000-2012 –> 2013-2019	2.0853791	702.2046
Atlantic	(-Inf,26]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	72	1.0760286	232.8790	2.653279	2000-2012 –> 2013-2019	2.0853791	702.2046
Atlantic	(-Inf,26]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	72	1.0760286	232.8790	2.653279	2000-2012 –> 2013-2019	2.0853791	702.2046
Atlantic	(-Inf,26]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	72	1.0760286	232.8790	2.653279	2000-2012 –> 2013-2019	2.0853791	702.2046
Atlantic	(-Inf,26]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	72	1.0760286	232.8790	2.653279	2000-2012 –> 2013-2019	2.0853791	702.2046
Atlantic	(-Inf,26]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	72	1.0760286	232.8790	2.653279	2000-2012 –> 2013-2019	2.0853791	702.2046
Atlantic	(-Inf,26]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	72	1.0760286	232.8790	2.653279	2000-2012 –> 2013-2019	2.0853791	702.2046
Atlantic	(-Inf,26]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	72	1.0760286	232.8790	2.653279	2000-2012 –> 2013-2019	2.0853791	702.2046
Atlantic	(-Inf,26]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	72	1.0760286	232.8790	2.653279	2000-2012 –> 2013-2019	2.0853791	702.2046
Atlantic	(-Inf,26]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	72	1.0760286	232.8790	2.653279	2000-2012 –> 2013-2019	2.0853791	702.2046
Atlantic	(26,26.5]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate + phosphate_star	1046	3.2125528	5427.9212	13.822807	1982-1999 –> 2000-2012	6.4070517	13316.0460
Atlantic	(26,26.5]	8	cstar_tref ~ sal + temp + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	1046	3.2125528	5427.9212	13.822807	1982-1999 –> 2000-2012	6.4070517	13316.0460
Atlantic	(26,26.5]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate_star + phosphate + phosphate_star	1046	3.2125528	5427.9212	13.822807	1982-1999 –> 2000-2012	6.4070517	13316.0460
Atlantic	(26,26.5]	7	cstar_tref ~ sal + temp + aou + silicate + silicate_star + phosphate + phosphate_star	1046	3.2125528	5427.9212	13.822807	1982-1999 –> 2000-2012	6.4070517	13316.0460
Atlantic	(26,26.5]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	1046	3.2125528	5427.9212	13.822807	1982-1999 –> 2000-2012	6.4070517	13316.0460
Atlantic	(26,26.5]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	1046	3.2125528	5427.9212	13.822807	1982-1999 –> 2000-2012	6.4070517	13316.0460
Atlantic	(26,26.5]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	1046	3.2125528	5427.9212	13.822807	1982-1999 –> 2000-2012	6.4070517	13316.0460
Atlantic	(26,26.5]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	1046	3.2125528	5427.9212	13.822807	1982-1999 –> 2000-2012	6.4070517	13316.0460
Atlantic	(26,26.5]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	1046	3.2125528	5427.9212	13.822807	1982-1999 –> 2000-2012	6.4070517	13316.0460
Atlantic	(26,26.5]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	1046	3.2125528	5427.9212	13.822807	1982-1999 –> 2000-2012	6.4070517	13316.0460
Atlantic	(26,26.5]	7	cstar_tref ~ sal + temp + aou + silicate + silicate_star + phosphate + phosphate_star	567	3.0166904	2879.1941	10.327831	2000-2012 –> 2013-2019	6.2292432	8307.1154
Atlantic	(26,26.5]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	567	3.0166904	2879.1941	10.327831	2000-2012 –> 2013-2019	6.2292432	8307.1154
Atlantic	(26,26.5]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	567	3.0166904	2879.1941	10.327831	2000-2012 –> 2013-2019	6.2292432	8307.1154
Atlantic	(26,26.5]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	567	3.0166904	2879.1941	10.327831	2000-2012 –> 2013-2019	6.2292432	8307.1154
Atlantic	(26,26.5]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	567	3.0166904	2879.1941	10.327831	2000-2012 –> 2013-2019	6.2292432	8307.1154
Atlantic	(26,26.5]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	567	3.0166904	2879.1941	10.327831	2000-2012 –> 2013-2019	6.2292432	8307.1154
Atlantic	(26,26.5]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	567	3.0166904	2879.1941	10.327831	2000-2012 –> 2013-2019	6.2292432	8307.1154
Atlantic	(26,26.5]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	567	3.0166904	2879.1941	10.327831	2000-2012 –> 2013-2019	6.2292432	8307.1154
Atlantic	(26,26.5]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	567	3.0166904	2879.1941	10.327831	2000-2012 –> 2013-2019	6.2292432	8307.1154
Atlantic	(26,26.5]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	567	3.0166904	2879.1941	10.327831	2000-2012 –> 2013-2019	6.2292432	8307.1154
Atlantic	(26.5,26.75]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate + phosphate_star	1397	3.1332896	7173.4954	8.794157	1982-1999 –> 2000-2012	6.0826758	16678.5275
Atlantic	(26.5,26.75]	8	cstar_tref ~ sal + temp + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	1397	3.1332896	7173.4954	8.794157	1982-1999 –> 2000-2012	6.0826758	16678.5275
Atlantic	(26.5,26.75]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate_star + phosphate + phosphate_star	1397	3.1332896	7173.4954	8.794157	1982-1999 –> 2000-2012	6.0826758	16678.5275
Atlantic	(26.5,26.75]	7	cstar_tref ~ sal + temp + aou + silicate + silicate_star + phosphate + phosphate_star	1397	3.1332896	7173.4954	8.794157	1982-1999 –> 2000-2012	6.0826758	16678.5275
Atlantic	(26.5,26.75]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	1397	3.1332896	7173.4954	8.794157	1982-1999 –> 2000-2012	6.0826758	16678.5275
Atlantic	(26.5,26.75]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	1397	3.1332896	7173.4954	8.794157	1982-1999 –> 2000-2012	6.0826758	16678.5275
Atlantic	(26.5,26.75]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	1397	3.1332896	7173.4954	8.794157	1982-1999 –> 2000-2012	6.0826758	16678.5275
Atlantic	(26.5,26.75]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	1397	3.1332896	7173.4954	8.794157	1982-1999 –> 2000-2012	6.0826758	16678.5275
Atlantic	(26.5,26.75]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	1397	3.1332896	7173.4954	8.794157	1982-1999 –> 2000-2012	6.0826758	16678.5275
Atlantic	(26.5,26.75]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	1397	3.1332896	7173.4954	8.794157	1982-1999 –> 2000-2012	6.0826758	16678.5275
Atlantic	(26.5,26.75]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	656	3.2676102	3433.1326	9.323441	2000-2012 –> 2013-2019	6.4008998	10606.6280
Atlantic	(26.5,26.75]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	656	3.2676102	3433.1326	9.323441	2000-2012 –> 2013-2019	6.4008998	10606.6280
Atlantic	(26.5,26.75]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	656	3.2676102	3433.1326	9.323441	2000-2012 –> 2013-2019	6.4008998	10606.6280
Atlantic	(26.5,26.75]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	656	3.2676102	3433.1326	9.323441	2000-2012 –> 2013-2019	6.4008998	10606.6280
Atlantic	(26.5,26.75]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	656	3.2676102	3433.1326	9.323441	2000-2012 –> 2013-2019	6.4008998	10606.6280
Atlantic	(26.5,26.75]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate_star	656	3.2676102	3433.1326	9.323441	2000-2012 –> 2013-2019	6.4008998	10606.6280
Atlantic	(26.5,26.75]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate_star	656	3.2676102	3433.1326	9.323441	2000-2012 –> 2013-2019	6.4008998	10606.6280
Atlantic	(26.5,26.75]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	656	3.2676102	3433.1326	9.323441	2000-2012 –> 2013-2019	6.4008998	10606.6280
Atlantic	(26.5,26.75]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	656	3.2676102	3433.1326	9.323441	2000-2012 –> 2013-2019	6.4008998	10606.6280
Atlantic	(26.5,26.75]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	656	3.2676102	3433.1326	9.323441	2000-2012 –> 2013-2019	6.4008998	10606.6280
Atlantic	(26.75,27]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate + phosphate_star	2398	1.9577047	10045.0510	12.852844	1982-1999 –> 2000-2012	3.4513206	21493.5208
Atlantic	(26.75,27]	8	cstar_tref ~ sal + temp + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	2398	1.9577047	10045.0510	12.852844	1982-1999 –> 2000-2012	3.4513206	21493.5208
Atlantic	(26.75,27]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate_star + phosphate + phosphate_star	2398	1.9577047	10045.0510	12.852844	1982-1999 –> 2000-2012	3.4513206	21493.5208
Atlantic	(26.75,27]	7	cstar_tref ~ sal + temp + aou + silicate + silicate_star + phosphate + phosphate_star	2398	1.9577047	10045.0510	12.852844	1982-1999 –> 2000-2012	3.4513206	21493.5208
Atlantic	(26.75,27]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	2398	1.9577047	10045.0510	12.852844	1982-1999 –> 2000-2012	3.4513206	21493.5208
Atlantic	(26.75,27]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	2398	1.9577047	10045.0510	12.852844	1982-1999 –> 2000-2012	3.4513206	21493.5208
Atlantic	(26.75,27]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	2398	1.9577047	10045.0510	12.852844	1982-1999 –> 2000-2012	3.4513206	21493.5208
Atlantic	(26.75,27]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	2398	1.9577047	10045.0510	12.852844	1982-1999 –> 2000-2012	3.4513206	21493.5208
Atlantic	(26.75,27]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	2398	1.9577047	10045.0510	12.852844	1982-1999 –> 2000-2012	3.4513206	21493.5208
Atlantic	(26.75,27]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	2398	1.9577047	10045.0510	12.852844	1982-1999 –> 2000-2012	3.4513206	21493.5208
Atlantic	(26.75,27]	7	cstar_tref ~ sal + temp + aou + silicate + silicate_star + phosphate + phosphate_star	1263	2.1936686	5586.6019	19.571794	2000-2012 –> 2013-2019	4.1513732	15631.6529
Atlantic	(26.75,27]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	1263	2.1936686	5586.6019	19.571794	2000-2012 –> 2013-2019	4.1513732	15631.6529
Atlantic	(26.75,27]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	1263	2.1936686	5586.6019	19.571794	2000-2012 –> 2013-2019	4.1513732	15631.6529
Atlantic	(26.75,27]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	1263	2.1936686	5586.6019	19.571794	2000-2012 –> 2013-2019	4.1513732	15631.6529
Atlantic	(26.75,27]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	1263	2.1936686	5586.6019	19.571794	2000-2012 –> 2013-2019	4.1513732	15631.6529
Atlantic	(26.75,27]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate_star	1263	2.1936686	5586.6019	19.571794	2000-2012 –> 2013-2019	4.1513732	15631.6529
Atlantic	(26.75,27]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate_star	1263	2.1936686	5586.6019	19.571794	2000-2012 –> 2013-2019	4.1513732	15631.6529
Atlantic	(26.75,27]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	1263	2.1936686	5586.6019	19.571794	2000-2012 –> 2013-2019	4.1513732	15631.6529
Atlantic	(26.75,27]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	1263	2.1936686	5586.6019	19.571794	2000-2012 –> 2013-2019	4.1513732	15631.6529
Atlantic	(26.75,27]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	1263	2.1936686	5586.6019	19.571794	2000-2012 –> 2013-2019	4.1513732	15631.6529
Atlantic	(27,27.25]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	2031	1.7317419	8012.2853	15.886046	1982-1999 –> 2000-2012	3.1780603	17367.0116
Atlantic	(27,27.25]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	2031	1.7317419	8012.2853	15.886046	1982-1999 –> 2000-2012	3.1780603	17367.0116
Atlantic	(27,27.25]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	2031	1.7317419	8012.2853	15.886046	1982-1999 –> 2000-2012	3.1780603	17367.0116
Atlantic	(27,27.25]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	2031	1.7317419	8012.2853	15.886046	1982-1999 –> 2000-2012	3.1780603	17367.0116
Atlantic	(27,27.25]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	2031	1.7317419	8012.2853	15.886046	1982-1999 –> 2000-2012	3.1780603	17367.0116
Atlantic	(27,27.25]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate_star	2031	1.7317419	8012.2853	15.886046	1982-1999 –> 2000-2012	3.1780603	17367.0116
Atlantic	(27,27.25]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	2031	1.7317419	8012.2853	15.886046	1982-1999 –> 2000-2012	3.1780603	17367.0116
Atlantic	(27,27.25]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	2031	1.7317419	8012.2853	15.886046	1982-1999 –> 2000-2012	3.1780603	17367.0116
Atlantic	(27,27.25]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	2031	1.7317419	8012.2853	15.886046	1982-1999 –> 2000-2012	3.1780603	17367.0116
Atlantic	(27,27.25]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	2031	1.7317419	8012.2853	15.886046	1982-1999 –> 2000-2012	3.1780603	17367.0116
Atlantic	(27,27.25]	7	cstar_tref ~ sal + temp + aou + silicate + silicate_star + phosphate + phosphate_star	1142	2.0775479	4928.8897	9.783184	2000-2012 –> 2013-2019	3.8092897	12941.1749
Atlantic	(27,27.25]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	1142	2.0775479	4928.8897	9.783184	2000-2012 –> 2013-2019	3.8092897	12941.1749
Atlantic	(27,27.25]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	1142	2.0775479	4928.8897	9.783184	2000-2012 –> 2013-2019	3.8092897	12941.1749
Atlantic	(27,27.25]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	1142	2.0775479	4928.8897	9.783184	2000-2012 –> 2013-2019	3.8092897	12941.1749
Atlantic	(27,27.25]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	1142	2.0775479	4928.8897	9.783184	2000-2012 –> 2013-2019	3.8092897	12941.1749
Atlantic	(27,27.25]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	1142	2.0775479	4928.8897	9.783184	2000-2012 –> 2013-2019	3.8092897	12941.1749
Atlantic	(27,27.25]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate_star	1142	2.0775479	4928.8897	9.783184	2000-2012 –> 2013-2019	3.8092897	12941.1749
Atlantic	(27,27.25]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	1142	2.0775479	4928.8897	9.783184	2000-2012 –> 2013-2019	3.8092897	12941.1749
Atlantic	(27,27.25]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	1142	2.0775479	4928.8897	9.783184	2000-2012 –> 2013-2019	3.8092897	12941.1749
Atlantic	(27,27.25]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	1142	2.0775479	4928.8897	9.783184	2000-2012 –> 2013-2019	3.8092897	12941.1749
Atlantic	(27.25,27.5]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate + phosphate_star	2025	2.1391730	8844.3992	17.846793	1982-1999 –> 2000-2012	3.7585043	19804.2779
Atlantic	(27.25,27.5]	8	cstar_tref ~ sal + temp + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	2025	2.1391730	8844.3992	17.846793	1982-1999 –> 2000-2012	3.7585043	19804.2779
Atlantic	(27.25,27.5]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate_star + phosphate + phosphate_star	2025	2.1391730	8844.3992	17.846793	1982-1999 –> 2000-2012	3.7585043	19804.2779
Atlantic	(27.25,27.5]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	2025	2.1391730	8844.3992	17.846793	1982-1999 –> 2000-2012	3.7585043	19804.2779
Atlantic	(27.25,27.5]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	2025	2.1391730	8844.3992	17.846793	1982-1999 –> 2000-2012	3.7585043	19804.2779
Atlantic	(27.25,27.5]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	2025	2.1391730	8844.3992	17.846793	1982-1999 –> 2000-2012	3.7585043	19804.2779
Atlantic	(27.25,27.5]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate_star	2025	2.1391730	8844.3992	17.846793	1982-1999 –> 2000-2012	3.7585043	19804.2779
Atlantic	(27.25,27.5]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	2025	2.1391730	8844.3992	17.846793	1982-1999 –> 2000-2012	3.7585043	19804.2779
Atlantic	(27.25,27.5]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	2025	2.1391730	8844.3992	17.846793	1982-1999 –> 2000-2012	3.7585043	19804.2779
Atlantic	(27.25,27.5]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate_star	2025	2.1391730	8844.3992	17.846793	1982-1999 –> 2000-2012	3.7585043	19804.2779
Atlantic	(27.25,27.5]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate + phosphate_star	1109	2.7053602	5372.6371	20.139108	2000-2012 –> 2013-2019	4.8445332	14217.0362
Atlantic	(27.25,27.5]	8	cstar_tref ~ sal + temp + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	1109	2.7053602	5372.6371	20.139108	2000-2012 –> 2013-2019	4.8445332	14217.0362
Atlantic	(27.25,27.5]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate_star + phosphate + phosphate_star	1109	2.7053602	5372.6371	20.139108	2000-2012 –> 2013-2019	4.8445332	14217.0362
Atlantic	(27.25,27.5]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	1109	2.7053602	5372.6371	20.139108	2000-2012 –> 2013-2019	4.8445332	14217.0362
Atlantic	(27.25,27.5]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	1109	2.7053602	5372.6371	20.139108	2000-2012 –> 2013-2019	4.8445332	14217.0362
Atlantic	(27.25,27.5]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	1109	2.7053602	5372.6371	20.139108	2000-2012 –> 2013-2019	4.8445332	14217.0362
Atlantic	(27.25,27.5]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate_star	1109	2.7053602	5372.6371	20.139108	2000-2012 –> 2013-2019	4.8445332	14217.0362
Atlantic	(27.25,27.5]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	1109	2.7053602	5372.6371	20.139108	2000-2012 –> 2013-2019	4.8445332	14217.0362
Atlantic	(27.25,27.5]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	1109	2.7053602	5372.6371	20.139108	2000-2012 –> 2013-2019	4.8445332	14217.0362
Atlantic	(27.25,27.5]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate_star	1109	2.7053602	5372.6371	20.139108	2000-2012 –> 2013-2019	4.8445332	14217.0362
Atlantic	(27.5,27.75]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate + phosphate_star	2804	2.3645727	12801.6371	15.906935	1982-1999 –> 2000-2012	4.4225209	28805.9963
Atlantic	(27.5,27.75]	8	cstar_tref ~ sal + temp + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	2804	2.3645727	12801.6371	15.906935	1982-1999 –> 2000-2012	4.4225209	28805.9963
Atlantic	(27.5,27.75]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate_star + phosphate + phosphate_star	2804	2.3645727	12801.6371	15.906935	1982-1999 –> 2000-2012	4.4225209	28805.9963
Atlantic	(27.5,27.75]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	2804	2.3645727	12801.6371	15.906935	1982-1999 –> 2000-2012	4.4225209	28805.9963
Atlantic	(27.5,27.75]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	2804	2.3645727	12801.6371	15.906935	1982-1999 –> 2000-2012	4.4225209	28805.9963
Atlantic	(27.5,27.75]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	2804	2.3645727	12801.6371	15.906935	1982-1999 –> 2000-2012	4.4225209	28805.9963
Atlantic	(27.5,27.75]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	2804	2.3645727	12801.6371	15.906935	1982-1999 –> 2000-2012	4.4225209	28805.9963
Atlantic	(27.5,27.75]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	2804	2.3645727	12801.6371	15.906935	1982-1999 –> 2000-2012	4.4225209	28805.9963
Atlantic	(27.5,27.75]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	2804	2.3645727	12801.6371	15.906935	1982-1999 –> 2000-2012	4.4225209	28805.9963
Atlantic	(27.5,27.75]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	2804	2.3645727	12801.6371	15.906935	1982-1999 –> 2000-2012	4.4225209	28805.9963
Atlantic	(27.5,27.75]	7	cstar_tref ~ sal + temp + aou + silicate + silicate_star + phosphate + phosphate_star	1515	2.6919437	7317.8822	19.885854	2000-2012 –> 2013-2019	5.0565164	20119.5192
Atlantic	(27.5,27.75]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	1515	2.6919437	7317.8822	19.885854	2000-2012 –> 2013-2019	5.0565164	20119.5192
Atlantic	(27.5,27.75]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	1515	2.6919437	7317.8822	19.885854	2000-2012 –> 2013-2019	5.0565164	20119.5192
Atlantic	(27.5,27.75]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	1515	2.6919437	7317.8822	19.885854	2000-2012 –> 2013-2019	5.0565164	20119.5192
Atlantic	(27.5,27.75]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	1515	2.6919437	7317.8822	19.885854	2000-2012 –> 2013-2019	5.0565164	20119.5192
Atlantic	(27.5,27.75]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	1515	2.6919437	7317.8822	19.885854	2000-2012 –> 2013-2019	5.0565164	20119.5192
Atlantic	(27.5,27.75]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	1515	2.6919437	7317.8822	19.885854	2000-2012 –> 2013-2019	5.0565164	20119.5192
Atlantic	(27.5,27.75]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	1515	2.6919437	7317.8822	19.885854	2000-2012 –> 2013-2019	5.0565164	20119.5192
Atlantic	(27.5,27.75]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	1515	2.6919437	7317.8822	19.885854	2000-2012 –> 2013-2019	5.0565164	20119.5192
Atlantic	(27.5,27.75]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	1515	2.6919437	7317.8822	19.885854	2000-2012 –> 2013-2019	5.0565164	20119.5192
Atlantic	(27.75,27.85]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate + phosphate_star	1036	1.4819371	3773.0620	12.346257	1982-1999 –> 2000-2012	2.9888702	8619.6617
Atlantic	(27.75,27.85]	8	cstar_tref ~ sal + temp + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	1036	1.4819371	3773.0620	12.346257	1982-1999 –> 2000-2012	2.9888702	8619.6617
Atlantic	(27.75,27.85]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate_star + phosphate + phosphate_star	1036	1.4819371	3773.0620	12.346257	1982-1999 –> 2000-2012	2.9888702	8619.6617
Atlantic	(27.75,27.85]	7	cstar_tref ~ sal + temp + aou + silicate + silicate_star + phosphate + phosphate_star	1036	1.4819371	3773.0620	12.346257	1982-1999 –> 2000-2012	2.9888702	8619.6617
Atlantic	(27.75,27.85]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	1036	1.4819371	3773.0620	12.346257	1982-1999 –> 2000-2012	2.9888702	8619.6617
Atlantic	(27.75,27.85]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	1036	1.4819371	3773.0620	12.346257	1982-1999 –> 2000-2012	2.9888702	8619.6617
Atlantic	(27.75,27.85]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	1036	1.4819371	3773.0620	12.346257	1982-1999 –> 2000-2012	2.9888702	8619.6617
Atlantic	(27.75,27.85]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	1036	1.4819371	3773.0620	12.346257	1982-1999 –> 2000-2012	2.9888702	8619.6617
Atlantic	(27.75,27.85]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	1036	1.4819371	3773.0620	12.346257	1982-1999 –> 2000-2012	2.9888702	8619.6617
Atlantic	(27.75,27.85]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	1036	1.4819371	3773.0620	12.346257	1982-1999 –> 2000-2012	2.9888702	8619.6617
Atlantic	(27.75,27.85]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	513	1.6130343	1964.3791	10.721751	2000-2012 –> 2013-2019	3.0949714	5737.4410
Atlantic	(27.75,27.85]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	513	1.6130343	1964.3791	10.721751	2000-2012 –> 2013-2019	3.0949714	5737.4410
Atlantic	(27.75,27.85]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	513	1.6130343	1964.3791	10.721751	2000-2012 –> 2013-2019	3.0949714	5737.4410
Atlantic	(27.75,27.85]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	513	1.6130343	1964.3791	10.721751	2000-2012 –> 2013-2019	3.0949714	5737.4410
Atlantic	(27.75,27.85]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	513	1.6130343	1964.3791	10.721751	2000-2012 –> 2013-2019	3.0949714	5737.4410
Atlantic	(27.75,27.85]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	513	1.6130343	1964.3791	10.721751	2000-2012 –> 2013-2019	3.0949714	5737.4410
Atlantic	(27.75,27.85]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	513	1.6130343	1964.3791	10.721751	2000-2012 –> 2013-2019	3.0949714	5737.4410
Atlantic	(27.75,27.85]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	513	1.6130343	1964.3791	10.721751	2000-2012 –> 2013-2019	3.0949714	5737.4410
Atlantic	(27.75,27.85]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	513	1.6130343	1964.3791	10.721751	2000-2012 –> 2013-2019	3.0949714	5737.4410
Atlantic	(27.75,27.85]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	513	1.6130343	1964.3791	10.721751	2000-2012 –> 2013-2019	3.0949714	5737.4410
Atlantic	(27.85,27.95]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate + phosphate_star	983	2.9167184	4912.1558	23.324808	1982-1999 –> 2000-2012	5.5922061	11351.1595
Atlantic	(27.85,27.95]	8	cstar_tref ~ sal + temp + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	983	2.9167184	4912.1558	23.324808	1982-1999 –> 2000-2012	5.5922061	11351.1595
Atlantic	(27.85,27.95]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate_star + phosphate + phosphate_star	983	2.9167184	4912.1558	23.324808	1982-1999 –> 2000-2012	5.5922061	11351.1595
Atlantic	(27.85,27.95]	7	cstar_tref ~ sal + temp + aou + silicate + silicate_star + phosphate + phosphate_star	983	2.9167184	4912.1558	23.324808	1982-1999 –> 2000-2012	5.5922061	11351.1595
Atlantic	(27.85,27.95]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	983	2.9167184	4912.1558	23.324808	1982-1999 –> 2000-2012	5.5922061	11351.1595
Atlantic	(27.85,27.95]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	983	2.9167184	4912.1558	23.324808	1982-1999 –> 2000-2012	5.5922061	11351.1595
Atlantic	(27.85,27.95]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	983	2.9167184	4912.1558	23.324808	1982-1999 –> 2000-2012	5.5922061	11351.1595
Atlantic	(27.85,27.95]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	983	2.9167184	4912.1558	23.324808	1982-1999 –> 2000-2012	5.5922061	11351.1595
Atlantic	(27.85,27.95]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	983	2.9167184	4912.1558	23.324808	1982-1999 –> 2000-2012	5.5922061	11351.1595
Atlantic	(27.85,27.95]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	983	2.9167184	4912.1558	23.324808	1982-1999 –> 2000-2012	5.5922061	11351.1595
Atlantic	(27.85,27.95]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate + phosphate_star	569	3.1989479	2956.0435	29.016215	2000-2012 –> 2013-2019	6.1156663	7868.1993
Atlantic	(27.85,27.95]	8	cstar_tref ~ sal + temp + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	569	3.1989479	2956.0435	29.016215	2000-2012 –> 2013-2019	6.1156663	7868.1993
Atlantic	(27.85,27.95]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate_star + phosphate + phosphate_star	569	3.1989479	2956.0435	29.016215	2000-2012 –> 2013-2019	6.1156663	7868.1993
Atlantic	(27.85,27.95]	7	cstar_tref ~ sal + temp + aou + silicate + silicate_star + phosphate + phosphate_star	569	3.1989479	2956.0435	29.016215	2000-2012 –> 2013-2019	6.1156663	7868.1993
Atlantic	(27.85,27.95]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	569	3.1989479	2956.0435	29.016215	2000-2012 –> 2013-2019	6.1156663	7868.1993
Atlantic	(27.85,27.95]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate_star	569	3.1989479	2956.0435	29.016215	2000-2012 –> 2013-2019	6.1156663	7868.1993
Atlantic	(27.85,27.95]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate_star	569	3.1989479	2956.0435	29.016215	2000-2012 –> 2013-2019	6.1156663	7868.1993
Atlantic	(27.85,27.95]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	569	3.1989479	2956.0435	29.016215	2000-2012 –> 2013-2019	6.1156663	7868.1993
Atlantic	(27.85,27.95]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	569	3.1989479	2956.0435	29.016215	2000-2012 –> 2013-2019	6.1156663	7868.1993
Atlantic	(27.85,27.95]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	569	3.1989479	2956.0435	29.016215	2000-2012 –> 2013-2019	6.1156663	7868.1993
Atlantic	(27.95,28.05]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate + phosphate_star	1204	3.5310176	6472.7033	22.926889	1982-1999 –> 2000-2012	6.7144051	15128.4199
Atlantic	(27.95,28.05]	8	cstar_tref ~ sal + temp + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	1204	3.5310176	6472.7033	22.926889	1982-1999 –> 2000-2012	6.7144051	15128.4199
Atlantic	(27.95,28.05]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate_star + phosphate + phosphate_star	1204	3.5310176	6472.7033	22.926889	1982-1999 –> 2000-2012	6.7144051	15128.4199
Atlantic	(27.95,28.05]	7	cstar_tref ~ sal + temp + aou + silicate + silicate_star + phosphate + phosphate_star	1204	3.5310176	6472.7033	22.926889	1982-1999 –> 2000-2012	6.7144051	15128.4199
Atlantic	(27.95,28.05]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	1204	3.5310176	6472.7033	22.926889	1982-1999 –> 2000-2012	6.7144051	15128.4199
Atlantic	(27.95,28.05]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	1204	3.5310176	6472.7033	22.926889	1982-1999 –> 2000-2012	6.7144051	15128.4199
Atlantic	(27.95,28.05]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate_star	1204	3.5310176	6472.7033	22.926889	1982-1999 –> 2000-2012	6.7144051	15128.4199
Atlantic	(27.95,28.05]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	1204	3.5310176	6472.7033	22.926889	1982-1999 –> 2000-2012	6.7144051	15128.4199
Atlantic	(27.95,28.05]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	1204	3.5310176	6472.7033	22.926889	1982-1999 –> 2000-2012	6.7144051	15128.4199
Atlantic	(27.95,28.05]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate_star	1204	3.5310176	6472.7033	22.926889	1982-1999 –> 2000-2012	6.7144051	15128.4199
Atlantic	(27.95,28.05]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate + phosphate_star	666	3.1959182	3455.6429	27.018389	2000-2012 –> 2013-2019	6.7269357	9928.3462
Atlantic	(27.95,28.05]	8	cstar_tref ~ sal + temp + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	666	3.1959182	3455.6429	27.018389	2000-2012 –> 2013-2019	6.7269357	9928.3462
Atlantic	(27.95,28.05]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate_star + phosphate + phosphate_star	666	3.1959182	3455.6429	27.018389	2000-2012 –> 2013-2019	6.7269357	9928.3462
Atlantic	(27.95,28.05]	7	cstar_tref ~ sal + temp + aou + silicate + silicate_star + phosphate + phosphate_star	666	3.1959182	3455.6429	27.018389	2000-2012 –> 2013-2019	6.7269357	9928.3462
Atlantic	(27.95,28.05]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	666	3.1959182	3455.6429	27.018389	2000-2012 –> 2013-2019	6.7269357	9928.3462
Atlantic	(27.95,28.05]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	666	3.1959182	3455.6429	27.018389	2000-2012 –> 2013-2019	6.7269357	9928.3462
Atlantic	(27.95,28.05]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate_star	666	3.1959182	3455.6429	27.018389	2000-2012 –> 2013-2019	6.7269357	9928.3462
Atlantic	(27.95,28.05]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	666	3.1959182	3455.6429	27.018389	2000-2012 –> 2013-2019	6.7269357	9928.3462
Atlantic	(27.95,28.05]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	666	3.1959182	3455.6429	27.018389	2000-2012 –> 2013-2019	6.7269357	9928.3462
Atlantic	(27.95,28.05]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate_star	666	3.1959182	3455.6429	27.018389	2000-2012 –> 2013-2019	6.7269357	9928.3462
Atlantic	(28.05,28.1]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	804	0.9868544	2278.3749	8.854339	1982-1999 –> 2000-2012	1.8168843	4949.0267
Atlantic	(28.05,28.1]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	804	0.9868544	2278.3749	8.854339	1982-1999 –> 2000-2012	1.8168843	4949.0267
Atlantic	(28.05,28.1]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	804	0.9868544	2278.3749	8.854339	1982-1999 –> 2000-2012	1.8168843	4949.0267
Atlantic	(28.05,28.1]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	804	0.9868544	2278.3749	8.854339	1982-1999 –> 2000-2012	1.8168843	4949.0267
Atlantic	(28.05,28.1]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate_star	804	0.9868544	2278.3749	8.854339	1982-1999 –> 2000-2012	1.8168843	4949.0267
Atlantic	(28.05,28.1]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	804	0.9868544	2278.3749	8.854339	1982-1999 –> 2000-2012	1.8168843	4949.0267
Atlantic	(28.05,28.1]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate_star	804	0.9868544	2278.3749	8.854339	1982-1999 –> 2000-2012	1.8168843	4949.0267
Atlantic	(28.05,28.1]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	804	0.9868544	2278.3749	8.854339	1982-1999 –> 2000-2012	1.8168843	4949.0267
Atlantic	(28.05,28.1]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	804	0.9868544	2278.3749	8.854339	1982-1999 –> 2000-2012	1.8168843	4949.0267
Atlantic	(28.05,28.1]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	804	0.9868544	2278.3749	8.854339	1982-1999 –> 2000-2012	1.8168843	4949.0267
Atlantic	(28.05,28.1]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	434	1.0166672	1263.9865	9.190678	2000-2012 –> 2013-2019	2.0035216	3542.3614
Atlantic	(28.05,28.1]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	434	1.0166672	1263.9865	9.190678	2000-2012 –> 2013-2019	2.0035216	3542.3614
Atlantic	(28.05,28.1]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	434	1.0166672	1263.9865	9.190678	2000-2012 –> 2013-2019	2.0035216	3542.3614
Atlantic	(28.05,28.1]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	434	1.0166672	1263.9865	9.190678	2000-2012 –> 2013-2019	2.0035216	3542.3614
Atlantic	(28.05,28.1]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	434	1.0166672	1263.9865	9.190678	2000-2012 –> 2013-2019	2.0035216	3542.3614
Atlantic	(28.05,28.1]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate_star	434	1.0166672	1263.9865	9.190678	2000-2012 –> 2013-2019	2.0035216	3542.3614
Atlantic	(28.05,28.1]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	434	1.0166672	1263.9865	9.190678	2000-2012 –> 2013-2019	2.0035216	3542.3614
Atlantic	(28.05,28.1]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	434	1.0166672	1263.9865	9.190678	2000-2012 –> 2013-2019	2.0035216	3542.3614
Atlantic	(28.05,28.1]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate_star	434	1.0166672	1263.9865	9.190678	2000-2012 –> 2013-2019	2.0035216	3542.3614
Atlantic	(28.05,28.1]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	434	1.0166672	1263.9865	9.190678	2000-2012 –> 2013-2019	2.0035216	3542.3614
Atlantic	(28.1,28.15]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate + phosphate_star	936	0.7332530	2093.4379	7.182115	1982-1999 –> 2000-2012	1.2960746	4297.7757
Atlantic	(28.1,28.15]	8	cstar_tref ~ sal + temp + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	936	0.7332530	2093.4379	7.182115	1982-1999 –> 2000-2012	1.2960746	4297.7757
Atlantic	(28.1,28.15]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate_star + phosphate + phosphate_star	936	0.7332530	2093.4379	7.182115	1982-1999 –> 2000-2012	1.2960746	4297.7757
Atlantic	(28.1,28.15]	7	cstar_tref ~ sal + temp + aou + silicate + silicate_star + phosphate + phosphate_star	936	0.7332530	2093.4379	7.182115	1982-1999 –> 2000-2012	1.2960746	4297.7757
Atlantic	(28.1,28.15]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	936	0.7332530	2093.4379	7.182115	1982-1999 –> 2000-2012	1.2960746	4297.7757
Atlantic	(28.1,28.15]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	936	0.7332530	2093.4379	7.182115	1982-1999 –> 2000-2012	1.2960746	4297.7757
Atlantic	(28.1,28.15]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate_star	936	0.7332530	2093.4379	7.182115	1982-1999 –> 2000-2012	1.2960746	4297.7757
Atlantic	(28.1,28.15]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	936	0.7332530	2093.4379	7.182115	1982-1999 –> 2000-2012	1.2960746	4297.7757
Atlantic	(28.1,28.15]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	936	0.7332530	2093.4379	7.182115	1982-1999 –> 2000-2012	1.2960746	4297.7757
Atlantic	(28.1,28.15]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	936	0.7332530	2093.4379	7.182115	1982-1999 –> 2000-2012	1.2960746	4297.7757
Atlantic	(28.1,28.15]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	497	0.8098062	1218.7303	8.335684	2000-2012 –> 2013-2019	1.5430592	3312.1683
Atlantic	(28.1,28.15]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	497	0.8098062	1218.7303	8.335684	2000-2012 –> 2013-2019	1.5430592	3312.1683
Atlantic	(28.1,28.15]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	497	0.8098062	1218.7303	8.335684	2000-2012 –> 2013-2019	1.5430592	3312.1683
Atlantic	(28.1,28.15]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate_star	497	0.8098062	1218.7303	8.335684	2000-2012 –> 2013-2019	1.5430592	3312.1683
Atlantic	(28.1,28.15]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	497	0.8098062	1218.7303	8.335684	2000-2012 –> 2013-2019	1.5430592	3312.1683
Atlantic	(28.1,28.15]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	497	0.8098062	1218.7303	8.335684	2000-2012 –> 2013-2019	1.5430592	3312.1683
Atlantic	(28.1,28.15]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate_star	497	0.8098062	1218.7303	8.335684	2000-2012 –> 2013-2019	1.5430592	3312.1683
Atlantic	(28.1,28.15]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	497	0.8098062	1218.7303	8.335684	2000-2012 –> 2013-2019	1.5430592	3312.1683
Atlantic	(28.1,28.15]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	497	0.8098062	1218.7303	8.335684	2000-2012 –> 2013-2019	1.5430592	3312.1683
Atlantic	(28.1,28.15]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	497	0.8098062	1218.7303	8.335684	2000-2012 –> 2013-2019	1.5430592	3312.1683
Atlantic	(28.15,28.2]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	1540	0.3383167	1050.3104	1.617842	1982-1999 –> 2000-2012	0.7244768	3147.4620
Atlantic	(28.15,28.2]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	1540	0.3383167	1050.3104	1.617842	1982-1999 –> 2000-2012	0.7244768	3147.4620
Atlantic	(28.15,28.2]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	1540	0.3383167	1050.3104	1.617842	1982-1999 –> 2000-2012	0.7244768	3147.4620
Atlantic	(28.15,28.2]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	1540	0.3383167	1050.3104	1.617842	1982-1999 –> 2000-2012	0.7244768	3147.4620
Atlantic	(28.15,28.2]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	1540	0.3383167	1050.3104	1.617842	1982-1999 –> 2000-2012	0.7244768	3147.4620
Atlantic	(28.15,28.2]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	1540	0.3383167	1050.3104	1.617842	1982-1999 –> 2000-2012	0.7244768	3147.4620
Atlantic	(28.15,28.2]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate_star	1540	0.3383167	1050.3104	1.617842	1982-1999 –> 2000-2012	0.7244768	3147.4620
Atlantic	(28.15,28.2]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	1540	0.3383167	1050.3104	1.617842	1982-1999 –> 2000-2012	0.7244768	3147.4620
Atlantic	(28.15,28.2]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	1540	0.3383167	1050.3104	1.617842	1982-1999 –> 2000-2012	0.7244768	3147.4620
Atlantic	(28.15,28.2]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	1540	0.3383167	1050.3104	1.617842	1982-1999 –> 2000-2012	0.7244768	3147.4620
Atlantic	(28.15,28.2]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	819	0.5837815	1460.6029	5.171558	2000-2012 –> 2013-2019	0.9220982	2510.9133
Atlantic	(28.15,28.2]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	819	0.5837815	1460.6029	5.171558	2000-2012 –> 2013-2019	0.9220982	2510.9133
Atlantic	(28.15,28.2]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	819	0.5837815	1460.6029	5.171558	2000-2012 –> 2013-2019	0.9220982	2510.9133
Atlantic	(28.15,28.2]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	819	0.5837815	1460.6029	5.171558	2000-2012 –> 2013-2019	0.9220982	2510.9133
Atlantic	(28.15,28.2]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	819	0.5837815	1460.6029	5.171558	2000-2012 –> 2013-2019	0.9220982	2510.9133
Atlantic	(28.15,28.2]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate_star	819	0.5837815	1460.6029	5.171558	2000-2012 –> 2013-2019	0.9220982	2510.9133
Atlantic	(28.15,28.2]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	819	0.5837815	1460.6029	5.171558	2000-2012 –> 2013-2019	0.9220982	2510.9133
Atlantic	(28.15,28.2]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	819	0.5837815	1460.6029	5.171558	2000-2012 –> 2013-2019	0.9220982	2510.9133
Atlantic	(28.15,28.2]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate_star	819	0.5837815	1460.6029	5.171558	2000-2012 –> 2013-2019	0.9220982	2510.9133
Atlantic	(28.15,28.2]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	819	0.5837815	1460.6029	5.171558	2000-2012 –> 2013-2019	0.9220982	2510.9133
Atlantic	(28.2, Inf]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate + phosphate_star	2693	0.3159903	1455.4947	2.104197	1982-1999 –> 2000-2012	0.5862118	2293.0711
Atlantic	(28.2, Inf]	8	cstar_tref ~ sal + temp + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	2693	0.3159903	1455.4947	2.104197	1982-1999 –> 2000-2012	0.5862118	2293.0711
Atlantic	(28.2, Inf]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate_star + phosphate + phosphate_star	2693	0.3159903	1455.4947	2.104197	1982-1999 –> 2000-2012	0.5862118	2293.0711
Atlantic	(28.2, Inf]	7	cstar_tref ~ sal + temp + aou + silicate + silicate_star + phosphate + phosphate_star	2693	0.3159903	1455.4947	2.104197	1982-1999 –> 2000-2012	0.5862118	2293.0711
Atlantic	(28.2, Inf]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	2693	0.3159903	1455.4947	2.104197	1982-1999 –> 2000-2012	0.5862118	2293.0711
Atlantic	(28.2, Inf]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	2693	0.3159903	1455.4947	2.104197	1982-1999 –> 2000-2012	0.5862118	2293.0711
Atlantic	(28.2, Inf]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	2693	0.3159903	1455.4947	2.104197	1982-1999 –> 2000-2012	0.5862118	2293.0711
Atlantic	(28.2, Inf]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	2693	0.3159903	1455.4947	2.104197	1982-1999 –> 2000-2012	0.5862118	2293.0711
Atlantic	(28.2, Inf]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	2693	0.3159903	1455.4947	2.104197	1982-1999 –> 2000-2012	0.5862118	2293.0711
Atlantic	(28.2, Inf]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	2693	0.3159903	1455.4947	2.104197	1982-1999 –> 2000-2012	0.5862118	2293.0711
Atlantic	(28.2, Inf]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	1557	0.3330850	1013.1749	2.720445	2000-2012 –> 2013-2019	0.6490752	2468.6696
Atlantic	(28.2, Inf]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	1557	0.3330850	1013.1749	2.720445	2000-2012 –> 2013-2019	0.6490752	2468.6696
Atlantic	(28.2, Inf]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	1557	0.3330850	1013.1749	2.720445	2000-2012 –> 2013-2019	0.6490752	2468.6696
Atlantic	(28.2, Inf]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	1557	0.3330850	1013.1749	2.720445	2000-2012 –> 2013-2019	0.6490752	2468.6696
Atlantic	(28.2, Inf]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	1557	0.3330850	1013.1749	2.720445	2000-2012 –> 2013-2019	0.6490752	2468.6696
Atlantic	(28.2, Inf]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate_star	1557	0.3330850	1013.1749	2.720445	2000-2012 –> 2013-2019	0.6490752	2468.6696
Atlantic	(28.2, Inf]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	1557	0.3330850	1013.1749	2.720445	2000-2012 –> 2013-2019	0.6490752	2468.6696
Atlantic	(28.2, Inf]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	1557	0.3330850	1013.1749	2.720445	2000-2012 –> 2013-2019	0.6490752	2468.6696
Atlantic	(28.2, Inf]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate_star	1557	0.3330850	1013.1749	2.720445	2000-2012 –> 2013-2019	0.6490752	2468.6696
Atlantic	(28.2, Inf]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	1557	0.3330850	1013.1749	2.720445	2000-2012 –> 2013-2019	0.6490752	2468.6696
Indo-Pacific	(-Inf,26]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate + phosphate_star	5288	5.8703920	33743.3827	32.845116	1982-1999 –> 2000-2012	12.0393795	80350.0591
Indo-Pacific	(-Inf,26]	8	cstar_tref ~ sal + temp + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	5288	5.8703920	33743.3827	32.845116	1982-1999 –> 2000-2012	12.0393795	80350.0591
Indo-Pacific	(-Inf,26]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate_star + phosphate + phosphate_star	5288	5.8703920	33743.3827	32.845116	1982-1999 –> 2000-2012	12.0393795	80350.0591
Indo-Pacific	(-Inf,26]	7	cstar_tref ~ sal + temp + aou + silicate + silicate_star + phosphate + phosphate_star	5288	5.8703920	33743.3827	32.845116	1982-1999 –> 2000-2012	12.0393795	80350.0591
Indo-Pacific	(-Inf,26]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	5288	5.8703920	33743.3827	32.845116	1982-1999 –> 2000-2012	12.0393795	80350.0591
Indo-Pacific	(-Inf,26]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate_star	5288	5.8703920	33743.3827	32.845116	1982-1999 –> 2000-2012	12.0393795	80350.0591
Indo-Pacific	(-Inf,26]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate_star	5288	5.8703920	33743.3827	32.845116	1982-1999 –> 2000-2012	12.0393795	80350.0591
Indo-Pacific	(-Inf,26]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	5288	5.8703920	33743.3827	32.845116	1982-1999 –> 2000-2012	12.0393795	80350.0591
Indo-Pacific	(-Inf,26]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	5288	5.8703920	33743.3827	32.845116	1982-1999 –> 2000-2012	12.0393795	80350.0591
Indo-Pacific	(-Inf,26]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	5288	5.8703920	33743.3827	32.845116	1982-1999 –> 2000-2012	12.0393795	80350.0591
Indo-Pacific	(-Inf,26]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate + phosphate_star	2816	5.5768716	17688.7746	25.837968	2000-2012 –> 2013-2019	11.4472636	51432.1573
Indo-Pacific	(-Inf,26]	8	cstar_tref ~ sal + temp + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	2816	5.5768716	17688.7746	25.837968	2000-2012 –> 2013-2019	11.4472636	51432.1573
Indo-Pacific	(-Inf,26]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate_star + phosphate + phosphate_star	2816	5.5768716	17688.7746	25.837968	2000-2012 –> 2013-2019	11.4472636	51432.1573
Indo-Pacific	(-Inf,26]	7	cstar_tref ~ sal + temp + aou + silicate + silicate_star + phosphate + phosphate_star	2816	5.5768716	17688.7746	25.837968	2000-2012 –> 2013-2019	11.4472636	51432.1573
Indo-Pacific	(-Inf,26]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	2816	5.5768716	17688.7746	25.837968	2000-2012 –> 2013-2019	11.4472636	51432.1573
Indo-Pacific	(-Inf,26]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	2816	5.5768716	17688.7746	25.837968	2000-2012 –> 2013-2019	11.4472636	51432.1573
Indo-Pacific	(-Inf,26]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	2816	5.5768716	17688.7746	25.837968	2000-2012 –> 2013-2019	11.4472636	51432.1573
Indo-Pacific	(-Inf,26]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	2816	5.5768716	17688.7746	25.837968	2000-2012 –> 2013-2019	11.4472636	51432.1573
Indo-Pacific	(-Inf,26]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	2816	5.5768716	17688.7746	25.837968	2000-2012 –> 2013-2019	11.4472636	51432.1573
Indo-Pacific	(-Inf,26]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	2816	5.5768716	17688.7746	25.837968	2000-2012 –> 2013-2019	11.4472636	51432.1573
Indo-Pacific	(26,26.5]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate + phosphate_star	5525	4.3695428	31992.2458	49.315429	1982-1999 –> 2000-2012	8.3981502	75388.0758
Indo-Pacific	(26,26.5]	8	cstar_tref ~ sal + temp + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	5525	4.3695428	31992.2458	49.315429	1982-1999 –> 2000-2012	8.3981502	75388.0758
Indo-Pacific	(26,26.5]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate_star + phosphate + phosphate_star	5525	4.3695428	31992.2458	49.315429	1982-1999 –> 2000-2012	8.3981502	75388.0758
Indo-Pacific	(26,26.5]	7	cstar_tref ~ sal + temp + aou + silicate + silicate_star + phosphate + phosphate_star	5525	4.3695428	31992.2458	49.315429	1982-1999 –> 2000-2012	8.3981502	75388.0758
Indo-Pacific	(26,26.5]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	5525	4.3695428	31992.2458	49.315429	1982-1999 –> 2000-2012	8.3981502	75388.0758
Indo-Pacific	(26,26.5]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	5525	4.3695428	31992.2458	49.315429	1982-1999 –> 2000-2012	8.3981502	75388.0758
Indo-Pacific	(26,26.5]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	5525	4.3695428	31992.2458	49.315429	1982-1999 –> 2000-2012	8.3981502	75388.0758
Indo-Pacific	(26,26.5]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	5525	4.3695428	31992.2458	49.315429	1982-1999 –> 2000-2012	8.3981502	75388.0758
Indo-Pacific	(26,26.5]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	5525	4.3695428	31992.2458	49.315429	1982-1999 –> 2000-2012	8.3981502	75388.0758
Indo-Pacific	(26,26.5]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	5525	4.3695428	31992.2458	49.315429	1982-1999 –> 2000-2012	8.3981502	75388.0758
Indo-Pacific	(26,26.5]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	2918	4.2073013	16684.2151	55.736884	2000-2012 –> 2013-2019	8.5768441	48676.4609
Indo-Pacific	(26,26.5]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	2918	4.2073013	16684.2151	55.736884	2000-2012 –> 2013-2019	8.5768441	48676.4609
Indo-Pacific	(26,26.5]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	2918	4.2073013	16684.2151	55.736884	2000-2012 –> 2013-2019	8.5768441	48676.4609
Indo-Pacific	(26,26.5]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	2918	4.2073013	16684.2151	55.736884	2000-2012 –> 2013-2019	8.5768441	48676.4609
Indo-Pacific	(26,26.5]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate_star	2918	4.2073013	16684.2151	55.736884	2000-2012 –> 2013-2019	8.5768441	48676.4609
Indo-Pacific	(26,26.5]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	2918	4.2073013	16684.2151	55.736884	2000-2012 –> 2013-2019	8.5768441	48676.4609
Indo-Pacific	(26,26.5]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate_star	2918	4.2073013	16684.2151	55.736884	2000-2012 –> 2013-2019	8.5768441	48676.4609
Indo-Pacific	(26,26.5]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	2918	4.2073013	16684.2151	55.736884	2000-2012 –> 2013-2019	8.5768441	48676.4609
Indo-Pacific	(26,26.5]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	2918	4.2073013	16684.2151	55.736884	2000-2012 –> 2013-2019	8.5768441	48676.4609
Indo-Pacific	(26,26.5]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	2918	4.2073013	16684.2151	55.736884	2000-2012 –> 2013-2019	8.5768441	48676.4609
Indo-Pacific	(26.5,26.75]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	4544	3.7883819	25017.9750	25.871154	1982-1999 –> 2000-2012	7.4920814	59346.7074
Indo-Pacific	(26.5,26.75]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	4544	3.7883819	25017.9750	25.871154	1982-1999 –> 2000-2012	7.4920814	59346.7074
Indo-Pacific	(26.5,26.75]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	4544	3.7883819	25017.9750	25.871154	1982-1999 –> 2000-2012	7.4920814	59346.7074
Indo-Pacific	(26.5,26.75]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	4544	3.7883819	25017.9750	25.871154	1982-1999 –> 2000-2012	7.4920814	59346.7074
Indo-Pacific	(26.5,26.75]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	4544	3.7883819	25017.9750	25.871154	1982-1999 –> 2000-2012	7.4920814	59346.7074
Indo-Pacific	(26.5,26.75]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate_star	4544	3.7883819	25017.9750	25.871154	1982-1999 –> 2000-2012	7.4920814	59346.7074
Indo-Pacific	(26.5,26.75]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate_star	4544	3.7883819	25017.9750	25.871154	1982-1999 –> 2000-2012	7.4920814	59346.7074
Indo-Pacific	(26.5,26.75]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	4544	3.7883819	25017.9750	25.871154	1982-1999 –> 2000-2012	7.4920814	59346.7074
Indo-Pacific	(26.5,26.75]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	4544	3.7883819	25017.9750	25.871154	1982-1999 –> 2000-2012	7.4920814	59346.7074
Indo-Pacific	(26.5,26.75]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	4544	3.7883819	25017.9750	25.871154	1982-1999 –> 2000-2012	7.4920814	59346.7074
Indo-Pacific	(26.5,26.75]	7	cstar_tref ~ sal + temp + aou + silicate + silicate_star + phosphate + phosphate_star	2447	4.1092818	13878.7222	29.174056	2000-2012 –> 2013-2019	7.8976637	38896.6972
Indo-Pacific	(26.5,26.75]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	2447	4.1092818	13878.7222	29.174056	2000-2012 –> 2013-2019	7.8976637	38896.6972
Indo-Pacific	(26.5,26.75]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	2447	4.1092818	13878.7222	29.174056	2000-2012 –> 2013-2019	7.8976637	38896.6972
Indo-Pacific	(26.5,26.75]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	2447	4.1092818	13878.7222	29.174056	2000-2012 –> 2013-2019	7.8976637	38896.6972
Indo-Pacific	(26.5,26.75]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	2447	4.1092818	13878.7222	29.174056	2000-2012 –> 2013-2019	7.8976637	38896.6972
Indo-Pacific	(26.5,26.75]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	2447	4.1092818	13878.7222	29.174056	2000-2012 –> 2013-2019	7.8976637	38896.6972
Indo-Pacific	(26.5,26.75]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	2447	4.1092818	13878.7222	29.174056	2000-2012 –> 2013-2019	7.8976637	38896.6972
Indo-Pacific	(26.5,26.75]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	2447	4.1092818	13878.7222	29.174056	2000-2012 –> 2013-2019	7.8976637	38896.6972
Indo-Pacific	(26.5,26.75]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	2447	4.1092818	13878.7222	29.174056	2000-2012 –> 2013-2019	7.8976637	38896.6972
Indo-Pacific	(26.5,26.75]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	2447	4.1092818	13878.7222	29.174056	2000-2012 –> 2013-2019	7.8976637	38896.6972
Indo-Pacific	(26.75,27]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	5354	4.0686072	30238.5379	20.211284	1982-1999 –> 2000-2012	7.5960628	70615.3985
Indo-Pacific	(26.75,27]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	5354	4.0686072	30238.5379	20.211284	1982-1999 –> 2000-2012	7.5960628	70615.3985
Indo-Pacific	(26.75,27]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	5354	4.0686072	30238.5379	20.211284	1982-1999 –> 2000-2012	7.5960628	70615.3985
Indo-Pacific	(26.75,27]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	5354	4.0686072	30238.5379	20.211284	1982-1999 –> 2000-2012	7.5960628	70615.3985
Indo-Pacific	(26.75,27]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	5354	4.0686072	30238.5379	20.211284	1982-1999 –> 2000-2012	7.5960628	70615.3985
Indo-Pacific	(26.75,27]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate_star	5354	4.0686072	30238.5379	20.211284	1982-1999 –> 2000-2012	7.5960628	70615.3985
Indo-Pacific	(26.75,27]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate_star	5354	4.0686072	30238.5379	20.211284	1982-1999 –> 2000-2012	7.5960628	70615.3985
Indo-Pacific	(26.75,27]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	5354	4.0686072	30238.5379	20.211284	1982-1999 –> 2000-2012	7.5960628	70615.3985
Indo-Pacific	(26.75,27]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	5354	4.0686072	30238.5379	20.211284	1982-1999 –> 2000-2012	7.5960628	70615.3985
Indo-Pacific	(26.75,27]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	5354	4.0686072	30238.5379	20.211284	1982-1999 –> 2000-2012	7.5960628	70615.3985
Indo-Pacific	(26.75,27]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	2868	4.4456859	16714.7657	18.413340	2000-2012 –> 2013-2019	8.5142930	46953.3037
Indo-Pacific	(26.75,27]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	2868	4.4456859	16714.7657	18.413340	2000-2012 –> 2013-2019	8.5142930	46953.3037
Indo-Pacific	(26.75,27]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	2868	4.4456859	16714.7657	18.413340	2000-2012 –> 2013-2019	8.5142930	46953.3037
Indo-Pacific	(26.75,27]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	2868	4.4456859	16714.7657	18.413340	2000-2012 –> 2013-2019	8.5142930	46953.3037
Indo-Pacific	(26.75,27]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	2868	4.4456859	16714.7657	18.413340	2000-2012 –> 2013-2019	8.5142930	46953.3037
Indo-Pacific	(26.75,27]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate_star	2868	4.4456859	16714.7657	18.413340	2000-2012 –> 2013-2019	8.5142930	46953.3037
Indo-Pacific	(26.75,27]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate_star	2868	4.4456859	16714.7657	18.413340	2000-2012 –> 2013-2019	8.5142930	46953.3037
Indo-Pacific	(26.75,27]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	2868	4.4456859	16714.7657	18.413340	2000-2012 –> 2013-2019	8.5142930	46953.3037
Indo-Pacific	(26.75,27]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	2868	4.4456859	16714.7657	18.413340	2000-2012 –> 2013-2019	8.5142930	46953.3037
Indo-Pacific	(26.75,27]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	2868	4.4456859	16714.7657	18.413340	2000-2012 –> 2013-2019	8.5142930	46953.3037
Indo-Pacific	(27,27.25]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate + phosphate_star	6466	3.6843471	35232.2482	42.550831	1982-1999 –> 2000-2012	6.7021130	79849.7973
Indo-Pacific	(27,27.25]	8	cstar_tref ~ sal + temp + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	6466	3.6843471	35232.2482	42.550831	1982-1999 –> 2000-2012	6.7021130	79849.7973
Indo-Pacific	(27,27.25]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate_star + phosphate + phosphate_star	6466	3.6843471	35232.2482	42.550831	1982-1999 –> 2000-2012	6.7021130	79849.7973
Indo-Pacific	(27,27.25]	7	cstar_tref ~ sal + temp + aou + silicate + silicate_star + phosphate + phosphate_star	6466	3.6843471	35232.2482	42.550831	1982-1999 –> 2000-2012	6.7021130	79849.7973
Indo-Pacific	(27,27.25]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	6466	3.6843471	35232.2482	42.550831	1982-1999 –> 2000-2012	6.7021130	79849.7973
Indo-Pacific	(27,27.25]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	6466	3.6843471	35232.2482	42.550831	1982-1999 –> 2000-2012	6.7021130	79849.7973
Indo-Pacific	(27,27.25]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	6466	3.6843471	35232.2482	42.550831	1982-1999 –> 2000-2012	6.7021130	79849.7973
Indo-Pacific	(27,27.25]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	6466	3.6843471	35232.2482	42.550831	1982-1999 –> 2000-2012	6.7021130	79849.7973
Indo-Pacific	(27,27.25]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	6466	3.6843471	35232.2482	42.550831	1982-1999 –> 2000-2012	6.7021130	79849.7973
Indo-Pacific	(27,27.25]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	6466	3.6843471	35232.2482	42.550831	1982-1999 –> 2000-2012	6.7021130	79849.7973
Indo-Pacific	(27,27.25]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	3470	4.3044979	19995.4772	31.639625	2000-2012 –> 2013-2019	7.9888450	55227.7255
Indo-Pacific	(27,27.25]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	3470	4.3044979	19995.4772	31.639625	2000-2012 –> 2013-2019	7.9888450	55227.7255
Indo-Pacific	(27,27.25]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	3470	4.3044979	19995.4772	31.639625	2000-2012 –> 2013-2019	7.9888450	55227.7255
Indo-Pacific	(27,27.25]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	3470	4.3044979	19995.4772	31.639625	2000-2012 –> 2013-2019	7.9888450	55227.7255
Indo-Pacific	(27,27.25]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	3470	4.3044979	19995.4772	31.639625	2000-2012 –> 2013-2019	7.9888450	55227.7255
Indo-Pacific	(27,27.25]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate_star	3470	4.3044979	19995.4772	31.639625	2000-2012 –> 2013-2019	7.9888450	55227.7255
Indo-Pacific	(27,27.25]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	3470	4.3044979	19995.4772	31.639625	2000-2012 –> 2013-2019	7.9888450	55227.7255
Indo-Pacific	(27,27.25]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate_star	3470	4.3044979	19995.4772	31.639625	2000-2012 –> 2013-2019	7.9888450	55227.7255
Indo-Pacific	(27,27.25]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	3470	4.3044979	19995.4772	31.639625	2000-2012 –> 2013-2019	7.9888450	55227.7255
Indo-Pacific	(27,27.25]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	3470	4.3044979	19995.4772	31.639625	2000-2012 –> 2013-2019	7.9888450	55227.7255
Indo-Pacific	(27.25,27.5]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate + phosphate_star	5847	3.1694873	30100.9131	32.456526	1982-1999 –> 2000-2012	5.9497930	68655.5399
Indo-Pacific	(27.25,27.5]	8	cstar_tref ~ sal + temp + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	5847	3.1694873	30100.9131	32.456526	1982-1999 –> 2000-2012	5.9497930	68655.5399
Indo-Pacific	(27.25,27.5]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate_star + phosphate + phosphate_star	5847	3.1694873	30100.9131	32.456526	1982-1999 –> 2000-2012	5.9497930	68655.5399
Indo-Pacific	(27.25,27.5]	7	cstar_tref ~ sal + temp + aou + silicate + silicate_star + phosphate + phosphate_star	5847	3.1694873	30100.9131	32.456526	1982-1999 –> 2000-2012	5.9497930	68655.5399
Indo-Pacific	(27.25,27.5]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	5847	3.1694873	30100.9131	32.456526	1982-1999 –> 2000-2012	5.9497930	68655.5399
Indo-Pacific	(27.25,27.5]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	5847	3.1694873	30100.9131	32.456526	1982-1999 –> 2000-2012	5.9497930	68655.5399
Indo-Pacific	(27.25,27.5]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate_star	5847	3.1694873	30100.9131	32.456526	1982-1999 –> 2000-2012	5.9497930	68655.5399
Indo-Pacific	(27.25,27.5]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	5847	3.1694873	30100.9131	32.456526	1982-1999 –> 2000-2012	5.9497930	68655.5399
Indo-Pacific	(27.25,27.5]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	5847	3.1694873	30100.9131	32.456526	1982-1999 –> 2000-2012	5.9497930	68655.5399
Indo-Pacific	(27.25,27.5]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	5847	3.1694873	30100.9131	32.456526	1982-1999 –> 2000-2012	5.9497930	68655.5399
Indo-Pacific	(27.25,27.5]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate + phosphate_star	3127	3.6527302	16993.9415	32.104464	2000-2012 –> 2013-2019	6.8222175	47094.8546
Indo-Pacific	(27.25,27.5]	8	cstar_tref ~ sal + temp + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	3127	3.6527302	16993.9415	32.104464	2000-2012 –> 2013-2019	6.8222175	47094.8546
Indo-Pacific	(27.25,27.5]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate_star + phosphate + phosphate_star	3127	3.6527302	16993.9415	32.104464	2000-2012 –> 2013-2019	6.8222175	47094.8546
Indo-Pacific	(27.25,27.5]	7	cstar_tref ~ sal + temp + aou + silicate + silicate_star + phosphate + phosphate_star	3127	3.6527302	16993.9415	32.104464	2000-2012 –> 2013-2019	6.8222175	47094.8546
Indo-Pacific	(27.25,27.5]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	3127	3.6527302	16993.9415	32.104464	2000-2012 –> 2013-2019	6.8222175	47094.8546
Indo-Pacific	(27.25,27.5]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	3127	3.6527302	16993.9415	32.104464	2000-2012 –> 2013-2019	6.8222175	47094.8546
Indo-Pacific	(27.25,27.5]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate_star	3127	3.6527302	16993.9415	32.104464	2000-2012 –> 2013-2019	6.8222175	47094.8546
Indo-Pacific	(27.25,27.5]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	3127	3.6527302	16993.9415	32.104464	2000-2012 –> 2013-2019	6.8222175	47094.8546
Indo-Pacific	(27.25,27.5]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	3127	3.6527302	16993.9415	32.104464	2000-2012 –> 2013-2019	6.8222175	47094.8546
Indo-Pacific	(27.25,27.5]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	3127	3.6527302	16993.9415	32.104464	2000-2012 –> 2013-2019	6.8222175	47094.8546
Indo-Pacific	(27.5,27.75]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	5836	2.9783170	29318.1854	12.869829	1982-1999 –> 2000-2012	5.1744982	65089.9113
Indo-Pacific	(27.5,27.75]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	5836	2.9783170	29318.1854	12.869829	1982-1999 –> 2000-2012	5.1744982	65089.9113
Indo-Pacific	(27.5,27.75]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	5836	2.9783170	29318.1854	12.869829	1982-1999 –> 2000-2012	5.1744982	65089.9113
Indo-Pacific	(27.5,27.75]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	5836	2.9783170	29318.1854	12.869829	1982-1999 –> 2000-2012	5.1744982	65089.9113
Indo-Pacific	(27.5,27.75]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate_star	5836	2.9783170	29318.1854	12.869829	1982-1999 –> 2000-2012	5.1744982	65089.9113
Indo-Pacific	(27.5,27.75]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	5836	2.9783170	29318.1854	12.869829	1982-1999 –> 2000-2012	5.1744982	65089.9113
Indo-Pacific	(27.5,27.75]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate_star	5836	2.9783170	29318.1854	12.869829	1982-1999 –> 2000-2012	5.1744982	65089.9113
Indo-Pacific	(27.5,27.75]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	5836	2.9783170	29318.1854	12.869829	1982-1999 –> 2000-2012	5.1744982	65089.9113
Indo-Pacific	(27.5,27.75]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	5836	2.9783170	29318.1854	12.869829	1982-1999 –> 2000-2012	5.1744982	65089.9113
Indo-Pacific	(27.5,27.75]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	5836	2.9783170	29318.1854	12.869829	1982-1999 –> 2000-2012	5.1744982	65089.9113
Indo-Pacific	(27.5,27.75]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	3104	3.4602176	16532.9562	13.845770	2000-2012 –> 2013-2019	6.4385346	45851.1416
Indo-Pacific	(27.5,27.75]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	3104	3.4602176	16532.9562	13.845770	2000-2012 –> 2013-2019	6.4385346	45851.1416
Indo-Pacific	(27.5,27.75]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	3104	3.4602176	16532.9562	13.845770	2000-2012 –> 2013-2019	6.4385346	45851.1416
Indo-Pacific	(27.5,27.75]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate_star	3104	3.4602176	16532.9562	13.845770	2000-2012 –> 2013-2019	6.4385346	45851.1416
Indo-Pacific	(27.5,27.75]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	3104	3.4602176	16532.9562	13.845770	2000-2012 –> 2013-2019	6.4385346	45851.1416
Indo-Pacific	(27.5,27.75]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	3104	3.4602176	16532.9562	13.845770	2000-2012 –> 2013-2019	6.4385346	45851.1416
Indo-Pacific	(27.5,27.75]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate_star	3104	3.4602176	16532.9562	13.845770	2000-2012 –> 2013-2019	6.4385346	45851.1416
Indo-Pacific	(27.5,27.75]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	3104	3.4602176	16532.9562	13.845770	2000-2012 –> 2013-2019	6.4385346	45851.1416
Indo-Pacific	(27.5,27.75]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	3104	3.4602176	16532.9562	13.845770	2000-2012 –> 2013-2019	6.4385346	45851.1416
Indo-Pacific	(27.5,27.75]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	3104	3.4602176	16532.9562	13.845770	2000-2012 –> 2013-2019	6.4385346	45851.1416
Indo-Pacific	(27.75,27.85]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate + phosphate_star	2317	2.6289513	11072.5162	21.934772	1982-1999 –> 2000-2012	4.5090898	23929.3942
Indo-Pacific	(27.75,27.85]	8	cstar_tref ~ sal + temp + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	2317	2.6289513	11072.5162	21.934772	1982-1999 –> 2000-2012	4.5090898	23929.3942
Indo-Pacific	(27.75,27.85]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate_star + phosphate + phosphate_star	2317	2.6289513	11072.5162	21.934772	1982-1999 –> 2000-2012	4.5090898	23929.3942
Indo-Pacific	(27.75,27.85]	7	cstar_tref ~ sal + temp + aou + silicate + silicate_star + phosphate + phosphate_star	2317	2.6289513	11072.5162	21.934772	1982-1999 –> 2000-2012	4.5090898	23929.3942
Indo-Pacific	(27.75,27.85]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	2317	2.6289513	11072.5162	21.934772	1982-1999 –> 2000-2012	4.5090898	23929.3942
Indo-Pacific	(27.75,27.85]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	2317	2.6289513	11072.5162	21.934772	1982-1999 –> 2000-2012	4.5090898	23929.3942
Indo-Pacific	(27.75,27.85]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	2317	2.6289513	11072.5162	21.934772	1982-1999 –> 2000-2012	4.5090898	23929.3942
Indo-Pacific	(27.75,27.85]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	2317	2.6289513	11072.5162	21.934772	1982-1999 –> 2000-2012	4.5090898	23929.3942
Indo-Pacific	(27.75,27.85]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	2317	2.6289513	11072.5162	21.934772	1982-1999 –> 2000-2012	4.5090898	23929.3942
Indo-Pacific	(27.75,27.85]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	2317	2.6289513	11072.5162	21.934772	1982-1999 –> 2000-2012	4.5090898	23929.3942
Indo-Pacific	(27.75,27.85]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate + phosphate_star	1264	2.9672019	6354.5784	33.274356	2000-2012 –> 2013-2019	5.5961533	17427.0947
Indo-Pacific	(27.75,27.85]	8	cstar_tref ~ sal + temp + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	1264	2.9672019	6354.5784	33.274356	2000-2012 –> 2013-2019	5.5961533	17427.0947
Indo-Pacific	(27.75,27.85]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate_star + phosphate + phosphate_star	1264	2.9672019	6354.5784	33.274356	2000-2012 –> 2013-2019	5.5961533	17427.0947
Indo-Pacific	(27.75,27.85]	7	cstar_tref ~ sal + temp + aou + silicate + silicate_star + phosphate + phosphate_star	1264	2.9672019	6354.5784	33.274356	2000-2012 –> 2013-2019	5.5961533	17427.0947
Indo-Pacific	(27.75,27.85]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	1264	2.9672019	6354.5784	33.274356	2000-2012 –> 2013-2019	5.5961533	17427.0947
Indo-Pacific	(27.75,27.85]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	1264	2.9672019	6354.5784	33.274356	2000-2012 –> 2013-2019	5.5961533	17427.0947
Indo-Pacific	(27.75,27.85]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	1264	2.9672019	6354.5784	33.274356	2000-2012 –> 2013-2019	5.5961533	17427.0947
Indo-Pacific	(27.75,27.85]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	1264	2.9672019	6354.5784	33.274356	2000-2012 –> 2013-2019	5.5961533	17427.0947
Indo-Pacific	(27.75,27.85]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	1264	2.9672019	6354.5784	33.274356	2000-2012 –> 2013-2019	5.5961533	17427.0947
Indo-Pacific	(27.75,27.85]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	1264	2.9672019	6354.5784	33.274356	2000-2012 –> 2013-2019	5.5961533	17427.0947
Indo-Pacific	(27.85,27.95]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate + phosphate_star	2907	2.7004857	14043.5202	28.455559	1982-1999 –> 2000-2012	4.7932875	31364.2073
Indo-Pacific	(27.85,27.95]	8	cstar_tref ~ sal + temp + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	2907	2.7004857	14043.5202	28.455559	1982-1999 –> 2000-2012	4.7932875	31364.2073
Indo-Pacific	(27.85,27.95]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate_star + phosphate + phosphate_star	2907	2.7004857	14043.5202	28.455559	1982-1999 –> 2000-2012	4.7932875	31364.2073
Indo-Pacific	(27.85,27.95]	7	cstar_tref ~ sal + temp + aou + silicate + silicate_star + phosphate + phosphate_star	2907	2.7004857	14043.5202	28.455559	1982-1999 –> 2000-2012	4.7932875	31364.2073
Indo-Pacific	(27.85,27.95]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	2907	2.7004857	14043.5202	28.455559	1982-1999 –> 2000-2012	4.7932875	31364.2073
Indo-Pacific	(27.85,27.95]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	2907	2.7004857	14043.5202	28.455559	1982-1999 –> 2000-2012	4.7932875	31364.2073
Indo-Pacific	(27.85,27.95]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	2907	2.7004857	14043.5202	28.455559	1982-1999 –> 2000-2012	4.7932875	31364.2073
Indo-Pacific	(27.85,27.95]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	2907	2.7004857	14043.5202	28.455559	1982-1999 –> 2000-2012	4.7932875	31364.2073
Indo-Pacific	(27.85,27.95]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	2907	2.7004857	14043.5202	28.455559	1982-1999 –> 2000-2012	4.7932875	31364.2073
Indo-Pacific	(27.85,27.95]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	2907	2.7004857	14043.5202	28.455559	1982-1999 –> 2000-2012	4.7932875	31364.2073
Indo-Pacific	(27.85,27.95]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate + phosphate_star	1542	3.1760747	7958.0189	31.994368	2000-2012 –> 2013-2019	5.8765604	22001.5392
Indo-Pacific	(27.85,27.95]	8	cstar_tref ~ sal + temp + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	1542	3.1760747	7958.0189	31.994368	2000-2012 –> 2013-2019	5.8765604	22001.5392
Indo-Pacific	(27.85,27.95]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate_star + phosphate + phosphate_star	1542	3.1760747	7958.0189	31.994368	2000-2012 –> 2013-2019	5.8765604	22001.5392
Indo-Pacific	(27.85,27.95]	7	cstar_tref ~ sal + temp + aou + silicate + silicate_star + phosphate + phosphate_star	1542	3.1760747	7958.0189	31.994368	2000-2012 –> 2013-2019	5.8765604	22001.5392
Indo-Pacific	(27.85,27.95]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	1542	3.1760747	7958.0189	31.994368	2000-2012 –> 2013-2019	5.8765604	22001.5392
Indo-Pacific	(27.85,27.95]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	1542	3.1760747	7958.0189	31.994368	2000-2012 –> 2013-2019	5.8765604	22001.5392
Indo-Pacific	(27.85,27.95]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	1542	3.1760747	7958.0189	31.994368	2000-2012 –> 2013-2019	5.8765604	22001.5392
Indo-Pacific	(27.85,27.95]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	1542	3.1760747	7958.0189	31.994368	2000-2012 –> 2013-2019	5.8765604	22001.5392
Indo-Pacific	(27.85,27.95]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	1542	3.1760747	7958.0189	31.994368	2000-2012 –> 2013-2019	5.8765604	22001.5392
Indo-Pacific	(27.85,27.95]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	1542	3.1760747	7958.0189	31.994368	2000-2012 –> 2013-2019	5.8765604	22001.5392
Indo-Pacific	(27.95,28.05]	7	cstar_tref ~ sal + temp + aou + silicate + silicate_star + phosphate + phosphate_star	2467	2.1513187	10798.8864	8.032399	1982-1999 –> 2000-2012	4.0117616	25270.5770
Indo-Pacific	(27.95,28.05]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	2467	2.1513187	10798.8864	8.032399	1982-1999 –> 2000-2012	4.0117616	25270.5770
Indo-Pacific	(27.95,28.05]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	2467	2.1513187	10798.8864	8.032399	1982-1999 –> 2000-2012	4.0117616	25270.5770
Indo-Pacific	(27.95,28.05]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	2467	2.1513187	10798.8864	8.032399	1982-1999 –> 2000-2012	4.0117616	25270.5770
Indo-Pacific	(27.95,28.05]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	2467	2.1513187	10798.8864	8.032399	1982-1999 –> 2000-2012	4.0117616	25270.5770
Indo-Pacific	(27.95,28.05]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	2467	2.1513187	10798.8864	8.032399	1982-1999 –> 2000-2012	4.0117616	25270.5770
Indo-Pacific	(27.95,28.05]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	2467	2.1513187	10798.8864	8.032399	1982-1999 –> 2000-2012	4.0117616	25270.5770
Indo-Pacific	(27.95,28.05]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	2467	2.1513187	10798.8864	8.032399	1982-1999 –> 2000-2012	4.0117616	25270.5770
Indo-Pacific	(27.95,28.05]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	2467	2.1513187	10798.8864	8.032399	1982-1999 –> 2000-2012	4.0117616	25270.5770
Indo-Pacific	(27.95,28.05]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	2467	2.1513187	10798.8864	8.032399	1982-1999 –> 2000-2012	4.0117616	25270.5770
Indo-Pacific	(27.95,28.05]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	1354	2.3752818	6203.2199	8.355865	2000-2012 –> 2013-2019	4.5266005	17002.1063
Indo-Pacific	(27.95,28.05]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	1354	2.3752818	6203.2199	8.355865	2000-2012 –> 2013-2019	4.5266005	17002.1063
Indo-Pacific	(27.95,28.05]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	1354	2.3752818	6203.2199	8.355865	2000-2012 –> 2013-2019	4.5266005	17002.1063
Indo-Pacific	(27.95,28.05]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	1354	2.3752818	6203.2199	8.355865	2000-2012 –> 2013-2019	4.5266005	17002.1063
Indo-Pacific	(27.95,28.05]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	1354	2.3752818	6203.2199	8.355865	2000-2012 –> 2013-2019	4.5266005	17002.1063
Indo-Pacific	(27.95,28.05]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	1354	2.3752818	6203.2199	8.355865	2000-2012 –> 2013-2019	4.5266005	17002.1063
Indo-Pacific	(27.95,28.05]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	1354	2.3752818	6203.2199	8.355865	2000-2012 –> 2013-2019	4.5266005	17002.1063
Indo-Pacific	(27.95,28.05]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	1354	2.3752818	6203.2199	8.355865	2000-2012 –> 2013-2019	4.5266005	17002.1063
Indo-Pacific	(27.95,28.05]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	1354	2.3752818	6203.2199	8.355865	2000-2012 –> 2013-2019	4.5266005	17002.1063
Indo-Pacific	(27.95,28.05]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	1354	2.3752818	6203.2199	8.355865	2000-2012 –> 2013-2019	4.5266005	17002.1063
Indo-Pacific	(28.05,28.1]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate_star + phosphate + phosphate_star	1718	1.6596412	6634.1553	10.112265	1982-1999 –> 2000-2012	3.0763941	15679.5548
Indo-Pacific	(28.05,28.1]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	1718	1.6596412	6634.1553	10.112265	1982-1999 –> 2000-2012	3.0763941	15679.5548
Indo-Pacific	(28.05,28.1]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	1718	1.6596412	6634.1553	10.112265	1982-1999 –> 2000-2012	3.0763941	15679.5548
Indo-Pacific	(28.05,28.1]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	1718	1.6596412	6634.1553	10.112265	1982-1999 –> 2000-2012	3.0763941	15679.5548
Indo-Pacific	(28.05,28.1]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	1718	1.6596412	6634.1553	10.112265	1982-1999 –> 2000-2012	3.0763941	15679.5548
Indo-Pacific	(28.05,28.1]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	1718	1.6596412	6634.1553	10.112265	1982-1999 –> 2000-2012	3.0763941	15679.5548
Indo-Pacific	(28.05,28.1]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate_star	1718	1.6596412	6634.1553	10.112265	1982-1999 –> 2000-2012	3.0763941	15679.5548
Indo-Pacific	(28.05,28.1]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	1718	1.6596412	6634.1553	10.112265	1982-1999 –> 2000-2012	3.0763941	15679.5548
Indo-Pacific	(28.05,28.1]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	1718	1.6596412	6634.1553	10.112265	1982-1999 –> 2000-2012	3.0763941	15679.5548
Indo-Pacific	(28.05,28.1]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate_star	1718	1.6596412	6634.1553	10.112265	1982-1999 –> 2000-2012	3.0763941	15679.5548
Indo-Pacific	(28.05,28.1]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	966	1.8625995	3961.0413	10.132291	2000-2012 –> 2013-2019	3.5222407	10595.1966
Indo-Pacific	(28.05,28.1]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	966	1.8625995	3961.0413	10.132291	2000-2012 –> 2013-2019	3.5222407	10595.1966
Indo-Pacific	(28.05,28.1]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	966	1.8625995	3961.0413	10.132291	2000-2012 –> 2013-2019	3.5222407	10595.1966
Indo-Pacific	(28.05,28.1]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	966	1.8625995	3961.0413	10.132291	2000-2012 –> 2013-2019	3.5222407	10595.1966
Indo-Pacific	(28.05,28.1]	9	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	966	1.8625995	3961.0413	10.132291	2000-2012 –> 2013-2019	3.5222407	10595.1966
Indo-Pacific	(28.05,28.1]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate_star	966	1.8625995	3961.0413	10.132291	2000-2012 –> 2013-2019	3.5222407	10595.1966
Indo-Pacific	(28.05,28.1]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	966	1.8625995	3961.0413	10.132291	2000-2012 –> 2013-2019	3.5222407	10595.1966
Indo-Pacific	(28.05,28.1]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate + phosphate_star	966	1.8625995	3961.0413	10.132291	2000-2012 –> 2013-2019	3.5222407	10595.1966
Indo-Pacific	(28.05,28.1]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate_star	966	1.8625995	3961.0413	10.132291	2000-2012 –> 2013-2019	3.5222407	10595.1966
Indo-Pacific	(28.05,28.1]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	966	1.8625995	3961.0413	10.132291	2000-2012 –> 2013-2019	3.5222407	10595.1966
Indo-Pacific	(28.1, Inf]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate	14627	1.1913295	46649.1221	13.074387	1982-1999 –> 2000-2012	2.2118753	105202.4830
Indo-Pacific	(28.1, Inf]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate + phosphate_star	14627	1.1913295	46649.1221	13.074387	1982-1999 –> 2000-2012	2.2118753	105202.4830
Indo-Pacific	(28.1, Inf]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	14627	1.1913295	46649.1221	13.074387	1982-1999 –> 2000-2012	2.2118753	105202.4830
Indo-Pacific	(28.1, Inf]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate	14627	1.1913295	46649.1221	13.074387	1982-1999 –> 2000-2012	2.2118753	105202.4830
Indo-Pacific	(28.1, Inf]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate_star	14627	1.1913295	46649.1221	13.074387	1982-1999 –> 2000-2012	2.2118753	105202.4830
Indo-Pacific	(28.1, Inf]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate	14627	1.1913295	46649.1221	13.074387	1982-1999 –> 2000-2012	2.2118753	105202.4830
Indo-Pacific	(28.1, Inf]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate_star	14627	1.1913295	46649.1221	13.074387	1982-1999 –> 2000-2012	2.2118753	105202.4830
Indo-Pacific	(28.1, Inf]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	14627	1.1913295	46649.1221	13.074387	1982-1999 –> 2000-2012	2.2118753	105202.4830
Indo-Pacific	(28.1, Inf]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	14627	1.1913295	46649.1221	13.074387	1982-1999 –> 2000-2012	2.2118753	105202.4830
Indo-Pacific	(28.1, Inf]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	14627	1.1913295	46649.1221	13.074387	1982-1999 –> 2000-2012	2.2118753	105202.4830
Indo-Pacific	(28.1, Inf]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate + phosphate_star	7804	1.3463098	26806.1024	10.367173	2000-2012 –> 2013-2019	2.5376393	73455.2245
Indo-Pacific	(28.1, Inf]	8	cstar_tref ~ sal + temp + aou + nitrate + silicate + silicate_star + phosphate + phosphate_star	7804	1.3463098	26806.1024	10.367173	2000-2012 –> 2013-2019	2.5376393	73455.2245
Indo-Pacific	(28.1, Inf]	7	cstar_tref ~ sal + temp + aou + nitrate + silicate_star + phosphate + phosphate_star	7804	1.3463098	26806.1024	10.367173	2000-2012 –> 2013-2019	2.5376393	73455.2245
Indo-Pacific	(28.1, Inf]	7	cstar_tref ~ sal + temp + aou + silicate + silicate_star + phosphate + phosphate_star	7804	1.3463098	26806.1024	10.367173	2000-2012 –> 2013-2019	2.5376393	73455.2245
Indo-Pacific	(28.1, Inf]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star	7804	1.3463098	26806.1024	10.367173	2000-2012 –> 2013-2019	2.5376393	73455.2245
Indo-Pacific	(28.1, Inf]	8	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + silicate_star + phosphate_star	7804	1.3463098	26806.1024	10.367173	2000-2012 –> 2013-2019	2.5376393	73455.2245
Indo-Pacific	(28.1, Inf]	7	cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate_star + phosphate_star	7804	1.3463098	26806.1024	10.367173	2000-2012 –> 2013-2019	2.5376393	73455.2245
Indo-Pacific	(28.1, Inf]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate	7804	1.3463098	26806.1024	10.367173	2000-2012 –> 2013-2019	2.5376393	73455.2245
Indo-Pacific	(28.1, Inf]	8	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate + phosphate_star	7804	1.3463098	26806.1024	10.367173	2000-2012 –> 2013-2019	2.5376393	73455.2245
Indo-Pacific	(28.1, Inf]	7	cstar_tref ~ sal + temp + oxygen + aou + silicate + silicate_star + phosphate_star	7804	1.3463098	26806.1024	10.367173	2000-2012 –> 2013-2019	2.5376393	73455.2245

5.3.2 per fitting unit

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

basin	gamma_slab	eras	rmse_sum_mean	ais_sum_mean
Atlantic	(-Inf,26]	1982-1999 –> 2000-2012	2.3449212	1197.9773
Atlantic	(-Inf,26]	2000-2012 –> 2013-2019	2.0853791	702.2046
Atlantic	(26,26.5]	1982-1999 –> 2000-2012	6.4070517	13316.0460
Atlantic	(26,26.5]	2000-2012 –> 2013-2019	6.2292432	8307.1154
Atlantic	(26.5,26.75]	1982-1999 –> 2000-2012	6.0826758	16678.5275
Atlantic	(26.5,26.75]	2000-2012 –> 2013-2019	6.4008998	10606.6280
Atlantic	(26.75,27]	1982-1999 –> 2000-2012	3.4513206	21493.5208
Atlantic	(26.75,27]	2000-2012 –> 2013-2019	4.1513732	15631.6529
Atlantic	(27,27.25]	1982-1999 –> 2000-2012	3.1780603	17367.0116
Atlantic	(27,27.25]	2000-2012 –> 2013-2019	3.8092897	12941.1749
Atlantic	(27.25,27.5]	1982-1999 –> 2000-2012	3.7585043	19804.2779
Atlantic	(27.25,27.5]	2000-2012 –> 2013-2019	4.8445332	14217.0362
Atlantic	(27.5,27.75]	1982-1999 –> 2000-2012	4.4225209	28805.9963
Atlantic	(27.5,27.75]	2000-2012 –> 2013-2019	5.0565164	20119.5192
Atlantic	(27.75,27.85]	1982-1999 –> 2000-2012	2.9888702	8619.6617
Atlantic	(27.75,27.85]	2000-2012 –> 2013-2019	3.0949714	5737.4410
Atlantic	(27.85,27.95]	1982-1999 –> 2000-2012	5.5922061	11351.1595
Atlantic	(27.85,27.95]	2000-2012 –> 2013-2019	6.1156663	7868.1993
Atlantic	(27.95,28.05]	1982-1999 –> 2000-2012	6.7144051	15128.4199
Atlantic	(27.95,28.05]	2000-2012 –> 2013-2019	6.7269357	9928.3462
Atlantic	(28.05,28.1]	1982-1999 –> 2000-2012	1.8168843	4949.0267
Atlantic	(28.05,28.1]	2000-2012 –> 2013-2019	2.0035216	3542.3614
Atlantic	(28.1,28.15]	1982-1999 –> 2000-2012	1.2960746	4297.7757
Atlantic	(28.1,28.15]	2000-2012 –> 2013-2019	1.5430592	3312.1683
Atlantic	(28.15,28.2]	1982-1999 –> 2000-2012	0.7244768	3147.4620
Atlantic	(28.15,28.2]	2000-2012 –> 2013-2019	0.9220982	2510.9133
Atlantic	(28.2, Inf]	1982-1999 –> 2000-2012	0.5862118	2293.0711
Atlantic	(28.2, Inf]	2000-2012 –> 2013-2019	0.6490752	2468.6696
Indo-Pacific	(-Inf,26]	1982-1999 –> 2000-2012	12.0393795	80350.0591
Indo-Pacific	(-Inf,26]	2000-2012 –> 2013-2019	11.4472636	51432.1573
Indo-Pacific	(26,26.5]	1982-1999 –> 2000-2012	8.3981502	75388.0758
Indo-Pacific	(26,26.5]	2000-2012 –> 2013-2019	8.5768441	48676.4609
Indo-Pacific	(26.5,26.75]	1982-1999 –> 2000-2012	7.4920814	59346.7074
Indo-Pacific	(26.5,26.75]	2000-2012 –> 2013-2019	7.8976637	38896.6972
Indo-Pacific	(26.75,27]	1982-1999 –> 2000-2012	7.5960628	70615.3985
Indo-Pacific	(26.75,27]	2000-2012 –> 2013-2019	8.5142930	46953.3037
Indo-Pacific	(27,27.25]	1982-1999 –> 2000-2012	6.7021130	79849.7973
Indo-Pacific	(27,27.25]	2000-2012 –> 2013-2019	7.9888450	55227.7255
Indo-Pacific	(27.25,27.5]	1982-1999 –> 2000-2012	5.9497930	68655.5399
Indo-Pacific	(27.25,27.5]	2000-2012 –> 2013-2019	6.8222175	47094.8546
Indo-Pacific	(27.5,27.75]	1982-1999 –> 2000-2012	5.1744982	65089.9113
Indo-Pacific	(27.5,27.75]	2000-2012 –> 2013-2019	6.4385346	45851.1416
Indo-Pacific	(27.75,27.85]	1982-1999 –> 2000-2012	4.5090898	23929.3942
Indo-Pacific	(27.75,27.85]	2000-2012 –> 2013-2019	5.5961533	17427.0947
Indo-Pacific	(27.85,27.95]	1982-1999 –> 2000-2012	4.7932875	31364.2073
Indo-Pacific	(27.85,27.95]	2000-2012 –> 2013-2019	5.8765604	22001.5392
Indo-Pacific	(27.95,28.05]	1982-1999 –> 2000-2012	4.0117616	25270.5770
Indo-Pacific	(27.95,28.05]	2000-2012 –> 2013-2019	4.5266005	17002.1063
Indo-Pacific	(28.05,28.1]	1982-1999 –> 2000-2012	3.0763941	15679.5548
Indo-Pacific	(28.05,28.1]	2000-2012 –> 2013-2019	3.5222407	10595.1966
Indo-Pacific	(28.1, Inf]	1982-1999 –> 2000-2012	2.2118753	105202.4830
Indo-Pacific	(28.1, Inf]	2000-2012 –> 2013-2019	2.5376393	73455.2245

5.4 Target variable coefficients

A data frame to map the target variable is prepared.

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

eras_backward <- lm_all_fitted_wide %>%
  arrange(era) %>% 
  group_by(basin, gamma_slab, model) %>% 
  mutate(eras = paste(lag(era), era, sep = " --> ")) %>% 
  ungroup() %>% 
  select(era, eras) %>% 
  unique()

eras_era <- full_join(eras_backward, eras_forward) %>% 
  filter(str_detect(eras, "NA") == FALSE)

# extend best model selection from eras to era
lm_best_target <- full_join(
  lm_best %>% select(basin, gamma_slab, model, eras),
  eras_era)

lm_best_target <- left_join(lm_best_target, lm_all_fitted_wide)

rm(eras_era, eras_forward, eras_backward,
   lm_all_fitted)

5.5 Plot selected model residuals

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

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

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

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

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

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

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

5.6 Cant coeffcients

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

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

# subtract coefficients of adjacent era
lm_best_long <- lm_best_long %>%
  arrange(era) %>%
  group_by(basin, gamma_slab, eras, model, term) %>%
  mutate(delta_coeff = estimate - lag(estimate)) %>%
  ungroup() %>%
  arrange(basin, gamma_slab, model, term, eras) %>%
  drop_na() %>%
  select(-c(era,estimate))

# pivot back to wide format
lm_best_cant <- lm_best_long %>%
  pivot_wider(values_from = delta_coeff,
              names_from = term,
              names_prefix = "delta_coeff_",
              values_fill = 0)

5.7 Write files

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

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

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

6 Model diagnotics

6.1 Selection criterion vs predictors

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

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
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
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
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
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	oxygen	phosphate	phosphate_star	sal	silicate	silicate_star	temp
Atlantic - 1982-1999 --> 2000-2012
(-Inf,26]	10	8	10	6	4	10	8	4	10
(26,26.5]	10	6	6	8	6	10	7	7	10
(26.5,26.75]	10	9	6	10	4	10	7	4	10
(26.75,27]	10	7	6	10	4	10	9	4	10
(27,27.25]	10	7	10	7	3	10	9	4	10
(27.25,27.5]	10	10	7	7	6	10	6	4	10
(27.5,27.75]	10	8	7	10	3	10	7	5	10
(27.75,27.85]	10	9	6	10	4	10	7	4	10
(27.85,27.95]	10	7	6	10	4	10	8	5	10
(27.95,28.05]	10	9	6	7	7	10	6	5	10
(28.05,28.1]	10	7	10	6	4	10	8	5	10
(28.1,28.15]	10	6	6	8	6	10	6	8	10
(28.15,28.2]	10	7	10	8	2	10	7	6	10
(28.2, Inf]	10	9	6	10	4	10	7	4	10
total	140.00	109.00	102.00	117.00	61.00	140.00	102.00	69.00	140.00
Atlantic - 2000-2012 --> 2013-2019
(-Inf,26]	10	7	10	8	2	10	8	5	10
(26,26.5]	10	6	9	9	2	10	8	6	10
(26.5,26.75]	10	7	10	6	4	10	9	4	10
(26.75,27]	10	6	9	6	5	10	9	5	10
(27,27.25]	10	6	9	7	4	10	10	4	10
(27.25,27.5]	10	10	7	7	6	10	6	4	10
(27.5,27.75]	10	6	9	9	2	10	8	6	10
(27.75,27.85]	10	7	10	8	2	10	8	5	10
(27.85,27.95]	10	6	6	6	8	10	8	6	10
(27.95,28.05]	10	9	6	7	7	10	6	5	10
(28.05,28.1]	10	9	10	6	4	10	7	4	10
(28.1,28.15]	10	7	10	6	4	10	7	6	10
(28.15,28.2]	10	9	10	7	3	10	7	4	10
(28.2, Inf]	10	9	10	6	4	10	7	4	10
total	140.00	104.00	125.00	98.00	57.00	140.00	108.00	68.00	140.00
Indo-Pacific - 1982-1999 --> 2000-2012
(-Inf,26]	10	6	6	6	8	10	8	6	10
(26,26.5]	10	6	6	9	5	10	9	5	10
(26.5,26.75]	10	7	10	6	4	10	9	4	10
(26.75,27]	10	7	10	6	4	10	9	4	10
(27,27.25]	10	9	6	10	4	10	7	4	10
(27.25,27.5]	10	6	6	8	6	10	6	8	10
(27.5,27.75]	10	7	10	6	4	10	8	5	10
(27.75,27.85]	10	8	6	10	4	10	8	4	10
(27.85,27.95]	10	9	6	10	4	10	7	4	10
(27.95,28.05]	10	6	9	9	2	10	8	6	10
(28.05,28.1]	10	10	9	7	4	10	6	4	10
(28.1, Inf]	10	7	10	6	4	10	8	5	10
total	120.00	88.00	94.00	93.00	53.00	120.00	93.00	59.00	120.00
Indo-Pacific - 2000-2012 --> 2013-2019
(-Inf,26]	10	9	6	10	4	10	7	4	10
(26,26.5]	10	7	10	6	4	10	8	5	10
(26.5,26.75]	10	6	9	9	2	10	8	6	10
(26.75,27]	10	7	10	6	4	10	9	4	10
(27,27.25]	10	8	10	6	4	10	8	4	10
(27.25,27.5]	10	6	6	8	6	10	6	8	10
(27.5,27.75]	10	7	10	6	4	10	7	6	10
(27.75,27.85]	10	8	6	8	6	10	9	3	10
(27.85,27.95]	10	9	6	10	4	10	7	4	10
(27.95,28.05]	10	7	10	8	2	10	8	5	10
(28.05,28.1]	10	9	10	6	4	10	7	4	10
(28.1, Inf]	10	6	6	6	8	10	8	6	10
total	120.00	89.00	99.00	89.00	52.00	120.00	92.00	59.00	120.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
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
cc2a956	Donghe-Zhu	2021-03-10
60689fb	Donghe-Zhu	2021-03-10
dba33c8	Donghe-Zhu	2021-03-09
17f1c4a	Donghe-Zhu	2021-03-09
c024d1a	Donghe-Zhu	2021-03-09
02f7242	Donghe-Zhu	2021-03-09
6f50bc6	Donghe-Zhu	2021-03-09
1691156	Donghe-Zhu	2021-03-08
c0ceaf8	Donghe-Zhu	2021-03-08
058e0a1	Donghe-Zhu	2021-03-08
112dea0	Donghe-Zhu	2021-03-08
1843412	Donghe-Zhu	2021-03-08
3fbbfa4	Donghe-Zhu	2021-03-07
627c8fb	Donghe-Zhu	2021-03-07
9ef3222	Donghe-Zhu	2021-03-05
8c1e978	Donghe-Zhu	2021-03-05
865f68c	Donghe-Zhu	2021-03-05
ee69bc1	Donghe-Zhu	2021-03-05
a79291f	Donghe-Zhu	2021-03-05
e8c6f30	Donghe-Zhu	2021-03-04
59288fe	Donghe-Zhu	2021-03-04
731abc8	Donghe-Zhu	2021-03-04
e2a5a33	Donghe-Zhu	2021-03-04
c7892c1	Donghe-Zhu	2021-03-04
924430b	Donghe-Zhu	2021-03-03
0d0bca1	Donghe-Zhu	2021-03-03
cb63c16	Donghe-Zhu	2021-03-03
ffda45a	Donghe-Zhu	2021-03-03
691ba81	Donghe-Zhu	2021-03-03
c5e45a2	Donghe-Zhu	2021-03-03
89c3e58	Donghe-Zhu	2021-03-03
c407a50	Donghe-Zhu	2021-03-03
c911669	Donghe-Zhu	2021-03-03
b71c719	Donghe-Zhu	2021-03-01
13666ca	Donghe-Zhu	2021-03-01
c6e60fe	Donghe-Zhu	2021-03-01
7a388f7	Donghe-Zhu	2021-03-01
799e913	Donghe-Zhu	2021-03-01
66ff99f	Donghe-Zhu	2021-03-01
ac9bb7a	Donghe-Zhu	2021-02-28
efdc047	Donghe-Zhu	2021-02-28
e9a7418	Donghe-Zhu	2021-02-28
e152917	Donghe-Zhu	2021-02-28
287123c	Donghe-Zhu	2021-02-27
54d5b5b	Donghe-Zhu	2021-02-27
330f064	Donghe-Zhu	2021-02-27
adbc9bc	Donghe-Zhu	2021-02-27
5937141	Donghe-Zhu	2021-02-27
4414bbf	Donghe-Zhu	2021-02-27
a265efb	Donghe-Zhu	2021-02-27
19edd1e	Donghe-Zhu	2021-02-27
f20483f	Donghe-Zhu	2021-02-26
6a2c7b3	Donghe-Zhu	2021-02-25
354c224	Donghe-Zhu	2021-02-24
1a0a88a	Donghe-Zhu	2021-02-24
57f701e	Donghe-Zhu	2021-02-24
06f3149	Donghe-Zhu	2021-02-16
401eab3	Donghe-Zhu	2021-02-15
e3bba84	Donghe-Zhu	2021-02-15
5dce4b1	Donghe-Zhu	2021-02-15
4469a0c	Donghe-Zhu	2021-02-13
5ae6a69	Donghe-Zhu	2021-02-10
05385dc	Donghe-Zhu	2021-02-10
f791ae4	Donghe-Zhu	2021-02-09
f71ae34	Donghe-Zhu	2021-02-09
c011832	Donghe-Zhu	2021-02-09
a145fa7	Donghe-Zhu	2021-02-09
c344e42	Donghe-Zhu	2021-02-08
2f095d7	Donghe-Zhu	2021-02-07
1fad5f1	Donghe-Zhu	2021-02-07
ca03c39	Donghe-Zhu	2021-02-07
e2ffc14	Donghe-Zhu	2021-02-05
cd7c52c	Donghe-Zhu	2021-02-04
bcf84f4	Donghe-Zhu	2021-02-02
a518739	Donghe-Zhu	2021-02-01
61666de	Donghe-Zhu	2021-01-31
865b582	Donghe-Zhu	2021-01-31
3e68089	Donghe-Zhu	2021-01-31
ecf335c	Donghe-Zhu	2021-01-31
a618965	Donghe-Zhu	2021-01-31
59e006e	Donghe-Zhu	2021-01-31
a1c8f87	Donghe-Zhu	2021-01-31
ae5c18f	Donghe-Zhu	2021-01-31
b50fe52	Donghe-Zhu	2021-01-31
ac99ae5	jens-daniel-mueller	2021-01-29
b5bdcaf	Donghe-Zhu	2021-01-29
442010d	Donghe-Zhu	2021-01-29
372adf5	Donghe-Zhu	2021-01-29
af8788e	Donghe-Zhu	2021-01-29
21c91c9	Donghe-Zhu	2021-01-29
eded038	Donghe-Zhu	2021-01-29
541d4dd	Donghe-Zhu	2021-01-29
6a75576	Donghe-Zhu	2021-01-28
16fba40	Donghe-Zhu	2021-01-28
12bc567	Donghe-Zhu	2021-01-27
ceed31b	Donghe-Zhu	2021-01-27
342402d	Donghe-Zhu	2021-01-27
5bad5c2	Donghe-Zhu	2021-01-27
61efb56	Donghe-Zhu	2021-01-25
48f638e	Donghe-Zhu	2021-01-25
c1cec47	Donghe-Zhu	2021-01-25
05ffb0c	Donghe-Zhu	2021-01-25
8b97165	Donghe-Zhu	2021-01-25
c569946	Donghe-Zhu	2021-01-24
a2f0d56	Donghe-Zhu	2021-01-23
28509fc	Donghe-Zhu	2021-01-23
4c28e4a	Donghe-Zhu	2021-01-22
24cc264	jens-daniel-mueller	2021-01-22
7891955	Donghe-Zhu	2021-01-21
d4cf1cb	Donghe-Zhu	2021-01-21
1f3e5b6	jens-daniel-mueller	2021-01-20
0e7bdf1	jens-daniel-mueller	2021-01-15
4571843	jens-daniel-mueller	2021-01-14
b3564aa	jens-daniel-mueller	2021-01-14
8d032c3	jens-daniel-mueller	2021-01-14
17dee1d	jens-daniel-mueller	2021-01-13
7cdea0c	jens-daniel-mueller	2021-01-06
fa85b93	jens-daniel-mueller	2021-01-06
e5cb81a	Donghe-Zhu	2021-01-05
a499f10	Donghe-Zhu	2021-01-05
fb8a752	Donghe-Zhu	2020-12-23
8fae0b2	Donghe-Zhu	2020-12-21
c8b76b3	jens-daniel-mueller	2020-12-19

Session information

sessionInfo()

R version 4.0.3 (2020-10-10)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: openSUSE Leap 15.2

Matrix products: default
BLAS:   /usr/local/R-4.0.3/lib64/R/lib/libRblas.so
LAPACK: /usr/local/R-4.0.3/lib64/R/lib/libRlapack.so

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
 [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8    
 [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8   
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C            
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
 [1] gt_0.2.2         corrr_0.4.3      broom_0.7.5      kableExtra_1.3.1
 [5] knitr_1.30       olsrr_0.5.3      GGally_2.0.0     lubridate_1.7.9 
 [9] metR_0.9.0       scico_1.2.0      patchwork_1.1.1  collapse_1.5.0  
[13] forcats_0.5.0    stringr_1.4.0    dplyr_1.0.2      purrr_0.3.4     
[17] readr_1.4.0      tidyr_1.1.2      tibble_3.0.4     ggplot2_3.3.3   
[21] tidyverse_1.3.0  workflowr_1.6.2 

loaded via a namespace (and not attached):
 [1] fs_1.5.0                 webshot_0.5.2            RColorBrewer_1.1-2      
 [4] httr_1.4.2               rprojroot_2.0.2          tools_4.0.3             
 [7] backports_1.1.10         R6_2.5.0                 nortest_1.0-4           
[10] DBI_1.1.0                colorspace_2.0-0         withr_2.3.0             
[13] gridExtra_2.3            tidyselect_1.1.0         curl_4.3                
[16] compiler_4.0.3           git2r_0.27.1             cli_2.2.0               
[19] rvest_0.3.6              xml2_1.3.2               sass_0.2.0              
[22] labeling_0.4.2           scales_1.1.1             checkmate_2.0.0         
[25] goftest_1.2-2            digest_0.6.27            foreign_0.8-80          
[28] rmarkdown_2.5            rio_0.5.16               pkgconfig_2.0.3         
[31] htmltools_0.5.0          highr_0.8                dbplyr_1.4.4            
[34] rlang_0.4.10             readxl_1.3.1             rstudioapi_0.13         
[37] farver_2.0.3             generics_0.1.0           jsonlite_1.7.2          
[40] zip_2.1.1                car_3.0-10               magrittr_2.0.1          
[43] Matrix_1.2-18            Rcpp_1.0.5               munsell_0.5.0           
[46] fansi_0.4.1              abind_1.4-5              lifecycle_0.2.0         
[49] stringi_1.5.3            whisker_0.4              yaml_2.2.1              
[52] carData_3.0-4            plyr_1.8.6               grid_4.0.3              
[55] blob_1.2.1               parallel_4.0.3           promises_1.1.1          
[58] crayon_1.3.4             lattice_0.20-41          haven_2.3.1             
[61] hms_0.5.3                pillar_1.4.7             reprex_0.3.0            
[64] glue_1.4.2               evaluate_0.14            RcppArmadillo_0.10.1.2.2
[67] data.table_1.13.6        modelr_0.1.8             vctrs_0.3.6             
[70] httpuv_1.5.4             cellranger_1.1.0         gtable_0.3.0            
[73] reshape_0.8.8            assertthat_0.2.1         xfun_0.20               
[76] openxlsx_4.2.3           RcppEigen_0.3.3.9.1      later_1.1.0.1           
[79] viridisLite_0.3.0        ellipsis_0.3.1           here_1.0.1