eMLR - model fitting

workflowr

• Summary
• Checks
• Past versions

Last updated: 2021-03-01

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: 61c083d

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 61c083d. 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:
    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	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, nitrate, silicate, phosphate, phosphate_star

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

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

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

i_gamma_slab <- unique(GLODAP_basin_era$gamma_slab)[1]
print(i_gamma_slab)

# subset one gamma slab
GLODAP_basin_era_slab <- GLODAP_basin_era %>%
  filter(gamma_slab == i_gamma_slab)

# fit the full linear model, i.e. all predictor combinations
lm_full <- lm(paste(
  params_local$MLR_target,
  paste(params_local$MLR_predictors, collapse = " + "),
  sep = " ~ "
),
data = GLODAP_basin_era_slab)

# fit linear models for all possible predictor combinations
# unfortunately, this functions does not provide model coefficients (yet)
lm_all <- ols_step_all_possible(lm_full)

# convert to tibble
lm_all <- as_tibble(lm_all)

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

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

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

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

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


#remove temp sal predictor combination
lm_all <- lm_all %>%
  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")
  )) %>%
  filter(!(
    str_detect(lm_coeff_filter, "temp") &
      str_detect(lm_coeff_filter, "phosphate_star")
  )) %>%
  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
)

3 Apply predictor threshold

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

• Minimum: 3
  
• 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:

• 39

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:

• aic

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, 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	model	nr_obs	rmse	aic	resid_max	eras	rmse_sum	aic_sum
Atlantic	(-Inf,26]	cstar_tref ~ sal + nitrate + phosphate_star	158	1.8666788	655.6194	7.225492	1982-1999 –> 2000-2012	3.9604888	1563.3157
Atlantic	(-Inf,26]	cstar_tref ~ sal + nitrate + silicate	158	1.8148164	646.7156	8.871634	1982-1999 –> 2000-2012	3.6615821	1502.1834
Atlantic	(-Inf,26]	cstar_tref ~ sal + nitrate + silicate + phosphate_star	158	1.6581469	620.1860	7.643496	1982-1999 –> 2000-2012	3.3615952	1444.0487
Atlantic	(-Inf,26]	cstar_tref ~ sal + oxygen + nitrate + phosphate_star	158	1.7528801	637.7428	7.185961	1982-1999 –> 2000-2012	3.5945001	1494.0498
Atlantic	(-Inf,26]	cstar_tref ~ sal + oxygen + nitrate + silicate	158	1.6268394	614.1625	8.329023	1982-1999 –> 2000-2012	3.2947666	1429.2589
Atlantic	(-Inf,26]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	158	1.6216413	615.1512	8.073641	1982-1999 –> 2000-2012	3.2856343	1431.2652
Atlantic	(-Inf,26]	cstar_tref ~ sal + oxygen + silicate + phosphate	158	1.9768012	675.7322	6.515418	1982-1999 –> 2000-2012	4.0711490	1585.5354
Atlantic	(-Inf,26]	cstar_tref ~ sal + oxygen + silicate + phosphate_star	158	1.9768012	675.7322	6.515418	1982-1999 –> 2000-2012	4.0711490	1585.5354
Atlantic	(-Inf,26]	cstar_tref ~ sal + silicate + phosphate	158	1.9790736	674.0953	6.354908	1982-1999 –> 2000-2012	4.0762175	1582.4535
Atlantic	(-Inf,26]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	158	1.9768012	675.7322	6.515418	1982-1999 –> 2000-2012	4.0711490	1585.5354
Atlantic	(-Inf,26]	cstar_tref ~ sal + nitrate + phosphate_star	72	2.0441204	317.2825	5.386407	2000-2012 –> 2013-2019	3.9107992	972.9019
Atlantic	(-Inf,26]	cstar_tref ~ sal + nitrate + silicate	72	1.8325789	301.5514	5.631280	2000-2012 –> 2013-2019	3.6473953	948.2671
Atlantic	(-Inf,26]	cstar_tref ~ sal + nitrate + silicate + phosphate_star	72	1.7638716	298.0488	5.052184	2000-2012 –> 2013-2019	3.4220185	918.2347
Atlantic	(-Inf,26]	cstar_tref ~ sal + oxygen + nitrate + phosphate_star	72	1.8628955	305.9142	4.679514	2000-2012 –> 2013-2019	3.6157756	943.6570
Atlantic	(-Inf,26]	cstar_tref ~ sal + oxygen + nitrate + silicate	72	1.7125308	293.7952	5.378554	2000-2012 –> 2013-2019	3.3393703	907.9577
Atlantic	(-Inf,26]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	72	1.7102518	295.6034	5.514006	2000-2012 –> 2013-2019	3.3318930	910.7546
Atlantic	(-Inf,26]	cstar_tref ~ sal + oxygen + silicate + phosphate	72	1.9854370	315.0880	4.406552	2000-2012 –> 2013-2019	3.9622382	990.8202
Atlantic	(-Inf,26]	cstar_tref ~ sal + oxygen + silicate + phosphate_star	72	1.9854370	315.0880	4.406552	2000-2012 –> 2013-2019	3.9622382	990.8202
Atlantic	(-Inf,26]	cstar_tref ~ sal + silicate + phosphate	72	1.9924896	313.5986	4.272959	2000-2012 –> 2013-2019	3.9715632	987.6939
Atlantic	(-Inf,26]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	72	1.9854370	315.0880	4.406552	2000-2012 –> 2013-2019	3.9622382	990.8202
Atlantic	(26,26.5]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	1046	4.5531745	6151.5246	12.264804	1982-1999 –> 2000-2012	9.0169247	15054.0544
Atlantic	(26,26.5]	cstar_tref ~ temp + oxygen + nitrate	1046	4.4337456	6093.9194	12.147787	1982-1999 –> 2000-2012	8.8915998	14990.4179
Atlantic	(26,26.5]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	1046	4.1683330	5966.7833	10.579872	1982-1999 –> 2000-2012	8.3007589	14634.0829
Atlantic	(26,26.5]	cstar_tref ~ temp + oxygen + nitrate + silicate	1046	3.9974267	5879.2009	14.350644	1982-1999 –> 2000-2012	8.0360096	14476.4398
Atlantic	(26,26.5]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	1046	3.9825944	5873.4242	15.972464	1982-1999 –> 2000-2012	8.0211559	14472.6469
Atlantic	(26,26.5]	cstar_tref ~ temp + oxygen + silicate	1046	5.1744624	6417.1140	15.577046	1982-1999 –> 2000-2012	10.3872223	15790.7603
Atlantic	(26,26.5]	cstar_tref ~ temp + oxygen + silicate + phosphate	1046	4.4898702	6122.2348	12.384079	1982-1999 –> 2000-2012	8.8545410	14956.3028
Atlantic	(26,26.5]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	1046	4.4898702	6122.2348	12.384079	1982-1999 –> 2000-2012	8.8545410	14956.3028
Atlantic	(26,26.5]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	1046	4.4898702	6122.2348	12.384079	1982-1999 –> 2000-2012	8.8545410	14956.3028
Atlantic	(26,26.5]	cstar_tref ~ temp + silicate + phosphate_star	1046	4.6317409	6185.3148	13.840219	1982-1999 –> 2000-2012	9.3965954	15284.9417
Atlantic	(26,26.5]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	567	4.4501173	3314.0594	11.610465	2000-2012 –> 2013-2019	9.0032917	9465.5840
Atlantic	(26,26.5]	cstar_tref ~ temp + oxygen + nitrate	567	4.2852908	3269.2599	12.306864	2000-2012 –> 2013-2019	8.7190363	9363.1793
Atlantic	(26,26.5]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	567	4.0757597	3214.4111	11.944274	2000-2012 –> 2013-2019	8.2440926	9181.1944
Atlantic	(26,26.5]	cstar_tref ~ temp + oxygen + nitrate + silicate	567	3.8537522	3150.8959	13.395954	2000-2012 –> 2013-2019	7.8511788	9030.0968
Atlantic	(26,26.5]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	567	3.8005103	3137.1198	14.416820	2000-2012 –> 2013-2019	7.7831048	9010.5440
Atlantic	(26,26.5]	cstar_tref ~ temp + oxygen + silicate + phosphate	567	4.4550925	3315.3265	11.906159	2000-2012 –> 2013-2019	8.9449627	9437.5613
Atlantic	(26,26.5]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	567	4.4550925	3315.3265	11.906159	2000-2012 –> 2013-2019	8.9449627	9437.5613
Atlantic	(26,26.5]	cstar_tref ~ temp + silicate + phosphate	567	4.8767953	3415.8860	13.895158	2000-2012 –> 2013-2019	10.0528007	9833.6237
Atlantic	(26,26.5]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	567	4.4550925	3315.3265	11.906159	2000-2012 –> 2013-2019	8.9449627	9437.5613
Atlantic	(26,26.5]	cstar_tref ~ temp + silicate + phosphate_star	567	4.4562951	3313.6326	11.744360	2000-2012 –> 2013-2019	9.0880360	9498.9474
Atlantic	(26.5,26.75]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	1397	3.5682017	7532.6548	10.073991	1982-1999 –> 2000-2012	6.9872721	17594.3324
Atlantic	(26.5,26.75]	cstar_tref ~ sal + oxygen + silicate + phosphate	1397	3.5707269	7532.6314	10.111585	1982-1999 –> 2000-2012	6.9903542	17592.9269
Atlantic	(26.5,26.75]	cstar_tref ~ sal + oxygen + silicate + phosphate_star	1397	3.5707269	7532.6314	10.111585	1982-1999 –> 2000-2012	6.9903542	17592.9269
Atlantic	(26.5,26.75]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	1397	3.5707269	7532.6314	10.111585	1982-1999 –> 2000-2012	6.9903542	17592.9269
Atlantic	(26.5,26.75]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	1397	3.7165763	7644.4856	12.547019	1982-1999 –> 2000-2012	7.2857185	17867.1404
Atlantic	(26.5,26.75]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	1397	3.7318580	7655.9503	26.389299	1982-1999 –> 2000-2012	7.2672629	17842.5719
Atlantic	(26.5,26.75]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	1397	3.3402414	7348.1986	9.940278	1982-1999 –> 2000-2012	6.4993586	17109.8616
Atlantic	(26.5,26.75]	cstar_tref ~ temp + oxygen + silicate + phosphate	1397	3.4413502	7429.5180	10.920943	1982-1999 –> 2000-2012	6.6456964	17243.1144
Atlantic	(26.5,26.75]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	1397	3.4413502	7429.5180	10.920943	1982-1999 –> 2000-2012	6.6456964	17243.1144
Atlantic	(26.5,26.75]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	1397	3.4413502	7429.5180	10.920943	1982-1999 –> 2000-2012	6.6456964	17243.1144
Atlantic	(26.5,26.75]	cstar_tref ~ sal + nitrate + silicate + phosphate_star	656	3.5914656	3551.1186	11.330901	2000-2012 –> 2013-2019	7.3317814	11213.3939
Atlantic	(26.5,26.75]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	656	3.5746856	3546.9743	10.593623	2000-2012 –> 2013-2019	7.1428873	11079.6291
Atlantic	(26.5,26.75]	cstar_tref ~ sal + oxygen + silicate + phosphate	656	3.5841303	3548.4361	11.020116	2000-2012 –> 2013-2019	7.1548572	11081.0675
Atlantic	(26.5,26.75]	cstar_tref ~ sal + oxygen + silicate + phosphate_star	656	3.5841303	3548.4361	11.020116	2000-2012 –> 2013-2019	7.1548572	11081.0675
Atlantic	(26.5,26.75]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	656	3.5841303	3548.4361	11.020116	2000-2012 –> 2013-2019	7.1548572	11081.0675
Atlantic	(26.5,26.75]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	656	3.6731778	3580.6344	12.617816	2000-2012 –> 2013-2019	7.3897541	11225.1199
Atlantic	(26.5,26.75]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	656	3.4254447	3491.0228	9.732422	2000-2012 –> 2013-2019	6.7656861	10839.2215
Atlantic	(26.5,26.75]	cstar_tref ~ temp + oxygen + silicate + phosphate	656	3.5288346	3528.0369	11.412513	2000-2012 –> 2013-2019	6.9701847	10957.5549
Atlantic	(26.5,26.75]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	656	3.5288346	3528.0369	11.412513	2000-2012 –> 2013-2019	6.9701847	10957.5549
Atlantic	(26.5,26.75]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	656	3.5288346	3528.0369	11.412513	2000-2012 –> 2013-2019	6.9701847	10957.5549
Atlantic	(26.75,27]	cstar_tref ~ sal + nitrate + silicate + phosphate_star	2398	2.6714987	11529.9688	38.212412	1982-1999 –> 2000-2012	5.1613576	26181.6825
Atlantic	(26.75,27]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	2398	2.5422316	11294.0999	39.132195	1982-1999 –> 2000-2012	4.8836368	25561.7519
Atlantic	(26.75,27]	cstar_tref ~ sal + oxygen + silicate + phosphate	2398	2.5470452	11301.1723	39.382843	1982-1999 –> 2000-2012	4.9052728	25611.7832
Atlantic	(26.75,27]	cstar_tref ~ sal + oxygen + silicate + phosphate_star	2398	2.5470452	11301.1723	39.382843	1982-1999 –> 2000-2012	4.9052728	25611.7832
Atlantic	(26.75,27]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	2398	2.5470452	11301.1723	39.382843	1982-1999 –> 2000-2012	4.9052728	25611.7832
Atlantic	(26.75,27]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	2398	2.7414356	11653.9075	22.253783	1982-1999 –> 2000-2012	5.1712500	26152.3190
Atlantic	(26.75,27]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	2398	2.2396447	10686.3266	20.763751	1982-1999 –> 2000-2012	4.0590677	23369.9626
Atlantic	(26.75,27]	cstar_tref ~ temp + oxygen + silicate + phosphate	2398	2.3001008	10812.0716	20.719299	1982-1999 –> 2000-2012	4.1580260	23625.2168
Atlantic	(26.75,27]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	2398	2.3001008	10812.0716	20.719299	1982-1999 –> 2000-2012	4.1580260	23625.2168
Atlantic	(26.75,27]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	2398	2.3001008	10812.0716	20.719299	1982-1999 –> 2000-2012	4.1580260	23625.2168
Atlantic	(26.75,27]	cstar_tref ~ sal + nitrate + silicate + phosphate_star	1263	2.8933893	6279.9334	27.695168	2000-2012 –> 2013-2019	5.5648880	17809.9022
Atlantic	(26.75,27]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	1263	2.7305419	6135.6060	22.303619	2000-2012 –> 2013-2019	5.2727735	17429.7059
Atlantic	(26.75,27]	cstar_tref ~ sal + oxygen + silicate + phosphate	1263	2.7359170	6138.5735	23.435321	2000-2012 –> 2013-2019	5.2829621	17439.7458
Atlantic	(26.75,27]	cstar_tref ~ sal + oxygen + silicate + phosphate_star	1263	2.7359170	6138.5735	23.435321	2000-2012 –> 2013-2019	5.2829621	17439.7458
Atlantic	(26.75,27]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	1263	2.7359170	6138.5735	23.435321	2000-2012 –> 2013-2019	5.2829621	17439.7458
Atlantic	(26.75,27]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	1263	3.1788514	6517.6081	32.113007	2000-2012 –> 2013-2019	5.9202870	18171.5156
Atlantic	(26.75,27]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	1263	2.6917313	6099.4450	22.955194	2000-2012 –> 2013-2019	4.9313760	16785.7717
Atlantic	(26.75,27]	cstar_tref ~ temp + oxygen + silicate + phosphate	1263	2.8252828	6219.7637	28.536665	2000-2012 –> 2013-2019	5.1253837	17031.8354
Atlantic	(26.75,27]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	1263	2.8252828	6219.7637	28.536665	2000-2012 –> 2013-2019	5.1253837	17031.8354
Atlantic	(26.75,27]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	1263	2.8252828	6219.7637	28.536665	2000-2012 –> 2013-2019	5.1253837	17031.8354
Atlantic	(27,27.25]	cstar_tref ~ sal + nitrate + silicate + phosphate_star	2031	2.2040749	8985.9589	10.013233	1982-1999 –> 2000-2012	4.5689583	20902.3812
Atlantic	(27,27.25]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	2031	2.1902432	8962.3874	9.942899	1982-1999 –> 2000-2012	4.5538173	20877.9179
Atlantic	(27,27.25]	cstar_tref ~ sal + oxygen + silicate + phosphate	2031	2.1910052	8961.8003	9.938424	1982-1999 –> 2000-2012	4.5832775	20938.3539
Atlantic	(27,27.25]	cstar_tref ~ sal + oxygen + silicate + phosphate_star	2031	2.1910052	8961.8003	9.938424	1982-1999 –> 2000-2012	4.5832775	20938.3539
Atlantic	(27,27.25]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	2031	2.1910052	8961.8003	9.938424	1982-1999 –> 2000-2012	4.5832775	20938.3539
Atlantic	(27,27.25]	cstar_tref ~ temp + oxygen + nitrate + silicate	2031	2.2533923	9075.8464	17.372391	1982-1999 –> 2000-2012	4.6876998	21143.3602
Atlantic	(27,27.25]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	2031	2.1067869	8804.5844	12.394761	1982-1999 –> 2000-2012	3.8168590	19030.0700
Atlantic	(27,27.25]	cstar_tref ~ temp + oxygen + silicate + phosphate	2031	2.6271564	9699.2224	12.195559	1982-1999 –> 2000-2012	4.7641127	21086.4149
Atlantic	(27,27.25]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	2031	2.6271564	9699.2224	12.195559	1982-1999 –> 2000-2012	4.7641127	21086.4149
Atlantic	(27,27.25]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	2031	2.6271564	9699.2224	12.195559	1982-1999 –> 2000-2012	4.7641127	21086.4149
Atlantic	(27,27.25]	cstar_tref ~ sal + nitrate + silicate + phosphate_star	1142	2.4064771	5258.5819	15.831426	2000-2012 –> 2013-2019	4.6105520	14244.5408
Atlantic	(27,27.25]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	1142	2.3332715	5190.0234	13.840333	2000-2012 –> 2013-2019	4.5235147	14152.4108
Atlantic	(27,27.25]	cstar_tref ~ sal + oxygen + silicate + phosphate	1142	2.3613801	5215.3740	15.167463	2000-2012 –> 2013-2019	4.5523853	14177.1744
Atlantic	(27,27.25]	cstar_tref ~ sal + oxygen + silicate + phosphate_star	1142	2.3613801	5215.3740	15.167463	2000-2012 –> 2013-2019	4.5523853	14177.1744
Atlantic	(27,27.25]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	1142	2.3613801	5215.3740	15.167463	2000-2012 –> 2013-2019	4.5523853	14177.1744
Atlantic	(27,27.25]	cstar_tref ~ sal + silicate + phosphate_star	1142	2.6802915	5502.7096	13.700856	2000-2012 –> 2013-2019	5.4265966	15380.0995
Atlantic	(27,27.25]	cstar_tref ~ temp + oxygen + nitrate + silicate	1142	2.6080258	5442.2833	24.208933	2000-2012 –> 2013-2019	4.8614181	14518.1297
Atlantic	(27,27.25]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	1142	2.5908455	5429.1876	20.743363	2000-2012 –> 2013-2019	4.6976323	14233.7721
Atlantic	(27,27.25]	cstar_tref ~ temp + oxygen + silicate + phosphate	1142	3.2173880	5921.8691	15.165532	2000-2012 –> 2013-2019	5.8445444	15621.0915
Atlantic	(27,27.25]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	1142	3.2173880	5921.8691	15.165532	2000-2012 –> 2013-2019	5.8445444	15621.0915
Atlantic	(27.25,27.5]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	2025	3.3018330	10598.3360	23.180595	1982-1999 –> 2000-2012	6.8189347	26018.7063
Atlantic	(27.25,27.5]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	2025	3.6658965	11019.9464	15.041394	1982-1999 –> 2000-2012	6.5633495	25322.7740
Atlantic	(27.25,27.5]	cstar_tref ~ temp + oxygen + nitrate	2025	3.3704618	10677.6526	10.689148	1982-1999 –> 2000-2012	6.0573060	24544.1053
Atlantic	(27.25,27.5]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	2025	3.3680346	10676.7350	10.384609	1982-1999 –> 2000-2012	5.9522077	24320.9237
Atlantic	(27.25,27.5]	cstar_tref ~ temp + oxygen + nitrate + silicate	2025	2.4916281	9456.0932	10.104738	1982-1999 –> 2000-2012	4.3239984	21121.3969
Atlantic	(27.25,27.5]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	2025	2.3326197	9191.0186	10.601623	1982-1999 –> 2000-2012	4.1616629	20847.8613
Atlantic	(27.25,27.5]	cstar_tref ~ temp + oxygen + silicate + phosphate	2025	3.6495556	11001.8530	15.721562	1982-1999 –> 2000-2012	6.4399972	25088.0695
Atlantic	(27.25,27.5]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	2025	3.6495556	11001.8530	15.721562	1982-1999 –> 2000-2012	6.4399972	25088.0695
Atlantic	(27.25,27.5]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	2025	3.6495556	11001.8530	15.721562	1982-1999 –> 2000-2012	6.4399972	25088.0695
Atlantic	(27.25,27.5]	cstar_tref ~ temp + silicate + phosphate_star	2025	3.6659043	11017.9550	15.059039	1982-1999 –> 2000-2012	6.7262296	25633.5733
Atlantic	(27.25,27.5]	cstar_tref ~ sal + nitrate + silicate + phosphate_star	1109	3.9462196	6203.9830	27.762167	2000-2012 –> 2013-2019	7.2956123	16858.2392
Atlantic	(27.25,27.5]	cstar_tref ~ sal + oxygen + nitrate + silicate	1109	3.6777016	6047.6804	19.783278	2000-2012 –> 2013-2019	7.1462844	16843.5532
Atlantic	(27.25,27.5]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	1109	3.6764640	6048.9339	20.304506	2000-2012 –> 2013-2019	6.9782970	16647.2699
Atlantic	(27.25,27.5]	cstar_tref ~ sal + oxygen + silicate + phosphate	1109	3.9750215	6220.1126	28.586368	2000-2012 –> 2013-2019	7.3205020	16869.6354
Atlantic	(27.25,27.5]	cstar_tref ~ sal + oxygen + silicate + phosphate_star	1109	3.9750215	6220.1126	28.586368	2000-2012 –> 2013-2019	7.3205020	16869.6354
Atlantic	(27.25,27.5]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	1109	3.9750215	6220.1126	28.586368	2000-2012 –> 2013-2019	7.3205020	16869.6354
Atlantic	(27.25,27.5]	cstar_tref ~ temp + oxygen + nitrate	1109	3.9761576	6218.7464	11.613376	2000-2012 –> 2013-2019	7.3466193	16896.3990
Atlantic	(27.25,27.5]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	1109	3.9721435	6218.5061	11.440074	2000-2012 –> 2013-2019	7.3401781	16895.2411
Atlantic	(27.25,27.5]	cstar_tref ~ temp + oxygen + nitrate + silicate	1109	3.2001964	5739.2103	12.589120	2000-2012 –> 2013-2019	5.6918245	15195.3035
Atlantic	(27.25,27.5]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	1109	2.8746758	5503.2804	14.300893	2000-2012 –> 2013-2019	5.2072955	14694.2990
Atlantic	(27.5,27.75]	cstar_tref ~ sal + oxygen + nitrate + phosphate_star	2804	2.9025448	13945.2240	19.212953	1982-1999 –> 2000-2012	6.1166042	33273.0065
Atlantic	(27.5,27.75]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	2804	2.8759174	13895.5398	18.881131	1982-1999 –> 2000-2012	6.0899038	33225.1527
Atlantic	(27.5,27.75]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	2804	3.7099257	15321.5617	17.533039	1982-1999 –> 2000-2012	6.6060295	33871.4143
Atlantic	(27.5,27.75]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	2804	3.9891085	15728.4554	15.599927	1982-1999 –> 2000-2012	6.8901927	34291.1396
Atlantic	(27.5,27.75]	cstar_tref ~ temp + oxygen + nitrate + silicate	2804	3.2454298	14571.4129	16.872661	1982-1999 –> 2000-2012	5.6444387	31714.9618
Atlantic	(27.5,27.75]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	2804	3.1572781	14418.9826	16.810973	1982-1999 –> 2000-2012	5.5460908	31532.7237
Atlantic	(27.5,27.75]	cstar_tref ~ temp + oxygen + silicate + phosphate	2804	3.7322940	15355.2726	18.019343	1982-1999 –> 2000-2012	6.6109791	33860.0731
Atlantic	(27.5,27.75]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	2804	3.7322940	15355.2726	18.019343	1982-1999 –> 2000-2012	6.6109791	33860.0731
Atlantic	(27.5,27.75]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	2804	3.7322940	15355.2726	18.019343	1982-1999 –> 2000-2012	6.6109791	33860.0731
Atlantic	(27.5,27.75]	cstar_tref ~ temp + silicate + phosphate_star	2804	3.7465465	15374.6470	18.562937	1982-1999 –> 2000-2012	6.6510914	33944.2347
Atlantic	(27.5,27.75]	cstar_tref ~ sal + nitrate + silicate + phosphate_star	1515	3.3077463	7936.0730	16.206508	2000-2012 –> 2013-2019	6.5222357	22453.7654
Atlantic	(27.5,27.75]	cstar_tref ~ sal + oxygen + nitrate + phosphate_star	1515	2.9543301	7593.6977	23.012182	2000-2012 –> 2013-2019	5.8568749	21538.9217
Atlantic	(27.5,27.75]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	1515	2.8464806	7483.0164	22.421750	2000-2012 –> 2013-2019	5.7223980	21378.5562
Atlantic	(27.5,27.75]	cstar_tref ~ sal + oxygen + phosphate	1515	3.7058316	8278.4041	14.544521	2000-2012 –> 2013-2019	7.1142988	23122.6925
Atlantic	(27.5,27.75]	cstar_tref ~ sal + oxygen + phosphate_star	1515	3.7058316	8278.4041	14.544521	2000-2012 –> 2013-2019	7.1142988	23122.6925
Atlantic	(27.5,27.75]	cstar_tref ~ sal + oxygen + silicate + phosphate	1515	3.2671643	7898.6687	16.135848	2000-2012 –> 2013-2019	6.4589075	22376.5370
Atlantic	(27.5,27.75]	cstar_tref ~ sal + oxygen + silicate + phosphate_star	1515	3.2671643	7898.6687	16.135848	2000-2012 –> 2013-2019	6.4589075	22376.5370
Atlantic	(27.5,27.75]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	1515	3.2671643	7898.6687	16.135848	2000-2012 –> 2013-2019	6.4589075	22376.5370
Atlantic	(27.5,27.75]	cstar_tref ~ sal + silicate + phosphate_star	1515	3.4628831	8072.9513	19.517533	2000-2012 –> 2013-2019	6.7421755	22700.5745
Atlantic	(27.5,27.75]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	1515	3.9981468	8512.4516	18.505654	2000-2012 –> 2013-2019	7.1554249	22931.4342
Atlantic	(27.75,27.85]	cstar_tref ~ sal + oxygen + nitrate + phosphate_star	1036	2.0673590	4456.8761	14.204780	1982-1999 –> 2000-2012	4.2394980	10262.7546
Atlantic	(27.75,27.85]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	1036	2.0673322	4458.8492	14.257932	1982-1999 –> 2000-2012	4.2060587	10225.8028
Atlantic	(27.75,27.85]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	1036	2.3495612	4722.0023	15.751424	1982-1999 –> 2000-2012	4.3651680	10330.4296
Atlantic	(27.75,27.85]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	1036	2.4544540	4812.4984	18.610143	1982-1999 –> 2000-2012	4.4779663	10431.2599
Atlantic	(27.75,27.85]	cstar_tref ~ temp + oxygen + nitrate + silicate	1036	2.3485017	4721.0678	16.287352	1982-1999 –> 2000-2012	4.4340791	10419.5862
Atlantic	(27.75,27.85]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	1036	2.3369261	4712.8298	15.853172	1982-1999 –> 2000-2012	4.3519055	10322.4352
Atlantic	(27.75,27.85]	cstar_tref ~ temp + oxygen + silicate + phosphate	1036	2.3475390	4720.2182	15.610251	1982-1999 –> 2000-2012	4.3638414	10329.5564
Atlantic	(27.75,27.85]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	1036	2.3475390	4720.2182	15.610251	1982-1999 –> 2000-2012	4.3638414	10329.5564
Atlantic	(27.75,27.85]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	1036	2.3475390	4720.2182	15.610251	1982-1999 –> 2000-2012	4.3638414	10329.5564
Atlantic	(27.75,27.85]	cstar_tref ~ temp + silicate + phosphate_star	1036	2.3804977	4747.1061	15.450631	1982-1999 –> 2000-2012	4.4466851	10418.9653
Atlantic	(27.75,27.85]	cstar_tref ~ sal + nitrate + silicate + phosphate_star	513	2.2842561	2315.3484	11.539783	2000-2012 –> 2013-2019	4.5828767	6991.9337
Atlantic	(27.75,27.85]	cstar_tref ~ sal + oxygen + nitrate + phosphate_star	513	2.1177551	2237.6968	14.216180	2000-2012 –> 2013-2019	4.1851141	6694.5729
Atlantic	(27.75,27.85]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	513	2.0908491	2226.5780	12.385110	2000-2012 –> 2013-2019	4.1581813	6685.4273
Atlantic	(27.75,27.85]	cstar_tref ~ sal + oxygen + silicate + phosphate	513	2.2850958	2315.7255	11.564796	2000-2012 –> 2013-2019	4.5875235	6995.7397
Atlantic	(27.75,27.85]	cstar_tref ~ sal + oxygen + silicate + phosphate_star	513	2.2850958	2315.7255	11.564796	2000-2012 –> 2013-2019	4.5875235	6995.7397
Atlantic	(27.75,27.85]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	513	2.2850958	2315.7255	11.564796	2000-2012 –> 2013-2019	4.5875235	6995.7397
Atlantic	(27.75,27.85]	cstar_tref ~ sal + silicate + phosphate_star	513	2.2851991	2313.7719	11.573863	2000-2012 –> 2013-2019	4.5932183	6996.8119
Atlantic	(27.75,27.85]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	513	2.7527803	2508.7702	15.319861	2000-2012 –> 2013-2019	5.0897064	7221.6000
Atlantic	(27.75,27.85]	cstar_tref ~ temp + oxygen + silicate + phosphate	513	2.7583236	2508.8342	15.146313	2000-2012 –> 2013-2019	5.1058626	7229.0524
Atlantic	(27.75,27.85]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	513	2.7583236	2508.8342	15.146313	2000-2012 –> 2013-2019	5.1058626	7229.0524
Atlantic	(27.85,27.95]	cstar_tref ~ temp + nitrate + phosphate_star	983	3.8389840	5444.3116	51.338778	1982-1999 –> 2000-2012	6.7922986	12139.2993
Atlantic	(27.85,27.95]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	983	3.8366506	5445.1162	50.582261	1982-1999 –> 2000-2012	6.7728918	12126.6121
Atlantic	(27.85,27.95]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	983	3.8383871	5446.0059	51.740798	1982-1999 –> 2000-2012	6.7907488	12142.1313
Atlantic	(27.85,27.95]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	983	3.8137613	5435.3521	51.408459	1982-1999 –> 2000-2012	6.7226058	12093.7996
Atlantic	(27.85,27.95]	cstar_tref ~ temp + oxygen + phosphate	983	3.8387704	5444.2022	52.175591	1982-1999 –> 2000-2012	6.7916245	12138.7733
Atlantic	(27.85,27.95]	cstar_tref ~ temp + oxygen + phosphate_star	983	3.8387704	5444.2022	52.175591	1982-1999 –> 2000-2012	6.7916245	12138.7733
Atlantic	(27.85,27.95]	cstar_tref ~ temp + oxygen + silicate + phosphate	983	3.8308048	5442.1184	50.328170	1982-1999 –> 2000-2012	6.7760129	12131.7618
Atlantic	(27.85,27.95]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	983	3.8308048	5442.1184	50.328170	1982-1999 –> 2000-2012	6.7760129	12131.7618
Atlantic	(27.85,27.95]	cstar_tref ~ temp + phosphate + phosphate_star	983	3.8387704	5444.2022	52.175591	1982-1999 –> 2000-2012	6.7916245	12138.7733
Atlantic	(27.85,27.95]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	983	3.8308048	5442.1184	50.328170	1982-1999 –> 2000-2012	6.7760129	12131.7618
Atlantic	(27.85,27.95]	cstar_tref ~ temp + nitrate + phosphate_star	569	3.8818373	3168.2312	24.594447	2000-2012 –> 2013-2019	7.7208213	8612.5428
Atlantic	(27.85,27.95]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	569	3.8561018	3162.6615	23.775108	2000-2012 –> 2013-2019	7.6927524	8607.7777
Atlantic	(27.85,27.95]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	569	3.8020643	3146.6013	23.264516	2000-2012 –> 2013-2019	7.6404514	8592.6072
Atlantic	(27.85,27.95]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	569	3.7705912	3139.1418	22.473597	2000-2012 –> 2013-2019	7.5843525	8574.4939
Atlantic	(27.85,27.95]	cstar_tref ~ temp + oxygen + phosphate	569	3.8865243	3169.6044	23.837557	2000-2012 –> 2013-2019	7.7252947	8613.8066
Atlantic	(27.85,27.95]	cstar_tref ~ temp + oxygen + phosphate_star	569	3.8865243	3169.6044	23.837557	2000-2012 –> 2013-2019	7.7252947	8613.8066
Atlantic	(27.85,27.95]	cstar_tref ~ temp + oxygen + silicate + phosphate	569	3.8608548	3164.0633	23.852541	2000-2012 –> 2013-2019	7.6916596	8606.1817
Atlantic	(27.85,27.95]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	569	3.8608548	3164.0633	23.852541	2000-2012 –> 2013-2019	7.6916596	8606.1817
Atlantic	(27.85,27.95]	cstar_tref ~ temp + phosphate + phosphate_star	569	3.8865243	3169.6044	23.837557	2000-2012 –> 2013-2019	7.7252947	8613.8066
Atlantic	(27.85,27.95]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	569	3.8608548	3164.0633	23.852541	2000-2012 –> 2013-2019	7.6916596	8606.1817
Atlantic	(27.95,28.05]	cstar_tref ~ oxygen + nitrate + silicate + phosphate_star	1204	4.7603831	7186.0742	34.069173	1982-1999 –> 2000-2012	8.7956135	16630.6086
Atlantic	(27.95,28.05]	cstar_tref ~ sal + oxygen + nitrate + phosphate_star	1204	4.4957529	7048.3486	26.854357	1982-1999 –> 2000-2012	8.6410505	16583.0894
Atlantic	(27.95,28.05]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	1204	4.1957717	6884.0621	30.223807	1982-1999 –> 2000-2012	7.7898469	15942.5131
Atlantic	(27.95,28.05]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	1204	4.6250887	7116.6453	24.067535	1982-1999 –> 2000-2012	9.5473302	17227.2242
Atlantic	(27.95,28.05]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	1204	4.1566877	6859.5262	23.454537	1982-1999 –> 2000-2012	8.2856315	16381.0167
Atlantic	(27.95,28.05]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	1204	3.9591041	6744.2548	25.125796	1982-1999 –> 2000-2012	7.5958976	15842.3119
Atlantic	(27.95,28.05]	cstar_tref ~ temp + oxygen + phosphate	1204	4.6374938	7121.0952	25.354382	1982-1999 –> 2000-2012	9.5744080	17239.6512
Atlantic	(27.95,28.05]	cstar_tref ~ temp + oxygen + silicate + phosphate	1204	4.5725557	7089.1381	24.432332	1982-1999 –> 2000-2012	9.4509550	17169.7271
Atlantic	(27.95,28.05]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	1204	4.5725557	7089.1381	24.432332	1982-1999 –> 2000-2012	9.4509550	17169.7271
Atlantic	(27.95,28.05]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	1204	4.5725557	7089.1381	24.432332	1982-1999 –> 2000-2012	9.4509550	17169.7271
Atlantic	(27.95,28.05]	cstar_tref ~ sal + oxygen + nitrate + phosphate_star	666	3.9013752	3715.3165	28.636553	2000-2012 –> 2013-2019	8.3971281	10763.6651
Atlantic	(27.95,28.05]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	666	3.4861342	3567.4190	32.154120	2000-2012 –> 2013-2019	7.6819059	10451.4811
Atlantic	(27.95,28.05]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	666	3.9030940	3715.9032	23.856002	2000-2012 –> 2013-2019	8.0597817	10575.4294
Atlantic	(27.95,28.05]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	666	3.5394291	3587.6281	27.891571	2000-2012 –> 2013-2019	7.4985332	10331.8828
Atlantic	(27.95,28.05]	cstar_tref ~ temp + oxygen + phosphate	666	4.4690013	3894.2499	24.065812	2000-2012 –> 2013-2019	9.1064951	11015.3451
Atlantic	(27.95,28.05]	cstar_tref ~ temp + oxygen + phosphate_star	666	4.4690013	3894.2499	24.065812	2000-2012 –> 2013-2019	9.1064951	11015.3451
Atlantic	(27.95,28.05]	cstar_tref ~ temp + oxygen + silicate + phosphate	666	4.4291749	3884.3263	23.917524	2000-2012 –> 2013-2019	9.0017306	10973.4643
Atlantic	(27.95,28.05]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	666	4.4291749	3884.3263	23.917524	2000-2012 –> 2013-2019	9.0017306	10973.4643
Atlantic	(27.95,28.05]	cstar_tref ~ temp + phosphate + phosphate_star	666	4.4690013	3894.2499	24.065812	2000-2012 –> 2013-2019	9.1064951	11015.3451
Atlantic	(27.95,28.05]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	666	4.4291749	3884.3263	23.917524	2000-2012 –> 2013-2019	9.0017306	10973.4643
Atlantic	(28.05,28.1]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	804	1.1861347	2568.1385	10.721163	1982-1999 –> 2000-2012	2.1695803	5597.7711
Atlantic	(28.05,28.1]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	804	1.1826495	2563.4068	11.075362	1982-1999 –> 2000-2012	2.1591614	5577.8130
Atlantic	(28.05,28.1]	cstar_tref ~ temp + oxygen + nitrate + silicate	804	1.1969610	2582.7489	10.502979	1982-1999 –> 2000-2012	2.2010569	5657.1009
Atlantic	(28.05,28.1]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	804	1.1823967	2565.0630	11.000864	1982-1999 –> 2000-2012	2.1548805	5572.5740
Atlantic	(28.05,28.1]	cstar_tref ~ temp + oxygen + phosphate	804	1.1853678	2565.0986	11.052937	1982-1999 –> 2000-2012	2.1623589	5578.5605
Atlantic	(28.05,28.1]	cstar_tref ~ temp + oxygen + phosphate_star	804	1.1853678	2565.0986	11.052937	1982-1999 –> 2000-2012	2.1623589	5578.5605
Atlantic	(28.05,28.1]	cstar_tref ~ temp + oxygen + silicate + phosphate	804	1.1824509	2563.1368	10.968185	1982-1999 –> 2000-2012	2.1574045	5574.1061
Atlantic	(28.05,28.1]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	804	1.1824509	2563.1368	10.968185	1982-1999 –> 2000-2012	2.1574045	5574.1061
Atlantic	(28.05,28.1]	cstar_tref ~ temp + phosphate + phosphate_star	804	1.1853678	2565.0986	11.052937	1982-1999 –> 2000-2012	2.1623589	5578.5605
Atlantic	(28.05,28.1]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	804	1.1824509	2563.1368	10.968185	1982-1999 –> 2000-2012	2.1574045	5574.1061
Atlantic	(28.05,28.1]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	434	1.2393117	1429.8734	11.911018	2000-2012 –> 2013-2019	2.4254464	3998.0119
Atlantic	(28.05,28.1]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	434	1.2345404	1426.5252	12.392191	2000-2012 –> 2013-2019	2.4171899	3989.9320
Atlantic	(28.05,28.1]	cstar_tref ~ temp + oxygen + nitrate + silicate	434	1.2485131	1436.2941	11.740228	2000-2012 –> 2013-2019	2.4454742	4019.0430
Atlantic	(28.05,28.1]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	434	1.2338376	1428.0309	12.265693	2000-2012 –> 2013-2019	2.4162343	3993.0939
Atlantic	(28.05,28.1]	cstar_tref ~ temp + oxygen + phosphate	434	1.2459665	1432.5219	12.315642	2000-2012 –> 2013-2019	2.4313344	3997.6205
Atlantic	(28.05,28.1]	cstar_tref ~ temp + oxygen + phosphate_star	434	1.2459665	1432.5219	12.315642	2000-2012 –> 2013-2019	2.4313344	3997.6205
Atlantic	(28.05,28.1]	cstar_tref ~ temp + oxygen + silicate + phosphate	434	1.2341739	1426.2674	12.180217	2000-2012 –> 2013-2019	2.4166248	3989.4043
Atlantic	(28.05,28.1]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	434	1.2341739	1426.2674	12.180217	2000-2012 –> 2013-2019	2.4166248	3989.4043
Atlantic	(28.05,28.1]	cstar_tref ~ temp + phosphate + phosphate_star	434	1.2459665	1432.5219	12.315642	2000-2012 –> 2013-2019	2.4313344	3997.6205
Atlantic	(28.05,28.1]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	434	1.2341739	1426.2674	12.180217	2000-2012 –> 2013-2019	2.4166248	3989.4043
Atlantic	(28.1,28.15]	cstar_tref ~ temp + nitrate + phosphate_star	936	0.9058253	2481.0956	10.374751	1982-1999 –> 2000-2012	1.6190427	5290.8067
Atlantic	(28.1,28.15]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	936	0.8931308	2456.6753	10.319681	1982-1999 –> 2000-2012	1.6061851	5267.7940
Atlantic	(28.1,28.15]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	936	0.8844217	2438.3315	10.086815	1982-1999 –> 2000-2012	1.5970909	5248.0509
Atlantic	(28.1,28.15]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	936	0.8807386	2432.5195	9.872456	1982-1999 –> 2000-2012	1.5856244	5215.7966
Atlantic	(28.1,28.15]	cstar_tref ~ temp + oxygen + phosphate	936	0.8846006	2436.7100	10.051092	1982-1999 –> 2000-2012	1.6121536	5297.9638
Atlantic	(28.1,28.15]	cstar_tref ~ temp + oxygen + phosphate_star	936	0.8846006	2436.7100	10.051092	1982-1999 –> 2000-2012	1.6121536	5297.9638
Atlantic	(28.1,28.15]	cstar_tref ~ temp + oxygen + silicate + phosphate	936	0.8830714	2435.4712	10.082608	1982-1999 –> 2000-2012	1.5884109	5218.4149
Atlantic	(28.1,28.15]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	936	0.8830714	2435.4712	10.082608	1982-1999 –> 2000-2012	1.5884109	5218.4149
Atlantic	(28.1,28.15]	cstar_tref ~ temp + phosphate + phosphate_star	936	0.8846006	2436.7100	10.051092	1982-1999 –> 2000-2012	1.6121536	5297.9638
Atlantic	(28.1,28.15]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	936	0.8830714	2435.4712	10.082608	1982-1999 –> 2000-2012	1.5884109	5218.4149
Atlantic	(28.1,28.15]	cstar_tref ~ temp + nitrate + phosphate_star	497	1.1083671	1522.6954	11.843631	2000-2012 –> 2013-2019	2.0141924	4003.7910
Atlantic	(28.1,28.15]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	497	1.0775245	1496.6431	11.777260	2000-2012 –> 2013-2019	1.9706553	3953.3185
Atlantic	(28.1,28.15]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	497	1.0559567	1476.5454	11.355356	2000-2012 –> 2013-2019	1.9403784	3914.8768
Atlantic	(28.1,28.15]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	497	1.0490759	1472.0472	11.035771	2000-2012 –> 2013-2019	1.9298146	3904.5667
Atlantic	(28.1,28.15]	cstar_tref ~ temp + oxygen + phosphate	497	1.0606641	1478.9667	11.544277	2000-2012 –> 2013-2019	1.9452647	3915.6768
Atlantic	(28.1,28.15]	cstar_tref ~ temp + oxygen + phosphate_star	497	1.0606641	1478.9667	11.544277	2000-2012 –> 2013-2019	1.9452647	3915.6768
Atlantic	(28.1,28.15]	cstar_tref ~ temp + oxygen + silicate + phosphate	497	1.0602381	1480.5675	11.532519	2000-2012 –> 2013-2019	1.9433095	3916.0387
Atlantic	(28.1,28.15]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	497	1.0602381	1480.5675	11.532519	2000-2012 –> 2013-2019	1.9433095	3916.0387
Atlantic	(28.1,28.15]	cstar_tref ~ temp + phosphate + phosphate_star	497	1.0606641	1478.9667	11.544277	2000-2012 –> 2013-2019	1.9452647	3915.6768
Atlantic	(28.1,28.15]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	497	1.0602381	1480.5675	11.532519	2000-2012 –> 2013-2019	1.9433095	3916.0387
Atlantic	(28.15,28.2]	cstar_tref ~ sal + oxygen + nitrate	1540	0.6939029	3254.8271	2.510912	1982-1999 –> 2000-2012	1.4529919	8350.2359
Atlantic	(28.15,28.2]	cstar_tref ~ sal + oxygen + nitrate + phosphate_star	1540	0.6650037	3125.8056	2.696128	1982-1999 –> 2000-2012	1.3251043	7601.7097
Atlantic	(28.15,28.2]	cstar_tref ~ sal + oxygen + nitrate + silicate	1540	0.6938712	3256.6862	2.523944	1982-1999 –> 2000-2012	1.4249928	8187.1205
Atlantic	(28.15,28.2]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	1540	0.6185754	2904.8953	2.211087	1982-1999 –> 2000-2012	1.2620075	7269.0379
Atlantic	(28.15,28.2]	cstar_tref ~ temp + nitrate + phosphate_star	1540	0.6352104	2982.6300	2.274003	1982-1999 –> 2000-2012	1.4034113	8131.1130
Atlantic	(28.15,28.2]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	1540	0.6283932	2951.3960	2.229212	1982-1999 –> 2000-2012	1.2745157	7332.0984
Atlantic	(28.15,28.2]	cstar_tref ~ temp + oxygen + nitrate	1540	0.6539635	3072.2429	2.506572	1982-1999 –> 2000-2012	1.2717050	7251.1464
Atlantic	(28.15,28.2]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	1540	0.6051750	2835.4386	2.198295	1982-1999 –> 2000-2012	1.1760176	6665.1429
Atlantic	(28.15,28.2]	cstar_tref ~ temp + oxygen + nitrate + silicate	1540	0.5772872	2690.1312	2.226077	1982-1999 –> 2000-2012	1.1364464	6427.8538
Atlantic	(28.15,28.2]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	1540	0.5669166	2636.2981	2.345127	1982-1999 –> 2000-2012	1.1259531	6375.0445
Atlantic	(28.15,28.2]	cstar_tref ~ sal + oxygen + nitrate	819	1.0269474	2377.7769	10.877467	2000-2012 –> 2013-2019	1.7208503	5632.6040
Atlantic	(28.15,28.2]	cstar_tref ~ sal + oxygen + nitrate + phosphate_star	819	0.8339184	2038.7282	5.301488	2000-2012 –> 2013-2019	1.4989220	5164.5337
Atlantic	(28.15,28.2]	cstar_tref ~ sal + oxygen + nitrate + silicate	819	0.9045133	2171.8345	5.875227	2000-2012 –> 2013-2019	1.5983844	5428.5206
Atlantic	(28.15,28.2]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	819	0.8297177	2032.4562	5.671364	2000-2012 –> 2013-2019	1.4482931	4937.3514
Atlantic	(28.15,28.2]	cstar_tref ~ temp + nitrate + phosphate_star	819	0.9292212	2213.9782	9.786344	2000-2012 –> 2013-2019	1.5644316	5196.6082
Atlantic	(28.15,28.2]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	819	0.8897687	2144.9132	7.849984	2000-2012 –> 2013-2019	1.5181619	5096.3091
Atlantic	(28.15,28.2]	cstar_tref ~ temp + oxygen + nitrate	819	0.9641793	2274.4703	9.030606	2000-2012 –> 2013-2019	1.6181427	5346.7132
Atlantic	(28.15,28.2]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	819	0.8724916	2112.7945	7.808630	2000-2012 –> 2013-2019	1.4776666	4948.2331
Atlantic	(28.15,28.2]	cstar_tref ~ temp + oxygen + nitrate + silicate	819	0.8696972	2107.5398	8.465903	2000-2012 –> 2013-2019	1.4469843	4797.6710
Atlantic	(28.15,28.2]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	819	0.8677435	2105.8561	8.181408	2000-2012 –> 2013-2019	1.4346601	4742.1542
Atlantic	(28.2, Inf]	cstar_tref ~ sal + nitrate + silicate + phosphate_star	2693	0.6073805	4968.9446	6.599511	1982-1999 –> 2000-2012	1.1735078	11289.6600
Atlantic	(28.2, Inf]	cstar_tref ~ sal + oxygen + nitrate	2693	0.4955127	3870.5572	8.085578	1982-1999 –> 2000-2012	1.0325571	9797.8500
Atlantic	(28.2, Inf]	cstar_tref ~ sal + oxygen + nitrate + phosphate_star	2693	0.4644953	3524.3978	6.704526	1982-1999 –> 2000-2012	0.9658160	8942.8000
Atlantic	(28.2, Inf]	cstar_tref ~ sal + oxygen + nitrate + silicate	2693	0.4778466	3677.0273	6.696866	1982-1999 –> 2000-2012	0.9678429	8925.8507
Atlantic	(28.2, Inf]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	2693	0.4618076	3495.1416	7.084808	1982-1999 –> 2000-2012	0.9515248	8741.7369
Atlantic	(28.2, Inf]	cstar_tref ~ sal + oxygen + silicate + phosphate_star	2693	0.7721739	6261.8769	7.611345	1982-1999 –> 2000-2012	1.4614046	14042.9241
Atlantic	(28.2, Inf]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	2693	0.7753565	6284.0299	5.358553	1982-1999 –> 2000-2012	1.4510348	13917.6842
Atlantic	(28.2, Inf]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	2693	0.7710158	6253.7932	4.326067	1982-1999 –> 2000-2012	1.3595153	12862.1642
Atlantic	(28.2, Inf]	cstar_tref ~ temp + oxygen + nitrate + silicate	2693	0.6190159	5071.1463	4.259725	1982-1999 –> 2000-2012	1.1183234	10459.6829
Atlantic	(28.2, Inf]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	2693	0.5451838	4389.0814	3.304918	1982-1999 –> 2000-2012	1.0135140	9304.2502
Atlantic	(28.2, Inf]	cstar_tref ~ sal + nitrate + phosphate_star	1557	0.8458975	3907.4245	14.110764	2000-2012 –> 2013-2019	1.5924911	9985.8538
Atlantic	(28.2, Inf]	cstar_tref ~ sal + nitrate + silicate + phosphate_star	1557	0.6309055	2996.2688	2.438558	2000-2012 –> 2013-2019	1.2382861	7965.2135
Atlantic	(28.2, Inf]	cstar_tref ~ sal + oxygen + nitrate	1557	0.5156492	2366.0835	7.293748	2000-2012 –> 2013-2019	1.0111620	6236.6407
Atlantic	(28.2, Inf]	cstar_tref ~ sal + oxygen + nitrate + phosphate_star	1557	0.4442799	1904.1849	3.889221	2000-2012 –> 2013-2019	0.9087752	5428.5827
Atlantic	(28.2, Inf]	cstar_tref ~ sal + oxygen + nitrate + silicate	1557	0.4718588	2091.7254	3.398074	2000-2012 –> 2013-2019	0.9497054	5768.7527
Atlantic	(28.2, Inf]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	1557	0.4385296	1865.6176	5.067855	2000-2012 –> 2013-2019	0.9003372	5360.7592
Atlantic	(28.2, Inf]	cstar_tref ~ sal + oxygen + silicate + phosphate_star	1557	0.8389904	3883.8932	2.971888	2000-2012 –> 2013-2019	1.6111643	10145.7701
Atlantic	(28.2, Inf]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	1557	0.8389904	3883.8932	2.971888	2000-2012 –> 2013-2019	1.6111643	10145.7701
Atlantic	(28.2, Inf]	cstar_tref ~ temp + oxygen + nitrate + silicate	1557	0.8180488	3805.1796	6.506855	2000-2012 –> 2013-2019	1.4370647	8876.3259
Atlantic	(28.2, Inf]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	1557	0.7225700	3420.7086	8.405924	2000-2012 –> 2013-2019	1.2677538	7809.7900
Indo-Pacific	(-Inf,26]	cstar_tref ~ temp + nitrate + phosphate_star	5288	7.8589632	36820.7539	32.936421	1982-1999 –> 2000-2012	16.2174312	87789.0400
Indo-Pacific	(-Inf,26]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	5288	7.7395948	36660.8846	30.024084	1982-1999 –> 2000-2012	15.9976756	87457.3454
Indo-Pacific	(-Inf,26]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	5288	7.5583116	36410.2175	34.882285	1982-1999 –> 2000-2012	15.4930477	86632.0709
Indo-Pacific	(-Inf,26]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	5288	7.1609299	35841.0290	31.554355	1982-1999 –> 2000-2012	14.6425732	85219.0198
Indo-Pacific	(-Inf,26]	cstar_tref ~ temp + oxygen + phosphate	5288	7.7813287	36715.7598	32.534467	1982-1999 –> 2000-2012	15.9375394	87331.6542
Indo-Pacific	(-Inf,26]	cstar_tref ~ temp + oxygen + phosphate_star	5288	7.7813287	36715.7598	32.534467	1982-1999 –> 2000-2012	15.9375394	87331.6542
Indo-Pacific	(-Inf,26]	cstar_tref ~ temp + oxygen + silicate + phosphate	5288	7.3907730	36173.1511	28.492599	1982-1999 –> 2000-2012	15.1050761	85989.6943
Indo-Pacific	(-Inf,26]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	5288	7.3907730	36173.1511	28.492599	1982-1999 –> 2000-2012	15.1050761	85989.6943
Indo-Pacific	(-Inf,26]	cstar_tref ~ temp + phosphate + phosphate_star	5288	7.7813287	36715.7598	32.534467	1982-1999 –> 2000-2012	15.9375394	87331.6542
Indo-Pacific	(-Inf,26]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	5288	7.3907730	36173.1511	28.492599	1982-1999 –> 2000-2012	15.1050761	85989.6943
Indo-Pacific	(-Inf,26]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	2816	7.3587257	19244.2962	26.521004	2000-2012 –> 2013-2019	15.0983206	55905.1807
Indo-Pacific	(-Inf,26]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	2816	7.4098204	19283.2663	28.866050	2000-2012 –> 2013-2019	14.9681320	55693.4838
Indo-Pacific	(-Inf,26]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	2816	6.9178353	18898.3293	26.006595	2000-2012 –> 2013-2019	14.0787652	54739.3584
Indo-Pacific	(-Inf,26]	cstar_tref ~ temp + oxygen + phosphate	2816	7.6384054	19452.3815	26.635654	2000-2012 –> 2013-2019	15.4197341	56168.1413
Indo-Pacific	(-Inf,26]	cstar_tref ~ temp + oxygen + phosphate_star	2816	7.6384054	19452.3815	26.635654	2000-2012 –> 2013-2019	15.4197341	56168.1413
Indo-Pacific	(-Inf,26]	cstar_tref ~ temp + oxygen + silicate + phosphate	2816	7.1344030	19069.9395	24.441528	2000-2012 –> 2013-2019	14.5251760	55243.0906
Indo-Pacific	(-Inf,26]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	2816	7.1344030	19069.9395	24.441528	2000-2012 –> 2013-2019	14.5251760	55243.0906
Indo-Pacific	(-Inf,26]	cstar_tref ~ temp + phosphate + phosphate_star	2816	7.6384054	19452.3815	26.635654	2000-2012 –> 2013-2019	15.4197341	56168.1413
Indo-Pacific	(-Inf,26]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	2816	7.1344030	19069.9395	24.441528	2000-2012 –> 2013-2019	14.5251760	55243.0906
Indo-Pacific	(-Inf,26]	cstar_tref ~ temp + silicate + phosphate_star	2816	7.3741321	19254.0751	27.379074	2000-2012 –> 2013-2019	15.1958698	56024.6147
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	5525	5.2631756	34042.3878	45.728160	1982-1999 –> 2000-2012	10.7299998	82140.8115
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	5525	5.2726312	34062.2219	53.955562	1982-1999 –> 2000-2012	10.1582823	80427.0605
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	5525	4.7637711	32942.7584	51.203206	1982-1999 –> 2000-2012	9.1789275	77747.7736
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + oxygen + phosphate	5525	5.4949921	34516.6714	53.690387	1982-1999 –> 2000-2012	10.7077534	81879.0953
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + oxygen + phosphate_star	5525	5.4949921	34516.6714	53.690387	1982-1999 –> 2000-2012	10.7077534	81879.0953
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + oxygen + silicate + phosphate	5525	4.9304894	33320.8635	50.805047	1982-1999 –> 2000-2012	9.5666218	78877.1440
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	5525	4.9304894	33320.8635	50.805047	1982-1999 –> 2000-2012	9.5666218	78877.1440
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + phosphate + phosphate_star	5525	5.4949921	34516.6714	53.690387	1982-1999 –> 2000-2012	10.7077534	81879.0953
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	5525	4.9304894	33320.8635	50.805047	1982-1999 –> 2000-2012	9.5666218	78877.1440
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + silicate + phosphate_star	5525	5.2809695	34077.6831	44.865042	1982-1999 –> 2000-2012	10.8686335	82511.3285
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + nitrate + phosphate_star	2918	5.2681921	17988.5320	39.285977	2000-2012 –> 2013-2019	10.7234846	52425.0809
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	2918	4.9958063	17680.7080	32.184396	2000-2012 –> 2013-2019	10.2589819	51723.0958
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	2918	5.1539157	17862.5457	47.713000	2000-2012 –> 2013-2019	10.4265469	51924.7676
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	2918	4.5779873	17172.9954	45.249093	2000-2012 –> 2013-2019	9.3417584	50115.7538
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + oxygen + phosphate_star	2918	5.4774950	18215.9063	35.116956	2000-2012 –> 2013-2019	10.9724872	52732.5777
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + oxygen + silicate + phosphate	2918	4.8517131	17509.9060	34.225344	2000-2012 –> 2013-2019	9.7822025	50830.7695
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	2918	4.8517131	17509.9060	34.225344	2000-2012 –> 2013-2019	9.7822025	50830.7695
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + phosphate + phosphate_star	2918	5.4774950	18215.9063	35.116956	2000-2012 –> 2013-2019	10.9724872	52732.5777
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	2918	4.8517131	17509.9060	34.225344	2000-2012 –> 2013-2019	9.7822025	50830.7695
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + silicate + phosphate_star	2918	5.0016139	17685.4884	31.434569	2000-2012 –> 2013-2019	10.2825834	51763.1714
Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + nitrate + phosphate_star	4544	4.1279653	25790.1399	36.509986	1982-1999 –> 2000-2012	8.2137297	61345.5422
Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	4544	4.0218809	25555.5345	32.963712	1982-1999 –> 2000-2012	8.0834114	61038.1226
Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	4544	4.1104968	25753.6003	34.578793	1982-1999 –> 2000-2012	8.1719835	61236.0528
Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	4544	4.0218771	25557.5260	32.942351	1982-1999 –> 2000-2012	8.0175184	60836.4236
Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + oxygen + phosphate	4544	4.4727413	26519.1509	35.203832	1982-1999 –> 2000-2012	8.7384525	62616.5797
Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + oxygen + phosphate_star	4544	4.4727413	26519.1509	35.203832	1982-1999 –> 2000-2012	8.7384525	62616.5797
Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + oxygen + silicate + phosphate	4544	4.2095240	25969.9463	32.334600	1982-1999 –> 2000-2012	8.3064308	61561.5980
Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	4544	4.2095240	25969.9463	32.334600	1982-1999 –> 2000-2012	8.3064308	61561.5980
Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	4544	4.2095240	25969.9463	32.334600	1982-1999 –> 2000-2012	8.3064308	61561.5980
Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + silicate + phosphate_star	4544	4.2542459	26063.9877	34.848227	1982-1999 –> 2000-2012	8.3702566	61712.1452
Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + nitrate + phosphate_star	2447	4.4875150	14301.6430	39.084814	2000-2012 –> 2013-2019	8.6154803	40091.7829
Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	2447	4.3388127	14138.7232	34.446250	2000-2012 –> 2013-2019	8.3606936	39694.2577
Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	2447	4.4113430	14219.8580	32.128406	2000-2012 –> 2013-2019	8.5218397	39973.4583
Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	2447	4.3205048	14120.0290	31.479242	2000-2012 –> 2013-2019	8.3423820	39677.5550
Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + oxygen + phosphate	2447	4.8558146	14687.6708	19.899362	2000-2012 –> 2013-2019	9.3285558	41206.8216
Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + oxygen + phosphate_star	2447	4.8558146	14687.6708	19.899362	2000-2012 –> 2013-2019	9.3285558	41206.8216
Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + oxygen + silicate + phosphate	2447	4.5556780	14377.4214	23.174250	2000-2012 –> 2013-2019	8.7652020	40347.3676
Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	2447	4.5556780	14377.4214	23.174250	2000-2012 –> 2013-2019	8.7652020	40347.3676
Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	2447	4.5556780	14377.4214	23.174250	2000-2012 –> 2013-2019	8.7652020	40347.3676
Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + silicate + phosphate_star	2447	4.7016885	14529.8140	29.527447	2000-2012 –> 2013-2019	8.9559344	40593.8017
Indo-Pacific	(26.75,27]	cstar_tref ~ temp + nitrate + phosphate_star	5354	5.7624413	33957.5698	32.774171	1982-1999 –> 2000-2012	10.3611936	78320.9959
Indo-Pacific	(26.75,27]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	5354	5.5953751	33644.5315	31.604100	1982-1999 –> 2000-2012	10.1571679	77888.4176
Indo-Pacific	(26.75,27]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	5354	4.7402141	31868.5232	21.890228	1982-1999 –> 2000-2012	9.0078158	75108.3206
Indo-Pacific	(26.75,27]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	5354	4.6674184	31704.8041	20.858437	1982-1999 –> 2000-2012	8.8798032	74750.4482
Indo-Pacific	(26.75,27]	cstar_tref ~ temp + oxygen + phosphate	5354	4.8452343	32101.1707	23.622430	1982-1999 –> 2000-2012	9.1905960	75610.9426
Indo-Pacific	(26.75,27]	cstar_tref ~ temp + oxygen + phosphate_star	5354	4.8452343	32101.1707	23.622430	1982-1999 –> 2000-2012	9.1905960	75610.9426
Indo-Pacific	(26.75,27]	cstar_tref ~ temp + oxygen + silicate + phosphate	5354	4.8218019	32051.2593	23.312826	1982-1999 –> 2000-2012	9.1486740	75498.8053
Indo-Pacific	(26.75,27]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	5354	4.8218019	32051.2593	23.312826	1982-1999 –> 2000-2012	9.1486740	75498.8053
Indo-Pacific	(26.75,27]	cstar_tref ~ temp + phosphate + phosphate_star	5354	4.8452343	32101.1707	23.622430	1982-1999 –> 2000-2012	9.1905960	75610.9426
Indo-Pacific	(26.75,27]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	5354	4.8218019	32051.2593	23.312826	1982-1999 –> 2000-2012	9.1486740	75498.8053
Indo-Pacific	(26.75,27]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	2868	5.7115592	18147.9658	21.260383	2000-2012 –> 2013-2019	11.0446157	51280.3434
Indo-Pacific	(26.75,27]	cstar_tref ~ sal + oxygen + silicate + phosphate	2868	5.8400264	18273.5532	23.549730	2000-2012 –> 2013-2019	11.2859513	51628.1892
Indo-Pacific	(26.75,27]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	2868	5.1160215	17514.3463	20.245716	2000-2012 –> 2013-2019	9.8562356	49382.8696
Indo-Pacific	(26.75,27]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	2868	5.0448252	17435.9616	18.619356	2000-2012 –> 2013-2019	9.7122436	49140.7657
Indo-Pacific	(26.75,27]	cstar_tref ~ temp + oxygen + phosphate	2868	5.2448998	17655.0527	17.217990	2000-2012 –> 2013-2019	10.0901341	49756.2234
Indo-Pacific	(26.75,27]	cstar_tref ~ temp + oxygen + phosphate_star	2868	5.2448998	17655.0527	17.217990	2000-2012 –> 2013-2019	10.0901341	49756.2234
Indo-Pacific	(26.75,27]	cstar_tref ~ temp + oxygen + silicate + phosphate	2868	5.2254515	17635.7438	17.046500	2000-2012 –> 2013-2019	10.0472534	49687.0031
Indo-Pacific	(26.75,27]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	2868	5.2254515	17635.7438	17.046500	2000-2012 –> 2013-2019	10.0472534	49687.0031
Indo-Pacific	(26.75,27]	cstar_tref ~ temp + phosphate + phosphate_star	2868	5.2448998	17655.0527	17.217990	2000-2012 –> 2013-2019	10.0901341	49756.2234
Indo-Pacific	(26.75,27]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	2868	5.2254515	17635.7438	17.046500	2000-2012 –> 2013-2019	10.0472534	49687.0031
Indo-Pacific	(27,27.25]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	6466	4.6876932	38342.8849	21.746845	1982-1999 –> 2000-2012	8.8474367	88628.6682
Indo-Pacific	(27,27.25]	cstar_tref ~ temp + nitrate + phosphate_star	6466	4.7823605	38597.4428	71.626560	1982-1999 –> 2000-2012	8.5941220	87335.1901
Indo-Pacific	(27,27.25]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	6466	4.6663213	38281.7913	72.764927	1982-1999 –> 2000-2012	8.1605463	85484.2572
Indo-Pacific	(27,27.25]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	6466	4.7286162	38453.2900	65.782175	1982-1999 –> 2000-2012	8.5031493	87019.5720
Indo-Pacific	(27,27.25]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	6466	4.1253304	36690.2467	48.506869	1982-1999 –> 2000-2012	7.4805153	83177.0650
Indo-Pacific	(27,27.25]	cstar_tref ~ temp + oxygen + phosphate	6466	5.1276743	39499.0354	52.755848	1982-1999 –> 2000-2012	9.3481294	90036.9072
Indo-Pacific	(27,27.25]	cstar_tref ~ temp + oxygen + phosphate_star	6466	5.1276743	39499.0354	52.755848	1982-1999 –> 2000-2012	9.3481294	90036.9072
Indo-Pacific	(27,27.25]	cstar_tref ~ temp + oxygen + silicate + phosphate	6466	4.6882058	38342.2992	36.052942	1982-1999 –> 2000-2012	8.5937887	87511.8009
Indo-Pacific	(27,27.25]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	6466	4.6882058	38342.2992	36.052942	1982-1999 –> 2000-2012	8.5937887	87511.8009
Indo-Pacific	(27,27.25]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	6466	4.6882058	38342.2992	36.052942	1982-1999 –> 2000-2012	8.5937887	87511.8009
Indo-Pacific	(27,27.25]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	3470	5.1816896	21278.6438	20.340514	2000-2012 –> 2013-2019	9.8693827	59621.5287
Indo-Pacific	(27,27.25]	cstar_tref ~ sal + oxygen + silicate + phosphate	3470	5.4833669	21669.3654	20.944689	2000-2012 –> 2013-2019	10.4723249	60815.7390
Indo-Pacific	(27,27.25]	cstar_tref ~ sal + oxygen + silicate + phosphate_star	3470	5.4833669	21669.3654	20.944689	2000-2012 –> 2013-2019	10.4723249	60815.7390
Indo-Pacific	(27,27.25]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	3470	5.4833669	21669.3654	20.944689	2000-2012 –> 2013-2019	10.4723249	60815.7390
Indo-Pacific	(27,27.25]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	3470	6.4677784	22815.2522	66.498629	2000-2012 –> 2013-2019	11.1340997	61097.0435
Indo-Pacific	(27,27.25]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	3470	5.9200231	22201.1155	51.106562	2000-2012 –> 2013-2019	10.6486393	60654.4055
Indo-Pacific	(27,27.25]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	3470	4.9160028	20913.3540	34.970533	2000-2012 –> 2013-2019	9.0413332	57603.6007
Indo-Pacific	(27,27.25]	cstar_tref ~ temp + oxygen + silicate + phosphate	3470	5.3226341	21462.8935	27.632792	2000-2012 –> 2013-2019	10.0108400	59805.1927
Indo-Pacific	(27,27.25]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	3470	5.3226341	21462.8935	27.632792	2000-2012 –> 2013-2019	10.0108400	59805.1927
Indo-Pacific	(27,27.25]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	3470	5.3226341	21462.8935	27.632792	2000-2012 –> 2013-2019	10.0108400	59805.1927
Indo-Pacific	(27.25,27.5]	cstar_tref ~ oxygen + nitrate + silicate + phosphate_star	5847	3.6972581	31896.0421	35.971452	1982-1999 –> 2000-2012	6.8980977	72667.8810
Indo-Pacific	(27.25,27.5]	cstar_tref ~ oxygen + silicate + phosphate	5847	3.9752652	32741.8567	34.799073	1982-1999 –> 2000-2012	7.3114607	74165.4363
Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	5847	3.4222130	30994.0484	31.732058	1982-1999 –> 2000-2012	6.3354476	70281.8372
Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + oxygen + silicate + phosphate	5847	3.7568967	32083.1670	30.801943	1982-1999 –> 2000-2012	6.8441839	72284.8816
Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + oxygen + silicate + phosphate_star	5847	3.7568967	32083.1670	30.801943	1982-1999 –> 2000-2012	6.8441839	72284.8816
Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	5847	3.7568967	32083.1670	30.801943	1982-1999 –> 2000-2012	6.8441839	72284.8816
Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	5847	3.6204935	31652.6884	33.074778	1982-1999 –> 2000-2012	6.6947648	71789.7182
Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + oxygen + silicate + phosphate	5847	3.9732569	32737.9473	34.309542	1982-1999 –> 2000-2012	7.2597842	73926.7734
Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	5847	3.9732569	32737.9473	34.309542	1982-1999 –> 2000-2012	7.2597842	73926.7734
Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	5847	3.9732569	32737.9473	34.309542	1982-1999 –> 2000-2012	7.2597842	73926.7734
Indo-Pacific	(27.25,27.5]	cstar_tref ~ oxygen + nitrate + silicate + phosphate_star	3127	4.3080384	18019.9002	34.741395	2000-2012 –> 2013-2019	8.0052965	49915.9424
Indo-Pacific	(27.25,27.5]	cstar_tref ~ oxygen + silicate + phosphate	3127	4.6569059	18504.8903	33.434540	2000-2012 –> 2013-2019	8.6321710	51246.7469
Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	3127	3.9789503	17524.9283	30.123196	2000-2012 –> 2013-2019	7.4011633	48518.9767
Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + oxygen + silicate + phosphate	3127	4.4121255	18169.2074	29.187384	2000-2012 –> 2013-2019	8.1690222	50252.3744
Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + oxygen + silicate + phosphate_star	3127	4.4121255	18169.2074	29.187384	2000-2012 –> 2013-2019	8.1690222	50252.3744
Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	3127	4.4121255	18169.2074	29.187384	2000-2012 –> 2013-2019	8.1690222	50252.3744
Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	3127	4.2264297	17902.2919	31.188098	2000-2012 –> 2013-2019	7.8469231	49554.9803
Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + oxygen + silicate + phosphate	3127	4.6568828	18506.8594	33.373407	2000-2012 –> 2013-2019	8.6301397	51244.8067
Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	3127	4.6568828	18506.8594	33.373407	2000-2012 –> 2013-2019	8.6301397	51244.8067
Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	3127	4.6568828	18506.8594	33.373407	2000-2012 –> 2013-2019	8.6301397	51244.8067
Indo-Pacific	(27.5,27.75]	cstar_tref ~ oxygen + nitrate + silicate + phosphate_star	5836	3.5543091	31375.8221	15.602716	1982-1999 –> 2000-2012	6.1053269	69569.3582
Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + nitrate + phosphate_star	5836	3.5608466	31395.2710	28.785023	1982-1999 –> 2000-2012	6.1945501	70103.8833
Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + nitrate + silicate + phosphate_star	5836	3.3957705	30843.2286	29.089521	1982-1999 –> 2000-2012	5.7808728	67946.6354
Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + oxygen + nitrate + phosphate_star	5836	3.5607118	31396.8290	28.693320	1982-1999 –> 2000-2012	6.1879632	70067.6811
Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	5836	3.2980912	30504.5604	26.307453	1982-1999 –> 2000-2012	5.6183359	67163.0675
Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + oxygen + silicate + phosphate	5836	3.4363673	30981.9413	25.357370	1982-1999 –> 2000-2012	5.8135896	68031.7039
Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + oxygen + silicate + phosphate_star	5836	3.4363673	30981.9413	25.357370	1982-1999 –> 2000-2012	5.8135896	68031.7039
Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	5836	3.4363673	30981.9413	25.357370	1982-1999 –> 2000-2012	5.8135896	68031.7039
Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + silicate + phosphate_star	5836	3.6584353	31710.8495	29.513815	1982-1999 –> 2000-2012	6.1698773	69648.9367
Indo-Pacific	(27.5,27.75]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	5836	3.5367123	31319.8925	11.083583	1982-1999 –> 2000-2012	6.0868856	69510.0613
Indo-Pacific	(27.5,27.75]	cstar_tref ~ oxygen + nitrate + silicate + phosphate_star	3104	4.2166617	17754.3540	17.626283	2000-2012 –> 2013-2019	7.7709707	49130.1761
Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + nitrate + phosphate_star	3104	4.3230640	17907.0615	32.442685	2000-2012 –> 2013-2019	7.8839106	49302.3325
Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + nitrate + silicate + phosphate_star	3104	4.1948886	17722.2154	32.318244	2000-2012 –> 2013-2019	7.5906591	48565.4440
Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + oxygen + nitrate + phosphate_star	3104	4.3053456	17883.5652	31.479458	2000-2012 –> 2013-2019	7.8660574	49280.3943
Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	3104	4.0075075	17440.5265	28.440018	2000-2012 –> 2013-2019	7.3055987	47945.0868
Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + oxygen + silicate + phosphate	3104	4.2088784	17742.8844	27.865993	2000-2012 –> 2013-2019	7.6452457	48724.8257
Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + oxygen + silicate + phosphate_star	3104	4.2088784	17742.8844	27.865993	2000-2012 –> 2013-2019	7.6452457	48724.8257
Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	3104	4.2088784	17742.8844	27.865993	2000-2012 –> 2013-2019	7.6452457	48724.8257
Indo-Pacific	(27.5,27.75]	cstar_tref ~ temp + oxygen + nitrate + silicate	3104	4.1769632	17695.6309	12.989269	2000-2012 –> 2013-2019	7.7660139	49184.9870
Indo-Pacific	(27.5,27.75]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	3104	4.1565100	17667.1576	11.836017	2000-2012 –> 2013-2019	7.6932222	48987.0501
Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	2317	2.7272370	11236.6024	21.622314	1982-1999 –> 2000-2012	4.7089141	24416.8498
Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	2317	2.7770530	11320.4838	22.304986	1982-1999 –> 2000-2012	4.7428380	24450.3103
Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	2317	2.7269514	11238.1171	21.350081	1982-1999 –> 2000-2012	4.6654113	24282.2891
Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + oxygen + phosphate	2317	2.7770643	11318.5027	22.348926	1982-1999 –> 2000-2012	4.7459420	24456.1732
Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + oxygen + phosphate_star	2317	2.7770643	11318.5027	22.348926	1982-1999 –> 2000-2012	4.7459420	24456.1732
Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + oxygen + silicate + phosphate	2317	2.7373185	11253.7009	20.385504	1982-1999 –> 2000-2012	4.6897567	24340.8662
Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	2317	2.7373185	11253.7009	20.385504	1982-1999 –> 2000-2012	4.6897567	24340.8662
Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + phosphate + phosphate_star	2317	2.7770643	11318.5027	22.348926	1982-1999 –> 2000-2012	4.7459420	24456.1732
Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	2317	2.7373185	11253.7009	20.385504	1982-1999 –> 2000-2012	4.6897567	24340.8662
Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + silicate + phosphate_star	2317	2.7507868	11274.4455	21.919660	1982-1999 –> 2000-2012	4.7347893	24460.0365
Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + nitrate + silicate + phosphate_star	1264	3.1488202	6498.7630	29.448400	2000-2012 –> 2013-2019	5.8762396	17735.6753
Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	1264	3.1390093	6492.8741	31.279189	2000-2012 –> 2013-2019	5.8653527	17729.9579
Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + oxygen + silicate + phosphate	1264	3.1768073	6521.1329	29.025171	2000-2012 –> 2013-2019	5.9117926	17770.8822
Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + oxygen + silicate + phosphate_star	1264	3.1768073	6521.1329	29.025171	2000-2012 –> 2013-2019	5.9117926	17770.8822
Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	1264	3.1768073	6521.1329	29.025171	2000-2012 –> 2013-2019	5.9117926	17770.8822
Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + silicate + phosphate_star	1264	3.1799375	6521.6226	29.845922	2000-2012 –> 2013-2019	5.9334554	17800.6665
Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	1264	3.2045095	6543.0819	27.185936	2000-2012 –> 2013-2019	5.9317465	17779.6843
Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	1264	3.1907808	6534.2282	30.135098	2000-2012 –> 2013-2019	5.9177322	17772.3453
Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + oxygen + silicate + phosphate	1264	3.2403449	6571.1951	27.030880	2000-2012 –> 2013-2019	5.9776634	17824.8960
Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	1264	3.2403449	6571.1951	27.030880	2000-2012 –> 2013-2019	5.9776634	17824.8960
Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + nitrate + phosphate_star	2907	2.7674113	14177.8510	31.750508	1982-1999 –> 2000-2012	4.8898936	31603.4801
Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	2907	2.7621671	14168.8231	31.138313	1982-1999 –> 2000-2012	4.8840444	31594.1657
Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + oxygen + nitrate	2907	2.7499312	14141.0108	28.547833	1982-1999 –> 2000-2012	4.9295241	31779.5857
Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	2907	2.7472864	14137.4165	28.371289	1982-1999 –> 2000-2012	4.8684199	31559.9473
Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + oxygen + nitrate + silicate	2907	2.7476621	14138.2115	28.325089	1982-1999 –> 2000-2012	4.9235034	31764.9698
Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	2907	2.7445406	14133.6026	28.109511	1982-1999 –> 2000-2012	4.8649180	31555.2745
Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + oxygen + phosphate	2907	2.7969089	14239.4940	33.225970	1982-1999 –> 2000-2012	4.9565475	31804.3068
Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + oxygen + silicate + phosphate	2907	2.7948101	14237.1294	33.015412	1982-1999 –> 2000-2012	4.9429857	31761.2604
Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	2907	2.7948101	14237.1294	33.015412	1982-1999 –> 2000-2012	4.9429857	31761.2604
Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	2907	2.7948101	14237.1294	33.015412	1982-1999 –> 2000-2012	4.9429857	31761.2604
Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + oxygen + nitrate + phosphate_star	1542	3.3109222	8080.2538	31.535597	2000-2012 –> 2013-2019	6.0927821	22290.3805
Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + oxygen + nitrate + silicate	1542	3.3366152	8104.0934	29.544640	2000-2012 –> 2013-2019	6.1199466	22317.2947
Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	1542	3.3033826	8075.2229	30.820497	2000-2012 –> 2013-2019	6.0846779	22286.1697
Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + nitrate + phosphate_star	1542	3.3221692	8088.7121	33.915337	2000-2012 –> 2013-2019	6.0895805	22266.5632
Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + nitrate + silicate	1542	3.3200594	8086.7530	34.851867	2000-2012 –> 2013-2019	6.0832309	22255.6898
Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	1542	3.2965108	8066.8008	32.540164	2000-2012 –> 2013-2019	6.0586779	22235.6239
Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + oxygen + nitrate	1542	3.3252514	8091.5720	31.716600	2000-2012 –> 2013-2019	6.0751825	22232.5829
Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	1542	3.3151344	8084.1748	31.858968	2000-2012 –> 2013-2019	6.0624208	22221.5913
Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + oxygen + nitrate + silicate	1542	3.3028816	8072.7552	31.234941	2000-2012 –> 2013-2019	6.0505437	22210.9666
Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	1542	3.2939987	8066.4498	31.381712	2000-2012 –> 2013-2019	6.0385393	22200.0524
Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	2467	2.2792684	11077.9408	8.417011	1982-1999 –> 2000-2012	4.2142805	25822.1041
Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	2467	2.2945444	11110.8990	7.469816	1982-1999 –> 2000-2012	4.3094440	26141.7320
Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + oxygen + nitrate + silicate	2467	2.2694946	11056.7378	8.186678	1982-1999 –> 2000-2012	4.2287407	25889.0947
Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	2467	2.2685749	11056.7379	8.088917	1982-1999 –> 2000-2012	4.2034596	25802.4346
Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + oxygen + silicate	2467	2.2889445	11096.8427	8.001514	1982-1999 –> 2000-2012	4.2995412	26110.5272
Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + oxygen + silicate + phosphate	2467	2.2789172	11077.1806	8.289996	1982-1999 –> 2000-2012	4.2138321	25820.9878
Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	2467	2.2789172	11077.1806	8.289996	1982-1999 –> 2000-2012	4.2138321	25820.9878
Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + silicate + phosphate	2467	2.2960807	11112.2014	8.241017	1982-1999 –> 2000-2012	4.3174106	26163.6125
Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	2467	2.2789172	11077.1806	8.289996	1982-1999 –> 2000-2012	4.2138321	25820.9878
Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + silicate + phosphate_star	2467	2.2794555	11076.3459	8.377885	1982-1999 –> 2000-2012	4.2146823	25819.2952
Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	1354	2.5167285	6355.8608	8.981549	2000-2012 –> 2013-2019	4.8166120	17480.2273
Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	1354	2.5592084	6399.1877	7.809156	2000-2012 –> 2013-2019	4.8384768	17477.1285
Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	1354	2.5449625	6384.0714	7.909962	2000-2012 –> 2013-2019	4.8395069	17494.9704
Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + oxygen + nitrate + silicate	1354	2.5600102	6400.0360	7.946110	2000-2012 –> 2013-2019	4.8295048	17456.7738
Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	1354	2.5387308	6379.4324	7.893914	2000-2012 –> 2013-2019	4.8073057	17436.1703
Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + oxygen + silicate	1354	2.5605769	6398.6353	7.837041	2000-2012 –> 2013-2019	4.8495214	17495.4780
Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + oxygen + silicate + phosphate	1354	2.5570958	6396.9513	7.768822	2000-2012 –> 2013-2019	4.8360130	17474.1319
Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	1354	2.5570958	6396.9513	7.768822	2000-2012 –> 2013-2019	4.8360130	17474.1319
Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	1354	2.5570958	6396.9513	7.768822	2000-2012 –> 2013-2019	4.8360130	17474.1319
Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + silicate + phosphate_star	1354	2.5592213	6397.2014	7.805707	2000-2012 –> 2013-2019	4.8386769	17473.5473
Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + nitrate + silicate + phosphate_star	1718	1.9562432	7193.1178	9.132845	1982-1999 –> 2000-2012	3.6029208	17000.7204
Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	1718	1.9538428	7190.8991	8.874628	1982-1999 –> 2000-2012	3.5863843	16956.4621
Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + oxygen + silicate + phosphate	1718	1.9539957	7189.1680	8.935869	1982-1999 –> 2000-2012	3.5950802	16979.3910
Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + oxygen + silicate + phosphate_star	1718	1.9539957	7189.1680	8.935869	1982-1999 –> 2000-2012	3.5950802	16979.3910
Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	1718	1.9539957	7189.1680	8.935869	1982-1999 –> 2000-2012	3.5950802	16979.3910
Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + silicate + phosphate_star	1718	1.9721486	7218.9416	9.184943	1982-1999 –> 2000-2012	3.6245406	17042.2394
Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	1718	1.8584907	7018.9844	10.244915	1982-1999 –> 2000-2012	3.4634280	16697.4390
Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + oxygen + silicate + phosphate	1718	1.8961482	7085.9101	9.843194	1982-1999 –> 2000-2012	3.6061027	17086.1201
Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	1718	1.8961482	7085.9101	9.843194	1982-1999 –> 2000-2012	3.6061027	17086.1201
Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	1718	1.8961482	7085.9101	9.843194	1982-1999 –> 2000-2012	3.6061027	17086.1201
Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + oxygen + silicate + phosphate	966	2.1032628	4189.8117	9.400527	2000-2012 –> 2013-2019	4.0572586	11378.9797
Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + oxygen + silicate + phosphate_star	966	2.1032628	4189.8117	9.400527	2000-2012 –> 2013-2019	4.0572586	11378.9797
Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	966	2.1032628	4189.8117	9.400527	2000-2012 –> 2013-2019	4.0572586	11378.9797
Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	966	2.1337577	4217.6223	10.080373	2000-2012 –> 2013-2019	4.0401750	11322.0907
Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	966	2.0835440	4173.6130	10.558303	2000-2012 –> 2013-2019	3.9420347	11192.5974
Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + oxygen + silicate + phosphate	966	2.1186866	4203.9279	10.028459	2000-2012 –> 2013-2019	4.0148349	11289.8379
Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	966	2.1186866	4203.9279	10.028459	2000-2012 –> 2013-2019	4.0148349	11289.8379
Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + silicate + phosphate	966	2.1199344	4203.0654	10.143714	2000-2012 –> 2013-2019	4.0665458	11377.2238
Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	966	2.1186866	4203.9279	10.028459	2000-2012 –> 2013-2019	4.0148349	11289.8379
Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + silicate + phosphate_star	966	2.1380335	4219.4899	10.237607	2000-2012 –> 2013-2019	4.0453228	11323.5296
Indo-Pacific	(28.1, Inf]	cstar_tref ~ oxygen + nitrate + silicate + phosphate_star	14627	1.4377139	52142.4183	8.775636	1982-1999 –> 2000-2012	2.6704554	118372.4597
Indo-Pacific	(28.1, Inf]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	14627	1.3133686	49498.1312	9.366747	1982-1999 –> 2000-2012	2.4634246	112906.4588
Indo-Pacific	(28.1, Inf]	cstar_tref ~ sal + oxygen + silicate + phosphate	14627	1.6072254	55402.9247	9.839246	1982-1999 –> 2000-2012	3.1030388	129499.7490
Indo-Pacific	(28.1, Inf]	cstar_tref ~ sal + oxygen + silicate + phosphate_star	14627	1.6072254	55402.9247	9.839246	1982-1999 –> 2000-2012	3.1030388	129499.7490
Indo-Pacific	(28.1, Inf]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	14627	1.6072254	55402.9247	9.839246	1982-1999 –> 2000-2012	3.1030388	129499.7490
Indo-Pacific	(28.1, Inf]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	14627	1.5669226	54659.9964	17.476131	1982-1999 –> 2000-2012	2.8961915	123956.0907
Indo-Pacific	(28.1, Inf]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	14627	1.3591254	50499.9672	12.907049	1982-1999 –> 2000-2012	2.4770940	112757.4410
Indo-Pacific	(28.1, Inf]	cstar_tref ~ temp + oxygen + silicate + phosphate	14627	1.6697977	56520.2259	18.221976	1982-1999 –> 2000-2012	3.1051732	128939.6726
Indo-Pacific	(28.1, Inf]	cstar_tref ~ temp + oxygen + silicate + phosphate_star	14627	1.6697977	56520.2259	18.221976	1982-1999 –> 2000-2012	3.1051732	128939.6726
Indo-Pacific	(28.1, Inf]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	14627	1.6697977	56520.2259	18.221976	1982-1999 –> 2000-2012	3.1051732	128939.6726
Indo-Pacific	(28.1, Inf]	cstar_tref ~ oxygen + nitrate + silicate + phosphate_star	7804	1.6277929	29763.4008	8.895830	2000-2012 –> 2013-2019	3.0655068	81905.8190
Indo-Pacific	(28.1, Inf]	cstar_tref ~ sal + nitrate + silicate + phosphate_star	7804	1.8010952	31342.4609	9.871618	2000-2012 –> 2013-2019	3.5052374	88458.2806
Indo-Pacific	(28.1, Inf]	cstar_tref ~ sal + oxygen + nitrate + phosphate_star	7804	1.8057962	31383.1458	9.458145	2000-2012 –> 2013-2019	3.5095128	88491.6590
Indo-Pacific	(28.1, Inf]	cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star	7804	1.4569278	28034.5504	8.140875	2000-2012 –> 2013-2019	2.7702964	77532.6816
Indo-Pacific	(28.1, Inf]	cstar_tref ~ sal + oxygen + silicate + phosphate	7804	1.7050926	30487.5244	9.446823	2000-2012 –> 2013-2019	3.3123180	85890.4491
Indo-Pacific	(28.1, Inf]	cstar_tref ~ sal + oxygen + silicate + phosphate_star	7804	1.7050926	30487.5244	9.446823	2000-2012 –> 2013-2019	3.3123180	85890.4491
Indo-Pacific	(28.1, Inf]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	7804	1.7050926	30487.5244	9.446823	2000-2012 –> 2013-2019	3.3123180	85890.4491
Indo-Pacific	(28.1, Inf]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	7804	1.7947239	31287.1497	15.290224	2000-2012 –> 2013-2019	3.3616465	85947.1460
Indo-Pacific	(28.1, Inf]	cstar_tref ~ temp + oxygen + nitrate + phosphate_star	7804	1.7788231	31148.2512	8.242671	2000-2012 –> 2013-2019	3.5042995	88628.0295
Indo-Pacific	(28.1, Inf]	cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star	7804	1.5743591	29244.4567	10.872611	2000-2012 –> 2013-2019	2.9334846	79744.4239

5.3.2 per fitting unit

lm_best %>%
  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	3.744823	1520.3182
Atlantic	(-Inf,26]	2000-2012 –> 2013-2019	3.712553	956.1927
Atlantic	(26,26.5]	1982-1999 –> 2000-2012	8.861389	14957.2252
Atlantic	(26,26.5]	2000-2012 –> 2013-2019	8.757643	9369.5854
Atlantic	(26.5,26.75]	1982-1999 –> 2000-2012	6.894776	17492.2030
Atlantic	(26.5,26.75]	2000-2012 –> 2013-2019	7.100524	11047.3232
Atlantic	(26.75,27]	1982-1999 –> 2000-2012	4.646521	24897.6716
Atlantic	(26.75,27]	2000-2012 –> 2013-2019	5.291436	17361.1639
Atlantic	(27,27.25]	1982-1999 –> 2000-2012	4.566950	20802.8036
Atlantic	(27,27.25]	2000-2012 –> 2013-2019	4.946596	14630.2659
Atlantic	(27.25,27.5]	1982-1999 –> 2000-2012	5.992368	24307.3549
Atlantic	(27.25,27.5]	2000-2012 –> 2013-2019	6.896762	16463.9211
Atlantic	(27.5,27.75]	1982-1999 –> 2000-2012	6.337729	33343.2852
Atlantic	(27.5,27.75]	2000-2012 –> 2013-2019	6.560443	22437.8248
Atlantic	(27.75,27.85]	1982-1999 –> 2000-2012	4.361288	10339.9903
Atlantic	(27.75,27.85]	2000-2012 –> 2013-2019	4.658339	7003.5670
Atlantic	(27.85,27.95]	1982-1999 –> 2000-2012	6.778146	12131.3448
Atlantic	(27.85,27.95]	2000-2012 –> 2013-2019	7.688924	8604.7387
Atlantic	(27.95,28.05]	1982-1999 –> 2000-2012	8.858264	16735.5597
Atlantic	(27.95,28.05]	2000-2012 –> 2013-2019	8.596203	10808.8887
Atlantic	(28.05,28.1]	1982-1999 –> 2000-2012	2.164397	5586.3259
Atlantic	(28.05,28.1]	2000-2012 –> 2013-2019	2.424822	3996.1155
Atlantic	(28.1,28.15]	1982-1999 –> 2000-2012	1.600964	5257.1584
Atlantic	(28.1,28.15]	2000-2012 –> 2013-2019	1.952076	3927.1699
Atlantic	(28.15,28.2]	1982-1999 –> 2000-2012	1.285315	7359.0503
Atlantic	(28.15,28.2]	2000-2012 –> 2013-2019	1.532650	5129.0699
Atlantic	(28.2, Inf]	1982-1999 –> 2000-2012	1.149504	10828.4603
Atlantic	(28.2, Inf]	2000-2012 –> 2013-2019	1.252790	7772.3459
Indo-Pacific	(-Inf,26]	1982-1999 –> 2000-2012	15.547857	86706.1521
Indo-Pacific	(-Inf,26]	2000-2012 –> 2013-2019	14.917582	55659.6333
Indo-Pacific	(26,26.5]	1982-1999 –> 2000-2012	10.175897	80509.5692
Indo-Pacific	(26,26.5]	2000-2012 –> 2013-2019	10.232494	51590.9333
Indo-Pacific	(26.5,26.75]	1982-1999 –> 2000-2012	8.325310	61608.6240
Indo-Pacific	(26.5,26.75]	2000-2012 –> 2013-2019	8.774905	40348.6602
Indo-Pacific	(26.75,27]	1982-1999 –> 2000-2012	9.342379	75939.7426
Indo-Pacific	(26.75,27]	2000-2012 –> 2013-2019	10.231121	49976.1847
Indo-Pacific	(27,27.25]	1982-1999 –> 2000-2012	8.606340	87425.3970
Indo-Pacific	(27,27.25]	2000-2012 –> 2013-2019	10.214295	60083.9374
Indo-Pacific	(27.25,27.5]	1982-1999 –> 2000-2012	6.955167	72753.9838
Indo-Pacific	(27.25,27.5]	2000-2012 –> 2013-2019	8.228304	50372.8190
Indo-Pacific	(27.5,27.75]	1982-1999 –> 2000-2012	5.958458	68810.4735
Indo-Pacific	(27.5,27.75]	2000-2012 –> 2013-2019	7.681217	48856.9948
Indo-Pacific	(27.75,27.85]	1982-1999 –> 2000-2012	4.715905	24400.0604
Indo-Pacific	(27.75,27.85]	2000-2012 –> 2013-2019	5.921523	17778.0768
Indo-Pacific	(27.85,27.95]	1982-1999 –> 2000-2012	4.914581	31694.5511
Indo-Pacific	(27.85,27.95]	2000-2012 –> 2013-2019	6.075558	22251.6915
Indo-Pacific	(27.95,28.05]	1982-1999 –> 2000-2012	4.242905	25921.1764
Indo-Pacific	(27.95,28.05]	2000-2012 –> 2013-2019	4.832764	17473.6691
Indo-Pacific	(28.05,28.1]	1982-1999 –> 2000-2012	3.588082	16989.3394
Indo-Pacific	(28.05,28.1]	2000-2012 –> 2013-2019	4.031036	11322.1894
Indo-Pacific	(28.1, Inf]	1982-1999 –> 2000-2012	2.913180	124331.0715
Indo-Pacific	(28.1, Inf]	2000-2012 –> 2013-2019	3.258694	84837.9387

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 = ""))

6 Model diagnotics

6.1 Selection criterion vs predictors

The selection criterion (aic) was plotted against the number of predictors (limited to 3 - 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
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
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
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
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")
  )

	nitrate	oxygen	phosphate	phosphate_star	sal	silicate	temp
Atlantic - 1982-1999 --> 2000-2012
(-Inf,26]	6	5	3	6	10	8	0
(26,26.5]	5	7	2	6	0	8	10
(26.5,26.75]	4	7	4	8	4	9	6
(26.75,27]	4	6	4	8	5	10	5
(27,27.25]	4	7	4	7	5	10	5
(27.25,27.5]	6	7	2	7	1	8	9
(27.5,27.75]	6	7	2	8	2	8	8
(27.75,27.85]	6	7	2	8	2	8	8
(27.85,27.95]	4	6	4	8	0	5	10
(27.95,28.05]	6	8	3	8	2	7	7
(28.05,28.1]	4	7	4	7	0	6	10
(28.1,28.15]	4	6	4	8	0	5	10
(28.15,28.2]	10	8	0	6	4	5	6
(28.2, Inf]	9	8	0	7	6	7	4
total	78.00	96.00	38.00	102.00	41.00	104.00	98.00
Atlantic - 2000-2012 --> 2013-2019
(-Inf,26]	6	5	3	6	10	8	0
(26,26.5]	5	6	3	6	0	8	10
(26.5,26.75]	4	6	4	8	5	10	5
(26.75,27]	4	6	4	8	5	10	5
(27,27.25]	4	6	4	7	6	10	4
(27.25,27.5]	7	8	2	6	6	8	4
(27.5,27.75]	4	7	3	8	9	7	1
(27.75,27.85]	4	7	3	8	7	9	3
(27.85,27.95]	4	6	4	8	0	5	10
(27.95,28.05]	4	8	4	8	2	5	8
(28.05,28.1]	4	7	4	7	0	6	10
(28.1,28.15]	4	6	4	8	0	5	10
(28.15,28.2]	10	8	0	6	4	5	6
(28.2, Inf]	8	7	1	7	8	7	2
total	72.00	93.00	43.00	101.00	62.00	103.00	78.00
Indo-Pacific - 1982-1999 --> 2000-2012
(-Inf,26]	4	6	4	8	0	5	10
(26,26.5]	3	6	4	8	0	6	10
(26.5,26.75]	4	6	3	8	0	6	10
(26.75,27]	4	6	4	8	0	5	10
(27,27.25]	5	7	3	8	1	6	9
(27.25,27.5]	3	8	5	7	4	10	4
(27.5,27.75]	6	6	2	9	8	8	1
(27.75,27.85]	3	6	4	8	0	6	10
(27.85,27.95]	6	7	3	6	0	6	10
(27.95,28.05]	4	6	3	6	0	9	10
(28.05,28.1]	3	6	4	8	6	10	4
(28.1, Inf]	4	7	4	8	4	10	5
total	49.00	77.00	43.00	92.00	23.00	87.00	93.00
Indo-Pacific - 2000-2012 --> 2013-2019
(-Inf,26]	3	6	4	8	0	6	10
(26,26.5]	4	5	3	9	0	6	10
(26.5,26.75]	4	6	3	8	0	6	10
(26.75,27]	3	8	5	7	2	6	8
(27,27.25]	4	7	4	8	4	9	6
(27.25,27.5]	3	8	5	7	4	10	4
(27.5,27.75]	7	7	2	8	7	8	2
(27.75,27.85]	4	5	4	8	6	10	4
(27.85,27.95]	10	7	0	6	3	6	7
(27.95,28.05]	5	7	2	7	1	9	9
(28.05,28.1]	2	5	5	7	3	10	7
(28.1, Inf]	7	7	2	9	6	8	3
total	56.00	78.00	39.00	92.00	36.00	94.00	80.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
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
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