Last updated: 2021-03-27

Checks: 7 0

Knit directory: emlr_mod_v_XXX/

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

R Markdown file: up-to-date

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

Environment: empty

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

Seed: set.seed(20200707)

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

Session information: recorded

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

Cache: none

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

File paths: relative

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

Repository version: a52ad87

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 a52ad87. 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/

Untracked files:
    Untracked:  2_Pac_zon_em.png
    Untracked:  so_Pac_zon_off.png

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

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

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

File Version Author Date Message
html e2704ca Donghe-Zhu 2021-03-27 Build site.
Rmd b1706de Donghe-Zhu 2021-03-27 local rebuild
html 3061a0b Donghe-Zhu 2021-03-27 Build site.
html b883157 Donghe-Zhu 2021-03-27 Build site.
html 0c20513 Donghe-Zhu 2021-03-26 Build site.
html d19654d Donghe-Zhu 2021-03-26 Build site.
Rmd 921c9db Donghe-Zhu 2021-03-26 local rebuild
html 6c53dbf Donghe-Zhu 2021-03-25 Build site.
Rmd 195db95 Donghe-Zhu 2021-03-24 local rebuild
html 1914a11 Donghe-Zhu 2021-03-24 Build site.
Rmd e32af14 Donghe-Zhu 2021-03-24 local rebuild
html 8be810e Donghe-Zhu 2021-03-23 Build site.
Rmd 85274db Donghe-Zhu 2021-03-23 local rebuild
html bf19764 Donghe-Zhu 2021-03-22 Build site.
html 3ec9d3d Donghe-Zhu 2021-03-22 Build site.
Rmd 073fa6a Donghe-Zhu 2021-03-22 local rebuild
html 134ace1 Donghe-Zhu 2021-03-22 Build site.
html f6d70a4 Donghe-Zhu 2021-03-22 Build site.
html 183443b Donghe-Zhu 2021-03-21 Build site.
Rmd adca794 Donghe-Zhu 2021-03-21 local rebuild
html 2e6976b Donghe-Zhu 2021-03-21 Build site.
Rmd 4e7de35 Donghe-Zhu 2021-03-21 local rebuild
html 51a42bd Donghe-Zhu 2021-03-16 Build site.
html f745381 Donghe-Zhu 2021-03-16 Build site.
html aecbf75 Donghe-Zhu 2021-03-14 Build site.
html 1b2a0c2 Donghe-Zhu 2021-03-14 Build site.
html 6733e48 Donghe-Zhu 2021-03-12 Build site.
html ba71e6a Donghe-Zhu 2021-03-12 Build site.
html 9dbf5bd Donghe-Zhu 2021-03-11 Build site.
html a49df30 Donghe-Zhu 2021-03-11 Build site.
html b3348a5 Donghe-Zhu 2021-03-11 Build site.
Rmd 026c061 Donghe-Zhu 2021-03-11 local rebuild
html 1c24ff7 Donghe-Zhu 2021-03-10 Build site.
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, aou, nitrate, silicate, phosphate, phosphate_star
# the following code is a workaround to find all predictor combinations
# using the olsrr package and fit all models for one era, slab, and basin

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

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

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

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

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

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

# convert to tibble
lm_all <- as_tibble(lm_all)

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

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

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

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

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

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

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

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

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

3 Apply predictor threshold

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

  • Minimum: 2
  • Maximum: 7
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)[100]
        
        # 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)
          
          # calculate max VIF
          if (is.na(max(i_lm_fit$coefficients))){
            i_vif_max <- as.double("NA")
          } else{
          i_vif_max <- max(vif(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,
              vif_max = i_vif_max,
              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_vif_max,
  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,
    vif_max,
    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)

5.3 VIF threshold

To avoid multicollinearity among predictors, models were excluded with a VIF above:

  • 10

After removing models affected by multicollinearity, the targeted number of MLRs (5) was undercut in following fitting units:

# remove models with predictors fitted as NA
lm_all_fitted_wide <- lm_all_fitted_wide %>%
  filter(vif_max <= params_local$vif_max)

lm_all_fitted_wide_check <- lm_all_fitted_wide %>%
  group_by(era, basin, gamma_slab) %>% 
  count() %>% 
  filter(n < params_local$MLR_number)
# 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() %>%
  filter(eras != "1982-1999 --> 2010-2019")

# 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.4 RMSE tables

5.4.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] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 123 1.5240065 466.70920 5.7925318 1982-1999 –> 2000-2009 3.2097663 1288.22963
Atlantic (-Inf,26] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 123 1.0373298 370.07477 4.0246479 1982-1999 –> 2000-2009 2.3736782 1092.96858
Atlantic (-Inf,26] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 123 1.0372836 372.06381 4.0109558 1982-1999 –> 2000-2009 2.3729762 1096.75340
Atlantic (-Inf,26] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 123 1.4327537 449.52006 6.5328965 1982-1999 –> 2000-2009 3.0950239 1263.20313
Atlantic (-Inf,26] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 123 1.4210216 449.49739 6.1762663 1982-1999 –> 2000-2009 3.0825101 1264.98479
Atlantic (-Inf,26] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 107 1.7005398 431.27523 3.9403910 2000-2009 –> 2010-2019 3.2245463 897.98443
Atlantic (-Inf,26] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 107 1.2078500 356.06502 3.5405245 2000-2009 –> 2010-2019 2.2451798 726.13979
Atlantic (-Inf,26] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 107 1.1371667 345.16037 3.1990511 2000-2009 –> 2010-2019 2.1744503 717.22418
Atlantic (-Inf,26] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 107 1.5675559 411.84963 6.2388419 2000-2009 –> 2010-2019 3.0003097 861.36969
Atlantic (-Inf,26] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 107 1.5143272 406.45670 5.5944269 2000-2009 –> 2010-2019 2.9353488 855.95409
Atlantic (26,26.5] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 815 3.5320379 4383.72607 13.5273360 1982-1999 –> 2000-2009 6.9872164 12437.97554
Atlantic (26,26.5] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 815 4.2289569 4675.25704 10.4983456 1982-1999 –> 2000-2009 8.3828315 13284.42881
Atlantic (26,26.5] 4 cstar_tref ~ temp + aou + nitrate + silicate 815 4.0306894 4596.98782 12.3130916 1982-1999 –> 2000-2009 8.0814731 13130.16288
Atlantic (26,26.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 815 4.0075189 4589.59069 13.2031304 1982-1999 –> 2000-2009 8.0582990 13124.76307
Atlantic (26,26.5] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 815 4.2984157 4703.81163 13.3440553 1982-1999 –> 2000-2009 8.5456468 13382.19396
Atlantic (26,26.5] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 787 3.6014952 4264.25272 11.3693942 2000-2009 –> 2010-2019 7.1335331 8647.97879
Atlantic (26,26.5] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 787 4.0998535 4466.24651 12.1009472 2000-2009 –> 2010-2019 8.3288104 9141.50355
Atlantic (26,26.5] 4 cstar_tref ~ temp + aou + nitrate + silicate 787 3.9028825 4388.74925 15.9584819 2000-2009 –> 2010-2019 7.9335719 8985.73707
Atlantic (26,26.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 787 3.8638909 4374.94518 18.4031563 2000-2009 –> 2010-2019 7.8714097 8964.53587
Atlantic (26,26.5] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 787 4.2226500 4514.69785 11.4755480 2000-2009 –> 2010-2019 8.5210657 9218.50948
Atlantic (26.5,26.75] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 1030 3.2179610 5342.63412 10.2321299 1982-1999 –> 2000-2009 6.3324329 14741.69730
Atlantic (26.5,26.75] 4 cstar_tref ~ sal + aou + silicate + phosphate 1030 3.2133679 5339.69170 10.0941448 1982-1999 –> 2000-2009 6.3191433 14728.48169
Atlantic (26.5,26.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1030 3.1434743 5296.39042 10.3581653 1982-1999 –> 2000-2009 6.1946657 14622.03568
Atlantic (26.5,26.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1030 3.2235806 5348.22840 9.8461989 1982-1999 –> 2000-2009 6.3766485 14794.54172
Atlantic (26.5,26.75] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1030 3.0858055 5258.24758 9.3381923 1982-1999 –> 2000-2009 6.0681722 14500.07075
Atlantic (26.5,26.75] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 976 3.6231389 5294.65718 11.5505370 2000-2009 –> 2010-2019 6.8410999 10637.29130
Atlantic (26.5,26.75] 4 cstar_tref ~ sal + aou + silicate + phosphate 976 3.5199757 5238.27039 11.2426091 2000-2009 –> 2010-2019 6.7333436 10577.96208
Atlantic (26.5,26.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 976 3.5053683 5232.15303 11.4141756 2000-2009 –> 2010-2019 6.6488427 10528.54345
Atlantic (26.5,26.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 976 3.5220570 5241.42422 9.2200905 2000-2009 –> 2010-2019 6.7456376 10589.65262
Atlantic (26.5,26.75] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 976 3.5013623 5229.92098 10.8565185 2000-2009 –> 2010-2019 6.5871678 10488.16856
Atlantic (26.75,27] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 1620 2.1299666 7059.14524 14.9056340 1982-1999 –> 2000-2009 3.8345830 17718.87300
Atlantic (26.75,27] 4 cstar_tref ~ sal + aou + silicate + phosphate 1620 2.1107392 7029.76464 16.9316002 1982-1999 –> 2000-2009 3.9534249 18114.27194
Atlantic (26.75,27] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1620 1.9837680 6830.75452 13.8624424 1982-1999 –> 2000-2009 3.6371854 17326.15807
Atlantic (26.75,27] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1620 2.2283003 7207.37553 15.0500811 1982-1999 –> 2000-2009 4.0604118 18262.49488
Atlantic (26.75,27] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1620 2.1293496 7060.20650 15.0247517 1982-1999 –> 2000-2009 3.6214756 16995.79883
Atlantic (26.75,27] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 1635 2.5945390 7769.57591 27.8116015 2000-2009 –> 2010-2019 4.7245056 14828.72114
Atlantic (26.75,27] 4 cstar_tref ~ sal + aou + silicate + phosphate 1635 2.4410726 7570.19974 24.6840365 2000-2009 –> 2010-2019 4.5518118 14599.96439
Atlantic (26.75,27] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1635 2.4271146 7553.44833 25.3806732 2000-2009 –> 2010-2019 4.4108826 14384.20285
Atlantic (26.75,27] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1635 2.6330315 7819.73326 23.4947616 2000-2009 –> 2010-2019 4.8613318 15027.10879
Atlantic (26.75,27] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1635 2.1699807 7187.25780 21.9087669 2000-2009 –> 2010-2019 4.2993303 14247.46429
Atlantic (27,27.25] 3 cstar_tref ~ sal + aou + phosphate 1310 2.3623438 5979.91316 11.5636335 1982-1999 –> 2000-2009 4.2272160 15150.91309
Atlantic (27,27.25] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 1310 2.3309138 5946.82116 12.3663950 1982-1999 –> 2000-2009 4.1833369 15089.77397
Atlantic (27,27.25] 4 cstar_tref ~ sal + aou + silicate + phosphate 1310 1.9895142 5531.89205 9.1174848 1982-1999 –> 2000-2009 3.7091375 14341.13433
Atlantic (27,27.25] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1310 1.9729970 5512.04966 9.4699382 1982-1999 –> 2000-2009 3.5841308 14030.95240
Atlantic (27,27.25] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1310 2.1082203 5685.73055 13.1497941 1982-1999 –> 2000-2009 3.8419338 14533.58065
Atlantic (27,27.25] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1389 2.2670798 6229.58358 10.9186762 2000-2009 –> 2010-2019 4.3556942 11890.83487
Atlantic (27,27.25] 4 cstar_tref ~ sal + aou + silicate + phosphate 1389 2.3218868 6293.94328 13.6907831 2000-2009 –> 2010-2019 4.3114010 11825.83533
Atlantic (27,27.25] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1389 2.3217985 6295.83758 13.7217883 2000-2009 –> 2010-2019 4.2947955 11807.88725
Atlantic (27,27.25] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 1389 2.3962856 6381.56059 14.8354885 2000-2009 –> 2010-2019 4.5794880 12156.85673
Atlantic (27,27.25] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1389 2.5106188 6513.04156 18.5880756 2000-2009 –> 2010-2019 4.6188391 12198.77211
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1275 2.2441459 5693.52200 13.1089562 1982-1999 –> 2000-2009 4.4555246 15659.59306
Atlantic (27.25,27.5] 3 cstar_tref ~ temp + aou + nitrate 1275 2.3897095 5849.78141 15.3256992 1982-1999 –> 2000-2009 4.2337940 14994.86625
Atlantic (27.25,27.5] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 1275 2.3648833 5825.15140 16.8005151 1982-1999 –> 2000-2009 4.1933453 14933.96814
Atlantic (27.25,27.5] 4 cstar_tref ~ temp + aou + nitrate + silicate 1275 2.3579158 5817.62738 14.1319042 1982-1999 –> 2000-2009 4.1053803 14722.64297
Atlantic (27.25,27.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1275 2.2544179 5705.16735 16.2664244 1982-1999 –> 2000-2009 3.9929150 14589.04135
Atlantic (27.25,27.5] 3 cstar_tref ~ sal + aou + nitrate 1250 2.4058874 5752.14338 17.2939253 2000-2009 –> 2010-2019 4.7270176 11527.67493
Atlantic (27.25,27.5] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 1250 2.4016139 5749.69876 17.8926773 2000-2009 –> 2010-2019 4.7210846 11525.40652
Atlantic (27.25,27.5] 4 cstar_tref ~ sal + aou + nitrate + silicate 1250 2.4058244 5754.07796 17.3496007 2000-2009 –> 2010-2019 4.6689870 11467.11735
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1250 2.3990587 5749.03751 17.8411943 2000-2009 –> 2010-2019 4.6432046 11442.55950
Atlantic (27.25,27.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1250 2.7359536 6077.54640 17.6328991 2000-2009 –> 2010-2019 4.9903714 11782.71374
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1631 2.5518331 7698.45813 15.0110445 1982-1999 –> 2000-2009 5.0736548 20796.21987
Atlantic (27.5,27.75] 3 cstar_tref ~ sal + aou + phosphate 1631 2.4792566 7600.33890 16.8090760 1982-1999 –> 2000-2009 4.8445891 20336.50151
Atlantic (27.5,27.75] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 1631 2.4727938 7593.82467 16.4805262 1982-1999 –> 2000-2009 4.7866808 20209.23993
Atlantic (27.5,27.75] 4 cstar_tref ~ sal + aou + silicate + phosphate 1631 2.4789089 7601.88145 16.6241717 1982-1999 –> 2000-2009 4.8065299 20250.33066
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1631 2.4721289 7594.94743 16.2149539 1982-1999 –> 2000-2009 4.7571593 20142.31168
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1563 2.7779397 7643.46587 16.3288604 2000-2009 –> 2010-2019 5.3297728 15341.92400
Atlantic (27.5,27.75] 3 cstar_tref ~ sal + aou + phosphate 1563 2.7066907 7558.24363 18.5776132 2000-2009 –> 2010-2019 5.1859472 15158.58253
Atlantic (27.5,27.75] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 1563 2.6535477 7498.25739 17.8405767 2000-2009 –> 2010-2019 5.1263415 15092.08207
Atlantic (27.5,27.75] 4 cstar_tref ~ sal + aou + silicate + phosphate 1563 2.6907008 7541.72194 19.9005117 2000-2009 –> 2010-2019 5.1696097 15143.60340
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1563 2.6430173 7487.82742 18.9469447 2000-2009 –> 2010-2019 5.1151462 15082.77485
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 548 1.6670153 2129.25079 12.7128400 1982-1999 –> 2000-2009 3.5147695 6038.30532
Atlantic (27.75,27.85] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 548 1.6298895 2102.56603 12.8120279 1982-1999 –> 2000-2009 3.3437122 5865.45203
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 548 1.6268177 2102.49849 13.5249128 1982-1999 –> 2000-2009 3.2966110 5817.51758
Atlantic (27.75,27.85] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 548 1.6451874 2112.80494 11.9279795 1982-1999 –> 2000-2009 3.3622383 5879.29657
Atlantic (27.75,27.85] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 548 1.6176283 2096.28995 12.8816205 1982-1999 –> 2000-2009 3.3011216 5826.96497
Atlantic (27.75,27.85] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 538 1.7555650 2144.34070 11.4664328 2000-2009 –> 2010-2019 3.4238619 4272.43370
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 538 1.7534279 2145.03006 11.3302496 2000-2009 –> 2010-2019 3.4204433 4274.28085
Atlantic (27.75,27.85] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 538 1.6741054 2093.21799 11.9174318 2000-2009 –> 2010-2019 3.3039949 4195.78402
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 538 1.6419754 2074.36632 14.0267093 2000-2009 –> 2010-2019 3.2687932 4176.86481
Atlantic (27.75,27.85] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 538 1.7404094 2137.01137 11.5109191 2000-2009 –> 2010-2019 3.3580377 4233.30132
Atlantic (27.85,27.95] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 508 2.5951281 2422.53560 18.4224106 1982-1999 –> 2000-2009 5.5537328 7251.57036
Atlantic (27.85,27.95] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 508 2.4911894 2383.00599 18.4030570 1982-1999 –> 2000-2009 5.4447613 7210.76506
Atlantic (27.85,27.95] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 508 2.6215708 2434.83560 19.8781217 1982-1999 –> 2000-2009 5.5563207 7250.29457
Atlantic (27.85,27.95] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 508 2.4714104 2372.90717 19.5024441 1982-1999 –> 2000-2009 5.5649150 7287.72702
Atlantic (27.85,27.95] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 508 2.4038991 2346.76709 18.8270916 1982-1999 –> 2000-2009 5.3572685 7174.39419
Atlantic (27.85,27.95] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 601 3.1147584 3083.21834 18.3655751 2000-2009 –> 2010-2019 5.7098865 5505.75395
Atlantic (27.85,27.95] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 601 3.0957711 3077.86860 18.3064684 2000-2009 –> 2010-2019 5.5869605 5460.87458
Atlantic (27.85,27.95] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 601 3.1756178 3108.47775 18.4104546 2000-2009 –> 2010-2019 5.7971886 5543.31335
Atlantic (27.85,27.95] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 601 3.2579695 3137.25131 19.0250062 2000-2009 –> 2010-2019 5.7293799 5510.15848
Atlantic (27.85,27.95] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 601 3.0792576 3071.43973 18.1625335 2000-2009 –> 2010-2019 5.4831568 5418.20682
Atlantic (27.95,28.05] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 561 2.6799343 2710.10799 13.7947032 1982-1999 –> 2000-2009 4.9571200 7174.47938
Atlantic (27.95,28.05] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 561 2.6709758 2708.35105 14.2072892 1982-1999 –> 2000-2009 4.9183744 7148.57293
Atlantic (27.95,28.05] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 561 3.1521128 2894.18688 17.2569229 1982-1999 –> 2000-2009 6.8999005 8350.03393
Atlantic (27.95,28.05] 3 cstar_tref ~ temp + nitrate + phosphate_star 561 3.1716693 2897.12655 16.3599089 1982-1999 –> 2000-2009 6.8603644 8317.41022
Atlantic (27.95,28.05] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 561 2.9074253 2801.52362 20.9187753 1982-1999 –> 2000-2009 6.1517445 7968.86994
Atlantic (27.95,28.05] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 548 2.0552036 2356.68755 14.6860391 2000-2009 –> 2010-2019 4.7351380 5066.79554
Atlantic (27.95,28.05] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 548 2.0410787 2351.12900 15.2708109 2000-2009 –> 2010-2019 4.7120545 5059.48005
Atlantic (27.95,28.05] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 548 2.8516516 2717.65320 13.8711633 2000-2009 –> 2010-2019 6.0037644 5611.84008
Atlantic (27.95,28.05] 3 cstar_tref ~ temp + nitrate + phosphate_star 548 2.7023344 2654.70775 16.9642960 2000-2009 –> 2010-2019 5.8740037 5551.83430
Atlantic (27.95,28.05] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 548 2.5076716 2574.76935 14.2753672 2000-2009 –> 2010-2019 5.4150969 5376.29297
Atlantic (28.05,28.1] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 341 0.2087418 -88.74406 1.2493465 1982-1999 –> 2000-2009 0.3245420 -956.66317
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 341 0.2012880 -111.54243 1.3589563 1982-1999 –> 2000-2009 0.3163031 -985.58446
Atlantic (28.05,28.1] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 341 0.1950798 -132.90819 1.1252386 1982-1999 –> 2000-2009 0.3108164 -999.48283
Atlantic (28.05,28.1] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 341 0.2081839 -88.56935 1.2646455 1982-1999 –> 2000-2009 0.3237401 -957.00675
Atlantic (28.05,28.1] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 341 0.2024689 -109.55325 1.2406897 1982-1999 –> 2000-2009 0.3207104 -952.56141
Atlantic (28.05,28.1] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 324 0.2306831 -18.95608 0.9145692 2000-2009 –> 2010-2019 0.4339468 -125.83719
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 324 0.2258683 -30.62419 1.3447007 2000-2009 –> 2010-2019 0.4271564 -142.16662
Atlantic (28.05,28.1] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 324 0.2197135 -48.52695 1.3007987 2000-2009 –> 2010-2019 0.4147934 -181.43514
Atlantic (28.05,28.1] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 324 0.2342013 -7.14807 1.6360117 2000-2009 –> 2010-2019 0.4423852 -95.71742
Atlantic (28.05,28.1] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 324 0.2280760 -26.32138 1.4090507 2000-2009 –> 2010-2019 0.4305449 -135.87463
Atlantic (28.1,28.15] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 415 0.1467305 -403.18174 1.2368741 1982-1999 –> 2000-2009 0.3094692 -1005.22829
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 415 0.1318290 -490.06797 1.0521829 1982-1999 –> 2000-2009 0.2667151 -1380.69790
Atlantic (28.1,28.15] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 415 0.1483940 -393.82501 1.3207171 1982-1999 –> 2000-2009 0.2766644 -1364.30368
Atlantic (28.1,28.15] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 415 0.1445007 -413.89189 1.3476647 1982-1999 –> 2000-2009 0.2643564 -1487.40494
Atlantic (28.1,28.15] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 415 0.1479727 -394.18496 1.1789236 1982-1999 –> 2000-2009 0.2922290 -1180.84885
Atlantic (28.1,28.15] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 414 0.1474216 -398.29069 0.7159833 2000-2009 –> 2010-2019 0.2941521 -801.47243
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 414 0.1400241 -438.91811 0.8769197 2000-2009 –> 2010-2019 0.2718531 -928.98608
Atlantic (28.1,28.15] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 414 0.1716549 -272.27781 1.1772440 2000-2009 –> 2010-2019 0.3200489 -666.10282
Atlantic (28.1,28.15] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 414 0.1664485 -295.78043 1.2753245 2000-2009 –> 2010-2019 0.3109492 -709.67231
Atlantic (28.1,28.15] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 414 0.1599298 -328.85974 0.8852545 2000-2009 –> 2010-2019 0.3079025 -723.04469
Atlantic (28.15,28.2] 3 cstar_tref ~ sal + aou + nitrate 894 0.1597045 -732.89925 0.9414553 1982-1999 –> 2000-2009 0.3658843 -1277.69482
Atlantic (28.15,28.2] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 894 0.1563762 -768.55513 0.8993689 1982-1999 –> 2000-2009 0.3609002 -1339.29809
Atlantic (28.15,28.2] 4 cstar_tref ~ sal + aou + nitrate + silicate 894 0.1592362 -736.15002 0.9284112 1982-1999 –> 2000-2009 0.3516730 -1518.03234
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 894 0.1348224 -1031.72684 0.6694043 1982-1999 –> 2000-2009 0.3091863 -2153.43935
Atlantic (28.15,28.2] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 894 0.1658694 -661.17736 0.6493076 1982-1999 –> 2000-2009 0.3881925 -940.69574
Atlantic (28.15,28.2] 3 cstar_tref ~ sal + aou + nitrate 897 0.1990963 -339.88087 2.1482948 2000-2009 –> 2010-2019 0.3588007 -1072.78012
Atlantic (28.15,28.2] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 897 0.1875675 -444.89188 1.6102019 2000-2009 –> 2010-2019 0.3439437 -1213.44702
Atlantic (28.15,28.2] 4 cstar_tref ~ sal + aou + nitrate + silicate 897 0.1789072 -529.69739 1.5098035 2000-2009 –> 2010-2019 0.3381434 -1265.84741
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 897 0.1758130 -558.99594 1.6224889 2000-2009 –> 2010-2019 0.3106354 -1590.72278
Atlantic (28.15,28.2] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 897 0.2340711 -45.54584 2.2218448 2000-2009 –> 2010-2019 0.3999405 -706.72321
Atlantic (28.2, Inf] 3 cstar_tref ~ sal + aou + nitrate 747 0.1843809 -396.08874 1.7379447 1982-1999 –> 2000-2009 0.4687906 14.35409
Atlantic (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 747 0.1502977 -699.44116 0.6477440 1982-1999 –> 2000-2009 0.3624924 -1012.83322
Atlantic (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + silicate 747 0.1686947 -526.92446 0.7415469 1982-1999 –> 2000-2009 0.4023233 -601.86291
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 747 0.1371243 -834.48596 0.5226156 1982-1999 –> 2000-2009 0.3450476 -1196.26882
Atlantic (28.2, Inf] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 747 0.1748600 -471.29725 0.7311327 1982-1999 –> 2000-2009 0.5290205 486.65569
Atlantic (28.2, Inf] 3 cstar_tref ~ sal + aou + nitrate 760 0.2767034 213.87669 3.6305716 2000-2009 –> 2010-2019 0.4610843 -182.21205
Atlantic (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 760 0.2084016 -215.01143 2.1396933 2000-2009 –> 2010-2019 0.3586993 -914.45259
Atlantic (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + silicate 760 0.2329760 -45.57945 2.3320636 2000-2009 –> 2010-2019 0.4016708 -572.50391
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 760 0.2013553 -265.29380 2.2977931 2000-2009 –> 2010-2019 0.3384795 -1099.77977
Atlantic (28.2, Inf] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 760 0.3463012 558.90795 5.8768655 2000-2009 –> 2010-2019 0.5211612 87.61071
Indo-Pacific (-Inf,26] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 4068 6.2543379 26472.01165 30.1837950 1982-1999 –> 2000-2009 12.7847142 73878.49896
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 4068 6.2120635 26418.83199 29.7389004 1982-1999 –> 2000-2009 12.6917293 73715.20265
Indo-Pacific (-Inf,26] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 4068 7.1874512 27605.41459 31.7177713 1982-1999 –> 2000-2009 14.7162042 77059.95522
Indo-Pacific (-Inf,26] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 4068 7.2871451 27717.48996 29.3839953 1982-1999 –> 2000-2009 14.8720894 77278.97281
Indo-Pacific (-Inf,26] 4 cstar_tref ~ temp + silicate + phosphate + phosphate_star 4068 7.3926659 27832.45762 28.6983368 1982-1999 –> 2000-2009 15.1076209 77636.36810
Indo-Pacific (-Inf,26] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 4035 6.0307092 25963.53139 24.7448809 2000-2009 –> 2010-2019 12.2850471 52435.54304
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 4035 5.9655736 25877.89601 24.0313097 2000-2009 –> 2010-2019 12.1776371 52296.72799
Indo-Pacific (-Inf,26] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 4035 7.0479983 27223.47524 25.7430663 2000-2009 –> 2010-2019 14.2354495 54828.88983
Indo-Pacific (-Inf,26] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 4035 7.1072570 27291.04313 25.2743149 2000-2009 –> 2010-2019 14.3944022 55008.53309
Indo-Pacific (-Inf,26] 4 cstar_tref ~ temp + silicate + phosphate + phosphate_star 4035 7.2226523 27419.01749 24.3078224 2000-2009 –> 2010-2019 14.6153182 55251.47511
Indo-Pacific (26,26.5] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 4157 4.7965533 24844.55567 41.2898774 1982-1999 –> 2000-2009 9.4286748 69702.63325
Indo-Pacific (26,26.5] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 4157 4.4613119 24244.16704 44.4049296 1982-1999 –> 2000-2009 8.7037624 67769.26409
Indo-Pacific (26,26.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 4157 4.6544602 24596.53985 49.2819809 1982-1999 –> 2000-2009 9.0525509 68668.99651
Indo-Pacific (26,26.5] 4 cstar_tref ~ temp + aou + silicate + phosphate 4157 4.8143407 24875.32995 39.4752610 1982-1999 –> 2000-2009 9.4241692 69660.11655
Indo-Pacific (26,26.5] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 4157 4.7928690 24840.16704 40.8786832 1982-1999 –> 2000-2009 9.4020789 69624.91492
Indo-Pacific (26,26.5] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 4147 4.9103659 24979.32040 41.8144330 2000-2009 –> 2010-2019 9.7069193 49823.87607
Indo-Pacific (26,26.5] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 4147 4.5645777 24375.67208 44.6175767 2000-2009 –> 2010-2019 9.0258896 48619.83912
Indo-Pacific (26,26.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 4147 4.7225019 24657.77357 47.4590040 2000-2009 –> 2010-2019 9.3769621 49254.31342
Indo-Pacific (26,26.5] 4 cstar_tref ~ temp + aou + silicate + phosphate 4147 4.9019527 24965.09752 40.5969285 2000-2009 –> 2010-2019 9.7162934 49840.42747
Indo-Pacific (26,26.5] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 4147 4.8970961 24958.87626 41.1659746 2000-2009 –> 2010-2019 9.6899651 49799.04330
Indo-Pacific (26.5,26.75] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 3304 4.0150866 18573.85511 21.5932210 1982-1999 –> 2000-2009 7.9621423 51686.72174
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 3304 3.9239979 18424.21536 19.7815762 1982-1999 –> 2000-2009 7.6795840 50949.53560
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 3304 3.9870041 18529.47473 33.0612751 1982-1999 –> 2000-2009 7.9776989 51774.70329
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 3304 3.9255794 18426.87811 19.9955038 1982-1999 –> 2000-2009 7.7903467 51291.95720
Indo-Pacific (26.5,26.75] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 3304 3.9874092 18528.14615 32.8662317 1982-1999 –> 2000-2009 8.0523685 51989.97992
Indo-Pacific (26.5,26.75] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 3289 4.3279970 18983.23737 25.1933400 2000-2009 –> 2010-2019 8.3430836 37557.09247
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 3289 4.1937765 18778.00926 25.3493782 2000-2009 –> 2010-2019 8.1177744 37202.22462
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 3289 4.2563980 18875.50584 32.7849538 2000-2009 –> 2010-2019 8.2434020 37404.98057
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 3289 4.1875904 18768.29906 23.3236715 2000-2009 –> 2010-2019 8.1131698 37195.17717
Indo-Pacific (26.5,26.75] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 3289 4.2579089 18875.84056 33.6994872 2000-2009 –> 2010-2019 8.2453181 37403.98671
Indo-Pacific (26.75,27] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 3405 4.5268950 19958.31846 18.6046621 1982-1999 –> 2000-2009 8.5054978 55971.03692
Indo-Pacific (26.75,27] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 3405 4.3576440 19698.82522 21.0579523 1982-1999 –> 2000-2009 8.0633970 54798.05734
Indo-Pacific (26.75,27] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 3405 4.3536840 19694.63386 20.7258047 1982-1999 –> 2000-2009 8.0547329 54779.53353
Indo-Pacific (26.75,27] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 3405 4.1357100 19342.84928 23.1363869 1982-1999 –> 2000-2009 7.7196231 54012.22259
Indo-Pacific (26.75,27] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 3405 4.1137006 19308.51105 21.9290321 1982-1999 –> 2000-2009 7.6794083 53914.40243
Indo-Pacific (26.75,27] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 3612 4.8788400 21711.78363 18.8780206 2000-2009 –> 2010-2019 9.4057350 41670.10209
Indo-Pacific (26.75,27] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 3612 4.7014331 21444.20582 20.4446770 2000-2009 –> 2010-2019 9.0590770 41143.03105
Indo-Pacific (26.75,27] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 3612 4.6989258 21442.35230 20.1493075 2000-2009 –> 2010-2019 9.0526098 41136.98616
Indo-Pacific (26.75,27] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 3612 4.4935599 21117.52118 22.8082610 2000-2009 –> 2010-2019 8.6292699 40460.37045
Indo-Pacific (26.75,27] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 3612 4.4764669 21091.98941 21.6916624 2000-2009 –> 2010-2019 8.5901674 40400.50047
Indo-Pacific (27,27.25] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 3305 4.1206902 18751.08040 42.1783348 1982-1999 –> 2000-2009 7.5215488 49247.05534
Indo-Pacific (27,27.25] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 3305 4.1676553 18825.99107 49.3300852 1982-1999 –> 2000-2009 7.5534851 49270.88292
Indo-Pacific (27,27.25] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 3305 4.0887569 18701.65670 43.8286562 1982-1999 –> 2000-2009 7.4252208 48979.14332
Indo-Pacific (27,27.25] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 3305 3.6228288 17899.94038 39.3747846 1982-1999 –> 2000-2009 6.6455441 47036.37513
Indo-Pacific (27,27.25] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 3305 3.3193717 17323.69994 28.2067888 1982-1999 –> 2000-2009 6.1167420 45568.52973
Indo-Pacific (27,27.25] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 3306 4.9389310 19954.37014 34.0026065 2000-2009 –> 2010-2019 9.0596212 38705.45054
Indo-Pacific (27,27.25] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 3306 4.9345637 19950.52090 34.6865331 2000-2009 –> 2010-2019 9.0233207 38652.17760
Indo-Pacific (27,27.25] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 3306 5.0949027 20161.94824 36.3638541 2000-2009 –> 2010-2019 9.3250022 39088.24262
Indo-Pacific (27,27.25] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 3306 4.4229353 19224.76673 36.2569818 2000-2009 –> 2010-2019 8.0457641 37124.70710
Indo-Pacific (27,27.25] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 3306 3.8572215 18321.87184 24.7883316 2000-2009 –> 2010-2019 7.1765932 35645.57178
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 2710 3.0383271 13727.93123 30.4385119 1982-1999 –> 2000-2009 5.9181494 37106.80865
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2710 3.1483319 13920.69726 29.0568899 1982-1999 –> 2000-2009 6.0624039 37411.11123
Indo-Pacific (27.25,27.5] 4 cstar_tref ~ sal + silicate + phosphate + phosphate_star 2710 3.1613680 13941.09308 29.6302147 1982-1999 –> 2000-2009 6.0877443 37469.25719
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2710 3.0326431 13717.78218 29.2804379 1982-1999 –> 2000-2009 5.9617813 37256.84591
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 2710 2.8791028 13436.18144 30.1616835 1982-1999 –> 2000-2009 5.6918465 36592.71713
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 2689 3.2832019 14038.52069 28.7570648 2000-2009 –> 2010-2019 6.3215290 27766.45192
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2689 3.4797226 14351.16236 27.2046175 2000-2009 –> 2010-2019 6.6280545 28271.85963
Indo-Pacific (27.25,27.5] 4 cstar_tref ~ sal + silicate + phosphate + phosphate_star 2689 3.4849721 14357.26944 27.7786076 2000-2009 –> 2010-2019 6.6463400 28298.36252
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2689 3.3276671 14110.86746 34.8561136 2000-2009 –> 2010-2019 6.3603102 27828.64964
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 2689 3.1715327 13852.41947 28.4926927 2000-2009 –> 2010-2019 6.0506354 27288.60092
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 2455 1.8024135 9873.59993 12.8005600 1982-1999 –> 2000-2009 3.4281450 27043.09658
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2455 1.8277780 9942.21414 9.5689598 1982-1999 –> 2000-2009 3.4782703 27247.84389
Indo-Pacific (27.5,27.75] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 2455 1.8474303 9992.72492 12.0296908 1982-1999 –> 2000-2009 3.4981424 27297.55350
Indo-Pacific (27.5,27.75] 4 cstar_tref ~ sal + silicate + phosphate + phosphate_star 2455 1.8699257 10052.15079 13.7510023 1982-1999 –> 2000-2009 3.5212956 27360.56755
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2455 1.9218175 10188.55062 12.3152024 1982-1999 –> 2000-2009 3.6188599 27744.66679
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 2398 2.1126340 10406.32796 14.0025218 2000-2009 –> 2010-2019 3.9150476 20279.92789
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2398 2.1203978 10423.92053 12.5218410 2000-2009 –> 2010-2019 3.9481757 20366.13467
Indo-Pacific (27.5,27.75] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 2398 2.1843959 10564.53256 13.2882456 2000-2009 –> 2010-2019 4.0318262 20557.25748
Indo-Pacific (27.5,27.75] 4 cstar_tref ~ sal + silicate + phosphate + phosphate_star 2398 2.2298224 10663.24681 15.1377772 2000-2009 –> 2010-2019 4.0997481 20715.39760
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2398 2.2481778 10704.56479 12.9184376 2000-2009 –> 2010-2019 4.1699953 20893.11542
Indo-Pacific (27.75,27.85] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 898 1.4765903 3260.37869 7.4235980 1982-1999 –> 2000-2009 2.8256624 9073.54699
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 898 1.4676250 3251.44088 7.1481337 1982-1999 –> 2000-2009 2.7959102 9014.18574
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 898 1.4796904 3266.14551 7.4705289 1982-1999 –> 2000-2009 2.8244257 9070.44348
Indo-Pacific (27.75,27.85] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 898 1.4030748 3168.65795 6.9627066 1982-1999 –> 2000-2009 2.7476322 8970.50958
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 898 1.3121619 3050.34386 6.8238078 1982-1999 –> 2000-2009 2.5615515 8606.36395
Indo-Pacific (27.75,27.85] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 1009 1.5363911 3742.02022 12.8415686 2000-2009 –> 2010-2019 3.0129813 7002.39890
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1009 1.5134906 3713.71471 12.2632116 2000-2009 –> 2010-2019 2.9811156 6965.15559
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1009 1.5364475 3744.09428 13.1533169 2000-2009 –> 2010-2019 3.0161379 7010.23979
Indo-Pacific (27.75,27.85] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 1009 1.5821543 3801.25096 10.9603455 2000-2009 –> 2010-2019 2.9852291 6969.90891
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1009 1.4557545 3635.22599 9.7195649 2000-2009 –> 2010-2019 2.7679164 6685.56985
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1359 1.4355555 4853.37294 9.5180780 1982-1999 –> 2000-2009 2.6825299 13167.32522
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1359 1.3919879 4769.60680 9.4266661 1982-1999 –> 2000-2009 2.6686657 13202.72433
Indo-Pacific (27.85,27.95] 4 cstar_tref ~ temp + aou + nitrate + silicate 1359 1.3746171 4733.47525 7.9782189 1982-1999 –> 2000-2009 2.6538372 13174.66266
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1359 1.3646222 4715.64023 8.0208967 1982-1999 –> 2000-2009 2.6053857 13004.31677
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1359 1.4028290 4790.69321 8.5722708 1982-1999 –> 2000-2009 2.6679488 13177.77459
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1425 1.5297374 5269.49859 14.5324837 2000-2009 –> 2010-2019 2.9652929 10122.87154
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1425 1.5031006 5219.43535 12.9206170 2000-2009 –> 2010-2019 2.8950884 9989.04215
Indo-Pacific (27.85,27.95] 4 cstar_tref ~ temp + aou + nitrate + silicate 1425 1.5189199 5247.27332 10.8619868 2000-2009 –> 2010-2019 2.8935370 9980.74857
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1425 1.5049755 5222.98820 11.4914425 2000-2009 –> 2010-2019 2.8695977 9938.62843
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1425 1.5431297 5294.34088 12.2442158 2000-2009 –> 2010-2019 2.9459587 10085.03409
Indo-Pacific (27.95,28.05] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 872 0.6553103 1749.53355 3.2233514 1982-1999 –> 2000-2009 1.2518201 4899.65926
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 872 0.5827310 1546.81710 3.0634322 1982-1999 –> 2000-2009 1.1248246 4366.24816
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 872 0.6592146 1761.89315 5.4562679 1982-1999 –> 2000-2009 1.2518727 4891.48920
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 872 0.5781727 1533.12154 2.9283577 1982-1999 –> 2000-2009 1.1106998 4290.62729
Indo-Pacific (27.95,28.05] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 872 0.6593438 1760.23498 5.6440662 1982-1999 –> 2000-2009 1.2522803 4889.46427
Indo-Pacific (27.95,28.05] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 938 0.6740407 1933.91287 3.3299230 2000-2009 –> 2010-2019 1.3293511 3683.44642
Indo-Pacific (27.95,28.05] 4 cstar_tref ~ sal + aou + silicate + phosphate 938 0.6454597 1852.63005 3.1698214 2000-2009 –> 2010-2019 1.2863574 3563.37860
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 938 0.5996316 1716.46766 2.7687038 2000-2009 –> 2010-2019 1.1823626 3263.28476
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 938 0.6719426 1930.06410 4.0306633 2000-2009 –> 2010-2019 1.3311571 3691.95725
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 938 0.5884165 1681.04778 2.8671423 2000-2009 –> 2010-2019 1.1665892 3214.16932
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 784 0.4821790 1095.13391 2.8225514 1982-1999 –> 2000-2009 0.9516485 3122.67928
Indo-Pacific (28.05,28.1] 4 cstar_tref ~ sal + aou + silicate + phosphate 784 0.4618292 1025.52122 2.4932032 1982-1999 –> 2000-2009 0.9187453 2968.72641
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 784 0.4565536 1009.50668 2.3886492 1982-1999 –> 2000-2009 0.9126696 2949.38666
Indo-Pacific (28.05,28.1] 4 cstar_tref ~ sal + nitrate + silicate + phosphate_star 784 0.4837745 1098.31374 2.5164479 1982-1999 –> 2000-2009 0.9540230 3128.89631
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 784 0.4535666 999.21403 2.3682555 1982-1999 –> 2000-2009 0.9084935 2931.16451
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 796 0.4726976 1080.06546 2.5144282 2000-2009 –> 2010-2019 0.9548766 2175.19937
Indo-Pacific (28.05,28.1] 4 cstar_tref ~ sal + aou + silicate + phosphate 796 0.4731804 1079.69075 2.7476265 2000-2009 –> 2010-2019 0.9350096 2105.21197
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 796 0.4695489 1069.42535 2.6500916 2000-2009 –> 2010-2019 0.9261025 2078.93203
Indo-Pacific (28.05,28.1] 4 cstar_tref ~ sal + nitrate + silicate + phosphate_star 796 0.4874817 1127.09410 2.7540394 2000-2009 –> 2010-2019 0.9712562 2225.40783
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 796 0.4698359 1070.39830 2.6322635 2000-2009 –> 2010-2019 0.9234025 2069.61233
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 7155 0.7920106 16982.19786 3.7498472 1982-1999 –> 2000-2009 1.5829272 47733.11891
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 7155 0.7748402 16668.55197 4.5631453 1982-1999 –> 2000-2009 1.5367473 46449.69794
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 7155 0.6768276 14733.26451 2.6889791 1982-1999 –> 2000-2009 1.3466469 41171.36135
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 7155 0.7739672 16652.41905 4.6223889 1982-1999 –> 2000-2009 1.5411749 46613.49283
Indo-Pacific (28.1, Inf] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 7155 0.7899356 16942.65802 3.2534357 1982-1999 –> 2000-2009 1.5800083 47663.87191
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 7256 0.7759252 16923.95341 4.2097623 2000-2009 –> 2010-2019 1.5679358 33906.15127
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 7256 0.7781670 16965.82233 4.9471578 2000-2009 –> 2010-2019 1.5530073 33634.37430
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 7256 0.6772317 14949.70221 3.0939910 2000-2009 –> 2010-2019 1.3540593 29682.96672
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 7256 0.7716749 16844.24241 5.0329559 2000-2009 –> 2010-2019 1.5456421 33496.66146
Indo-Pacific (28.1, Inf] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 7256 0.7927594 17233.43409 3.1123519 2000-2009 –> 2010-2019 1.5826950 34176.09211
SO (26,26.5] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 78 2.8562219 399.07637 11.9154298 1982-1999 –> 2000-2009 6.4038503 1105.86258
SO (26,26.5] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 78 1.0591712 242.32233 3.8762980 1982-1999 –> 2000-2009 2.3805664 692.31003
SO (26,26.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 78 1.0226157 238.84316 2.7948800 1982-1999 –> 2000-2009 2.2403945 669.76267
SO (26,26.5] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 78 2.8976746 399.32414 11.4913466 1982-1999 –> 2000-2009 6.3856765 1099.73719
SO (26,26.5] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 78 2.7631466 393.90814 11.4529696 1982-1999 –> 2000-2009 6.2013853 1092.61379
SO (26,26.5] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 72 3.1071690 381.58168 11.1922962 2000-2009 –> 2010-2019 6.0170175 783.55985
SO (26,26.5] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 72 3.1319558 382.72585 9.3286768 2000-2009 –> 2010-2019 5.9881777 781.80222
SO (26,26.5] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 72 1.0484115 223.13492 3.0304865 2000-2009 –> 2010-2019 2.1075827 465.45725
SO (26,26.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 72 0.9948997 217.59082 2.4104351 2000-2009 –> 2010-2019 2.0175154 456.43398
SO (26,26.5] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 72 3.2253225 386.95589 9.2938943 2000-2009 –> 2010-2019 5.9884690 780.86404
SO (26.5,26.75] 4 cstar_tref ~ aou + nitrate + silicate + phosphate_star 257 3.0441942 1313.53783 14.4571480 1982-1999 –> 2000-2009 6.6390420 3592.22759
SO (26.5,26.75] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 257 3.0854718 1320.46056 13.2108213 1982-1999 –> 2000-2009 6.5716114 3573.35672
SO (26.5,26.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 257 3.0310386 1313.31175 13.9891916 1982-1999 –> 2000-2009 6.5147995 3567.63456
SO (26.5,26.75] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 257 1.7228381 1020.93653 7.4551373 1982-1999 –> 2000-2009 3.6377498 2770.56894
SO (26.5,26.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 257 1.7071493 1018.23439 7.0687379 1982-1999 –> 2000-2009 3.5862772 2754.02132
SO (26.5,26.75] 4 cstar_tref ~ aou + nitrate + silicate + phosphate_star 188 2.9416066 951.20831 10.8740706 2000-2009 –> 2010-2019 5.9858008 2264.74614
SO (26.5,26.75] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 188 3.0096139 959.80213 10.0071951 2000-2009 –> 2010-2019 6.0950858 2280.26269
SO (26.5,26.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 188 2.9044679 948.43096 10.4652547 2000-2009 –> 2010-2019 5.9355065 2261.74271
SO (26.5,26.75] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 188 1.6334978 730.03297 5.1951588 2000-2009 –> 2010-2019 3.3563360 1750.96950
SO (26.5,26.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 188 1.6334957 732.03247 5.1944445 2000-2009 –> 2010-2019 3.3406449 1750.26686
SO (26.75,27] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 815 2.4345797 3775.20163 12.3077325 1982-1999 –> 2000-2009 5.1293572 11090.27178
SO (26.75,27] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 815 2.4345397 3777.17485 12.3694531 1982-1999 –> 2000-2009 5.1291966 11094.10938
SO (26.75,27] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 815 2.3521796 3719.07793 11.1711976 1982-1999 –> 2000-2009 4.8399717 10791.99033
SO (26.75,27] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 815 2.3096708 3691.35094 9.2239862 1982-1999 –> 2000-2009 4.7448287 10701.46987
SO (26.75,27] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 815 2.8971151 4058.72596 10.3378660 1982-1999 –> 2000-2009 5.9533660 11755.19232
SO (26.75,27] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 796 2.6059668 3795.77370 16.9226168 2000-2009 –> 2010-2019 5.0405465 7570.97533
SO (26.75,27] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 796 2.5987395 3793.35234 16.2150824 2000-2009 –> 2010-2019 5.0332791 7570.52719
SO (26.75,27] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 796 2.4271619 3682.61259 15.0082632 2000-2009 –> 2010-2019 4.7793416 7401.69052
SO (26.75,27] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 796 2.3280611 3618.24709 12.4194086 2000-2009 –> 2010-2019 4.6377318 7309.59803
SO (26.75,27] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 796 3.0819710 4062.85652 11.8274140 2000-2009 –> 2010-2019 5.9790861 8121.58248
SO (27,27.25] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 2024 1.8640463 8276.75332 10.4025131 1982-1999 –> 2000-2009 3.6875499 22176.19834
SO (27,27.25] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2024 1.8185527 8178.73298 9.8310092 1982-1999 –> 2000-2009 3.6359697 22057.18852
SO (27,27.25] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 2024 1.6683604 7827.79693 6.9020159 1982-1999 –> 2000-2009 3.2781094 20869.93254
SO (27,27.25] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 2024 1.6158565 7700.35738 7.0899316 1982-1999 –> 2000-2009 3.2230355 20733.50623
SO (27,27.25] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 2024 1.9135450 8382.84341 10.2754627 1982-1999 –> 2000-2009 3.8253606 22607.47954
SO (27,27.25] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 1775 2.0933918 7671.92067 18.4851689 2000-2009 –> 2010-2019 3.9574380 15948.67399
SO (27,27.25] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1775 2.0297548 7564.33000 18.1743043 2000-2009 –> 2010-2019 3.8483075 15743.06298
SO (27,27.25] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 1775 1.9244383 7373.18304 15.5965902 2000-2009 –> 2010-2019 3.5927987 15200.97998
SO (27,27.25] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1775 1.8576089 7249.71164 15.2959730 2000-2009 –> 2010-2019 3.4734654 14950.06902
SO (27,27.25] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 1775 2.1180905 7713.55998 19.2129986 2000-2009 –> 2010-2019 4.0316356 16096.40339
SO (27.25,27.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2079 1.6845878 8082.43041 8.4323284 1982-1999 –> 2000-2009 3.1232189 21658.68380
SO (27.25,27.5] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 2079 1.6927112 8102.43285 10.3868706 1982-1999 –> 2000-2009 3.1309881 21676.81266
SO (27.25,27.5] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 2079 1.6927467 8100.52005 10.3740389 1982-1999 –> 2000-2009 3.1340475 21688.87880
SO (27.25,27.5] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 2079 1.6849131 8081.23318 8.0351732 1982-1999 –> 2000-2009 3.1263541 21670.33178
SO (27.25,27.5] 4 cstar_tref ~ temp + silicate + phosphate + phosphate_star 2079 1.6931654 8101.54828 10.3381690 1982-1999 –> 2000-2009 3.1361189 21698.62624
SO (27.25,27.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2105 1.7988627 8459.65233 9.5107422 2000-2009 –> 2010-2019 3.4834506 16542.08274
SO (27.25,27.5] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 2105 1.7989215 8459.78984 9.4176821 2000-2009 –> 2010-2019 3.4916327 16562.22269
SO (27.25,27.5] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 2105 1.8006974 8461.94397 9.3693810 2000-2009 –> 2010-2019 3.4934441 16562.46402
SO (27.25,27.5] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 2105 1.8053542 8472.81744 9.7286635 2000-2009 –> 2010-2019 3.4902673 16554.05063
SO (27.25,27.5] 4 cstar_tref ~ temp + silicate + phosphate + phosphate_star 2105 1.7995239 8459.19935 9.3914773 2000-2009 –> 2010-2019 3.4926892 16560.74763
SO (27.5,27.75] 4 cstar_tref ~ aou + nitrate + silicate + phosphate_star 2601 1.4783874 9427.05000 8.3660068 1982-1999 –> 2000-2009 2.7376758 24432.87189
SO (27.5,27.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 2601 1.4606551 9366.27809 8.1666672 1982-1999 –> 2000-2009 2.7198717 24373.58141
SO (27.5,27.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2601 1.5124249 9547.45967 9.0019142 1982-1999 –> 2000-2009 2.7951205 24722.69231
SO (27.5,27.75] 4 cstar_tref ~ temp + aou + nitrate + silicate 2601 1.5110951 9540.88384 8.8320439 1982-1999 –> 2000-2009 2.7694588 24540.02878
SO (27.5,27.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2601 1.4667692 9388.00757 8.0728106 1982-1999 –> 2000-2009 2.7225665 24370.59430
SO (27.5,27.75] 4 cstar_tref ~ aou + nitrate + silicate + phosphate_star 2611 1.6584349 10063.37284 9.9303479 2000-2009 –> 2010-2019 3.1368223 19490.42284
SO (27.5,27.75] 4 cstar_tref ~ sal + aou + nitrate + silicate 2611 1.6352467 9989.84364 9.5769674 2000-2009 –> 2010-2019 3.1218506 19445.72513
SO (27.5,27.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 2611 1.6184017 9937.77162 9.5868091 2000-2009 –> 2010-2019 3.0790567 19304.04971
SO (27.5,27.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2611 1.6787747 10129.02804 9.9982436 2000-2009 –> 2010-2019 3.1911995 19676.48772
SO (27.5,27.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2611 1.6077308 9903.22660 9.5387931 2000-2009 –> 2010-2019 3.0745000 19291.23418
SO (27.75,27.85] 4 cstar_tref ~ aou + nitrate + silicate + phosphate_star 1063 1.4106588 3760.12813 11.1958760 1982-1999 –> 2000-2009 2.6398578 9639.46578
SO (27.75,27.85] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1063 1.3875379 3726.99398 10.5261950 1982-1999 –> 2000-2009 2.6078977 9582.27824
SO (27.75,27.85] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1063 1.4972227 3888.74214 13.2525930 1982-1999 –> 2000-2009 2.7865423 9942.46403
SO (27.75,27.85] 4 cstar_tref ~ temp + aou + nitrate + silicate 1063 1.5277925 3929.71285 13.2544747 1982-1999 –> 2000-2009 2.7871606 9896.58338
SO (27.75,27.85] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1063 1.2966591 3582.97905 8.5237654 1982-1999 –> 2000-2009 2.5222516 9453.70939
SO (27.75,27.85] 4 cstar_tref ~ aou + nitrate + silicate + phosphate_star 1074 1.4514212 3860.10274 8.4602018 2000-2009 –> 2010-2019 2.8620800 7620.23087
SO (27.75,27.85] 4 cstar_tref ~ sal + aou + nitrate + silicate 1074 1.4815502 3904.23514 8.5439266 2000-2009 –> 2010-2019 2.9523445 7753.11477
SO (27.75,27.85] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1074 1.4142801 3806.42107 8.3084410 2000-2009 –> 2010-2019 2.8018180 7533.41505
SO (27.75,27.85] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1074 1.4934056 3923.35503 8.6482567 2000-2009 –> 2010-2019 2.9906283 7812.09717
SO (27.75,27.85] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1074 1.3396681 3690.00224 8.5248280 2000-2009 –> 2010-2019 2.6363272 7272.98129
SO (27.85,27.95] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1051 1.3668670 3653.52844 7.2992266 1982-1999 –> 2000-2009 2.7013196 9995.37405
SO (27.85,27.95] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 1051 1.5567256 3924.92165 8.0265059 1982-1999 –> 2000-2009 3.0015535 10559.27980
SO (27.85,27.95] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1051 1.4608005 3793.23432 9.4318823 1982-1999 –> 2000-2009 2.8757886 10352.25143
SO (27.85,27.95] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 1051 1.4680755 3801.67669 8.4781742 1982-1999 –> 2000-2009 2.9785872 10600.79688
SO (27.85,27.95] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1051 1.4355684 3756.60970 9.3165775 1982-1999 –> 2000-2009 2.9398192 10542.33703
SO (27.85,27.95] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1057 1.4691394 3826.84276 8.4021962 2000-2009 –> 2010-2019 2.8360063 7480.37120
SO (27.85,27.95] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1057 1.5471283 3936.18670 10.5380499 2000-2009 –> 2010-2019 3.0079287 7729.42101
SO (27.85,27.95] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1057 1.5155053 3892.52921 8.8397525 2000-2009 –> 2010-2019 3.0389422 7774.01471
SO (27.85,27.95] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 1057 1.5643012 3957.52258 8.0089666 2000-2009 –> 2010-2019 3.0323767 7759.19927
SO (27.85,27.95] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1057 1.5031488 3875.22234 9.0444242 2000-2009 –> 2010-2019 2.9387172 7631.83203
SO (27.95,28.05] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 1365 1.6950458 5326.34991 9.3382116 1982-1999 –> 2000-2009 3.1778494 14355.58088
SO (27.95,28.05] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1365 1.5697405 5118.68738 7.5465351 1982-1999 –> 2000-2009 2.9075387 13638.04654
SO (27.95,28.05] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 1365 1.7865409 5469.87008 10.1040513 1982-1999 –> 2000-2009 3.2737828 14513.96717
SO (27.95,28.05] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1365 1.7813097 5463.86463 9.6564982 1982-1999 –> 2000-2009 3.2530261 14457.76448
SO (27.95,28.05] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1365 1.8422859 5555.75168 10.8449566 1982-1999 –> 2000-2009 3.5255062 15217.55746
SO (27.95,28.05] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 1407 1.8393145 5719.72485 9.0348462 2000-2009 –> 2010-2019 3.5343604 11046.07476
SO (27.95,28.05] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1407 1.5901971 5312.18930 7.0616860 2000-2009 –> 2010-2019 3.1599376 10430.87667
SO (27.95,28.05] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 1407 1.7690176 5610.06727 9.3101180 2000-2009 –> 2010-2019 3.5555585 11079.93735
SO (27.95,28.05] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1407 1.7137009 5522.66907 8.7326955 2000-2009 –> 2010-2019 3.4950106 10986.53370
SO (27.95,28.05] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1407 1.8215961 5694.48574 8.6865634 2000-2009 –> 2010-2019 3.6382856 11212.04128
SO (28.05,28.1] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 821 1.2027264 2644.99549 6.8514612 1982-1999 –> 2000-2009 2.3036546 7244.69453
SO (28.05,28.1] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 821 1.1301940 2544.86034 7.9400594 1982-1999 –> 2000-2009 2.1671589 6965.31663
SO (28.05,28.1] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 821 1.2219116 2670.98100 10.1026904 1982-1999 –> 2000-2009 2.3077347 7228.84723
SO (28.05,28.1] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 821 1.2215427 2672.48525 10.0765460 1982-1999 –> 2000-2009 2.3047202 7224.96466
SO (28.05,28.1] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 821 1.2408077 2698.17920 8.1202016 1982-1999 –> 2000-2009 2.4203107 7508.62718
SO (28.05,28.1] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 856 1.3785601 2990.84251 7.0256157 2000-2009 –> 2010-2019 2.5812866 5635.83800
SO (28.05,28.1] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 856 1.2592408 2837.85424 8.9348982 2000-2009 –> 2010-2019 2.3894349 5382.71458
SO (28.05,28.1] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 856 1.3172094 2912.90520 11.1410627 2000-2009 –> 2010-2019 2.5391210 5583.88619
SO (28.05,28.1] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 856 1.3161849 2913.57312 11.2565409 2000-2009 –> 2010-2019 2.5377277 5586.05838
SO (28.05,28.1] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 856 1.3250652 2925.08520 8.7792074 2000-2009 –> 2010-2019 2.5658730 5623.26440
SO (28.1, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 4069 1.0834975 12213.94262 5.9964967 1982-1999 –> 2000-2009 2.0726643 32951.56620
SO (28.1, Inf] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 4069 1.0886453 12252.51507 6.0780739 1982-1999 –> 2000-2009 2.0677264 32839.30271
SO (28.1, Inf] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 4069 1.0886459 12250.51991 6.0758449 1982-1999 –> 2000-2009 2.0780008 32988.94303
SO (28.1, Inf] 4 cstar_tref ~ sal + nitrate + silicate + phosphate_star 4069 1.0834982 12211.94793 5.9950112 1982-1999 –> 2000-2009 2.0766196 33006.29616
SO (28.1, Inf] 4 cstar_tref ~ sal + silicate + phosphate + phosphate_star 4069 1.0937618 12288.67298 5.8346568 1982-1999 –> 2000-2009 2.1035317 33327.70584
SO (28.1, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 3951 1.1694332 12463.26688 6.2409335 2000-2009 –> 2010-2019 2.2529307 24677.20949
SO (28.1, Inf] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 3951 1.1776558 12518.63386 6.1193967 2000-2009 –> 2010-2019 2.2663011 24771.14893
SO (28.1, Inf] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 3951 1.1785273 12522.47882 6.1357060 2000-2009 –> 2010-2019 2.2671732 24772.99873
SO (28.1, Inf] 4 cstar_tref ~ sal + nitrate + silicate + phosphate_star 3951 1.1702010 12466.45330 6.2291769 2000-2009 –> 2010-2019 2.2536992 24678.40123
SO (28.1, Inf] 4 cstar_tref ~ sal + silicate + phosphate + phosphate_star 3951 1.1791229 12526.47187 6.0728150 2000-2009 –> 2010-2019 2.2728847 24815.14486
SO (28.1,28.15] 4 cstar_tref ~ sal + aou + nitrate + silicate 317 0.5074183 481.48897 1.8520118 1982-1999 –> 2000-2009 0.9999260 1234.03153
SO (28.1,28.15] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 317 0.4660177 429.52794 1.4620719 1982-1999 –> 2000-2009 0.9347512 1132.51668
SO (28.1,28.15] 4 cstar_tref ~ sal + aou + silicate + phosphate 317 0.5107772 485.67200 1.7536115 1982-1999 –> 2000-2009 1.0062340 1244.43536
SO (28.1,28.15] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 317 0.4547179 413.96550 1.7741564 1982-1999 –> 2000-2009 0.8957964 1053.58875
SO (28.1,28.15] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 317 0.4894693 458.65621 1.8204796 1982-1999 –> 2000-2009 0.9743042 1194.83743
SO (28.1,28.15] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 287 0.5019925 432.88706 2.0864768 2000-2009 –> 2010-2019 0.9680102 862.41500
SO (28.1,28.15] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 287 0.4979308 428.22387 2.3615801 2000-2009 –> 2010-2019 0.9526487 842.18937
SO (28.1,28.15] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 287 0.5563859 489.93843 2.1103649 2000-2009 –> 2010-2019 1.0458552 948.59464
SO (28.1,28.15] 4 cstar_tref ~ sal + nitrate + silicate + phosphate_star 287 0.5751917 509.01891 2.0875759 2000-2009 –> 2010-2019 1.0842619 992.56860
SO (28.1,28.15] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 287 0.5587642 494.38676 2.6541867 2000-2009 –> 2010-2019 1.0907799 1007.88763
SO (28.15,28.2] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 293 0.2980246 136.09844 1.1340674 1982-1999 –> 2000-2009 0.6050153 383.81517
SO (28.15,28.2] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 293 0.2999722 137.91558 1.2075725 1982-1999 –> 2000-2009 0.6120788 399.86231
SO (28.15,28.2] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 293 0.2998770 139.72954 1.2196852 1982-1999 –> 2000-2009 0.5698869 261.39816
SO (28.15,28.2] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 293 0.2988984 135.81426 1.0216917 1982-1999 –> 2000-2009 0.6301386 456.18872
SO (28.15,28.2] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 293 0.2988106 137.64208 1.0327212 1982-1999 –> 2000-2009 0.6292472 457.63138
SO (28.15,28.2] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 275 0.3282835 181.78341 0.7568624 2000-2009 –> 2010-2019 0.6263080 317.88185
SO (28.15,28.2] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 275 0.3418232 202.01235 0.9409124 2000-2009 –> 2010-2019 0.6417954 339.92793
SO (28.15,28.2] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 275 0.3349910 192.90778 0.9285407 2000-2009 –> 2010-2019 0.6348680 332.63732
SO (28.15,28.2] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 275 0.3324954 186.79516 1.2839449 2000-2009 –> 2010-2019 0.6313939 322.60942
SO (28.15,28.2] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 275 0.3293587 183.58191 1.1533303 2000-2009 –> 2010-2019 0.6281694 321.22400
SO (28.2, Inf] 3 cstar_tref ~ sal + aou + nitrate 1351 0.3716747 1169.70459 1.9520150 1982-1999 –> 2000-2009 0.6754315 2304.02223
SO (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 1351 0.3627143 1105.76586 1.8413473 1982-1999 –> 2000-2009 0.6664474 2241.69906
SO (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + silicate 1351 0.3635557 1112.02715 1.9253716 1982-1999 –> 2000-2009 0.6673032 2248.19301
SO (28.2, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1351 0.3626909 1107.59174 1.8488744 1982-1999 –> 2000-2009 0.6664033 2245.18679
SO (28.2, Inf] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1351 0.3902633 1305.56908 2.0860609 1982-1999 –> 2000-2009 0.7153295 2779.09780
SO (28.2, Inf] 3 cstar_tref ~ sal + aou + nitrate 1392 0.4050057 1443.99501 2.4996960 2000-2009 –> 2010-2019 0.7766804 2613.69961
SO (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 1392 0.3912704 1349.94051 2.3406736 2000-2009 –> 2010-2019 0.7539846 2455.70638
SO (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + silicate 1392 0.3890765 1334.28618 2.4375060 2000-2009 –> 2010-2019 0.7526322 2446.31333
SO (28.2, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1392 0.3890463 1336.07053 2.4521701 2000-2009 –> 2010-2019 0.7517372 2443.66227
SO (28.2, Inf] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1392 0.4153073 1517.92209 2.5374216 2000-2009 –> 2010-2019 0.8055705 2823.49117

5.4.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-2009 2.8267909 1201.2279
Atlantic (-Inf,26] 2000-2009 –> 2010-2019 2.7159670 811.7344
Atlantic (26,26.5] 1982-1999 –> 2000-2009 8.0110933 13071.9049
Atlantic (26,26.5] 2000-2009 –> 2010-2019 7.9576782 8991.6530
Atlantic (26.5,26.75] 1982-1999 –> 2000-2009 6.2582126 14677.3654
Atlantic (26.5,26.75] 2000-2009 –> 2010-2019 6.7112183 10564.3236
Atlantic (26.75,27] 1982-1999 –> 2000-2009 3.8214161 17683.5193
Atlantic (26.75,27] 2000-2009 –> 2010-2019 4.5695724 14617.4923
Atlantic (27,27.25] 1982-1999 –> 2000-2009 3.9091510 14629.2709
Atlantic (27,27.25] 2000-2009 –> 2010-2019 4.4320436 11976.0373
Atlantic (27.25,27.5] 1982-1999 –> 2000-2009 4.1961918 14980.0224
Atlantic (27.25,27.5] 2000-2009 –> 2010-2019 4.7501331 11549.0944
Atlantic (27.5,27.75] 1982-1999 –> 2000-2009 4.8537228 20346.9207
Atlantic (27.5,27.75] 2000-2009 –> 2010-2019 5.1853635 15163.7934
Atlantic (27.75,27.85] 1982-1999 –> 2000-2009 3.3636905 5885.5073
Atlantic (27.75,27.85] 2000-2009 –> 2010-2019 3.3550262 4230.5329
Atlantic (27.85,27.95] 1982-1999 –> 2000-2009 5.4953997 7234.9502
Atlantic (27.85,27.95] 2000-2009 –> 2010-2019 5.6613144 5487.6614
Atlantic (27.95,28.05] 1982-1999 –> 2000-2009 5.9575008 7791.8733
Atlantic (27.95,28.05] 2000-2009 –> 2010-2019 5.3480115 5333.2486
Atlantic (28.05,28.1] 1982-1999 –> 2000-2009 0.3192224 -970.2597
Atlantic (28.05,28.1] 2000-2009 –> 2010-2019 0.4297653 -136.2062
Atlantic (28.1,28.15] 1982-1999 –> 2000-2009 0.2818868 -1283.6967
Atlantic (28.1,28.15] 2000-2009 –> 2010-2019 0.3009812 -765.8557
Atlantic (28.15,28.2] 1982-1999 –> 2000-2009 0.3551673 -1445.8321
Atlantic (28.15,28.2] 2000-2009 –> 2010-2019 0.3502928 -1169.9041
Atlantic (28.2, Inf] 1982-1999 –> 2000-2009 0.4215349 -461.9910
Atlantic (28.2, Inf] 2000-2009 –> 2010-2019 0.4162190 -536.2675
Indo-Pacific (-Inf,26] 1982-1999 –> 2000-2009 14.0344716 75913.7995
Indo-Pacific (-Inf,26] 2000-2009 –> 2010-2019 13.5415708 53964.2338
Indo-Pacific (26,26.5] 1982-1999 –> 2000-2009 9.2022473 69085.1851
Indo-Pacific (26,26.5] 2000-2009 –> 2010-2019 9.5032059 49467.4999
Indo-Pacific (26.5,26.75] 1982-1999 –> 2000-2009 7.8924281 51538.5796
Indo-Pacific (26.5,26.75] 2000-2009 –> 2010-2019 8.2125496 37352.6923
Indo-Pacific (26.75,27] 1982-1999 –> 2000-2009 8.0045318 54695.0506
Indo-Pacific (26.75,27] 2000-2009 –> 2010-2019 8.9473718 40962.1980
Indo-Pacific (27,27.25] 1982-1999 –> 2000-2009 7.0525082 48020.3973
Indo-Pacific (27,27.25] 2000-2009 –> 2010-2019 8.5260603 37843.2299
Indo-Pacific (27.25,27.5] 1982-1999 –> 2000-2009 5.9443851 37167.3480
Indo-Pacific (27.25,27.5] 2000-2009 –> 2010-2019 6.4013738 27890.7849
Indo-Pacific (27.5,27.75] 1982-1999 –> 2000-2009 3.5089426 27338.7457
Indo-Pacific (27.5,27.75] 2000-2009 –> 2010-2019 4.0329586 20562.3666
Indo-Pacific (27.75,27.85] 1982-1999 –> 2000-2009 2.7510364 8947.0099
Indo-Pacific (27.75,27.85] 2000-2009 –> 2010-2019 2.9526761 6926.6546
Indo-Pacific (27.85,27.95] 1982-1999 –> 2000-2009 2.6556735 13145.3607
Indo-Pacific (27.85,27.95] 2000-2009 –> 2010-2019 2.9138949 10023.2650
Indo-Pacific (27.95,28.05] 1982-1999 –> 2000-2009 1.1982995 4667.4976
Indo-Pacific (27.95,28.05] 2000-2009 –> 2010-2019 1.2591635 3483.2473
Indo-Pacific (28.05,28.1] 1982-1999 –> 2000-2009 0.9291160 3020.1706
Indo-Pacific (28.05,28.1] 2000-2009 –> 2010-2019 0.9421295 2130.8727
Indo-Pacific (28.1, Inf] 1982-1999 –> 2000-2009 1.5175009 45926.3086
Indo-Pacific (28.1, Inf] 2000-2009 –> 2010-2019 1.5206679 32979.2492
SO (26,26.5] 1982-1999 –> 2000-2009 4.7223746 932.0573
SO (26,26.5] 2000-2009 –> 2010-2019 4.4237525 653.6235
SO (26.5,26.75] 1982-1999 –> 2000-2009 5.3898960 3251.5618
SO (26.5,26.75] 2000-2009 –> 2010-2019 4.9426748 2061.5976
SO (26.75,27] 1982-1999 –> 2000-2009 5.1593440 11086.6067
SO (26.75,27] 2000-2009 –> 2010-2019 5.0939970 7594.8747
SO (27,27.25] 1982-1999 –> 2000-2009 3.5300050 21688.8610
SO (27,27.25] 2000-2009 –> 2010-2019 3.7807290 15587.8379
SO (27.25,27.5] 1982-1999 –> 2000-2009 3.1301455 21678.6667
SO (27.25,27.5] 2000-2009 –> 2010-2019 3.4902968 16556.3135
SO (27.5,27.75] 1982-1999 –> 2000-2009 2.7489387 24487.9537
SO (27.5,27.75] 2000-2009 –> 2010-2019 3.1206858 19441.5839
SO (27.75,27.85] 1982-1999 –> 2000-2009 2.6687420 9702.9002
SO (27.75,27.85] 2000-2009 –> 2010-2019 2.8486396 7598.3678
SO (27.85,27.95] 1982-1999 –> 2000-2009 2.8994136 10410.0078
SO (27.85,27.95] 2000-2009 –> 2010-2019 2.9707942 7674.9676
SO (27.95,28.05] 1982-1999 –> 2000-2009 3.2275406 14436.5833
SO (27.95,28.05] 2000-2009 –> 2010-2019 3.4766305 10951.0928
SO (28.05,28.1] 1982-1999 –> 2000-2009 2.3007158 7234.4900
SO (28.05,28.1] 2000-2009 –> 2010-2019 2.5226886 5562.3523
SO (28.1, Inf] 1982-1999 –> 2000-2009 2.0797086 33022.7628
SO (28.1, Inf] 2000-2009 –> 2010-2019 2.2625978 24742.9806
SO (28.1,28.15] 1982-1999 –> 2000-2009 0.9622024 1171.8820
SO (28.1,28.15] 2000-2009 –> 2010-2019 1.0283112 930.7310
SO (28.15,28.2] 1982-1999 –> 2000-2009 0.6092734 391.7792
SO (28.15,28.2] 2000-2009 –> 2010-2019 0.6325069 326.8561
SO (28.2, Inf] 1982-1999 –> 2000-2009 0.6781830 2363.6398
SO (28.2, Inf] 2000-2009 –> 2010-2019 0.7681210 2556.5746

5.5 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.6 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.7 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.8 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 - 7).

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
e2704ca Donghe-Zhu 2021-03-27
3061a0b Donghe-Zhu 2021-03-27
b883157 Donghe-Zhu 2021-03-27
d19654d Donghe-Zhu 2021-03-26
6c53dbf Donghe-Zhu 2021-03-25
1914a11 Donghe-Zhu 2021-03-24
8be810e Donghe-Zhu 2021-03-23
bf19764 Donghe-Zhu 2021-03-22
3ec9d3d Donghe-Zhu 2021-03-22
134ace1 Donghe-Zhu 2021-03-22
f6d70a4 Donghe-Zhu 2021-03-22
183443b Donghe-Zhu 2021-03-21
2e6976b Donghe-Zhu 2021-03-21
51a42bd Donghe-Zhu 2021-03-16
f745381 Donghe-Zhu 2021-03-16
aecbf75 Donghe-Zhu 2021-03-14
1b2a0c2 Donghe-Zhu 2021-03-14
6733e48 Donghe-Zhu 2021-03-12
ba71e6a Donghe-Zhu 2021-03-12
9dbf5bd Donghe-Zhu 2021-03-11
a49df30 Donghe-Zhu 2021-03-11
b3348a5 Donghe-Zhu 2021-03-11
1c24ff7 Donghe-Zhu 2021-03-10
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
e2704ca Donghe-Zhu 2021-03-27
3061a0b Donghe-Zhu 2021-03-27
b883157 Donghe-Zhu 2021-03-27
d19654d Donghe-Zhu 2021-03-26
6c53dbf Donghe-Zhu 2021-03-25
1914a11 Donghe-Zhu 2021-03-24
8be810e Donghe-Zhu 2021-03-23
bf19764 Donghe-Zhu 2021-03-22
3ec9d3d Donghe-Zhu 2021-03-22
134ace1 Donghe-Zhu 2021-03-22
f6d70a4 Donghe-Zhu 2021-03-22
183443b Donghe-Zhu 2021-03-21
2e6976b Donghe-Zhu 2021-03-21
51a42bd Donghe-Zhu 2021-03-16
f745381 Donghe-Zhu 2021-03-16
aecbf75 Donghe-Zhu 2021-03-14
1b2a0c2 Donghe-Zhu 2021-03-14
6733e48 Donghe-Zhu 2021-03-12
ba71e6a Donghe-Zhu 2021-03-12
9dbf5bd Donghe-Zhu 2021-03-11
a49df30 Donghe-Zhu 2021-03-11
b3348a5 Donghe-Zhu 2021-03-11
1c24ff7 Donghe-Zhu 2021-03-10
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
e2704ca Donghe-Zhu 2021-03-27
3061a0b Donghe-Zhu 2021-03-27
b883157 Donghe-Zhu 2021-03-27
d19654d Donghe-Zhu 2021-03-26
6c53dbf Donghe-Zhu 2021-03-25
1914a11 Donghe-Zhu 2021-03-24
8be810e Donghe-Zhu 2021-03-23
bf19764 Donghe-Zhu 2021-03-22
3ec9d3d Donghe-Zhu 2021-03-22
134ace1 Donghe-Zhu 2021-03-22
f6d70a4 Donghe-Zhu 2021-03-22
183443b Donghe-Zhu 2021-03-21
2e6976b Donghe-Zhu 2021-03-21
51a42bd Donghe-Zhu 2021-03-16
f745381 Donghe-Zhu 2021-03-16
aecbf75 Donghe-Zhu 2021-03-14
1b2a0c2 Donghe-Zhu 2021-03-14
6733e48 Donghe-Zhu 2021-03-12
ba71e6a Donghe-Zhu 2021-03-12
9dbf5bd Donghe-Zhu 2021-03-11
a49df30 Donghe-Zhu 2021-03-11
b3348a5 Donghe-Zhu 2021-03-11
1c24ff7 Donghe-Zhu 2021-03-10
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
e2704ca Donghe-Zhu 2021-03-27
3061a0b Donghe-Zhu 2021-03-27
b883157 Donghe-Zhu 2021-03-27
d19654d Donghe-Zhu 2021-03-26
6c53dbf Donghe-Zhu 2021-03-25
1914a11 Donghe-Zhu 2021-03-24
8be810e Donghe-Zhu 2021-03-23
bf19764 Donghe-Zhu 2021-03-22
3ec9d3d Donghe-Zhu 2021-03-22
134ace1 Donghe-Zhu 2021-03-22
f6d70a4 Donghe-Zhu 2021-03-22
183443b Donghe-Zhu 2021-03-21
2e6976b Donghe-Zhu 2021-03-21
51a42bd Donghe-Zhu 2021-03-16
f745381 Donghe-Zhu 2021-03-16
aecbf75 Donghe-Zhu 2021-03-14
1b2a0c2 Donghe-Zhu 2021-03-14
6733e48 Donghe-Zhu 2021-03-12
ba71e6a Donghe-Zhu 2021-03-12
9dbf5bd Donghe-Zhu 2021-03-11
a49df30 Donghe-Zhu 2021-03-11
b3348a5 Donghe-Zhu 2021-03-11
1c24ff7 Donghe-Zhu 2021-03-10
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 phosphate phosphate_star sal silicate temp
Atlantic - 1982-1999 --> 2000-2009
(-Inf,26] 5 1 4 5 3 3 2
(26,26.5] 5 3 2 4 1 4 4
(26.5,26.75] 5 1 4 4 2 5 2
(26.75,27] 5 1 4 4 2 5 2
(27,27.25] 5 1 4 3 4 3 1
(27.25,27.5] 5 5 0 3 1 3 4
(27.5,27.75] 5 1 4 3 5 3 0
(27.75,27.85] 5 1 4 5 3 3 2
(27.85,27.95] 5 1 4 5 1 4 3
(27.95,28.05] 3 4 1 5 0 3 5
(28.05,28.1] 4 2 3 5 1 5 3
(28.1,28.15] 5 2 3 5 2 3 3
(28.15,28.2] 5 5 0 3 4 3 1
(28.2, Inf] 5 5 0 3 4 3 1
total 67.00 33.00 37.00 57.00 33.00 50.00 33.00
Atlantic - 2000-2009 --> 2010-2019
(-Inf,26] 5 1 4 5 3 3 2
(26,26.5] 5 3 2 4 1 4 4
(26.5,26.75] 5 1 4 4 2 5 2
(26.75,27] 5 1 4 4 2 5 2
(27,27.25] 5 2 2 4 4 5 1
(27.25,27.5] 5 5 0 3 4 3 1
(27.5,27.75] 5 1 4 3 5 3 0
(27.75,27.85] 5 2 3 5 4 3 1
(27.85,27.95] 5 1 4 5 1 4 3
(27.95,28.05] 3 4 1 5 0 3 5
(28.05,28.1] 4 2 3 5 2 4 3
(28.1,28.15] 5 2 3 5 2 3 3
(28.15,28.2] 5 5 0 3 4 3 1
(28.2, Inf] 5 5 0 3 4 3 1
total 67.00 35.00 34.00 58.00 38.00 51.00 29.00
Indo-Pacific - 1982-1999 --> 2000-2009
(-Inf,26] 4 1 4 5 2 4 3
(26,26.5] 5 1 4 4 1 5 3
(26.5,26.75] 4 2 3 5 1 5 3
(26.75,27] 5 0 5 5 2 3 2
(27,27.25] 5 0 5 5 2 3 2
(27.25,27.5] 4 2 3 5 3 5 2
(27.5,27.75] 4 2 2 5 4 5 1
(27.75,27.85] 5 1 4 5 1 4 3
(27.85,27.95] 5 3 2 4 2 5 3
(27.95,28.05] 4 2 3 5 1 5 3
(28.05,28.1] 4 2 3 4 4 5 1
(28.1, Inf] 4 3 2 5 2 5 3
total 53.00 19.00 40.00 57.00 25.00 54.00 29.00
Indo-Pacific - 2000-2009 --> 2010-2019
(-Inf,26] 4 1 4 5 2 4 3
(26,26.5] 5 1 4 4 1 5 3
(26.5,26.75] 4 2 3 5 1 5 3
(26.75,27] 5 0 5 5 2 3 2
(27,27.25] 5 1 4 5 1 4 3
(27.25,27.5] 4 2 3 5 3 5 2
(27.5,27.75] 4 2 2 5 4 5 1
(27.75,27.85] 5 1 4 5 1 4 3
(27.85,27.95] 5 3 2 4 2 5 3
(27.95,28.05] 5 1 4 4 2 5 2
(28.05,28.1] 4 2 3 4 4 5 1
(28.1, Inf] 4 3 2 5 2 5 3
total 54.00 19.00 40.00 56.00 25.00 55.00 29.00
SO - 1982-1999 --> 2000-2009
(26,26.5] 5 2 3 5 1 3 4
(26.5,26.75] 5 5 0 5 2 3 2
(26.75,27] 5 2 2 5 0 3 5
(27,27.25] 5 2 2 5 0 3 5
(27.25,27.5] 3 2 2 5 0 5 5
(27.5,27.75] 5 4 1 4 2 5 2
(27.75,27.85] 5 4 1 4 2 5 2
(27.85,27.95] 5 1 4 5 3 3 2
(27.95,28.05] 5 2 3 5 4 3 1
(28.05,28.1] 5 2 3 5 4 3 1
(28.1, Inf] 3 2 2 5 5 5 0
(28.1,28.15] 5 2 2 3 5 5 0
(28.15,28.2] 5 2 3 5 1 3 4
(28.2, Inf] 5 4 1 3 5 3 0
total 66.00 36.00 29.00 64.00 34.00 52.00 33.00
SO - 2000-2009 --> 2010-2019
(26,26.5] 5 3 2 5 2 4 3
(26.5,26.75] 5 5 0 5 2 3 2
(26.75,27] 5 2 2 5 0 3 5
(27,27.25] 5 2 2 5 0 3 5
(27.25,27.5] 3 2 2 5 0 5 5
(27.5,27.75] 5 4 1 4 3 5 1
(27.75,27.85] 5 4 1 4 3 5 1
(27.85,27.95] 5 2 3 5 2 4 3
(27.95,28.05] 5 3 2 5 4 3 1
(28.05,28.1] 5 2 3 5 4 3 1
(28.1, Inf] 3 2 2 5 5 5 0
(28.1,28.15] 4 3 1 5 4 5 1
(28.15,28.2] 5 2 3 5 1 3 4
(28.2, Inf] 5 4 1 3 5 3 0
total 65.00 40.00 25.00 66.00 35.00 54.00 32.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
e2704ca Donghe-Zhu 2021-03-27
3061a0b Donghe-Zhu 2021-03-27
b883157 Donghe-Zhu 2021-03-27
d19654d Donghe-Zhu 2021-03-26
6c53dbf Donghe-Zhu 2021-03-25
1914a11 Donghe-Zhu 2021-03-24
8be810e Donghe-Zhu 2021-03-23
bf19764 Donghe-Zhu 2021-03-22
3ec9d3d Donghe-Zhu 2021-03-22
134ace1 Donghe-Zhu 2021-03-22
f6d70a4 Donghe-Zhu 2021-03-22
183443b Donghe-Zhu 2021-03-21
2e6976b Donghe-Zhu 2021-03-21
51a42bd Donghe-Zhu 2021-03-16
f745381 Donghe-Zhu 2021-03-16
aecbf75 Donghe-Zhu 2021-03-14
1b2a0c2 Donghe-Zhu 2021-03-14
6733e48 Donghe-Zhu 2021-03-12
ba71e6a Donghe-Zhu 2021-03-12
9dbf5bd Donghe-Zhu 2021-03-11
a49df30 Donghe-Zhu 2021-03-11
b3348a5 Donghe-Zhu 2021-03-11
1c24ff7 Donghe-Zhu 2021-03-10
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
e2704ca Donghe-Zhu 2021-03-27
3061a0b Donghe-Zhu 2021-03-27
b883157 Donghe-Zhu 2021-03-27
d19654d Donghe-Zhu 2021-03-26
6c53dbf Donghe-Zhu 2021-03-25
1914a11 Donghe-Zhu 2021-03-24
8be810e Donghe-Zhu 2021-03-23
bf19764 Donghe-Zhu 2021-03-22
3ec9d3d Donghe-Zhu 2021-03-22
134ace1 Donghe-Zhu 2021-03-22
f6d70a4 Donghe-Zhu 2021-03-22
183443b Donghe-Zhu 2021-03-21
2e6976b Donghe-Zhu 2021-03-21
51a42bd Donghe-Zhu 2021-03-16
f745381 Donghe-Zhu 2021-03-16
aecbf75 Donghe-Zhu 2021-03-14
1b2a0c2 Donghe-Zhu 2021-03-14
6733e48 Donghe-Zhu 2021-03-12
ba71e6a Donghe-Zhu 2021-03-12
9dbf5bd Donghe-Zhu 2021-03-11
a49df30 Donghe-Zhu 2021-03-11
b3348a5 Donghe-Zhu 2021-03-11
1c24ff7 Donghe-Zhu 2021-03-10
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] car_3.0-10       carData_3.0-4    gt_0.2.2         corrr_0.4.3     
 [5] broom_0.7.5      kableExtra_1.3.1 knitr_1.30       olsrr_0.5.3     
 [9] GGally_2.0.0     lubridate_1.7.9  metR_0.9.0       scico_1.2.0     
[13] patchwork_1.1.1  collapse_1.5.0   forcats_0.5.0    stringr_1.4.0   
[17] dplyr_1.0.5      purrr_0.3.4      readr_1.4.0      tidyr_1.1.2     
[21] tibble_3.0.4     ggplot2_3.3.3    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                magrittr_2.0.1           Matrix_1.2-18           
[43] Rcpp_1.0.5               munsell_0.5.0            fansi_0.4.1             
[46] abind_1.4-5              lifecycle_1.0.0          stringi_1.5.3           
[49] whisker_0.4              yaml_2.2.1               plyr_1.8.6              
[52] grid_4.0.3               blob_1.2.1               parallel_4.0.3          
[55] promises_1.1.1           crayon_1.3.4             lattice_0.20-41         
[58] haven_2.3.1              hms_0.5.3                pillar_1.4.7            
[61] reprex_0.3.0             glue_1.4.2               evaluate_0.14           
[64] RcppArmadillo_0.10.1.2.2 data.table_1.13.6        modelr_0.1.8            
[67] vctrs_0.3.6              httpuv_1.5.4             cellranger_1.1.0        
[70] gtable_0.3.0             reshape_0.8.8            assertthat_0.2.1        
[73] xfun_0.20                openxlsx_4.2.3           RcppEigen_0.3.3.9.1     
[76] later_1.1.0.1            viridisLite_0.3.0        ellipsis_0.3.1          
[79] here_1.0.1