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

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

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


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

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

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

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

File Version Author Date Message
Rmd dbd8ff8 Donghe-Zhu 2021-03-10 local rebuild
html 94ce5a8 Donghe-Zhu 2021-03-10 Build site.
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_star
# the following code is a workaround to find all predictor combinations
# using the olsrr package and fit all models for one era, slab, and basin

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

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

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

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

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

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

# convert to tibble
lm_all <- as_tibble(lm_all)

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

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

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

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

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

# lm_rm <- lm_rm_ph %>%
#   select(-lm_coeff_filter) %>%
#   unique()
# 
# lm_all <- anti_join(lm_all, lm_rm)

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

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

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

3 Apply predictor threshold

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

  • Minimum: 2
  • Maximum: 9
lm_all <- lm_all %>% 
  filter(n >= params_local$MLR_predictors_min,
         n <= params_local$MLR_predictors_max)

This results in a total number of MLR models of:

  • 64

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 + aou + nitrate + phosphate_star 158 1.6001372 608.9328252 6.6397414 1982-1999 –> 2000-2012 3.2898388 1429.42487
Atlantic (-Inf,26] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 158 1.5921923 609.3599345 6.1184078 1982-1999 –> 2000-2012 3.2779521 1430.88036
Atlantic (-Inf,26] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 158 1.5933338 607.5863993 5.8962936 1982-1999 –> 2000-2012 3.3151005 1435.89880
Atlantic (-Inf,26] 4 cstar_tref ~ sal + oxygen + nitrate + silicate 158 1.6268394 614.1625449 8.3290234 1982-1999 –> 2000-2012 3.2947666 1429.25892
Atlantic (-Inf,26] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 158 1.6216413 615.1512255 8.0736408 1982-1999 –> 2000-2012 3.2856343 1431.26523
Atlantic (-Inf,26] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 72 1.7037742 293.0569596 4.3374446 2000-2012 –> 2013-2019 3.3039114 901.98978
Atlantic (-Inf,26] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 72 1.6825475 293.2516466 3.6000125 2000-2012 –> 2013-2019 3.2747398 902.61158
Atlantic (-Inf,26] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 72 1.6918467 292.0453252 4.1335038 2000-2012 –> 2013-2019 3.2851805 899.63172
Atlantic (-Inf,26] 4 cstar_tref ~ sal + oxygen + nitrate + silicate 72 1.7125308 293.7951588 5.3785542 2000-2012 –> 2013-2019 3.3393703 907.95770
Atlantic (-Inf,26] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 72 1.7102518 295.6033948 5.5140059 2000-2012 –> 2013-2019 3.3318930 910.75462
Atlantic (26,26.5] 4 cstar_tref ~ temp + aou + nitrate + silicate 1040 4.0312783 5863.0858536 14.5182831 1982-1999 –> 2000-2012 8.0820619 14396.26092
Atlantic (26,26.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1040 4.0161192 5857.2495844 16.2265490 1982-1999 –> 2000-2012 8.0668993 14392.42196
Atlantic (26,26.5] 4 cstar_tref ~ temp + oxygen + nitrate + phosphate_star 1040 4.1746063 5935.7538704 10.5920650 1982-1999 –> 2000-2012 8.3053987 14528.07501
Atlantic (26,26.5] 4 cstar_tref ~ temp + oxygen + nitrate + silicate 1040 3.9945290 5844.0375649 14.3000743 1982-1999 –> 2000-2012 8.0094945 14350.35460
Atlantic (26,26.5] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 1040 3.9714887 5834.0054376 16.3663159 1982-1999 –> 2000-2012 7.9849218 14341.16812
Atlantic (26,26.5] 4 cstar_tref ~ temp + aou + nitrate + silicate 562 3.8522413 3122.7752745 13.6121936 2000-2012 –> 2013-2019 7.8835195 8985.86113
Atlantic (26,26.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 562 3.7832546 3104.4640637 14.8635127 2000-2012 –> 2013-2019 7.7993739 8961.71365
Atlantic (26,26.5] 4 cstar_tref ~ temp + oxygen + nitrate + phosphate_star 562 4.0692661 3184.3789542 11.9278024 2000-2012 –> 2013-2019 8.2438724 9120.13282
Atlantic (26,26.5] 4 cstar_tref ~ temp + oxygen + nitrate + silicate 562 3.8236516 3114.4023119 13.3534455 2000-2012 –> 2013-2019 7.8181806 8958.43988
Atlantic (26,26.5] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 562 3.7365212 3090.4931588 14.6731544 2000-2012 –> 2013-2019 7.7080099 8924.49860
Atlantic (26.5,26.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1360 3.4609857 7250.5381279 9.7064368 1982-1999 –> 2000-2012 6.7283877 16827.71654
Atlantic (26.5,26.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1360 3.3435115 7156.6115389 11.0105664 1982-1999 –> 2000-2012 6.4965794 16602.92486
Atlantic (26.5,26.75] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 1360 3.4651325 7251.7951708 10.8919683 1982-1999 –> 2000-2012 6.6787016 16765.93747
Atlantic (26.5,26.75] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 1360 3.3591700 7169.3202672 10.6626989 1982-1999 –> 2000-2012 6.5272991 16633.14132
Atlantic (26.5,26.75] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 1360 3.4729302 7257.9091697 10.9054150 1982-1999 –> 2000-2012 6.6968185 16783.83041
Atlantic (26.5,26.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 646 3.4853031 3460.4016522 9.9449249 2000-2012 –> 2013-2019 6.9462887 10710.93978
Atlantic (26.5,26.75] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 646 3.4893493 3459.9007006 10.2599092 2000-2012 –> 2013-2019 6.9792579 10731.07502
Atlantic (26.5,26.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 646 3.4054258 3430.4466004 9.5018855 2000-2012 –> 2013-2019 6.7489373 10587.05814
Atlantic (26.5,26.75] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 646 3.5365520 3477.2612981 11.4151893 2000-2012 –> 2013-2019 7.0016845 10729.05647
Atlantic (26.5,26.75] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 646 3.4224157 3436.8764420 9.6021654 2000-2012 –> 2013-2019 6.7815857 10606.19671
Atlantic (26.75,27] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 2134 2.3649570 9743.7525764 33.1153986 1982-1999 –> 2000-2012 4.4849721 21594.90562
Atlantic (26.75,27] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2134 2.2560885 9542.6134942 15.8146255 1982-1999 –> 2000-2012 4.0882000 20597.73284
Atlantic (26.75,27] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 2134 2.3393246 9695.2417178 16.9658625 1982-1999 –> 2000-2012 4.2284899 20915.61339
Atlantic (26.75,27] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 2134 2.2995434 9624.0385119 20.3717919 1982-1999 –> 2000-2012 4.1579971 20757.01735
Atlantic (26.75,27] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 2134 2.3692574 9749.5063883 20.0133380 1982-1999 –> 2000-2012 4.2790073 21028.98430
Atlantic (26.75,27] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1121 2.6413897 5372.9263797 21.1206682 2000-2012 –> 2013-2019 5.0063466 15116.67896
Atlantic (26.75,27] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 1121 2.6767952 5400.7788164 23.7846217 2000-2012 –> 2013-2019 5.0479467 15153.69603
Atlantic (26.75,27] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 1121 2.6938814 5417.0442244 21.5824879 2000-2012 –> 2013-2019 5.1593071 15338.36540
Atlantic (26.75,27] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1121 2.7759227 5484.3045017 22.8492521 2000-2012 –> 2013-2019 5.0320112 15026.91800
Atlantic (26.75,27] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 1121 2.8075535 5509.7070069 23.2297681 2000-2012 –> 2013-2019 5.1070970 15133.74552
Atlantic (27,27.25] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1717 2.1245891 7474.4231860 10.8920661 1982-1999 –> 2000-2012 4.3279600 17397.66951
Atlantic (27,27.25] 4 cstar_tref ~ temp + aou + nitrate + silicate 1717 2.1782321 7558.0506248 13.4210426 1982-1999 –> 2000-2012 4.2804948 17268.57061
Atlantic (27,27.25] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1717 2.1502220 7515.6061171 12.8861292 1982-1999 –> 2000-2012 3.8839356 16363.45622
Atlantic (27,27.25] 4 cstar_tref ~ temp + oxygen + nitrate + silicate 1717 2.2153691 7616.1037936 12.9543064 1982-1999 –> 2000-2012 4.3344817 17362.43661
Atlantic (27,27.25] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 1717 2.1930270 7583.2959569 11.9341409 1982-1999 –> 2000-2012 3.9617129 16520.73662
Atlantic (27,27.25] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 982 2.3112528 4446.2142611 10.6624255 2000-2012 –> 2013-2019 4.4358419 11920.63745
Atlantic (27,27.25] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 982 2.4537027 4561.6780993 14.9835525 2000-2012 –> 2013-2019 4.6724715 12183.04779
Atlantic (27,27.25] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 982 2.4101764 4528.5259682 12.9297873 2000-2012 –> 2013-2019 4.6700084 12214.86844
Atlantic (27,27.25] 4 cstar_tref ~ sal + oxygen + silicate + phosphate_star 982 2.4701167 4574.7725503 15.1052250 2000-2012 –> 2013-2019 4.7371689 12270.06903
Atlantic (27,27.25] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 982 2.6296834 4699.7151170 21.1977049 2000-2012 –> 2013-2019 4.7799054 12215.32123
Atlantic (27.25,27.5] 3 cstar_tref ~ temp + aou + nitrate 1621 2.4518761 7517.7976652 15.0667238 1982-1999 –> 2000-2012 4.2959606 16662.88251
Atlantic (27.25,27.5] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 1621 2.4080396 7461.3103609 17.0152980 1982-1999 –> 2000-2012 4.2365015 16570.12710
Atlantic (27.25,27.5] 4 cstar_tref ~ temp + aou + nitrate + silicate 1621 2.4361730 7498.9674242 14.2153496 1982-1999 –> 2000-2012 4.1836374 16403.98301
Atlantic (27.25,27.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1621 2.3121875 7331.6238130 16.5099918 1982-1999 –> 2000-2012 4.0506846 16215.49782
Atlantic (27.25,27.5] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 1621 2.3464706 7379.3404680 13.7613741 1982-1999 –> 2000-2012 4.1649109 16465.43622
Atlantic (27.25,27.5] 3 cstar_tref ~ sal + aou + nitrate 904 2.5315398 4254.7614075 17.3212618 2000-2012 –> 2013-2019 4.8109801 11536.14096
Atlantic (27.25,27.5] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 904 2.5244498 4251.6906907 18.1263712 2000-2012 –> 2013-2019 4.8027876 11533.50176
Atlantic (27.25,27.5] 4 cstar_tref ~ sal + aou + nitrate + silicate 904 2.5311276 4256.4669738 17.4690433 2000-2012 –> 2013-2019 4.7699467 11481.55093
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 904 2.5215176 4251.5894594 18.0824185 2000-2012 –> 2013-2019 4.7464288 11458.47071
Atlantic (27.25,27.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 904 2.8515321 4473.9652800 17.6060112 2000-2012 –> 2013-2019 5.1637196 11805.58909
Atlantic (27.5,27.75] 3 cstar_tref ~ sal + aou + nitrate 2079 2.6053240 9891.4686200 16.4290502 1982-1999 –> 2000-2012 5.2359965 23221.11451
Atlantic (27.5,27.75] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 2079 2.5549227 9812.2417630 15.0971834 1982-1999 –> 2000-2012 5.1855234 23143.73526
Atlantic (27.5,27.75] 4 cstar_tref ~ sal + aou + nitrate + silicate 2079 2.5935623 9874.6549095 16.1196094 1982-1999 –> 2000-2012 5.1154176 22970.49110
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 2079 2.5437146 9795.9611173 14.8081133 1982-1999 –> 2000-2012 5.0655363 22893.72285
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 2079 2.7248162 10081.9296654 14.9156207 1982-1999 –> 2000-2012 5.4634166 23640.01360
Atlantic (27.5,27.75] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 1115 2.8838250 5538.0650607 17.2399006 2000-2012 –> 2013-2019 5.4387477 15350.30682
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1115 2.8834250 5539.7557331 17.3414728 2000-2012 –> 2013-2019 5.4271397 15335.71685
Atlantic (27.5,27.75] 4 cstar_tref ~ sal + nitrate + silicate + phosphate_star 1115 2.9432352 5583.5388798 17.4024408 2000-2012 –> 2013-2019 5.6870507 15692.36032
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 1115 2.8827347 5539.2217910 18.3857778 2000-2012 –> 2013-2019 5.6075509 15621.15146
Atlantic (27.5,27.75] 4 cstar_tref ~ sal + oxygen + silicate + phosphate_star 1115 2.9330221 5575.7872625 17.4417334 2000-2012 –> 2013-2019 5.6734040 15679.40222
Atlantic (27.75,27.85] 3 cstar_tref ~ sal + aou + nitrate 722 1.7460691 2863.7852719 11.8036580 1982-1999 –> 2000-2012 3.6413019 6817.45014
Atlantic (27.75,27.85] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 722 1.6905278 2819.1061822 12.5552253 1982-1999 –> 2000-2012 3.5387296 6726.62484
Atlantic (27.75,27.85] 4 cstar_tref ~ sal + aou + nitrate + silicate 722 1.7417024 2862.1694615 12.0409779 1982-1999 –> 2000-2012 3.6349084 6815.78420
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 722 1.6904824 2821.0673683 12.5742519 1982-1999 –> 2000-2012 3.5382366 6730.12190
Atlantic (27.75,27.85] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 722 1.8502959 2951.5062510 14.7695936 1982-1999 –> 2000-2012 3.6806578 6842.44070
Atlantic (27.75,27.85] 3 cstar_tref ~ sal + aou + nitrate 364 1.8197291 1478.8318604 10.9791463 2000-2012 –> 2013-2019 3.5657983 4342.61713
Atlantic (27.75,27.85] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 364 1.8007426 1473.1962392 11.3084900 2000-2012 –> 2013-2019 3.4912705 4292.30242
Atlantic (27.75,27.85] 4 cstar_tref ~ sal + aou + nitrate + silicate 364 1.8196215 1480.7888133 11.0138981 2000-2012 –> 2013-2019 3.5613239 4342.95827
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 364 1.8006139 1475.1441738 11.2744869 2000-2012 –> 2013-2019 3.4910963 4296.21154
Atlantic (27.75,27.85] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 364 2.1139028 1591.9213770 14.5080776 2000-2012 –> 2013-2019 3.9641987 4543.42763
Atlantic (27.85,27.95] 4 cstar_tref ~ temp + aou + nitrate + silicate 694 2.9129781 3465.5029364 21.0165066 1982-1999 –> 2000-2012 5.9626466 8352.86394
Atlantic (27.85,27.95] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 694 2.9086574 3465.4426035 21.2167232 1982-1999 –> 2000-2012 5.8434073 8280.90157
Atlantic (27.85,27.95] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 694 3.0757988 3540.9944070 21.2352560 1982-1999 –> 2000-2012 6.0292854 8366.69790
Atlantic (27.85,27.95] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 694 3.0238411 3519.3474037 21.2372375 1982-1999 –> 2000-2012 5.9556282 8332.86297
Atlantic (27.85,27.95] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 694 3.0992996 3551.5591643 21.6586223 1982-1999 –> 2000-2012 6.0530208 8377.41553
Atlantic (27.85,27.95] 4 cstar_tref ~ temp + aou + nitrate + silicate 415 3.1804173 2150.0392779 19.0957860 2000-2012 –> 2013-2019 6.0933954 5615.54221
Atlantic (27.85,27.95] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 415 3.1182188 2135.6463927 18.1103786 2000-2012 –> 2013-2019 6.0268762 5601.08900
Atlantic (27.85,27.95] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 415 3.1188968 2133.8268487 18.2013225 2000-2012 –> 2013-2019 6.1946956 5674.82126
Atlantic (27.85,27.95] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 415 3.1150672 2134.8070717 17.7698175 2000-2012 –> 2013-2019 6.1389082 5654.15448
Atlantic (27.85,27.95] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 415 3.1278954 2136.2180882 18.2362943 2000-2012 –> 2013-2019 6.2271949 5687.77725
Atlantic (27.95,28.05] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 723 2.5928883 3441.4940712 13.9368647 1982-1999 –> 2000-2012 4.8700740 7905.86546
Atlantic (27.95,28.05] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 723 2.5849087 3439.0371527 13.5046374 1982-1999 –> 2000-2012 4.8323074 7879.25904
Atlantic (27.95,28.05] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 723 2.8215315 3563.6917628 20.6751203 1982-1999 –> 2000-2012 6.0658508 8731.03808
Atlantic (27.95,28.05] 4 cstar_tref ~ temp + oxygen + nitrate + phosphate_star 723 2.6355853 3465.1113735 14.2708620 1982-1999 –> 2000-2012 4.9203983 7936.12367
Atlantic (27.95,28.05] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 723 2.6146774 3455.5945865 15.0371265 1982-1999 –> 2000-2012 4.8850580 7916.02226
Atlantic (27.95,28.05] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 386 1.9424018 1619.9708427 14.5879054 2000-2012 –> 2013-2019 4.5352901 5061.46491
Atlantic (27.95,28.05] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 386 1.9183328 1612.3449546 15.4301667 2000-2012 –> 2013-2019 4.5032416 5051.38211
Atlantic (27.95,28.05] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 386 2.4606514 1802.5495198 13.8311274 2000-2012 –> 2013-2019 5.2821830 5366.24128
Atlantic (27.95,28.05] 4 cstar_tref ~ temp + oxygen + nitrate + phosphate_star 386 1.9514658 1623.5648945 15.4815166 2000-2012 –> 2013-2019 4.5870511 5088.67627
Atlantic (27.95,28.05] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 386 1.9374096 1619.9841602 16.1697987 2000-2012 –> 2013-2019 4.5520870 5075.57875
Atlantic (28.05,28.1] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 436 0.2026977 -140.4320462 1.1989259 1982-1999 –> 2000-2012 0.3184343 -1007.00668
Atlantic (28.05,28.1] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 436 0.2223078 -61.9053227 1.4225064 1982-1999 –> 2000-2012 0.3430222 -880.20070
Atlantic (28.05,28.1] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 436 0.2117907 -104.1663554 1.3461511 1982-1999 –> 2000-2012 0.3300322 -947.17451
Atlantic (28.05,28.1] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 436 0.2012900 -146.5090680 1.1821042 1982-1999 –> 2000-2012 0.3161426 -1022.23922
Atlantic (28.05,28.1] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 436 0.2230988 -58.8082625 1.4266656 1982-1999 –> 2000-2012 0.3442268 -873.01874
Atlantic (28.05,28.1] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 229 0.2594994 44.0314618 0.8986441 2000-2012 –> 2013-2019 0.4771556 -36.31332
Atlantic (28.05,28.1] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 229 0.2412331 12.6016075 1.0576928 2000-2012 –> 2013-2019 0.4439308 -127.83044
Atlantic (28.05,28.1] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 229 0.2498995 26.7669599 1.1561502 2000-2012 –> 2013-2019 0.4616902 -77.39940
Atlantic (28.05,28.1] 4 cstar_tref ~ temp + oxygen + nitrate + phosphate_star 229 0.2596292 44.2603895 0.9023449 2000-2012 –> 2013-2019 0.4773394 -35.86801
Atlantic (28.05,28.1] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 229 0.2403679 10.9560254 1.0657336 2000-2012 –> 2013-2019 0.4416579 -135.55304
Atlantic (28.1,28.15] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 541 0.1576027 -451.8958141 1.4966383 1982-1999 –> 2000-2012 0.2858731 -1422.37448
Atlantic (28.1,28.15] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 541 0.1549911 -467.9756571 1.5279244 1982-1999 –> 2000-2012 0.2748469 -1541.48871
Atlantic (28.1,28.15] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 541 0.1759643 -332.6556168 1.4526178 1982-1999 –> 2000-2012 0.3270778 -1049.43197
Atlantic (28.1,28.15] 4 cstar_tref ~ temp + oxygen + nitrate + phosphate_star 541 0.1561830 -461.6869484 1.4893757 1982-1999 –> 2000-2012 0.2843957 -1432.86267
Atlantic (28.1,28.15] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 541 0.1540499 -474.5665659 1.5195644 1982-1999 –> 2000-2012 0.2741208 -1545.30295
Atlantic (28.1,28.15] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 288 0.1706821 -189.0321719 0.9575421 2000-2012 –> 2013-2019 0.3282848 -640.92799
Atlantic (28.1,28.15] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 288 0.1624875 -215.3720860 1.0609876 2000-2012 –> 2013-2019 0.3174787 -683.34774
Atlantic (28.1,28.15] 4 cstar_tref ~ temp + oxygen + nitrate + phosphate_star 288 0.1684916 -196.4720512 0.9439858 2000-2012 –> 2013-2019 0.3246746 -658.15900
Atlantic (28.1,28.15] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 288 0.1613135 -219.5489446 1.0473782 2000-2012 –> 2013-2019 0.3153634 -694.11551
Atlantic (28.1,28.15] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 288 0.1720813 -184.3295683 0.9678522 2000-2012 –> 2013-2019 0.3456807 -531.62529
Atlantic (28.15,28.2] 3 cstar_tref ~ sal + aou + nitrate 1155 0.1574471 -982.6698392 0.9830553 1982-1999 –> 2000-2012 0.3636269 -1527.46540
Atlantic (28.15,28.2] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 1155 0.1550542 -1016.0477142 0.9483180 1982-1999 –> 2000-2012 0.3595782 -1586.79067
Atlantic (28.15,28.2] 4 cstar_tref ~ sal + aou + nitrate + silicate 1155 0.1567453 -990.9892289 0.9661276 1982-1999 –> 2000-2012 0.3491822 -1772.87154
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1155 0.1338912 -1353.0327195 0.6859787 1982-1999 –> 2000-2012 0.3082551 -2474.74523
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 1155 0.1609252 -928.1953795 0.7964617 1982-1999 –> 2000-2012 0.3834485 -1204.59395
Atlantic (28.15,28.2] 3 cstar_tref ~ sal + aou + nitrate 636 0.2156021 -136.7663848 1.9725799 2000-2012 –> 2013-2019 0.3730492 -1119.43622
Atlantic (28.15,28.2] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 636 0.2000269 -230.1440140 1.4705162 2000-2012 –> 2013-2019 0.3550811 -1246.19173
Atlantic (28.15,28.2] 4 cstar_tref ~ sal + aou + nitrate + silicate 636 0.1903257 -293.3813460 1.3866890 2000-2012 –> 2013-2019 0.3470711 -1284.37057
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 636 0.1879213 -307.5531446 1.4575343 2000-2012 –> 2013-2019 0.3218125 -1660.58586
Atlantic (28.15,28.2] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 636 0.2394319 0.5830451 1.9671533 2000-2012 –> 2013-2019 0.4098310 -795.47123
Atlantic (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 955 0.1541931 -848.6675542 0.7135736 1982-1999 –> 2000-2012 0.3663878 -1162.05962
Atlantic (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + silicate 955 0.1744896 -612.4775343 0.9944731 1982-1999 –> 2000-2012 0.4081182 -687.41598
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 955 0.1427395 -994.0896712 0.8634543 1982-1999 –> 2000-2012 0.3506627 -1355.87253
Atlantic (28.2, Inf] 4 cstar_tref ~ sal + oxygen + nitrate + phosphate_star 955 0.2140864 -221.8551418 2.5899588 1982-1999 –> 2000-2012 0.4568162 -202.09333
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 955 0.2134016 -225.9742680 2.5091980 1982-1999 –> 2000-2012 0.4561179 -204.35048
Atlantic (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 552 0.2005924 -195.0462133 1.3544749 2000-2012 –> 2013-2019 0.3547854 -1043.71377
Atlantic (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + silicate 552 0.2343110 -23.5130276 1.6064160 2000-2012 –> 2013-2019 0.4088006 -635.99056
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 552 0.1875683 -267.1597642 0.8704042 2000-2012 –> 2013-2019 0.3303078 -1261.24944
Atlantic (28.2, Inf] 4 cstar_tref ~ sal + oxygen + nitrate + phosphate_star 552 0.1756256 -341.7901968 1.3611522 2000-2012 –> 2013-2019 0.3897120 -563.64534
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 552 0.1636030 -418.0766403 0.9738777 2000-2012 –> 2013-2019 0.3770046 -644.05091
Indo-Pacific (-Inf,26] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5288 7.1987438 35896.7294426 31.3718077 1982-1999 –> 2000-2012 14.7274968 85351.27007
Indo-Pacific (-Inf,26] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 5288 7.4060853 36195.0399372 28.5153217 1982-1999 –> 2000-2012 15.1440696 86041.81693
Indo-Pacific (-Inf,26] 4 cstar_tref ~ temp + oxygen + nitrate + phosphate_star 5288 7.5583116 36410.2174511 34.8822848 1982-1999 –> 2000-2012 15.4934811 86618.89571
Indo-Pacific (-Inf,26] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 5288 7.1609299 35841.0290108 31.5543548 1982-1999 –> 2000-2012 14.6429836 85206.08317
Indo-Pacific (-Inf,26] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 5288 7.3907730 36173.1510635 28.4925990 1982-1999 –> 2000-2012 15.1057280 85977.06154
Indo-Pacific (-Inf,26] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2815 6.9518302 18919.2219176 25.9076786 2000-2012 –> 2013-2019 14.1505740 54815.95136
Indo-Pacific (-Inf,26] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 2815 7.1476758 19073.6361248 24.4727304 2000-2012 –> 2013-2019 14.5537612 55268.67606
Indo-Pacific (-Inf,26] 4 cstar_tref ~ temp + oxygen + nitrate + phosphate_star 2815 7.4095774 19276.2382450 28.8913693 2000-2012 –> 2013-2019 14.9678891 55686.45570
Indo-Pacific (-Inf,26] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 2815 6.9175391 18891.3821880 26.0064859 2000-2012 –> 2013-2019 14.0784690 54732.41120
Indo-Pacific (-Inf,26] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 2815 7.1346964 19063.4033807 24.4397815 2000-2012 –> 2013-2019 14.5254694 55236.55444
Indo-Pacific (26,26.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5428 4.7468521 32326.0177685 49.2620011 1982-1999 –> 2000-2012 9.1449428 76398.47442
Indo-Pacific (26,26.5] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 5428 4.9237293 32721.1797544 49.7471800 1982-1999 –> 2000-2012 9.5483116 77554.51099
Indo-Pacific (26,26.5] 4 cstar_tref ~ temp + oxygen + nitrate + phosphate_star 5428 5.2629233 33444.4108717 54.3987847 1982-1999 –> 2000-2012 10.1328347 79063.01160
Indo-Pacific (26,26.5] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 5428 4.7562153 32347.4102632 51.4538473 1982-1999 –> 2000-2012 9.1587692 76435.27599
Indo-Pacific (26,26.5] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 5428 4.9303396 32735.7446235 51.0420312 1982-1999 –> 2000-2012 9.5603521 77586.90379
Indo-Pacific (26,26.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2876 4.5772124 16925.0450303 45.7195416 2000-2012 –> 2013-2019 9.3240644 49251.06280
Indo-Pacific (26,26.5] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 2876 4.8581955 17265.7314706 34.3531399 2000-2012 –> 2013-2019 9.7819248 49986.91123
Indo-Pacific (26,26.5] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 2876 5.0114519 17444.3805309 32.0664466 2000-2012 –> 2013-2019 10.2914304 50923.91483
Indo-Pacific (26,26.5] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 2876 4.5840145 16933.5866150 45.2479329 2000-2012 –> 2013-2019 9.3402297 49280.99688
Indo-Pacific (26,26.5] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 2876 4.8618053 17270.0037363 34.1438489 2000-2012 –> 2013-2019 9.7921449 50005.74836
Indo-Pacific (26.5,26.75] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 4269 4.1159083 24206.9679477 35.7079690 1982-1999 –> 2000-2012 8.1757605 57653.89704
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 4269 4.0277915 24024.1942408 33.2108036 1982-1999 –> 2000-2012 8.0184863 57269.42280
Indo-Pacific (26.5,26.75] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 4269 4.0285176 24023.7333576 32.9352042 1982-1999 –> 2000-2012 8.0934769 57485.56713
Indo-Pacific (26.5,26.75] 4 cstar_tref ~ temp + oxygen + nitrate + phosphate_star 4269 4.1144340 24203.9091240 35.4385661 1982-1999 –> 2000-2012 8.1771095 57659.08006
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 4269 4.0282209 24025.1045514 33.1411767 1982-1999 –> 2000-2012 8.0241088 57285.75097
Indo-Pacific (26.5,26.75] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 2324 4.4297229 13525.0168465 33.1064604 2000-2012 –> 2013-2019 8.5456312 37731.98479
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2324 4.3368286 13428.5085435 32.2595908 2000-2012 –> 2013-2019 8.3646201 37452.70278
Indo-Pacific (26.5,26.75] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 2324 4.3478487 13438.3043881 34.5430589 2000-2012 –> 2013-2019 8.3763663 37462.03775
Indo-Pacific (26.5,26.75] 4 cstar_tref ~ temp + oxygen + nitrate + phosphate_star 2324 4.4268240 13521.9740380 32.9280194 2000-2012 –> 2013-2019 8.5412580 37725.88316
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 2324 4.3351224 13426.6795295 32.0750065 2000-2012 –> 2013-2019 8.3633433 37451.78408
Indo-Pacific (26.75,27] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 4557 4.8429822 27321.8204666 21.3675951 1982-1999 –> 2000-2012 9.1459669 64342.32595
Indo-Pacific (26.75,27] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 4557 4.8279313 27295.4520360 21.1638765 1982-1999 –> 2000-2012 9.1211296 64288.68115
Indo-Pacific (26.75,27] 4 cstar_tref ~ temp + oxygen + nitrate + phosphate_star 4557 4.7938814 27228.9461788 20.4474901 1982-1999 –> 2000-2012 9.0679743 64162.82732
Indo-Pacific (26.75,27] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 4557 4.7734283 27191.9781616 20.1346476 1982-1999 –> 2000-2012 9.0333117 64085.04076
Indo-Pacific (26.75,27] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 4557 4.8874824 27405.1828608 22.0215384 1982-1999 –> 2000-2012 9.2350911 64558.34546
Indo-Pacific (26.75,27] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 2460 5.2579457 15159.1003849 19.1176723 2000-2012 –> 2013-2019 10.1009279 42480.92085
Indo-Pacific (26.75,27] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2460 5.2437406 15147.7902440 18.5642516 2000-2012 –> 2013-2019 10.0716718 42443.24228
Indo-Pacific (26.75,27] 4 cstar_tref ~ temp + oxygen + nitrate + phosphate_star 2460 5.1910578 15096.1100164 18.7070307 2000-2012 –> 2013-2019 9.9849392 42325.05620
Indo-Pacific (26.75,27] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 2460 5.1714478 15079.4888523 18.0525304 2000-2012 –> 2013-2019 9.9448762 42271.46701
Indo-Pacific (26.75,27] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 2460 5.3165788 15213.6612807 15.8748678 2000-2012 –> 2013-2019 10.2040612 42618.84414
Indo-Pacific (27,27.25] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 4267 4.5843810 25115.5599684 57.7149040 1982-1999 –> 2000-2012 8.1980626 56311.63924
Indo-Pacific (27,27.25] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 4267 4.3243017 24619.1367362 50.1571966 1982-1999 –> 2000-2012 7.7734614 55279.77227
Indo-Pacific (27,27.25] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 4267 4.5313368 25016.2405221 62.9179070 1982-1999 –> 2000-2012 7.9914714 55711.51811
Indo-Pacific (27,27.25] 4 cstar_tref ~ temp + oxygen + nitrate + phosphate_star 4267 4.5777806 25103.2642633 57.1639566 1982-1999 –> 2000-2012 8.1990179 56323.43415
Indo-Pacific (27,27.25] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 4267 4.2969843 24565.0548215 49.1433506 1982-1999 –> 2000-2012 7.7403218 55206.20483
Indo-Pacific (27,27.25] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 2344 5.3519510 14529.9218387 20.7727658 2000-2012 –> 2013-2019 10.1903141 40107.64678
Indo-Pacific (27,27.25] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2344 5.2740156 14461.1529911 34.1511699 2000-2012 –> 2013-2019 9.5983174 39080.28973
Indo-Pacific (27,27.25] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 2344 5.4504530 14613.4194305 30.1690654 2000-2012 –> 2013-2019 10.1259553 39896.94151
Indo-Pacific (27,27.25] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 2344 5.2055927 14399.9348141 33.2928122 2000-2012 –> 2013-2019 9.5025770 38964.98964
Indo-Pacific (27,27.25] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 2344 5.3701277 14543.8165582 29.3952372 2000-2012 –> 2013-2019 10.0017331 39746.83681
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 3520 3.0170713 17777.5047944 30.6141183 1982-1999 –> 2000-2012 5.8968935 41156.38222
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 3520 2.9700972 17667.0338834 30.3495438 1982-1999 –> 2000-2012 5.8222458 40954.81698
Indo-Pacific (27.25,27.5] 4 cstar_tref ~ sal + oxygen + silicate + phosphate_star 3520 3.1493898 18077.6766494 29.8455275 1982-1999 –> 2000-2012 6.0757661 41605.84076
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 3520 3.0094328 17759.6586184 29.1482995 1982-1999 –> 2000-2012 5.9385710 41298.72235
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 3520 3.0036251 17746.0594371 29.2061733 1982-1999 –> 2000-2012 5.9288624 41272.55108
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1879 3.4636269 10014.9955686 28.8274843 2000-2012 –> 2013-2019 6.4806982 27792.50036
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 1879 3.4196778 9967.0060932 28.3215760 2000-2012 –> 2013-2019 6.3897750 27634.03998
Indo-Pacific (27.25,27.5] 4 cstar_tref ~ sal + oxygen + silicate + phosphate_star 1879 3.6875866 10248.4565706 27.7368313 2000-2012 –> 2013-2019 6.8369764 28326.13322
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1879 3.5497395 10107.2843807 33.8094442 2000-2012 –> 2013-2019 6.5591723 27866.94300
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 1879 3.5422402 10099.3366498 33.5571287 2000-2012 –> 2013-2019 6.5458652 27845.39609
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 3170 1.8542474 12924.8866415 13.2838580 1982-1999 –> 2000-2012 3.4799789 30094.38329
Indo-Pacific (27.5,27.75] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 3170 1.9027986 13086.7555124 12.4473296 1982-1999 –> 2000-2012 3.5535107 30391.58410
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 3170 1.8700935 12978.8369461 14.8905892 1982-1999 –> 2000-2012 3.4940074 30138.25969
Indo-Pacific (27.5,27.75] 4 cstar_tref ~ sal + oxygen + silicate + phosphate_star 3170 1.9247016 13159.3181099 14.1925496 1982-1999 –> 2000-2012 3.5760715 30467.73487
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 3170 1.9756992 13327.1182198 12.8757708 1982-1999 –> 2000-2012 3.6722342 30880.54150
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1683 2.1594104 7381.4124538 13.7085041 2000-2012 –> 2013-2019 4.0136578 20306.29910
Indo-Pacific (27.5,27.75] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 1683 2.2370458 7498.3027878 13.7718692 2000-2012 –> 2013-2019 4.1398444 20585.05830
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 1683 2.2078715 7456.1164551 15.6847995 2000-2012 –> 2013-2019 4.0779650 20434.95340
Indo-Pacific (27.5,27.75] 4 cstar_tref ~ sal + oxygen + silicate + phosphate_star 1683 2.2955622 7585.2182760 14.9967179 2000-2012 –> 2013-2019 4.2202638 20744.53639
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 1683 2.2884244 7576.7357604 13.2590663 2000-2012 –> 2013-2019 4.2641236 20903.85398
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1216 1.4724348 4405.8414666 9.7647262 1982-1999 –> 2000-2012 2.8012636 10169.96785
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1216 1.4429788 4356.6962523 8.1238690 1982-1999 –> 2000-2012 2.7877140 10160.99422
Indo-Pacific (27.75,27.85] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 1216 1.4433627 4355.3431320 8.1927273 1982-1999 –> 2000-2012 2.7893920 10160.88841
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 1216 1.4440743 4358.5418520 8.1544572 1982-1999 –> 2000-2012 2.7889516 10163.19670
Indo-Pacific (27.75,27.85] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 1216 1.4444374 4357.1532702 8.2212238 1982-1999 –> 2000-2012 2.7905233 10162.84055
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 691 1.6308082 2650.8757184 15.5062349 2000-2012 –> 2013-2019 3.1032430 7056.71719
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 691 1.6406123 2659.1591558 13.9456772 2000-2012 –> 2013-2019 3.0835911 7015.85541
Indo-Pacific (27.75,27.85] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 691 1.6406255 2657.1702968 13.9682355 2000-2012 –> 2013-2019 3.0839882 7012.51343
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 691 1.6412704 2659.7133533 13.9902086 2000-2012 –> 2013-2019 3.0853446 7018.25521
Indo-Pacific (27.75,27.85] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 691 1.6412785 2657.7201870 14.0078568 2000-2012 –> 2013-2019 3.0857158 7014.87346
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1797 1.4125952 6355.1355196 11.3038103 1982-1999 –> 2000-2012 2.6595697 14669.08780
Indo-Pacific (27.85,27.95] 4 cstar_tref ~ temp + aou + nitrate + silicate 1797 1.3683349 6238.7241576 9.2082326 1982-1999 –> 2000-2012 2.6475550 14679.91157
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1797 1.3568501 6210.4314395 9.3793503 1982-1999 –> 2000-2012 2.5976136 14499.10798
Indo-Pacific (27.85,27.95] 4 cstar_tref ~ temp + oxygen + nitrate + silicate 1797 1.3690177 6240.5169444 9.2508865 1982-1999 –> 2000-2012 2.6487710 14683.81411
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 1797 1.3580391 6213.5794139 9.4298260 1982-1999 –> 2000-2012 2.6001429 14507.72119
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 987 1.6494324 3802.8359632 12.1825689 2000-2012 –> 2013-2019 3.0620277 10157.97148
Indo-Pacific (27.85,27.95] 4 cstar_tref ~ temp + aou + nitrate + silicate 987 1.6038431 3745.5076106 8.9600325 2000-2012 –> 2013-2019 2.9721781 9984.23177
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 987 1.5960835 3737.9339228 9.2774322 2000-2012 –> 2013-2019 2.9529336 9948.36536
Indo-Pacific (27.85,27.95] 4 cstar_tref ~ temp + oxygen + nitrate + silicate 987 1.6044265 3746.2255310 9.0053111 2000-2012 –> 2013-2019 2.9734442 9986.74248
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 987 1.5970590 3739.1400061 9.3241425 2000-2012 –> 2013-2019 2.9550981 9952.71942
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1162 0.6596354 2344.6710734 5.9778740 1982-1999 –> 2000-2012 1.2522936 5474.26713
Indo-Pacific (27.95,28.05] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 1162 0.6867923 2436.4322045 3.4262432 1982-1999 –> 2000-2012 1.3076557 5725.73146
Indo-Pacific (27.95,28.05] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 1162 0.6598787 2343.5280283 5.7094931 1982-1999 –> 2000-2012 1.2528152 5472.75732
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 1162 0.6593779 2343.7636371 6.0893651 1982-1999 –> 2000-2012 1.2517940 5471.93904
Indo-Pacific (27.95,28.05] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 1162 0.6891088 2444.2574253 3.4403173 1982-1999 –> 2000-2012 1.3122791 5746.45639
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 648 0.6715663 1336.9516306 3.3640626 2000-2012 –> 2013-2019 1.3312017 3681.62270
Indo-Pacific (27.95,28.05] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 648 0.6863576 1363.1863850 3.2348373 2000-2012 –> 2013-2019 1.3731500 3799.61859
Indo-Pacific (27.95,28.05] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 648 0.6715817 1334.9813500 3.4211288 2000-2012 –> 2013-2019 1.3314604 3678.50938
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 648 0.6715777 1336.9737100 3.4499367 2000-2012 –> 2013-2019 1.3309556 3680.73735
Indo-Pacific (27.95,28.05] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 648 0.6880082 1366.2992959 3.2459241 2000-2012 –> 2013-2019 1.3771170 3810.55672
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1015 0.4795310 1402.5034846 3.1382674 1982-1999 –> 2000-2012 0.9490005 3430.04886
Indo-Pacific (28.05,28.1] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 1015 0.4952435 1465.9524385 2.6603949 1982-1999 –> 2000-2012 0.9763616 3565.95826
Indo-Pacific (28.05,28.1] 4 cstar_tref ~ sal + nitrate + silicate + phosphate_star 1015 0.4832957 1416.3783333 2.5369659 1982-1999 –> 2000-2012 0.9535442 3446.96091
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 1015 0.4720649 1370.6486070 3.6896280 1982-1999 –> 2000-2012 0.9379775 3375.08919
Indo-Pacific (28.05,28.1] 4 cstar_tref ~ sal + oxygen + silicate + phosphate_star 1015 0.5024246 1495.1766922 2.7042536 1982-1999 –> 2000-2012 0.9908390 3640.90839
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 565 0.4790009 785.6609647 2.5225697 2000-2012 –> 2013-2019 0.9585320 2188.16445
Indo-Pacific (28.05,28.1] 4 cstar_tref ~ sal + nitrate + silicate + phosphate_star 565 0.4904495 810.3512236 2.7325905 2000-2012 –> 2013-2019 0.9737452 2226.72956
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 565 0.4665997 756.0202249 2.6691654 2000-2012 –> 2013-2019 0.9386647 2126.66883
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 565 0.4925610 817.2057193 2.9003756 2000-2012 –> 2013-2019 0.9977401 2325.48118
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 565 0.4914534 814.6618814 2.9243530 2000-2012 –> 2013-2019 0.9962342 2321.33609
Indo-Pacific (28.1, Inf] 4 cstar_tref ~ oxygen + nitrate + silicate + phosphate_star 9318 0.7121604 20129.3080619 3.2377817 1982-1999 –> 2000-2012 1.4288133 48319.33964
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 9318 0.7861340 21972.9954474 3.7460043 1982-1999 –> 2000-2012 1.5770507 52723.91650
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 9318 0.7048635 19939.3772377 3.1418437 1982-1999 –> 2000-2012 1.4079498 47635.41142
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 9318 0.6753077 19141.0900046 2.6646187 1982-1999 –> 2000-2012 1.3451269 45579.18684
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 9318 0.6732113 19083.1461670 2.6753591 1982-1999 –> 2000-2012 1.3406036 45427.04217
Indo-Pacific (28.1, Inf] 4 cstar_tref ~ oxygen + nitrate + silicate + phosphate_star 5093 0.7126114 11014.0973107 3.4303306 2000-2012 –> 2013-2019 1.4247718 31143.40537
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 5093 0.7791371 11925.2057059 4.2437545 2000-2012 –> 2013-2019 1.5652711 33898.20115
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 5093 0.7039983 10892.2318125 3.3161171 2000-2012 –> 2013-2019 1.4088618 30831.60905
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5093 0.6800660 10539.9378299 3.0769840 2000-2012 –> 2013-2019 1.3553737 29681.02783
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 5093 0.6779216 10507.7692399 3.0793019 2000-2012 –> 2013-2019 1.3511329 29590.91541
SO (26,26.5] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 103 1.9289815 441.6417224 8.5446739 1982-1999 –> 2000-2012 4.5881505 1074.05539
SO (26,26.5] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 103 1.0023937 304.7938596 4.5747847 1982-1999 –> 2000-2012 2.3237890 754.78155
SO (26,26.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 103 0.9768158 301.4691489 3.5663094 1982-1999 –> 2000-2012 2.1945946 732.38866
SO (26,26.5] 4 cstar_tref ~ temp + oxygen + nitrate + phosphate_star 103 1.1061695 325.0873780 5.3380492 1982-1999 –> 2000-2012 2.5894313 804.88838
SO (26,26.5] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 103 1.0511102 316.5698067 3.8353803 1982-1999 –> 2000-2012 2.3888479 771.72875
SO (26,26.5] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 47 2.1638927 219.9396465 6.6773028 2000-2012 –> 2013-2019 4.0928742 661.58137
SO (26,26.5] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 47 1.1019644 154.5070979 2.3246978 2000-2012 –> 2013-2019 2.1043582 459.30096
SO (26,26.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 47 1.0381247 150.8973184 2.2629736 2000-2012 –> 2013-2019 2.0149405 452.36647
SO (26,26.5] 4 cstar_tref ~ temp + oxygen + nitrate + phosphate_star 47 1.1971200 162.2925751 2.3579846 2000-2012 –> 2013-2019 2.3032894 487.37995
SO (26,26.5] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 47 1.1040262 156.6828085 2.2663783 2000-2012 –> 2013-2019 2.1551364 473.25262
SO (26.5,26.75] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 312 2.3977151 1445.1157772 13.1971265 1982-1999 –> 2000-2012 5.1969906 3515.68707
SO (26.5,26.75] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 312 1.6818786 1221.8423319 7.4624349 1982-1999 –> 2000-2012 3.5967902 2971.47474
SO (26.5,26.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 312 1.6697360 1219.3209394 7.1133689 1982-1999 –> 2000-2012 3.5488640 2955.10787
SO (26.5,26.75] 4 cstar_tref ~ temp + oxygen + nitrate + phosphate_star 312 1.8198506 1271.0403995 7.7301109 1982-1999 –> 2000-2012 3.8762708 3080.56079
SO (26.5,26.75] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 312 1.8045776 1267.7814048 7.3425552 1982-1999 –> 2000-2012 3.8145175 3060.09785
SO (26.5,26.75] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 133 2.2050757 601.7803004 7.1755301 2000-2012 –> 2013-2019 4.6027908 2046.89608
SO (26.5,26.75] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 133 1.6548865 523.4304747 4.9197792 2000-2012 –> 2013-2019 3.3367651 1745.27281
SO (26.5,26.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 133 1.6547816 525.4136173 4.9257339 2000-2012 –> 2013-2019 3.3245176 1744.73456
SO (26.5,26.75] 4 cstar_tref ~ temp + oxygen + nitrate + phosphate_star 133 1.7625663 540.1986949 4.9732020 2000-2012 –> 2013-2019 3.5824169 1811.23909
SO (26.5,26.75] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 133 1.7623504 542.1661107 4.9816205 2000-2012 –> 2013-2019 3.5669279 1809.94752
SO (26.75,27] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 1061 2.4677656 4939.8180549 13.5857513 1982-1999 –> 2000-2012 5.1625431 12254.88820
SO (26.75,27] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1061 2.4672141 4941.3437400 13.4599518 1982-1999 –> 2000-2012 5.1618710 12258.27828
SO (26.75,27] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 1061 2.9578329 5324.2050702 10.9471230 1982-1999 –> 2000-2012 6.0140839 13020.67143
SO (26.75,27] 4 cstar_tref ~ temp + oxygen + nitrate + phosphate_star 1061 2.5622442 5019.5424051 12.9363272 1982-1999 –> 2000-2012 5.3328558 12418.70243
SO (26.75,27] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 1061 2.5621875 5021.4954099 12.8935474 1982-1999 –> 2000-2012 5.3326619 12422.50543
SO (26.75,27] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 550 2.5663729 2609.5753272 17.6995339 2000-2012 –> 2013-2019 5.0341385 7549.39338
SO (26.75,27] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 550 2.5605544 2609.0785756 17.0705169 2000-2012 –> 2013-2019 5.0277685 7550.42232
SO (26.75,27] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 550 3.0086847 2784.4856997 12.4224158 2000-2012 –> 2013-2019 5.9665176 8108.69077
SO (26.75,27] 4 cstar_tref ~ temp + oxygen + nitrate + phosphate_star 550 2.6716534 2653.7996514 16.9241953 2000-2012 –> 2013-2019 5.2338976 7673.34206
SO (26.75,27] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 550 2.6687763 2654.6144635 16.4534147 2000-2012 –> 2013-2019 5.2309638 7676.10987
SO (27,27.25] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 2513 1.9002603 10370.2312203 19.2728362 1982-1999 –> 2000-2012 3.7237639 24269.67624
SO (27,27.25] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2513 1.8530732 10245.8502101 18.9619811 1982-1999 –> 2000-2012 3.6704902 24124.30576
SO (27,27.25] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 2513 1.9463574 10490.6980852 19.9217964 1982-1999 –> 2000-2012 3.8581729 24715.33421
SO (27,27.25] 4 cstar_tref ~ temp + oxygen + nitrate + phosphate_star 2513 1.9401112 10474.5429330 19.6469841 1982-1999 –> 2000-2012 3.8009118 24513.20726
SO (27,27.25] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 2513 1.8909124 10347.4459479 19.3045392 1982-1999 –> 2000-2012 3.7443899 24360.99662
SO (27,27.25] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 1286 2.1008436 5570.8058048 13.9602984 2000-2012 –> 2013-2019 4.0011039 15941.03703
SO (27,27.25] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1286 2.0262434 5479.8139370 13.8567921 2000-2012 –> 2013-2019 3.8793165 15725.66415
SO (27,27.25] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 1286 2.1184198 5592.2342269 14.6526181 2000-2012 –> 2013-2019 4.0647771 16082.93231
SO (27,27.25] 4 cstar_tref ~ temp + oxygen + nitrate + phosphate_star 1286 2.1435716 5622.5914916 14.3344432 2000-2012 –> 2013-2019 4.0836828 16097.13442
SO (27,27.25] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 1286 2.0660156 5529.8094651 14.1995949 2000-2012 –> 2013-2019 3.9569280 15877.25541
SO (27.25,27.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2731 1.6825482 10606.1722440 8.7665619 1982-1999 –> 2000-2012 3.1211793 24182.42563
SO (27.25,27.5] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 2731 1.6893649 10626.2566473 10.5431309 1982-1999 –> 2000-2012 3.1306658 24214.61540
SO (27.25,27.5] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 2731 1.6834191 10606.9989239 8.1155062 1982-1999 –> 2000-2012 3.1248601 24196.09752
SO (27.25,27.5] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 2731 1.6827611 10606.8633951 8.6884540 1982-1999 –> 2000-2012 3.1222924 24187.87591
SO (27.25,27.5] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 2731 1.6897063 10627.3602422 10.5063777 1982-1999 –> 2000-2012 3.1326599 24224.43821
SO (27.25,27.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1453 1.8604040 5941.4617956 9.5764858 2000-2012 –> 2013-2019 3.5429522 16547.63404
SO (27.25,27.5] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 1453 1.8608548 5940.1658157 9.5050972 2000-2012 –> 2013-2019 3.5502197 16566.42246
SO (27.25,27.5] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 1453 1.8710766 5956.0850177 9.8630288 2000-2012 –> 2013-2019 3.5544957 16563.08394
SO (27.25,27.5] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 1453 1.8591019 5939.4271407 9.5678324 2000-2012 –> 2013-2019 3.5418630 16546.29054
SO (27.25,27.5] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 1453 1.8592539 5937.6646737 9.5274690 2000-2012 –> 2013-2019 3.5489602 16565.02492
SO (27.5,27.75] 4 cstar_tref ~ aou + nitrate + silicate + phosphate_star 3391 1.4983094 12377.4572107 8.4631334 1982-1999 –> 2000-2012 2.7575978 27383.27911
SO (27.5,27.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 3391 1.4809706 12300.5166921 8.2778025 1982-1999 –> 2000-2012 2.7401872 27307.82001
SO (27.5,27.75] 4 cstar_tref ~ temp + aou + nitrate + silicate 3391 1.5314112 12525.6595815 8.9583230 1982-1999 –> 2000-2012 2.7897750 27524.80452
SO (27.5,27.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 3391 1.4865171 12325.8693341 8.1453800 1982-1999 –> 2000-2012 2.7423144 27308.45606
SO (27.5,27.75] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 3391 1.4728687 12263.3130481 7.7650825 1982-1999 –> 2000-2012 2.7255384 27223.23217
SO (27.5,27.75] 4 cstar_tref ~ aou + nitrate + silicate + phosphate_star 1821 1.7424835 7202.2181577 10.3871064 2000-2012 –> 2013-2019 3.2407928 19579.67537
SO (27.5,27.75] 4 cstar_tref ~ sal + aou + nitrate + silicate 1821 1.7135136 7141.1585933 9.9489584 2000-2012 –> 2013-2019 3.2205821 19558.14808
SO (27.5,27.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1821 1.6982942 7110.6659963 9.9919662 2000-2012 –> 2013-2019 3.1792648 19411.18269
SO (27.5,27.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1821 1.6931774 7099.6764854 9.9569404 2000-2012 –> 2013-2019 3.1796946 19425.54582
SO (27.5,27.75] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 1821 1.6789480 7068.9398858 9.6484521 2000-2012 –> 2013-2019 3.1518168 19332.25293
SO (27.75,27.85] 4 cstar_tref ~ aou + nitrate + silicate + phosphate_star 1415 1.4036591 4987.1994177 11.2501595 1982-1999 –> 2000-2012 2.6328581 10866.53706
SO (27.75,27.85] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1415 1.3799763 4941.0437367 10.6092986 1982-1999 –> 2000-2012 2.6003361 10796.32800
SO (27.75,27.85] 4 cstar_tref ~ temp + aou + nitrate + silicate 1415 1.5056329 5185.6697430 13.4580336 1982-1999 –> 2000-2012 2.7650010 11152.54027
SO (27.75,27.85] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1415 1.3062806 4785.7263346 8.6276320 1982-1999 –> 2000-2012 2.5318731 10656.45667
SO (27.75,27.85] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 1415 1.2832711 4735.4330180 8.9326433 1982-1999 –> 2000-2012 2.5044896 10593.25645
SO (27.75,27.85] 4 cstar_tref ~ aou + nitrate + silicate + phosphate_star 722 1.5404965 2684.9065454 8.3689795 2000-2012 –> 2013-2019 2.9441556 7672.10596
SO (27.75,27.85] 4 cstar_tref ~ sal + aou + nitrate + silicate 722 1.5764019 2718.1765705 8.5554720 2000-2012 –> 2013-2019 3.0324771 7809.12992
SO (27.75,27.85] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 722 1.5043391 2652.6099150 8.2718276 2000-2012 –> 2013-2019 2.8843154 7593.65365
SO (27.75,27.85] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 722 1.4200481 2569.3447020 8.5861197 2000-2012 –> 2013-2019 2.7263287 7355.07104
SO (27.75,27.85] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 722 1.3849071 2533.1613767 8.2492850 2000-2012 –> 2013-2019 2.6681782 7268.59439
SO (27.85,27.95] 4 cstar_tref ~ aou + nitrate + silicate + phosphate_star 1399 1.5026265 5121.5763882 9.8099973 1982-1999 –> 2000-2012 3.0891268 12102.59481
SO (27.85,27.95] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1399 1.3661540 4857.1645441 8.5867618 1982-1999 –> 2000-2012 2.7006066 11199.01016
SO (27.85,27.95] 4 cstar_tref ~ temp + aou + nitrate + silicate 1399 1.5200790 5153.8870333 9.6302273 1982-1999 –> 2000-2012 3.0894478 12094.66913
SO (27.85,27.95] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1399 1.5025970 5123.5214924 9.8177631 1982-1999 –> 2000-2012 3.0719532 12066.27393
SO (27.85,27.95] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 1399 1.5146295 5145.8381005 9.8001556 1982-1999 –> 2000-2012 3.0896307 12101.89705
SO (27.85,27.95] 4 cstar_tref ~ aou + nitrate + silicate + phosphate_star 709 1.5653688 2659.4910609 8.8715048 2000-2012 –> 2013-2019 3.0679953 7781.06745
SO (27.85,27.95] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 709 1.5223861 2622.0102956 7.8651764 2000-2012 –> 2013-2019 2.8885400 7479.17484
SO (27.85,27.95] 4 cstar_tref ~ temp + aou + nitrate + silicate 709 1.6394437 2725.0530335 9.1707714 2000-2012 –> 2013-2019 3.1595227 7878.94007
SO (27.85,27.95] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 709 1.5516698 2649.0271053 9.0966674 2000-2012 –> 2013-2019 3.0542668 7772.54860
SO (27.85,27.95] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 709 1.5648116 2660.9862234 9.0697766 2000-2012 –> 2013-2019 3.0794411 7806.82432
SO (27.95,28.05] 4 cstar_tref ~ aou + nitrate + silicate + phosphate_star 1786 1.7878376 7155.8048596 9.9981828 1982-1999 –> 2000-2012 3.5042943 16912.86903
SO (27.95,28.05] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 1786 1.6882407 6951.0583598 9.4364874 1982-1999 –> 2000-2012 3.1710443 15980.28933
SO (27.95,28.05] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1786 1.5483054 6643.9876245 7.5549297 1982-1999 –> 2000-2012 2.8861036 15163.34679
SO (27.95,28.05] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 1786 1.8559804 7291.4200292 8.9354048 1982-1999 –> 2000-2012 3.5003230 16836.99264
SO (27.95,28.05] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1786 1.7862290 7154.5895927 10.1233100 1982-1999 –> 2000-2012 3.4995003 16904.41411
SO (27.95,28.05] 4 cstar_tref ~ aou + nitrate + silicate + phosphate_star 986 1.9008717 4076.7871557 8.6940193 2000-2012 –> 2013-2019 3.6887093 11232.59202
SO (27.95,28.05] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 986 1.9327722 4109.6067463 9.1048682 2000-2012 –> 2013-2019 3.6210129 11060.66511
SO (27.95,28.05] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 986 1.6648161 3817.3040710 7.0125027 2000-2012 –> 2013-2019 3.2131214 10461.29170
SO (27.95,28.05] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 986 1.8908856 4068.4000974 8.6762070 2000-2012 –> 2013-2019 3.6771146 11222.98969
SO (27.95,28.05] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 986 1.8962543 4073.9911786 8.8729848 2000-2012 –> 2013-2019 3.6863263 11236.25753
SO (28.05,28.1] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 1071 1.2084510 3456.9312497 6.8893057 1982-1999 –> 2000-2012 2.3093792 8056.63029
SO (28.05,28.1] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1071 1.1155347 3287.5593365 8.1328311 1982-1999 –> 2000-2012 2.1524995 7708.01563
SO (28.05,28.1] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 1071 1.2795987 3579.4689177 10.2387822 1982-1999 –> 2000-2012 2.4750477 8428.57908
SO (28.05,28.1] 4 cstar_tref ~ sal + nitrate + silicate + phosphate_star 1071 1.2890735 3595.2709915 9.8929533 1982-1999 –> 2000-2012 2.4900398 8458.32396
SO (28.05,28.1] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 1071 1.2416885 3517.0497248 7.9412274 1982-1999 –> 2000-2012 2.4228048 8331.63650
SO (28.05,28.1] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 606 1.4409135 2174.4695378 7.1949964 2000-2012 –> 2013-2019 2.6493644 5631.40079
SO (28.05,28.1] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 606 1.3483022 2095.9550940 8.8409161 2000-2012 –> 2013-2019 2.4638369 5383.51443
SO (28.05,28.1] 4 cstar_tref ~ sal + nitrate + silicate + phosphate_star 606 1.5082442 2229.8202625 11.4315183 2000-2012 –> 2013-2019 2.7973177 5825.09125
SO (28.05,28.1] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 606 1.4874983 2215.0335236 9.9011290 2000-2012 –> 2013-2019 2.7291868 5732.08325
SO (28.05,28.1] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 606 1.5143347 2236.7046774 9.5724758 2000-2012 –> 2013-2019 2.8015603 5830.90291
SO (28.1, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 5309 1.0966504 16059.9105945 6.0094365 1982-1999 –> 2000-2012 2.0858171 36797.53418
SO (28.1, Inf] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 5309 1.1018770 16108.3952774 6.0859562 1982-1999 –> 2000-2012 2.0912319 36846.81840
SO (28.1, Inf] 4 cstar_tref ~ sal + nitrate + silicate + phosphate_star 5309 1.0966518 16057.9241425 6.0074321 1982-1999 –> 2000-2012 2.0897732 36852.27237
SO (28.1, Inf] 4 cstar_tref ~ sal + oxygen + nitrate + silicate 5309 1.0921155 16013.9119813 6.1402836 1982-1999 –> 2000-2012 2.0772669 36689.66808
SO (28.1, Inf] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 5309 1.0921105 16015.8628575 6.1368799 1982-1999 –> 2000-2012 2.0724357 36621.34083
SO (28.1, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 2711 1.1833450 8620.2523271 6.1959706 2000-2012 –> 2013-2019 2.2799954 24680.16292
SO (28.1, Inf] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 2711 1.1946469 8669.7908607 6.0803943 2000-2012 –> 2013-2019 2.2965238 24778.18614
SO (28.1, Inf] 4 cstar_tref ~ sal + nitrate + silicate + phosphate_star 2711 1.1852384 8626.9207957 6.1789071 2000-2012 –> 2013-2019 2.2818902 24684.84494
SO (28.1, Inf] 4 cstar_tref ~ sal + oxygen + nitrate + silicate 2711 1.1836919 8619.8415600 6.3248519 2000-2012 –> 2013-2019 2.2758074 24633.75354
SO (28.1, Inf] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 2711 1.1827113 8617.3478342 6.2791456 2000-2012 –> 2013-2019 2.2748217 24633.21069
SO (28.1,28.15] 4 cstar_tref ~ sal + aou + nitrate + silicate 395 0.5288565 629.7012772 1.8285447 1982-1999 –> 2000-2012 1.0213642 1382.24384
SO (28.1,28.15] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 395 0.4829853 560.0239392 1.5318770 1982-1999 –> 2000-2012 0.9517188 1263.01268
SO (28.1,28.15] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 395 0.5103244 601.5215519 1.7907538 1982-1999 –> 2000-2012 0.9951592 1337.70278
SO (28.1,28.15] 4 cstar_tref ~ sal + oxygen + nitrate + phosphate_star 395 0.5650355 681.9767823 3.4016063 1982-1999 –> 2000-2012 1.0255225 1364.47026
SO (28.1,28.15] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 395 0.5157574 611.8875562 2.4216515 1982-1999 –> 2000-2012 0.9669854 1275.21630
SO (28.1,28.15] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 209 0.5031920 320.0407883 2.0966657 2000-2012 –> 2013-2019 0.9861773 880.06473
SO (28.1,28.15] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 209 0.5573440 360.7649404 2.1406620 2000-2012 –> 2013-2019 1.0676684 962.28649
SO (28.1,28.15] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 209 0.5670742 369.9995113 2.4309188 2000-2012 –> 2013-2019 1.0828316 981.88707
SO (28.1,28.15] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 209 0.5580602 363.3017219 2.6731375 2000-2012 –> 2013-2019 1.0983891 1011.95719
SO (28.1,28.15] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 209 0.5556540 361.4955364 2.6910012 2000-2012 –> 2013-2019 1.0930170 1005.80272
SO (28.15,28.2] 4 cstar_tref ~ oxygen + nitrate + silicate + phosphate_star 385 0.3253916 240.0837885 1.3423359 1982-1999 –> 2000-2012 0.6085395 406.40761
SO (28.15,28.2] 4 cstar_tref ~ sal + oxygen + nitrate + phosphate_star 385 0.3164134 218.5393499 1.2843448 1982-1999 –> 2000-2012 0.6055513 405.42080
SO (28.15,28.2] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 385 0.3164122 220.5363856 1.2824638 1982-1999 –> 2000-2012 0.5891766 352.17177
SO (28.15,28.2] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 385 0.3073698 198.2108824 1.1858788 1982-1999 –> 2000-2012 0.5773798 319.87950
SO (28.15,28.2] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 385 0.3038928 189.4509187 1.1855544 1982-1999 –> 2000-2012 0.5687874 292.33674
SO (28.15,28.2] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 183 0.3470133 145.9599817 0.8256096 2000-2012 –> 2013-2019 0.6543831 344.17086
SO (28.15,28.2] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 183 0.3477757 144.7632517 0.8263895 2000-2012 –> 2013-2019 0.6679819 372.47753
SO (28.15,28.2] 4 cstar_tref ~ temp + oxygen + nitrate + phosphate_star 183 0.3582778 155.6520285 0.8521815 2000-2012 –> 2013-2019 0.6638681 347.39203
SO (28.15,28.2] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 183 0.3456801 144.5510795 0.8275075 2000-2012 –> 2013-2019 0.6495729 334.00200
SO (28.15,28.2] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 183 0.3471746 144.1300076 0.8242149 2000-2012 –> 2013-2019 0.6668994 370.68584
SO (28.2, Inf] 3 cstar_tref ~ sal + aou + nitrate 1738 0.3752615 1535.2908035 2.0011297 1982-1999 –> 2000-2012 0.6790183 2669.60844
SO (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 1738 0.3665194 1455.3553954 1.8883521 1982-1999 –> 2000-2012 0.6702525 2591.28860
SO (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + silicate 1738 0.3672295 1462.0838414 1.9673041 1982-1999 –> 2000-2012 0.6709769 2598.24970
SO (28.2, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1738 0.3664835 1457.0151862 1.8977416 1982-1999 –> 2000-2012 0.6701959 2594.61023
SO (28.2, Inf] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1738 0.4262613 1982.2357759 1.9959611 1982-1999 –> 2000-2012 0.7640018 3644.86728
SO (28.2, Inf] 3 cstar_tref ~ sal + aou + nitrate 1005 0.4164558 1101.3567107 2.5953775 2000-2012 –> 2013-2019 0.7917173 2636.64751
SO (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 1005 0.3988626 1016.5985475 2.4052531 2000-2012 –> 2013-2019 0.7653820 2471.95394
SO (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + silicate 1005 0.3952080 998.0966826 2.5161713 2000-2012 –> 2013-2019 0.7624375 2460.18052
SO (28.2, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1005 0.3950006 999.0416755 2.5557103 2000-2012 –> 2013-2019 0.7614841 2456.05686
SO (28.2, Inf] 5 cstar_tref ~ sal + oxygen + nitrate + silicate + phosphate_star 1005 0.4448838 1238.0826928 3.0898823 2000-2012 –> 2013-2019 0.8786152 3280.70645

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 3.2926584 1431.34564
Atlantic (-Inf,26] 2000-2012 –> 2013-2019 3.3070190 904.58908
Atlantic (26,26.5] 1982-1999 –> 2000-2012 8.0897552 14401.65612
Atlantic (26,26.5] 2000-2012 –> 2013-2019 7.8905913 8990.12921
Atlantic (26.5,26.75] 1982-1999 –> 2000-2012 6.6255573 16722.71012
Atlantic (26.5,26.75] 2000-2012 –> 2013-2019 6.8915508 10672.86522
Atlantic (26.75,27] 1982-1999 –> 2000-2012 4.2477333 20978.85070
Atlantic (26.75,27] 2000-2012 –> 2013-2019 5.0705417 15153.88078
Atlantic (27,27.25] 1982-1999 –> 2000-2012 4.1577170 16982.57391
Atlantic (27,27.25] 2000-2012 –> 2013-2019 4.6590792 12160.78879
Atlantic (27.25,27.5] 1982-1999 –> 2000-2012 4.1863390 16463.58533
Atlantic (27.25,27.5] 2000-2012 –> 2013-2019 4.8587725 11563.05069
Atlantic (27.5,27.75] 1982-1999 –> 2000-2012 5.2131781 23173.81547
Atlantic (27.5,27.75] 2000-2012 –> 2013-2019 5.5667786 15535.78753
Atlantic (27.75,27.85] 1982-1999 –> 2000-2012 3.6067669 6786.48436
Atlantic (27.75,27.85] 2000-2012 –> 2013-2019 3.6147375 4363.50340
Atlantic (27.85,27.95] 1982-1999 –> 2000-2012 5.9687977 8342.14838
Atlantic (27.85,27.95] 2000-2012 –> 2013-2019 6.1362141 5646.67684
Atlantic (27.95,28.05] 1982-1999 –> 2000-2012 5.1147377 8073.66170
Atlantic (27.95,28.05] 2000-2012 –> 2013-2019 4.6919705 5128.66866
Atlantic (28.05,28.1] 1982-1999 –> 2000-2012 0.3303716 -945.92797
Atlantic (28.05,28.1] 2000-2012 –> 2013-2019 0.4603548 -82.59284
Atlantic (28.1,28.15] 1982-1999 –> 2000-2012 0.2892629 -1398.29216
Atlantic (28.1,28.15] 2000-2012 –> 2013-2019 0.3262964 -641.63511
Atlantic (28.15,28.2] 1982-1999 –> 2000-2012 0.3528182 -1713.29336
Atlantic (28.15,28.2] 2000-2012 –> 2013-2019 0.3613690 -1221.21112
Atlantic (28.2, Inf] 1982-1999 –> 2000-2012 0.4076206 -722.35839
Atlantic (28.2, Inf] 2000-2012 –> 2013-2019 0.3721221 -829.73000
Indo-Pacific (-Inf,26] 1982-1999 –> 2000-2012 15.0227518 85839.02548
Indo-Pacific (-Inf,26] 2000-2012 –> 2013-2019 14.4552325 55148.00975
Indo-Pacific (26,26.5] 1982-1999 –> 2000-2012 9.5090421 77407.63536
Indo-Pacific (26,26.5] 2000-2012 –> 2013-2019 9.7059589 49889.72682
Indo-Pacific (26.5,26.75] 1982-1999 –> 2000-2012 8.0977884 57470.74360
Indo-Pacific (26.5,26.75] 2000-2012 –> 2013-2019 8.4382438 37564.87851
Indo-Pacific (26.75,27] 1982-1999 –> 2000-2012 9.1206948 64287.44413
Indo-Pacific (26.75,27] 2000-2012 –> 2013-2019 10.0612953 42427.90610
Indo-Pacific (27,27.25] 1982-1999 –> 2000-2012 7.9804670 55766.51372
Indo-Pacific (27,27.25] 2000-2012 –> 2013-2019 9.8837794 39559.34089
Indo-Pacific (27.25,27.5] 1982-1999 –> 2000-2012 5.9324677 41257.66268
Indo-Pacific (27.25,27.5] 2000-2012 –> 2013-2019 6.5624974 27893.00253
Indo-Pacific (27.5,27.75] 1982-1999 –> 2000-2012 3.5551605 30394.50069
Indo-Pacific (27.5,27.75] 2000-2012 –> 2013-2019 4.1431709 20594.94023
Indo-Pacific (27.75,27.85] 1982-1999 –> 2000-2012 2.7915689 10163.57755
Indo-Pacific (27.75,27.85] 2000-2012 –> 2013-2019 3.0883766 7023.64294
Indo-Pacific (27.85,27.95] 1982-1999 –> 2000-2012 2.6307305 14607.92853
Indo-Pacific (27.85,27.95] 2000-2012 –> 2013-2019 2.9831363 10006.00610
Indo-Pacific (27.95,28.05] 1982-1999 –> 2000-2012 1.2753675 5578.23027
Indo-Pacific (27.95,28.05] 2000-2012 –> 2013-2019 1.3487769 3730.20895
Indo-Pacific (28.05,28.1] 1982-1999 –> 2000-2012 0.9615446 3491.79312
Indo-Pacific (28.05,28.1] 2000-2012 –> 2013-2019 0.9729832 2237.67602
Indo-Pacific (28.1, Inf] 1982-1999 –> 2000-2012 1.4199089 47936.97932
Indo-Pacific (28.1, Inf] 2000-2012 –> 2013-2019 1.4210822 31029.03176
SO (26,26.5] 1982-1999 –> 2000-2012 2.8169627 827.56855
SO (26,26.5] 2000-2012 –> 2013-2019 2.5341197 506.77627
SO (26.5,26.75] 1982-1999 –> 2000-2012 4.0066866 3116.58566
SO (26.5,26.75] 2000-2012 –> 2013-2019 3.6826837 1831.61801
SO (26.75,27] 1982-1999 –> 2000-2012 5.4008031 12475.00915
SO (26.75,27] 2000-2012 –> 2013-2019 5.2986572 7711.59168
SO (27,27.25] 1982-1999 –> 2000-2012 3.7595457 24396.70402
SO (27,27.25] 2000-2012 –> 2013-2019 3.9971617 15944.80466
SO (27.25,27.5] 1982-1999 –> 2000-2012 3.1263315 24201.09053
SO (27.25,27.5] 2000-2012 –> 2013-2019 3.5476982 16557.69118
SO (27.5,27.75] 1982-1999 –> 2000-2012 2.7510826 27349.51837
SO (27.5,27.75] 2000-2012 –> 2013-2019 3.1944302 19461.36098
SO (27.75,27.85] 1982-1999 –> 2000-2012 2.6069116 10813.02369
SO (27.75,27.85] 2000-2012 –> 2013-2019 2.8510910 7539.71099
SO (27.85,27.95] 1982-1999 –> 2000-2012 3.0081530 11912.88902
SO (27.85,27.95] 2000-2012 –> 2013-2019 3.0499532 7743.71106
SO (27.95,28.05] 1982-1999 –> 2000-2012 3.3122531 16359.58238
SO (27.95,28.05] 2000-2012 –> 2013-2019 3.5772569 11042.75921
SO (28.05,28.1] 1982-1999 –> 2000-2012 2.3699542 8196.63709
SO (28.05,28.1] 2000-2012 –> 2013-2019 2.6882532 5680.59853
SO (28.1, Inf] 1982-1999 –> 2000-2012 2.0833050 36761.52677
SO (28.1, Inf] 2000-2012 –> 2013-2019 2.2818077 24682.03165
SO (28.1,28.15] 1982-1999 –> 2000-2012 0.9921500 1324.52917
SO (28.1,28.15] 2000-2012 –> 2013-2019 1.0656167 968.39964
SO (28.15,28.2] 1982-1999 –> 2000-2012 0.5898869 355.24329
SO (28.15,28.2] 2000-2012 –> 2013-2019 0.6605411 353.74565
SO (28.2, Inf] 1982-1999 –> 2000-2012 0.6908891 2819.72485
SO (28.2, Inf] 2000-2012 –> 2013-2019 0.7919272 2661.10906

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
94ce5a8 Donghe-Zhu 2021-03-10
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
94ce5a8 Donghe-Zhu 2021-03-10
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
94ce5a8 Donghe-Zhu 2021-03-10
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
94ce5a8 Donghe-Zhu 2021-03-10
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 nitrate oxygen phosphate_star sal silicate temp
Atlantic - 1982-1999 --> 2000-2012
(-Inf,26] 3 4 2 4 5 4 0
(26,26.5] 2 5 3 3 0 4 5
(26.5,26.75] 3 3 2 5 1 5 4
(26.75,27] 3 3 2 5 1 5 4
(27,27.25] 3 5 2 3 1 5 4
(27.25,27.5] 4 5 1 3 0 3 5
(27.5,27.75] 4 5 1 3 5 3 0
(27.75,27.85] 5 5 0 3 4 3 1
(27.85,27.95] 3 3 2 4 0 5 5
(27.95,28.05] 2 5 2 5 0 3 5
(28.05,28.1] 2 3 2 5 0 5 5
(28.1,28.15] 2 5 2 5 0 3 5
(28.15,28.2] 4 5 1 3 5 3 0
(28.2, Inf] 3 5 2 4 5 3 0
total 43.00 61.00 24.00 55.00 27.00 54.00 43.00
Atlantic - 2000-2012 --> 2013-2019
(-Inf,26] 3 4 2 4 5 4 0
(26,26.5] 2 5 3 3 0 4 5
(26.5,26.75] 4 3 1 5 2 5 3
(26.75,27] 3 4 2 5 3 5 2
(27,27.25] 3 3 2 5 4 5 1
(27.25,27.5] 5 5 0 3 4 3 1
(27.5,27.75] 2 4 2 5 5 4 0
(27.75,27.85] 5 5 0 3 4 3 1
(27.85,27.95] 3 3 2 4 0 5 5
(27.95,28.05] 2 5 2 5 0 3 5
(28.05,28.1] 2 5 2 5 0 3 5
(28.1,28.15] 2 4 3 5 0 3 5
(28.15,28.2] 5 5 0 3 4 3 1
(28.2, Inf] 3 5 2 4 5 3 0
total 44.00 60.00 23.00 59.00 36.00 53.00 34.00
Indo-Pacific - 1982-1999 --> 2000-2012
(-Inf,26] 2 3 3 5 0 4 5
(26,26.5] 2 3 3 5 0 4 5
(26.5,26.75] 2 5 2 5 0 3 5
(26.75,27] 2 4 3 5 0 3 5
(27,27.25] 2 5 2 5 0 3 5
(27.25,27.5] 2 4 3 5 3 5 2
(27.5,27.75] 2 3 3 5 4 5 1
(27.75,27.85] 3 3 2 5 1 5 4
(27.85,27.95] 3 5 2 3 1 5 4
(27.95,28.05] 2 3 2 5 0 5 5
(28.05,28.1] 2 3 2 5 5 5 0
(28.1, Inf] 2 5 3 5 2 5 2
total 26.00 46.00 30.00 58.00 16.00 52.00 43.00
Indo-Pacific - 2000-2012 --> 2013-2019
(-Inf,26] 2 3 3 5 0 4 5
(26,26.5] 2 3 2 5 0 5 5
(26.5,26.75] 2 5 2 5 0 3 5
(26.75,27] 2 4 3 5 0 3 5
(27,27.25] 2 3 3 5 1 5 4
(27.25,27.5] 2 4 3 5 3 5 2
(27.5,27.75] 2 3 3 5 4 5 1
(27.75,27.85] 3 3 2 5 1 5 4
(27.85,27.95] 3 5 2 3 1 5 4
(27.95,28.05] 2 3 2 5 0 5 5
(28.05,28.1] 2 5 2 5 3 5 2
(28.1, Inf] 2 5 3 5 2 5 2
total 26.00 46.00 30.00 58.00 15.00 55.00 44.00
SO - 1982-1999 --> 2000-2012
(26,26.5] 2 5 3 5 1 3 4
(26.5,26.75] 2 5 3 5 1 3 4
(26.75,27] 3 4 2 5 0 3 5
(27,27.25] 3 4 2 5 0 3 5
(27.25,27.5] 2 3 2 5 0 5 5
(27.5,27.75] 4 5 1 4 1 5 3
(27.75,27.85] 4 5 1 4 1 5 3
(27.85,27.95] 4 5 1 4 1 5 3
(27.95,28.05] 4 5 1 5 3 4 1
(28.05,28.1] 3 4 1 5 5 4 0
(28.1, Inf] 2 4 2 4 5 5 0
(28.1,28.15] 3 4 2 4 5 4 0
(28.15,28.2] 1 5 4 5 2 4 2
(28.2, Inf] 5 5 0 3 4 3 1
total 42.00 63.00 25.00 63.00 29.00 56.00 36.00
SO - 2000-2012 --> 2013-2019
(26,26.5] 2 5 3 5 1 3 4
(26.5,26.75] 2 5 3 5 1 3 4
(26.75,27] 3 4 2 5 0 3 5
(27,27.25] 3 4 2 5 0 3 5
(27.25,27.5] 2 3 2 5 0 5 5
(27.5,27.75] 4 5 1 4 2 5 2
(27.75,27.85] 4 5 1 4 2 5 2
(27.85,27.95] 4 5 1 4 1 5 3
(27.95,28.05] 4 5 1 5 2 4 2
(28.05,28.1] 2 5 2 5 4 4 1
(28.1, Inf] 2 4 2 4 5 5 0
(28.1,28.15] 3 4 2 5 3 5 2
(28.15,28.2] 2 3 3 5 0 4 5
(28.2, Inf] 4 5 1 3 5 3 0
total 41.00 62.00 26.00 64.00 26.00 57.00 40.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
94ce5a8 Donghe-Zhu 2021-03-10
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
94ce5a8 Donghe-Zhu 2021-03-10
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