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: 05657ed

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 05657ed. 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 db33928 Donghe-Zhu 2021-03-10 Build site.
html 3d1f470 Donghe-Zhu 2021-03-10 Build site.
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, nitrate, silicate, phosphate
# 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:

  • 88

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 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 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 158 1.0538228 478.950654 4.4968455 1982-1999 –> 2000-2012 2.3907671 1204.02994
Atlantic (-Inf,26] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 158 1.0447877 476.229718 3.9913380 1982-1999 –> 2000-2012 2.3809752 1201.07344
Atlantic (-Inf,26] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 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 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 72 1.1990801 244.471022 3.3796520 2000-2012 –> 2013-2019 2.2529029 723.42168
Atlantic (-Inf,26] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 72 1.1767614 243.765461 2.7882798 2000-2012 –> 2013-2019 2.2626619 734.19147
Atlantic (-Inf,26] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 72 1.1875291 245.077119 3.1942008 2000-2012 –> 2013-2019 2.2307669 722.83769
Atlantic (-Inf,26] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 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 + nitrate + silicate 1040 3.2353488 5407.596546 14.5591803 1982-1999 –> 2000-2012 6.4866697 13277.94818
Atlantic (26,26.5] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 1040 3.3140799 5457.606621 14.9918101 1982-1999 –> 2000-2012 6.6603260 13414.98229
Atlantic (26,26.5] 5 cstar_tref ~ sal + temp + oxygen + aou + nitrate 1040 3.4152862 5520.175649 15.2363993 1982-1999 –> 2000-2012 6.8200055 13529.93713
Atlantic (26,26.5] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 1040 3.1872705 5378.455168 13.6598787 1982-1999 –> 2000-2012 6.3907896 13206.01701
Atlantic (26,26.5] 5 cstar_tref ~ sal + temp + oxygen + nitrate + silicate 1040 3.3405660 5474.163895 15.2055707 1982-1999 –> 2000-2012 6.7195852 13461.01249
Atlantic (26,26.5] 5 cstar_tref ~ sal + oxygen + aou + nitrate + silicate 562 3.0957021 2879.023463 9.4138668 2000-2012 –> 2013-2019 6.3310508 8286.62001
Atlantic (26,26.5] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 562 3.1825680 2910.128705 9.4627336 2000-2012 –> 2013-2019 6.4966480 8367.73533
Atlantic (26,26.5] 5 cstar_tref ~ sal + temp + oxygen + aou + nitrate 562 3.3333181 2962.147207 11.1442415 2000-2012 –> 2013-2019 6.7486043 8482.32286
Atlantic (26,26.5] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 562 3.0399853 2860.609307 10.8756225 2000-2012 –> 2013-2019 6.2272558 8239.06447
Atlantic (26,26.5] 5 cstar_tref ~ sal + temp + oxygen + nitrate + silicate 562 3.2090047 2919.426883 9.5659666 2000-2012 –> 2013-2019 6.5495707 8393.59078
Atlantic (26.5,26.75] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 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 1360 3.2964180 7118.027900 10.5109060 1982-1999 –> 2000-2012 6.3652079 16464.80301
Atlantic (26.5,26.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1360 3.2276124 7062.652885 10.1005384 1982-1999 –> 2000-2012 6.2017135 16296.27934
Atlantic (26.5,26.75] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 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 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 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 646 3.3445545 3407.143458 10.7277964 2000-2012 –> 2013-2019 6.6409725 10525.17136
Atlantic (26.5,26.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 646 3.2958735 3390.199829 10.5343935 2000-2012 –> 2013-2019 6.5234860 10452.85271
Atlantic (26.5,26.75] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 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 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 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 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 2134 2.0721249 9179.586005 14.1852345 1982-1999 –> 2000-2012 3.7872525 19874.84163
Atlantic (26.75,27] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 2134 1.9756584 8978.118057 13.5842460 1982-1999 –> 2000-2012 3.4606942 18889.73246
Atlantic (26.75,27] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 2134 2.1557967 9348.538081 19.2601021 1982-1999 –> 2000-2012 3.6479228 19284.13041
Atlantic (26.75,27] 4 cstar_tref ~ sal + aou + silicate + phosphate 1121 2.6063943 5341.023922 24.8434412 2000-2012 –> 2013-2019 4.7269587 14617.23359
Atlantic (26.75,27] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 1121 2.5950428 5333.238180 25.4741707 2000-2012 –> 2013-2019 4.6101370 14393.70999
Atlantic (26.75,27] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 1121 2.6063657 5342.999353 24.8816347 2000-2012 –> 2013-2019 4.6784906 14522.58536
Atlantic (26.75,27] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1121 2.2698743 5035.082409 21.4501194 2000-2012 –> 2013-2019 4.2455326 14013.20047
Atlantic (26.75,27] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 1121 2.3110503 5073.388314 22.0043753 2000-2012 –> 2013-2019 4.4668470 14421.92640
Atlantic (27,27.25] 4 cstar_tref ~ sal + aou + silicate + phosphate 1717 2.0250912 7307.715917 9.3807341 1982-1999 –> 2000-2012 3.7447144 16116.95820
Atlantic (27,27.25] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 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 1717 2.0148260 7292.264615 9.3883458 1982-1999 –> 2000-2012 3.6145946 15779.41047
Atlantic (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1717 1.9602218 7199.914902 13.1630682 1982-1999 –> 2000-2012 3.5598163 15688.57230
Atlantic (27,27.25] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 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 982 2.4264509 4539.743124 13.9625358 2000-2012 –> 2013-2019 4.4515421 11847.45904
Atlantic (27,27.25] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 982 2.4198641 4536.404447 14.2305742 2000-2012 –> 2013-2019 4.4271125 11815.72975
Atlantic (27,27.25] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 982 2.4120191 4530.026992 14.3657953 2000-2012 –> 2013-2019 4.4268451 11822.29161
Atlantic (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 982 2.3584087 4487.881942 19.9575086 2000-2012 –> 2013-2019 4.4550936 11918.90451
Atlantic (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 982 2.3482249 4479.382903 14.3974130 2000-2012 –> 2013-2019 4.3084468 11679.29781
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + oxygen + aou + nitrate + silicate 1621 2.0544117 6948.404835 17.9117945 1982-1999 –> 2000-2012 3.7116240 15616.89592
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 1621 2.0187330 6891.606684 14.5715494 1982-1999 –> 2000-2012 3.6660378 15533.12636
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + temp + oxygen + aou + nitrate 1621 2.0269950 6904.848048 11.4882386 1982-1999 –> 2000-2012 3.7841406 15836.71427
Atlantic (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 1621 2.0072942 6875.184353 10.8817923 1982-1999 –> 2000-2012 3.6512482 15509.54506
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + temp + oxygen + nitrate + silicate 1621 2.0195251 6892.878564 10.5563778 1982-1999 –> 2000-2012 3.7197074 15676.51218
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 904 2.4753224 4218.159019 17.7494577 2000-2012 –> 2013-2019 4.4940554 11109.76570
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + temp + oxygen + aou + nitrate 904 2.4587292 4205.998352 13.0983715 2000-2012 –> 2013-2019 4.4857242 11110.84640
Atlantic (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 904 2.4552021 4205.402918 13.6176582 2000-2012 –> 2013-2019 4.4624964 11080.58727
Atlantic (27.25,27.5] 4 cstar_tref ~ sal + temp + oxygen + nitrate 904 2.4682151 4210.960297 13.0763333 2000-2012 –> 2013-2019 4.5035365 11127.09846
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + temp + oxygen + nitrate + silicate 904 2.4623754 4208.677565 13.7371184 2000-2012 –> 2013-2019 4.4819005 11101.55613
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + temp + oxygen + aou + nitrate 2079 2.1493013 9095.410309 13.2353104 1982-1999 –> 2000-2012 4.0590191 20641.23057
Atlantic (27.5,27.75] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 2079 2.1421629 9083.577489 13.9206011 1982-1999 –> 2000-2012 4.0510182 20628.87626
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 2079 2.3075871 9390.876085 16.0914994 1982-1999 –> 2000-2012 4.4239676 21510.25739
Atlantic (27.5,27.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 2079 2.3067859 9391.432216 16.3811627 1982-1999 –> 2000-2012 4.4091130 21475.62356
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + temp + oxygen + nitrate + silicate 2079 2.3723587 9505.978976 14.8799689 1982-1999 –> 2000-2012 4.4937874 21638.65892
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + temp + oxygen + aou + nitrate 1115 2.5428241 5259.436901 15.6212383 2000-2012 –> 2013-2019 4.6921255 14354.84721
Atlantic (27.5,27.75] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 1115 2.5158530 5237.657523 16.3872161 2000-2012 –> 2013-2019 4.6580160 14321.23501
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 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 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 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 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 722 1.6360802 2773.833134 14.5982784 1982-1999 –> 2000-2012 3.3135078 6497.59228
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 722 1.5992193 2740.927704 12.1339144 1982-1999 –> 2000-2012 3.2778246 6466.03163
Atlantic (27.75,27.85] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 722 1.5916813 2736.105232 13.2839690 1982-1999 –> 2000-2012 3.2452848 6434.45688
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 722 1.6535346 2789.156784 14.7399711 1982-1999 –> 2000-2012 3.3439506 6527.69445
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 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 + oxygen + aou + nitrate 364 1.7150839 1439.715595 11.2024479 2000-2012 –> 2013-2019 3.3442322 4207.41768
Atlantic (27.75,27.85] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 364 1.7133169 1440.965149 11.0768950 2000-2012 –> 2013-2019 3.3409938 4209.36249
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 364 1.6780668 1423.830913 10.6183742 2000-2012 –> 2013-2019 3.2772860 4164.75862
Atlantic (27.75,27.85] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 364 1.6649520 1420.118948 10.6920877 2000-2012 –> 2013-2019 3.2566332 4156.22418
Atlantic (27.85,27.95] 4 cstar_tref ~ oxygen + aou + silicate + phosphate 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 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 + aou + silicate + phosphate 694 2.8772505 3450.373867 20.1802195 1982-1999 –> 2000-2012 5.8300293 8277.61627
Atlantic (27.85,27.95] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 694 2.6858576 3356.830803 20.1208556 1982-1999 –> 2000-2012 5.6386261 8186.06643
Atlantic (27.85,27.95] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 694 2.6878487 3355.859387 20.1028452 1982-1999 –> 2000-2012 5.6412180 8183.48649
Atlantic (27.85,27.95] 4 cstar_tref ~ oxygen + aou + silicate + phosphate 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 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 + aou + silicate + phosphate 415 3.0407759 2114.772546 18.0277960 2000-2012 –> 2013-2019 5.9180264 5565.14641
Atlantic (27.85,27.95] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 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 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 723 3.1275877 3714.603974 22.5137593 1982-1999 –> 2000-2012 6.8809068 9173.38001
Atlantic (27.95,28.05] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 723 3.0495655 3680.073829 12.6814197 1982-1999 –> 2000-2012 6.7958653 9137.13227
Atlantic (27.95,28.05] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 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 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 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 386 2.8428108 1916.000961 14.0211453 2000-2012 –> 2013-2019 5.9703985 5630.60493
Atlantic (27.95,28.05] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 386 2.8070603 1908.230904 12.9707292 2000-2012 –> 2013-2019 5.8566257 5588.30473
Atlantic (27.95,28.05] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 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 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 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 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 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 436 0.2175423 -78.801383 1.4425099 1982-1999 –> 2000-2012 0.3337872 -940.14322
Atlantic (28.05,28.1] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 436 0.1931541 -180.486583 1.6293925 1982-1999 –> 2000-2012 0.3078192 -1056.16703
Atlantic (28.05,28.1] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 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 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 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 + aou + silicate + phosphate 229 0.2552708 38.506822 1.1720676 2000-2012 –> 2013-2019 0.4728131 -40.29456
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 229 0.2362004 2.945556 1.0819569 2000-2012 –> 2013-2019 0.4378010 -142.21867
Atlantic (28.05,28.1] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 229 0.2182829 -31.185279 0.9320609 2000-2012 –> 2013-2019 0.4114370 -211.67186
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 541 0.1394443 -582.345612 1.0824605 1982-1999 –> 2000-2012 0.2743304 -1472.97554
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 541 0.1413974 -567.296224 1.1077546 1982-1999 –> 2000-2012 0.2772035 -1447.40398
Atlantic (28.1,28.15] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 541 0.1354765 -611.579729 0.9957761 1982-1999 –> 2000-2012 0.2695767 -1509.25535
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 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 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 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 288 0.1305243 -341.540085 0.4627303 2000-2012 –> 2013-2019 0.2719217 -908.83631
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 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 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 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 + nitrate + silicate 1155 0.1374276 -1292.812100 0.8111592 1982-1999 –> 2000-2012 0.3180622 -2292.06363
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 1155 0.1372783 -1295.322840 0.8080558 1982-1999 –> 2000-2012 0.3178075 -2296.59800
Atlantic (28.15,28.2] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 1155 0.1354897 -1323.618192 0.7466368 1982-1999 –> 2000-2012 0.3145638 -2350.94297
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + temp + oxygen + nitrate + silicate 1155 0.1376435 -1289.186490 0.8138382 1982-1999 –> 2000-2012 0.3185413 -2283.39147
Atlantic (28.15,28.2] 5 cstar_tref ~ temp + oxygen + aou + nitrate + silicate 1155 0.1369421 -1300.986910 0.7217059 1982-1999 –> 2000-2012 0.3188855 -2275.21659
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + oxygen + aou + nitrate + silicate 636 0.1898096 -294.835912 1.4147466 2000-2012 –> 2013-2019 0.3272372 -1587.64801
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 636 0.1898163 -294.790951 1.4151881 2000-2012 –> 2013-2019 0.3270946 -1590.11379
Atlantic (28.15,28.2] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 636 0.1896121 -294.159817 1.3775091 2000-2012 –> 2013-2019 0.3251018 -1617.77801
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + temp + oxygen + nitrate + silicate 636 0.1897111 -295.495890 1.4036411 2000-2012 –> 2013-2019 0.3273546 -1584.68238
Atlantic (28.15,28.2] 5 cstar_tref ~ temp + oxygen + aou + nitrate + silicate 636 0.2008104 -223.171584 1.7076436 2000-2012 –> 2013-2019 0.3377525 -1524.15849
Atlantic (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + silicate 955 0.1744896 -612.477534 0.9944731 1982-1999 –> 2000-2012 0.4081182 -687.41598
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + oxygen + aou + nitrate + silicate 955 0.1695779 -665.013678 0.9467358 1982-1999 –> 2000-2012 0.4007085 -764.58995
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 955 0.1698189 -662.301029 0.9333870 1982-1999 –> 2000-2012 0.4014157 -756.88416
Atlantic (28.2, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 955 0.1683740 -676.621511 1.1006316 1982-1999 –> 2000-2012 0.3865669 -916.93887
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + temp + oxygen + nitrate + silicate 955 0.1708923 -650.265552 0.8843184 1982-1999 –> 2000-2012 0.3993587 -778.57131
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + oxygen + aou + nitrate + silicate 552 0.2329105 -28.131217 1.5004320 2000-2012 –> 2013-2019 0.4024884 -693.14489
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 552 0.2324457 -30.336709 1.4787637 2000-2012 –> 2013-2019 0.4022646 -692.63774
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + temp + oxygen + aou + nitrate 552 0.2354907 -15.968438 1.1910259 2000-2012 –> 2013-2019 0.4039971 -693.08900
Atlantic (28.2, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 552 0.2241162 -68.624020 1.3072927 2000-2012 –> 2013-2019 0.3924902 -745.24553
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + temp + oxygen + nitrate + silicate 552 0.2299633 -42.190459 1.4503171 2000-2012 –> 2013-2019 0.4008556 -692.45601
Indo-Pacific (-Inf,26] 4 cstar_tref ~ sal + oxygen + aou + phosphate 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 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 5288 6.3160008 34513.173890 29.8497761 1982-1999 –> 2000-2012 12.9040860 82048.19652
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 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 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 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 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 2815 5.9886169 18079.538454 24.9064180 2000-2012 –> 2013-2019 12.3046177 52592.71234
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 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 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 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 5428 4.5972029 31978.260876 39.9466263 1982-1999 –> 2000-2012 8.8916713 75688.49856
Indo-Pacific (26,26.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 5428 4.5340876 31830.185580 43.0233728 1982-1999 –> 2000-2012 8.7741402 75348.69317
Indo-Pacific (26,26.5] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 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 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 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 2876 4.5713284 16917.646172 39.7335019 2000-2012 –> 2013-2019 9.1685314 48895.90705
Indo-Pacific (26,26.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 2876 4.4815444 16805.548776 46.3106843 2000-2012 –> 2013-2019 9.0156321 48635.73436
Indo-Pacific (26,26.5] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 2876 4.6557423 17022.893429 36.3745282 2000-2012 –> 2013-2019 9.3343493 49191.70271
Indo-Pacific (26,26.5] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 2876 4.8127496 17213.671221 38.8375973 2000-2012 –> 2013-2019 9.6831832 49818.70223
Indo-Pacific (26.5,26.75] 4 cstar_tref ~ oxygen + aou + silicate + phosphate 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 4269 3.9438120 23844.294920 20.4830181 1982-1999 –> 2000-2012 7.6993981 56369.61516
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 4269 4.0637389 24100.056558 20.7793361 1982-1999 –> 2000-2012 7.9220345 56945.26579
Indo-Pacific (26.5,26.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 4269 3.8963414 23742.901689 26.7023232 1982-1999 –> 2000-2012 7.6360524 56219.99995
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 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 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 2324 4.3199526 13410.386398 24.6916724 2000-2012 –> 2013-2019 8.2637646 37254.68132
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 2324 4.4417007 13539.567831 24.4985770 2000-2012 –> 2013-2019 8.5054396 37639.62439
Indo-Pacific (26.5,26.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 2324 4.2544954 13341.419522 23.3613336 2000-2012 –> 2013-2019 8.1508368 37084.32121
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 2324 4.3137865 13403.747333 23.1641121 2000-2012 –> 2013-2019 8.2596379 37252.45617
Indo-Pacific (26.75,27] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 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 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 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 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 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 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 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 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 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 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 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 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 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 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 4267 3.4041802 22577.406832 30.9255112 1982-1999 –> 2000-2012 6.2015505 50822.23662
Indo-Pacific (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 2344 4.7398557 13962.543204 41.3080758 2000-2012 –> 2013-2019 9.1373896 38726.99338
Indo-Pacific (27,27.25] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 2344 4.5554949 13774.558428 35.2893459 2000-2012 –> 2013-2019 8.2904863 37143.42005
Indo-Pacific (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 2344 3.9472749 13104.727079 23.9342640 2000-2012 –> 2013-2019 7.3509079 35682.76198
Indo-Pacific (27,27.25] 4 cstar_tref ~ temp + oxygen + aou + phosphate 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 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 3520 3.1346725 18046.701221 29.1767786 1982-1999 –> 2000-2012 6.0487445 41537.11519
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 3520 3.1202129 18014.152053 29.4518972 1982-1999 –> 2000-2012 6.0220192 41464.77333
Indo-Pacific (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 3520 2.8617081 17407.314891 30.1089756 1982-1999 –> 2000-2012 5.6400680 40449.81559
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 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 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 1879 3.6867091 10249.562204 27.5141850 2000-2012 –> 2013-2019 6.8213816 28296.26342
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 1879 3.6579924 10220.175560 27.9920522 2000-2012 –> 2013-2019 6.7782053 28234.32761
Indo-Pacific (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1879 3.3074238 9843.576259 29.1273995 2000-2012 –> 2013-2019 6.1691319 27250.89115
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 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 1879 3.3526141 9892.575280 28.9794819 2000-2012 –> 2013-2019 6.2300248 27336.41416
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 3170 1.8834035 13023.800742 9.5549910 1982-1999 –> 2000-2012 3.5338959 30329.43049
Indo-Pacific (27.5,27.75] 4 cstar_tref ~ sal + oxygen + silicate + phosphate 3170 1.9247016 13159.318110 14.1925496 1982-1999 –> 2000-2012 3.5760715 30467.73487
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 3170 1.8808520 13015.205979 9.5198488 1982-1999 –> 2000-2012 3.5317875 30323.25341
Indo-Pacific (27.5,27.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 3170 1.8800469 13014.491400 9.6798100 1982-1999 –> 2000-2012 3.5305329 30322.08642
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 3170 1.8817537 13018.244743 9.5293548 1982-1999 –> 2000-2012 3.5322649 30323.97696
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 1683 2.1483716 7364.161462 12.8217994 2000-2012 –> 2013-2019 4.0317751 20387.96220
Indo-Pacific (27.5,27.75] 4 cstar_tref ~ sal + oxygen + silicate + phosphate 1683 2.2955622 7585.218276 14.9967179 2000-2012 –> 2013-2019 4.2202638 20744.53639
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 1683 2.1449046 7358.725091 12.7298567 2000-2012 –> 2013-2019 4.0257566 20373.93107
Indo-Pacific (27.5,27.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1683 2.1441647 7359.563742 12.6747027 2000-2012 –> 2013-2019 4.0242115 20374.05514
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 1683 2.1454909 7359.645130 12.7476542 2000-2012 –> 2013-2019 4.0272446 20377.88987
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 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 + oxygen + aou + phosphate 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 1216 1.1484687 3803.519113 5.5653347 1982-1999 –> 2000-2012 2.2742935 9009.96521
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 1216 1.4365283 4345.800244 8.0082261 1982-1999 –> 2000-2012 2.7661797 10112.01580
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 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 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 + oxygen + aou + phosphate 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 691 1.3067937 2346.763836 7.0585944 2000-2012 –> 2013-2019 2.4552624 6150.28295
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 691 1.6200931 2641.765389 13.6362004 2000-2012 –> 2013-2019 3.0566214 6987.56563
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 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 + nitrate + silicate 1797 1.3367652 6156.833290 9.1081321 1982-1999 –> 2000-2012 2.5868821 14483.52596
Indo-Pacific (27.85,27.95] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 1797 1.2314246 5863.834129 6.7020452 1982-1999 –> 2000-2012 2.3640740 13693.02108
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 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 1797 1.1268453 5544.867661 6.8163957 1982-1999 –> 2000-2012 2.0904038 12555.62250
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ sal + oxygen + aou + nitrate + silicate 987 1.5418683 3669.716560 8.7959041 2000-2012 –> 2013-2019 2.8786336 9826.54985
Indo-Pacific (27.85,27.95] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 987 1.4141378 3501.015168 6.9780896 2000-2012 –> 2013-2019 2.6455624 9364.84930
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 987 1.4361554 3529.512839 8.9979380 2000-2012 –> 2013-2019 2.7232012 9550.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 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 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 1162 0.5807005 2048.472523 3.1107181 1982-1999 –> 2000-2012 1.1238891 4874.92133
Indo-Pacific (27.95,28.05] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1162 0.5749179 2027.213849 2.9262058 1982-1999 –> 2000-2012 1.1061869 4778.49345
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 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 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 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 648 0.6059601 1203.724657 2.6446851 2000-2012 –> 2013-2019 1.1866607 3252.19718
Indo-Pacific (27.95,28.05] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 648 0.5887663 1168.419405 2.7201223 2000-2012 –> 2013-2019 1.1636841 3195.63325
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 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 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 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 1015 0.4568373 1304.086554 2.4288325 1982-1999 –> 2000-2012 0.9129949 3244.24376
Indo-Pacific (28.05,28.1] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1015 0.4542785 1294.684324 2.4180424 1982-1999 –> 2000-2012 0.9081081 3221.29740
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 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 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 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 565 0.4757604 777.990361 2.6230255 2000-2012 –> 2013-2019 0.9325977 2082.07691
Indo-Pacific (28.05,28.1] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 565 0.4745332 777.071782 2.6197171 2000-2012 –> 2013-2019 0.9288117 2071.75611
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 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 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 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 9318 0.7771207 21758.092479 4.8370707 1982-1999 –> 2000-2012 1.5392513 51546.85223
Indo-Pacific (28.1, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 9318 0.7739529 21683.971466 4.9061100 1982-1999 –> 2000-2012 1.5356708 51460.67513
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 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 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 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 5093 0.7773139 11901.342823 4.2794415 2000-2012 –> 2013-2019 1.5544347 33659.43530
Indo-Pacific (28.1, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 5093 0.7692293 11796.846452 4.3310811 2000-2012 –> 2013-2019 1.5431823 33480.81792
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 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 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 103 2.7385358 513.830562 11.4817285 1982-1999 –> 2000-2012 6.2369968 1217.01608
SO (26,26.5] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 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 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 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 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 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 47 3.0856619 253.296243 8.0888582 2000-2012 –> 2013-2019 5.8241977 767.12680
SO (26,26.5] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 47 3.1990190 256.687576 7.9977611 2000-2012 –> 2013-2019 5.9461226 771.16162
SO (26,26.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 47 3.0856552 255.296039 8.0922982 2000-2012 –> 2013-2019 5.7584324 766.11970
SO (26,26.5] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 47 3.0875443 253.353570 8.0376410 2000-2012 –> 2013-2019 5.8122599 766.14191
SO (26.5,26.75] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 312 3.2153502 1628.209882 10.5226323 1982-1999 –> 2000-2012 6.7106987 3885.32203
SO (26.5,26.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 312 3.1290661 1613.236019 9.9544018 1982-1999 –> 2000-2012 6.6028383 3867.14691
SO (26.5,26.75] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 312 3.2879009 1642.133236 10.0473415 1982-1999 –> 2000-2012 6.9083141 3928.77566
SO (26.5,26.75] 4 cstar_tref ~ temp + oxygen + aou + phosphate 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 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 133 3.2950584 708.622405 10.9719273 2000-2012 –> 2013-2019 6.6038186 2354.70195
SO (26.5,26.75] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 133 3.3899758 716.176510 10.0861110 2000-2012 –> 2013-2019 6.6053260 2344.38639
SO (26.5,26.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 133 3.2742828 708.939942 10.4020030 2000-2012 –> 2013-2019 6.4033488 2322.17596
SO (26.5,26.75] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 133 3.2903253 708.240039 10.8891344 2000-2012 –> 2013-2019 6.5782262 2350.37328
SO (26.5,26.75] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 133 3.3148729 710.217175 10.3205968 2000-2012 –> 2013-2019 6.5618670 2344.53815
SO (26.75,27] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 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 1061 2.2880755 4783.390414 9.3213970 1982-1999 –> 2000-2012 4.7000755 11766.55673
SO (26.75,27] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 1061 2.8433092 5242.410902 9.0743310 1982-1999 –> 2000-2012 5.8088340 12849.56830
SO (26.75,27] 4 cstar_tref ~ temp + oxygen + aou + phosphate 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 1061 2.3191394 4810.005845 10.6756905 1982-1999 –> 2000-2012 4.7542973 11820.12478
SO (26.75,27] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 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 550 2.2521270 2469.895006 11.5622287 2000-2012 –> 2013-2019 4.5402025 7253.28542
SO (26.75,27] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 550 2.8956977 2744.381089 9.1121669 2000-2012 –> 2013-2019 5.7390069 7986.79199
SO (26.75,27] 4 cstar_tref ~ temp + oxygen + aou + phosphate 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 550 2.2748295 2478.928005 12.6508529 2000-2012 –> 2013-2019 4.5939689 7288.93385
SO (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 2513 1.8266872 10175.770410 7.3417975 1982-1999 –> 2000-2012 3.5919249 23855.92174
SO (27,27.25] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 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 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 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 2513 1.6615541 9697.553575 16.3323879 1982-1999 –> 2000-2012 3.2687331 22730.70242
SO (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 1286 1.9225746 5346.736882 7.0614429 2000-2012 –> 2013-2019 3.7492618 15522.50729
SO (27,27.25] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 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 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 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 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 + nitrate + silicate 2731 1.6653543 10550.069231 7.6152449 1982-1999 –> 2000-2012 3.0994883 24102.50297
SO (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 2731 1.6617833 10540.344591 7.3898985 1982-1999 –> 2000-2012 3.0863767 24043.99706
SO (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 2731 1.6887194 10628.169039 10.6094178 1982-1999 –> 2000-2012 3.1238256 24187.75864
SO (27.25,27.5] 5 cstar_tref ~ sal + temp + oxygen + nitrate + silicate 2731 1.6623528 10540.215983 7.4913929 1982-1999 –> 2000-2012 3.0903673 24060.11706
SO (27.25,27.5] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 2731 1.6893569 10628.230775 10.5492429 1982-1999 –> 2000-2012 3.1276338 24202.61058
SO (27.25,27.5] 5 cstar_tref ~ sal + oxygen + aou + nitrate + silicate 1453 1.8281321 5890.609837 8.3046413 2000-2012 –> 2013-2019 3.5085857 16489.97855
SO (27.25,27.5] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 1453 1.8209928 5879.239046 8.3707548 2000-2012 –> 2013-2019 3.4863471 16429.30828
SO (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 1453 1.8209112 5881.108844 8.3718926 2000-2012 –> 2013-2019 3.4826945 16421.45343
SO (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 1453 1.8445309 5918.561145 8.9141838 2000-2012 –> 2013-2019 3.5332503 16546.73018
SO (27.25,27.5] 5 cstar_tref ~ sal + temp + oxygen + nitrate + silicate 1453 1.8228962 5882.274894 8.3460548 2000-2012 –> 2013-2019 3.4852489 16422.49088
SO (27.5,27.75] 5 cstar_tref ~ sal + oxygen + aou + nitrate + silicate 3391 1.4915886 12348.967798 8.4899800 1982-1999 –> 2000-2012 2.7507785 27356.07857
SO (27.5,27.75] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 3391 1.4899251 12341.399836 8.4595019 1982-1999 –> 2000-2012 2.7459402 27325.56271
SO (27.5,27.75] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 3391 1.4885691 12337.224464 8.4034326 1982-1999 –> 2000-2012 2.7436284 27316.46812
SO (27.5,27.75] 5 cstar_tref ~ sal + temp + oxygen + nitrate + silicate 3391 1.4885692 12335.224989 8.4029779 1982-1999 –> 2000-2012 2.7438968 27314.41162
SO (27.5,27.75] 5 cstar_tref ~ temp + oxygen + aou + nitrate + silicate 3391 1.5313544 12527.407763 8.9657450 1982-1999 –> 2000-2012 2.7893155 27525.64374
SO (27.5,27.75] 4 cstar_tref ~ sal + aou + nitrate + silicate 1821 1.7135136 7141.158593 9.9489584 2000-2012 –> 2013-2019 3.2205821 19558.14808
SO (27.5,27.75] 5 cstar_tref ~ sal + oxygen + aou + nitrate + silicate 1821 1.7088024 7133.131447 10.0111865 2000-2012 –> 2013-2019 3.2003910 19482.09925
SO (27.5,27.75] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 1821 1.7084237 7132.324177 9.9997974 2000-2012 –> 2013-2019 3.1983488 19473.72401
SO (27.5,27.75] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 1821 1.7082430 7133.938865 9.9811291 2000-2012 –> 2013-2019 3.1968120 19471.16333
SO (27.5,27.75] 5 cstar_tref ~ sal + temp + oxygen + nitrate + silicate 1821 1.7085924 7132.683905 9.9566587 2000-2012 –> 2013-2019 3.1971616 19467.90889
SO (27.75,27.85] 5 cstar_tref ~ sal + oxygen + aou + nitrate + silicate 1415 1.4208930 5023.734079 11.3449694 1982-1999 –> 2000-2012 2.6760780 10980.59382
SO (27.75,27.85] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 1415 1.4221268 5026.190477 11.3474457 1982-1999 –> 2000-2012 2.6814715 10994.99388
SO (27.75,27.85] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 1415 1.4184928 5020.949559 11.3324638 1982-1999 –> 2000-2012 2.6688137 10965.79268
SO (27.75,27.85] 5 cstar_tref ~ sal + temp + oxygen + nitrate + silicate 1415 1.4249706 5031.843850 11.3778380 1982-1999 –> 2000-2012 2.6893537 11015.06147
SO (27.75,27.85] 5 cstar_tref ~ temp + oxygen + aou + nitrate + silicate 1415 1.4767020 5132.761719 12.9090205 1982-1999 –> 2000-2012 2.7281763 11078.93338
SO (27.75,27.85] 4 cstar_tref ~ sal + aou + nitrate + silicate 722 1.5764019 2718.176570 8.5554720 2000-2012 –> 2013-2019 3.0324771 7809.12992
SO (27.75,27.85] 5 cstar_tref ~ sal + oxygen + aou + nitrate + silicate 722 1.5447113 2690.851952 8.3631592 2000-2012 –> 2013-2019 2.9656043 7714.58603
SO (27.75,27.85] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 722 1.5451324 2691.245552 8.3303789 2000-2012 –> 2013-2019 2.9672593 7717.43603
SO (27.75,27.85] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 722 1.5445744 2692.723951 8.3971448 2000-2012 –> 2013-2019 2.9630672 7713.67351
SO (27.75,27.85] 5 cstar_tref ~ sal + temp + oxygen + nitrate + silicate 722 1.5462635 2692.302152 8.2814577 2000-2012 –> 2013-2019 2.9712340 7724.14600
SO (27.85,27.95] 5 cstar_tref ~ sal + oxygen + aou + nitrate + silicate 1399 1.4344988 4993.751876 9.0306421 1982-1999 –> 2000-2012 2.8289679 11498.63400
SO (27.85,27.95] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 1399 1.2870785 4692.335133 8.5164180 1982-1999 –> 2000-2012 2.5480100 10826.15996
SO (27.85,27.95] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 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 1399 1.2836451 4684.861093 8.5215752 1982-1999 –> 2000-2012 2.5448741 10819.56053
SO (27.85,27.95] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 709 1.4977451 2600.871146 8.8681972 2000-2012 –> 2013-2019 2.7848237 7293.20628
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 709 1.4880087 2591.623013 9.0787780 2000-2012 –> 2013-2019 2.7716538 7276.48411
SO (27.85,27.95] 5 cstar_tref ~ temp + oxygen + aou + nitrate + silicate 709 1.5376429 2636.150288 8.7215503 2000-2012 –> 2013-2019 2.9798013 7644.80250
SO (27.85,27.95] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 709 1.5409988 2639.241700 8.7822837 2000-2012 –> 2013-2019 2.9762393 7634.43986
SO (27.95,28.05] 5 cstar_tref ~ sal + temp + oxygen + aou + nitrate 1786 1.5469380 6640.831597 8.9903259 1982-1999 –> 2000-2012 2.8545132 15046.53171
SO (27.95,28.05] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 1786 1.5258734 6593.857755 8.3812420 1982-1999 –> 2000-2012 2.8092581 14908.67514
SO (27.95,28.05] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 1786 1.5620816 6675.629386 9.1134438 1982-1999 –> 2000-2012 2.8752317 15102.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 1786 1.5465396 6641.911675 8.3764032 1982-1999 –> 2000-2012 2.8450920 15015.16971
SO (27.95,28.05] 5 cstar_tref ~ sal + temp + oxygen + aou + nitrate 986 1.7028887 3861.893740 8.3557835 2000-2012 –> 2013-2019 3.2498267 10502.72534
SO (27.95,28.05] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 986 1.6261350 3772.945116 8.0450596 2000-2012 –> 2013-2019 3.1520085 10366.80287
SO (27.95,28.05] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 986 1.6936583 3851.175591 8.4129789 2000-2012 –> 2013-2019 3.2557399 10526.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 986 1.6280632 3775.281999 8.1157057 2000-2012 –> 2013-2019 3.1746028 10417.19367
SO (28.05,28.1] 5 cstar_tref ~ sal + temp + oxygen + aou + nitrate 1071 1.0450829 3147.820495 8.4427617 1982-1999 –> 2000-2012 1.9983465 7313.43451
SO (28.05,28.1] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 1071 1.0389462 3137.205569 8.6555860 1982-1999 –> 2000-2012 1.9898140 7297.19998
SO (28.05,28.1] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 1071 1.0487379 3155.298539 9.0027805 1982-1999 –> 2000-2012 2.0069805 7336.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 1071 1.0473850 3154.533537 9.0409967 1982-1999 –> 2000-2012 2.0055492 7337.67458
SO (28.05,28.1] 5 cstar_tref ~ sal + temp + oxygen + aou + nitrate 606 1.2264599 1981.161326 9.4690189 2000-2012 –> 2013-2019 2.2715428 5128.98182
SO (28.05,28.1] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 606 1.2236721 1980.403317 9.7409069 2000-2012 –> 2013-2019 2.2626183 5117.60889
SO (28.05,28.1] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 606 1.2283070 1982.985334 10.1816140 2000-2012 –> 2013-2019 2.2770449 5138.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 606 1.2282847 1984.963330 10.1968303 2000-2012 –> 2013-2019 2.2756697 5139.49687
SO (28.1, Inf] 5 cstar_tref ~ sal + oxygen + aou + nitrate + silicate 5309 1.0911988 16006.995295 6.2254904 1982-1999 –> 2000-2012 2.0608315 36451.06021
SO (28.1, Inf] 4 cstar_tref ~ sal + oxygen + nitrate + silicate 5309 1.0921155 16013.911981 6.1402836 1982-1999 –> 2000-2012 2.0772669 36689.66808
SO (28.1, Inf] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 5309 1.0913601 16008.564990 6.2232166 1982-1999 –> 2000-2012 2.0610962 36454.19956
SO (28.1, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 5309 1.0896709 15994.117646 6.2333995 1982-1999 –> 2000-2012 2.0589958 36435.50978
SO (28.1, Inf] 5 cstar_tref ~ sal + temp + oxygen + nitrate + silicate 5309 1.0912471 16007.465144 6.2230771 1982-1999 –> 2000-2012 2.0609801 36453.05279
SO (28.1, Inf] 5 cstar_tref ~ sal + oxygen + aou + nitrate + silicate 2711 1.1826059 8616.864791 6.3094698 2000-2012 –> 2013-2019 2.2738047 24623.86009
SO (28.1, Inf] 4 cstar_tref ~ sal + oxygen + nitrate + silicate 2711 1.1836919 8619.841560 6.3248519 2000-2012 –> 2013-2019 2.2758074 24633.75354
SO (28.1, Inf] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 2711 1.1826938 8617.267740 6.3100661 2000-2012 –> 2013-2019 2.2740539 24625.83273
SO (28.1, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 2711 1.1818836 8615.551963 6.2901289 2000-2012 –> 2013-2019 2.2715545 24609.66961
SO (28.1, Inf] 5 cstar_tref ~ sal + temp + oxygen + nitrate + silicate 2711 1.1825727 8616.712428 6.3087855 2000-2012 –> 2013-2019 2.2738197 24624.17757
SO (28.1,28.15] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 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 395 0.4707828 539.808215 1.7513525 1982-1999 –> 2000-2012 0.9131138 1182.38618
SO (28.1,28.15] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 395 0.4483007 503.151413 1.3556998 1982-1999 –> 2000-2012 0.8827308 1128.94896
SO (28.1,28.15] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 395 0.4688350 536.532945 1.2066978 1982-1999 –> 2000-2012 0.8955957 1141.77071
SO (28.1,28.15] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 395 0.4412700 490.663828 1.3656130 1982-1999 –> 2000-2012 0.8604524 1079.23170
SO (28.1,28.15] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 209 0.5007130 317.976471 2.3778112 2000-2012 –> 2013-2019 0.9699983 855.26781
SO (28.1,28.15] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 209 0.4655560 289.545781 1.9653861 2000-2012 –> 2013-2019 0.9138567 792.69719
SO (28.1,28.15] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 209 0.4863454 305.806779 1.9335059 2000-2012 –> 2013-2019 0.9551804 842.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 209 0.4627191 286.990868 2.0315187 2000-2012 –> 2013-2019 0.9039892 777.65470
SO (28.15,28.2] 5 cstar_tref ~ sal + oxygen + aou + silicate + phosphate 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 385 0.3068525 196.913871 1.0656804 1982-1999 –> 2000-2012 0.6144770 446.65572
SO (28.15,28.2] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 385 0.2958284 168.741280 1.2331036 1982-1999 –> 2000-2012 0.5978993 400.59322
SO (28.15,28.2] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 385 0.2908445 157.658441 1.0402887 1982-1999 –> 2000-2012 0.5915367 387.01785
SO (28.15,28.2] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate 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 183 0.3423788 141.038937 0.8031784 2000-2012 –> 2013-2019 0.6484427 335.97138
SO (28.15,28.2] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 183 0.3255062 124.542672 0.7719731 2000-2012 –> 2013-2019 0.6272571 310.54715
SO (28.15,28.2] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate 183 0.3300672 127.635476 0.8682035 2000-2012 –> 2013-2019 0.6258955 296.37676
SO (28.15,28.2] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate 183 0.3212016 119.670270 0.7924252 2000-2012 –> 2013-2019 0.6120461 277.32871
SO (28.15,28.2] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate 183 0.3337812 131.730824 0.8025278 2000-2012 –> 2013-2019 0.6459444 341.85710
SO (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + silicate 1738 0.3672295 1462.083841 1.9673041 1982-1999 –> 2000-2012 0.6709769 2598.24970
SO (28.2, Inf] 5 cstar_tref ~ sal + oxygen + aou + nitrate + silicate 1738 0.3666071 1458.187160 1.8944599 1982-1999 –> 2000-2012 0.6703163 2595.73028
SO (28.2, Inf] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 1738 0.3665086 1457.253822 1.8900397 1982-1999 –> 2000-2012 0.6701913 2594.36525
SO (28.2, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 1738 0.3508194 1307.177393 1.8206271 1982-1999 –> 2000-2012 0.6420929 2240.02105
SO (28.2, Inf] 5 cstar_tref ~ sal + temp + oxygen + nitrate + silicate 1738 0.3651559 1444.400770 1.8866124 1982-1999 –> 2000-2012 0.6674325 2558.56726
SO (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + silicate 1005 0.3952080 998.096683 2.5161713 2000-2012 –> 2013-2019 0.7624375 2460.18052
SO (28.2, Inf] 5 cstar_tref ~ sal + oxygen + aou + nitrate + silicate 1005 0.3946691 997.354221 2.5805895 2000-2012 –> 2013-2019 0.7612762 2455.54138
SO (28.2, Inf] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 1005 0.3947488 997.759814 2.5756848 2000-2012 –> 2013-2019 0.7612574 2455.01364
SO (28.2, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 1005 0.3808359 927.639329 2.5626963 2000-2012 –> 2013-2019 0.7316554 2234.81672
SO (28.2, Inf] 5 cstar_tref ~ sal + temp + oxygen + nitrate + silicate 1005 0.3935678 991.737484 2.5748446 2000-2012 –> 2013-2019 0.7587237 2436.13825

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.3839867 1201.9807
Atlantic (-Inf,26] 2000-2012 –> 2013-2019 2.2495187 725.2599
Atlantic (26,26.5] 1982-1999 –> 2000-2012 6.6154752 13377.9794
Atlantic (26,26.5] 2000-2012 –> 2013-2019 6.4706259 8353.8667
Atlantic (26.5,26.75] 1982-1999 –> 2000-2012 6.3080465 16404.8404
Atlantic (26.5,26.75] 2000-2012 –> 2013-2019 6.6070583 10502.5059
Atlantic (26.75,27] 1982-1999 –> 2000-2012 3.6859779 19524.1842
Atlantic (26.75,27] 2000-2012 –> 2013-2019 4.5455932 14393.7312
Atlantic (27,27.25] 1982-1999 –> 2000-2012 3.6441909 15865.0674
Atlantic (27,27.25] 2000-2012 –> 2013-2019 4.4138080 11816.7365
Atlantic (27.25,27.5] 1982-1999 –> 2000-2012 3.7065516 15634.5588
Atlantic (27.25,27.5] 2000-2012 –> 2013-2019 4.4855426 11105.9708
Atlantic (27.5,27.75] 1982-1999 –> 2000-2012 4.2873811 21178.9293
Atlantic (27.5,27.75] 2000-2012 –> 2013-2019 4.8264336 14574.4593
Atlantic (27.75,27.85] 1982-1999 –> 2000-2012 3.2954246 6481.2759
Atlantic (27.75,27.85] 2000-2012 –> 2013-2019 3.3145474 4189.6065
Atlantic (27.85,27.95] 1982-1999 –> 2000-2012 5.7268448 8226.5940
Atlantic (27.85,27.95] 2000-2012 –> 2013-2019 5.8062728 5511.2311
Atlantic (27.95,28.05] 1982-1999 –> 2000-2012 6.8529001 9160.1322
Atlantic (27.95,28.05] 2000-2012 –> 2013-2019 5.9329860 5615.5119
Atlantic (28.05,28.1] 1982-1999 –> 2000-2012 0.3268920 -974.0309
Atlantic (28.05,28.1] 2000-2012 –> 2013-2019 0.4484627 -110.4223
Atlantic (28.1,28.15] 1982-1999 –> 2000-2012 0.2803841 -1419.1502
Atlantic (28.1,28.15] 2000-2012 –> 2013-2019 0.2728189 -904.9017
Atlantic (28.15,28.2] 1982-1999 –> 2000-2012 0.3175721 -2299.6425
Atlantic (28.15,28.2] 2000-2012 –> 2013-2019 0.3289081 -1580.8761
Atlantic (28.2, Inf] 1982-1999 –> 2000-2012 0.3992336 -780.8801
Atlantic (28.2, Inf] 2000-2012 –> 2013-2019 0.4004192 -703.3146
Indo-Pacific (-Inf,26] 1982-1999 –> 2000-2012 12.7547190 81761.8787
Indo-Pacific (-Inf,26] 2000-2012 –> 2013-2019 12.1633766 52397.7109
Indo-Pacific (26,26.5] 1982-1999 –> 2000-2012 8.9932512 75980.8040
Indo-Pacific (26,26.5] 2000-2012 –> 2013-2019 9.2444026 49036.9874
Indo-Pacific (26.5,26.75] 1982-1999 –> 2000-2012 7.8139551 56688.0807
Indo-Pacific (26.5,26.75] 2000-2012 –> 2013-2019 8.3319414 37366.0884
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.0861844 36877.7863
Indo-Pacific (27.25,27.5] 1982-1999 –> 2000-2012 5.8937909 41127.9330
Indo-Pacific (27.25,27.5] 2000-2012 –> 2013-2019 6.5660867 27887.8141
Indo-Pacific (27.5,27.75] 1982-1999 –> 2000-2012 3.5409105 30353.2964
Indo-Pacific (27.5,27.75] 2000-2012 –> 2013-2019 4.0658503 20451.6749
Indo-Pacific (27.75,27.85] 1982-1999 –> 2000-2012 2.5561318 9653.7398
Indo-Pacific (27.75,27.85] 2000-2012 –> 2013-2019 2.7918277 6625.7747
Indo-Pacific (27.85,27.95] 1982-1999 –> 2000-2012 2.3598156 13627.6634
Indo-Pacific (27.85,27.95] 2000-2012 –> 2013-2019 2.6687072 9405.1626
Indo-Pacific (27.95,28.05] 1982-1999 –> 2000-2012 1.1170991 4837.3706
Indo-Pacific (27.95,28.05] 2000-2012 –> 2013-2019 1.1782283 3230.4877
Indo-Pacific (28.05,28.1] 1982-1999 –> 2000-2012 0.9114053 3236.4769
Indo-Pacific (28.05,28.1] 2000-2012 –> 2013-2019 0.9314292 2078.1849
Indo-Pacific (28.1, Inf] 1982-1999 –> 2000-2012 1.5389779 51550.2437
Indo-Pacific (28.1, Inf] 2000-2012 –> 2013-2019 1.5499468 33588.1122
SO (26,26.5] 1982-1999 –> 2000-2012 6.2045880 1215.3819
SO (26,26.5] 2000-2012 –> 2013-2019 5.8436789 768.2439
SO (26.5,26.75] 1982-1999 –> 2000-2012 6.7728314 3899.2219
SO (26.5,26.75] 2000-2012 –> 2013-2019 6.5505173 2343.2351
SO (26.75,27] 1982-1999 –> 2000-2012 4.9939347 12057.1929
SO (26.75,27] 2000-2012 –> 2013-2019 4.8787943 7460.0025
SO (27,27.25] 1982-1999 –> 2000-2012 3.3296917 22922.6595
SO (27,27.25] 2000-2012 –> 2013-2019 3.5827504 15108.1956
SO (27.25,27.5] 1982-1999 –> 2000-2012 3.1055383 24119.3973
SO (27.25,27.5] 2000-2012 –> 2013-2019 3.4992253 16461.9923
SO (27.5,27.75] 1982-1999 –> 2000-2012 2.7547119 27367.6330
SO (27.5,27.75] 2000-2012 –> 2013-2019 3.2026591 19490.6087
SO (27.75,27.85] 1982-1999 –> 2000-2012 2.6887786 11007.0750
SO (27.75,27.85] 2000-2012 –> 2013-2019 2.9799284 7735.7943
SO (27.85,27.95] 1982-1999 –> 2000-2012 2.6473591 11060.6427
SO (27.85,27.95] 2000-2012 –> 2013-2019 2.8606990 7429.8217
SO (27.95,28.05] 1982-1999 –> 2000-2012 2.8518679 15035.1071
SO (27.95,28.05] 2000-2012 –> 2013-2019 3.2044757 10452.4828
SO (28.05,28.1] 1982-1999 –> 2000-2012 2.0144010 7359.5544
SO (28.05,28.1] 2000-2012 –> 2013-2019 2.2833341 5147.2736
SO (28.1, Inf] 1982-1999 –> 2000-2012 2.0638341 36496.6981
SO (28.1, Inf] 2000-2012 –> 2013-2019 2.2738080 24623.4587
SO (28.1,28.15] 1982-1999 –> 2000-2012 0.8924513 1141.8504
SO (28.1,28.15] 2000-2012 –> 2013-2019 0.9407138 822.7854
SO (28.15,28.2] 1982-1999 –> 2000-2012 0.6064628 425.2148
SO (28.15,28.2] 2000-2012 –> 2013-2019 0.6319172 312.4162
SO (28.2, Inf] 1982-1999 –> 2000-2012 0.6642020 2517.3867
SO (28.2, Inf] 2000-2012 –> 2013-2019 0.7550700 2408.3381

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
db33928 Donghe-Zhu 2021-03-10
3d1f470 Donghe-Zhu 2021-03-10
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
db33928 Donghe-Zhu 2021-03-10
3d1f470 Donghe-Zhu 2021-03-10
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
db33928 Donghe-Zhu 2021-03-10
3d1f470 Donghe-Zhu 2021-03-10
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
db33928 Donghe-Zhu 2021-03-10
3d1f470 Donghe-Zhu 2021-03-10
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 sal silicate temp nitrate
Atlantic - 1982-1999 --> 2000-2012
(-Inf,26] 5 4 5 5 3 3 0
(26,26.5] 4 4 0 5 4 4 5
(26.5,26.75] 4 4 5 4 5 4 0
(26.75,27] 5 4 5 3 5 3 0
(27,27.25] 4 3 5 5 5 3 0
(27.25,27.5] 4 4 0 5 4 4 5
(27.5,27.75] 4 5 2 5 3 5 3
(27.75,27.85] 4 4 5 5 4 4 0
(27.85,27.95] 5 4 5 3 5 3 0
(27.95,28.05] 3 4 5 3 5 5 0
(28.05,28.1] 5 4 5 3 5 3 0
(28.1,28.15] 4 4 5 4 5 4 0
(28.15,28.2] 4 4 0 4 5 4 5
(28.2, Inf] 4 3 0 5 5 3 5
total 59.00 55.00 47.00 59.00 63.00 52.00 23.00
Atlantic - 2000-2012 --> 2013-2019
(-Inf,26] 4 4 4 5 5 4 1
(26,26.5] 4 4 0 5 4 4 5
(26.5,26.75] 4 4 5 4 5 4 0
(26.75,27] 5 3 5 4 5 3 0
(27,27.25] 5 3 4 5 5 3 1
(27.25,27.5] 3 4 0 5 3 5 5
(27.5,27.75] 4 5 3 5 3 5 2
(27.75,27.85] 5 5 3 5 3 4 2
(27.85,27.95] 5 4 5 3 5 3 0
(27.95,28.05] 3 4 5 3 5 5 0
(28.05,28.1] 5 4 5 5 3 3 0
(28.1,28.15] 4 4 5 5 4 4 0
(28.15,28.2] 4 4 0 4 5 4 5
(28.2, Inf] 4 4 0 5 4 4 5
total 59.00 56.00 44.00 63.00 59.00 55.00 26.00
Indo-Pacific - 1982-1999 --> 2000-2012
(-Inf,26] 5 4 5 5 3 3 0
(26,26.5] 4 4 5 4 5 4 0
(26.5,26.75] 5 4 5 3 5 3 0
(26.75,27] 5 5 5 3 3 4 0
(27,27.25] 5 5 5 3 3 4 0
(27.25,27.5] 4 4 5 4 5 4 0
(27.5,27.75] 3 4 5 5 5 3 0
(27.75,27.85] 4 5 5 4 4 4 0
(27.85,27.95] 5 5 2 5 4 4 2
(27.95,28.05] 4 4 5 4 5 4 0
(28.05,28.1] 4 4 5 4 5 4 0
(28.1, Inf] 4 4 5 4 5 4 0
total 52.00 52.00 57.00 48.00 52.00 45.00 2.00
Indo-Pacific - 2000-2012 --> 2013-2019
(-Inf,26] 5 4 5 5 3 3 0
(26,26.5] 4 4 5 4 5 4 0
(26.5,26.75] 5 4 5 3 5 3 0
(26.75,27] 5 5 5 3 3 4 0
(27,27.25] 5 5 4 3 3 5 1
(27.25,27.5] 4 4 5 4 5 4 0
(27.5,27.75] 3 4 5 5 5 3 0
(27.75,27.85] 4 5 5 4 4 4 0
(27.85,27.95] 5 5 2 5 4 4 2
(27.95,28.05] 4 4 5 4 5 4 0
(28.05,28.1] 4 4 5 4 5 4 0
(28.1, Inf] 4 4 5 4 5 4 0
total 52.00 52.00 56.00 48.00 52.00 46.00 3.00
SO - 1982-1999 --> 2000-2012
(26,26.5] 4 4 5 4 4 5 0
(26.5,26.75] 4 5 5 3 3 5 0
(26.75,27] 4 5 5 3 3 5 0
(27,27.25] 5 5 4 3 3 5 1
(27.25,27.5] 4 4 2 4 5 5 3
(27.5,27.75] 4 4 0 4 5 4 5
(27.75,27.85] 4 4 0 4 5 4 5
(27.85,27.95] 5 5 2 5 4 4 2
(27.95,28.05] 5 5 2 5 3 5 2
(28.05,28.1] 5 5 2 5 3 5 2
(28.1, Inf] 3 4 0 5 5 3 5
(28.1,28.15] 5 4 4 5 4 4 1
(28.15,28.2] 4 4 5 5 4 4 0
(28.2, Inf] 4 3 0 5 5 3 5
total 60.00 61.00 36.00 60.00 56.00 61.00 31.00
SO - 2000-2012 --> 2013-2019
(26,26.5] 4 4 5 5 4 4 0
(26.5,26.75] 4 4 5 4 4 5 0
(26.75,27] 4 5 5 3 3 5 0
(27,27.25] 5 5 4 3 3 5 1
(27.25,27.5] 4 4 1 5 5 4 4
(27.5,27.75] 4 3 0 5 5 3 5
(27.75,27.85] 4 3 0 5 5 3 5
(27.85,27.95] 5 5 2 3 5 5 2
(27.95,28.05] 5 5 2 5 3 5 2
(28.05,28.1] 5 5 2 5 3 5 2
(28.1, Inf] 3 4 0 5 5 3 5
(28.1,28.15] 5 5 3 5 4 4 1
(28.15,28.2] 5 5 4 4 4 4 1
(28.2, Inf] 4 3 0 5 5 3 5
total 61.00 60.00 33.00 62.00 58.00 58.00 33.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
db33928 Donghe-Zhu 2021-03-10
3d1f470 Donghe-Zhu 2021-03-10
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
db33928 Donghe-Zhu 2021-03-10
3d1f470 Donghe-Zhu 2021-03-10
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