eMLR - model fitting

workflowr

• Summary
• Checks
• Past versions

Last updated: 2021-02-27

Checks: 7 0

Knit directory: emlr_mod_v_XXX/

This reproducible R Markdown analysis was created with workflowr (version 1.6.2). The Checks tab describes the reproducibility checks that were applied when the results were created. The Past versions tab lists the development history.

---

R Markdown file: up-to-date

Great! Since the R Markdown file has been committed to the Git repository, you know the exact version of the code that produced these results.

Environment: empty

Great job! The global environment was empty. Objects defined in the global environment can affect the analysis in your R Markdown file in unknown ways. For reproduciblity it’s best to always run the code in an empty environment.

Seed: set.seed(20200707)

The command set.seed(20200707) was run prior to running the code in the R Markdown file. Setting a seed ensures that any results that rely on randomness, e.g. subsampling or permutations, are reproducible.

Session information: recorded

Great job! Recording the operating system, R version, and package versions is critical for reproducibility.

Cache: none

Nice! There were no cached chunks for this analysis, so you can be confident that you successfully produced the results during this run.

File paths: relative

Great job! Using relative paths to the files within your workflowr project makes it easier to run your code on other machines.

Repository version: 764a9c2

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 764a9c2. 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
Rmd	764a9c2	Donghe-Zhu	2021-02-27	local rebuild
html	54d5b5b	Donghe-Zhu	2021-02-27	Build site.
Rmd	2e836bf	Donghe-Zhu	2021-02-27	local rebuild
html	330f064	Donghe-Zhu	2021-02-27	Build site.
Rmd	f360282	Donghe-Zhu	2021-02-27	local rebuild
html	adbc9bc	Donghe-Zhu	2021-02-27	Build site.
Rmd	a8a8f9c	Donghe-Zhu	2021-02-27	local rebuild
html	5937141	Donghe-Zhu	2021-02-27	Build site.
Rmd	5f0bfca	Donghe-Zhu	2021-02-27	local rebuild
html	4414bbf	Donghe-Zhu	2021-02-27	Build site.
html	a265efb	Donghe-Zhu	2021-02-27	Build site.
html	19edd1e	Donghe-Zhu	2021-02-27	Build site.
Rmd	1ae5bb4	Donghe-Zhu	2021-02-27	local rebuild
html	f20483f	Donghe-Zhu	2021-02-26	Build site.
html	6a2c7b3	Donghe-Zhu	2021-02-25	Build site.
html	02b976d	Donghe-Zhu	2021-02-24	Build site.
html	354c224	Donghe-Zhu	2021-02-24	Build site.
Rmd	d910d39	Donghe-Zhu	2021-02-24	local rebuild
html	1a0a88a	Donghe-Zhu	2021-02-24	Build site.
Rmd	663e891	Donghe-Zhu	2021-02-24	local rebuild
html	57f701e	Donghe-Zhu	2021-02-24	Build site.
html	06f3149	Donghe-Zhu	2021-02-16	Build site.
html	401eab3	Donghe-Zhu	2021-02-15	Build site.
html	e3bba84	Donghe-Zhu	2021-02-15	Build site.
html	5dce4b1	Donghe-Zhu	2021-02-15	Build site.
Rmd	35b0f2e	Donghe-Zhu	2021-02-15	local rebuild
html	4469a0c	Donghe-Zhu	2021-02-13	Build site.
Rmd	8186d57	Donghe-Zhu	2021-02-10	apply nr_obs
html	5ae6a69	Donghe-Zhu	2021-02-10	Build site.
Rmd	eeb6557	Donghe-Zhu	2021-02-10	local rebuild
html	05385dc	Donghe-Zhu	2021-02-10	Build site.
html	f791ae4	Donghe-Zhu	2021-02-09	Build site.
html	f71ae34	Donghe-Zhu	2021-02-09	Build site.
html	c011832	Donghe-Zhu	2021-02-09	Build site.
html	a145fa7	Donghe-Zhu	2021-02-09	Build site.
Rmd	37a41c8	Donghe-Zhu	2021-02-09	local rebuild
html	c344e42	Donghe-Zhu	2021-02-08	Build site.
Rmd	8f5fa79	Donghe-Zhu	2021-02-08	local rebuild
html	2f095d7	Donghe-Zhu	2021-02-07	Build site.
html	2305044	Donghe-Zhu	2021-02-07	Build site.
Rmd	c3638c1	Donghe-Zhu	2021-02-07	local rebuild
html	1fad5f1	Donghe-Zhu	2021-02-07	Build site.
html	ca03c39	Donghe-Zhu	2021-02-07	Build site.
html	e2ffc14	Donghe-Zhu	2021-02-05	Build site.
html	cd7c52c	Donghe-Zhu	2021-02-04	Build site.
Rmd	bcf91a8	jens-daniel-mueller	2021-02-04	separate MLR tables, added plot residuals vs location
html	bcf84f4	Donghe-Zhu	2021-02-02	Build site.
html	a518739	Donghe-Zhu	2021-02-01	Build site.
html	61666de	Donghe-Zhu	2021-01-31	Build site.
html	865b582	Donghe-Zhu	2021-01-31	Build site.
html	3e68089	Donghe-Zhu	2021-01-31	Build site.
html	ecf335c	Donghe-Zhu	2021-01-31	Build site.
html	a618965	Donghe-Zhu	2021-01-31	Build site.
html	59e006e	Donghe-Zhu	2021-01-31	Build site.
html	a1c8f87	Donghe-Zhu	2021-01-31	Build site.
html	ae5c18f	Donghe-Zhu	2021-01-31	Build site.
html	b50fe52	Donghe-Zhu	2021-01-31	Build site.
Rmd	ac99ae5	jens-daniel-mueller	2021-01-29	code review
html	ac99ae5	jens-daniel-mueller	2021-01-29	code review
html	b5bdcaf	Donghe-Zhu	2021-01-29	Build site.
Rmd	b234505	Donghe-Zhu	2021-01-29	MLR approach across all depth
html	442010d	Donghe-Zhu	2021-01-29	Build site.
Rmd	e67e7dd	Donghe-Zhu	2021-01-29	surface equilibrium approach across all latitudes irrespective of gamma
html	372adf5	Donghe-Zhu	2021-01-29	Build site.
html	af8788e	Donghe-Zhu	2021-01-29	Build site.
html	21c91c9	Donghe-Zhu	2021-01-29	Build site.
html	eded038	Donghe-Zhu	2021-01-29	Build site.
html	541d4dd	Donghe-Zhu	2021-01-29	Build site.
html	6a75576	Donghe-Zhu	2021-01-28	Build site.
html	16fba40	Donghe-Zhu	2021-01-28	Build site.
Rmd	aecf0c6	Donghe-Zhu	2021-01-28	diagnostic ploting
Rmd	a02684e	Donghe-Zhu	2021-01-28	error
Rmd	4584be9	Donghe-Zhu	2021-01-28	latest
html	12bc567	Donghe-Zhu	2021-01-27	Build site.
html	ceed31b	Donghe-Zhu	2021-01-27	Build site.
html	342402d	Donghe-Zhu	2021-01-27	Build site.
html	5bad5c2	Donghe-Zhu	2021-01-27	Build site.
Rmd	c2c9529	Donghe-Zhu	2021-01-27	random subsetting based on lat
html	61efb56	Donghe-Zhu	2021-01-25	Build site.
html	48f638e	Donghe-Zhu	2021-01-25	Build site.
html	c1cec47	Donghe-Zhu	2021-01-25	Build site.
html	05ffb0c	Donghe-Zhu	2021-01-25	Build site.
html	8b97165	Donghe-Zhu	2021-01-25	Build site.
html	c569946	Donghe-Zhu	2021-01-24	Build site.
html	a2f0d56	Donghe-Zhu	2021-01-23	Build site.
html	28509fc	Donghe-Zhu	2021-01-23	Build site.
html	4c28e4a	Donghe-Zhu	2021-01-22	Build site.
html	24cc264	jens-daniel-mueller	2021-01-22	cleaned /docs before creating copies
html	88eb28f	Donghe-Zhu	2021-01-21	Build site.
html	2679490	Donghe-Zhu	2021-01-21	Build site.
html	7891955	Donghe-Zhu	2021-01-21	Build site.
html	d4cf1cb	Donghe-Zhu	2021-01-21	Build site.
Rmd	167eeec	Donghe-Zhu	2021-01-21	surface DIC calculation with atmospheric equilibrium option
html	1f3e5b6	jens-daniel-mueller	2021-01-20	Build site.
html	0e7bdf1	jens-daniel-mueller	2021-01-15	cleaning template repository
html	73cbef3	jens-daniel-mueller	2021-01-15	Build site.
html	4571843	jens-daniel-mueller	2021-01-14	revision and html deleted for template copying
html	23151cd	jens-daniel-mueller	2021-01-14	Build site.
html	b3564aa	jens-daniel-mueller	2021-01-14	Build site.
html	8d032c3	jens-daniel-mueller	2021-01-14	Build site.
html	022871c	Donghe-Zhu	2021-01-13	Build site.
Rmd	d44f36f	Donghe-Zhu	2021-01-13	reorder analysis final
html	17dee1d	jens-daniel-mueller	2021-01-13	Build site.
Rmd	9e04fd7	jens-daniel-mueller	2021-01-13	local rebuild after revision
html	a076226	Donghe-Zhu	2021-01-11	Build site.
Rmd	52eff18	Donghe-Zhu	2021-01-09	Implemet model_run and subsetting
html	7cdea0c	jens-daniel-mueller	2021-01-06	Build site.
Rmd	b5934dd	jens-daniel-mueller	2021-01-06	local rebuild after revision
html	fa85b93	jens-daniel-mueller	2021-01-06	Build site.
html	e5cb81a	Donghe-Zhu	2021-01-05	Build site.
Rmd	608cc45	Donghe-Zhu	2021-01-05	modification of analysis
html	a499f10	Donghe-Zhu	2021-01-05	Build site.
Rmd	715bdb4	Donghe-Zhu	2021-01-02	model modification
html	fb8a752	Donghe-Zhu	2020-12-23	Build site.
Rmd	82e3c9c	Donghe-Zhu	2020-12-23	first build after creating model template
html	8fae0b2	Donghe-Zhu	2020-12-21	Build site.
Rmd	00a1322	Donghe-Zhu	2020-12-21	first build after creating model template
Rmd	d73ae35	Donghe-Zhu	2020-12-21	first version with lm error
html	c8b76b3	jens-daniel-mueller	2020-12-19	Build site.
Rmd	b5fedce	jens-daniel-mueller	2020-12-19	first build after creating model template
Rmd	8e8abf5	Jens Müller	2020-12-18	Initial commit

---

1 Required data

Required are:

• cleaned and prepared GLODAP-based synthetic model subsetting file

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

2 Predictor combinations

Find all possible combinations of following considered predictor variables:

• sal, temp, oxygen, aou, nitrate, silicate, silicate_star, phosphate, phosphate_star

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

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

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

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

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

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

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

# convert to tibble
lm_all <- as_tibble(lm_all)

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

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

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

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

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


#remove temp sal predictor combination
lm_all <- lm_all %>%
  mutate(lm_coeff_filter = str_remove(lm_coeff, "phosphate_star")) %>%
  filter(!(
    str_detect(lm_coeff_filter, "nitrate") &
      str_detect(lm_coeff_filter, "phosphate")
  )) %>%
  filter(!(
    str_detect(lm_coeff_filter, "temp") &
      str_detect(lm_coeff_filter, "sal")
  )) %>%
  filter(!(
    str_detect(lm_coeff_filter, "oxygen") &
      str_detect(lm_coeff_filter, "aou")
  )) %>%
  select(-lm_coeff_filter)

lm_rm <- bind_rows(lm_rm_ph, lm_rm_si) %>%
  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:

• 111

4 Fit all models

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

• cstar_tref

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

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

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

5 Prepare coeffcients

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

5.1 Formatting

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

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

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

5.2 Predictor selection

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

• 10

The criterion used to select the best models was:

• rmse

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

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

# remove models with predictors fitted as NA

lm_all_fitted_wide <- lm_all_fitted_wide %>%
  filter(na_predictor == FALSE)

# calculate RMSE sum for adjacent eras
lm_all_fitted_wide_eras <- lm_all_fitted_wide  %>%
  select(basin, gamma_slab, 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 + aou + nitrate + silicate + phosphate_star	158	1.5921923	609.3599	6.118408	1982-1999 –> 2000-2012	3.2779521	1430.8804
Atlantic	(-Inf,26]	cstar_tref ~ sal + aou + nitrate + silicate_star + phosphate_star	158	1.5921923	609.3599	6.118408	1982-1999 –> 2000-2012	3.2779521	1430.8804
Atlantic	(-Inf,26]	cstar_tref ~ sal + aou + phosphate + phosphate_star	158	1.0492038	475.5626	4.357844	1982-1999 –> 2000-2012	2.3855522	1198.4564
Atlantic	(-Inf,26]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	158	1.0480832	477.2249	4.299021	1982-1999 –> 2000-2012	2.3837757	1201.9145
Atlantic	(-Inf,26]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	158	1.0475856	477.0748	4.332795	1982-1999 –> 2000-2012	2.3832111	1201.7435
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 + nitrate + silicate_star + phosphate_star	158	1.6216413	615.1512	8.073641	1982-1999 –> 2000-2012	3.2856343	1431.2652
Atlantic	(-Inf,26]	cstar_tref ~ temp + aou + phosphate + phosphate_star	158	1.4654213	581.1417	7.417048	1982-1999 –> 2000-2012	3.1276914	1394.8248
Atlantic	(-Inf,26]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	158	1.4545749	580.7941	7.074044	1982-1999 –> 2000-2012	3.1160634	1396.2815
Atlantic	(-Inf,26]	cstar_tref ~ temp + aou + silicate_star + phosphate + phosphate_star	158	1.4627166	582.5579	7.363430	1982-1999 –> 2000-2012	3.1249791	1398.2391
Atlantic	(-Inf,26]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	72	1.6825475	293.2516	3.600013	2000-2012 –> 2013-2019	3.2747398	902.6116
Atlantic	(-Inf,26]	cstar_tref ~ sal + aou + nitrate + silicate_star + phosphate_star	72	1.6825475	293.2516	3.600013	2000-2012 –> 2013-2019	3.2747398	902.6116
Atlantic	(-Inf,26]	cstar_tref ~ sal + aou + phosphate + phosphate_star	72	1.2725075	251.0296	3.658367	2000-2012 –> 2013-2019	2.3217113	726.5922
Atlantic	(-Inf,26]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	72	1.1936365	243.8158	3.237974	2000-2012 –> 2013-2019	2.2417197	721.0407
Atlantic	(-Inf,26]	cstar_tref ~ sal + aou + silicate + phosphate_star	72	1.6918467	292.0453	4.133504	2000-2012 –> 2013-2019	3.2851805	899.6317
Atlantic	(-Inf,26]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	72	1.2456464	249.9574	3.695592	2000-2012 –> 2013-2019	2.2932320	727.0322
Atlantic	(-Inf,26]	cstar_tref ~ sal + aou + silicate_star + phosphate_star	72	1.6838402	291.3622	3.802365	2000-2012 –> 2013-2019	3.3019542	903.8254
Atlantic	(-Inf,26]	cstar_tref ~ temp + aou + phosphate + phosphate_star	72	1.5761474	281.8448	5.713265	2000-2012 –> 2013-2019	3.0415687	862.9865
Atlantic	(-Inf,26]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	72	1.5044066	277.1365	5.003414	2000-2012 –> 2013-2019	2.9589816	857.9306
Atlantic	(-Inf,26]	cstar_tref ~ temp + aou + silicate_star + phosphate + phosphate_star	72	1.5517609	281.5994	5.508599	2000-2012 –> 2013-2019	3.0144775	864.1573
Atlantic	(26,26.5]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1046	3.5518213	5633.9468	13.477555	1982-1999 –> 2000-2012	6.9985663	13749.8718
Atlantic	(26,26.5]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	1046	3.8736280	5815.3881	13.955824	1982-1999 –> 2000-2012	7.6071748	14175.0946
Atlantic	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate	1046	4.0405761	5901.6616	14.579174	1982-1999 –> 2000-2012	8.1249524	14533.2897
Atlantic	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	1046	4.0338370	5900.1695	15.685804	1982-1999 –> 2000-2012	8.1154929	14531.7655
Atlantic	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate_star	1046	4.0405761	5901.6616	14.579174	1982-1999 –> 2000-2012	8.1249524	14533.2897
Atlantic	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	1046	4.0338370	5900.1695	15.685804	1982-1999 –> 2000-2012	8.1154929	14531.7655
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 + nitrate + silicate_star	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_star + phosphate_star	1046	3.9825944	5873.4242	15.972464	1982-1999 –> 2000-2012	8.0211559	14472.6469
Atlantic	(26,26.5]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	567	3.6070302	3077.8676	10.902628	2000-2012 –> 2013-2019	7.1588515	8711.8144
Atlantic	(26,26.5]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	567	3.9979703	3194.5585	12.046814	2000-2012 –> 2013-2019	7.8715984	9009.9467
Atlantic	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate	567	3.8959809	3163.2545	13.662757	2000-2012 –> 2013-2019	7.9365570	9064.9161
Atlantic	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	567	3.8631919	3155.6702	14.511046	2000-2012 –> 2013-2019	7.8970289	9055.8397
Atlantic	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate_star	567	3.8959809	3163.2545	13.662757	2000-2012 –> 2013-2019	7.9365570	9064.9161
Atlantic	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	567	3.8631919	3155.6702	14.511046	2000-2012 –> 2013-2019	7.8970289	9055.8397
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 + nitrate + silicate_star	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_star + phosphate_star	567	3.8005103	3137.1198	14.416820	2000-2012 –> 2013-2019	7.7831048	9010.5440
Atlantic	(26.5,26.75]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	1397	3.3429532	7348.4661	10.697715	1982-1999 –> 2000-2012	6.4415667	17034.7612
Atlantic	(26.5,26.75]	cstar_tref ~ sal + aou + silicate + phosphate	1397	3.3316510	7339.0038	10.073920	1982-1999 –> 2000-2012	6.4785991	17084.0238
Atlantic	(26.5,26.75]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1397	3.2700730	7288.8799	10.507891	1982-1999 –> 2000-2012	6.3292231	16928.5449
Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	1397	3.3250708	7335.4800	10.297171	1982-1999 –> 2000-2012	6.4698137	17079.8406
Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	1397	3.3250708	7335.4800	10.297171	1982-1999 –> 2000-2012	6.4698137	17079.8406
Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + silicate + phosphate	1397	3.4027533	7398.0044	10.818893	1982-1999 –> 2000-2012	6.5639402	17160.1523
Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1397	3.2249399	7250.0489	10.037456	1982-1999 –> 2000-2012	6.2079736	16794.1190
Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + silicate + phosphate_star	1397	3.4336491	7423.2585	10.906572	1982-1999 –> 2000-2012	6.6279234	17224.9109
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 + nitrate + silicate_star + phosphate_star	1397	3.3402414	7348.1986	9.940278	1982-1999 –> 2000-2012	6.4993586	17109.8616
Atlantic	(26.5,26.75]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	656	3.4323365	3491.6598	11.340269	2000-2012 –> 2013-2019	6.7752897	10840.1259
Atlantic	(26.5,26.75]	cstar_tref ~ sal + aou + silicate + phosphate	656	3.3850399	3473.4551	10.375376	2000-2012 –> 2013-2019	6.7166909	10812.4589
Atlantic	(26.5,26.75]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	656	3.3455096	3460.0435	10.735554	2000-2012 –> 2013-2019	6.6155826	10748.9234
Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	656	3.4093443	3484.8416	9.603997	2000-2012 –> 2013-2019	6.7344151	10820.3216
Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	656	3.4093443	3484.8416	9.603997	2000-2012 –> 2013-2019	6.7344151	10820.3216
Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + silicate + phosphate	656	3.4910702	3513.9207	11.381386	2000-2012 –> 2013-2019	6.8938235	10911.9252
Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	656	3.3238240	3451.5114	11.113656	2000-2012 –> 2013-2019	6.5487639	10701.5603
Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + silicate + phosphate_star	656	3.5219283	3525.4667	11.414557	2000-2012 –> 2013-2019	6.9555774	10948.7252
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 + nitrate + silicate_star + phosphate_star	656	3.4254447	3491.0228	9.732422	2000-2012 –> 2013-2019	6.7656861	10839.2215
Atlantic	(26.75,27]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	2398	2.0935628	10360.8368	14.739529	1982-1999 –> 2000-2012	3.7661002	22513.8376
Atlantic	(26.75,27]	cstar_tref ~ sal + aou + silicate + phosphate	2398	2.1776482	10549.6945	20.348896	1982-1999 –> 2000-2012	4.1329846	23683.7595
Atlantic	(26.75,27]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2398	2.0001820	10143.9994	14.081577	1982-1999 –> 2000-2012	3.6531709	22225.1681
Atlantic	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	2398	2.1989701	10598.4250	16.309198	1982-1999 –> 2000-2012	3.9946362	23199.5206
Atlantic	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	2398	2.1989701	10598.4250	16.309198	1982-1999 –> 2000-2012	3.9946362	23199.5206
Atlantic	(26.75,27]	cstar_tref ~ temp + aou + silicate + phosphate	2398	2.2004206	10599.5875	15.280820	1982-1999 –> 2000-2012	3.9880121	23170.3800
Atlantic	(26.75,27]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	2398	2.0860637	10345.6268	19.975942	1982-1999 –> 2000-2012	3.5995129	21872.9381
Atlantic	(26.75,27]	cstar_tref ~ temp + aou + silicate + phosphate_star	2398	2.2716861	10752.4544	17.733985	1982-1999 –> 2000-2012	4.1106467	23501.1674
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 + nitrate + silicate_star + phosphate_star	2398	2.2396447	10686.3266	20.763751	1982-1999 –> 2000-2012	4.0590677	23369.9626
Atlantic	(26.75,27]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	1263	2.6364762	6045.0524	27.471250	2000-2012 –> 2013-2019	4.7300390	16405.8891
Atlantic	(26.75,27]	cstar_tref ~ sal + aou + silicate + phosphate	1263	2.5846959	5994.9482	24.744636	2000-2012 –> 2013-2019	4.7623441	16544.6427
Atlantic	(26.75,27]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1263	2.5330762	5945.9902	25.891764	2000-2012 –> 2013-2019	4.5332582	16089.9896
Atlantic	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	1263	2.6623762	6071.7460	22.596317	2000-2012 –> 2013-2019	4.8613463	16670.1710
Atlantic	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	1263	2.6623762	6071.7460	22.596317	2000-2012 –> 2013-2019	4.8613463	16670.1710
Atlantic	(26.75,27]	cstar_tref ~ temp + aou + silicate + phosphate	1263	2.7633578	6163.7827	28.240304	2000-2012 –> 2013-2019	4.9637784	16763.3702
Atlantic	(26.75,27]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1263	2.2552518	5652.5376	22.583567	2000-2012 –> 2013-2019	4.3413155	15998.1644
Atlantic	(26.75,27]	cstar_tref ~ temp + aou + silicate + phosphate_star	1263	2.8094508	6205.5690	28.490393	2000-2012 –> 2013-2019	5.0811370	16958.0234
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 + nitrate + silicate_star + phosphate_star	1263	2.6917313	6099.4450	22.955194	2000-2012 –> 2013-2019	4.9313760	16785.7717
Atlantic	(27,27.25]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	2031	2.0909522	8773.9389	9.715574	1982-1999 –> 2000-2012	4.3252336	20395.7044
Atlantic	(27,27.25]	cstar_tref ~ sal + aou + nitrate + silicate_star + phosphate_star	2031	2.0909522	8773.9389	9.715574	1982-1999 –> 2000-2012	4.3252336	20395.7044
Atlantic	(27,27.25]	cstar_tref ~ sal + aou + silicate + phosphate	2031	1.9148200	8414.4995	9.235730	1982-1999 –> 2000-2012	3.6311952	18657.1971
Atlantic	(27,27.25]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2031	1.8904789	8364.5325	9.689213	1982-1999 –> 2000-2012	3.4604645	18143.6984
Atlantic	(27,27.25]	cstar_tref ~ sal + aou + silicate_star + phosphate	2031	2.1705611	8923.7203	10.608147	1982-1999 –> 2000-2012	4.0349174	19598.2839
Atlantic	(27,27.25]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	2031	2.1677153	8920.3911	10.441195	1982-1999 –> 2000-2012	3.9875628	19470.7749
Atlantic	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	2031	2.0594140	8712.2044	13.384239	1982-1999 –> 2000-2012	3.7359927	18834.3973
Atlantic	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	2031	2.0594140	8712.2044	13.384239	1982-1999 –> 2000-2012	3.7359927	18834.3973
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 + nitrate + silicate_star + phosphate_star	2031	2.1067869	8804.5844	12.394761	1982-1999 –> 2000-2012	3.8168590	19030.0700
Atlantic	(27,27.25]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	1142	2.2408412	5097.7041	11.555324	2000-2012 –> 2013-2019	4.3317934	13871.6430
Atlantic	(27,27.25]	cstar_tref ~ sal + aou + nitrate + silicate_star + phosphate_star	1142	2.2408412	5097.7041	11.555324	2000-2012 –> 2013-2019	4.3317934	13871.6430
Atlantic	(27,27.25]	cstar_tref ~ sal + aou + silicate + phosphate	1142	2.3003827	5155.6000	13.777336	2000-2012 –> 2013-2019	4.2152027	13570.0995
Atlantic	(27,27.25]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1142	2.2946948	5151.9457	14.080782	2000-2012 –> 2013-2019	4.1851737	13516.4782
Atlantic	(27,27.25]	cstar_tref ~ sal + aou + silicate + phosphate_star	1142	2.3424794	5197.0192	15.033313	2000-2012 –> 2013-2019	4.4776528	14053.9689
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 + nitrate + silicate_star + 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.25,27.5]	cstar_tref ~ sal + aou + nitrate + silicate_star + phosphate_star	2025	2.4514644	9392.2776	12.813379	1982-1999 –> 2000-2012	4.9864330	22927.8111
Atlantic	(27.25,27.5]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2025	2.6905754	9769.2091	21.422922	1982-1999 –> 2000-2012	4.8576234	22402.1133
Atlantic	(27.25,27.5]	cstar_tref ~ temp + aou + nitrate + silicate	2025	2.4063497	9315.0504	11.614926	1982-1999 –> 2000-2012	4.1264292	20616.3452
Atlantic	(27.25,27.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	2025	2.2547324	9053.4778	15.509020	1982-1999 –> 2000-2012	3.9706310	20342.7650
Atlantic	(27.25,27.5]	cstar_tref ~ temp + aou + nitrate + silicate_star	2025	2.4063497	9315.0504	11.614926	1982-1999 –> 2000-2012	4.1264292	20616.3452
Atlantic	(27.25,27.5]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	2025	2.2547324	9053.4778	15.509020	1982-1999 –> 2000-2012	3.9706310	20342.7650
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 + nitrate + silicate_star	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_star + phosphate_star	2025	2.3326197	9191.0186	10.601623	1982-1999 –> 2000-2012	4.1616629	20847.8613
Atlantic	(27.25,27.5]	cstar_tref ~ sal + aou + nitrate + silicate	1109	3.0025580	5597.8181	16.245244	2000-2012 –> 2013-2019	5.5311321	15113.5242
Atlantic	(27.25,27.5]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	1109	2.7952620	5441.1451	16.992150	2000-2012 –> 2013-2019	5.2467264	14833.4228
Atlantic	(27.25,27.5]	cstar_tref ~ sal + aou + nitrate + silicate_star	1109	3.0025580	5597.8181	16.245244	2000-2012 –> 2013-2019	5.5311321	15113.5242
Atlantic	(27.25,27.5]	cstar_tref ~ sal + aou + nitrate + silicate_star + phosphate_star	1109	2.7952620	5441.1451	16.992150	2000-2012 –> 2013-2019	5.2467264	14833.4228
Atlantic	(27.25,27.5]	cstar_tref ~ temp + aou + nitrate + silicate	1109	3.1081847	5674.5038	12.222568	2000-2012 –> 2013-2019	5.5145344	14989.5543
Atlantic	(27.25,27.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	1109	2.7923328	5438.8197	16.610907	2000-2012 –> 2013-2019	5.0470652	14492.2975
Atlantic	(27.25,27.5]	cstar_tref ~ temp + aou + nitrate + silicate_star	1109	3.1081847	5674.5038	12.222568	2000-2012 –> 2013-2019	5.5145344	14989.5543
Atlantic	(27.25,27.5]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	1109	2.7923328	5438.8197	16.610907	2000-2012 –> 2013-2019	5.0470652	14492.2975
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.25,27.5]	cstar_tref ~ temp + oxygen + nitrate + silicate_star + 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 + aou + nitrate + silicate + phosphate_star	2804	2.8239152	13793.2083	15.894239	1982-1999 –> 2000-2012	5.5291200	31835.8277
Atlantic	(27.5,27.75]	cstar_tref ~ sal + aou + nitrate + silicate_star + phosphate_star	2804	2.8239152	13793.2083	15.894239	1982-1999 –> 2000-2012	5.5291200	31835.8277
Atlantic	(27.5,27.75]	cstar_tref ~ sal + aou + phosphate	2804	2.7836924	13708.7557	13.816311	1982-1999 –> 2000-2012	5.2587422	31083.3422
Atlantic	(27.5,27.75]	cstar_tref ~ sal + aou + phosphate + phosphate_star	2804	2.7656318	13674.2524	14.269310	1982-1999 –> 2000-2012	5.0314996	30391.3946
Atlantic	(27.5,27.75]	cstar_tref ~ sal + aou + silicate + phosphate	2804	2.4538915	13003.5694	13.901382	1982-1999 –> 2000-2012	4.8210764	30047.3887
Atlantic	(27.5,27.75]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2804	2.3944064	12867.9502	14.685866	1982-1999 –> 2000-2012	4.4949743	29021.3907
Atlantic	(27.5,27.75]	cstar_tref ~ sal + aou + silicate_star + phosphate	2804	2.4397971	12971.2660	14.134608	1982-1999 –> 2000-2012	4.7984006	29987.9633
Atlantic	(27.5,27.75]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	2804	2.3858262	12847.8183	14.889428	1982-1999 –> 2000-2012	4.4857148	28998.8431
Atlantic	(27.5,27.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	2804	3.1165165	14346.1100	16.244315	1982-1999 –> 2000-2012	5.5026775	31451.5564
Atlantic	(27.5,27.75]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	2804	3.1165165	14346.1100	16.244315	1982-1999 –> 2000-2012	5.5026775	31451.5564
Atlantic	(27.5,27.75]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	1515	3.2267256	7862.9314	14.600864	2000-2012 –> 2013-2019	6.0506408	21656.1398
Atlantic	(27.5,27.75]	cstar_tref ~ sal + aou + nitrate + silicate_star	1515	3.2789491	7909.5783	16.534959	2000-2012 –> 2013-2019	6.1058258	21706.6649
Atlantic	(27.5,27.75]	cstar_tref ~ sal + aou + nitrate + silicate_star + phosphate_star	1515	3.2267256	7862.9314	14.600864	2000-2012 –> 2013-2019	6.0506408	21656.1398
Atlantic	(27.5,27.75]	cstar_tref ~ sal + aou + silicate + phosphate	1515	2.8745087	7510.7055	15.545879	2000-2012 –> 2013-2019	5.3284002	20514.2749
Atlantic	(27.5,27.75]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1515	2.8711798	7509.1946	15.461950	2000-2012 –> 2013-2019	5.2655862	20377.1448
Atlantic	(27.5,27.75]	cstar_tref ~ sal + aou + silicate_star + phosphate	1515	2.8460603	7480.5689	15.954890	2000-2012 –> 2013-2019	5.2858574	20451.8349
Atlantic	(27.5,27.75]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	1515	2.8440453	7480.4229	15.889857	2000-2012 –> 2013-2019	5.2298715	20328.2412
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 + nitrate + silicate_star + phosphate_star	1515	2.8464806	7483.0164	22.421750	2000-2012 –> 2013-2019	5.7223980	21378.5562
Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + nitrate + silicate	1036	1.9259602	4310.0805	13.267819	1982-1999 –> 2000-2012	3.9306085	9904.1154
Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	1036	1.8800580	4262.0997	14.057027	1982-1999 –> 2000-2012	3.8317057	9787.3967
Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + nitrate + silicate_star	1036	1.9259602	4310.0805	13.267819	1982-1999 –> 2000-2012	3.9306085	9904.1154
Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + nitrate + silicate_star + phosphate_star	1036	1.8800580	4262.0997	14.057027	1982-1999 –> 2000-2012	3.8317057	9787.3967
Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + phosphate + phosphate_star	1036	1.9943230	4382.3519	12.367472	1982-1999 –> 2000-2012	3.6776199	9515.1410
Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + silicate + phosphate	1036	1.7389083	4098.3902	12.226980	1982-1999 –> 2000-2012	3.5475856	9420.8408
Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1036	1.5641536	3880.9392	12.798383	1982-1999 –> 2000-2012	3.1611985	8876.8661
Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + silicate_star + phosphate	1036	1.7292672	4086.8704	12.169194	1982-1999 –> 2000-2012	3.5327060	9401.6636
Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	1036	1.5578540	3872.5773	12.740813	1982-1999 –> 2000-2012	3.1520469	8863.7855
Atlantic	(27.75,27.85]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1036	2.2923345	4672.9112	15.392219	1982-1999 –> 2000-2012	4.1974560	10134.5099
Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + nitrate + silicate	513	2.0314562	2195.0114	12.644220	2000-2012 –> 2013-2019	3.9574164	6505.0919
Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	513	2.0209724	2191.7028	12.956555	2000-2012 –> 2013-2019	3.9010304	6453.8025
Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + nitrate + silicate_star	513	2.0314562	2195.0114	12.644220	2000-2012 –> 2013-2019	3.9574164	6505.0919
Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + nitrate + silicate_star + phosphate_star	513	2.0209724	2191.7028	12.956555	2000-2012 –> 2013-2019	3.9010304	6453.8025
Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + silicate + phosphate	513	1.8570079	2102.8906	11.721525	2000-2012 –> 2013-2019	3.5959162	6201.2808
Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	513	1.7475521	2042.5606	12.039326	2000-2012 –> 2013-2019	3.3117057	5923.4997
Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + silicate_star + phosphate	513	1.8449695	2096.2177	11.652747	2000-2012 –> 2013-2019	3.5742368	6183.0881
Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	513	1.7382180	2037.0658	11.972402	2000-2012 –> 2013-2019	3.2960719	5909.6431
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 + nitrate + silicate_star + phosphate_star	513	2.0908491	2226.5780	12.385110	2000-2012 –> 2013-2019	4.1581813	6685.4273
Atlantic	(27.85,27.95]	cstar_tref ~ aou + phosphate + phosphate_star	983	3.4444208	5231.0950	40.406696	1982-1999 –> 2000-2012	6.3613830	11892.9889
Atlantic	(27.85,27.95]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	983	3.4332878	5226.7302	38.381356	1982-1999 –> 2000-2012	6.3108859	11854.3203
Atlantic	(27.85,27.95]	cstar_tref ~ aou + silicate_star + phosphate + phosphate_star	983	3.4335770	5226.8958	38.391899	1982-1999 –> 2000-2012	6.3120415	11855.2903
Atlantic	(27.85,27.95]	cstar_tref ~ sal + aou + phosphate + phosphate_star	983	3.4396213	5230.3536	36.781212	1982-1999 –> 2000-2012	6.1696286	11717.2584
Atlantic	(27.85,27.95]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	983	3.4272905	5225.2930	34.214054	1982-1999 –> 2000-2012	6.1571686	11714.0714
Atlantic	(27.85,27.95]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	983	3.4275839	5225.4613	34.223662	1982-1999 –> 2000-2012	6.1575457	11714.3215
Atlantic	(27.85,27.95]	cstar_tref ~ temp + aou + phosphate + phosphate_star	983	3.0946626	5022.5819	23.332862	1982-1999 –> 2000-2012	6.0114682	11686.3324
Atlantic	(27.85,27.95]	cstar_tref ~ temp + aou + silicate + phosphate	983	3.6929852	5370.0846	46.749916	1982-1999 –> 2000-2012	6.6056230	12030.0144
Atlantic	(27.85,27.95]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	983	3.0757549	5012.5332	23.130965	1982-1999 –> 2000-2012	5.9400103	11629.7050
Atlantic	(27.85,27.95]	cstar_tref ~ temp + aou + silicate_star + phosphate + phosphate_star	983	3.0758688	5012.6060	23.140031	1982-1999 –> 2000-2012	5.9407219	11630.3353
Atlantic	(27.85,27.95]	cstar_tref ~ aou + phosphate + phosphate_star	569	3.7400627	3125.8905	23.720140	2000-2012 –> 2013-2019	7.1844835	8356.9855
Atlantic	(27.85,27.95]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	569	3.7372699	3127.0405	23.268863	2000-2012 –> 2013-2019	7.1705577	8353.7707
Atlantic	(27.85,27.95]	cstar_tref ~ aou + silicate_star + phosphate + phosphate_star	569	3.7368389	3126.9092	23.232522	2000-2012 –> 2013-2019	7.1704159	8353.8050
Atlantic	(27.85,27.95]	cstar_tref ~ sal + aou + nitrate + silicate_star + phosphate_star	569	3.6673760	3107.5562	28.923043	2000-2012 –> 2013-2019	7.3896551	8495.1743
Atlantic	(27.85,27.95]	cstar_tref ~ sal + aou + phosphate + phosphate_star	569	3.4103854	3022.8793	30.447244	2000-2012 –> 2013-2019	6.8500067	8253.2329
Atlantic	(27.85,27.95]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	569	3.3215129	2994.8304	29.207068	2000-2012 –> 2013-2019	6.7488034	8220.1234
Atlantic	(27.85,27.95]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	569	3.3174363	2993.4329	29.173144	2000-2012 –> 2013-2019	6.7450201	8218.8941
Atlantic	(27.85,27.95]	cstar_tref ~ temp + aou + phosphate + phosphate_star	569	3.6730463	3107.3144	22.693963	2000-2012 –> 2013-2019	6.7677089	8129.8962
Atlantic	(27.85,27.95]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	569	3.6600265	3105.2733	23.151675	2000-2012 –> 2013-2019	6.7357814	8117.8065
Atlantic	(27.85,27.95]	cstar_tref ~ temp + aou + silicate_star + phosphate + phosphate_star	569	3.6603862	3105.3852	23.147178	2000-2012 –> 2013-2019	6.7362550	8117.9912
Atlantic	(27.95,28.05]	cstar_tref ~ oxygen + nitrate + silicate_star + 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 ~ sal + oxygen + nitrate + silicate_star + 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 + aou + nitrate + phosphate_star	1204	4.0815088	6815.5758	22.528945	1982-1999 –> 2000-2012	8.1122385	16256.3696
Atlantic	(27.95,28.05]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	1204	3.8590720	6682.6316	25.099675	1982-1999 –> 2000-2012	7.3908813	15682.5017
Atlantic	(27.95,28.05]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	1204	3.8590720	6682.6316	25.099675	1982-1999 –> 2000-2012	7.3908813	15682.5017
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 + nitrate + silicate_star + phosphate_star	1204	3.9591041	6744.2548	25.125796	1982-1999 –> 2000-2012	7.5958976	15842.3119
Atlantic	(27.95,28.05]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	666	4.2650178	3836.0207	24.914131	2000-2012 –> 2013-2019	8.7958040	10905.0610
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 ~ sal + oxygen + nitrate + silicate_star + 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 + aou + nitrate + phosphate_star	666	3.8465745	3696.4739	22.933854	2000-2012 –> 2013-2019	7.9280833	10512.0498
Atlantic	(27.95,28.05]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	666	3.5060655	3575.0128	26.986021	2000-2012 –> 2013-2019	7.3651375	10257.6443
Atlantic	(27.95,28.05]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	666	3.5060655	3575.0128	26.986021	2000-2012 –> 2013-2019	7.3651375	10257.6443
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 + nitrate + silicate_star + phosphate_star	666	3.5394291	3587.6281	27.891571	2000-2012 –> 2013-2019	7.4985332	10331.8828
Atlantic	(28.05,28.1]	cstar_tref ~ aou + phosphate + phosphate_star	804	1.1608853	2531.5393	11.005507	1982-1999 –> 2000-2012	2.1206840	5506.7947
Atlantic	(28.05,28.1]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	804	1.1604167	2532.8901	10.979649	1982-1999 –> 2000-2012	2.1134849	5495.0017
Atlantic	(28.05,28.1]	cstar_tref ~ aou + silicate_star + phosphate + phosphate_star	804	1.1604195	2532.8940	10.981201	1982-1999 –> 2000-2012	2.1135700	5495.1914
Atlantic	(28.05,28.1]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	804	1.1335334	2497.1994	6.397009	1982-1999 –> 2000-2012	2.1255963	5547.6063
Atlantic	(28.05,28.1]	cstar_tref ~ sal + aou + phosphate + phosphate_star	804	1.0598757	2387.1610	7.838670	1982-1999 –> 2000-2012	1.9527650	5208.9110
Atlantic	(28.05,28.1]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	804	1.0324411	2346.9902	6.945256	1982-1999 –> 2000-2012	1.9247967	5169.4535
Atlantic	(28.05,28.1]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	804	1.0319223	2346.1819	6.940660	1982-1999 –> 2000-2012	1.9242150	5168.4935
Atlantic	(28.05,28.1]	cstar_tref ~ temp + aou + phosphate + phosphate_star	804	1.1065403	2456.4443	11.331523	1982-1999 –> 2000-2012	2.0473492	5390.6954
Atlantic	(28.05,28.1]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	804	1.0436142	2364.2985	12.188672	1982-1999 –> 2000-2012	1.8828911	5054.7923
Atlantic	(28.05,28.1]	cstar_tref ~ temp + aou + silicate_star + phosphate + phosphate_star	804	1.0427660	2362.9910	12.172235	1982-1999 –> 2000-2012	1.8811182	5051.1125
Atlantic	(28.05,28.1]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	434	1.2132881	1411.4526	12.215619	2000-2012 –> 2013-2019	2.3737048	3944.3427
Atlantic	(28.05,28.1]	cstar_tref ~ aou + silicate_star + phosphate + phosphate_star	434	1.2132804	1411.4471	12.220901	2000-2012 –> 2013-2019	2.3736999	3944.3411
Atlantic	(28.05,28.1]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	434	1.1676572	1380.1780	6.819178	2000-2012 –> 2013-2019	2.3011906	3877.3774
Atlantic	(28.05,28.1]	cstar_tref ~ sal + aou + nitrate + silicate_star + phosphate_star	434	1.1676572	1380.1780	6.819178	2000-2012 –> 2013-2019	2.3011906	3877.3774
Atlantic	(28.05,28.1]	cstar_tref ~ sal + aou + phosphate + phosphate_star	434	1.1212849	1343.0032	8.929761	2000-2012 –> 2013-2019	2.1811606	3730.1642
Atlantic	(28.05,28.1]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	434	1.0586072	1295.0748	7.429722	2000-2012 –> 2013-2019	2.0910483	3642.0650
Atlantic	(28.05,28.1]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	434	1.0577152	1294.3431	7.425265	2000-2012 –> 2013-2019	2.0896374	3640.5250
Atlantic	(28.05,28.1]	cstar_tref ~ temp + aou + phosphate + phosphate_star	434	1.1480057	1363.4454	12.833810	2000-2012 –> 2013-2019	2.2545460	3819.8898
Atlantic	(28.05,28.1]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	434	1.1075442	1334.3006	13.578379	2000-2012 –> 2013-2019	2.1511584	3698.5991
Atlantic	(28.05,28.1]	cstar_tref ~ temp + aou + silicate_star + phosphate + phosphate_star	434	1.1066639	1333.6104	13.570865	2000-2012 –> 2013-2019	2.1494299	3696.6015
Atlantic	(28.1,28.15]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	936	0.8652007	2397.1990	9.973159	1982-1999 –> 2000-2012	1.5501093	5104.0130
Atlantic	(28.1,28.15]	cstar_tref ~ aou + silicate_star + phosphate + phosphate_star	936	0.8653236	2397.4650	9.979480	1982-1999 –> 2000-2012	1.5502922	5104.5055
Atlantic	(28.1,28.15]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	936	0.8394339	2342.6013	6.649159	1982-1999 –> 2000-2012	1.5477300	5138.3794
Atlantic	(28.1,28.15]	cstar_tref ~ sal + aou + nitrate + silicate_star + phosphate_star	936	0.8394339	2342.6013	6.649159	1982-1999 –> 2000-2012	1.5477300	5138.3794
Atlantic	(28.1,28.15]	cstar_tref ~ sal + aou + phosphate + phosphate_star	936	0.7732900	2186.9596	6.039973	1982-1999 –> 2000-2012	1.4451653	4844.0128
Atlantic	(28.1,28.15]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	936	0.7598609	2156.1646	5.927820	1982-1999 –> 2000-2012	1.3893138	4646.2932
Atlantic	(28.1,28.15]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	936	0.7598340	2156.0982	5.933231	1982-1999 –> 2000-2012	1.3890162	4645.1130
Atlantic	(28.1,28.15]	cstar_tref ~ temp + aou + phosphate + phosphate_star	936	0.8570872	2379.5613	9.929510	1982-1999 –> 2000-2012	1.5754781	5209.9919
Atlantic	(28.1,28.15]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	936	0.7856962	2218.7544	10.025704	1982-1999 –> 2000-2012	1.3744859	4535.9182
Atlantic	(28.1,28.15]	cstar_tref ~ temp + aou + silicate_star + phosphate + phosphate_star	936	0.7848370	2216.7061	10.063082	1982-1999 –> 2000-2012	1.3705842	4520.4516
Atlantic	(28.1,28.15]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	497	1.0349299	1456.5526	11.397933	2000-2012 –> 2013-2019	1.9001307	3853.7516
Atlantic	(28.1,28.15]	cstar_tref ~ aou + silicate_star + phosphate + phosphate_star	497	1.0349664	1456.5876	11.398469	2000-2012 –> 2013-2019	1.9002900	3854.0526
Atlantic	(28.1,28.15]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	497	0.9633503	1387.3108	6.933462	2000-2012 –> 2013-2019	1.8027842	3729.9121
Atlantic	(28.1,28.15]	cstar_tref ~ sal + aou + nitrate + silicate_star + phosphate_star	497	0.9633503	1387.3108	6.933462	2000-2012 –> 2013-2019	1.8027842	3729.9121
Atlantic	(28.1,28.15]	cstar_tref ~ sal + aou + phosphate + phosphate_star	497	0.8714406	1285.6430	6.181914	2000-2012 –> 2013-2019	1.6447306	3472.6026
Atlantic	(28.1,28.15]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	497	0.8712585	1287.4353	6.182326	2000-2012 –> 2013-2019	1.6311195	3443.5998
Atlantic	(28.1,28.15]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	497	0.8712841	1287.4645	6.182969	2000-2012 –> 2013-2019	1.6311181	3443.5627
Atlantic	(28.1,28.15]	cstar_tref ~ temp + aou + phosphate + phosphate_star	497	0.9653756	1387.3983	11.632299	2000-2012 –> 2013-2019	1.8224628	3766.9596
Atlantic	(28.1,28.15]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	497	0.8791505	1296.3985	11.933687	2000-2012 –> 2013-2019	1.6648467	3515.1529
Atlantic	(28.1,28.15]	cstar_tref ~ temp + aou + silicate_star + phosphate + phosphate_star	497	0.8784638	1295.6218	11.978025	2000-2012 –> 2013-2019	1.6633008	3512.3279
Atlantic	(28.15,28.2]	cstar_tref ~ sal + aou + nitrate + phosphate_star	1540	0.5939756	2777.9061	2.353787	1982-1999 –> 2000-2012	1.1267120	6300.3123
Atlantic	(28.15,28.2]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	1540	0.5471615	2527.0561	2.040783	1982-1999 –> 2000-2012	1.0648603	5924.1022
Atlantic	(28.15,28.2]	cstar_tref ~ sal + aou + nitrate + silicate_star + phosphate_star	1540	0.5471615	2527.0561	2.040783	1982-1999 –> 2000-2012	1.0648603	5924.1022
Atlantic	(28.15,28.2]	cstar_tref ~ temp + aou + nitrate + silicate	1540	0.5718023	2660.7281	2.192349	1982-1999 –> 2000-2012	1.1224532	6330.2485
Atlantic	(28.15,28.2]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	1540	0.5576121	2585.3284	2.314250	1982-1999 –> 2000-2012	1.1076613	6251.9858
Atlantic	(28.15,28.2]	cstar_tref ~ temp + aou + nitrate + silicate_star	1540	0.5718023	2660.7281	2.192349	1982-1999 –> 2000-2012	1.1224532	6330.2485
Atlantic	(28.15,28.2]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	1540	0.5576121	2585.3284	2.314250	1982-1999 –> 2000-2012	1.1076613	6251.9858
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 ~ temp + oxygen + nitrate + silicate_star + 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 + aou + nitrate + phosphate_star	819	0.7431790	1850.0329	4.262643	2000-2012 –> 2013-2019	1.3371546	4627.9389
Atlantic	(28.15,28.2]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	819	0.7372762	1838.9708	4.731734	2000-2012 –> 2013-2019	1.2844377	4366.0269
Atlantic	(28.15,28.2]	cstar_tref ~ sal + aou + nitrate + silicate_star + phosphate_star	819	0.7372762	1838.9708	4.731734	2000-2012 –> 2013-2019	1.2844377	4366.0269
Atlantic	(28.15,28.2]	cstar_tref ~ temp + aou + nitrate + silicate	819	0.8619250	2092.8357	8.363302	2000-2012 –> 2013-2019	1.4337273	4753.5638
Atlantic	(28.15,28.2]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	819	0.8613203	2093.6861	8.208278	2000-2012 –> 2013-2019	1.4189324	4679.0145
Atlantic	(28.15,28.2]	cstar_tref ~ temp + aou + nitrate + silicate_star	819	0.8619250	2092.8357	8.363302	2000-2012 –> 2013-2019	1.4337273	4753.5638
Atlantic	(28.15,28.2]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	819	0.8613203	2093.6861	8.208278	2000-2012 –> 2013-2019	1.4189324	4679.0145
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.15,28.2]	cstar_tref ~ temp + oxygen + nitrate + silicate_star + phosphate_star	819	0.8677435	2105.8561	8.181408	2000-2012 –> 2013-2019	1.4346601	4742.1542
Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate	2693	0.3593168	2139.5581	1.966552	1982-1999 –> 2000-2012	0.6655852	3898.5165
Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate + phosphate_star	2693	0.3540505	2062.0339	1.937681	1982-1999 –> 2000-2012	0.6560011	3717.6111
Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate + silicate	2693	0.3588037	2133.8610	1.963165	1982-1999 –> 2000-2012	0.6635911	3858.8422
Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	2693	0.3463114	1944.9979	1.900529	1982-1999 –> 2000-2012	0.6453988	3531.8596
Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate + silicate_star	2693	0.3588037	2133.8610	1.963165	1982-1999 –> 2000-2012	0.6635911	3858.8422
Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate + silicate_star + phosphate_star	2693	0.3463114	1944.9979	1.900529	1982-1999 –> 2000-2012	0.6453988	3531.8596
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 + nitrate + silicate_star + phosphate_star	2693	0.4618076	3495.1416	7.084808	1982-1999 –> 2000-2012	0.9515248	8741.7369
Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate	1557	0.4186021	1716.7958	2.479100	2000-2012 –> 2013-2019	0.7779188	3856.3538
Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate + phosphate_star	1557	0.3991609	1570.7062	2.411345	2000-2012 –> 2013-2019	0.7532114	3632.7401
Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate + silicate	1557	0.4110188	1661.8663	2.465092	2000-2012 –> 2013-2019	0.7698225	3795.7273
Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	1557	0.3901368	1501.4981	2.457959	2000-2012 –> 2013-2019	0.7364483	3446.4960
Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate + silicate_star	1557	0.4110188	1661.8663	2.465092	2000-2012 –> 2013-2019	0.7698225	3795.7273
Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate + silicate_star + phosphate_star	1557	0.3901368	1501.4981	2.457959	2000-2012 –> 2013-2019	0.7364483	3446.4960
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 + nitrate + silicate_star + phosphate_star	1557	0.4385296	1865.6176	5.067855	2000-2012 –> 2013-2019	0.9003372	5360.7592
Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + aou + phosphate + phosphate_star	5288	6.2360939	34376.5179	29.968567	1982-1999 –> 2000-2012	12.7656352	81794.3473
Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	5288	6.1928478	34304.9200	29.527354	1982-1999 –> 2000-2012	12.6716991	81612.6324
Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	5288	6.2330691	34373.3868	30.214602	1982-1999 –> 2000-2012	12.7604122	81788.3722
Indo-Pacific	(-Inf,26]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	5288	7.1987438	35896.7294	31.371808	1982-1999 –> 2000-2012	14.7271008	85364.2640
Indo-Pacific	(-Inf,26]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	5288	7.1987438	35896.7294	31.371808	1982-1999 –> 2000-2012	14.7271008	85364.2640
Indo-Pacific	(-Inf,26]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	5288	7.2933052	36034.7494	29.166512	1982-1999 –> 2000-2012	14.8774849	85608.5624
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 + nitrate + silicate_star + 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 + 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 ~ sal + aou + phosphate + phosphate_star	2816	5.9520247	18049.4373	24.894024	2000-2012 –> 2013-2019	12.1881186	52425.9552
Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2816	5.8778725	17980.8313	24.128023	2000-2012 –> 2013-2019	12.0707203	52285.7513
Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	2816	5.9519636	18051.3795	24.924420	2000-2012 –> 2013-2019	12.1850327	52424.7664
Indo-Pacific	(-Inf,26]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	2816	6.9521542	18926.2002	25.907885	2000-2012 –> 2013-2019	14.1508980	54822.9297
Indo-Pacific	(-Inf,26]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	2816	6.9521542	18926.2002	25.907885	2000-2012 –> 2013-2019	14.1508980	54822.9297
Indo-Pacific	(-Inf,26]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	2816	7.0001768	18964.9700	25.465195	2000-2012 –> 2013-2019	14.2934820	54999.7194
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 + nitrate + silicate_star + 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 + 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 + silicate + phosphate + phosphate_star	2816	7.1344030	19069.9395	24.441528	2000-2012 –> 2013-2019	14.5251760	55243.0906
Indo-Pacific	(26,26.5]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	5525	4.8905258	33230.9339	41.688307	1982-1999 –> 2000-2012	9.5377100	78823.9388
Indo-Pacific	(26,26.5]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	5525	4.5360650	32401.5327	43.772780	1982-1999 –> 2000-2012	8.7824377	76605.3494
Indo-Pacific	(26,26.5]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	5525	4.9090739	33274.7635	40.067718	1982-1999 –> 2000-2012	9.5566695	78871.1338
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	5525	4.7555301	32923.6260	49.041150	1982-1999 –> 2000-2012	9.1671720	77716.3587
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	5525	4.7555301	32923.6260	49.041150	1982-1999 –> 2000-2012	9.1671720	77716.3587
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + aou + silicate + phosphate	5525	4.8893893	33228.3657	40.409229	1982-1999 –> 2000-2012	9.5091497	78730.0815
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	5525	4.8803203	33209.8508	41.154556	1982-1999 –> 2000-2012	9.4999894	78713.2617
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + aou + silicate + phosphate_star	5525	4.9245814	33307.6149	49.529819	1982-1999 –> 2000-2012	9.5559722	78848.1124
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 + nitrate + silicate_star + phosphate_star	5525	4.7637711	32942.7584	51.203206	1982-1999 –> 2000-2012	9.1789275	77747.7736
Indo-Pacific	(26,26.5]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	2918	4.8175903	17468.7156	40.382988	2000-2012 –> 2013-2019	9.7081161	50699.6495
Indo-Pacific	(26,26.5]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2918	4.4724496	17036.8813	46.310491	2000-2012 –> 2013-2019	9.0085146	49438.4140
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	2918	4.5715543	17164.7889	45.715284	2000-2012 –> 2013-2019	9.3270844	50088.4149
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	2918	4.5715543	17164.7889	45.715284	2000-2012 –> 2013-2019	9.3270844	50088.4149
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + aou + silicate + phosphate	2918	4.8142859	17464.7112	37.043514	2000-2012 –> 2013-2019	9.7036752	50693.0769
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	2918	4.8065118	17457.2797	38.739984	2000-2012 –> 2013-2019	9.6868321	50667.1305
Indo-Pacific	(26,26.5]	cstar_tref ~ temp + aou + silicate + phosphate_star	2918	4.8483037	17505.8034	34.431113	2000-2012 –> 2013-2019	9.7728851	50813.4183
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 + nitrate + silicate_star + 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 + silicate + phosphate_star	2918	4.8517131	17509.9060	34.225344	2000-2012 –> 2013-2019	9.7822025	50830.7695
Indo-Pacific	(26.5,26.75]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	4544	4.0473152	25612.8259	22.127132	1982-1999 –> 2000-2012	7.9963826	60742.2745
Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	4544	3.9265162	25339.4488	20.226531	1982-1999 –> 2000-2012	7.6839109	59845.1956
Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	4544	4.1898965	25929.4731	16.421495	1982-1999 –> 2000-2012	8.0951171	60920.5086
Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	4544	4.0218289	25557.4171	33.030413	1982-1999 –> 2000-2012	8.0127086	60821.3190
Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	4544	4.0218289	25557.4171	33.030413	1982-1999 –> 2000-2012	8.0127086	60821.3190
Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	4544	3.9477481	25388.4579	23.151958	1982-1999 –> 2000-2012	7.8249401	60288.9073
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 + nitrate + silicate_star + 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 + 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 + nitrate + silicate_star + 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 ~ aou + silicate + phosphate + phosphate_star	2447	4.4033896	14211.0265	24.538450	2000-2012 –> 2013-2019	8.4507048	39823.8524
Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2447	4.2938981	14089.7973	24.434470	2000-2012 –> 2013-2019	8.2204143	39429.2462
Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	2447	4.3225383	14122.3318	31.664753	2000-2012 –> 2013-2019	8.3443672	39679.7489
Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	2447	4.3225383	14122.3318	31.664753	2000-2012 –> 2013-2019	8.3443672	39679.7489
Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	2447	4.2923135	14087.9909	22.815590	2000-2012 –> 2013-2019	8.2400615	39476.4489
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 + nitrate + silicate_star + 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 + nitrate + silicate_star + phosphate_star	2447	4.3205048	14120.0290	31.479242	2000-2012 –> 2013-2019	8.3423820	39677.5550
Indo-Pacific	(26.75,27]	cstar_tref ~ aou + phosphate + phosphate_star	5354	4.5403612	31405.2681	19.876020	1982-1999 –> 2000-2012	8.5493390	73701.4539
Indo-Pacific	(26.75,27]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	5354	4.5232238	31366.7748	18.928405	1982-1999 –> 2000-2012	8.5202080	73619.8323
Indo-Pacific	(26.75,27]	cstar_tref ~ aou + silicate_star + phosphate + phosphate_star	5354	4.4364272	31159.3008	18.611622	1982-1999 –> 2000-2012	8.3694630	73169.4305
Indo-Pacific	(26.75,27]	cstar_tref ~ sal + aou + phosphate + phosphate_star	5354	4.3077986	30844.2463	21.980567	1982-1999 –> 2000-2012	7.9956057	71884.6902
Indo-Pacific	(26.75,27]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	5354	4.2906105	30803.4359	21.052395	1982-1999 –> 2000-2012	7.9645082	71788.9629
Indo-Pacific	(26.75,27]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	5354	4.2213053	30629.0600	20.756473	1982-1999 –> 2000-2012	7.8498644	71427.5545
Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + phosphate + phosphate_star	5354	4.2057274	30587.4712	23.104196	1982-1999 –> 2000-2012	7.8859232	71596.7963
Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	5354	4.1702723	30498.8181	21.453131	1982-1999 –> 2000-2012	7.8227361	71396.2141
Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + silicate_star + phosphate	5354	4.6086416	31567.1021	20.206504	1982-1999 –> 2000-2012	8.7190287	74241.5503
Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + silicate_star + phosphate + phosphate_star	5354	4.1165993	30360.1074	20.285252	1982-1999 –> 2000-2012	7.7366826	71123.3780
Indo-Pacific	(26.75,27]	cstar_tref ~ aou + phosphate + phosphate_star	2868	4.9134802	17280.6429	19.343238	2000-2012 –> 2013-2019	9.4538414	48685.9111
Indo-Pacific	(26.75,27]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	2868	4.9018613	17269.0630	18.494849	2000-2012 –> 2013-2019	9.4250851	48635.8378
Indo-Pacific	(26.75,27]	cstar_tref ~ aou + silicate_star + phosphate + phosphate_star	2868	4.8115172	17162.3589	18.671411	2000-2012 –> 2013-2019	9.2479444	48321.6597
Indo-Pacific	(26.75,27]	cstar_tref ~ sal + aou + phosphate + phosphate_star	2868	4.6612052	16980.3073	20.095371	2000-2012 –> 2013-2019	8.9690038	47824.5536
Indo-Pacific	(26.75,27]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2868	4.6499950	16968.4956	19.282213	2000-2012 –> 2013-2019	8.9406055	47771.9315
Indo-Pacific	(26.75,27]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	2868	4.5834469	16885.8122	19.621984	2000-2012 –> 2013-2019	8.8047522	47514.8721
Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + phosphate + phosphate_star	2868	4.5774010	16876.2410	16.522089	2000-2012 –> 2013-2019	8.7831284	47463.7121
Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	2868	4.5537010	16848.4650	17.066955	2000-2012 –> 2013-2019	8.7239733	47347.2831
Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + silicate_star + phosphate	2868	5.0050395	17388.5457	19.071061	2000-2012 –> 2013-2019	9.6136811	48955.6478
Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + silicate_star + phosphate + phosphate_star	2868	4.5010869	16781.8046	18.240062	2000-2012 –> 2013-2019	8.6176862	47141.9121
Indo-Pacific	(27,27.25]	cstar_tref ~ aou + silicate_star + phosphate + phosphate_star	6466	4.0461089	36437.4893	40.815494	1982-1999 –> 2000-2012	7.3572918	82688.9862
Indo-Pacific	(27,27.25]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	6466	4.0166063	36344.8487	42.576745	1982-1999 –> 2000-2012	7.2664602	82267.9222
Indo-Pacific	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	6466	4.1535561	36778.4265	49.359080	1982-1999 –> 2000-2012	7.5202091	83325.5519
Indo-Pacific	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	6466	4.1535561	36778.4265	49.359080	1982-1999 –> 2000-2012	7.5202091	83325.5519
Indo-Pacific	(27,27.25]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	6466	3.9500266	36128.6900	31.512269	1982-1999 –> 2000-2012	7.2528573	82337.5496
Indo-Pacific	(27,27.25]	cstar_tref ~ temp + aou + silicate_star + phosphate	6466	4.1697945	36826.8859	39.883706	1982-1999 –> 2000-2012	7.5915871	83659.1340
Indo-Pacific	(27,27.25]	cstar_tref ~ temp + aou + silicate_star + phosphate + phosphate_star	6466	3.7382856	35416.1990	36.297381	1982-1999 –> 2000-2012	6.8506238	80575.1233
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 + nitrate + silicate_star + 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 + silicate_star + phosphate	6466	4.2699704	37133.8935	39.103729	1982-1999 –> 2000-2012	7.7978077	84505.5602
Indo-Pacific	(27,27.25]	cstar_tref ~ aou + silicate_star + phosphate + phosphate_star	3470	4.7690854	20700.7860	33.838511	2000-2012 –> 2013-2019	8.8151944	57138.2754
Indo-Pacific	(27,27.25]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	3470	4.7342458	20651.9011	35.629548	2000-2012 –> 2013-2019	8.7508521	56996.7498
Indo-Pacific	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	3470	4.9716940	20991.5322	35.708885	2000-2012 –> 2013-2019	9.1252501	57769.9588
Indo-Pacific	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	3470	4.9716940	20991.5322	35.708885	2000-2012 –> 2013-2019	9.1252501	57769.9588
Indo-Pacific	(27,27.25]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	3470	4.4723380	20256.9380	25.616236	2000-2012 –> 2013-2019	8.4223646	56385.6280
Indo-Pacific	(27,27.25]	cstar_tref ~ temp + aou + silicate_star + phosphate	3470	4.8755660	20854.0327	32.501494	2000-2012 –> 2013-2019	9.0453605	57680.9186
Indo-Pacific	(27,27.25]	cstar_tref ~ temp + aou + silicate_star + phosphate + phosphate_star	3470	4.3257326	20025.6292	30.341767	2000-2012 –> 2013-2019	8.0640183	55441.8282
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 + nitrate + silicate_star + 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_star + phosphate_star	3470	4.9514942	20961.2778	31.586747	2000-2012 –> 2013-2019	9.2214646	58095.1713
Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	5847	3.4261572	31007.5182	32.305003	1982-1999 –> 2000-2012	6.3475428	70339.4077
Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + aou + nitrate + silicate_star + phosphate_star	5847	3.4261572	31007.5182	32.305003	1982-1999 –> 2000-2012	6.3475428	70339.4077
Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	5847	3.4479927	31081.8095	32.787008	1982-1999 –> 2000-2012	6.3639079	70384.1153
Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + aou + silicate_star + phosphate_star	5847	3.4481958	31080.4983	32.737786	1982-1999 –> 2000-2012	6.3699577	70412.4213
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 + nitrate + silicate_star + 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_star + phosphate	5847	3.4677649	31146.6762	32.201955	1982-1999 –> 2000-2012	6.3842745	70450.1996
Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + silicate_star + phosphate + phosphate_star	5847	3.4677649	31146.6762	32.201955	1982-1999 –> 2000-2012	6.3842745	70450.1996
Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	5847	3.3461863	30731.3292	33.390589	1982-1999 –> 2000-2012	6.2704796	70078.9213
Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + aou + silicate_star + phosphate + phosphate_star	5847	3.2436714	30367.4652	32.825325	1982-1999 –> 2000-2012	6.1083715	69390.0778
Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	3127	3.9776245	17522.8442	31.036835	2000-2012 –> 2013-2019	7.4037817	48530.3624
Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + aou + nitrate + silicate_star + phosphate_star	3127	3.9776245	17522.8442	31.036835	2000-2012 –> 2013-2019	7.4037817	48530.3624
Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	3127	4.0149436	17581.2473	32.065458	2000-2012 –> 2013-2019	7.4629363	48663.0568
Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + aou + silicate_star + phosphate_star	3127	4.0191802	17585.8431	31.731513	2000-2012 –> 2013-2019	7.4673759	48666.3413
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 + nitrate + silicate_star + 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_star + phosphate	3127	4.0587225	17647.0717	30.845627	2000-2012 –> 2013-2019	7.5264874	48793.7479
Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + silicate_star + phosphate + phosphate_star	3127	4.0587225	17647.0717	30.845627	2000-2012 –> 2013-2019	7.5264874	48793.7479
Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	3127	3.8437243	17308.6883	33.232474	2000-2012 –> 2013-2019	7.1899105	48040.0175
Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + aou + silicate_star + phosphate + phosphate_star	3127	3.7456570	17147.0555	32.042107	2000-2012 –> 2013-2019	6.9893284	47514.5206
Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	5836	3.2683124	30398.6938	24.119856	1982-1999 –> 2000-2012	5.5822735	67013.2424
Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + nitrate + silicate_star + phosphate_star	5836	3.2683124	30398.6938	24.119856	1982-1999 –> 2000-2012	5.5822735	67013.2424
Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	5836	3.2280103	30253.8699	9.417658	1982-1999 –> 2000-2012	5.5715676	67074.4342
Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	5836	3.1845769	30095.7550	10.967937	1982-1999 –> 2000-2012	5.4967039	66697.4500
Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + silicate_star + phosphate_star	5836	3.3188975	30575.9628	22.448148	1982-1999 –> 2000-2012	5.6411441	67246.4502
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 + nitrate + silicate_star + 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_star + phosphate	5836	3.3745368	30770.0148	25.183989	1982-1999 –> 2000-2012	5.7083873	67521.2992
Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + silicate_star + phosphate + phosphate_star	5836	3.3745368	30770.0148	25.183989	1982-1999 –> 2000-2012	5.7083873	67521.2992
Indo-Pacific	(27.5,27.75]	cstar_tref ~ temp + aou + silicate_star + phosphate + phosphate_star	5836	3.3018702	30517.9267	14.570342	1982-1999 –> 2000-2012	5.6994111	67707.6509
Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	3104	3.9417327	17337.7899	25.386634	2000-2012 –> 2013-2019	7.2100451	47736.4837
Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + nitrate + silicate_star + phosphate_star	3104	3.9417327	17337.7899	25.386634	2000-2012 –> 2013-2019	7.2100451	47736.4837
Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	3104	3.7804110	17078.3721	11.020697	2000-2012 –> 2013-2019	7.0084214	47332.2420
Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + silicate + phosphate_star	3104	4.0989818	17578.6356	23.944569	2000-2012 –> 2013-2019	7.4797062	48370.0328
Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	3104	3.7296814	16994.5025	10.660379	2000-2012 –> 2013-2019	6.9142583	47090.2575
Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + silicate_star + phosphate_star	3104	4.0195344	17457.1293	23.663404	2000-2012 –> 2013-2019	7.3384319	48033.0922
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 + nitrate + silicate_star + 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 ~ temp + aou + silicate + phosphate + phosphate_star	3104	4.0127640	17448.6640	19.751773	2000-2012 –> 2013-2019	7.3819468	48202.1453
Indo-Pacific	(27.5,27.75]	cstar_tref ~ temp + aou + silicate_star + phosphate + phosphate_star	3104	3.9060342	17281.3106	18.885415	2000-2012 –> 2013-2019	7.2079044	47799.2372
Indo-Pacific	(27.75,27.85]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	2317	2.7304690	11242.0909	21.046900	1982-1999 –> 2000-2012	4.6623579	24262.9999
Indo-Pacific	(27.75,27.85]	cstar_tref ~ aou + silicate_star + phosphate + phosphate_star	2317	2.7260850	11234.6446	21.214199	1982-1999 –> 2000-2012	4.6543997	24243.9574
Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2317	2.7290365	11241.6591	21.325228	1982-1999 –> 2000-2012	4.6507603	24231.5321
Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	2317	2.7249081	11234.6436	21.471767	1982-1999 –> 2000-2012	4.6432603	24213.5203
Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	2317	2.7271269	11238.4154	21.453348	1982-1999 –> 2000-2012	4.6635446	24275.9866
Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	2317	2.7271269	11238.4154	21.453348	1982-1999 –> 2000-2012	4.6635446	24275.9866
Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	2317	2.6406178	11089.0349	20.968772	1982-1999 –> 2000-2012	4.5482001	24032.6569
Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + aou + silicate_star + phosphate + phosphate_star	2317	2.6387258	11085.7135	21.377698	1982-1999 –> 2000-2012	4.5432730	24019.3645
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 + nitrate + silicate_star + 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 ~ sal + aou + nitrate + silicate + phosphate_star	1264	3.1358059	6490.2930	31.989646	2000-2012 –> 2013-2019	5.8629052	17728.6613
Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + aou + nitrate + silicate_star + phosphate_star	1264	3.1358059	6490.2930	31.989646	2000-2012 –> 2013-2019	5.8629052	17728.6613
Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	1264	3.1665018	6514.9188	30.688074	2000-2012 –> 2013-2019	5.8914099	17749.5624
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 + nitrate + silicate_star + 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 + nitrate + silicate_star + 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 ~ temp + aou + phosphate + phosphate_star	1264	3.1128492	6469.7178	31.766441	2000-2012 –> 2013-2019	5.8536894	17729.3767
Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1264	3.0665742	6433.8549	28.861392	2000-2012 –> 2013-2019	5.7071919	17522.8898
Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + aou + silicate_star + phosphate + phosphate_star	1264	3.0580694	6426.8341	29.049935	2000-2012 –> 2013-2019	5.6967952	17512.5476
Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + aou + nitrate + phosphate_star	2907	2.7487099	14140.4282	28.492785	1982-1999 –> 2000-2012	4.8701028	31563.9400
Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	2907	2.7457460	14136.1557	28.208152	1982-1999 –> 2000-2012	4.8664196	31558.9474
Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	2907	2.7457460	14136.1557	28.208152	1982-1999 –> 2000-2012	4.8664196	31558.9474
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 + nitrate + silicate_star + 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 + 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 + nitrate + silicate_star + 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 + aou + nitrate + silicate	1542	3.3046795	8074.4335	31.424170	2000-2012 –> 2013-2019	6.0531124	22214.2756
Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	1542	3.2946985	8067.1049	31.554855	2000-2012 –> 2013-2019	6.0404446	22203.2606
Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + aou + nitrate + silicate_star	1542	3.3046795	8074.4335	31.424170	2000-2012 –> 2013-2019	6.0531124	22214.2756
Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	1542	3.2946985	8067.1049	31.554855	2000-2012 –> 2013-2019	6.0404446	22203.2606
Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1542	3.2637446	8037.9935	37.426960	2000-2012 –> 2013-2019	6.0494318	22258.1136
Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + aou + silicate_star + phosphate + phosphate_star	1542	3.2661935	8040.3067	37.431137	2000-2012 –> 2013-2019	6.0559877	22268.9922
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.85,27.95]	cstar_tref ~ temp + oxygen + nitrate + silicate_star	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_star + 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 ~ aou + silicate + phosphate + phosphate_star	2467	2.2624239	11041.3418	8.237248	1982-1999 –> 2000-2012	4.1845877	25738.2975
Indo-Pacific	(27.95,28.05]	cstar_tref ~ aou + silicate_star + phosphate + phosphate_star	2467	2.2603950	11036.9149	8.232769	1982-1999 –> 2000-2012	4.1832993	25736.6002
Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	2467	2.1936153	10890.9515	8.036190	1982-1999 –> 2000-2012	4.1649943	25769.0539
Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + aou + nitrate + silicate_star + phosphate_star	2467	2.1936153	10890.9515	8.036190	1982-1999 –> 2000-2012	4.1649943	25769.0539
Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	2467	2.2115502	10931.1275	8.468880	1982-1999 –> 2000-2012	4.1022689	25513.2034
Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	2467	2.2063914	10919.6046	8.494541	1982-1999 –> 2000-2012	4.0944790	25491.8129
Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	2467	2.2668883	11053.0683	8.062293	1982-1999 –> 2000-2012	4.2018862	25799.1794
Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	2467	2.2668883	11053.0683	8.062293	1982-1999 –> 2000-2012	4.2018862	25799.1794
Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	2467	2.1886299	10879.7253	7.878228	1982-1999 –> 2000-2012	4.0642895	25405.1370
Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + aou + silicate_star + phosphate + phosphate_star	2467	2.1895879	10881.8845	7.849347	1982-1999 –> 2000-2012	4.0701805	25425.9079
Indo-Pacific	(27.95,28.05]	cstar_tref ~ aou + silicate_star + phosphate + phosphate_star	1354	2.5367274	6375.2946	7.703273	2000-2012 –> 2013-2019	4.7971224	17412.2095
Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + aou + nitrate + phosphate_star	1354	2.4060015	6232.0181	8.104133	2000-2012 –> 2013-2019	4.6173063	17160.5981
Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	1354	2.4007452	6228.0956	8.259744	2000-2012 –> 2013-2019	4.5943605	17119.0470
Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + aou + nitrate + silicate_star + phosphate_star	1354	2.4007452	6228.0956	8.259744	2000-2012 –> 2013-2019	4.5943605	17119.0470
Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + aou + phosphate + phosphate_star	1354	2.4939482	6329.2376	8.352020	2000-2012 –> 2013-2019	4.7451070	17345.9505
Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1354	2.4779038	6313.7598	8.329809	2000-2012 –> 2013-2019	4.6894540	17244.8873
Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	1354	2.4735198	6308.9645	8.317297	2000-2012 –> 2013-2019	4.6799112	17228.5692
Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + aou + phosphate + phosphate_star	1354	2.4648617	6297.4690	8.721893	2000-2012 –> 2013-2019	4.7165003	17315.2334
Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1354	2.4266912	6257.2052	8.523746	2000-2012 –> 2013-2019	4.6153211	17136.9304
Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + aou + silicate_star + phosphate + phosphate_star	1354	2.4267117	6257.2281	8.510320	2000-2012 –> 2013-2019	4.6162997	17139.1126
Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + aou + silicate + phosphate	1718	1.9503440	7182.7407	8.816710	1982-1999 –> 2000-2012	3.5709134	16908.7065
Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	1718	1.8450332	6994.0133	9.467590	1982-1999 –> 2000-2012	3.4414783	16645.3685
Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + aou + silicate_star + phosphate	1718	1.9495859	7181.4047	8.849010	1982-1999 –> 2000-2012	3.5737884	16918.8092
Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	1718	1.8482228	6999.9483	9.397411	1982-1999 –> 2000-2012	3.4507256	16670.6487
Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	1718	1.8612616	7024.1034	10.244079	1982-1999 –> 2000-2012	3.4707315	16716.9636
Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	1718	1.8612616	7024.1034	10.244079	1982-1999 –> 2000-2012	3.4707315	16716.9636
Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	1718	1.8032691	6915.3424	9.697645	1982-1999 –> 2000-2012	3.4298781	16662.3098
Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + aou + silicate_star + phosphate + phosphate_star	1718	1.8210070	6948.9757	9.802522	1982-1999 –> 2000-2012	3.4729607	16774.9183
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 + nitrate + silicate_star + 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 ~ aou + silicate + phosphate + phosphate_star	966	2.1066397	4192.9110	9.987464	2000-2012 –> 2013-2019	3.9914816	11258.2717
Indo-Pacific	(28.05,28.1]	cstar_tref ~ aou + silicate_star + phosphate + phosphate_star	966	2.1150545	4200.6129	9.941901	2000-2012 –> 2013-2019	4.0093242	11283.1172
Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	966	2.0390556	4131.9137	9.747515	2000-2012 –> 2013-2019	3.8840888	11125.9270
Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	966	2.0395415	4132.3740	9.670797	2000-2012 –> 2013-2019	3.8877643	11132.3223
Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	966	2.0868098	4176.6389	10.548800	2000-2012 –> 2013-2019	3.9480714	11200.7423
Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	966	2.0868098	4176.6389	10.548800	2000-2012 –> 2013-2019	3.9480714	11200.7423
Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	966	2.0001288	4094.6740	9.111786	2000-2012 –> 2013-2019	3.8033978	11010.0164
Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + aou + silicate_star + phosphate + phosphate_star	966	2.0180718	4111.9285	9.104081	2000-2012 –> 2013-2019	3.8390788	11060.9042
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 + nitrate + silicate_star + phosphate_star	966	2.0835440	4173.6130	10.558303	2000-2012 –> 2013-2019	3.9420347	11192.5974
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 + aou + silicate + phosphate + phosphate_star	14627	1.3378772	50039.0036	16.225048	1982-1999 –> 2000-2012	2.5090135	114186.0580
Indo-Pacific	(28.1, Inf]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	14627	1.3532531	50373.2971	16.386271	1982-1999 –> 2000-2012	2.5431232	115165.7673
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 + nitrate + silicate_star + phosphate_star	14627	1.3133686	49498.1312	9.366747	1982-1999 –> 2000-2012	2.4634246	112906.4588
Indo-Pacific	(28.1, Inf]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	14627	1.3692488	50717.0570	13.004762	1982-1999 –> 2000-2012	2.4951048	113260.4541
Indo-Pacific	(28.1, Inf]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	14627	1.3692488	50717.0570	13.004762	1982-1999 –> 2000-2012	2.4951048	113260.4541
Indo-Pacific	(28.1, Inf]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	14627	1.3940533	51242.2624	17.034744	1982-1999 –> 2000-2012	2.6340769	117713.8450
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 + nitrate + silicate_star + phosphate_star	14627	1.3591254	50499.9672	12.907049	1982-1999 –> 2000-2012	2.4770940	112757.4410
Indo-Pacific	(28.1, Inf]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	7804	1.4923380	28409.3615	12.644005	2000-2012 –> 2013-2019	2.8302151	78448.3652
Indo-Pacific	(28.1, Inf]	cstar_tref ~ sal + aou + silicate_star + phosphate + phosphate_star	7804	1.5067868	28559.7515	12.779785	2000-2012 –> 2013-2019	2.8600399	78933.0486
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 + nitrate + silicate_star + phosphate_star	7804	1.4569278	28034.5504	8.140875	2000-2012 –> 2013-2019	2.7702964	77532.6816
Indo-Pacific	(28.1, Inf]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	7804	1.5858456	29357.9189	10.996188	2000-2012 –> 2013-2019	2.9550945	80074.9759
Indo-Pacific	(28.1, Inf]	cstar_tref ~ temp + aou + nitrate + silicate_star + phosphate_star	7804	1.5858456	29357.9189	10.996188	2000-2012 –> 2013-2019	2.9550945	80074.9759
Indo-Pacific	(28.1, Inf]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	7804	1.5521124	29022.3323	13.774789	2000-2012 –> 2013-2019	2.9461657	80264.5947
Indo-Pacific	(28.1, Inf]	cstar_tref ~ temp + aou + silicate_star + phosphate + phosphate_star	7804	1.5749885	29250.6951	13.940540	2000-2012 –> 2013-2019	2.9909908	80949.9639
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
Indo-Pacific	(28.1, Inf]	cstar_tref ~ temp + oxygen + nitrate + silicate_star + 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	2.9648446	1351.575
Atlantic	(-Inf,26]	2000-2012 –> 2013-2019	2.9008305	836.842
Atlantic	(26,26.5]	1982-1999 –> 2000-2012	7.9200963	14395.325
Atlantic	(26,26.5]	2000-2012 –> 2013-2019	7.7966189	9004.455
Atlantic	(26.5,26.75]	1982-1999 –> 2000-2012	6.4587571	17060.592
Atlantic	(26.5,26.75]	2000-2012 –> 2013-2019	6.7505930	10828.281
Atlantic	(26.75,27]	1982-1999 –> 2000-2012	3.9357835	23010.622
Atlantic	(26.75,27]	2000-2012 –> 2013-2019	4.7997317	16567.196
Atlantic	(27,27.25]	1982-1999 –> 2000-2012	3.8870311	19239.030
Atlantic	(27,27.25]	2000-2012 –> 2013-2019	4.4245801	13972.018
Atlantic	(27.25,27.5]	1982-1999 –> 2000-2012	4.3009500	21118.666
Atlantic	(27.25,27.5]	2000-2012 –> 2013-2019	5.3093507	14824.620
Atlantic	(27.5,27.75]	1982-1999 –> 2000-2012	5.0954003	30610.509
Atlantic	(27.5,27.75]	2000-2012 –> 2013-2019	5.6618494	21098.647
Atlantic	(27.75,27.85]	1982-1999 –> 2000-2012	3.6793241	9559.583
Atlantic	(27.75,27.85]	2000-2012 –> 2013-2019	3.7811187	6350.616
Atlantic	(27.85,27.95]	1982-1999 –> 2000-2012	6.1966477	11772.464
Atlantic	(27.85,27.95]	2000-2012 –> 2013-2019	6.9498688	8261.768
Atlantic	(27.95,28.05]	1982-1999 –> 2000-2012	7.9387786	16078.574
Atlantic	(27.95,28.05]	2000-2012 –> 2013-2019	7.8271950	10483.822
Atlantic	(28.05,28.1]	1982-1999 –> 2000-2012	2.0086471	5308.805
Atlantic	(28.05,28.1]	2000-2012 –> 2013-2019	2.2266766	3787.128
Atlantic	(28.1,28.15]	1982-1999 –> 2000-2012	1.4739905	4888.706
Atlantic	(28.1,28.15]	2000-2012 –> 2013-2019	1.7463568	3632.183
Atlantic	(28.15,28.2]	1982-1999 –> 2000-2012	1.1105014	6249.093
Atlantic	(28.15,28.2]	2000-2012 –> 2013-2019	1.3927654	4650.713
Atlantic	(28.2, Inf]	1982-1999 –> 2000-2012	0.7776275	5774.966
Atlantic	(28.2, Inf]	2000-2012 –> 2013-2019	0.8202827	4389.239
Indo-Pacific	(-Inf,26]	1982-1999 –> 2000-2012	14.2024731	84394.987
Indo-Pacific	(-Inf,26]	2000-2012 –> 2013-2019	13.6247032	54174.695
Indo-Pacific	(26,26.5]	1982-1999 –> 2000-2012	9.3134127	78152.014
Indo-Pacific	(26,26.5]	2000-2012 –> 2013-2019	9.4999911	50355.080
Indo-Pacific	(26.5,26.75]	1982-1999 –> 2000-2012	7.9827628	60718.862
Indo-Pacific	(26.5,26.75]	2000-2012 –> 2013-2019	8.3527906	39680.613
Indo-Pacific	(26.75,27]	1982-1999 –> 2000-2012	8.1413359	72394.986
Indo-Pacific	(26.75,27]	2000-2012 –> 2013-2019	9.0579701	47966.332
Indo-Pacific	(27,27.25]	1982-1999 –> 2000-2012	7.4118077	82903.951
Indo-Pacific	(27,27.25]	2000-2012 –> 2013-2019	8.8652421	57248.569
Indo-Pacific	(27.25,27.5]	1982-1999 –> 2000-2012	6.3247247	70240.842
Indo-Pacific	(27.25,27.5]	2000-2012 –> 2013-2019	7.3772416	48457.011
Indo-Pacific	(27.5,27.75]	1982-1999 –> 2000-2012	5.6226820	67212.120
Indo-Pacific	(27.5,27.75]	2000-2012 –> 2013-2019	7.2361956	47819.015
Indo-Pacific	(27.75,27.85]	1982-1999 –> 2000-2012	4.6360163	24212.058
Indo-Pacific	(27.75,27.85]	2000-2012 –> 2013-2019	5.8358081	17690.297
Indo-Pacific	(27.85,27.95]	1982-1999 –> 2000-2012	4.8782684	31590.911
Indo-Pacific	(27.85,27.95]	2000-2012 –> 2013-2019	6.0470699	22218.422
Indo-Pacific	(27.95,28.05]	1982-1999 –> 2000-2012	4.1432866	25644.743
Indo-Pacific	(27.95,28.05]	2000-2012 –> 2013-2019	4.6665743	17222.159
Indo-Pacific	(28.05,28.1]	1982-1999 –> 2000-2012	3.4808064	16740.957
Indo-Pacific	(28.05,28.1]	2000-2012 –> 2013-2019	3.9195348	11165.724
Indo-Pacific	(28.1, Inf]	1982-1999 –> 2000-2012	2.5227916	114328.684
Indo-Pacific	(28.1, Inf]	2000-2012 –> 2013-2019	2.8945162	79330.014

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 (rmse) 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
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
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
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
54d5b5b	Donghe-Zhu	2021-02-27
330f064	Donghe-Zhu	2021-02-27
adbc9bc	Donghe-Zhu	2021-02-27
5937141	Donghe-Zhu	2021-02-27
4414bbf	Donghe-Zhu	2021-02-27
a265efb	Donghe-Zhu	2021-02-27
19edd1e	Donghe-Zhu	2021-02-27
f20483f	Donghe-Zhu	2021-02-26
6a2c7b3	Donghe-Zhu	2021-02-25
354c224	Donghe-Zhu	2021-02-24
1a0a88a	Donghe-Zhu	2021-02-24
57f701e	Donghe-Zhu	2021-02-24
06f3149	Donghe-Zhu	2021-02-16
401eab3	Donghe-Zhu	2021-02-15
e3bba84	Donghe-Zhu	2021-02-15
5dce4b1	Donghe-Zhu	2021-02-15
4469a0c	Donghe-Zhu	2021-02-13
5ae6a69	Donghe-Zhu	2021-02-10
05385dc	Donghe-Zhu	2021-02-10
f791ae4	Donghe-Zhu	2021-02-09
f71ae34	Donghe-Zhu	2021-02-09
c011832	Donghe-Zhu	2021-02-09
a145fa7	Donghe-Zhu	2021-02-09
c344e42	Donghe-Zhu	2021-02-08
2f095d7	Donghe-Zhu	2021-02-07
1fad5f1	Donghe-Zhu	2021-02-07
ca03c39	Donghe-Zhu	2021-02-07
e2ffc14	Donghe-Zhu	2021-02-05
cd7c52c	Donghe-Zhu	2021-02-04
bcf84f4	Donghe-Zhu	2021-02-02
a518739	Donghe-Zhu	2021-02-01
61666de	Donghe-Zhu	2021-01-31
865b582	Donghe-Zhu	2021-01-31
3e68089	Donghe-Zhu	2021-01-31
ecf335c	Donghe-Zhu	2021-01-31
a618965	Donghe-Zhu	2021-01-31
59e006e	Donghe-Zhu	2021-01-31
a1c8f87	Donghe-Zhu	2021-01-31
ae5c18f	Donghe-Zhu	2021-01-31
b50fe52	Donghe-Zhu	2021-01-31
ac99ae5	jens-daniel-mueller	2021-01-29
b5bdcaf	Donghe-Zhu	2021-01-29
442010d	Donghe-Zhu	2021-01-29
372adf5	Donghe-Zhu	2021-01-29
af8788e	Donghe-Zhu	2021-01-29
21c91c9	Donghe-Zhu	2021-01-29
eded038	Donghe-Zhu	2021-01-29
541d4dd	Donghe-Zhu	2021-01-29
6a75576	Donghe-Zhu	2021-01-28
16fba40	Donghe-Zhu	2021-01-28
12bc567	Donghe-Zhu	2021-01-27
ceed31b	Donghe-Zhu	2021-01-27
342402d	Donghe-Zhu	2021-01-27
5bad5c2	Donghe-Zhu	2021-01-27
61efb56	Donghe-Zhu	2021-01-25
48f638e	Donghe-Zhu	2021-01-25
c1cec47	Donghe-Zhu	2021-01-25
05ffb0c	Donghe-Zhu	2021-01-25
8b97165	Donghe-Zhu	2021-01-25
c569946	Donghe-Zhu	2021-01-24
a2f0d56	Donghe-Zhu	2021-01-23
28509fc	Donghe-Zhu	2021-01-23
4c28e4a	Donghe-Zhu	2021-01-22
24cc264	jens-daniel-mueller	2021-01-22
7891955	Donghe-Zhu	2021-01-21
d4cf1cb	Donghe-Zhu	2021-01-21
1f3e5b6	jens-daniel-mueller	2021-01-20
0e7bdf1	jens-daniel-mueller	2021-01-15
4571843	jens-daniel-mueller	2021-01-14
b3564aa	jens-daniel-mueller	2021-01-14
8d032c3	jens-daniel-mueller	2021-01-14
17dee1d	jens-daniel-mueller	2021-01-13
7cdea0c	jens-daniel-mueller	2021-01-06
fa85b93	jens-daniel-mueller	2021-01-06
e5cb81a	Donghe-Zhu	2021-01-05
a499f10	Donghe-Zhu	2021-01-05
fb8a752	Donghe-Zhu	2020-12-23
8fae0b2	Donghe-Zhu	2020-12-21
c8b76b3	jens-daniel-mueller	2020-12-19

rm(max_rmse)

6.3 Predictor counts

The number of models where a particular predictor was included were counted for each basin, density slab and compared eras

# calculate cases of predictor used
lm_all_stats <- lm_best_long %>% 
  filter(term != "(Intercept)",
         delta_coeff != 0) %>% 
  group_by(basin, eras, gamma_slab) %>% 
  count(term) %>% 
  ungroup() %>% 
  pivot_wider(values_from = n,
              names_from = term,
              values_fill = 0)

# print table
lm_all_stats %>%
  gt(rowname_col = "gamma_slab",
     groupname_col = c("basin", "eras")) %>% 
  summary_rows(
    groups = TRUE,
    fns = list(total = "sum")
  )

	aou	nitrate	oxygen	phosphate	phosphate_star	sal	silicate	silicate_star	temp
Atlantic - 1982-1999 --> 2000-2012
(-Inf,26]	8	4	2	6	10	7	4	4	3
(26,26.5]	6	8	4	2	6	2	5	5	8
(26.5,26.75]	8	4	2	5	8	2	8	2	7
(26.75,27]	8	4	2	5	8	2	8	2	7
(27,27.25]	8	6	2	4	8	6	5	5	4
(27.25,27.5]	6	9	4	1	6	2	5	5	8
(27.5,27.75]	10	4	0	6	7	8	4	4	2
(27.75,27.85]	10	4	0	6	6	9	5	4	1
(27.85,27.95]	10	0	0	10	9	3	4	3	4
(27.95,28.05]	3	10	7	0	10	3	3	4	6
(28.05,28.1]	10	1	0	9	10	4	4	3	3
(28.1,28.15]	10	2	0	8	10	5	4	4	3
(28.15,28.2]	7	10	3	0	7	3	5	4	7
(28.2, Inf]	6	10	4	0	6	10	4	3	0
total	110.00	76.00	30.00	62.00	111.00	66.00	68.00	52.00	63.00
Atlantic - 2000-2012 --> 2013-2019
(-Inf,26]	10	2	0	6	10	7	4	4	3
(26,26.5]	6	8	4	2	6	2	5	5	8
(26.5,26.75]	8	4	2	5	8	2	8	2	7
(26.75,27]	8	4	2	5	8	2	8	2	7
(27,27.25]	5	4	4	4	8	10	8	2	0
(27.25,27.5]	8	10	2	0	6	4	5	5	6
(27.5,27.75]	7	6	3	4	7	10	4	5	0
(27.75,27.85]	8	6	2	4	6	10	5	5	0
(27.85,27.95]	10	1	0	9	10	4	3	4	3
(27.95,28.05]	4	9	6	1	10	4	3	4	6
(28.05,28.1]	10	2	0	8	10	5	4	4	3
(28.1,28.15]	10	2	0	8	10	5	4	4	3
(28.15,28.2]	7	10	3	0	7	3	5	4	7
(28.2, Inf]	6	10	4	0	6	10	4	3	0
total	107.00	78.00	32.00	56.00	112.00	78.00	70.00	53.00	53.00
Indo-Pacific - 1982-1999 --> 2000-2012
(-Inf,26]	6	4	4	5	9	3	6	3	7
(26,26.5]	8	4	2	5	9	2	7	3	7
(26.5,26.75]	6	6	2	4	10	2	6	4	7
(26.75,27]	10	0	0	10	9	3	3	4	4
(27,27.25]	7	4	3	6	8	1	3	7	8
(27.25,27.5]	6	4	3	5	9	8	3	7	2
(27.5,27.75]	6	4	3	5	9	9	3	7	1
(27.75,27.85]	8	4	2	6	10	2	5	5	6
(27.85,27.95]	3	10	4	0	9	0	4	3	10
(27.95,28.05]	10	4	0	6	10	4	5	5	4
(28.05,28.1]	8	4	2	6	8	4	5	5	6
(28.1, Inf]	5	7	5	3	10	4	6	4	5
total	83.00	55.00	30.00	61.00	110.00	42.00	56.00	57.00	67.00
Indo-Pacific - 2000-2012 --> 2013-2019
(-Inf,26]	6	4	3	5	10	3	6	3	7
(26,26.5]	7	4	3	4	9	1	8	2	8
(26.5,26.75]	5	7	3	3	10	1	6	3	8
(26.75,27]	10	0	0	10	9	3	3	4	4
(27,27.25]	7	4	3	5	9	1	3	7	8
(27.25,27.5]	6	4	3	5	9	8	3	7	2
(27.5,27.75]	8	4	2	4	10	8	5	5	2
(27.75,27.85]	6	6	2	4	10	7	4	5	3
(27.85,27.95]	6	8	4	2	6	0	5	5	10
(27.95,28.05]	10	3	0	7	10	6	3	4	3
(28.05,28.1]	8	4	2	6	10	2	5	5	6
(28.1, Inf]	6	6	4	4	10	4	5	5	6
total	85.00	54.00	29.00	59.00	112.00	44.00	56.00	55.00	67.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
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
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