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: abf383e

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

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


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

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

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

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

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

1 Required data

Required are:

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

2 Predictor combinations

Find all possible combinations of following considered predictor variables:

  • sal, temp, oxygen, aou, silicate, 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:

  • 57

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:

  • rmse

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

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

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

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

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

5.3 RMSE tables

5.3.1 per model

lm_best %>%
  kable() %>%
  add_header_above() %>%
  kable_styling() %>%
  scroll_box(width = "100%", height = "400px")
basin gamma_slab n_predictors model nr_obs rmse aic resid_max eras rmse_sum aic_sum
Atlantic (-Inf,26] 4 cstar_tref ~ sal + oxygen + aou + phosphate_star 158 1.0492038 475.562558 4.3578441 1982-1999 –> 2000-2012 2.3855522 1198.45637
Atlantic (-Inf,26] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 158 1.0480832 477.224868 4.2990206 1982-1999 –> 2000-2012 2.3837757 1201.91446
Atlantic (-Inf,26] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 158 1.0553928 479.421090 4.5407576 1982-1999 –> 2000-2012 2.3925975 1204.58138
Atlantic (-Inf,26] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 158 1.0447877 476.229718 3.9913380 1982-1999 –> 2000-2012 2.3809752 1201.07344
Atlantic (-Inf,26] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 158 1.0432378 477.760573 3.9034959 1982-1999 –> 2000-2012 2.3788630 1204.42921
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.04067
Atlantic (-Inf,26] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 72 1.2012446 244.730723 3.4250804 2000-2012 –> 2013-2019 2.2566374 724.15181
Atlantic (-Inf,26] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 72 1.2718853 252.959191 3.5478791 2000-2012 –> 2013-2019 2.3166730 729.18891
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.83769
Atlantic (-Inf,26] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 72 1.2024324 244.873045 3.4452953 2000-2012 –> 2013-2019 2.2595421 724.80780
Atlantic (26,26.5] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 1040 3.5443539 5597.332405 13.4662233 1982-1999 –> 2000-2012 6.9995324 13651.58188
Atlantic (26,26.5] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1040 3.6046520 5632.420615 14.1850927 1982-1999 –> 2000-2012 7.1565388 13770.14714
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.75927
Atlantic (26,26.5] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 1040 3.6059532 5633.171340 14.2183254 1982-1999 –> 2000-2012 7.1620346 13774.46693
Atlantic (26,26.5] 4 cstar_tref ~ sal + temp + silicate + phosphate_star 1040 3.6722318 5669.055330 16.3679712 1982-1999 –> 2000-2012 7.5024485 14032.91968
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.70102
Atlantic (26,26.5] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 562 3.6408932 3061.352359 10.4893496 2000-2012 –> 2013-2019 7.2455452 8693.77297
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.64326
Atlantic (26,26.5] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 562 3.6406124 3061.265647 10.4898621 2000-2012 –> 2013-2019 7.2465656 8694.43699
Atlantic (26,26.5] 4 cstar_tref ~ sal + temp + silicate + phosphate_star 562 3.6468894 3061.201945 10.4063646 2000-2012 –> 2013-2019 7.3191212 8730.25727
Atlantic (26.5,26.75] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 1360 3.2843587 7108.059025 10.5247223 1982-1999 –> 2000-2012 6.3355501 16433.70428
Atlantic (26.5,26.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1360 3.3125149 7131.277720 10.5523314 1982-1999 –> 2000-2012 6.4028241 16503.72606
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.27934
Atlantic (26.5,26.75] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 1360 3.3168229 7134.812862 10.5593163 1982-1999 –> 2000-2012 6.4138058 16515.18682
Atlantic (26.5,26.75] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 1360 3.2415884 7072.405398 10.0491820 1982-1999 –> 2000-2012 6.2239552 16314.22857
Atlantic (26.5,26.75] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 646 3.3363780 3403.981011 10.7144118 2000-2012 –> 2013-2019 6.6207366 10512.04004
Atlantic (26.5,26.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 646 3.3594523 3412.885718 10.8878591 2000-2012 –> 2013-2019 6.6719672 10544.16344
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.85271
Atlantic (26.5,26.75] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 646 3.3631670 3414.313520 10.9361612 2000-2012 –> 2013-2019 6.6799899 10549.12638
Atlantic (26.5,26.75] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 646 3.3285183 3400.933773 11.1260032 2000-2012 –> 2013-2019 6.5701066 10473.33917
Atlantic (26.75,27] 4 cstar_tref ~ oxygen + aou + silicate + phosphate_star 2134 2.1608919 9356.613416 15.2322922 1982-1999 –> 2000-2012 3.8655083 20016.34118
Atlantic (26.75,27] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 2134 2.0150942 9060.471809 14.1929877 1982-1999 –> 2000-2012 3.6685116 19555.87536
Atlantic (26.75,27] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2134 2.1507424 9338.519878 22.8281773 1982-1999 –> 2000-2012 3.9108945 20175.10339
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.73246
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.13041
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.70999
Atlantic (26.75,27] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1121 2.6436309 5374.827940 25.0893644 2000-2012 –> 2013-2019 4.7943733 14713.34782
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.20047
Atlantic (26.75,27] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 1121 2.6552629 5384.671078 25.1149369 2000-2012 –> 2013-2019 4.8356436 14781.60466
Atlantic (26.75,27] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 1121 2.3110503 5073.388314 22.0043753 2000-2012 –> 2013-2019 4.4668470 14421.92640
Atlantic (27,27.25] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 1717 2.0072484 7279.325305 9.7017253 1982-1999 –> 2000-2012 3.6183821 15798.22804
Atlantic (27,27.25] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1717 2.0359536 7328.086282 8.6521412 1982-1999 –> 2000-2012 3.6652699 15897.33254
Atlantic (27,27.25] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 1717 2.3210685 7778.157816 14.1099732 1982-1999 –> 2000-2012 4.1730415 16922.02063
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.57230
Atlantic (27,27.25] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 1717 2.0406890 7336.064154 8.5013493 1982-1999 –> 2000-2012 3.6834472 15942.16805
Atlantic (27,27.25] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 982 2.4537027 4561.678099 14.9835525 2000-2012 –> 2013-2019 4.6724715 12183.04779
Atlantic (27,27.25] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 982 2.4198641 4536.404447 14.2305742 2000-2012 –> 2013-2019 4.4271125 11815.72975
Atlantic (27,27.25] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 982 2.4324626 4546.603049 14.3792540 2000-2012 –> 2013-2019 4.4684162 11874.68933
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.29781
Atlantic (27,27.25] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 982 2.4360162 4549.470195 14.3742943 2000-2012 –> 2013-2019 4.4767052 11885.53435
Atlantic (27.25,27.5] 4 cstar_tref ~ sal + oxygen + aou + phosphate_star 1621 2.8128853 7965.109998 20.4458923 1982-1999 –> 2000-2012 5.1118884 18103.97065
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 1621 2.8109649 7964.895905 20.2408438 1982-1999 –> 2000-2012 5.1068355 18099.62362
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1621 2.8948304 8060.206531 21.8649642 1982-1999 –> 2000-2012 5.3345237 18468.23380
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.36215
Atlantic (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1621 2.8001886 7954.443216 19.6927272 1982-1999 –> 2000-2012 5.0677704 18035.40432
Atlantic (27.25,27.5] 4 cstar_tref ~ sal + oxygen + aou + phosphate_star 904 3.6637613 4925.111334 23.1796634 2000-2012 –> 2013-2019 6.4766466 12890.22133
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 904 3.6351396 4912.931565 22.4385896 2000-2012 –> 2013-2019 6.4461045 12877.82747
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.04468
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.20889
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 904 3.7450377 4966.781357 24.4908360 2000-2012 –> 2013-2019 6.6390002 13026.01573
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2079 2.3539368 9473.565176 17.1082120 1982-1999 –> 2000-2012 4.4979050 21665.23929
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.25739
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.62356
Atlantic (27.5,27.75] 4 cstar_tref ~ sal + temp + oxygen + phosphate_star 2079 2.3560092 9475.224293 16.7638042 1982-1999 –> 2000-2012 4.5113248 21694.36440
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 2079 2.3546463 9474.818186 17.1405940 1982-1999 –> 2000-2012 4.4968909 21662.00303
Atlantic (27.5,27.75] 4 cstar_tref ~ sal + temp + aou + phosphate_star 1115 2.6930488 5385.435813 18.3772087 2000-2012 –> 2013-2019 5.0483685 14859.44294
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1115 2.6676325 5366.289734 19.9843851 2000-2012 –> 2013-2019 5.0215694 14839.85491
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.01984
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.80519
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 1115 2.6715610 5369.571260 20.0496977 2000-2012 –> 2013-2019 5.0262072 14844.38945
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 722 1.6267617 2765.585146 14.5734851 1982-1999 –> 2000-2012 3.2965550 6480.60424
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 722 1.6517965 2787.638078 14.7352390 1982-1999 –> 2000-2012 3.3409598 6524.75540
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.03163
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.45688
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 722 1.6535346 2789.156784 14.7399711 1982-1999 –> 2000-2012 3.3439506 6527.69445
Atlantic (27.75,27.85] 4 cstar_tref ~ sal + oxygen + aou + phosphate_star 364 1.7438081 1449.807133 12.2711012 2000-2012 –> 2013-2019 3.3792145 4221.04551
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 364 1.7268302 1444.684538 13.7518672 2000-2012 –> 2013-2019 3.3535919 4210.26968
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 364 1.7547046 1456.342042 14.1928028 2000-2012 –> 2013-2019 3.4065011 4243.98012
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.75862
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.22418
Atlantic (27.85,27.95] 4 cstar_tref ~ oxygen + aou + silicate + phosphate_star 694 2.8412740 3430.909237 20.0203261 1982-1999 –> 2000-2012 5.7998788 8259.94399
Atlantic (27.85,27.95] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 694 2.7709000 3398.097680 20.1418519 1982-1999 –> 2000-2012 5.7244719 8225.85675
Atlantic (27.85,27.95] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 694 2.7808350 3403.065411 20.6486441 1982-1999 –> 2000-2012 5.8720788 8318.47864
Atlantic (27.85,27.95] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 694 2.6858576 3356.830803 20.1208556 1982-1999 –> 2000-2012 5.6386261 8186.06643
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.48649
Atlantic (27.85,27.95] 4 cstar_tref ~ oxygen + aou + silicate + phosphate_star 415 3.0542817 2116.450869 18.0269148 2000-2012 –> 2013-2019 5.8955557 5547.36011
Atlantic (27.85,27.95] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 415 3.0317519 2112.305705 17.9822303 2000-2012 –> 2013-2019 5.8026519 5510.40338
Atlantic (27.85,27.95] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 415 3.1817733 2152.393097 18.7816098 2000-2012 –> 2013-2019 5.9626083 5555.45851
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.10749
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.13824
Atlantic (27.95,28.05] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 723 3.1223151 3712.164200 22.1400867 1982-1999 –> 2000-2012 6.8744392 9170.30781
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.13227
Atlantic (27.95,28.05] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 723 3.1176154 3709.986044 21.8324808 1982-1999 –> 2000-2012 6.8686487 9167.55224
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.00761
Atlantic (27.95,28.05] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 723 3.1228707 3710.421472 22.8193044 1982-1999 –> 2000-2012 6.8769297 9167.58898
Atlantic (27.95,28.05] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 386 2.8383959 1914.801104 13.9583534 2000-2012 –> 2013-2019 5.9607110 5626.96530
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.30473
Atlantic (27.95,28.05] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 386 2.8350282 1913.884600 13.8967441 2000-2012 –> 2013-2019 5.9526436 5623.87064
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.83694
Atlantic (27.95,28.05] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 386 2.8484126 1915.520690 14.0128878 2000-2012 –> 2013-2019 5.9712832 5625.94216
Atlantic (28.05,28.1] 4 cstar_tref ~ oxygen + aou + silicate + phosphate_star 436 0.2179047 -79.349662 1.3640876 1982-1999 –> 2000-2012 0.3337049 -947.26878
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.71285
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 436 0.2223054 -59.914873 1.4246008 1982-1999 –> 2000-2012 0.3404537 -901.86471
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.16703
Atlantic (28.05,28.1] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 436 0.2176361 -78.425344 1.3751833 1982-1999 –> 2000-2012 0.3331923 -946.86274
Atlantic (28.05,28.1] 4 cstar_tref ~ sal + oxygen + aou + phosphate_star 229 0.2500674 27.074573 0.9325065 2000-2012 –> 2013-2019 0.4621040 -76.08002
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 229 0.2472173 23.824553 1.0717180 2000-2012 –> 2013-2019 0.4581585 -81.84627
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.21867
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.67186
Atlantic (28.05,28.1] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 229 0.2589000 44.972401 1.3338851 2000-2012 –> 2013-2019 0.4765361 -33.45294
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 541 0.1394443 -582.345612 1.0824605 1982-1999 –> 2000-2012 0.2743304 -1472.97554
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 541 0.1428962 -555.887490 1.1218747 1982-1999 –> 2000-2012 0.2796776 -1424.91708
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.25535
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 541 0.1430848 -554.460424 1.1249335 1982-1999 –> 2000-2012 0.2799199 -1422.88225
Atlantic (28.1,28.15] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 541 0.1566336 -456.569795 1.2740811 1982-1999 –> 2000-2012 0.3008900 -1243.23369
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.97959
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 288 0.1322538 -333.958174 0.4590119 2000-2012 –> 2013-2019 0.2751500 -889.84566
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 288 0.1409148 -297.420897 0.6451259 2000-2012 –> 2013-2019 0.2886767 -817.07924
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.40216
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 288 0.1325312 -332.751023 0.4591722 2000-2012 –> 2013-2019 0.2756160 -887.21145
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 1155 0.2973070 489.741021 1.1555289 1982-1999 –> 2000-2012 0.6778125 2072.74363
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1155 0.2982610 497.141683 1.1508210 1982-1999 –> 2000-2012 0.6794854 2086.68593
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.91861
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 1155 0.2982387 496.968549 1.1503200 1982-1999 –> 2000-2012 0.6794721 2086.59522
Atlantic (28.15,28.2] 4 cstar_tref ~ sal + temp + silicate + phosphate_star 1155 0.3004695 512.183308 1.2576325 1982-1999 –> 2000-2012 0.6828614 2110.32575
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 636 0.3580823 512.555409 2.3528717 2000-2012 –> 2013-2019 0.6553893 1002.29643
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 636 0.3589840 515.754532 2.3318667 2000-2012 –> 2013-2019 0.6572450 1012.89622
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.43517
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 636 0.3589968 515.799831 2.3317557 2000-2012 –> 2013-2019 0.6572355 1012.76838
Atlantic (28.15,28.2] 4 cstar_tref ~ sal + temp + silicate + phosphate_star 636 0.3605425 519.264825 2.3812408 2000-2012 –> 2013-2019 0.6610120 1031.44813
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.25281
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 955 0.3377975 651.233205 1.9994147 1982-1999 –> 2000-2012 0.8860519 2692.04536
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.96723
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 955 0.3377542 650.988114 2.0093171 1982-1999 –> 2000-2012 0.8860409 2691.94645
Atlantic (28.2, Inf] 4 cstar_tref ~ sal + temp + silicate + phosphate_star 955 0.3378894 649.752409 1.9553433 1982-1999 –> 2000-2012 0.8890811 2701.80551
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 552 0.5028733 821.599741 3.7448893 2000-2012 –> 2013-2019 0.8330790 1429.41714
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 552 0.5001058 815.507256 3.5900054 2000-2012 –> 2013-2019 0.8379033 1466.74046
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.46638
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.88155
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 552 0.4999646 815.195549 3.5780560 2000-2012 –> 2013-2019 0.8377188 1466.18366
Indo-Pacific (-Inf,26] 4 cstar_tref ~ sal + oxygen + aou + phosphate_star 5288 6.2360939 34376.517925 29.9685673 1982-1999 –> 2000-2012 12.7664702 81783.00523
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.29062
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5288 6.3788089 34617.825136 29.9652791 1982-1999 –> 2000-2012 13.0199305 82268.16435
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 5288 6.2323348 34372.140776 29.8694812 1982-1999 –> 2000-2012 12.7612356 81777.37841
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.52295
Indo-Pacific (-Inf,26] 4 cstar_tref ~ sal + oxygen + aou + phosphate_star 2815 5.9530030 18043.957291 24.8937268 2000-2012 –> 2013-2019 12.1890969 52420.47522
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 2815 5.8787943 17975.333992 24.1274152 2000-2012 –> 2013-2019 12.0716422 52280.25395
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2815 6.0477578 18134.865093 25.1292575 2000-2012 –> 2013-2019 12.4265667 52752.69023
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 2815 5.9511813 18044.234154 24.7705850 2000-2012 –> 2013-2019 12.1835161 52416.37493
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.73820
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.37494
Indo-Pacific (26,26.5] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5428 4.6702501 32149.400985 38.5622127 1982-1999 –> 2000-2012 9.0197595 76053.11295
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.69317
Indo-Pacific (26,26.5] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 5428 4.6786070 32168.809282 39.5746597 1982-1999 –> 2000-2012 9.0396276 76112.67451
Indo-Pacific (26,26.5] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 5428 4.8704336 32605.031006 41.2998501 1982-1999 –> 2000-2012 9.4796435 77389.77888
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.89050
Indo-Pacific (26,26.5] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2876 4.6511214 17017.181667 36.5566719 2000-2012 –> 2013-2019 9.3213715 49166.58265
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.73436
Indo-Pacific (26,26.5] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 2876 4.6557423 17022.893429 36.3745282 2000-2012 –> 2013-2019 9.3343493 49191.70271
Indo-Pacific (26,26.5] 4 cstar_tref ~ sal + temp + silicate + phosphate_star 2876 4.7308237 17112.913778 34.8668443 2000-2012 –> 2013-2019 9.6494996 49822.94581
Indo-Pacific (26.5,26.75] 4 cstar_tref ~ oxygen + aou + silicate + phosphate_star 4269 4.0546161 24078.867852 22.3581891 1982-1999 –> 2000-2012 8.0016719 57191.73449
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.61516
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 56219.99995
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 4269 4.1376444 24253.938495 25.6713855 1982-1999 –> 2000-2012 8.0635838 57305.20648
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 4269 3.9458514 23848.708841 23.7921959 1982-1999 –> 2000-2012 7.8106187 56713.78794
Indo-Pacific (26.5,26.75] 4 cstar_tref ~ oxygen + aou + silicate + phosphate_star 2324 4.4254119 13520.491182 24.9112112 2000-2012 –> 2013-2019 8.4800280 37599.35903
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 2324 4.3199526 13410.386398 24.6916724 2000-2012 –> 2013-2019 8.2637646 37254.68132
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.32121
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 2324 4.5083111 13608.754592 23.7952691 2000-2012 –> 2013-2019 8.6459555 37862.69309
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 2324 4.3137865 13403.747333 23.1641121 2000-2012 –> 2013-2019 8.2596379 37252.45617
Indo-Pacific (26.75,27] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 4557 4.3862266 26420.975208 21.0605286 1982-1999 –> 2000-2012 8.0872756 61505.87489
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.17863
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.49880
Indo-Pacific (26.75,27] 4 cstar_tref ~ temp + oxygen + aou + phosphate_star 4557 4.1598161 25935.947216 23.1614170 1982-1999 –> 2000-2012 7.7437292 60605.32053
Indo-Pacific (26.75,27] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 4557 4.1362281 25886.119686 21.9180570 1982-1999 –> 2000-2012 7.7019358 60492.01106
Indo-Pacific (26.75,27] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 2460 4.7910605 14703.596395 19.4544536 2000-2012 –> 2013-2019 9.1772872 41124.57160
Indo-Pacific (26.75,27] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 2460 4.6048662 14508.576616 17.1326056 2000-2012 –> 2013-2019 8.7604987 40437.35307
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.87403
Indo-Pacific (26.75,27] 4 cstar_tref ~ temp + oxygen + aou + phosphate_star 2460 4.6128201 14515.067487 16.6134065 2000-2012 –> 2013-2019 8.7726362 40451.01470
Indo-Pacific (26.75,27] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 2460 4.5994232 14502.757680 16.6240534 2000-2012 –> 2013-2019 8.7356513 40388.87737
Indo-Pacific (27,27.25] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 4267 4.2143064 24399.252678 44.7114638 1982-1999 –> 2000-2012 7.5507703 54676.73929
Indo-Pacific (27,27.25] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 4267 3.7349914 23368.861623 42.9143849 1982-1999 –> 2000-2012 6.7282481 52394.33781
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.99091
Indo-Pacific (27,27.25] 4 cstar_tref ~ temp + oxygen + aou + phosphate_star 4267 3.7387791 23375.511647 42.2525000 1982-1999 –> 2000-2012 6.7614943 52511.94640
Indo-Pacific (27,27.25] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 4267 3.4041802 22577.406832 30.9255112 1982-1999 –> 2000-2012 6.2015505 50822.23662
Indo-Pacific (27,27.25] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 2344 4.5554949 13774.558428 35.2893459 2000-2012 –> 2013-2019 8.2904863 37143.42005
Indo-Pacific (27,27.25] 5 cstar_tref ~ sal + temp + oxygen + aou + silicate 2344 4.7694276 13989.700735 41.9566499 2000-2012 –> 2013-2019 9.2207679 38855.93567
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.76198
Indo-Pacific (27,27.25] 4 cstar_tref ~ temp + oxygen + aou + phosphate_star 2344 4.5558378 13772.911312 35.4736376 2000-2012 –> 2013-2019 8.2946169 37148.42296
Indo-Pacific (27,27.25] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 2344 3.9533412 13109.926205 24.4569335 2000-2012 –> 2013-2019 7.3575214 35687.33304
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 3520 3.1346725 18046.701221 29.1767786 1982-1999 –> 2000-2012 6.0487445 41537.11519
Indo-Pacific (27.25,27.5] 4 cstar_tref ~ sal + oxygen + silicate + phosphate_star 3520 3.1493898 18077.676649 29.8455275 1982-1999 –> 2000-2012 6.0757661 41605.84076
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.81559
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 3520 3.1464564 18073.116463 29.5363811 1982-1999 –> 2000-2012 6.0679684 41587.58622
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.37457
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 1879 3.6867091 10249.562204 27.5141850 2000-2012 –> 2013-2019 6.8213816 28296.26342
Indo-Pacific (27.25,27.5] 4 cstar_tref ~ sal + oxygen + silicate + phosphate_star 1879 3.6875866 10248.456571 27.7368313 2000-2012 –> 2013-2019 6.8369764 28326.13322
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.89115
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 1879 3.6852335 10248.057760 28.1152441 2000-2012 –> 2013-2019 6.8316899 28321.17422
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.41416
Indo-Pacific (27.5,27.75] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 3170 1.9027986 13086.755512 12.4473296 1982-1999 –> 2000-2012 3.5535107 30391.58410
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 3170 1.8834035 13023.800742 9.5549910 1982-1999 –> 2000-2012 3.5338959 30329.43049
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 3170 1.8822945 13020.066462 9.5359358 1982-1999 –> 2000-2012 3.5327824 30325.67165
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.08642
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 3170 1.8817537 13018.244743 9.5293548 1982-1999 –> 2000-2012 3.5322649 30323.97696
Indo-Pacific (27.5,27.75] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 1683 2.2370458 7498.302788 13.7718692 2000-2012 –> 2013-2019 4.1398444 20585.05830
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 1683 2.1483716 7364.161462 12.8217994 2000-2012 –> 2013-2019 4.0317751 20387.96220
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1683 2.1461494 7360.678050 12.7660610 2000-2012 –> 2013-2019 4.0284439 20380.74451
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.05514
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 1683 2.1454909 7359.645130 12.7476542 2000-2012 –> 2013-2019 4.0272446 20377.88987
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 1216 1.4268879 4329.424189 7.6562558 1982-1999 –> 2000-2012 2.7551730 10092.16905
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1216 1.4351721 4343.503102 7.9724211 1982-1999 –> 2000-2012 2.7646432 10109.26081
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.33983
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.96521
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 1216 1.3124934 4126.189160 6.8457618 1982-1999 –> 2000-2012 2.5618830 9682.20926
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 691 1.6140386 2636.590993 13.0991822 2000-2012 –> 2013-2019 3.0409264 6966.01518
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 691 1.6190756 2640.897202 13.5827675 2000-2012 –> 2013-2019 3.0542477 6984.40030
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.68176
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.28295
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 691 1.5393786 2571.138596 10.1972442 2000-2012 –> 2013-2019 2.8518720 6697.32776
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 1797 1.3767808 6262.839408 10.8206192 1982-1999 –> 2000-2012 2.6534586 14695.95693
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.63724
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ sal + temp + oxygen + aou + silicate 1797 1.2373876 5879.195630 6.8169027 1982-1999 –> 2000-2012 2.3851572 13773.51013
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.62250
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 1797 1.3914252 6300.866006 9.9533251 1982-1999 –> 2000-2012 2.6565451 14687.94739
Indo-Pacific (27.85,27.95] 4 cstar_tref ~ sal + oxygen + aou + silicate 987 1.6167964 3761.386363 12.6522924 2000-2012 –> 2013-2019 3.0074315 10058.21095
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 987 1.5919784 3732.850266 10.5660006 2000-2012 –> 2013-2019 2.9687592 9995.68967
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 987 1.4361554 3529.512839 8.9979380 2000-2012 –> 2013-2019 2.7232012 9550.12217
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.41753
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.87429
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 1162 0.5801738 2046.363310 3.0668291 1982-1999 –> 2000-2012 1.1222674 4865.79437
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1162 0.5812399 2050.630181 3.1553323 1982-1999 –> 2000-2012 1.1253410 4882.91731
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.49345
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 1162 0.5813748 2051.169436 3.1662919 1982-1999 –> 2000-2012 1.1257026 4884.90525
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.73873
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 648 0.6051146 1201.915076 2.6426945 2000-2012 –> 2013-2019 1.1852884 3248.27839
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 648 0.6068029 1205.525847 2.6420479 2000-2012 –> 2013-2019 1.1880428 3256.15603
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.63325
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 648 0.6069777 1205.899106 2.6420805 2000-2012 –> 2013-2019 1.1883525 3257.06854
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.26128
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 1015 0.4567175 1303.553975 2.4275171 1982-1999 –> 2000-2012 0.9128334 3243.43396
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1015 0.4570675 1305.109306 2.4292266 1982-1999 –> 2000-2012 0.9134348 3246.66276
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.29740
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 1015 0.4571160 1305.324491 2.4294579 1982-1999 –> 2000-2012 0.9135298 3247.18817
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.22099
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 565 0.4757060 777.861154 2.6220126 2000-2012 –> 2013-2019 0.9324235 2081.41513
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 565 0.4759048 778.333277 2.6231400 2000-2012 –> 2013-2019 0.9329724 2083.44258
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.75611
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 565 0.4759340 778.402535 2.6232713 2000-2012 –> 2013-2019 0.9330500 2083.72703
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 565 0.4756293 777.678976 2.6025948 2000-2012 –> 2013-2019 0.9302630 2071.94948
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.14490
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 9318 0.7775186 21767.632646 4.8314997 1982-1999 –> 2000-2012 1.5398228 51562.30395
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.67513
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 9318 0.7775458 21768.284693 4.8316142 1982-1999 –> 2000-2012 1.5398616 51563.35068
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 9318 0.7744993 21695.121826 4.9051251 1982-1999 –> 2000-2012 1.5417070 51656.19560
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 5093 0.7763424 11888.603170 4.2785787 2000-2012 –> 2013-2019 1.5528339 33631.60210
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5093 0.7779155 11909.222793 4.2755260 2000-2012 –> 2013-2019 1.5554342 33676.85544
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.81792
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 5093 0.7779515 11909.693562 4.2760690 2000-2012 –> 2013-2019 1.5554973 33677.97826
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 5093 0.7692867 11795.605640 4.3360410 2000-2012 –> 2013-2019 1.5437859 33490.72747
SO (26,26.5] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 103 2.7254296 512.842306 11.4100082 1982-1999 –> 2000-2012 6.2111286 1215.08495
SO (26,26.5] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 103 2.7471036 514.474048 11.2693383 1982-1999 –> 2000-2012 6.2350014 1216.87938
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.97576
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.84184
SO (26,26.5] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 103 2.8011962 518.490922 11.6880130 1982-1999 –> 2000-2012 6.2394349 1217.19657
SO (26,26.5] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 47 3.0953135 253.589806 8.2440313 2000-2012 –> 2013-2019 5.8773819 770.66925
SO (26,26.5] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 47 3.0868892 253.333626 8.0486275 2000-2012 –> 2013-2019 5.8123188 766.17593
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.11970
SO (26,26.5] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 47 3.0875443 253.353570 8.0376410 2000-2012 –> 2013-2019 5.8122599 766.14191
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.84725
SO (26.5,26.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 312 3.2867067 1641.906542 10.0386801 1982-1999 –> 2000-2012 6.9067437 3928.46167
SO (26.5,26.75] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 312 3.2153502 1628.209882 10.5226323 1982-1999 –> 2000-2012 6.7106987 3885.32203
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.14691
SO (26.5,26.75] 4 cstar_tref ~ temp + oxygen + aou + phosphate_star 312 3.2799039 1638.613652 10.5986533 1982-1999 –> 2000-2012 6.8388507 3908.87235
SO (26.5,26.75] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 312 3.2469942 1634.320979 10.2848760 1982-1999 –> 2000-2012 6.8034551 3905.99273
SO (26.5,26.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 133 3.2899907 708.212988 10.8830055 2000-2012 –> 2013-2019 6.5766974 2350.11953
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.17596
SO (26.5,26.75] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 133 3.2903253 708.240039 10.8891344 2000-2012 –> 2013-2019 6.5782262 2350.37328
SO (26.5,26.75] 4 cstar_tref ~ sal + temp + silicate + phosphate_star 133 3.2932626 706.477397 10.9049735 2000-2012 –> 2013-2019 6.5827544 2346.91248
SO (26.5,26.75] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 133 3.3148729 710.217175 10.3205968 2000-2012 –> 2013-2019 6.5618670 2344.53815
SO (26.75,27] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1061 2.8283940 5231.250170 9.1586744 1982-1999 –> 2000-2012 5.7833809 12827.62141
SO (26.75,27] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 1061 2.3650556 4851.608393 11.6713433 1982-1999 –> 2000-2012 4.8527893 11926.44971
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.55673
SO (26.75,27] 4 cstar_tref ~ temp + oxygen + aou + phosphate_star 1061 2.3658852 4850.352635 11.7790894 1982-1999 –> 2000-2012 4.8536773 11923.26504
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.12478
SO (26.75,27] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 550 2.8763685 2737.013799 9.5868252 2000-2012 –> 2013-2019 5.7047625 7968.26397
SO (26.75,27] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 550 2.3941154 2535.147577 15.5973947 2000-2012 –> 2013-2019 4.7591710 7386.75597
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.28542
SO (26.75,27] 4 cstar_tref ~ temp + oxygen + aou + phosphate_star 550 2.3957372 2533.892462 15.4998931 2000-2012 –> 2013-2019 4.7616224 7384.24510
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.93385
SO (27,27.25] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2513 1.8324114 10189.495572 14.0129707 1982-1999 –> 2000-2012 3.6402091 24031.46072
SO (27,27.25] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 2513 1.7156657 9858.625769 16.4970589 1982-1999 –> 2000-2012 3.3035748 22808.83415
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.35801
SO (27,27.25] 4 cstar_tref ~ temp + oxygen + aou + phosphate_star 2513 1.7159114 9857.345457 16.6672621 1982-1999 –> 2000-2012 3.3256604 22899.48106
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.70242
SO (27,27.25] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1286 1.9306015 5355.452856 8.6195738 2000-2012 –> 2013-2019 3.7630129 15544.94843
SO (27,27.25] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 1286 1.9259139 5349.200280 11.1221020 2000-2012 –> 2013-2019 3.6415796 15207.82605
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.70269
SO (27,27.25] 4 cstar_tref ~ temp + oxygen + aou + phosphate_star 1286 1.9277877 5349.701445 11.5137595 2000-2012 –> 2013-2019 3.6436991 15207.04690
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.89494
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.64538
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.75864
SO (27.25,27.5] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 2731 1.6889628 10626.956131 10.5782629 1982-1999 –> 2000-2012 3.1276979 24203.75969
SO (27.25,27.5] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 2731 1.6893649 10626.256647 10.5431309 1982-1999 –> 2000-2012 3.1306658 24214.61540
SO (27.25,27.5] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 2731 1.6893569 10628.230775 10.5492429 1982-1999 –> 2000-2012 3.1276338 24202.61058
SO (27.25,27.5] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1453 1.8486094 5922.979657 8.8767744 2000-2012 –> 2013-2019 3.5373288 16549.14871
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.73018
SO (27.25,27.5] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 1453 1.8467306 5920.024690 8.8926590 2000-2012 –> 2013-2019 3.5356934 16546.98082
SO (27.25,27.5] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 1453 1.8579855 5937.681599 9.5635309 2000-2012 –> 2013-2019 3.5473425 16565.91237
SO (27.25,27.5] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 1453 1.8592539 5937.664674 9.5274690 2000-2012 –> 2013-2019 3.5489602 16565.02492
SO (27.5,27.75] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 3391 1.5320138 12530.327695 9.1542418 1982-1999 –> 2000-2012 2.8147095 27705.56033
SO (27.5,27.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 3391 1.5327196 12533.451505 9.1052349 1982-1999 –> 2000-2012 2.8175679 27723.92592
SO (27.5,27.75] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 3391 1.5347080 12542.244087 8.6514924 1982-1999 –> 2000-2012 2.8178293 27720.49199
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.16027
SO (27.5,27.75] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 3391 1.5319189 12529.907533 9.1108749 1982-1999 –> 2000-2012 2.8155627 27711.85689
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.49316
SO (27.5,27.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1821 1.7574181 7235.300257 10.3457829 2000-2012 –> 2013-2019 3.2901377 19768.75176
SO (27.5,27.75] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 1821 1.7645761 7250.104197 9.5282977 2000-2012 –> 2013-2019 3.2992842 19792.34828
SO (27.5,27.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1821 1.7554202 7233.157544 10.3879446 2000-2012 –> 2013-2019 3.2869838 19763.49181
SO (27.5,27.75] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 1821 1.7560943 7232.555962 10.3552020 2000-2012 –> 2013-2019 3.2880132 19762.46349
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.39033
SO (27.75,27.85] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1415 1.4814190 5141.787239 13.5369619 1982-1999 –> 2000-2012 2.7844305 11233.64330
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.56712
SO (27.75,27.85] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 1415 1.4825470 5143.941155 13.5477979 1982-1999 –> 2000-2012 2.7869620 11239.68368
SO (27.75,27.85] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 1415 1.5002561 5177.545528 13.7566640 1982-1999 –> 2000-2012 2.7734180 11185.74101
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.43122
SO (27.75,27.85] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 722 1.5955525 2737.613017 8.5704585 2000-2012 –> 2013-2019 3.0769715 7879.40026
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.55165
SO (27.75,27.85] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 722 1.5956695 2737.718886 8.5534574 2000-2012 –> 2013-2019 3.0782164 7881.66004
SO (27.75,27.85] 4 cstar_tref ~ sal + temp + silicate + phosphate_star 722 1.5960379 2736.052251 8.4492259 2000-2012 –> 2013-2019 3.0932979 7905.94044
SO (27.85,27.95] 4 cstar_tref ~ sal + oxygen + aou + silicate 1399 1.4396856 5001.850417 9.7403080 1982-1999 –> 2000-2012 2.8655377 11587.21281
SO (27.85,27.95] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 1399 1.4396129 5003.709245 9.8030845 1982-1999 –> 2000-2012 2.8546010 11562.72636
SO (27.85,27.95] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 1399 1.4232287 4971.682542 7.2306924 1982-1999 –> 2000-2012 2.7624199 11326.66498
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.19404
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.56053
SO (27.85,27.95] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 709 1.6432910 2730.376750 7.3480083 2000-2012 –> 2013-2019 3.0665196 7702.05929
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.17557
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.48411
SO (27.85,27.95] 4 cstar_tref ~ temp + oxygen + aou + silicate 709 1.5425225 2638.643090 8.9509076 2000-2012 –> 2013-2019 3.0269454 7726.11609
SO (27.85,27.95] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 709 1.5409988 2639.241700 8.7822837 2000-2012 –> 2013-2019 2.9762393 7634.43986
SO (27.95,28.05] 4 cstar_tref ~ sal + oxygen + aou + phosphate_star 1786 1.7323419 7043.170193 10.1736604 1982-1999 –> 2000-2012 3.2195838 16087.26729
SO (27.95,28.05] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 1786 1.7234592 7026.807513 9.6018392 1982-1999 –> 2000-2012 3.1951756 16020.70736
SO (27.95,28.05] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 1786 1.5620816 6675.629386 9.1134438 1982-1999 –> 2000-2012 2.8752317 15102.49093
SO (27.95,28.05] 5 cstar_tref ~ sal + temp + oxygen + aou + silicate 1786 1.5619964 6675.434498 7.6946622 1982-1999 –> 2000-2012 2.8752447 15102.66796
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.16971
SO (27.95,28.05] 4 cstar_tref ~ sal + oxygen + aou + silicate 986 1.8090331 3979.133623 8.8118450 2000-2012 –> 2013-2019 3.5496875 11039.40270
SO (27.95,28.05] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 986 1.8071065 3979.032373 8.8048453 2000-2012 –> 2013-2019 3.5305658 11005.83989
SO (27.95,28.05] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 986 1.6936583 3851.175591 8.4129789 2000-2012 –> 2013-2019 3.2557399 10526.80498
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.88714
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.19367
SO (28.05,28.1] 4 cstar_tref ~ sal + oxygen + aou + phosphate_star 1071 1.1929018 3429.191146 10.1579725 1982-1999 –> 2000-2012 2.2787248 7987.05738
SO (28.05,28.1] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 1071 1.1927699 3430.954357 10.1729214 1982-1999 –> 2000-2012 2.2759473 7983.43376
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.68729
SO (28.05,28.1] 5 cstar_tref ~ sal + temp + oxygen + aou + silicate 1071 1.0700400 3198.371237 9.2962436 1982-1999 –> 2000-2012 2.0713151 7512.77554
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.67458
SO (28.05,28.1] 4 cstar_tref ~ sal + oxygen + aou + phosphate_star 606 1.4121339 2150.017064 11.3700248 2000-2012 –> 2013-2019 2.6050356 5579.20821
SO (28.05,28.1] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 606 1.4121202 2152.005351 11.3570568 2000-2012 –> 2013-2019 2.6048901 5582.95971
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.28387
SO (28.05,28.1] 5 cstar_tref ~ sal + temp + oxygen + aou + silicate 606 1.2597549 2013.625136 10.8095062 2000-2012 –> 2013-2019 2.3297949 5211.99637
SO (28.05,28.1] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 606 1.2282847 1984.963330 10.1968303 2000-2012 –> 2013-2019 2.2756697 5139.49687
SO (28.1, Inf] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 5309 1.1018770 16108.395277 6.0859562 1982-1999 –> 2000-2012 2.0912319 36846.81840
SO (28.1, Inf] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 5309 1.1018498 16110.133861 6.1005720 1982-1999 –> 2000-2012 2.0809310 36696.92150
SO (28.1, Inf] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5309 1.1018543 16110.177014 6.0992752 1982-1999 –> 2000-2012 2.0810102 36698.08779
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.92485
SO (28.1, Inf] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 5309 1.1019213 16110.822324 6.0983887 1982-1999 –> 2000-2012 2.0811392 36699.66564
SO (28.1, Inf] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 2711 1.1946469 8669.790861 6.0803943 2000-2012 –> 2013-2019 2.2965238 24778.18614
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.31567
SO (28.1, Inf] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2711 1.1918665 8659.157262 6.0517786 2000-2012 –> 2013-2019 2.2937208 24769.33428
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.47146
SO (28.1, Inf] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 2711 1.1918732 8659.187802 6.0524883 2000-2012 –> 2013-2019 2.2937945 24770.01013
SO (28.1,28.15] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 395 0.4692853 537.291343 1.7413747 1982-1999 –> 2000-2012 0.9103638 1176.91460
SO (28.1,28.15] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 395 0.4712779 540.638668 1.7559543 1982-1999 –> 2000-2012 0.9142886 1184.81639
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.77071
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.23170
SO (28.1,28.15] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 395 0.4711414 540.409822 1.7546095 1982-1999 –> 2000-2012 0.9140250 1184.28863
SO (28.1,28.15] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 209 0.5007130 317.976471 2.3778112 2000-2012 –> 2013-2019 0.9699983 855.26781
SO (28.1,28.15] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 209 0.4863454 305.806779 1.9335059 2000-2012 –> 2013-2019 0.9551804 842.33972
SO (28.1,28.15] 5 cstar_tref ~ sal + temp + oxygen + aou + silicate 209 0.4928520 311.361923 2.0181406 2000-2012 –> 2013-2019 0.9605446 845.96764
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.65470
SO (28.1,28.15] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 209 0.5029754 319.860825 2.3841073 2000-2012 –> 2013-2019 0.9741168 860.27065
SO (28.15,28.2] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 385 0.3060639 194.932446 1.0624911 1982-1999 –> 2000-2012 0.6130546 442.64918
SO (28.15,28.2] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 385 0.3077233 199.095828 1.0703532 1982-1999 –> 2000-2012 0.6161479 451.38865
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.59322
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.01785
SO (28.15,28.2] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 385 0.3072492 197.908741 1.0678544 1982-1999 –> 2000-2012 0.6153464 449.15828
SO (28.15,28.2] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 183 0.3423788 141.038937 0.8031784 2000-2012 –> 2013-2019 0.6484427 335.97138
SO (28.15,28.2] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 183 0.3300672 127.635476 0.8682035 2000-2012 –> 2013-2019 0.6258955 296.37676
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.32871
SO (28.15,28.2] 4 cstar_tref ~ temp + oxygen + aou + phosphate_star 183 0.3366745 132.889698 0.8481580 2000-2012 –> 2013-2019 0.6489057 341.18383
SO (28.15,28.2] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 183 0.3337812 131.730824 0.8025278 2000-2012 –> 2013-2019 0.6459444 341.85710
SO (28.2, Inf] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 1738 0.3918565 1689.706304 2.1137168 1982-1999 –> 2000-2012 0.7169227 3163.23502
SO (28.2, Inf] 4 cstar_tref ~ sal + oxygen + silicate + phosphate_star 1738 0.4338139 2041.285083 2.8272800 1982-1999 –> 2000-2012 0.8013633 4120.07926
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.77632
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.29544
SO (28.2, Inf] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 1738 0.3914217 1685.847473 2.1087884 1982-1999 –> 2000-2012 0.7160679 3152.98326
SO (28.2, Inf] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate_star 1005 0.4243307 1143.009619 2.6684674 2000-2012 –> 2013-2019 0.8161872 2832.71592
SO (28.2, Inf] 4 cstar_tref ~ sal + oxygen + silicate + phosphate_star 1005 0.4467609 1244.545547 3.0434845 2000-2012 –> 2013-2019 0.8805748 3285.83063
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.52355
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.77141
SO (28.2, Inf] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 1005 0.4240694 1141.771704 2.6663305 2000-2012 –> 2013-2019 0.8154911 2827.61918

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.3843527 1202.0910
Atlantic (-Inf,26] 2000-2012 –> 2013-2019 2.2610678 724.4054
Atlantic (26,26.5] 1982-1999 –> 2000-2012 7.1497958 13765.5750
Atlantic (26,26.5] 2000-2012 –> 2013-2019 7.2173645 8680.9623
Atlantic (26.5,26.75] 1982-1999 –> 2000-2012 6.3155697 16412.6250
Atlantic (26.5,26.75] 2000-2012 –> 2013-2019 6.6132573 10506.3043
Atlantic (26.75,27] 1982-1999 –> 2000-2012 3.7107062 19584.2366
Atlantic (26.75,27] 2000-2012 –> 2013-2019 4.5905067 14464.7579
Atlantic (27,27.25] 1982-1999 –> 2000-2012 3.7399914 16049.6643
Atlantic (27,27.25] 2000-2012 –> 2013-2019 4.4706304 11887.6598
Atlantic (27.25,27.5] 1982-1999 –> 2000-2012 5.1381403 18148.3189
Atlantic (27.25,27.5] 2000-2012 –> 2013-2019 6.4732144 12899.2636
Atlantic (27.5,27.75] 1982-1999 –> 2000-2012 4.4678403 21601.4975
Atlantic (27.5,27.75] 2000-2012 –> 2013-2019 4.9703929 14779.1025
Atlantic (27.75,27.85] 1982-1999 –> 2000-2012 3.3009150 6486.7085
Atlantic (27.75,27.85] 2000-2012 –> 2013-2019 3.3346454 4199.2556
Atlantic (27.85,27.95] 1982-1999 –> 2000-2012 5.7352547 8234.7665
Atlantic (27.85,27.95] 2000-2012 –> 2013-2019 5.8151892 5509.2935
Atlantic (27.95,28.05] 1982-1999 –> 2000-2012 6.8516066 9159.5178
Atlantic (27.95,28.05] 2000-2012 –> 2013-2019 5.9310485 5614.7840
Atlantic (28.05,28.1] 1982-1999 –> 2000-2012 0.3282253 -966.3752
Atlantic (28.05,28.1] 2000-2012 –> 2013-2019 0.4492073 -109.0540
Atlantic (28.1,28.15] 1982-1999 –> 2000-2012 0.2808789 -1414.6528
Atlantic (28.1,28.15] 2000-2012 –> 2013-2019 0.2734646 -901.1036
Atlantic (28.15,28.2] 1982-1999 –> 2000-2012 0.6775554 2070.8538
Atlantic (28.15,28.2] 2000-2012 –> 2013-2019 0.6551443 999.1689
Atlantic (28.2, Inf] 1982-1999 –> 2000-2012 0.8775077 2640.2035
Atlantic (28.2, Inf] 2000-2012 –> 2013-2019 0.8305790 1388.9378
Indo-Pacific (-Inf,26] 1982-1999 –> 2000-2012 12.7778879 81805.8723
Indo-Pacific (-Inf,26] 2000-2012 –> 2013-2019 12.1877664 52429.7065
Indo-Pacific (26,26.5] 1982-1999 –> 2000-2012 9.0188688 76053.7269
Indo-Pacific (26,26.5] 2000-2012 –> 2013-2019 9.2682339 49091.9712
Indo-Pacific (26.5,26.75] 1982-1999 –> 2000-2012 7.8422650 56760.0688
Indo-Pacific (26.5,26.75] 2000-2012 –> 2013-2019 8.3600446 37410.7022
Indo-Pacific (26.75,27] 1982-1999 –> 2000-2012 7.7775168 60678.1768
Indo-Pacific (26.75,27] 2000-2012 –> 2013-2019 8.8340206 40555.5382
Indo-Pacific (27,27.25] 1982-1999 –> 2000-2012 6.6865375 52237.4502
Indo-Pacific (27,27.25] 2000-2012 –> 2013-2019 8.1028600 36903.5747
Indo-Pacific (27.25,27.5] 1982-1999 –> 2000-2012 5.9045403 41156.1465
Indo-Pacific (27.25,27.5] 2000-2012 –> 2013-2019 6.5778409 27906.1752
Indo-Pacific (27.5,27.75] 1982-1999 –> 2000-2012 3.5365973 30338.5499
Indo-Pacific (27.5,27.75] 2000-2012 –> 2013-2019 4.0503039 20421.1420
Indo-Pacific (27.75,27.85] 1982-1999 –> 2000-2012 2.5558244 9653.1888
Indo-Pacific (27.75,27.85] 2000-2012 –> 2013-2019 2.7913529 6625.1416
Indo-Pacific (27.85,27.95] 1982-1999 –> 2000-2012 2.4316251 13869.1348
Indo-Pacific (27.85,27.95] 2000-2012 –> 2013-2019 2.7591061 9577.6629
Indo-Pacific (27.95,28.05] 1982-1999 –> 2000-2012 1.1173895 4838.9698
Indo-Pacific (27.95,28.05] 2000-2012 –> 2013-2019 1.1785047 3231.2795
Indo-Pacific (28.05,28.1] 1982-1999 –> 2000-2012 0.9114933 3236.9607
Indo-Pacific (28.05,28.1] 2000-2012 –> 2013-2019 0.9315041 2078.4581
Indo-Pacific (28.1, Inf] 1982-1999 –> 2000-2012 1.5390922 51553.3341
Indo-Pacific (28.1, Inf] 2000-2012 –> 2013-2019 1.5501467 33591.5962
SO (26,26.5] 1982-1999 –> 2000-2012 6.1994144 1214.9957
SO (26,26.5] 2000-2012 –> 2013-2019 5.8375387 766.9908
SO (26.5,26.75] 1982-1999 –> 2000-2012 6.7725173 3899.1591
SO (26.5,26.75] 2000-2012 –> 2013-2019 6.5405788 2342.8239
SO (26.75,27] 1982-1999 –> 2000-2012 4.9888440 12052.8035
SO (26.75,27] 2000-2012 –> 2013-2019 4.8719454 7456.2969
SO (27,27.25] 1982-1999 –> 2000-2012 3.3393486 22957.7673
SO (27,27.25] 2000-2012 –> 2013-2019 3.5855006 15112.6838
SO (27.25,27.5] 1982-1999 –> 2000-2012 3.1271785 24200.8779
SO (27.25,27.5] 2000-2012 –> 2013-2019 3.5405150 16554.7594
SO (27.5,27.75] 1982-1999 –> 2000-2012 2.8158332 27712.7991
SO (27.5,27.75] 2000-2012 –> 2013-2019 3.2903713 19769.7097
SO (27.75,27.85] 1982-1999 –> 2000-2012 2.7676177 11189.2051
SO (27.75,27.85] 2000-2012 –> 2013-2019 3.0716064 7871.5967
SO (27.85,27.95] 1982-1999 –> 2000-2012 2.7159913 11225.6717
SO (27.85,27.95] 2000-2012 –> 2013-2019 2.9264670 7527.8550
SO (27.95,28.05] 1982-1999 –> 2000-2012 3.0020656 15465.6606
SO (27.95,28.05] 2000-2012 –> 2013-2019 3.3401593 10687.6257
SO (28.05,28.1] 1982-1999 –> 2000-2012 2.1277034 7631.5257
SO (28.05,28.1] 2000-2012 –> 2013-2019 2.4184871 5330.3890
SO (28.1, Inf] 1982-1999 –> 2000-2012 2.0828593 36726.2836
SO (28.1, Inf] 2000-2012 –> 2013-2019 2.2941740 24770.8635
SO (28.1,28.15] 1982-1999 –> 2000-2012 0.8989451 1153.4044
SO (28.1,28.15] 2000-2012 –> 2013-2019 0.9527659 836.3001
SO (28.15,28.2] 1982-1999 –> 2000-2012 0.6067970 426.1614
SO (28.15,28.2] 2000-2012 –> 2013-2019 0.6362469 318.5436
SO (28.2, Inf] 1982-1999 –> 2000-2012 0.7321128 3335.4739
SO (28.2, Inf] 2000-2012 –> 2013-2019 0.8274764 2912.6921

5.4 Target variable coefficients

A data frame to map the target variable is prepared.

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

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

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

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

lm_best_target <- left_join(lm_best_target, lm_all_fitted_wide)

rm(eras_era, eras_forward, eras_backward,
   lm_all_fitted)

5.5 Plot selected model residuals

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

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

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

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

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

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

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

5.6 Cant coeffcients

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

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

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

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

5.7 Write files

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

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

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

6 Model diagnotics

6.1 Selection criterion vs predictors

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

6.1.1 All models

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

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

6.1.2 Best models

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

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

6.2 RMSE correlation between eras

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

6.2.1 All models

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

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

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

6.2.2 Best models

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

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

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

6.3 Predictor counts

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

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

# print table
lm_all_stats %>%
  gt(rowname_col = "gamma_slab",
     groupname_col = c("basin", "eras")) %>% 
  summary_rows(
    groups = TRUE,
    fns = list(total = "sum")
  )
aou oxygen phosphate_star sal silicate temp
Atlantic - 1982-1999 --> 2000-2012
(-Inf,26] 5 4 5 5 3 3
(26,26.5] 3 3 5 5 5 4
(26.5,26.75] 4 4 5 4 5 4
(26.75,27] 5 4 5 3 5 3
(27,27.25] 4 4 5 5 4 4
(27.25,27.5] 5 4 5 5 3 3
(27.5,27.75] 3 4 5 5 3 5
(27.75,27.85] 4 4 5 5 4 4
(27.85,27.95] 5 5 5 3 4 3
(27.95,28.05] 3 4 5 3 5 5
(28.05,28.1] 5 4 5 3 5 3
(28.1,28.15] 4 4 5 4 5 4
(28.15,28.2] 3 3 5 5 5 4
(28.2, Inf] 3 3 5 5 5 4
total 56.00 54.00 70.00 60.00 61.00 53.00
Atlantic - 2000-2012 --> 2013-2019
(-Inf,26] 4 4 5 5 4 4
(26,26.5] 3 3 5 5 5 4
(26.5,26.75] 4 4 5 4 5 4
(26.75,27] 4 4 5 4 5 4
(27,27.25] 4 3 5 5 5 3
(27.25,27.5] 4 5 5 5 3 3
(27.5,27.75] 4 3 5 5 3 5
(27.75,27.85] 5 4 5 5 3 3
(27.85,27.95] 5 5 5 3 4 3
(27.95,28.05] 3 4 5 3 5 5
(28.05,28.1] 5 5 5 4 3 3
(28.1,28.15] 4 4 5 5 4 4
(28.15,28.2] 3 3 5 5 5 4
(28.2, Inf] 4 4 5 5 4 4
total 56.00 55.00 70.00 63.00 58.00 53.00
Indo-Pacific - 1982-1999 --> 2000-2012
(-Inf,26] 5 4 5 5 3 3
(26,26.5] 4 4 5 4 5 4
(26.5,26.75] 4 5 5 3 5 3
(26.75,27] 5 5 5 3 3 4
(27,27.25] 5 5 5 3 3 4
(27.25,27.5] 3 5 5 4 5 3
(27.5,27.75] 4 3 5 5 5 3
(27.75,27.85] 5 4 5 4 4 4
(27.85,27.95] 5 5 4 4 4 4
(27.95,28.05] 4 4 5 4 5 4
(28.05,28.1] 4 4 5 4 5 4
(28.1, Inf] 4 4 5 4 5 4
total 52.00 52.00 59.00 47.00 52.00 44.00
Indo-Pacific - 2000-2012 --> 2013-2019
(-Inf,26] 5 4 5 5 3 3
(26,26.5] 3 3 5 5 5 4
(26.5,26.75] 4 5 5 3 5 3
(26.75,27] 5 5 5 3 3 4
(27,27.25] 5 5 4 3 3 5
(27.25,27.5] 3 5 5 4 5 3
(27.5,27.75] 4 3 5 5 5 3
(27.75,27.85] 5 4 5 4 4 4
(27.85,27.95] 5 5 3 5 4 3
(27.95,28.05] 4 4 5 4 5 4
(28.05,28.1] 4 4 5 4 5 4
(28.1, Inf] 4 4 5 4 5 4
total 51.00 51.00 57.00 49.00 52.00 44.00
SO - 1982-1999 --> 2000-2012
(26,26.5] 4 4 5 4 4 5
(26.5,26.75] 5 4 5 3 3 5
(26.75,27] 5 4 5 3 3 5
(27,27.25] 5 4 5 3 3 5
(27.25,27.5] 4 3 5 3 5 5
(27.5,27.75] 4 4 5 5 4 4
(27.75,27.85] 4 4 5 4 5 4
(27.85,27.95] 5 5 3 5 4 3
(27.95,28.05] 5 5 4 5 3 3
(28.05,28.1] 5 5 4 5 3 3
(28.1, Inf] 4 3 5 5 5 3
(28.1,28.15] 4 4 5 5 4 4
(28.15,28.2] 4 4 5 5 4 4
(28.2, Inf] 3 4 5 5 5 3
total 61.00 57.00 66.00 60.00 55.00 56.00
SO - 2000-2012 --> 2013-2019
(26,26.5] 3 3 5 5 5 4
(26.5,26.75] 3 3 5 4 5 5
(26.75,27] 5 4 5 3 3 5
(27,27.25] 5 4 5 3 3 5
(27.25,27.5] 3 4 5 3 5 5
(27.5,27.75] 4 4 5 5 4 4
(27.75,27.85] 3 3 5 5 5 4
(27.85,27.95] 5 5 3 3 4 5
(27.95,28.05] 5 5 3 5 4 3
(28.05,28.1] 5 5 4 5 3 3
(28.1, Inf] 4 3 5 5 5 3
(28.1,28.15] 4 5 4 5 4 4
(28.15,28.2] 5 5 5 3 3 4
(28.2, Inf] 3 4 5 5 5 3
total 57.00 57.00 64.00 59.00 58.00 57.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
9b7bc66 Donghe-Zhu 2021-03-10
5365f80 Donghe-Zhu 2021-03-10
2b6c392 Donghe-Zhu 2021-03-10
9f58753 Donghe-Zhu 2021-03-10
9d5a62c Donghe-Zhu 2021-03-10
5d1e70b Donghe-Zhu 2021-03-10
2093979 Donghe-Zhu 2021-03-10
b865899 Donghe-Zhu 2021-03-10
cc2a956 Donghe-Zhu 2021-03-10
dba33c8 Donghe-Zhu 2021-03-09
17f1c4a Donghe-Zhu 2021-03-09
c024d1a Donghe-Zhu 2021-03-09
02f7242 Donghe-Zhu 2021-03-09
1691156 Donghe-Zhu 2021-03-08
c0ceaf8 Donghe-Zhu 2021-03-08
058e0a1 Donghe-Zhu 2021-03-08
112dea0 Donghe-Zhu 2021-03-08
1843412 Donghe-Zhu 2021-03-08
3fbbfa4 Donghe-Zhu 2021-03-07
627c8fb Donghe-Zhu 2021-03-07
9ef3222 Donghe-Zhu 2021-03-05
8c1e978 Donghe-Zhu 2021-03-05
865f68c Donghe-Zhu 2021-03-05
ee69bc1 Donghe-Zhu 2021-03-05
a79291f Donghe-Zhu 2021-03-05
e8c6f30 Donghe-Zhu 2021-03-04
59288fe Donghe-Zhu 2021-03-04
731abc8 Donghe-Zhu 2021-03-04
e2a5a33 Donghe-Zhu 2021-03-04
c7892c1 Donghe-Zhu 2021-03-04
924430b Donghe-Zhu 2021-03-03
0d0bca1 Donghe-Zhu 2021-03-03
cb63c16 Donghe-Zhu 2021-03-03
ffda45a Donghe-Zhu 2021-03-03
691ba81 Donghe-Zhu 2021-03-03
c5e45a2 Donghe-Zhu 2021-03-03
89c3e58 Donghe-Zhu 2021-03-03
c911669 Donghe-Zhu 2021-03-03
13666ca Donghe-Zhu 2021-03-01
c6e60fe Donghe-Zhu 2021-03-01
7a388f7 Donghe-Zhu 2021-03-01
799e913 Donghe-Zhu 2021-03-01
66ff99f Donghe-Zhu 2021-03-01
ac9bb7a Donghe-Zhu 2021-02-28
efdc047 Donghe-Zhu 2021-02-28
e9a7418 Donghe-Zhu 2021-02-28
287123c Donghe-Zhu 2021-02-27
54d5b5b Donghe-Zhu 2021-02-27
330f064 Donghe-Zhu 2021-02-27
adbc9bc Donghe-Zhu 2021-02-27
5937141 Donghe-Zhu 2021-02-27
4414bbf Donghe-Zhu 2021-02-27
a265efb Donghe-Zhu 2021-02-27
19edd1e Donghe-Zhu 2021-02-27
f20483f Donghe-Zhu 2021-02-26
6a2c7b3 Donghe-Zhu 2021-02-25
02b976d Donghe-Zhu 2021-02-24
354c224 Donghe-Zhu 2021-02-24
1a0a88a Donghe-Zhu 2021-02-24
57f701e Donghe-Zhu 2021-02-24
06f3149 Donghe-Zhu 2021-02-16
e3bba84 Donghe-Zhu 2021-02-15
5dce4b1 Donghe-Zhu 2021-02-15
4469a0c Donghe-Zhu 2021-02-13
5ae6a69 Donghe-Zhu 2021-02-10
05385dc Donghe-Zhu 2021-02-10
f791ae4 Donghe-Zhu 2021-02-09
f71ae34 Donghe-Zhu 2021-02-09
c011832 Donghe-Zhu 2021-02-09
a145fa7 Donghe-Zhu 2021-02-09
c344e42 Donghe-Zhu 2021-02-08
2f095d7 Donghe-Zhu 2021-02-07
1fad5f1 Donghe-Zhu 2021-02-07
ca03c39 Donghe-Zhu 2021-02-07
cd7c52c Donghe-Zhu 2021-02-04
bcf84f4 Donghe-Zhu 2021-02-02
61666de Donghe-Zhu 2021-01-31
865b582 Donghe-Zhu 2021-01-31
3e68089 Donghe-Zhu 2021-01-31
ecf335c Donghe-Zhu 2021-01-31
a618965 Donghe-Zhu 2021-01-31
59e006e Donghe-Zhu 2021-01-31
a1c8f87 Donghe-Zhu 2021-01-31
ae5c18f Donghe-Zhu 2021-01-31
b50fe52 Donghe-Zhu 2021-01-31
ac99ae5 jens-daniel-mueller 2021-01-29
b5bdcaf Donghe-Zhu 2021-01-29
442010d Donghe-Zhu 2021-01-29
372adf5 Donghe-Zhu 2021-01-29
af8788e Donghe-Zhu 2021-01-29
21c91c9 Donghe-Zhu 2021-01-29
eded038 Donghe-Zhu 2021-01-29
541d4dd Donghe-Zhu 2021-01-29
6a75576 Donghe-Zhu 2021-01-28
16fba40 Donghe-Zhu 2021-01-28
12bc567 Donghe-Zhu 2021-01-27
ceed31b Donghe-Zhu 2021-01-27
342402d Donghe-Zhu 2021-01-27
5bad5c2 Donghe-Zhu 2021-01-27
61efb56 Donghe-Zhu 2021-01-25
48f638e Donghe-Zhu 2021-01-25
c1cec47 Donghe-Zhu 2021-01-25
05ffb0c Donghe-Zhu 2021-01-25
8b97165 Donghe-Zhu 2021-01-25
c569946 Donghe-Zhu 2021-01-24
a2f0d56 Donghe-Zhu 2021-01-23
28509fc Donghe-Zhu 2021-01-23
4c28e4a Donghe-Zhu 2021-01-22
24cc264 jens-daniel-mueller 2021-01-22
7891955 Donghe-Zhu 2021-01-21
d4cf1cb Donghe-Zhu 2021-01-21
1f3e5b6 jens-daniel-mueller 2021-01-20
0e7bdf1 jens-daniel-mueller 2021-01-15
4571843 jens-daniel-mueller 2021-01-14
b3564aa jens-daniel-mueller 2021-01-14
8d032c3 jens-daniel-mueller 2021-01-14
17dee1d jens-daniel-mueller 2021-01-13
7cdea0c jens-daniel-mueller 2021-01-06
fa85b93 jens-daniel-mueller 2021-01-06
e5cb81a Donghe-Zhu 2021-01-05
a499f10 Donghe-Zhu 2021-01-05
fb8a752 Donghe-Zhu 2020-12-23
8fae0b2 Donghe-Zhu 2020-12-21
c8b76b3 jens-daniel-mueller 2020-12-19 
lm_best %>% 
  ggplot(aes(rmse, aic, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  facet_grid(eras~basin)

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

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