Last updated: 2021-03-10

Checks: 7 0

Knit directory: emlr_mod_v_XXX/

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

R Markdown file: up-to-date

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

Environment: empty

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

Seed: set.seed(20200707)

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

Session information: recorded

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

Cache: none

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

File paths: relative

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

Repository version: 498be1d

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

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


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

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

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

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

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

1 Required data

Required are:

  • cleaned and prepared GLODAP-based synthetic model subsetting file
GLODAP <-
  read_csv(paste(path_version_data,
                 "GLODAPv2.2020_MLR_fitting_ready.csv",
                 sep = ""))

2 Predictor combinations

Find all possible combinations of following considered predictor variables:

  • sal, temp, oxygen, aou, silicate, phosphate, phosphate_star
# the following code is a workaround to find all predictor combinations
# using the olsrr package and fit all models for one era, slab, and basin

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

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

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

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

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

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

# convert to tibble
lm_all <- as_tibble(lm_all)

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

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

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

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

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

lm_rm <- lm_rm_ph %>%
  select(-lm_coeff_filter) %>%
  unique()

lm_all <- anti_join(lm_all, lm_rm)

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

# remove helper objects
rm(
  i_gamma_slab,
  i_era,
  i_basin,
  GLODAP_basin_era,
  GLODAP_basin_era_slab,
  lm_full,
  lm_rm_ph,
  lm_rm_si,
  lm_rm_o2,
  lm_rm
)

# lm_full <- lm(paste(
#   params_local$MLR_target,
#   paste(params_local$MLR_predictors, collapse = " + "),
#   sep = " ~ "
# ),
# data = GLODAP_basin_era_slab)
# 
# lm_full <-
#   loess(
#     cstar_tref ~ sal + temp + aou + phosphate,
#     span = 0.1,
#     degree = 1,
#     data = GLODAP_basin_era_slab
#   )

3 Apply predictor threshold

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

  • Minimum: 2
  • Maximum: 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:

  • 104

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:

  • 5

The criterion used to select the best models was:

  • aic

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

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

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

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

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

5.3 RMSE tables

5.3.1 per model

lm_best %>%
  kable() %>%
  add_header_above() %>%
  kable_styling() %>%
  scroll_box(width = "100%", height = "400px")
basin gamma_slab n_predictors model nr_obs rmse aic resid_max eras rmse_sum aic_sum
Atlantic (-Inf,26] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 158 1.0492038 475.562558 4.3578441 1982-1999 –> 2000-2012 2.3855522 1198.4564
Atlantic (-Inf,26] 4 cstar_tref ~ sal + oxygen + aou + phosphate 158 1.0492038 475.562558 4.3578441 1982-1999 –> 2000-2012 2.3855522 1198.4564
Atlantic (-Inf,26] 4 cstar_tref ~ sal + oxygen + aou + phosphate_star 158 1.0492038 475.562558 4.3578441 1982-1999 –> 2000-2012 2.3855522 1198.4564
Atlantic (-Inf,26] 4 cstar_tref ~ sal + temp + aou + phosphate 158 1.0547002 477.213652 4.5473106 1982-1999 –> 2000-2012 2.3926132 1200.5942
Atlantic (-Inf,26] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 158 1.0447877 476.229718 3.9913380 1982-1999 –> 2000-2012 2.3809752 1201.0734
Atlantic (-Inf,26] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 72 1.1936365 243.815800 3.2379740 2000-2012 –> 2013-2019 2.2417197 721.0407
Atlantic (-Inf,26] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 72 1.1936365 243.815800 3.2379740 2000-2012 –> 2013-2019 2.2417197 721.0407
Atlantic (-Inf,26] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 72 1.1936365 243.815800 3.2379740 2000-2012 –> 2013-2019 2.2417197 721.0407
Atlantic (-Inf,26] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 72 1.1875291 245.077119 3.1942008 2000-2012 –> 2013-2019 2.2307669 722.8377
Atlantic (-Inf,26] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 72 1.1875291 245.077119 3.1942008 2000-2012 –> 2013-2019 2.2307669 722.8377
Atlantic (26,26.5] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1040 3.5443539 5597.332405 13.4662233 1982-1999 –> 2000-2012 6.9995324 13651.5819
Atlantic (26,26.5] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 1040 3.5443539 5597.332405 13.4662233 1982-1999 –> 2000-2012 6.9995324 13651.5819
Atlantic (26,26.5] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 1040 3.5243146 5587.538970 12.6794726 1982-1999 –> 2000-2012 6.9284246 13598.7593
Atlantic (26,26.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1040 3.5243146 5587.538970 12.6794726 1982-1999 –> 2000-2012 6.9284246 13598.7593
Atlantic (26,26.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1040 3.5243146 5587.538970 12.6794726 1982-1999 –> 2000-2012 6.9284246 13598.7593
Atlantic (26,26.5] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 562 3.6054876 3050.368613 11.1243922 2000-2012 –> 2013-2019 7.1498416 8647.7010
Atlantic (26,26.5] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 562 3.6054876 3050.368613 11.1243922 2000-2012 –> 2013-2019 7.1498416 8647.7010
Atlantic (26,26.5] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 562 3.6014343 3051.104289 11.6005294 2000-2012 –> 2013-2019 7.1257489 8638.6433
Atlantic (26,26.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 562 3.6014343 3051.104289 11.6005294 2000-2012 –> 2013-2019 7.1257489 8638.6433
Atlantic (26,26.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 562 3.6014343 3051.104289 11.6005294 2000-2012 –> 2013-2019 7.1257489 8638.6433
Atlantic (26.5,26.75] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 1360 3.2276124 7062.652885 10.1005384 1982-1999 –> 2000-2012 6.2017135 16296.2793
Atlantic (26.5,26.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1360 3.2276124 7062.652885 10.1005384 1982-1999 –> 2000-2012 6.2017135 16296.2793
Atlantic (26.5,26.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1360 3.2276124 7062.652885 10.1005384 1982-1999 –> 2000-2012 6.2017135 16296.2793
Atlantic (26.5,26.75] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1360 3.2415884 7072.405398 10.0491820 1982-1999 –> 2000-2012 6.2239552 16314.2286
Atlantic (26.5,26.75] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 1360 3.2415884 7072.405398 10.0491820 1982-1999 –> 2000-2012 6.2239552 16314.2286
Atlantic (26.5,26.75] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 646 3.2958735 3390.199829 10.5343935 2000-2012 –> 2013-2019 6.5234860 10452.8527
Atlantic (26.5,26.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 646 3.2958735 3390.199829 10.5343935 2000-2012 –> 2013-2019 6.5234860 10452.8527
Atlantic (26.5,26.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 646 3.2958735 3390.199829 10.5343935 2000-2012 –> 2013-2019 6.5234860 10452.8527
Atlantic (26.5,26.75] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 646 3.3285183 3400.933773 11.1260032 2000-2012 –> 2013-2019 6.5701066 10473.3392
Atlantic (26.5,26.75] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 646 3.3285183 3400.933773 11.1260032 2000-2012 –> 2013-2019 6.5701066 10473.3392
Atlantic (26.75,27] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 2134 1.9756584 8978.118057 13.5842460 1982-1999 –> 2000-2012 3.4606942 18889.7325
Atlantic (26.75,27] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 2134 1.9756584 8978.118057 13.5842460 1982-1999 –> 2000-2012 3.4606942 18889.7325
Atlantic (26.75,27] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 2134 1.9756584 8978.118057 13.5842460 1982-1999 –> 2000-2012 3.4606942 18889.7325
Atlantic (26.75,27] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 2134 2.1557967 9348.538081 19.2601021 1982-1999 –> 2000-2012 3.6479228 19284.1304
Atlantic (26.75,27] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 2134 2.1557967 9348.538081 19.2601021 1982-1999 –> 2000-2012 3.6479228 19284.1304
Atlantic (26.75,27] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 1121 2.5950428 5333.238180 25.4741707 2000-2012 –> 2013-2019 4.6101370 14393.7100
Atlantic (26.75,27] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 1121 2.5950428 5333.238180 25.4741707 2000-2012 –> 2013-2019 4.6101370 14393.7100
Atlantic (26.75,27] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 1121 2.2698743 5035.082409 21.4501194 2000-2012 –> 2013-2019 4.2455326 14013.2005
Atlantic (26.75,27] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1121 2.2698743 5035.082409 21.4501194 2000-2012 –> 2013-2019 4.2455326 14013.2005
Atlantic (26.75,27] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1121 2.2698743 5035.082409 21.4501194 2000-2012 –> 2013-2019 4.2455326 14013.2005
Atlantic (27,27.25] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1717 2.0072484 7279.325305 9.7017253 1982-1999 –> 2000-2012 3.6183821 15798.2280
Atlantic (27,27.25] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 1717 2.0148260 7292.264615 9.3883458 1982-1999 –> 2000-2012 3.6145946 15779.4105
Atlantic (27,27.25] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 1717 1.9602218 7199.914902 13.1630682 1982-1999 –> 2000-2012 3.5598163 15688.5723
Atlantic (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1717 1.9602218 7199.914902 13.1630682 1982-1999 –> 2000-2012 3.5598163 15688.5723
Atlantic (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1717 1.9602218 7199.914902 13.1630682 1982-1999 –> 2000-2012 3.5598163 15688.5723
Atlantic (27,27.25] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 982 2.4198641 4536.404447 14.2305742 2000-2012 –> 2013-2019 4.4271125 11815.7298
Atlantic (27,27.25] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 982 2.4198641 4536.404447 14.2305742 2000-2012 –> 2013-2019 4.4271125 11815.7298
Atlantic (27,27.25] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 982 2.3482249 4479.382903 14.3974130 2000-2012 –> 2013-2019 4.3084468 11679.2978
Atlantic (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 982 2.3482249 4479.382903 14.3974130 2000-2012 –> 2013-2019 4.3084468 11679.2978
Atlantic (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 982 2.3482249 4479.382903 14.3974130 2000-2012 –> 2013-2019 4.3084468 11679.2978
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 1621 2.8012018 7953.616078 19.8213460 1982-1999 –> 2000-2012 5.0696837 18034.3621
Atlantic (27.25,27.5] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 1621 2.8001886 7954.443216 19.6927272 1982-1999 –> 2000-2012 5.0677704 18035.4043
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 1621 2.8012018 7953.616078 19.8213460 1982-1999 –> 2000-2012 5.0696837 18034.3621
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 1621 2.8012018 7953.616078 19.8213460 1982-1999 –> 2000-2012 5.0696837 18034.3621
Atlantic (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1621 2.8001886 7954.443216 19.6927272 1982-1999 –> 2000-2012 5.0677704 18035.4043
Atlantic (27.25,27.5] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 904 3.5908452 4892.765676 20.9092338 2000-2012 –> 2013-2019 6.3910338 12847.2089
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 904 3.6120853 4901.428601 21.4215855 2000-2012 –> 2013-2019 6.4132871 12855.0447
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 904 3.6120853 4901.428601 21.4215855 2000-2012 –> 2013-2019 6.4132871 12855.0447
Atlantic (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 904 3.5908452 4892.765676 20.9092338 2000-2012 –> 2013-2019 6.3910338 12847.2089
Atlantic (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 904 3.5908452 4892.765676 20.9092338 2000-2012 –> 2013-2019 6.3910338 12847.2089
Atlantic (27.5,27.75] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 2079 2.3067859 9391.432216 16.3811627 1982-1999 –> 2000-2012 4.4091130 21475.6236
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 2079 2.3075871 9390.876085 16.0914994 1982-1999 –> 2000-2012 4.4239676 21510.2574
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 2079 2.3075871 9390.876085 16.0914994 1982-1999 –> 2000-2012 4.4239676 21510.2574
Atlantic (27.5,27.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 2079 2.3067859 9391.432216 16.3811627 1982-1999 –> 2000-2012 4.4091130 21475.6236
Atlantic (27.5,27.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 2079 2.3067859 9391.432216 16.3811627 1982-1999 –> 2000-2012 4.4091130 21475.6236
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 1115 2.5792371 5291.143752 17.7062520 2000-2012 –> 2013-2019 4.8868241 14682.0198
Atlantic (27.5,27.75] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 1115 2.5622095 5278.372978 18.9011113 2000-2012 –> 2013-2019 4.8689954 14669.8052
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 1115 2.5792371 5291.143752 17.7062520 2000-2012 –> 2013-2019 4.8868241 14682.0198
Atlantic (27.5,27.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1115 2.5622095 5278.372978 18.9011113 2000-2012 –> 2013-2019 4.8689954 14669.8052
Atlantic (27.5,27.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1115 2.5622095 5278.372978 18.9011113 2000-2012 –> 2013-2019 4.8689954 14669.8052
Atlantic (27.75,27.85] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 722 1.5916813 2736.105232 13.2839690 1982-1999 –> 2000-2012 3.2452848 6434.4569
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 722 1.5992193 2740.927704 12.1339144 1982-1999 –> 2000-2012 3.2778246 6466.0316
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 722 1.5992193 2740.927704 12.1339144 1982-1999 –> 2000-2012 3.2778246 6466.0316
Atlantic (27.75,27.85] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 722 1.5916813 2736.105232 13.2839690 1982-1999 –> 2000-2012 3.2452848 6434.4569
Atlantic (27.75,27.85] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 722 1.5916813 2736.105232 13.2839690 1982-1999 –> 2000-2012 3.2452848 6434.4569
Atlantic (27.75,27.85] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 364 1.6649520 1420.118948 10.6920877 2000-2012 –> 2013-2019 3.2566332 4156.2242
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 364 1.6780668 1423.830913 10.6183742 2000-2012 –> 2013-2019 3.2772860 4164.7586
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 364 1.6780668 1423.830913 10.6183742 2000-2012 –> 2013-2019 3.2772860 4164.7586
Atlantic (27.75,27.85] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 364 1.6649520 1420.118948 10.6920877 2000-2012 –> 2013-2019 3.2566332 4156.2242
Atlantic (27.75,27.85] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 364 1.6649520 1420.118948 10.6920877 2000-2012 –> 2013-2019 3.2566332 4156.2242
Atlantic (27.85,27.95] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 694 2.6858576 3356.830803 20.1208556 1982-1999 –> 2000-2012 5.6386261 8186.0664
Atlantic (27.85,27.95] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 694 2.6858576 3356.830803 20.1208556 1982-1999 –> 2000-2012 5.6386261 8186.0664
Atlantic (27.85,27.95] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 694 2.6878487 3355.859387 20.1028452 1982-1999 –> 2000-2012 5.6412180 8183.4865
Atlantic (27.85,27.95] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 694 2.6878487 3355.859387 20.1028452 1982-1999 –> 2000-2012 5.6412180 8183.4865
Atlantic (27.85,27.95] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 694 2.6878487 3355.859387 20.1028452 1982-1999 –> 2000-2012 5.6412180 8183.4865
Atlantic (27.85,27.95] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 415 3.0207079 2111.276685 17.8599673 2000-2012 –> 2013-2019 5.7065655 5468.1075
Atlantic (27.85,27.95] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 415 3.0207079 2111.276685 17.8599673 2000-2012 –> 2013-2019 5.7065655 5468.1075
Atlantic (27.85,27.95] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 415 3.0207158 2109.278858 17.8578045 2000-2012 –> 2013-2019 5.7085645 5465.1382
Atlantic (27.85,27.95] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 415 3.0207158 2109.278858 17.8578045 2000-2012 –> 2013-2019 5.7085645 5465.1382
Atlantic (27.85,27.95] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 415 3.0207158 2109.278858 17.8578045 2000-2012 –> 2013-2019 5.7085645 5465.1382
Atlantic (27.95,28.05] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 723 3.0495655 3680.073829 12.6814197 1982-1999 –> 2000-2012 6.7958653 9137.1323
Atlantic (27.95,28.05] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 723 3.0495655 3680.073829 12.6814197 1982-1999 –> 2000-2012 6.7958653 9137.1323
Atlantic (27.95,28.05] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 723 3.0495655 3680.073829 12.6814197 1982-1999 –> 2000-2012 6.7958653 9137.1323
Atlantic (27.95,28.05] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 723 3.0943625 3699.160555 18.3736306 1982-1999 –> 2000-2012 6.8421502 9155.0076
Atlantic (27.95,28.05] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 723 3.0943625 3699.160555 18.3736306 1982-1999 –> 2000-2012 6.8421502 9155.0076
Atlantic (27.95,28.05] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 386 2.8070603 1908.230904 12.9707292 2000-2012 –> 2013-2019 5.8566257 5588.3047
Atlantic (27.95,28.05] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 386 2.8070603 1908.230904 12.9707292 2000-2012 –> 2013-2019 5.8566257 5588.3047
Atlantic (27.95,28.05] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 386 2.8070603 1908.230904 12.9707292 2000-2012 –> 2013-2019 5.8566257 5588.3047
Atlantic (27.95,28.05] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 386 2.8196164 1909.676389 13.5039640 2000-2012 –> 2013-2019 5.9139789 5608.8369
Atlantic (27.95,28.05] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 386 2.8196164 1909.676389 13.5039640 2000-2012 –> 2013-2019 5.9139789 5608.8369
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 436 0.2109412 -105.670821 1.4763185 1982-1999 –> 2000-2012 0.3259563 -979.7129
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 436 0.2109412 -105.670821 1.4763185 1982-1999 –> 2000-2012 0.3259563 -979.7129
Atlantic (28.05,28.1] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 436 0.1931541 -180.486583 1.6293925 1982-1999 –> 2000-2012 0.3078192 -1056.1670
Atlantic (28.05,28.1] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 436 0.1931541 -180.486583 1.6293925 1982-1999 –> 2000-2012 0.3078192 -1056.1670
Atlantic (28.05,28.1] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 436 0.1931541 -180.486583 1.6293925 1982-1999 –> 2000-2012 0.3078192 -1056.1670
Atlantic (28.05,28.1] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 229 0.2182829 -31.185279 0.9320609 2000-2012 –> 2013-2019 0.4114370 -211.6719
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 229 0.2362004 2.945556 1.0819569 2000-2012 –> 2013-2019 0.4378010 -142.2187
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 229 0.2362004 2.945556 1.0819569 2000-2012 –> 2013-2019 0.4378010 -142.2187
Atlantic (28.05,28.1] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 229 0.2182829 -31.185279 0.9320609 2000-2012 –> 2013-2019 0.4114370 -211.6719
Atlantic (28.05,28.1] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 229 0.2182829 -31.185279 0.9320609 2000-2012 –> 2013-2019 0.4114370 -211.6719
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 541 0.1394443 -582.345612 1.0824605 1982-1999 –> 2000-2012 0.2743304 -1472.9755
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 541 0.1394443 -582.345612 1.0824605 1982-1999 –> 2000-2012 0.2743304 -1472.9755
Atlantic (28.1,28.15] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 541 0.1354765 -611.579729 0.9957761 1982-1999 –> 2000-2012 0.2695767 -1509.2554
Atlantic (28.1,28.15] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 541 0.1354765 -611.579729 0.9957761 1982-1999 –> 2000-2012 0.2695767 -1509.2554
Atlantic (28.1,28.15] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 541 0.1354765 -611.579729 0.9957761 1982-1999 –> 2000-2012 0.2695767 -1509.2554
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 288 0.1282569 -351.633975 0.4661753 2000-2012 –> 2013-2019 0.2677013 -933.9796
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 288 0.1282569 -351.633975 0.4661753 2000-2012 –> 2013-2019 0.2677013 -933.9796
Atlantic (28.1,28.15] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 288 0.1247025 -365.822431 0.4871925 2000-2012 –> 2013-2019 0.2601790 -977.4022
Atlantic (28.1,28.15] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 288 0.1247025 -365.822431 0.4871925 2000-2012 –> 2013-2019 0.2601790 -977.4022
Atlantic (28.1,28.15] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 288 0.1247025 -365.822431 0.4871925 2000-2012 –> 2013-2019 0.2601790 -977.4022
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1155 0.2973070 489.741021 1.1555289 1982-1999 –> 2000-2012 0.6778125 2072.7436
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 1155 0.2973070 489.741021 1.1555289 1982-1999 –> 2000-2012 0.6778125 2072.7436
Atlantic (28.15,28.2] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 1155 0.2898158 432.790267 1.1774606 1982-1999 –> 2000-2012 0.6681457 1997.9186
Atlantic (28.15,28.2] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1155 0.2898158 432.790267 1.1774606 1982-1999 –> 2000-2012 0.6681457 1997.9186
Atlantic (28.15,28.2] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1155 0.2898158 432.790267 1.1774606 1982-1999 –> 2000-2012 0.6681457 1997.9186
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 636 0.3580823 512.555409 2.3528717 2000-2012 –> 2013-2019 0.6553893 1002.2964
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 636 0.3580823 512.555409 2.3528717 2000-2012 –> 2013-2019 0.6553893 1002.2964
Atlantic (28.15,28.2] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 636 0.3550240 503.644899 2.5256649 2000-2012 –> 2013-2019 0.6448398 936.4352
Atlantic (28.15,28.2] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 636 0.3550240 503.644899 2.5256649 2000-2012 –> 2013-2019 0.6448398 936.4352
Atlantic (28.15,28.2] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 636 0.3550240 503.644899 2.5256649 2000-2012 –> 2013-2019 0.6448398 936.4352
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 955 0.3302057 607.817396 2.0435994 1982-1999 –> 2000-2012 0.8794839 2653.2528
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 955 0.3302057 607.817396 2.0435994 1982-1999 –> 2000-2012 0.8794839 2653.2528
Atlantic (28.2, Inf] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 955 0.2993420 422.390811 2.7766346 1982-1999 –> 2000-2012 0.8468809 2461.9672
Atlantic (28.2, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 955 0.2993420 422.390811 2.7766346 1982-1999 –> 2000-2012 0.8468809 2461.9672
Atlantic (28.2, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 955 0.2993420 422.390811 2.7766346 1982-1999 –> 2000-2012 0.8468809 2461.9672
Atlantic (28.2, Inf] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 552 0.4982897 813.490736 3.3670347 2000-2012 –> 2013-2019 0.7976317 1235.8815
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 552 0.5446404 909.685228 4.6301732 2000-2012 –> 2013-2019 0.8465623 1346.4664
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 552 0.5446404 909.685228 4.6301732 2000-2012 –> 2013-2019 0.8465623 1346.4664
Atlantic (28.2, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 552 0.4982897 813.490736 3.3670347 2000-2012 –> 2013-2019 0.7976317 1235.8815
Atlantic (28.2, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 552 0.4982897 813.490736 3.3670347 2000-2012 –> 2013-2019 0.7976317 1235.8815
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 5288 6.1928478 34304.919955 29.5273536 1982-1999 –> 2000-2012 12.6725136 81601.2906
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 5288 6.1928478 34304.919955 29.5273536 1982-1999 –> 2000-2012 12.6725136 81601.2906
Indo-Pacific (-Inf,26] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 5288 6.1901333 34302.283112 29.4483781 1982-1999 –> 2000-2012 12.6692896 81599.5230
Indo-Pacific (-Inf,26] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 5288 6.1901333 34302.283112 29.4483781 1982-1999 –> 2000-2012 12.6692896 81599.5230
Indo-Pacific (-Inf,26] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 5288 6.1901333 34302.283112 29.4483781 1982-1999 –> 2000-2012 12.6692896 81599.5230
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2815 5.8787943 17975.333992 24.1274152 2000-2012 –> 2013-2019 12.0716422 52280.2539
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 2815 5.8787943 17975.333992 24.1274152 2000-2012 –> 2013-2019 12.0716422 52280.2539
Indo-Pacific (-Inf,26] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 2815 5.8778767 17976.455087 24.0444134 2000-2012 –> 2013-2019 12.0680099 52278.7382
Indo-Pacific (-Inf,26] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 2815 5.8778767 17976.455087 24.0444134 2000-2012 –> 2013-2019 12.0680099 52278.7382
Indo-Pacific (-Inf,26] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 2815 5.8778767 17976.455087 24.0444134 2000-2012 –> 2013-2019 12.0680099 52278.7382
Indo-Pacific (26,26.5] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 5428 4.5387228 31839.277890 43.7366220 1982-1999 –> 2000-2012 8.7811733 75364.3749
Indo-Pacific (26,26.5] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 5428 4.5387228 31839.277890 43.7366220 1982-1999 –> 2000-2012 8.7811733 75364.3749
Indo-Pacific (26,26.5] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 5428 4.5340876 31830.185580 43.0233728 1982-1999 –> 2000-2012 8.7741402 75348.6932
Indo-Pacific (26,26.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 5428 4.5340876 31830.185580 43.0233728 1982-1999 –> 2000-2012 8.7741402 75348.6932
Indo-Pacific (26,26.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 5428 4.5340876 31830.185580 43.0233728 1982-1999 –> 2000-2012 8.7741402 75348.6932
Indo-Pacific (26,26.5] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 2876 4.4815941 16803.612608 46.1900332 2000-2012 –> 2013-2019 9.0203169 48642.8905
Indo-Pacific (26,26.5] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 2876 4.4815941 16803.612608 46.1900332 2000-2012 –> 2013-2019 9.0203169 48642.8905
Indo-Pacific (26,26.5] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 2876 4.4815444 16805.548776 46.3106843 2000-2012 –> 2013-2019 9.0156321 48635.7344
Indo-Pacific (26,26.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 2876 4.4815444 16805.548776 46.3106843 2000-2012 –> 2013-2019 9.0156321 48635.7344
Indo-Pacific (26,26.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 2876 4.4815444 16805.548776 46.3106843 2000-2012 –> 2013-2019 9.0156321 48635.7344
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 4269 3.9438120 23844.294920 20.4830181 1982-1999 –> 2000-2012 7.6993981 56369.6152
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 4269 3.9438120 23844.294920 20.4830181 1982-1999 –> 2000-2012 7.6993981 56369.6152
Indo-Pacific (26.5,26.75] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 4269 3.8963414 23742.901689 26.7023232 1982-1999 –> 2000-2012 7.6360524 56220.0000
Indo-Pacific (26.5,26.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 4269 3.8963414 23742.901689 26.7023232 1982-1999 –> 2000-2012 7.6360524 56220.0000
Indo-Pacific (26.5,26.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 4269 3.8963414 23742.901689 26.7023232 1982-1999 –> 2000-2012 7.6360524 56220.0000
Indo-Pacific (26.5,26.75] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 2324 4.2544954 13341.419522 23.3613336 2000-2012 –> 2013-2019 8.1508368 37084.3212
Indo-Pacific (26.5,26.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 2324 4.2544954 13341.419522 23.3613336 2000-2012 –> 2013-2019 8.1508368 37084.3212
Indo-Pacific (26.5,26.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 2324 4.2544954 13341.419522 23.3613336 2000-2012 –> 2013-2019 8.1508368 37084.3212
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 2324 4.3137865 13403.747333 23.1641121 2000-2012 –> 2013-2019 8.2596379 37252.4562
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 2324 4.3137865 13403.747333 23.1641121 2000-2012 –> 2013-2019 8.2596379 37252.4562
Indo-Pacific (26.75,27] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 4557 4.1322056 25879.252031 21.7852477 1982-1999 –> 2000-2012 7.6563812 60336.4988
Indo-Pacific (26.75,27] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 4557 4.1556325 25928.776454 23.0206730 1982-1999 –> 2000-2012 7.6982624 60451.1786
Indo-Pacific (26.75,27] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 4557 4.1556325 25928.776454 23.0206730 1982-1999 –> 2000-2012 7.6982624 60451.1786
Indo-Pacific (26.75,27] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 4557 4.1322056 25879.252031 21.7852477 1982-1999 –> 2000-2012 7.6563812 60336.4988
Indo-Pacific (26.75,27] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 4557 4.1322056 25879.252031 21.7852477 1982-1999 –> 2000-2012 7.6563812 60336.4988
Indo-Pacific (26.75,27] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 2460 4.5918239 14496.621995 16.9140190 2000-2012 –> 2013-2019 8.7240295 40375.8740
Indo-Pacific (26.75,27] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 2460 4.5918239 14496.621995 16.9140190 2000-2012 –> 2013-2019 8.7240295 40375.8740
Indo-Pacific (26.75,27] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 2460 4.5918239 14496.621995 16.9140190 2000-2012 –> 2013-2019 8.7240295 40375.8740
Indo-Pacific (26.75,27] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 2460 4.5994232 14502.757680 16.6240534 2000-2012 –> 2013-2019 8.7356513 40388.8774
Indo-Pacific (26.75,27] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 2460 4.5994232 14502.757680 16.6240534 2000-2012 –> 2013-2019 8.7356513 40388.8774
Indo-Pacific (27,27.25] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 4267 3.4036330 22578.034898 30.6158677 1982-1999 –> 2000-2012 6.1906243 50781.9909
Indo-Pacific (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 4267 3.4036330 22578.034898 30.6158677 1982-1999 –> 2000-2012 6.1906243 50781.9909
Indo-Pacific (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 4267 3.4036330 22578.034898 30.6158677 1982-1999 –> 2000-2012 6.1906243 50781.9909
Indo-Pacific (27,27.25] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 4267 3.4041802 22577.406832 30.9255112 1982-1999 –> 2000-2012 6.2015505 50822.2366
Indo-Pacific (27,27.25] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 4267 3.4041802 22577.406832 30.9255112 1982-1999 –> 2000-2012 6.2015505 50822.2366
Indo-Pacific (27,27.25] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 2344 3.9472749 13104.727079 23.9342640 2000-2012 –> 2013-2019 7.3509079 35682.7620
Indo-Pacific (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 2344 3.9472749 13104.727079 23.9342640 2000-2012 –> 2013-2019 7.3509079 35682.7620
Indo-Pacific (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 2344 3.9472749 13104.727079 23.9342640 2000-2012 –> 2013-2019 7.3509079 35682.7620
Indo-Pacific (27,27.25] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 2344 3.9533412 13109.926205 24.4569335 2000-2012 –> 2013-2019 7.3575214 35687.3330
Indo-Pacific (27,27.25] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 2344 3.9533412 13109.926205 24.4569335 2000-2012 –> 2013-2019 7.3575214 35687.3330
Indo-Pacific (27.25,27.5] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 3520 2.8617081 17407.314891 30.1089756 1982-1999 –> 2000-2012 5.6400680 40449.8156
Indo-Pacific (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 3520 2.8617081 17407.314891 30.1089756 1982-1999 –> 2000-2012 5.6400680 40449.8156
Indo-Pacific (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 3520 2.8617081 17407.314891 30.1089756 1982-1999 –> 2000-2012 5.6400680 40449.8156
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 3520 2.8774108 17443.838876 30.1555714 1982-1999 –> 2000-2012 5.6901545 40600.3746
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 3520 2.8774108 17443.838876 30.1555714 1982-1999 –> 2000-2012 5.6901545 40600.3746
Indo-Pacific (27.25,27.5] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 1879 3.3074238 9843.576259 29.1273995 2000-2012 –> 2013-2019 6.1691319 27250.8912
Indo-Pacific (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1879 3.3074238 9843.576259 29.1273995 2000-2012 –> 2013-2019 6.1691319 27250.8912
Indo-Pacific (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1879 3.3074238 9843.576259 29.1273995 2000-2012 –> 2013-2019 6.1691319 27250.8912
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 1879 3.3526141 9892.575280 28.9794819 2000-2012 –> 2013-2019 6.2300248 27336.4142
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 1879 3.3526141 9892.575280 28.9794819 2000-2012 –> 2013-2019 6.2300248 27336.4142
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 3170 1.8808520 13015.205979 9.5198488 1982-1999 –> 2000-2012 3.5317875 30323.2534
Indo-Pacific (27.5,27.75] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 3170 1.8800469 13014.491400 9.6798100 1982-1999 –> 2000-2012 3.5305329 30322.0864
Indo-Pacific (27.5,27.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 3170 1.8800469 13014.491400 9.6798100 1982-1999 –> 2000-2012 3.5305329 30322.0864
Indo-Pacific (27.5,27.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 3170 1.8800469 13014.491400 9.6798100 1982-1999 –> 2000-2012 3.5305329 30322.0864
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 3170 1.8817537 13018.244743 9.5293548 1982-1999 –> 2000-2012 3.5322649 30323.9770
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 1683 2.1449046 7358.725091 12.7298567 2000-2012 –> 2013-2019 4.0257566 20373.9311
Indo-Pacific (27.5,27.75] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 1683 2.1441647 7359.563742 12.6747027 2000-2012 –> 2013-2019 4.0242115 20374.0551
Indo-Pacific (27.5,27.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1683 2.1441647 7359.563742 12.6747027 2000-2012 –> 2013-2019 4.0242115 20374.0551
Indo-Pacific (27.5,27.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1683 2.1441647 7359.563742 12.6747027 2000-2012 –> 2013-2019 4.0242115 20374.0551
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 1683 2.1454909 7359.645130 12.7476542 2000-2012 –> 2013-2019 4.0272446 20377.8899
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 1216 1.2257896 3959.977795 5.7133287 1982-1999 –> 2000-2012 2.4231295 9372.3398
Indo-Pacific (27.75,27.85] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 1216 1.1484687 3803.519113 5.5653347 1982-1999 –> 2000-2012 2.2742935 9009.9652
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 1216 1.2257896 3959.977795 5.7133287 1982-1999 –> 2000-2012 2.4231295 9372.3398
Indo-Pacific (27.75,27.85] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1216 1.1484687 3803.519113 5.5653347 1982-1999 –> 2000-2012 2.2742935 9009.9652
Indo-Pacific (27.75,27.85] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1216 1.1484687 3803.519113 5.5653347 1982-1999 –> 2000-2012 2.2742935 9009.9652
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 691 1.3286664 2367.703966 8.2042737 2000-2012 –> 2013-2019 2.5544561 6327.6818
Indo-Pacific (27.75,27.85] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 691 1.3067937 2346.763836 7.0585944 2000-2012 –> 2013-2019 2.4552624 6150.2829
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 691 1.3286664 2367.703966 8.2042737 2000-2012 –> 2013-2019 2.5544561 6327.6818
Indo-Pacific (27.75,27.85] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 691 1.3067937 2346.763836 7.0585944 2000-2012 –> 2013-2019 2.4552624 6150.2829
Indo-Pacific (27.75,27.85] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 691 1.3067937 2346.763836 7.0585944 2000-2012 –> 2013-2019 2.4552624 6150.2829
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 1797 1.2870458 6020.609331 8.7933519 1982-1999 –> 2000-2012 2.3725611 13632.6372
Indo-Pacific (27.85,27.95] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 1797 1.1268453 5544.867661 6.8163957 1982-1999 –> 2000-2012 2.0904038 12555.6225
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 1797 1.2870458 6020.609331 8.7933519 1982-1999 –> 2000-2012 2.3725611 13632.6372
Indo-Pacific (27.85,27.95] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1797 1.1268453 5544.867661 6.8163957 1982-1999 –> 2000-2012 2.0904038 12555.6225
Indo-Pacific (27.85,27.95] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1797 1.1268453 5544.867661 6.8163957 1982-1999 –> 2000-2012 2.0904038 12555.6225
Indo-Pacific (27.85,27.95] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 987 1.3133146 3355.006629 7.3269651 2000-2012 –> 2013-2019 2.4401599 8899.8743
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 987 1.4361554 3529.512839 8.9979380 2000-2012 –> 2013-2019 2.7232012 9550.1222
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ sal + temp + oxygen + aou + silicate 987 1.4185912 3505.221899 7.0604103 2000-2012 –> 2013-2019 2.6559788 9384.4175
Indo-Pacific (27.85,27.95] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 987 1.3133146 3355.006629 7.3269651 2000-2012 –> 2013-2019 2.4401599 8899.8743
Indo-Pacific (27.85,27.95] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 987 1.3133146 3355.006629 7.3269651 2000-2012 –> 2013-2019 2.4401599 8899.8743
Indo-Pacific (27.95,28.05] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 1162 0.5749179 2027.213849 2.9262058 1982-1999 –> 2000-2012 1.1061869 4778.4935
Indo-Pacific (27.95,28.05] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1162 0.5749179 2027.213849 2.9262058 1982-1999 –> 2000-2012 1.1061869 4778.4935
Indo-Pacific (27.95,28.05] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1162 0.5749179 2027.213849 2.9262058 1982-1999 –> 2000-2012 1.1061869 4778.4935
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 1162 0.5749226 2025.232979 2.9259808 1982-1999 –> 2000-2012 1.1074497 4782.7387
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 1162 0.5749226 2025.232979 2.9259808 1982-1999 –> 2000-2012 1.1074497 4782.7387
Indo-Pacific (27.95,28.05] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 648 0.5887663 1168.419405 2.7201223 2000-2012 –> 2013-2019 1.1636841 3195.6333
Indo-Pacific (27.95,28.05] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 648 0.5887663 1168.419405 2.7201223 2000-2012 –> 2013-2019 1.1636841 3195.6333
Indo-Pacific (27.95,28.05] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 648 0.5887663 1168.419405 2.7201223 2000-2012 –> 2013-2019 1.1636841 3195.6333
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 648 0.5922331 1174.028301 2.7169734 2000-2012 –> 2013-2019 1.1671557 3199.2613
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 648 0.5922331 1174.028301 2.7169734 2000-2012 –> 2013-2019 1.1671557 3199.2613
Indo-Pacific (28.05,28.1] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 1015 0.4542785 1294.684324 2.4180424 1982-1999 –> 2000-2012 0.9081081 3221.2974
Indo-Pacific (28.05,28.1] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1015 0.4542785 1294.684324 2.4180424 1982-1999 –> 2000-2012 0.9081081 3221.2974
Indo-Pacific (28.05,28.1] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1015 0.4542785 1294.684324 2.4180424 1982-1999 –> 2000-2012 0.9081081 3221.2974
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 1015 0.4546336 1294.270504 2.4086220 1982-1999 –> 2000-2012 0.9095606 3226.2210
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 1015 0.4546336 1294.270504 2.4086220 1982-1999 –> 2000-2012 0.9095606 3226.2210
Indo-Pacific (28.05,28.1] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 565 0.4745332 777.071782 2.6197171 2000-2012 –> 2013-2019 0.9288117 2071.7561
Indo-Pacific (28.05,28.1] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 565 0.4745332 777.071782 2.6197171 2000-2012 –> 2013-2019 0.9288117 2071.7561
Indo-Pacific (28.05,28.1] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 565 0.4745332 777.071782 2.6197171 2000-2012 –> 2013-2019 0.9288117 2071.7561
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 565 0.4756293 777.678976 2.6025948 2000-2012 –> 2013-2019 0.9302630 2071.9495
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 565 0.4756293 777.678976 2.6025948 2000-2012 –> 2013-2019 0.9302630 2071.9495
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 9318 0.7764916 21742.998933 4.8420893 1982-1999 –> 2000-2012 1.5383986 51524.1449
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 9318 0.7764916 21742.998933 4.8420893 1982-1999 –> 2000-2012 1.5383986 51524.1449
Indo-Pacific (28.1, Inf] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 9318 0.7739529 21683.971466 4.9061100 1982-1999 –> 2000-2012 1.5356708 51460.6751
Indo-Pacific (28.1, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 9318 0.7739529 21683.971466 4.9061100 1982-1999 –> 2000-2012 1.5356708 51460.6751
Indo-Pacific (28.1, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 9318 0.7739529 21683.971466 4.9061100 1982-1999 –> 2000-2012 1.5356708 51460.6751
Indo-Pacific (28.1, Inf] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 5093 0.7692293 11796.846452 4.3310811 2000-2012 –> 2013-2019 1.5431823 33480.8179
Indo-Pacific (28.1, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 5093 0.7692293 11796.846452 4.3310811 2000-2012 –> 2013-2019 1.5431823 33480.8179
Indo-Pacific (28.1, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 5093 0.7692293 11796.846452 4.3310811 2000-2012 –> 2013-2019 1.5431823 33480.8179
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 5093 0.7692867 11795.605640 4.3360410 2000-2012 –> 2013-2019 1.5437859 33490.7275
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 5093 0.7692867 11795.605640 4.3360410 2000-2012 –> 2013-2019 1.5437859 33490.7275
SO (26,26.5] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 103 2.6727772 510.823665 11.1616695 1982-1999 –> 2000-2012 6.1036469 1210.9758
SO (26,26.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 103 2.6727772 510.823665 11.1616695 1982-1999 –> 2000-2012 6.1036469 1210.9758
SO (26,26.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 103 2.6727772 510.823665 11.1616695 1982-1999 –> 2000-2012 6.1036469 1210.9758
SO (26,26.5] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 103 2.7247157 512.788337 11.4146613 1982-1999 –> 2000-2012 6.2078603 1214.8418
SO (26,26.5] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 103 2.7247157 512.788337 11.4146613 1982-1999 –> 2000-2012 6.2078603 1214.8418
SO (26,26.5] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 47 3.0856552 255.296039 8.0922982 2000-2012 –> 2013-2019 5.7584324 766.1197
SO (26,26.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 47 3.0856552 255.296039 8.0922982 2000-2012 –> 2013-2019 5.7584324 766.1197
SO (26,26.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 47 3.0856552 255.296039 8.0922982 2000-2012 –> 2013-2019 5.7584324 766.1197
SO (26,26.5] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 47 3.0875443 253.353570 8.0376410 2000-2012 –> 2013-2019 5.8122599 766.1419
SO (26,26.5] 4 cstar_tref ~ sal + temp + silicate + phosphate_star 47 3.1968068 254.622550 7.8513989 2000-2012 –> 2013-2019 5.9273003 765.8473
SO (26.5,26.75] 5 cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 312 3.2153502 1628.209882 10.5226323 1982-1999 –> 2000-2012 6.7106987 3885.3220
SO (26.5,26.75] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 312 3.1290661 1613.236019 9.9544018 1982-1999 –> 2000-2012 6.6028383 3867.1469
SO (26.5,26.75] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 312 3.2153502 1628.209882 10.5226323 1982-1999 –> 2000-2012 6.7106987 3885.3220
SO (26.5,26.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 312 3.1290661 1613.236019 9.9544018 1982-1999 –> 2000-2012 6.6028383 3867.1469
SO (26.5,26.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 312 3.1290661 1613.236019 9.9544018 1982-1999 –> 2000-2012 6.6028383 3867.1469
SO (26.5,26.75] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 133 3.2742828 708.939942 10.4020030 2000-2012 –> 2013-2019 6.4033488 2322.1760
SO (26.5,26.75] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 133 3.3899758 716.176510 10.0861110 2000-2012 –> 2013-2019 6.6053260 2344.3864
SO (26.5,26.75] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 133 3.3899758 716.176510 10.0861110 2000-2012 –> 2013-2019 6.6053260 2344.3864
SO (26.5,26.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 133 3.2742828 708.939942 10.4020030 2000-2012 –> 2013-2019 6.4033488 2322.1760
SO (26.5,26.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 133 3.2742828 708.939942 10.4020030 2000-2012 –> 2013-2019 6.4033488 2322.1760
SO (26.75,27] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 1061 2.2880755 4783.390414 9.3213970 1982-1999 –> 2000-2012 4.7000755 11766.5567
SO (26.75,27] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1061 2.2880755 4783.390414 9.3213970 1982-1999 –> 2000-2012 4.7000755 11766.5567
SO (26.75,27] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1061 2.2880755 4783.390414 9.3213970 1982-1999 –> 2000-2012 4.7000755 11766.5567
SO (26.75,27] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1061 2.3191394 4810.005845 10.6756905 1982-1999 –> 2000-2012 4.7542973 11820.1248
SO (26.75,27] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 1061 2.3191394 4810.005845 10.6756905 1982-1999 –> 2000-2012 4.7542973 11820.1248
SO (26.75,27] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 550 2.2521270 2469.895006 11.5622287 2000-2012 –> 2013-2019 4.5402025 7253.2854
SO (26.75,27] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 550 2.2521270 2469.895006 11.5622287 2000-2012 –> 2013-2019 4.5402025 7253.2854
SO (26.75,27] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 550 2.2521270 2469.895006 11.5622287 2000-2012 –> 2013-2019 4.5402025 7253.2854
SO (26.75,27] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 550 2.2748295 2478.928005 12.6508529 2000-2012 –> 2013-2019 4.5939689 7288.9339
SO (26.75,27] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 550 2.2748295 2478.928005 12.6508529 2000-2012 –> 2013-2019 4.5939689 7288.9339
SO (27,27.25] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 2513 1.6118710 9546.975590 12.9579746 1982-1999 –> 2000-2012 3.1585654 22318.3580
SO (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 2513 1.6118710 9546.975590 12.9579746 1982-1999 –> 2000-2012 3.1585654 22318.3580
SO (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 2513 1.6118710 9546.975590 12.9579746 1982-1999 –> 2000-2012 3.1585654 22318.3580
SO (27,27.25] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 2513 1.6615541 9697.553575 16.3323879 1982-1999 –> 2000-2012 3.2687331 22730.7024
SO (27,27.25] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 2513 1.6615541 9697.553575 16.3323879 1982-1999 –> 2000-2012 3.2687331 22730.7024
SO (27,27.25] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 1286 1.7553735 5112.727102 7.7602872 2000-2012 –> 2013-2019 3.3672445 14659.7027
SO (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1286 1.7553735 5112.727102 7.7602872 2000-2012 –> 2013-2019 3.3672445 14659.7027
SO (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1286 1.7553735 5112.727102 7.7602872 2000-2012 –> 2013-2019 3.3672445 14659.7027
SO (27,27.25] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 1286 1.8504129 5246.341361 11.4263441 2000-2012 –> 2013-2019 3.5119670 14943.8949
SO (27,27.25] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 1286 1.8504129 5246.341361 11.4263441 2000-2012 –> 2013-2019 3.5119670 14943.8949
SO (27.25,27.5] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 2731 1.6887194 10628.169039 10.6094178 1982-1999 –> 2000-2012 3.1238256 24187.7586
SO (27.25,27.5] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2731 1.6887194 10626.169048 10.6093166 1982-1999 –> 2000-2012 3.1260696 24195.6454
SO (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 2731 1.6887194 10628.169039 10.6094178 1982-1999 –> 2000-2012 3.1238256 24187.7586
SO (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 2731 1.6887194 10628.169039 10.6094178 1982-1999 –> 2000-2012 3.1238256 24187.7586
SO (27.25,27.5] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 2731 1.6893569 10628.230775 10.5492429 1982-1999 –> 2000-2012 3.1276338 24202.6106
SO (27.25,27.5] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 1453 1.8445309 5918.561145 8.9141838 2000-2012 –> 2013-2019 3.5332503 16546.7302
SO (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1453 1.8445309 5918.561145 8.9141838 2000-2012 –> 2013-2019 3.5332503 16546.7302
SO (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1453 1.8445309 5918.561145 8.9141838 2000-2012 –> 2013-2019 3.5332503 16546.7302
SO (27.25,27.5] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 1453 1.8467306 5920.024690 8.8926590 2000-2012 –> 2013-2019 3.5356934 16546.9808
SO (27.25,27.5] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 1453 1.8467306 5920.024690 8.8926590 2000-2012 –> 2013-2019 3.5356934 16546.9808
SO (27.5,27.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 3391 1.5320138 12530.327695 9.1542418 1982-1999 –> 2000-2012 2.8147095 27705.5603
SO (27.5,27.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 3391 1.5317818 12529.300408 9.1490232 1982-1999 –> 2000-2012 2.8137341 27699.2635
SO (27.5,27.75] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 3391 1.5315636 12530.334268 9.1313823 1982-1999 –> 2000-2012 2.8134965 27702.1603
SO (27.5,27.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 3391 1.5315636 12530.334268 9.1313823 1982-1999 –> 2000-2012 2.8134965 27702.1603
SO (27.5,27.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 3391 1.5315636 12530.334268 9.1313823 1982-1999 –> 2000-2012 2.8134965 27702.1603
SO (27.5,27.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1821 1.7554240 7231.165465 10.3901005 2000-2012 –> 2013-2019 3.2874378 19761.4932
SO (27.5,27.75] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 1821 1.7554240 7231.165465 10.3901005 2000-2012 –> 2013-2019 3.2874378 19761.4932
SO (27.5,27.75] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 1821 1.7554240 7231.165465 10.3901005 2000-2012 –> 2013-2019 3.2874378 19761.4932
SO (27.5,27.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 1821 1.7554279 7231.173500 10.3915870 2000-2012 –> 2013-2019 3.2872096 19760.4739
SO (27.5,27.75] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 1821 1.7560943 7232.555962 10.3552020 2000-2012 –> 2013-2019 3.2880132 19762.4635
SO (27.75,27.85] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1415 1.4740467 5127.668440 13.4027185 1982-1999 –> 2000-2012 2.7633662 11181.3903
SO (27.75,27.85] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 1415 1.4740467 5127.668440 13.4027185 1982-1999 –> 2000-2012 2.7633662 11181.3903
SO (27.75,27.85] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 1415 1.4577478 5098.202260 12.9139862 1982-1999 –> 2000-2012 2.7299117 11105.5671
SO (27.75,27.85] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1415 1.4577478 5098.202260 12.9139862 1982-1999 –> 2000-2012 2.7299117 11105.5671
SO (27.75,27.85] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1415 1.4577478 5098.202260 12.9139862 1982-1999 –> 2000-2012 2.7299117 11105.5671
SO (27.75,27.85] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 722 1.5913038 2733.762782 8.6142494 2000-2012 –> 2013-2019 3.0653505 7861.4312
SO (27.75,27.85] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 722 1.5913038 2733.762782 8.6142494 2000-2012 –> 2013-2019 3.0653505 7861.4312
SO (27.75,27.85] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 722 1.5864476 2731.349388 8.7949376 2000-2012 –> 2013-2019 3.0441955 7829.5516
SO (27.75,27.85] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 722 1.5864476 2731.349388 8.7949376 2000-2012 –> 2013-2019 3.0441955 7829.5516
SO (27.75,27.85] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 722 1.5864476 2731.349388 8.7949376 2000-2012 –> 2013-2019 3.0441955 7829.5516
SO (27.85,27.95] 5 cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 1399 1.4232287 4971.682542 7.2306924 1982-1999 –> 2000-2012 2.7624199 11326.6650
SO (27.85,27.95] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 1399 1.2836451 4684.861093 8.5215752 1982-1999 –> 2000-2012 2.5448741 10819.5605
SO (27.85,27.95] 5 cstar_tref ~ sal + temp + oxygen + aou + silicate 1399 1.2912832 4699.460807 8.0313297 1982-1999 –> 2000-2012 2.5525237 10832.1940
SO (27.85,27.95] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1399 1.2836451 4684.861093 8.5215752 1982-1999 –> 2000-2012 2.5448741 10819.5605
SO (27.85,27.95] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1399 1.2836451 4684.861093 8.5215752 1982-1999 –> 2000-2012 2.5448741 10819.5605
SO (27.85,27.95] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 709 1.4880087 2591.623013 9.0787780 2000-2012 –> 2013-2019 2.7716538 7276.4841
SO (27.85,27.95] 5 cstar_tref ~ sal + temp + oxygen + aou + silicate 709 1.4996937 2600.714760 8.5827164 2000-2012 –> 2013-2019 2.7909769 7300.1756
SO (27.85,27.95] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 709 1.4880087 2591.623013 9.0787780 2000-2012 –> 2013-2019 2.7716538 7276.4841
SO (27.85,27.95] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 709 1.4880087 2591.623013 9.0787780 2000-2012 –> 2013-2019 2.7716538 7276.4841
SO (27.85,27.95] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 709 1.5409988 2639.241700 8.7822837 2000-2012 –> 2013-2019 2.9762393 7634.4399
SO (27.95,28.05] 5 cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 1786 1.5620816 6675.629386 9.1134438 1982-1999 –> 2000-2012 2.8752317 15102.4909
SO (27.95,28.05] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 1786 1.5465396 6641.911675 8.3764032 1982-1999 –> 2000-2012 2.8450920 15015.1697
SO (27.95,28.05] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 1786 1.5620816 6675.629386 9.1134438 1982-1999 –> 2000-2012 2.8752317 15102.4909
SO (27.95,28.05] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1786 1.5465396 6641.911675 8.3764032 1982-1999 –> 2000-2012 2.8450920 15015.1697
SO (27.95,28.05] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1786 1.5465396 6641.911675 8.3764032 1982-1999 –> 2000-2012 2.8450920 15015.1697
SO (27.95,28.05] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 986 1.6280632 3775.281999 8.1157057 2000-2012 –> 2013-2019 3.1746028 10417.1937
SO (27.95,28.05] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 986 1.6936583 3851.175591 8.4129789 2000-2012 –> 2013-2019 3.2557399 10526.8050
SO (27.95,28.05] 5 cstar_tref ~ sal + temp + oxygen + aou + silicate 986 1.6282041 3773.452640 8.1160135 2000-2012 –> 2013-2019 3.1902005 10448.8871
SO (27.95,28.05] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 986 1.6280632 3775.281999 8.1157057 2000-2012 –> 2013-2019 3.1746028 10417.1937
SO (27.95,28.05] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 986 1.6280632 3775.281999 8.1157057 2000-2012 –> 2013-2019 3.1746028 10417.1937
SO (28.05,28.1] 5 cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 1071 1.0487379 3155.298539 9.0027805 1982-1999 –> 2000-2012 2.0069805 7336.6873
SO (28.05,28.1] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 1071 1.0487379 3155.298539 9.0027805 1982-1999 –> 2000-2012 2.0069805 7336.6873
SO (28.05,28.1] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 1071 1.0487379 3155.298539 9.0027805 1982-1999 –> 2000-2012 2.0069805 7336.6873
SO (28.05,28.1] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1071 1.0473850 3154.533537 9.0409967 1982-1999 –> 2000-2012 2.0055492 7337.6746
SO (28.05,28.1] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1071 1.0473850 3154.533537 9.0409967 1982-1999 –> 2000-2012 2.0055492 7337.6746
SO (28.05,28.1] 5 cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 606 1.2283070 1982.985334 10.1816140 2000-2012 –> 2013-2019 2.2770449 5138.2839
SO (28.05,28.1] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 606 1.2282847 1984.963330 10.1968303 2000-2012 –> 2013-2019 2.2756697 5139.4969
SO (28.05,28.1] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 606 1.2283070 1982.985334 10.1816140 2000-2012 –> 2013-2019 2.2770449 5138.2839
SO (28.05,28.1] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 606 1.2283070 1982.985334 10.1816140 2000-2012 –> 2013-2019 2.2770449 5138.2839
SO (28.05,28.1] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 606 1.2282847 1984.963330 10.1968303 2000-2012 –> 2013-2019 2.2756697 5139.4969
SO (28.1, Inf] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 5309 1.1018498 16110.133861 6.1005720 1982-1999 –> 2000-2012 2.0809310 36696.9215
SO (28.1, Inf] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 5309 1.1018498 16110.133861 6.1005720 1982-1999 –> 2000-2012 2.0809310 36696.9215
SO (28.1, Inf] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 5309 1.1012492 16106.344648 6.1101131 1982-1999 –> 2000-2012 2.0799840 36689.9248
SO (28.1, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 5309 1.1012492 16106.344648 6.1101131 1982-1999 –> 2000-2012 2.0799840 36689.9248
SO (28.1, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 5309 1.1012492 16106.344648 6.1101131 1982-1999 –> 2000-2012 2.0799840 36689.9248
SO (28.1, Inf] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 2711 1.1918719 8659.181813 6.0524338 2000-2012 –> 2013-2019 2.2937217 24769.3157
SO (28.1, Inf] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 2711 1.1918719 8659.181813 6.0524338 2000-2012 –> 2013-2019 2.2937217 24769.3157
SO (28.1, Inf] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 2711 1.1918598 8661.126810 6.0499653 2000-2012 –> 2013-2019 2.2931091 24767.4715
SO (28.1, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 2711 1.1918598 8661.126810 6.0499653 2000-2012 –> 2013-2019 2.2931091 24767.4715
SO (28.1, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 2711 1.1918598 8661.126810 6.0499653 2000-2012 –> 2013-2019 2.2931091 24767.4715
SO (28.1,28.15] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 395 0.4412700 490.663828 1.3656130 1982-1999 –> 2000-2012 0.8604524 1079.2317
SO (28.1,28.15] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 395 0.4688350 536.532945 1.2066978 1982-1999 –> 2000-2012 0.8955957 1141.7707
SO (28.1,28.15] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 395 0.4688350 536.532945 1.2066978 1982-1999 –> 2000-2012 0.8955957 1141.7707
SO (28.1,28.15] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 395 0.4412700 490.663828 1.3656130 1982-1999 –> 2000-2012 0.8604524 1079.2317
SO (28.1,28.15] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 395 0.4412700 490.663828 1.3656130 1982-1999 –> 2000-2012 0.8604524 1079.2317
SO (28.1,28.15] 5 cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 209 0.4863454 305.806779 1.9335059 2000-2012 –> 2013-2019 0.9551804 842.3397
SO (28.1,28.15] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 209 0.4627191 286.990868 2.0315187 2000-2012 –> 2013-2019 0.9039892 777.6547
SO (28.1,28.15] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 209 0.4863454 305.806779 1.9335059 2000-2012 –> 2013-2019 0.9551804 842.3397
SO (28.1,28.15] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 209 0.4627191 286.990868 2.0315187 2000-2012 –> 2013-2019 0.9039892 777.6547
SO (28.1,28.15] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 209 0.4627191 286.990868 2.0315187 2000-2012 –> 2013-2019 0.9039892 777.6547
SO (28.15,28.2] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 385 0.2908445 157.658441 1.0402887 1982-1999 –> 2000-2012 0.5915367 387.0178
SO (28.15,28.2] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 385 0.2958284 168.741280 1.2331036 1982-1999 –> 2000-2012 0.5978993 400.5932
SO (28.15,28.2] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 385 0.2958284 168.741280 1.2331036 1982-1999 –> 2000-2012 0.5978993 400.5932
SO (28.15,28.2] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 385 0.2908445 157.658441 1.0402887 1982-1999 –> 2000-2012 0.5915367 387.0178
SO (28.15,28.2] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 385 0.2908445 157.658441 1.0402887 1982-1999 –> 2000-2012 0.5915367 387.0178
SO (28.15,28.2] 5 cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 183 0.3300672 127.635476 0.8682035 2000-2012 –> 2013-2019 0.6258955 296.3768
SO (28.15,28.2] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 183 0.3212016 119.670270 0.7924252 2000-2012 –> 2013-2019 0.6120461 277.3287
SO (28.15,28.2] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 183 0.3300672 127.635476 0.8682035 2000-2012 –> 2013-2019 0.6258955 296.3768
SO (28.15,28.2] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 183 0.3212016 119.670270 0.7924252 2000-2012 –> 2013-2019 0.6120461 277.3287
SO (28.15,28.2] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 183 0.3212016 119.670270 0.7924252 2000-2012 –> 2013-2019 0.6120461 277.3287
SO (28.2, Inf] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 1738 0.3884149 1661.042850 2.0901512 1982-1999 –> 2000-2012 0.7108438 3096.2954
SO (28.2, Inf] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1738 0.3910307 1682.373485 2.1070558 1982-1999 –> 2000-2012 0.7153662 3144.7763
SO (28.2, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1738 0.3884149 1661.042850 2.0901512 1982-1999 –> 2000-2012 0.7108438 3096.2954
SO (28.2, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1738 0.3884149 1661.042850 2.0901512 1982-1999 –> 2000-2012 0.7108438 3096.2954
SO (28.2, Inf] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 1738 0.3914217 1685.847473 2.1087884 1982-1999 –> 2000-2012 0.7160679 3152.9833
SO (28.2, Inf] 6 cstar_tref ~ sal + temp + aou + silicate + phosphate + phosphate_star 1005 0.4219558 1133.728559 2.6670990 2000-2012 –> 2013-2019 0.8103707 2794.7714
SO (28.2, Inf] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1005 0.4237274 1140.150063 2.6659440 2000-2012 –> 2013-2019 0.8147581 2822.5235
SO (28.2, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1005 0.4219558 1133.728559 2.6670990 2000-2012 –> 2013-2019 0.8103707 2794.7714
SO (28.2, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1005 0.4219558 1133.728559 2.6670990 2000-2012 –> 2013-2019 0.8103707 2794.7714
SO (28.2, Inf] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 1005 0.4240694 1141.771704 2.6663305 2000-2012 –> 2013-2019 0.8154911 2827.6192

5.3.2 per fitting unit

lm_best %>%
  select(-n_predictors) %>%
  group_by(basin, gamma_slab, eras) %>% 
  summarise(rmse_sum_mean = mean(rmse_sum),
            ais_sum_mean = mean(aic_sum)) %>% 
  ungroup() %>% 
  kable() %>%
  add_header_above() %>%
  kable_styling() %>%
  scroll_box(width = "100%", height = "400px")
basin gamma_slab eras rmse_sum_mean ais_sum_mean
Atlantic (-Inf,26] 1982-1999 –> 2000-2012 2.3860490 1199.4074
Atlantic (-Inf,26] 2000-2012 –> 2013-2019 2.2373386 721.7595
Atlantic (26,26.5] 1982-1999 –> 2000-2012 6.9568677 13619.8883
Atlantic (26,26.5] 2000-2012 –> 2013-2019 7.1353860 8642.2664
Atlantic (26.5,26.75] 1982-1999 –> 2000-2012 6.2106101 16303.4590
Atlantic (26.5,26.75] 2000-2012 –> 2013-2019 6.5421342 10461.0473
Atlantic (26.75,27] 1982-1999 –> 2000-2012 3.5355856 19047.4916
Atlantic (26.75,27] 2000-2012 –> 2013-2019 4.3913744 14165.4043
Atlantic (27,27.25] 1982-1999 –> 2000-2012 3.5824851 15728.6711
Atlantic (27,27.25] 2000-2012 –> 2013-2019 4.3559131 11733.8706
Atlantic (27.25,27.5] 1982-1999 –> 2000-2012 5.0689184 18034.7790
Atlantic (27.25,27.5] 2000-2012 –> 2013-2019 6.3999351 12850.3432
Atlantic (27.5,27.75] 1982-1999 –> 2000-2012 4.4150548 21489.4771
Atlantic (27.5,27.75] 2000-2012 –> 2013-2019 4.8761269 14674.6911
Atlantic (27.75,27.85] 1982-1999 –> 2000-2012 3.2583008 6447.0868
Atlantic (27.75,27.85] 2000-2012 –> 2013-2019 3.2648944 4159.6380
Atlantic (27.85,27.95] 1982-1999 –> 2000-2012 5.6401812 8184.5185
Atlantic (27.85,27.95] 2000-2012 –> 2013-2019 5.7077649 5466.3259
Atlantic (27.95,28.05] 1982-1999 –> 2000-2012 6.8143793 9144.2824
Atlantic (27.95,28.05] 2000-2012 –> 2013-2019 5.8795670 5596.5176
Atlantic (28.05,28.1] 1982-1999 –> 2000-2012 0.3150740 -1025.5854
Atlantic (28.05,28.1] 2000-2012 –> 2013-2019 0.4219826 -183.8906
Atlantic (28.1,28.15] 1982-1999 –> 2000-2012 0.2714782 -1494.7434
Atlantic (28.1,28.15] 2000-2012 –> 2013-2019 0.2631879 -960.0331
Atlantic (28.15,28.2] 1982-1999 –> 2000-2012 0.6720124 2027.8486
Atlantic (28.15,28.2] 2000-2012 –> 2013-2019 0.6490596 962.7797
Atlantic (28.2, Inf] 1982-1999 –> 2000-2012 0.8599221 2538.4815
Atlantic (28.2, Inf] 2000-2012 –> 2013-2019 0.8172039 1280.1155
Indo-Pacific (-Inf,26] 1982-1999 –> 2000-2012 12.6705792 81600.2300
Indo-Pacific (-Inf,26] 2000-2012 –> 2013-2019 12.0694628 52279.3445
Indo-Pacific (26,26.5] 1982-1999 –> 2000-2012 8.7769534 75354.9659
Indo-Pacific (26,26.5] 2000-2012 –> 2013-2019 9.0175060 48638.5968
Indo-Pacific (26.5,26.75] 1982-1999 –> 2000-2012 7.6613907 56279.8460
Indo-Pacific (26.5,26.75] 2000-2012 –> 2013-2019 8.1943572 37151.5752
Indo-Pacific (26.75,27] 1982-1999 –> 2000-2012 7.6731337 60382.3707
Indo-Pacific (26.75,27] 2000-2012 –> 2013-2019 8.7286782 40381.0754
Indo-Pacific (27,27.25] 1982-1999 –> 2000-2012 6.1949948 50798.0892
Indo-Pacific (27,27.25] 2000-2012 –> 2013-2019 7.3535533 35684.5904
Indo-Pacific (27.25,27.5] 1982-1999 –> 2000-2012 5.6601026 40510.0392
Indo-Pacific (27.25,27.5] 2000-2012 –> 2013-2019 6.1934891 27285.1004
Indo-Pacific (27.5,27.75] 1982-1999 –> 2000-2012 3.5311302 30322.6979
Indo-Pacific (27.5,27.75] 2000-2012 –> 2013-2019 4.0251272 20374.7973
Indo-Pacific (27.75,27.85] 1982-1999 –> 2000-2012 2.3338279 9154.9151
Indo-Pacific (27.75,27.85] 2000-2012 –> 2013-2019 2.4949399 6221.2425
Indo-Pacific (27.85,27.95] 1982-1999 –> 2000-2012 2.2032667 12986.4284
Indo-Pacific (27.85,27.95] 2000-2012 –> 2013-2019 2.5399320 9126.8325
Indo-Pacific (27.95,28.05] 1982-1999 –> 2000-2012 1.1066920 4780.1916
Indo-Pacific (27.95,28.05] 2000-2012 –> 2013-2019 1.1650727 3197.0845
Indo-Pacific (28.05,28.1] 1982-1999 –> 2000-2012 0.9086891 3223.2668
Indo-Pacific (28.05,28.1] 2000-2012 –> 2013-2019 0.9293922 2071.8335
Indo-Pacific (28.1, Inf] 1982-1999 –> 2000-2012 1.5367619 51486.0630
Indo-Pacific (28.1, Inf] 2000-2012 –> 2013-2019 1.5434237 33484.7817
SO (26,26.5] 1982-1999 –> 2000-2012 6.1453322 1212.5222
SO (26,26.5] 2000-2012 –> 2013-2019 5.8029715 766.0697
SO (26.5,26.75] 1982-1999 –> 2000-2012 6.6459824 3874.4170
SO (26.5,26.75] 2000-2012 –> 2013-2019 6.4841397 2331.0601
SO (26.75,27] 1982-1999 –> 2000-2012 4.7217642 11787.9840
SO (26.75,27] 2000-2012 –> 2013-2019 4.5617090 7267.5448
SO (27,27.25] 1982-1999 –> 2000-2012 3.2026325 22483.2958
SO (27,27.25] 2000-2012 –> 2013-2019 3.4251335 14773.3796
SO (27.25,27.5] 1982-1999 –> 2000-2012 3.1250360 24192.3064
SO (27.25,27.5] 2000-2012 –> 2013-2019 3.5342275 16546.8304
SO (27.5,27.75] 1982-1999 –> 2000-2012 2.8137866 27702.2609
SO (27.5,27.75] 2000-2012 –> 2013-2019 3.2875073 19761.4834
SO (27.75,27.85] 1982-1999 –> 2000-2012 2.7432935 11135.8964
SO (27.75,27.85] 2000-2012 –> 2013-2019 3.0526575 7842.3035
SO (27.85,27.95] 1982-1999 –> 2000-2012 2.5899132 10923.5081
SO (27.85,27.95] 2000-2012 –> 2013-2019 2.8164355 7352.8135
SO (27.95,28.05] 1982-1999 –> 2000-2012 2.8571479 15050.0982
SO (27.95,28.05] 2000-2012 –> 2013-2019 3.1939498 10445.4546
SO (28.05,28.1] 1982-1999 –> 2000-2012 2.0064080 7337.0822
SO (28.05,28.1] 2000-2012 –> 2013-2019 2.2764948 5138.7691
SO (28.1, Inf] 1982-1999 –> 2000-2012 2.0803628 36692.7235
SO (28.1, Inf] 2000-2012 –> 2013-2019 2.2933541 24768.2091
SO (28.1,28.15] 1982-1999 –> 2000-2012 0.8745097 1104.2473
SO (28.1,28.15] 2000-2012 –> 2013-2019 0.9244657 803.5287
SO (28.15,28.2] 1982-1999 –> 2000-2012 0.5940817 392.4480
SO (28.15,28.2] 2000-2012 –> 2013-2019 0.6175859 284.9479
SO (28.2, Inf] 1982-1999 –> 2000-2012 0.7127931 3117.3292
SO (28.2, Inf] 2000-2012 –> 2013-2019 0.8122723 2806.8914

5.4 Target variable coefficients

A data frame to map the target variable is prepared.

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

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

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

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

lm_best_target <- left_join(lm_best_target, lm_all_fitted_wide)

rm(eras_era, eras_forward, eras_backward,
   lm_all_fitted)

5.5 Plot selected model residuals

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

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

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

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

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

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

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

5.6 Cant coeffcients

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

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

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

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

5.7 Write files

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

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

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

6 Model diagnotics

6.1 Selection criterion vs predictors

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

6.1.1 All models

lm_all_fitted_wide %>%
  ggplot(aes(as.factor(n_predictors),
             !!sym(params_local$MLR_criterion),
             col = basin)) +
  geom_hline(yintercept = 10) +
  geom_boxplot() +
  facet_grid(gamma_slab~era) +
  scale_color_brewer(palette = "Set1") +
  labs(x="Number of predictors")

Version Author Date
4168b43 Donghe-Zhu 2021-03-10
5365f80 Donghe-Zhu 2021-03-10
2b6c392 Donghe-Zhu 2021-03-10
9f58753 Donghe-Zhu 2021-03-10
9d5a62c Donghe-Zhu 2021-03-10
5d1e70b Donghe-Zhu 2021-03-10
2093979 Donghe-Zhu 2021-03-10
b865899 Donghe-Zhu 2021-03-10
cc2a956 Donghe-Zhu 2021-03-10
60689fb Donghe-Zhu 2021-03-10
dba33c8 Donghe-Zhu 2021-03-09
17f1c4a Donghe-Zhu 2021-03-09
c024d1a Donghe-Zhu 2021-03-09
02f7242 Donghe-Zhu 2021-03-09
6f50bc6 Donghe-Zhu 2021-03-09
1691156 Donghe-Zhu 2021-03-08
c0ceaf8 Donghe-Zhu 2021-03-08
058e0a1 Donghe-Zhu 2021-03-08
112dea0 Donghe-Zhu 2021-03-08
1843412 Donghe-Zhu 2021-03-08
3fbbfa4 Donghe-Zhu 2021-03-07
627c8fb Donghe-Zhu 2021-03-07
9ef3222 Donghe-Zhu 2021-03-05
8c1e978 Donghe-Zhu 2021-03-05
865f68c Donghe-Zhu 2021-03-05
ee69bc1 Donghe-Zhu 2021-03-05
a79291f Donghe-Zhu 2021-03-05
e8c6f30 Donghe-Zhu 2021-03-04
59288fe Donghe-Zhu 2021-03-04
731abc8 Donghe-Zhu 2021-03-04
e2a5a33 Donghe-Zhu 2021-03-04
c7892c1 Donghe-Zhu 2021-03-04
924430b Donghe-Zhu 2021-03-03
0d0bca1 Donghe-Zhu 2021-03-03
cb63c16 Donghe-Zhu 2021-03-03
ffda45a Donghe-Zhu 2021-03-03
691ba81 Donghe-Zhu 2021-03-03
c5e45a2 Donghe-Zhu 2021-03-03
89c3e58 Donghe-Zhu 2021-03-03
c911669 Donghe-Zhu 2021-03-03
b71c719 Donghe-Zhu 2021-03-01
13666ca Donghe-Zhu 2021-03-01
c6e60fe Donghe-Zhu 2021-03-01
7a388f7 Donghe-Zhu 2021-03-01
799e913 Donghe-Zhu 2021-03-01
66ff99f Donghe-Zhu 2021-03-01
ac9bb7a Donghe-Zhu 2021-02-28
efdc047 Donghe-Zhu 2021-02-28
e9a7418 Donghe-Zhu 2021-02-28
e152917 Donghe-Zhu 2021-02-28
287123c Donghe-Zhu 2021-02-27
54d5b5b Donghe-Zhu 2021-02-27
330f064 Donghe-Zhu 2021-02-27
adbc9bc Donghe-Zhu 2021-02-27
5937141 Donghe-Zhu 2021-02-27
4414bbf Donghe-Zhu 2021-02-27
a265efb Donghe-Zhu 2021-02-27
19edd1e Donghe-Zhu 2021-02-27
f20483f Donghe-Zhu 2021-02-26
6a2c7b3 Donghe-Zhu 2021-02-25
02b976d Donghe-Zhu 2021-02-24
354c224 Donghe-Zhu 2021-02-24
1a0a88a Donghe-Zhu 2021-02-24
57f701e Donghe-Zhu 2021-02-24
06f3149 Donghe-Zhu 2021-02-16
401eab3 Donghe-Zhu 2021-02-15
e3bba84 Donghe-Zhu 2021-02-15
5dce4b1 Donghe-Zhu 2021-02-15
4469a0c Donghe-Zhu 2021-02-13
5ae6a69 Donghe-Zhu 2021-02-10
05385dc Donghe-Zhu 2021-02-10
f791ae4 Donghe-Zhu 2021-02-09
f71ae34 Donghe-Zhu 2021-02-09
c011832 Donghe-Zhu 2021-02-09
a145fa7 Donghe-Zhu 2021-02-09
c344e42 Donghe-Zhu 2021-02-08
2f095d7 Donghe-Zhu 2021-02-07
1fad5f1 Donghe-Zhu 2021-02-07
ca03c39 Donghe-Zhu 2021-02-07
cd7c52c Donghe-Zhu 2021-02-04
bcf84f4 Donghe-Zhu 2021-02-02
a518739 Donghe-Zhu 2021-02-01
61666de Donghe-Zhu 2021-01-31
865b582 Donghe-Zhu 2021-01-31
3e68089 Donghe-Zhu 2021-01-31
ecf335c Donghe-Zhu 2021-01-31
a618965 Donghe-Zhu 2021-01-31
59e006e Donghe-Zhu 2021-01-31
a1c8f87 Donghe-Zhu 2021-01-31
ae5c18f Donghe-Zhu 2021-01-31
b50fe52 Donghe-Zhu 2021-01-31
ac99ae5 jens-daniel-mueller 2021-01-29
b5bdcaf Donghe-Zhu 2021-01-29
442010d Donghe-Zhu 2021-01-29
372adf5 Donghe-Zhu 2021-01-29
af8788e Donghe-Zhu 2021-01-29
21c91c9 Donghe-Zhu 2021-01-29
eded038 Donghe-Zhu 2021-01-29
541d4dd Donghe-Zhu 2021-01-29
6a75576 Donghe-Zhu 2021-01-28
16fba40 Donghe-Zhu 2021-01-28
12bc567 Donghe-Zhu 2021-01-27
ceed31b Donghe-Zhu 2021-01-27
342402d Donghe-Zhu 2021-01-27
5bad5c2 Donghe-Zhu 2021-01-27
61efb56 Donghe-Zhu 2021-01-25
48f638e Donghe-Zhu 2021-01-25
c1cec47 Donghe-Zhu 2021-01-25
05ffb0c Donghe-Zhu 2021-01-25
8b97165 Donghe-Zhu 2021-01-25
c569946 Donghe-Zhu 2021-01-24
a2f0d56 Donghe-Zhu 2021-01-23
28509fc Donghe-Zhu 2021-01-23
4c28e4a Donghe-Zhu 2021-01-22
24cc264 jens-daniel-mueller 2021-01-22
7891955 Donghe-Zhu 2021-01-21
d4cf1cb Donghe-Zhu 2021-01-21
1f3e5b6 jens-daniel-mueller 2021-01-20
0e7bdf1 jens-daniel-mueller 2021-01-15
4571843 jens-daniel-mueller 2021-01-14
b3564aa jens-daniel-mueller 2021-01-14
8d032c3 jens-daniel-mueller 2021-01-14
17dee1d jens-daniel-mueller 2021-01-13
7cdea0c jens-daniel-mueller 2021-01-06
fa85b93 jens-daniel-mueller 2021-01-06
e5cb81a Donghe-Zhu 2021-01-05
a499f10 Donghe-Zhu 2021-01-05
fb8a752 Donghe-Zhu 2020-12-23
8fae0b2 Donghe-Zhu 2020-12-21
c8b76b3 jens-daniel-mueller 2020-12-19

6.1.2 Best models

lm_best_target %>%
  ggplot(aes("",
             !!sym(params_local$MLR_criterion),
             col = basin)) +
  geom_hline(yintercept = 10) +
  geom_boxplot() +
  facet_grid(gamma_slab~era) +
  scale_color_brewer(palette = "Set1") +
  labs(x="Number of predictors pooled")

Version Author Date
4168b43 Donghe-Zhu 2021-03-10
5365f80 Donghe-Zhu 2021-03-10
2b6c392 Donghe-Zhu 2021-03-10
9f58753 Donghe-Zhu 2021-03-10
fd528ed Donghe-Zhu 2021-03-10
9d5a62c Donghe-Zhu 2021-03-10
5d1e70b Donghe-Zhu 2021-03-10
2093979 Donghe-Zhu 2021-03-10
b865899 Donghe-Zhu 2021-03-10
cc2a956 Donghe-Zhu 2021-03-10
60689fb Donghe-Zhu 2021-03-10
dba33c8 Donghe-Zhu 2021-03-09
17f1c4a Donghe-Zhu 2021-03-09
c024d1a Donghe-Zhu 2021-03-09
02f7242 Donghe-Zhu 2021-03-09
6f50bc6 Donghe-Zhu 2021-03-09
1691156 Donghe-Zhu 2021-03-08
c0ceaf8 Donghe-Zhu 2021-03-08
058e0a1 Donghe-Zhu 2021-03-08
112dea0 Donghe-Zhu 2021-03-08
1843412 Donghe-Zhu 2021-03-08
3fbbfa4 Donghe-Zhu 2021-03-07
627c8fb Donghe-Zhu 2021-03-07
9ef3222 Donghe-Zhu 2021-03-05
8c1e978 Donghe-Zhu 2021-03-05
865f68c Donghe-Zhu 2021-03-05
ee69bc1 Donghe-Zhu 2021-03-05
a79291f Donghe-Zhu 2021-03-05
e8c6f30 Donghe-Zhu 2021-03-04
59288fe Donghe-Zhu 2021-03-04
731abc8 Donghe-Zhu 2021-03-04
e2a5a33 Donghe-Zhu 2021-03-04
c7892c1 Donghe-Zhu 2021-03-04
924430b Donghe-Zhu 2021-03-03
0d0bca1 Donghe-Zhu 2021-03-03
cb63c16 Donghe-Zhu 2021-03-03
ffda45a Donghe-Zhu 2021-03-03
691ba81 Donghe-Zhu 2021-03-03
c5e45a2 Donghe-Zhu 2021-03-03
89c3e58 Donghe-Zhu 2021-03-03
c407a50 Donghe-Zhu 2021-03-03
c911669 Donghe-Zhu 2021-03-03
b71c719 Donghe-Zhu 2021-03-01
13666ca Donghe-Zhu 2021-03-01
c6e60fe Donghe-Zhu 2021-03-01
7a388f7 Donghe-Zhu 2021-03-01
799e913 Donghe-Zhu 2021-03-01
66ff99f Donghe-Zhu 2021-03-01
ac9bb7a Donghe-Zhu 2021-02-28
efdc047 Donghe-Zhu 2021-02-28
e9a7418 Donghe-Zhu 2021-02-28
e152917 Donghe-Zhu 2021-02-28
287123c Donghe-Zhu 2021-02-27
54d5b5b Donghe-Zhu 2021-02-27
330f064 Donghe-Zhu 2021-02-27
adbc9bc Donghe-Zhu 2021-02-27
5937141 Donghe-Zhu 2021-02-27
4414bbf Donghe-Zhu 2021-02-27
a265efb Donghe-Zhu 2021-02-27
19edd1e Donghe-Zhu 2021-02-27
f20483f Donghe-Zhu 2021-02-26
6a2c7b3 Donghe-Zhu 2021-02-25
354c224 Donghe-Zhu 2021-02-24
1a0a88a Donghe-Zhu 2021-02-24
57f701e Donghe-Zhu 2021-02-24
06f3149 Donghe-Zhu 2021-02-16
401eab3 Donghe-Zhu 2021-02-15
e3bba84 Donghe-Zhu 2021-02-15
5dce4b1 Donghe-Zhu 2021-02-15
4469a0c Donghe-Zhu 2021-02-13
5ae6a69 Donghe-Zhu 2021-02-10
05385dc Donghe-Zhu 2021-02-10
f791ae4 Donghe-Zhu 2021-02-09
f71ae34 Donghe-Zhu 2021-02-09
c011832 Donghe-Zhu 2021-02-09
a145fa7 Donghe-Zhu 2021-02-09
c344e42 Donghe-Zhu 2021-02-08
2f095d7 Donghe-Zhu 2021-02-07
1fad5f1 Donghe-Zhu 2021-02-07
ca03c39 Donghe-Zhu 2021-02-07
e2ffc14 Donghe-Zhu 2021-02-05
cd7c52c Donghe-Zhu 2021-02-04
bcf84f4 Donghe-Zhu 2021-02-02
a518739 Donghe-Zhu 2021-02-01
61666de Donghe-Zhu 2021-01-31
865b582 Donghe-Zhu 2021-01-31
3e68089 Donghe-Zhu 2021-01-31
ecf335c Donghe-Zhu 2021-01-31
a618965 Donghe-Zhu 2021-01-31
59e006e Donghe-Zhu 2021-01-31
a1c8f87 Donghe-Zhu 2021-01-31
ae5c18f Donghe-Zhu 2021-01-31
b50fe52 Donghe-Zhu 2021-01-31
ac99ae5 jens-daniel-mueller 2021-01-29
b5bdcaf Donghe-Zhu 2021-01-29
442010d Donghe-Zhu 2021-01-29
372adf5 Donghe-Zhu 2021-01-29
af8788e Donghe-Zhu 2021-01-29
21c91c9 Donghe-Zhu 2021-01-29
eded038 Donghe-Zhu 2021-01-29
541d4dd Donghe-Zhu 2021-01-29
6a75576 Donghe-Zhu 2021-01-28
16fba40 Donghe-Zhu 2021-01-28
12bc567 Donghe-Zhu 2021-01-27
ceed31b Donghe-Zhu 2021-01-27
342402d Donghe-Zhu 2021-01-27
5bad5c2 Donghe-Zhu 2021-01-27
61efb56 Donghe-Zhu 2021-01-25
48f638e Donghe-Zhu 2021-01-25
c1cec47 Donghe-Zhu 2021-01-25
05ffb0c Donghe-Zhu 2021-01-25
8b97165 Donghe-Zhu 2021-01-25
c569946 Donghe-Zhu 2021-01-24
a2f0d56 Donghe-Zhu 2021-01-23
28509fc Donghe-Zhu 2021-01-23
4c28e4a Donghe-Zhu 2021-01-22
24cc264 jens-daniel-mueller 2021-01-22
7891955 Donghe-Zhu 2021-01-21
d4cf1cb Donghe-Zhu 2021-01-21
1f3e5b6 jens-daniel-mueller 2021-01-20
0e7bdf1 jens-daniel-mueller 2021-01-15
4571843 jens-daniel-mueller 2021-01-14
b3564aa jens-daniel-mueller 2021-01-14
8d032c3 jens-daniel-mueller 2021-01-14
17dee1d jens-daniel-mueller 2021-01-13
7cdea0c jens-daniel-mueller 2021-01-06
fa85b93 jens-daniel-mueller 2021-01-06
e5cb81a Donghe-Zhu 2021-01-05
a499f10 Donghe-Zhu 2021-01-05
fb8a752 Donghe-Zhu 2020-12-23
8fae0b2 Donghe-Zhu 2020-12-21
c8b76b3 jens-daniel-mueller 2020-12-19

6.2 RMSE correlation between eras

RMSE was plotted to compare the agreement for one model applied to two adjacent eras (ie check whether the same predictor combination performs equal in both eras).

6.2.1 All models

# find max rmse to scale axis
max_rmse <-
  max(c(lm_all_fitted_wide_eras$rmse,
        lm_all_fitted_wide_eras$rmse_sum - lm_all_fitted_wide_eras$rmse))

lm_all_fitted_wide_eras %>%
  ggplot(aes(rmse, rmse_sum - rmse, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  coord_equal(xlim = c(0,max_rmse),
              ylim = c(0,max_rmse)) +
  facet_grid(eras ~ basin)

Version Author Date
5365f80 Donghe-Zhu 2021-03-10
2b6c392 Donghe-Zhu 2021-03-10
9f58753 Donghe-Zhu 2021-03-10
9d5a62c Donghe-Zhu 2021-03-10
5d1e70b Donghe-Zhu 2021-03-10
2093979 Donghe-Zhu 2021-03-10
b865899 Donghe-Zhu 2021-03-10
cc2a956 Donghe-Zhu 2021-03-10
dba33c8 Donghe-Zhu 2021-03-09
17f1c4a Donghe-Zhu 2021-03-09
c024d1a Donghe-Zhu 2021-03-09
02f7242 Donghe-Zhu 2021-03-09
1691156 Donghe-Zhu 2021-03-08
c0ceaf8 Donghe-Zhu 2021-03-08
058e0a1 Donghe-Zhu 2021-03-08
112dea0 Donghe-Zhu 2021-03-08
1843412 Donghe-Zhu 2021-03-08
3fbbfa4 Donghe-Zhu 2021-03-07
627c8fb Donghe-Zhu 2021-03-07
9ef3222 Donghe-Zhu 2021-03-05
8c1e978 Donghe-Zhu 2021-03-05
865f68c Donghe-Zhu 2021-03-05
ee69bc1 Donghe-Zhu 2021-03-05
a79291f Donghe-Zhu 2021-03-05
e8c6f30 Donghe-Zhu 2021-03-04
59288fe Donghe-Zhu 2021-03-04
731abc8 Donghe-Zhu 2021-03-04
e2a5a33 Donghe-Zhu 2021-03-04
c7892c1 Donghe-Zhu 2021-03-04
924430b Donghe-Zhu 2021-03-03
0d0bca1 Donghe-Zhu 2021-03-03
cb63c16 Donghe-Zhu 2021-03-03
ffda45a Donghe-Zhu 2021-03-03
691ba81 Donghe-Zhu 2021-03-03
c5e45a2 Donghe-Zhu 2021-03-03
89c3e58 Donghe-Zhu 2021-03-03
c911669 Donghe-Zhu 2021-03-03
13666ca Donghe-Zhu 2021-03-01
c6e60fe Donghe-Zhu 2021-03-01
7a388f7 Donghe-Zhu 2021-03-01
799e913 Donghe-Zhu 2021-03-01
66ff99f Donghe-Zhu 2021-03-01
ac9bb7a Donghe-Zhu 2021-02-28
efdc047 Donghe-Zhu 2021-02-28
e9a7418 Donghe-Zhu 2021-02-28
287123c Donghe-Zhu 2021-02-27
54d5b5b Donghe-Zhu 2021-02-27
330f064 Donghe-Zhu 2021-02-27
adbc9bc Donghe-Zhu 2021-02-27
5937141 Donghe-Zhu 2021-02-27
4414bbf Donghe-Zhu 2021-02-27
a265efb Donghe-Zhu 2021-02-27
19edd1e Donghe-Zhu 2021-02-27
f20483f Donghe-Zhu 2021-02-26
6a2c7b3 Donghe-Zhu 2021-02-25
02b976d Donghe-Zhu 2021-02-24
354c224 Donghe-Zhu 2021-02-24
1a0a88a Donghe-Zhu 2021-02-24
57f701e Donghe-Zhu 2021-02-24
06f3149 Donghe-Zhu 2021-02-16
e3bba84 Donghe-Zhu 2021-02-15
5dce4b1 Donghe-Zhu 2021-02-15
4469a0c Donghe-Zhu 2021-02-13
5ae6a69 Donghe-Zhu 2021-02-10
05385dc Donghe-Zhu 2021-02-10
f791ae4 Donghe-Zhu 2021-02-09
f71ae34 Donghe-Zhu 2021-02-09
c011832 Donghe-Zhu 2021-02-09
a145fa7 Donghe-Zhu 2021-02-09
c344e42 Donghe-Zhu 2021-02-08
2f095d7 Donghe-Zhu 2021-02-07
1fad5f1 Donghe-Zhu 2021-02-07
ca03c39 Donghe-Zhu 2021-02-07
cd7c52c Donghe-Zhu 2021-02-04
bcf84f4 Donghe-Zhu 2021-02-02
61666de Donghe-Zhu 2021-01-31
865b582 Donghe-Zhu 2021-01-31
3e68089 Donghe-Zhu 2021-01-31
ecf335c Donghe-Zhu 2021-01-31
a618965 Donghe-Zhu 2021-01-31
59e006e Donghe-Zhu 2021-01-31
a1c8f87 Donghe-Zhu 2021-01-31
ae5c18f Donghe-Zhu 2021-01-31
b50fe52 Donghe-Zhu 2021-01-31
ac99ae5 jens-daniel-mueller 2021-01-29
b5bdcaf Donghe-Zhu 2021-01-29
442010d Donghe-Zhu 2021-01-29
372adf5 Donghe-Zhu 2021-01-29
af8788e Donghe-Zhu 2021-01-29
21c91c9 Donghe-Zhu 2021-01-29
eded038 Donghe-Zhu 2021-01-29
541d4dd Donghe-Zhu 2021-01-29
6a75576 Donghe-Zhu 2021-01-28
16fba40 Donghe-Zhu 2021-01-28
12bc567 Donghe-Zhu 2021-01-27
ceed31b Donghe-Zhu 2021-01-27
342402d Donghe-Zhu 2021-01-27
5bad5c2 Donghe-Zhu 2021-01-27
61efb56 Donghe-Zhu 2021-01-25
48f638e Donghe-Zhu 2021-01-25
c1cec47 Donghe-Zhu 2021-01-25
05ffb0c Donghe-Zhu 2021-01-25
8b97165 Donghe-Zhu 2021-01-25
c569946 Donghe-Zhu 2021-01-24
a2f0d56 Donghe-Zhu 2021-01-23
28509fc Donghe-Zhu 2021-01-23
4c28e4a Donghe-Zhu 2021-01-22
24cc264 jens-daniel-mueller 2021-01-22
7891955 Donghe-Zhu 2021-01-21
d4cf1cb Donghe-Zhu 2021-01-21
1f3e5b6 jens-daniel-mueller 2021-01-20
0e7bdf1 jens-daniel-mueller 2021-01-15
4571843 jens-daniel-mueller 2021-01-14
b3564aa jens-daniel-mueller 2021-01-14
8d032c3 jens-daniel-mueller 2021-01-14
17dee1d jens-daniel-mueller 2021-01-13
7cdea0c jens-daniel-mueller 2021-01-06
fa85b93 jens-daniel-mueller 2021-01-06
e5cb81a Donghe-Zhu 2021-01-05
a499f10 Donghe-Zhu 2021-01-05
fb8a752 Donghe-Zhu 2020-12-23
8fae0b2 Donghe-Zhu 2020-12-21
c8b76b3 jens-daniel-mueller 2020-12-19 
rm(max_rmse)

6.2.2 Best models

# find max rmse to scale axis
max_rmse <-
  max(c(lm_best$rmse,
        lm_best$rmse_sum - lm_best$rmse))

lm_best %>%
  ggplot(aes(rmse, rmse_sum - rmse, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  coord_equal(xlim = c(0,max_rmse),
              ylim = c(0,max_rmse)) +
  facet_grid(eras ~ basin)

Version Author Date
4168b43 Donghe-Zhu 2021-03-10
5365f80 Donghe-Zhu 2021-03-10
2b6c392 Donghe-Zhu 2021-03-10
9f58753 Donghe-Zhu 2021-03-10
fd528ed Donghe-Zhu 2021-03-10
9d5a62c Donghe-Zhu 2021-03-10
5d1e70b Donghe-Zhu 2021-03-10
2093979 Donghe-Zhu 2021-03-10
b865899 Donghe-Zhu 2021-03-10
cc2a956 Donghe-Zhu 2021-03-10
60689fb Donghe-Zhu 2021-03-10
dba33c8 Donghe-Zhu 2021-03-09
17f1c4a Donghe-Zhu 2021-03-09
c024d1a Donghe-Zhu 2021-03-09
02f7242 Donghe-Zhu 2021-03-09
6f50bc6 Donghe-Zhu 2021-03-09
1691156 Donghe-Zhu 2021-03-08
c0ceaf8 Donghe-Zhu 2021-03-08
058e0a1 Donghe-Zhu 2021-03-08
112dea0 Donghe-Zhu 2021-03-08
1843412 Donghe-Zhu 2021-03-08
3fbbfa4 Donghe-Zhu 2021-03-07
627c8fb Donghe-Zhu 2021-03-07
9ef3222 Donghe-Zhu 2021-03-05
8c1e978 Donghe-Zhu 2021-03-05
865f68c Donghe-Zhu 2021-03-05
ee69bc1 Donghe-Zhu 2021-03-05
a79291f Donghe-Zhu 2021-03-05
e8c6f30 Donghe-Zhu 2021-03-04
59288fe Donghe-Zhu 2021-03-04
731abc8 Donghe-Zhu 2021-03-04
e2a5a33 Donghe-Zhu 2021-03-04
c7892c1 Donghe-Zhu 2021-03-04
924430b Donghe-Zhu 2021-03-03
0d0bca1 Donghe-Zhu 2021-03-03
cb63c16 Donghe-Zhu 2021-03-03
ffda45a Donghe-Zhu 2021-03-03
691ba81 Donghe-Zhu 2021-03-03
c5e45a2 Donghe-Zhu 2021-03-03
89c3e58 Donghe-Zhu 2021-03-03
c407a50 Donghe-Zhu 2021-03-03
c911669 Donghe-Zhu 2021-03-03
b71c719 Donghe-Zhu 2021-03-01
13666ca Donghe-Zhu 2021-03-01
c6e60fe Donghe-Zhu 2021-03-01
7a388f7 Donghe-Zhu 2021-03-01
799e913 Donghe-Zhu 2021-03-01
66ff99f Donghe-Zhu 2021-03-01
ac9bb7a Donghe-Zhu 2021-02-28
efdc047 Donghe-Zhu 2021-02-28
e9a7418 Donghe-Zhu 2021-02-28
e152917 Donghe-Zhu 2021-02-28
287123c Donghe-Zhu 2021-02-27
54d5b5b Donghe-Zhu 2021-02-27
330f064 Donghe-Zhu 2021-02-27
adbc9bc Donghe-Zhu 2021-02-27
5937141 Donghe-Zhu 2021-02-27
4414bbf Donghe-Zhu 2021-02-27
a265efb Donghe-Zhu 2021-02-27
19edd1e Donghe-Zhu 2021-02-27
f20483f Donghe-Zhu 2021-02-26
6a2c7b3 Donghe-Zhu 2021-02-25
354c224 Donghe-Zhu 2021-02-24
1a0a88a Donghe-Zhu 2021-02-24
57f701e Donghe-Zhu 2021-02-24
06f3149 Donghe-Zhu 2021-02-16
401eab3 Donghe-Zhu 2021-02-15
e3bba84 Donghe-Zhu 2021-02-15
5dce4b1 Donghe-Zhu 2021-02-15
4469a0c Donghe-Zhu 2021-02-13
5ae6a69 Donghe-Zhu 2021-02-10
05385dc Donghe-Zhu 2021-02-10
f791ae4 Donghe-Zhu 2021-02-09
f71ae34 Donghe-Zhu 2021-02-09
c011832 Donghe-Zhu 2021-02-09
a145fa7 Donghe-Zhu 2021-02-09
c344e42 Donghe-Zhu 2021-02-08
2f095d7 Donghe-Zhu 2021-02-07
1fad5f1 Donghe-Zhu 2021-02-07
ca03c39 Donghe-Zhu 2021-02-07
e2ffc14 Donghe-Zhu 2021-02-05
cd7c52c Donghe-Zhu 2021-02-04
bcf84f4 Donghe-Zhu 2021-02-02
a518739 Donghe-Zhu 2021-02-01
61666de Donghe-Zhu 2021-01-31
865b582 Donghe-Zhu 2021-01-31
3e68089 Donghe-Zhu 2021-01-31
ecf335c Donghe-Zhu 2021-01-31
a618965 Donghe-Zhu 2021-01-31
59e006e Donghe-Zhu 2021-01-31
a1c8f87 Donghe-Zhu 2021-01-31
ae5c18f Donghe-Zhu 2021-01-31
b50fe52 Donghe-Zhu 2021-01-31
ac99ae5 jens-daniel-mueller 2021-01-29
b5bdcaf Donghe-Zhu 2021-01-29
442010d Donghe-Zhu 2021-01-29
372adf5 Donghe-Zhu 2021-01-29
af8788e Donghe-Zhu 2021-01-29
21c91c9 Donghe-Zhu 2021-01-29
eded038 Donghe-Zhu 2021-01-29
541d4dd Donghe-Zhu 2021-01-29
6a75576 Donghe-Zhu 2021-01-28
16fba40 Donghe-Zhu 2021-01-28
12bc567 Donghe-Zhu 2021-01-27
ceed31b Donghe-Zhu 2021-01-27
342402d Donghe-Zhu 2021-01-27
5bad5c2 Donghe-Zhu 2021-01-27
61efb56 Donghe-Zhu 2021-01-25
48f638e Donghe-Zhu 2021-01-25
c1cec47 Donghe-Zhu 2021-01-25
05ffb0c Donghe-Zhu 2021-01-25
8b97165 Donghe-Zhu 2021-01-25
c569946 Donghe-Zhu 2021-01-24
a2f0d56 Donghe-Zhu 2021-01-23
28509fc Donghe-Zhu 2021-01-23
4c28e4a Donghe-Zhu 2021-01-22
24cc264 jens-daniel-mueller 2021-01-22
7891955 Donghe-Zhu 2021-01-21
d4cf1cb Donghe-Zhu 2021-01-21
1f3e5b6 jens-daniel-mueller 2021-01-20
0e7bdf1 jens-daniel-mueller 2021-01-15
4571843 jens-daniel-mueller 2021-01-14
b3564aa jens-daniel-mueller 2021-01-14
8d032c3 jens-daniel-mueller 2021-01-14
17dee1d jens-daniel-mueller 2021-01-13
7cdea0c jens-daniel-mueller 2021-01-06
fa85b93 jens-daniel-mueller 2021-01-06
e5cb81a Donghe-Zhu 2021-01-05
a499f10 Donghe-Zhu 2021-01-05
fb8a752 Donghe-Zhu 2020-12-23
8fae0b2 Donghe-Zhu 2020-12-21
c8b76b3 jens-daniel-mueller 2020-12-19 
rm(max_rmse)

6.3 Predictor counts

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

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

# print table
lm_all_stats %>%
  gt(rowname_col = "gamma_slab",
     groupname_col = c("basin", "eras")) %>% 
  summary_rows(
    groups = TRUE,
    fns = list(total = "sum")
  )
aou oxygen phosphate phosphate_star sal temp silicate
Atlantic - 1982-1999 --> 2000-2012
(-Inf,26] 5 3 4 2 5 2 0
(26,26.5] 5 3 4 3 5 3 5
(26.5,26.75] 5 3 4 3 3 5 5
(26.75,27] 5 4 3 3 3 5 5
(27,27.25] 5 2 4 3 5 4 5
(27.25,27.5] 5 3 4 3 5 5 2
(27.5,27.75] 5 4 3 3 5 5 3
(27.75,27.85] 5 4 3 3 5 5 3
(27.85,27.95] 5 3 4 3 2 5 5
(27.95,28.05] 5 4 3 3 3 5 5
(28.05,28.1] 5 4 3 3 5 3 5
(28.1,28.15] 5 3 4 3 5 3 5
(28.15,28.2] 5 3 4 3 5 3 5
(28.2, Inf] 5 4 3 3 5 3 5
total 70.00 47.00 50.00 41.00 61.00 56.00 58.00
Atlantic - 2000-2012 --> 2013-2019
(-Inf,26] 5 4 3 3 5 2 5
(26,26.5] 5 3 3 4 5 3 5
(26.5,26.75] 5 3 4 3 3 5 5
(26.75,27] 5 4 3 3 5 3 5
(27,27.25] 5 3 4 3 5 3 5
(27.25,27.5] 5 4 3 3 5 5 3
(27.5,27.75] 5 3 3 4 5 5 3
(27.75,27.85] 5 4 3 3 5 5 3
(27.85,27.95] 5 3 3 4 2 5 5
(27.95,28.05] 5 3 3 4 3 5 5
(28.05,28.1] 5 4 3 3 5 5 3
(28.1,28.15] 5 3 3 4 5 3 5
(28.15,28.2] 5 3 4 3 5 3 5
(28.2, Inf] 5 4 3 3 5 5 3
total 70.00 48.00 45.00 47.00 63.00 57.00 60.00
Indo-Pacific - 1982-1999 --> 2000-2012
(-Inf,26] 5 3 3 4 5 3 5
(26,26.5] 5 3 3 4 5 3 5
(26.5,26.75] 5 4 3 3 5 3 5
(26.75,27] 5 4 3 3 5 5 3
(27,27.25] 5 3 4 3 3 5 5
(27.25,27.5] 5 4 3 3 3 5 5
(27.5,27.75] 4 2 4 3 5 5 5
(27.75,27.85] 5 3 3 4 5 5 3
(27.85,27.95] 5 3 3 4 5 5 3
(27.95,28.05] 5 4 3 3 3 5 5
(28.05,28.1] 5 4 3 3 3 5 5
(28.1, Inf] 5 4 3 3 5 3 5
total 59.00 41.00 38.00 40.00 52.00 52.00 54.00
Indo-Pacific - 2000-2012 --> 2013-2019
(-Inf,26] 5 3 4 3 5 3 5
(26,26.5] 5 4 3 3 5 3 5
(26.5,26.75] 5 3 4 3 3 5 5
(26.75,27] 5 3 4 3 3 5 5
(27,27.25] 5 3 4 3 3 5 5
(27.25,27.5] 5 4 3 3 3 5 5
(27.5,27.75] 4 3 4 2 5 5 5
(27.75,27.85] 5 3 3 4 5 5 3
(27.85,27.95] 5 4 3 2 5 5 4
(27.95,28.05] 5 4 3 3 3 5 5
(28.05,28.1] 5 3 4 3 3 5 5
(28.1, Inf] 5 3 4 3 3 5 5
total 59.00 40.00 43.00 35.00 46.00 56.00 57.00
SO - 1982-1999 --> 2000-2012
(26,26.5] 3 3 3 4 5 5 5
(26.5,26.75] 5 3 4 3 5 5 3
(26.75,27] 5 3 3 4 3 5 5
(27,27.25] 5 4 3 3 3 5 5
(27.25,27.5] 5 3 3 3 4 5 5
(27.5,27.75] 5 2 4 3 5 4 5
(27.75,27.85] 5 3 3 4 5 3 5
(27.85,27.95] 5 3 3 3 5 5 4
(27.95,28.05] 5 3 4 3 5 5 3
(28.05,28.1] 5 4 3 3 5 5 2
(28.1, Inf] 5 3 4 3 5 3 5
(28.1,28.15] 5 4 3 3 5 5 3
(28.15,28.2] 5 4 3 3 5 5 3
(28.2, Inf] 4 2 3 4 5 5 5
total 67.00 44.00 46.00 46.00 65.00 65.00 58.00
SO - 2000-2012 --> 2013-2019
(26,26.5] 3 3 3 3 5 5 5
(26.5,26.75] 5 4 3 3 5 5 3
(26.75,27] 5 3 3 4 3 5 5
(27,27.25] 5 4 3 3 3 5 5
(27.25,27.5] 3 3 4 3 5 5 5
(27.5,27.75] 4 2 4 3 5 2 5
(27.75,27.85] 5 4 3 3 5 3 5
(27.85,27.95] 5 3 3 3 4 5 5
(27.95,28.05] 5 4 3 2 5 5 4
(28.05,28.1] 5 3 4 3 5 5 2
(28.1, Inf] 5 4 3 3 5 3 5
(28.1,28.15] 5 3 4 3 5 5 3
(28.15,28.2] 5 3 4 3 5 5 3
(28.2, Inf] 4 3 3 3 5 5 5
total 64.00 46.00 47.00 42.00 65.00 63.00 60.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)

Version Author Date
5365f80 Donghe-Zhu 2021-03-10
2b6c392 Donghe-Zhu 2021-03-10
9f58753 Donghe-Zhu 2021-03-10
9d5a62c Donghe-Zhu 2021-03-10
5d1e70b Donghe-Zhu 2021-03-10
2093979 Donghe-Zhu 2021-03-10
b865899 Donghe-Zhu 2021-03-10
cc2a956 Donghe-Zhu 2021-03-10
dba33c8 Donghe-Zhu 2021-03-09
17f1c4a Donghe-Zhu 2021-03-09
c024d1a Donghe-Zhu 2021-03-09
02f7242 Donghe-Zhu 2021-03-09
1691156 Donghe-Zhu 2021-03-08
c0ceaf8 Donghe-Zhu 2021-03-08
058e0a1 Donghe-Zhu 2021-03-08
112dea0 Donghe-Zhu 2021-03-08
1843412 Donghe-Zhu 2021-03-08
3fbbfa4 Donghe-Zhu 2021-03-07
627c8fb Donghe-Zhu 2021-03-07
9ef3222 Donghe-Zhu 2021-03-05
8c1e978 Donghe-Zhu 2021-03-05
865f68c Donghe-Zhu 2021-03-05
ee69bc1 Donghe-Zhu 2021-03-05
a79291f Donghe-Zhu 2021-03-05
e8c6f30 Donghe-Zhu 2021-03-04
59288fe Donghe-Zhu 2021-03-04
731abc8 Donghe-Zhu 2021-03-04
e2a5a33 Donghe-Zhu 2021-03-04
c7892c1 Donghe-Zhu 2021-03-04
924430b Donghe-Zhu 2021-03-03
0d0bca1 Donghe-Zhu 2021-03-03
cb63c16 Donghe-Zhu 2021-03-03
ffda45a Donghe-Zhu 2021-03-03
691ba81 Donghe-Zhu 2021-03-03
c5e45a2 Donghe-Zhu 2021-03-03
89c3e58 Donghe-Zhu 2021-03-03
c911669 Donghe-Zhu 2021-03-03
13666ca Donghe-Zhu 2021-03-01
c6e60fe Donghe-Zhu 2021-03-01
7a388f7 Donghe-Zhu 2021-03-01
799e913 Donghe-Zhu 2021-03-01
66ff99f Donghe-Zhu 2021-03-01
ac9bb7a Donghe-Zhu 2021-02-28
efdc047 Donghe-Zhu 2021-02-28
e9a7418 Donghe-Zhu 2021-02-28
287123c Donghe-Zhu 2021-02-27
54d5b5b Donghe-Zhu 2021-02-27
330f064 Donghe-Zhu 2021-02-27
adbc9bc Donghe-Zhu 2021-02-27
5937141 Donghe-Zhu 2021-02-27
4414bbf Donghe-Zhu 2021-02-27
a265efb Donghe-Zhu 2021-02-27
19edd1e Donghe-Zhu 2021-02-27
f20483f Donghe-Zhu 2021-02-26
6a2c7b3 Donghe-Zhu 2021-02-25
02b976d Donghe-Zhu 2021-02-24
354c224 Donghe-Zhu 2021-02-24
1a0a88a Donghe-Zhu 2021-02-24
57f701e Donghe-Zhu 2021-02-24
06f3149 Donghe-Zhu 2021-02-16
e3bba84 Donghe-Zhu 2021-02-15
5dce4b1 Donghe-Zhu 2021-02-15
4469a0c Donghe-Zhu 2021-02-13
5ae6a69 Donghe-Zhu 2021-02-10
05385dc Donghe-Zhu 2021-02-10
f791ae4 Donghe-Zhu 2021-02-09
f71ae34 Donghe-Zhu 2021-02-09
c011832 Donghe-Zhu 2021-02-09
a145fa7 Donghe-Zhu 2021-02-09
c344e42 Donghe-Zhu 2021-02-08
2f095d7 Donghe-Zhu 2021-02-07
1fad5f1 Donghe-Zhu 2021-02-07
ca03c39 Donghe-Zhu 2021-02-07
cd7c52c Donghe-Zhu 2021-02-04
bcf84f4 Donghe-Zhu 2021-02-02
61666de Donghe-Zhu 2021-01-31
865b582 Donghe-Zhu 2021-01-31
3e68089 Donghe-Zhu 2021-01-31
ecf335c Donghe-Zhu 2021-01-31
a618965 Donghe-Zhu 2021-01-31
59e006e Donghe-Zhu 2021-01-31
a1c8f87 Donghe-Zhu 2021-01-31
ae5c18f Donghe-Zhu 2021-01-31
b50fe52 Donghe-Zhu 2021-01-31
ac99ae5 jens-daniel-mueller 2021-01-29
b5bdcaf Donghe-Zhu 2021-01-29
442010d Donghe-Zhu 2021-01-29
372adf5 Donghe-Zhu 2021-01-29
af8788e Donghe-Zhu 2021-01-29
21c91c9 Donghe-Zhu 2021-01-29
eded038 Donghe-Zhu 2021-01-29
541d4dd Donghe-Zhu 2021-01-29
6a75576 Donghe-Zhu 2021-01-28
16fba40 Donghe-Zhu 2021-01-28
12bc567 Donghe-Zhu 2021-01-27
ceed31b Donghe-Zhu 2021-01-27
342402d Donghe-Zhu 2021-01-27
5bad5c2 Donghe-Zhu 2021-01-27
61efb56 Donghe-Zhu 2021-01-25
48f638e Donghe-Zhu 2021-01-25
c1cec47 Donghe-Zhu 2021-01-25
05ffb0c Donghe-Zhu 2021-01-25
8b97165 Donghe-Zhu 2021-01-25
c569946 Donghe-Zhu 2021-01-24
a2f0d56 Donghe-Zhu 2021-01-23
28509fc Donghe-Zhu 2021-01-23
4c28e4a Donghe-Zhu 2021-01-22
24cc264 jens-daniel-mueller 2021-01-22
7891955 Donghe-Zhu 2021-01-21
d4cf1cb Donghe-Zhu 2021-01-21
1f3e5b6 jens-daniel-mueller 2021-01-20
0e7bdf1 jens-daniel-mueller 2021-01-15
4571843 jens-daniel-mueller 2021-01-14
b3564aa jens-daniel-mueller 2021-01-14
8d032c3 jens-daniel-mueller 2021-01-14
17dee1d jens-daniel-mueller 2021-01-13
7cdea0c jens-daniel-mueller 2021-01-06
fa85b93 jens-daniel-mueller 2021-01-06
e5cb81a Donghe-Zhu 2021-01-05
a499f10 Donghe-Zhu 2021-01-05
fb8a752 Donghe-Zhu 2020-12-23
8fae0b2 Donghe-Zhu 2020-12-21
c8b76b3 jens-daniel-mueller 2020-12-19 
lm_best %>% 
  ggplot(aes(rmse, aic, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  facet_grid(eras~basin)

Version Author Date
4168b43 Donghe-Zhu 2021-03-10
5365f80 Donghe-Zhu 2021-03-10
2b6c392 Donghe-Zhu 2021-03-10
9f58753 Donghe-Zhu 2021-03-10
fd528ed Donghe-Zhu 2021-03-10
9d5a62c Donghe-Zhu 2021-03-10
5d1e70b Donghe-Zhu 2021-03-10
2093979 Donghe-Zhu 2021-03-10
b865899 Donghe-Zhu 2021-03-10
cc2a956 Donghe-Zhu 2021-03-10
60689fb Donghe-Zhu 2021-03-10
dba33c8 Donghe-Zhu 2021-03-09
17f1c4a Donghe-Zhu 2021-03-09
c024d1a Donghe-Zhu 2021-03-09
02f7242 Donghe-Zhu 2021-03-09
6f50bc6 Donghe-Zhu 2021-03-09
1691156 Donghe-Zhu 2021-03-08
c0ceaf8 Donghe-Zhu 2021-03-08
058e0a1 Donghe-Zhu 2021-03-08
112dea0 Donghe-Zhu 2021-03-08
1843412 Donghe-Zhu 2021-03-08
3fbbfa4 Donghe-Zhu 2021-03-07
627c8fb Donghe-Zhu 2021-03-07
9ef3222 Donghe-Zhu 2021-03-05
8c1e978 Donghe-Zhu 2021-03-05
865f68c Donghe-Zhu 2021-03-05
ee69bc1 Donghe-Zhu 2021-03-05
a79291f Donghe-Zhu 2021-03-05
e8c6f30 Donghe-Zhu 2021-03-04
59288fe Donghe-Zhu 2021-03-04
731abc8 Donghe-Zhu 2021-03-04
e2a5a33 Donghe-Zhu 2021-03-04
c7892c1 Donghe-Zhu 2021-03-04
924430b Donghe-Zhu 2021-03-03
0d0bca1 Donghe-Zhu 2021-03-03
cb63c16 Donghe-Zhu 2021-03-03
ffda45a Donghe-Zhu 2021-03-03
691ba81 Donghe-Zhu 2021-03-03
c5e45a2 Donghe-Zhu 2021-03-03
89c3e58 Donghe-Zhu 2021-03-03
c407a50 Donghe-Zhu 2021-03-03
c911669 Donghe-Zhu 2021-03-03
b71c719 Donghe-Zhu 2021-03-01
13666ca Donghe-Zhu 2021-03-01
c6e60fe Donghe-Zhu 2021-03-01
7a388f7 Donghe-Zhu 2021-03-01
799e913 Donghe-Zhu 2021-03-01
66ff99f Donghe-Zhu 2021-03-01
ac9bb7a Donghe-Zhu 2021-02-28
efdc047 Donghe-Zhu 2021-02-28
e9a7418 Donghe-Zhu 2021-02-28
e152917 Donghe-Zhu 2021-02-28
287123c Donghe-Zhu 2021-02-27
54d5b5b Donghe-Zhu 2021-02-27
330f064 Donghe-Zhu 2021-02-27
adbc9bc Donghe-Zhu 2021-02-27
5937141 Donghe-Zhu 2021-02-27
4414bbf Donghe-Zhu 2021-02-27
a265efb Donghe-Zhu 2021-02-27
19edd1e Donghe-Zhu 2021-02-27
f20483f Donghe-Zhu 2021-02-26
6a2c7b3 Donghe-Zhu 2021-02-25
354c224 Donghe-Zhu 2021-02-24
1a0a88a Donghe-Zhu 2021-02-24
57f701e Donghe-Zhu 2021-02-24
06f3149 Donghe-Zhu 2021-02-16
401eab3 Donghe-Zhu 2021-02-15
e3bba84 Donghe-Zhu 2021-02-15
5dce4b1 Donghe-Zhu 2021-02-15
4469a0c Donghe-Zhu 2021-02-13
5ae6a69 Donghe-Zhu 2021-02-10
05385dc Donghe-Zhu 2021-02-10
f791ae4 Donghe-Zhu 2021-02-09
f71ae34 Donghe-Zhu 2021-02-09
c011832 Donghe-Zhu 2021-02-09
a145fa7 Donghe-Zhu 2021-02-09
c344e42 Donghe-Zhu 2021-02-08
2f095d7 Donghe-Zhu 2021-02-07
1fad5f1 Donghe-Zhu 2021-02-07
ca03c39 Donghe-Zhu 2021-02-07
e2ffc14 Donghe-Zhu 2021-02-05
cd7c52c Donghe-Zhu 2021-02-04
bcf84f4 Donghe-Zhu 2021-02-02
a518739 Donghe-Zhu 2021-02-01
61666de Donghe-Zhu 2021-01-31
865b582 Donghe-Zhu 2021-01-31
3e68089 Donghe-Zhu 2021-01-31
ecf335c Donghe-Zhu 2021-01-31
a618965 Donghe-Zhu 2021-01-31
59e006e Donghe-Zhu 2021-01-31
a1c8f87 Donghe-Zhu 2021-01-31
ae5c18f Donghe-Zhu 2021-01-31
b50fe52 Donghe-Zhu 2021-01-31
ac99ae5 jens-daniel-mueller 2021-01-29
b5bdcaf Donghe-Zhu 2021-01-29
442010d Donghe-Zhu 2021-01-29
372adf5 Donghe-Zhu 2021-01-29
af8788e Donghe-Zhu 2021-01-29
21c91c9 Donghe-Zhu 2021-01-29
eded038 Donghe-Zhu 2021-01-29
541d4dd Donghe-Zhu 2021-01-29
6a75576 Donghe-Zhu 2021-01-28
16fba40 Donghe-Zhu 2021-01-28
12bc567 Donghe-Zhu 2021-01-27
ceed31b Donghe-Zhu 2021-01-27
342402d Donghe-Zhu 2021-01-27
5bad5c2 Donghe-Zhu 2021-01-27
61efb56 Donghe-Zhu 2021-01-25
48f638e Donghe-Zhu 2021-01-25
c1cec47 Donghe-Zhu 2021-01-25
05ffb0c Donghe-Zhu 2021-01-25
8b97165 Donghe-Zhu 2021-01-25
c569946 Donghe-Zhu 2021-01-24
a2f0d56 Donghe-Zhu 2021-01-23
28509fc Donghe-Zhu 2021-01-23
4c28e4a Donghe-Zhu 2021-01-22
24cc264 jens-daniel-mueller 2021-01-22
7891955 Donghe-Zhu 2021-01-21
d4cf1cb Donghe-Zhu 2021-01-21
1f3e5b6 jens-daniel-mueller 2021-01-20
0e7bdf1 jens-daniel-mueller 2021-01-15
4571843 jens-daniel-mueller 2021-01-14
b3564aa jens-daniel-mueller 2021-01-14
8d032c3 jens-daniel-mueller 2021-01-14
17dee1d jens-daniel-mueller 2021-01-13
7cdea0c jens-daniel-mueller 2021-01-06
fa85b93 jens-daniel-mueller 2021-01-06
e5cb81a Donghe-Zhu 2021-01-05
a499f10 Donghe-Zhu 2021-01-05
fb8a752 Donghe-Zhu 2020-12-23
8fae0b2 Donghe-Zhu 2020-12-21
c8b76b3 jens-daniel-mueller 2020-12-19 

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