Last updated: 2021-03-25

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: 195db95

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

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


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

Unstaged changes:
    Deleted:    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
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, oxygen, aou, silicate, phosphate, phosphate_star
# the following code is a workaround to find all predictor combinations
# using the olsrr package and fit all models for one era, slab, and basin

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

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

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

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

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

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

# convert to tibble
lm_all <- as_tibble(lm_all)

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

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

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

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

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

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

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

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

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

3 Apply predictor threshold

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

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

  • 112

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:

  • 10

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 (10) 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] 3 cstar_tref ~ sal + temp + aou 123 1.1497327 393.3831 3.863522 1982-1999 –> 2000-2009 2.543046 1131.6424
Atlantic (-Inf,26] 4 cstar_tref ~ sal + temp + aou + phosphate_star 123 1.0447717 371.8333 4.212269 1982-1999 –> 2000-2009 2.383036 1095.3232
Atlantic (-Inf,26] 4 cstar_tref ~ sal + temp + aou + silicate 123 1.1041694 385.4359 4.330897 1982-1999 –> 2000-2009 2.461988 1114.9602
Atlantic (-Inf,26] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 123 1.0447623 373.8311 4.206456 1982-1999 –> 2000-2009 2.381967 1098.9914
Atlantic (-Inf,26] 3 cstar_tref ~ sal + temp + oxygen 123 1.1496665 393.3690 3.912193 1982-1999 –> 2000-2009 2.542540 1131.4970
Atlantic (-Inf,26] 4 cstar_tref ~ sal + temp + oxygen + phosphate_star 123 1.0466870 372.2838 4.248311 1982-1999 –> 2000-2009 2.385867 1096.0581
Atlantic (-Inf,26] 4 cstar_tref ~ sal + temp + oxygen + silicate 123 1.1049811 385.6166 4.375130 1982-1999 –> 2000-2009 2.462633 1115.0898
Atlantic (-Inf,26] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 123 1.0466758 374.2812 4.241859 1982-1999 –> 2000-2009 2.384687 1099.6923
Atlantic (-Inf,26] 4 cstar_tref ~ sal + temp + phosphate + phosphate_star 123 1.0466870 372.2838 4.248311 1982-1999 –> 2000-2009 2.385867 1096.0581
Atlantic (-Inf,26] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 123 1.0466758 374.2812 4.241859 1982-1999 –> 2000-2009 2.384687 1099.6923
Atlantic (-Inf,26] 3 cstar_tref ~ sal + temp + aou 107 1.3371422 375.8273 3.358010 2000-2009 –> 2010-2019 2.486875 769.2104
Atlantic (-Inf,26] 4 cstar_tref ~ sal + temp + aou + phosphate_star 107 1.2091891 356.3021 3.722170 2000-2009 –> 2010-2019 2.253961 728.1354
Atlantic (-Inf,26] 4 cstar_tref ~ sal + temp + aou + silicate 107 1.3232430 375.5911 3.424852 2000-2009 –> 2010-2019 2.427412 761.0270
Atlantic (-Inf,26] 3 cstar_tref ~ sal + temp + oxygen 107 1.3360164 375.6470 3.360126 2000-2009 –> 2010-2019 2.485683 769.0160
Atlantic (-Inf,26] 4 cstar_tref ~ sal + temp + oxygen + phosphate_star 107 1.2098022 356.4106 3.758556 2000-2009 –> 2010-2019 2.256489 728.6945
Atlantic (-Inf,26] 4 cstar_tref ~ sal + temp + oxygen + silicate 107 1.3222101 375.4240 3.475512 2000-2009 –> 2010-2019 2.427191 761.0407
Atlantic (-Inf,26] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 107 1.1463671 346.8848 3.464326 2000-2009 –> 2010-2019 2.193043 721.1660
Atlantic (-Inf,26] 4 cstar_tref ~ sal + temp + phosphate + phosphate_star 107 1.2098022 356.4106 3.758556 2000-2009 –> 2010-2019 2.256489 728.6945
Atlantic (-Inf,26] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 107 1.1463671 346.8848 3.464326 2000-2009 –> 2010-2019 2.193043 721.1660
Atlantic (-Inf,26] 4 cstar_tref ~ sal + temp + silicate + phosphate_star 107 1.3634664 381.9993 4.209507 2000-2009 –> 2010-2019 2.972035 859.9930
Atlantic (26,26.5] 3 cstar_tref ~ aou + silicate + phosphate 821 5.1945928 5045.2862 16.170517 1982-1999 –> 2000-2009 10.184744 14285.8205
Atlantic (26,26.5] 3 cstar_tref ~ aou + silicate + phosphate_star 821 5.3344024 5088.8955 16.789997 1982-1999 –> 2000-2009 10.489154 14428.4105
Atlantic (26,26.5] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 821 5.2776207 5073.3236 15.936910 1982-1999 –> 2000-2009 10.363762 14373.9703
Atlantic (26,26.5] 4 cstar_tref ~ sal + oxygen + silicate + phosphate_star 821 5.3496503 5095.5823 15.926798 1982-1999 –> 2000-2009 10.561302 14470.5803
Atlantic (26,26.5] 3 cstar_tref ~ temp + aou + silicate 821 5.2575413 5065.0645 15.499679 1982-1999 –> 2000-2009 10.481147 14445.0502
Atlantic (26,26.5] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 821 4.5374656 4825.2063 12.323404 1982-1999 –> 2000-2009 8.908995 13664.0630
Atlantic (26,26.5] 3 cstar_tref ~ temp + oxygen + silicate 821 5.2546632 5064.1654 15.530795 1982-1999 –> 2000-2009 10.467423 14437.8117
Atlantic (26,26.5] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 821 4.5329123 4823.5578 12.329884 1982-1999 –> 2000-2009 8.897583 13657.6258
Atlantic (26,26.5] 4 cstar_tref ~ temp + silicate + phosphate + phosphate_star 821 4.5329123 4823.5578 12.329884 1982-1999 –> 2000-2009 8.897583 13657.6258
Atlantic (26,26.5] 3 cstar_tref ~ temp + silicate + phosphate_star 821 4.6582673 4866.3498 12.531600 1982-1999 –> 2000-2009 9.423122 13965.9767
Atlantic (26,26.5] 3 cstar_tref ~ aou + silicate + phosphate 792 4.9672752 4796.5470 15.371700 2000-2009 –> 2010-2019 10.161868 9841.8332
Atlantic (26,26.5] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 792 5.0353902 4820.1204 15.184840 2000-2009 –> 2010-2019 10.313011 9893.4440
Atlantic (26,26.5] 3 cstar_tref ~ temp + aou + silicate 792 5.0241506 4814.5808 15.233427 2000-2009 –> 2010-2019 10.281692 9879.6453
Atlantic (26,26.5] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 792 4.4278377 4616.4502 12.143959 2000-2009 –> 2010-2019 8.965303 9441.6566
Atlantic (26,26.5] 3 cstar_tref ~ temp + oxygen + silicate 792 5.0247838 4814.7804 15.243114 2000-2009 –> 2010-2019 10.279447 9878.9458
Atlantic (26,26.5] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 792 4.4258792 4615.7495 12.162944 2000-2009 –> 2010-2019 8.958792 9439.3073
Atlantic (26,26.5] 2 cstar_tref ~ temp + silicate 792 5.0252108 4812.9150 15.257003 2000-2009 –> 2010-2019 10.294538 9879.6564
Atlantic (26,26.5] 3 cstar_tref ~ temp + silicate + phosphate 792 4.9299870 4784.6114 16.804661 2000-2009 –> 2010-2019 10.155756 9839.7228
Atlantic (26,26.5] 4 cstar_tref ~ temp + silicate + phosphate + phosphate_star 792 4.4258792 4615.7495 12.162944 2000-2009 –> 2010-2019 8.958792 9439.3073
Atlantic (26,26.5] 3 cstar_tref ~ temp + silicate + phosphate_star 792 4.4546547 4624.0147 12.125009 2000-2009 –> 2010-2019 9.112922 9490.3645
Atlantic (26.5,26.75] 3 cstar_tref ~ sal + aou + phosphate_star 1062 4.3974032 6169.4996 13.278488 1982-1999 –> 2000-2009 8.356222 16783.2899
Atlantic (26.5,26.75] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 1062 3.4176049 5636.0940 9.697426 1982-1999 –> 2000-2009 6.757867 15607.2974
Atlantic (26.5,26.75] 3 cstar_tref ~ sal + oxygen + phosphate_star 1062 4.4237693 6182.1967 13.487430 1982-1999 –> 2000-2009 8.419835 16831.5165
Atlantic (26.5,26.75] 4 cstar_tref ~ sal + oxygen + silicate + phosphate_star 1062 3.4901335 5680.6980 9.973189 1982-1999 –> 2000-2009 6.909761 15740.9935
Atlantic (26.5,26.75] 3 cstar_tref ~ sal + phosphate + phosphate_star 1062 4.4237693 6182.1967 13.487430 1982-1999 –> 2000-2009 8.419835 16831.5165
Atlantic (26.5,26.75] 4 cstar_tref ~ sal + silicate + phosphate + phosphate_star 1062 3.4901335 5680.6980 9.973189 1982-1999 –> 2000-2009 6.909761 15740.9935
Atlantic (26.5,26.75] 3 cstar_tref ~ temp + aou + phosphate_star 1062 4.5962890 6263.4548 12.965493 1982-1999 –> 2000-2009 8.756120 17065.1577
Atlantic (26.5,26.75] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 1062 3.2809089 5549.3934 10.744534 1982-1999 –> 2000-2009 6.475183 15351.0457
Atlantic (26.5,26.75] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 1062 3.2889038 5554.5628 10.782412 1982-1999 –> 2000-2009 6.493250 15368.1593
Atlantic (26.5,26.75] 4 cstar_tref ~ temp + silicate + phosphate + phosphate_star 1062 3.2889038 5554.5628 10.782412 1982-1999 –> 2000-2009 6.493250 15368.1593
Atlantic (26.5,26.75] 3 cstar_tref ~ sal + aou + phosphate_star 991 4.9676380 5999.3721 17.576055 2000-2009 –> 2010-2019 9.365041 12168.8717
Atlantic (26.5,26.75] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 991 3.6595917 5395.6870 10.764565 2000-2009 –> 2010-2019 7.077197 11031.7810
Atlantic (26.5,26.75] 3 cstar_tref ~ sal + oxygen + phosphate_star 991 4.9778260 6003.4328 17.834646 2000-2009 –> 2010-2019 9.401595 12185.6295
Atlantic (26.5,26.75] 4 cstar_tref ~ sal + oxygen + silicate + phosphate_star 991 3.7067219 5421.0493 11.179730 2000-2009 –> 2010-2019 7.196855 11101.7473
Atlantic (26.5,26.75] 3 cstar_tref ~ sal + phosphate + phosphate_star 991 4.9778260 6003.4328 17.834646 2000-2009 –> 2010-2019 9.401595 12185.6295
Atlantic (26.5,26.75] 4 cstar_tref ~ sal + silicate + phosphate + phosphate_star 991 3.7067219 5421.0493 11.179730 2000-2009 –> 2010-2019 7.196855 11101.7473
Atlantic (26.5,26.75] 3 cstar_tref ~ temp + aou + phosphate_star 991 5.2643531 6114.3554 14.876368 2000-2009 –> 2010-2019 9.860642 12377.8103
Atlantic (26.5,26.75] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 991 3.6944860 5414.4959 11.677687 2000-2009 –> 2010-2019 6.975395 10963.8893
Atlantic (26.5,26.75] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 991 3.7013236 5418.1607 11.679905 2000-2009 –> 2010-2019 6.990227 10972.7235
Atlantic (26.5,26.75] 4 cstar_tref ~ temp + silicate + phosphate + phosphate_star 991 3.7013236 5418.1607 11.679905 2000-2009 –> 2010-2019 6.990227 10972.7235
Atlantic (26.75,27] 3 cstar_tref ~ sal + aou + phosphate_star 1823 4.0543880 10287.1198 33.076058 1982-1999 –> 2000-2009 7.390442 26774.1270
Atlantic (26.75,27] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 1823 2.4866609 8506.7400 34.126676 1982-1999 –> 2000-2009 4.737788 22525.4598
Atlantic (26.75,27] 3 cstar_tref ~ sal + oxygen + phosphate_star 1823 4.0612365 10293.2732 34.597274 1982-1999 –> 2000-2009 7.412792 26809.3936
Atlantic (26.75,27] 4 cstar_tref ~ sal + oxygen + silicate + phosphate_star 1823 2.5970905 8665.1623 38.485667 1982-1999 –> 2000-2009 4.955318 22975.7732
Atlantic (26.75,27] 3 cstar_tref ~ sal + phosphate + phosphate_star 1823 4.0612365 10293.2732 34.597274 1982-1999 –> 2000-2009 7.412792 26809.3936
Atlantic (26.75,27] 4 cstar_tref ~ sal + silicate + phosphate + phosphate_star 1823 2.5970905 8665.1623 38.485667 1982-1999 –> 2000-2009 4.955318 22975.7732
Atlantic (26.75,27] 3 cstar_tref ~ temp + aou + phosphate_star 1823 4.6423383 10780.8554 23.266594 1982-1999 –> 2000-2009 8.384491 27989.2601
Atlantic (26.75,27] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 1823 2.2467361 8136.8086 16.994135 1982-1999 –> 2000-2009 4.085697 20885.5216
Atlantic (26.75,27] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 1823 2.2781392 8187.4167 19.889087 1982-1999 –> 2000-2009 4.136064 21000.5619
Atlantic (26.75,27] 4 cstar_tref ~ temp + silicate + phosphate + phosphate_star 1823 2.2781392 8187.4167 19.889087 1982-1999 –> 2000-2009 4.136064 21000.5619
Atlantic (26.75,27] 3 cstar_tref ~ sal + aou + phosphate_star 1838 4.6423331 10869.4760 12.826802 2000-2009 –> 2010-2019 8.696721 21156.5958
Atlantic (26.75,27] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 1838 2.5464232 8663.9374 23.624288 2000-2009 –> 2010-2019 5.033084 17170.6774
Atlantic (26.75,27] 3 cstar_tref ~ sal + oxygen + phosphate_star 1838 4.6403865 10867.9343 12.906049 2000-2009 –> 2010-2019 8.701623 21161.2075
Atlantic (26.75,27] 4 cstar_tref ~ sal + oxygen + silicate + phosphate_star 1838 2.5954926 8734.0998 23.199581 2000-2009 –> 2010-2019 5.192583 17399.2620
Atlantic (26.75,27] 3 cstar_tref ~ sal + phosphate + phosphate_star 1838 4.6403865 10867.9343 12.906049 2000-2009 –> 2010-2019 8.701623 21161.2075
Atlantic (26.75,27] 4 cstar_tref ~ sal + silicate + phosphate + phosphate_star 1838 2.5954926 8734.0998 23.199581 2000-2009 –> 2010-2019 5.192583 17399.2620
Atlantic (26.75,27] 3 cstar_tref ~ sal + silicate + phosphate_star 1838 4.3522653 10632.2983 14.511306 2000-2009 –> 2010-2019 8.882745 21324.2265
Atlantic (26.75,27] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 1838 2.6580468 8821.6446 28.396389 2000-2009 –> 2010-2019 4.904783 16958.4533
Atlantic (26.75,27] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 1838 2.6744528 8844.2639 28.451846 2000-2009 –> 2010-2019 4.952592 17031.6806
Atlantic (26.75,27] 4 cstar_tref ~ temp + silicate + phosphate + phosphate_star 1838 2.6744528 8844.2639 28.451846 2000-2009 –> 2010-2019 4.952592 17031.6806
Atlantic (27,27.25] 3 cstar_tref ~ sal + aou + phosphate_star 1559 2.9909389 7850.2917 12.853728 1982-1999 –> 2000-2009 5.728793 20529.4529
Atlantic (27,27.25] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 1559 2.1052588 6757.4093 9.766882 1982-1999 –> 2000-2009 4.368811 18445.1436
Atlantic (27,27.25] 3 cstar_tref ~ sal + oxygen + phosphate_star 1559 2.9894022 7848.6893 12.800944 1982-1999 –> 2000-2009 5.763133 20595.8352
Atlantic (27,27.25] 4 cstar_tref ~ sal + oxygen + silicate + phosphate_star 1559 2.1608863 6838.7269 10.174447 1982-1999 –> 2000-2009 4.553159 18815.2804
Atlantic (27,27.25] 3 cstar_tref ~ sal + phosphate + phosphate_star 1559 2.9894022 7848.6893 12.800944 1982-1999 –> 2000-2009 5.763133 20595.8352
Atlantic (27,27.25] 4 cstar_tref ~ sal + silicate + phosphate + phosphate_star 1559 2.1608863 6838.7269 10.174447 1982-1999 –> 2000-2009 4.553159 18815.2804
Atlantic (27,27.25] 3 cstar_tref ~ sal + silicate + phosphate_star 1559 2.7013086 7532.7202 12.000646 1982-1999 –> 2000-2009 5.842569 20929.6444
Atlantic (27,27.25] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 1559 2.5907112 7404.3757 10.397666 1982-1999 –> 2000-2009 4.712393 18754.1089
Atlantic (27,27.25] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 1559 2.6024440 7418.4645 10.395423 1982-1999 –> 2000-2009 4.739400 18805.6571
Atlantic (27,27.25] 4 cstar_tref ~ temp + silicate + phosphate + phosphate_star 1559 2.6024440 7418.4645 10.395423 1982-1999 –> 2000-2009 4.739400 18805.6571
Atlantic (27,27.25] 3 cstar_tref ~ sal + aou + phosphate_star 1614 3.4086962 8548.9264 13.892510 2000-2009 –> 2010-2019 6.399635 16399.2181
Atlantic (27,27.25] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 1614 2.2884650 7264.7344 14.940040 2000-2009 –> 2010-2019 4.393724 14022.1437
Atlantic (27,27.25] 3 cstar_tref ~ sal + oxygen + phosphate_star 1614 3.4100421 8550.2007 13.797278 2000-2009 –> 2010-2019 6.399444 16398.8901
Atlantic (27,27.25] 4 cstar_tref ~ sal + oxygen + silicate + phosphate_star 1614 2.3175306 7305.4747 15.091983 2000-2009 –> 2010-2019 4.478417 14144.2016
Atlantic (27,27.25] 3 cstar_tref ~ sal + phosphate + phosphate_star 1614 3.4100421 8550.2007 13.797278 2000-2009 –> 2010-2019 6.399444 16398.8901
Atlantic (27,27.25] 4 cstar_tref ~ sal + silicate + phosphate + phosphate_star 1614 2.3175306 7305.4747 15.091983 2000-2009 –> 2010-2019 4.478417 14144.2016
Atlantic (27,27.25] 3 cstar_tref ~ sal + silicate + phosphate_star 1614 2.7224187 7823.2425 13.516236 2000-2009 –> 2010-2019 5.423727 15355.9627
Atlantic (27,27.25] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 1614 3.0365539 8177.7483 14.563534 2000-2009 –> 2010-2019 5.627265 15582.1240
Atlantic (27,27.25] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 1614 3.0472502 8189.0991 14.580528 2000-2009 –> 2010-2019 5.649694 15607.5636
Atlantic (27,27.25] 4 cstar_tref ~ temp + silicate + phosphate + phosphate_star 1614 3.0472502 8189.0991 14.580528 2000-2009 –> 2010-2019 5.649694 15607.5636
Atlantic (27.25,27.5] 2 cstar_tref ~ aou + phosphate 1594 7.8584666 11103.9297 58.627415 1982-1999 –> 2000-2009 14.733733 30376.5073
Atlantic (27.25,27.5] 3 cstar_tref ~ sal + aou + phosphate_star 1594 3.9672368 8926.8627 25.819532 1982-1999 –> 2000-2009 7.820120 24868.0902
Atlantic (27.25,27.5] 3 cstar_tref ~ sal + oxygen + phosphate_star 1594 3.8976689 8870.4632 25.177551 1982-1999 –> 2000-2009 7.734994 24788.4008
Atlantic (27.25,27.5] 3 cstar_tref ~ sal + phosphate + phosphate_star 1594 3.8976689 8870.4632 25.177551 1982-1999 –> 2000-2009 7.734994 24788.4008
Atlantic (27.25,27.5] 2 cstar_tref ~ sal + phosphate_star 1594 4.4216101 9270.5505 30.989773 1982-1999 –> 2000-2009 8.286016 25226.9679
Atlantic (27.25,27.5] 3 cstar_tref ~ temp + aou + phosphate 1594 3.9237702 8891.7409 14.417675 1982-1999 –> 2000-2009 6.927036 23399.0268
Atlantic (27.25,27.5] 3 cstar_tref ~ temp + aou + phosphate_star 1594 4.0402761 8985.0220 14.433386 1982-1999 –> 2000-2009 7.160894 23712.9380
Atlantic (27.25,27.5] 3 cstar_tref ~ temp + oxygen + phosphate_star 1594 4.0317694 8978.3027 14.436277 1982-1999 –> 2000-2009 7.153495 23708.2623
Atlantic (27.25,27.5] 3 cstar_tref ~ temp + phosphate + phosphate_star 1594 4.0317694 8978.3027 14.436277 1982-1999 –> 2000-2009 7.153495 23708.2623
Atlantic (27.25,27.5] 2 cstar_tref ~ temp + phosphate_star 1594 4.4074198 9260.3027 17.320916 1982-1999 –> 2000-2009 7.537261 24003.2079
Atlantic (27.25,27.5] 3 cstar_tref ~ sal + aou + phosphate 1540 3.8327915 8518.5983 24.636192 2000-2009 –> 2010-2019 7.038028 16765.5113
Atlantic (27.25,27.5] 3 cstar_tref ~ sal + aou + phosphate_star 1540 4.4005717 8944.0728 28.470202 2000-2009 –> 2010-2019 8.367809 17870.9355
Atlantic (27.25,27.5] 3 cstar_tref ~ sal + oxygen + phosphate_star 1540 4.3056986 8876.9440 27.659517 2000-2009 –> 2010-2019 8.203367 17747.4072
Atlantic (27.25,27.5] 3 cstar_tref ~ sal + phosphate + phosphate_star 1540 4.3056986 8876.9440 27.659517 2000-2009 –> 2010-2019 8.203367 17747.4072
Atlantic (27.25,27.5] 4 cstar_tref ~ sal + silicate + phosphate + phosphate_star 1540 3.7958826 8490.7949 28.606538 2000-2009 –> 2010-2019 7.148711 16883.2269
Atlantic (27.25,27.5] 3 cstar_tref ~ sal + silicate + phosphate_star 1540 3.8176963 8506.4440 29.476780 2000-2009 –> 2010-2019 7.187308 16912.7944
Atlantic (27.25,27.5] 3 cstar_tref ~ temp + aou + phosphate 1540 4.5656232 9057.4802 16.030129 2000-2009 –> 2010-2019 8.489393 17949.2211
Atlantic (27.25,27.5] 3 cstar_tref ~ temp + oxygen + phosphate 1540 4.6812383 9134.5037 16.681978 2000-2009 –> 2010-2019 8.713008 18112.8064
Atlantic (27.25,27.5] 4 cstar_tref ~ temp + silicate + phosphate + phosphate_star 1540 4.2290773 8823.6409 18.391600 2000-2009 –> 2010-2019 7.760756 17381.7501
Atlantic (27.25,27.5] 3 cstar_tref ~ temp + silicate + phosphate_star 1540 4.2407045 8830.0972 19.031159 2000-2009 –> 2010-2019 7.800389 17411.3878
Atlantic (27.5,27.75] 3 cstar_tref ~ sal + aou + phosphate 2189 2.7630085 10671.5623 13.627639 1982-1999 –> 2000-2009 5.238058 28046.1489
Atlantic (27.5,27.75] 3 cstar_tref ~ sal + aou + phosphate_star 2189 3.4919682 11696.6510 13.181354 1982-1999 –> 2000-2009 7.165867 32021.0449
Atlantic (27.5,27.75] 3 cstar_tref ~ sal + oxygen + phosphate 2189 3.3798413 11553.7672 12.897114 1982-1999 –> 2000-2009 6.986569 31740.3570
Atlantic (27.5,27.75] 3 cstar_tref ~ sal + oxygen + phosphate_star 2189 3.3798413 11553.7672 12.897114 1982-1999 –> 2000-2009 6.986569 31740.3570
Atlantic (27.5,27.75] 3 cstar_tref ~ sal + phosphate + phosphate_star 2189 3.3798413 11553.7672 12.897114 1982-1999 –> 2000-2009 6.986569 31740.3570
Atlantic (27.5,27.75] 3 cstar_tref ~ sal + silicate + phosphate_star 2189 3.2670012 11405.1067 14.068179 1982-1999 –> 2000-2009 6.803351 31444.5336
Atlantic (27.5,27.75] 3 cstar_tref ~ temp + aou + phosphate 2189 4.0561151 12352.3008 17.450143 1982-1999 –> 2000-2009 7.104849 31283.7099
Atlantic (27.5,27.75] 3 cstar_tref ~ temp + oxygen + phosphate_star 2189 4.1397026 12441.6045 17.120635 1982-1999 –> 2000-2009 7.265273 31558.8946
Atlantic (27.5,27.75] 3 cstar_tref ~ temp + phosphate + phosphate_star 2189 4.1397026 12441.6045 17.120635 1982-1999 –> 2000-2009 7.265273 31558.8946
Atlantic (27.5,27.75] 3 cstar_tref ~ temp + silicate + phosphate_star 2189 3.6966241 11945.9976 18.229750 1982-1999 –> 2000-2009 6.601169 30515.5854
Atlantic (27.5,27.75] 3 cstar_tref ~ sal + aou + phosphate 2130 3.2316574 11051.6374 13.890500 2000-2009 –> 2010-2019 5.994666 21723.1998
Atlantic (27.5,27.75] 3 cstar_tref ~ sal + aou + phosphate_star 2130 3.7574472 11693.8097 14.526143 2000-2009 –> 2010-2019 7.249415 23390.4607
Atlantic (27.5,27.75] 3 cstar_tref ~ sal + oxygen + phosphate 2130 3.6208038 11536.0033 14.181762 2000-2009 –> 2010-2019 7.000645 23089.7705
Atlantic (27.5,27.75] 3 cstar_tref ~ sal + oxygen + phosphate_star 2130 3.6208038 11536.0033 14.181762 2000-2009 –> 2010-2019 7.000645 23089.7705
Atlantic (27.5,27.75] 3 cstar_tref ~ sal + phosphate + phosphate_star 2130 3.6208038 11536.0033 14.181762 2000-2009 –> 2010-2019 7.000645 23089.7705
Atlantic (27.5,27.75] 2 cstar_tref ~ sal + phosphate_star 2130 4.7270661 12669.7563 22.187512 2000-2009 –> 2010-2019 8.725610 24957.4719
Atlantic (27.5,27.75] 4 cstar_tref ~ sal + silicate + phosphate + phosphate_star 2130 3.2350841 11058.1521 15.479156 2000-2009 –> 2010-2019 6.418275 22351.4821
Atlantic (27.5,27.75] 3 cstar_tref ~ sal + silicate + phosphate_star 2130 3.3925320 11258.5942 18.351277 2000-2009 –> 2010-2019 6.659533 22663.7010
Atlantic (27.5,27.75] 4 cstar_tref ~ temp + silicate + phosphate + phosphate_star 2130 4.2910702 12261.5222 19.836394 2000-2009 –> 2010-2019 7.974643 24194.0359
Atlantic (27.5,27.75] 3 cstar_tref ~ temp + silicate + phosphate_star 2130 4.3503079 12317.9287 21.259027 2000-2009 –> 2010-2019 8.046932 24263.9264
Atlantic (27.75,27.85] 2 cstar_tref ~ aou + phosphate 779 5.1790958 4781.0405 21.338901 1982-1999 –> 2000-2009 9.413988 12345.5033
Atlantic (27.75,27.85] 2 cstar_tref ~ aou + phosphate_star 779 6.5119890 5137.8450 25.195117 1982-1999 –> 2000-2009 12.003079 13388.0976
Atlantic (27.75,27.85] 2 cstar_tref ~ oxygen + phosphate 779 5.7010391 4930.6365 21.604655 1982-1999 –> 2000-2009 10.603747 12881.6737
Atlantic (27.75,27.85] 2 cstar_tref ~ oxygen + phosphate_star 779 5.7010391 4930.6365 21.604655 1982-1999 –> 2000-2009 10.603747 12881.6737
Atlantic (27.75,27.85] 2 cstar_tref ~ phosphate + phosphate_star 779 5.7010391 4930.6365 21.604655 1982-1999 –> 2000-2009 10.603747 12881.6737
Atlantic (27.75,27.85] 3 cstar_tref ~ sal + aou + phosphate 779 2.2543220 3487.1254 11.487811 1982-1999 –> 2000-2009 4.228110 9038.2025
Atlantic (27.75,27.85] 2 cstar_tref ~ sal + phosphate_star 779 3.0722686 3967.4307 11.746784 1982-1999 –> 2000-2009 5.781594 10352.7152
Atlantic (27.75,27.85] 2 cstar_tref ~ silicate + phosphate 779 6.4392248 5120.3381 19.853324 1982-1999 –> 2000-2009 12.941463 13816.8027
Atlantic (27.75,27.85] 3 cstar_tref ~ temp + aou + phosphate 779 2.4686857 3628.6489 14.823817 1982-1999 –> 2000-2009 4.479538 9228.8414
Atlantic (27.75,27.85] 2 cstar_tref ~ temp + phosphate_star 779 2.5483968 3676.1598 14.886774 1982-1999 –> 2000-2009 4.648977 9389.6011
Atlantic (27.75,27.85] 2 cstar_tref ~ aou + phosphate 770 5.5628985 4835.9890 26.633822 2000-2009 –> 2010-2019 10.741994 9617.0296
Atlantic (27.75,27.85] 2 cstar_tref ~ aou + phosphate_star 770 6.7389900 5131.3468 31.176281 2000-2009 –> 2010-2019 13.250979 10269.1919
Atlantic (27.75,27.85] 2 cstar_tref ~ oxygen + phosphate 770 5.9380521 4936.4923 27.114589 2000-2009 –> 2010-2019 11.639091 9867.1288
Atlantic (27.75,27.85] 2 cstar_tref ~ oxygen + phosphate_star 770 5.9380521 4936.4923 27.114589 2000-2009 –> 2010-2019 11.639091 9867.1288
Atlantic (27.75,27.85] 2 cstar_tref ~ phosphate + phosphate_star 770 5.9380521 4936.4923 27.114589 2000-2009 –> 2010-2019 11.639091 9867.1288
Atlantic (27.75,27.85] 3 cstar_tref ~ sal + aou + phosphate 770 2.5533119 3638.7479 10.680820 2000-2009 –> 2010-2019 4.807634 7125.8733
Atlantic (27.75,27.85] 2 cstar_tref ~ sal + phosphate_star 770 3.3509189 4055.3866 10.634616 2000-2009 –> 2010-2019 6.423188 8022.8174
Atlantic (27.75,27.85] 2 cstar_tref ~ silicate + phosphate 770 6.2474442 5014.7109 26.041238 2000-2009 –> 2010-2019 12.686669 10135.0491
Atlantic (27.75,27.85] 2 cstar_tref ~ temp + phosphate_star 770 2.9975532 3883.7717 14.333073 2000-2009 –> 2010-2019 5.545950 7559.9316
Atlantic (27.75,27.85] 3 cstar_tref ~ temp + silicate + phosphate 770 6.2444059 5015.9618 26.357040 2000-2009 –> 2010-2019 12.627398 10124.6343
Atlantic (27.85,27.95] 2 cstar_tref ~ aou + phosphate 738 4.9607995 4466.2660 26.713923 1982-1999 –> 2000-2009 9.839694 12500.5723
Atlantic (27.85,27.95] 2 cstar_tref ~ aou + phosphate_star 738 5.8835133 4718.0527 28.376497 1982-1999 –> 2000-2009 11.611860 13181.2372
Atlantic (27.85,27.95] 2 cstar_tref ~ oxygen + phosphate 738 5.0599664 4495.4804 26.366275 1982-1999 –> 2000-2009 10.113294 12623.6495
Atlantic (27.85,27.95] 2 cstar_tref ~ oxygen + phosphate_star 738 5.0599664 4495.4804 26.366275 1982-1999 –> 2000-2009 10.113294 12623.6495
Atlantic (27.85,27.95] 2 cstar_tref ~ phosphate + phosphate_star 738 5.0599664 4495.4804 26.366275 1982-1999 –> 2000-2009 10.113294 12623.6495
Atlantic (27.85,27.95] 2 cstar_tref ~ sal + phosphate 738 5.0189787 4483.4755 24.856453 1982-1999 –> 2000-2009 10.197946 12677.2655
Atlantic (27.85,27.95] 2 cstar_tref ~ sal + phosphate_star 738 5.8311952 4704.8689 98.149901 1982-1999 –> 2000-2009 9.463385 11950.7140
Atlantic (27.85,27.95] 2 cstar_tref ~ silicate + phosphate 738 5.0797621 4501.2436 26.251893 1982-1999 –> 2000-2009 10.235835 12683.1956
Atlantic (27.85,27.95] 2 cstar_tref ~ temp + phosphate 738 5.0488960 4492.2476 25.145689 1982-1999 –> 2000-2009 10.230903 12687.6054
Atlantic (27.85,27.95] 2 cstar_tref ~ temp + phosphate_star 738 3.8524384 4093.0438 61.562426 1982-1999 –> 2000-2009 6.901410 10871.2052
Atlantic (27.85,27.95] 2 cstar_tref ~ aou + phosphate 814 6.1156373 5266.0941 29.554772 2000-2009 –> 2010-2019 11.076437 9732.3601
Atlantic (27.85,27.95] 2 cstar_tref ~ aou + phosphate_star 814 6.8856339 5459.1557 30.767144 2000-2009 –> 2010-2019 12.769147 10177.2084
Atlantic (27.85,27.95] 2 cstar_tref ~ oxygen + phosphate 814 6.1006190 5262.0912 29.183001 2000-2009 –> 2010-2019 11.160585 9757.5717
Atlantic (27.85,27.95] 2 cstar_tref ~ oxygen + phosphate_star 814 6.1006190 5262.0912 29.183001 2000-2009 –> 2010-2019 11.160585 9757.5717
Atlantic (27.85,27.95] 2 cstar_tref ~ phosphate + phosphate_star 814 6.1006190 5262.0912 29.183001 2000-2009 –> 2010-2019 11.160585 9757.5717
Atlantic (27.85,27.95] 2 cstar_tref ~ sal + phosphate 814 5.8061161 5181.5404 27.139019 2000-2009 –> 2010-2019 10.825095 9665.0160
Atlantic (27.85,27.95] 2 cstar_tref ~ sal + phosphate_star 814 4.2196120 4661.9338 35.964033 2000-2009 –> 2010-2019 10.050807 9366.8028
Atlantic (27.85,27.95] 2 cstar_tref ~ silicate + phosphate 814 6.0768848 5255.7452 28.861665 2000-2009 –> 2010-2019 11.156647 9756.9888
Atlantic (27.85,27.95] 2 cstar_tref ~ temp + phosphate 814 5.8781258 5201.6073 27.909955 2000-2009 –> 2010-2019 10.927022 9693.8550
Atlantic (27.85,27.95] 2 cstar_tref ~ temp + phosphate_star 814 3.8956603 4531.8892 23.958127 2000-2009 –> 2010-2019 7.748099 8624.9330
Atlantic (27.95,28.05] 2 cstar_tref ~ aou + phosphate 924 6.2743271 6023.9880 36.996618 1982-1999 –> 2000-2009 12.479755 16907.0923
Atlantic (27.95,28.05] 2 cstar_tref ~ aou + phosphate_star 924 6.9139313 6203.3774 39.264370 1982-1999 –> 2000-2009 13.744737 17408.3349
Atlantic (27.95,28.05] 2 cstar_tref ~ oxygen + phosphate 924 6.2775881 6024.9483 36.567293 1982-1999 –> 2000-2009 12.599096 16970.1764
Atlantic (27.95,28.05] 2 cstar_tref ~ oxygen + phosphate_star 924 6.2775881 6024.9483 36.567293 1982-1999 –> 2000-2009 12.599096 16970.1764
Atlantic (27.95,28.05] 2 cstar_tref ~ phosphate + phosphate_star 924 6.2775881 6024.9483 36.567293 1982-1999 –> 2000-2009 12.599096 16970.1764
Atlantic (27.95,28.05] 3 cstar_tref ~ sal + aou + phosphate 924 4.8981822 5568.4194 27.990154 1982-1999 –> 2000-2009 9.820542 15677.0790
Atlantic (27.95,28.05] 2 cstar_tref ~ sal + phosphate 924 5.7567977 5864.9032 31.615177 1982-1999 –> 2000-2009 12.018371 16778.1995
Atlantic (27.95,28.05] 2 cstar_tref ~ sal + silicate 924 8.0377240 6481.7001 34.344975 1982-1999 –> 2000-2009 16.236414 18298.5043
Atlantic (27.95,28.05] 2 cstar_tref ~ silicate + phosphate_star 924 8.4054667 6564.3729 41.317363 1982-1999 –> 2000-2009 16.333300 18268.5677
Atlantic (27.95,28.05] 2 cstar_tref ~ temp + phosphate 924 5.9733795 5933.1525 32.923149 1982-1999 –> 2000-2009 12.370692 16918.3376
Atlantic (27.95,28.05] 2 cstar_tref ~ aou + phosphate 946 6.2629535 6163.7931 36.037517 2000-2009 –> 2010-2019 12.537281 12187.7811
Atlantic (27.95,28.05] 2 cstar_tref ~ oxygen + phosphate 946 6.1851014 6140.1270 35.219503 2000-2009 –> 2010-2019 12.462689 12165.0753
Atlantic (27.95,28.05] 2 cstar_tref ~ oxygen + phosphate_star 946 6.1851014 6140.1270 35.219503 2000-2009 –> 2010-2019 12.462689 12165.0753
Atlantic (27.95,28.05] 2 cstar_tref ~ phosphate + phosphate_star 946 6.1851014 6140.1270 35.219503 2000-2009 –> 2010-2019 12.462689 12165.0753
Atlantic (27.95,28.05] 3 cstar_tref ~ sal + aou + phosphate 946 4.5757533 5571.9311 26.530895 2000-2009 –> 2010-2019 9.473935 11140.3505
Atlantic (27.95,28.05] 3 cstar_tref ~ sal + aou + phosphate_star 946 4.8654512 5688.0774 27.864577 2000-2009 –> 2010-2019 10.046029 11360.0817
Atlantic (27.95,28.05] 3 cstar_tref ~ sal + oxygen + phosphate 946 4.6602564 5606.5530 27.040106 2000-2009 –> 2010-2019 9.656702 11211.6785
Atlantic (27.95,28.05] 3 cstar_tref ~ sal + oxygen + phosphate_star 946 4.6602564 5606.5530 27.040106 2000-2009 –> 2010-2019 9.656702 11211.6785
Atlantic (27.95,28.05] 2 cstar_tref ~ sal + phosphate 946 5.2562325 5832.2440 29.159702 2000-2009 –> 2010-2019 11.013030 11697.1471
Atlantic (27.95,28.05] 2 cstar_tref ~ temp + phosphate 946 5.5962486 5950.8382 30.367975 2000-2009 –> 2010-2019 11.569628 11883.9908
Atlantic (28.05,28.1] 2 cstar_tref ~ aou + phosphate 622 2.3959938 2860.1644 25.743477 1982-1999 –> 2000-2009 4.186544 7175.3100
Atlantic (28.05,28.1] 2 cstar_tref ~ aou + phosphate_star 622 3.0188673 3147.6324 31.150128 1982-1999 –> 2000-2009 5.435739 8108.2717
Atlantic (28.05,28.1] 2 cstar_tref ~ aou + silicate 622 4.3215042 3593.8823 33.438659 1982-1999 –> 2000-2009 8.665656 9816.3615
Atlantic (28.05,28.1] 2 cstar_tref ~ oxygen + phosphate 622 2.4368732 2881.2099 25.878079 1982-1999 –> 2000-2009 4.378103 7370.2344
Atlantic (28.05,28.1] 2 cstar_tref ~ oxygen + phosphate_star 622 2.4368732 2881.2099 25.878079 1982-1999 –> 2000-2009 4.378103 7370.2344
Atlantic (28.05,28.1] 2 cstar_tref ~ oxygen + silicate 622 4.3392134 3598.9697 33.902428 1982-1999 –> 2000-2009 8.669438 9814.5383
Atlantic (28.05,28.1] 2 cstar_tref ~ sal + phosphate 622 2.2376316 2775.0995 19.566535 1982-1999 –> 2000-2009 4.457994 7553.2422
Atlantic (28.05,28.1] 2 cstar_tref ~ sal + silicate 622 4.3812585 3610.9655 33.241155 1982-1999 –> 2000-2009 8.713883 9827.7265
Atlantic (28.05,28.1] 2 cstar_tref ~ silicate + phosphate_star 622 3.8499317 3450.1405 33.080202 1982-1999 –> 2000-2009 7.199874 9113.3417
Atlantic (28.05,28.1] 2 cstar_tref ~ temp + phosphate 622 2.2971245 2807.7423 21.916115 1982-1999 –> 2000-2009 4.510173 7578.7838
Atlantic (28.05,28.1] 2 cstar_tref ~ aou + phosphate 616 2.3228498 2794.4555 28.621772 2000-2009 –> 2010-2019 4.718844 5654.6199
Atlantic (28.05,28.1] 2 cstar_tref ~ aou + phosphate_star 616 2.8586132 3050.1470 33.641083 2000-2009 –> 2010-2019 5.877480 6197.7793
Atlantic (28.05,28.1] 2 cstar_tref ~ aou + silicate 616 3.8747005 3424.8373 36.020863 2000-2009 –> 2010-2019 8.196205 7018.7196
Atlantic (28.05,28.1] 2 cstar_tref ~ oxygen + phosphate 616 2.2693259 2765.7351 27.950305 2000-2009 –> 2010-2019 4.706199 5646.9450
Atlantic (28.05,28.1] 2 cstar_tref ~ oxygen + phosphate_star 616 2.2693259 2765.7351 27.950305 2000-2009 –> 2010-2019 4.706199 5646.9450
Atlantic (28.05,28.1] 2 cstar_tref ~ oxygen + silicate 616 3.8792442 3426.2812 36.121073 2000-2009 –> 2010-2019 8.218458 7025.2509
Atlantic (28.05,28.1] 2 cstar_tref ~ phosphate + phosphate_star 616 2.2693259 2765.7351 27.950305 2000-2009 –> 2010-2019 4.706199 5646.9450
Atlantic (28.05,28.1] 2 cstar_tref ~ sal + phosphate 616 1.8694828 2526.9476 19.340104 2000-2009 –> 2010-2019 4.107114 5302.0471
Atlantic (28.05,28.1] 2 cstar_tref ~ silicate + phosphate_star 616 3.6401408 3347.9038 35.794680 2000-2009 –> 2010-2019 7.490073 6798.0443
Atlantic (28.05,28.1] 2 cstar_tref ~ temp + phosphate 616 1.8959767 2544.2847 22.349264 2000-2009 –> 2010-2019 4.193101 5352.0270
Atlantic (28.1,28.15] 2 cstar_tref ~ aou + phosphate 732 1.8240938 2965.3120 25.423564 1982-1999 –> 2000-2009 3.204293 7482.9324
Atlantic (28.1,28.15] 2 cstar_tref ~ aou + phosphate_star 732 2.2744001 3288.3187 30.733835 1982-1999 –> 2000-2009 4.106151 8539.0240
Atlantic (28.1,28.15] 2 cstar_tref ~ aou + silicate 732 2.9746830 3681.2873 26.442024 1982-1999 –> 2000-2009 6.424631 10571.6882
Atlantic (28.1,28.15] 2 cstar_tref ~ oxygen + phosphate 732 1.8530281 2988.3521 25.870497 1982-1999 –> 2000-2009 3.347456 7711.9192
Atlantic (28.1,28.15] 2 cstar_tref ~ oxygen + phosphate_star 732 1.8530281 2988.3521 25.870497 1982-1999 –> 2000-2009 3.347456 7711.9192
Atlantic (28.1,28.15] 2 cstar_tref ~ oxygen + silicate 732 2.9516059 3669.8855 28.607684 1982-1999 –> 2000-2009 6.280124 10467.4814
Atlantic (28.1,28.15] 2 cstar_tref ~ sal + silicate 732 2.9409317 3664.5815 34.826386 1982-1999 –> 2000-2009 5.909276 10165.5637
Atlantic (28.1,28.15] 2 cstar_tref ~ silicate + phosphate_star 732 2.4184126 3378.2011 27.803745 1982-1999 –> 2000-2009 4.708434 9207.2253
Atlantic (28.1,28.15] 2 cstar_tref ~ temp + phosphate 732 1.7045759 2866.1011 22.052391 1982-1999 –> 2000-2009 3.322595 7795.4674
Atlantic (28.1,28.15] 2 cstar_tref ~ temp + silicate 732 2.9983926 3692.9098 32.079187 1982-1999 –> 2000-2009 6.192943 10384.1063
Atlantic (28.1,28.15] 2 cstar_tref ~ aou + phosphate 701 1.7913674 2814.6887 28.221794 2000-2009 –> 2010-2019 3.615461 5780.0007
Atlantic (28.1,28.15] 2 cstar_tref ~ aou + phosphate_star 701 2.1704697 3083.8227 33.161223 2000-2009 –> 2010-2019 4.444870 6372.1414
Atlantic (28.1,28.15] 2 cstar_tref ~ aou + silicate 701 2.8526982 3467.0197 29.478976 2000-2009 –> 2010-2019 5.827381 7148.3070
Atlantic (28.1,28.15] 2 cstar_tref ~ oxygen + phosphate 701 1.7430372 2776.3438 27.654958 2000-2009 –> 2010-2019 3.596065 5764.6959
Atlantic (28.1,28.15] 2 cstar_tref ~ oxygen + phosphate_star 701 1.7430372 2776.3438 27.654958 2000-2009 –> 2010-2019 3.596065 5764.6959
Atlantic (28.1,28.15] 2 cstar_tref ~ oxygen + silicate 701 2.8413641 3461.4383 30.557526 2000-2009 –> 2010-2019 5.792970 7131.3239
Atlantic (28.1,28.15] 2 cstar_tref ~ sal + silicate 701 2.7865426 3434.1237 35.063331 2000-2009 –> 2010-2019 5.727474 7098.7052
Atlantic (28.1,28.15] 2 cstar_tref ~ silicate + phosphate_star 701 2.4718887 3266.1373 30.083182 2000-2009 –> 2010-2019 4.890301 6644.3384
Atlantic (28.1,28.15] 2 cstar_tref ~ temp + phosphate 701 1.4666638 2534.3043 22.314498 2000-2009 –> 2010-2019 3.171240 5400.4054
Atlantic (28.1,28.15] 2 cstar_tref ~ temp + silicate 701 2.8443904 3462.9308 32.688640 2000-2009 –> 2010-2019 5.842783 7155.8406
Atlantic (28.15,28.2] 2 cstar_tref ~ aou + phosphate_star 1187 1.0017753 3380.7709 3.133811 1982-1999 –> 2000-2009 2.203598 10517.9282
Atlantic (28.15,28.2] 2 cstar_tref ~ oxygen + phosphate 1187 0.9302641 3204.9513 2.603778 1982-1999 –> 2000-2009 2.056323 10052.4725
Atlantic (28.15,28.2] 2 cstar_tref ~ oxygen + phosphate_star 1187 0.9302641 3204.9513 2.603778 1982-1999 –> 2000-2009 2.056323 10052.4725
Atlantic (28.15,28.2] 2 cstar_tref ~ sal + aou 1187 1.6214142 4523.9113 6.272678 1982-1999 –> 2000-2009 3.544207 13751.3568
Atlantic (28.15,28.2] 3 cstar_tref ~ sal + oxygen + silicate 1187 1.0661275 3530.5742 5.222547 1982-1999 –> 2000-2009 2.484369 11406.2224
Atlantic (28.15,28.2] 2 cstar_tref ~ sal + phosphate 1187 0.8692799 3043.9861 3.200506 1982-1999 –> 2000-2009 1.954114 9725.6104
Atlantic (28.15,28.2] 2 cstar_tref ~ sal + silicate 1187 1.1882531 3786.0377 5.351390 1982-1999 –> 2000-2009 2.655565 11810.9837
Atlantic (28.15,28.2] 2 cstar_tref ~ silicate + phosphate 1187 0.8685613 3042.0230 2.881704 1982-1999 –> 2000-2009 1.933874 9642.8784
Atlantic (28.15,28.2] 2 cstar_tref ~ silicate + phosphate_star 1187 1.6657623 4587.9716 4.906122 1982-1999 –> 2000-2009 3.324608 13158.6411
Atlantic (28.15,28.2] 2 cstar_tref ~ temp + phosphate 1187 0.9561402 3270.0845 2.763329 1982-1999 –> 2000-2009 2.192823 10534.4252
Atlantic (28.15,28.2] 2 cstar_tref ~ aou + phosphate 1172 1.3622391 4058.5920 24.555971 2000-2009 –> 2010-2019 2.282828 7238.7229
Atlantic (28.15,28.2] 2 cstar_tref ~ aou + phosphate_star 1172 1.6013554 4437.6652 30.460753 2000-2009 –> 2010-2019 2.603131 7818.4361
Atlantic (28.15,28.2] 2 cstar_tref ~ oxygen + phosphate 1172 1.3785623 4086.5124 25.282962 2000-2009 –> 2010-2019 2.308826 7291.4636
Atlantic (28.15,28.2] 2 cstar_tref ~ oxygen + phosphate_star 1172 1.3785623 4086.5124 25.282962 2000-2009 –> 2010-2019 2.308826 7291.4636
Atlantic (28.15,28.2] 3 cstar_tref ~ sal + oxygen + silicate 1172 1.9782366 4935.0824 40.533152 2000-2009 –> 2010-2019 3.044364 8465.6565
Atlantic (28.15,28.2] 2 cstar_tref ~ sal + phosphate 1172 1.3680237 4068.5246 23.861522 2000-2009 –> 2010-2019 2.237304 7112.5107
Atlantic (28.15,28.2] 2 cstar_tref ~ sal + silicate 1172 2.0274651 4990.6990 39.359255 2000-2009 –> 2010-2019 3.215718 8776.7367
Atlantic (28.15,28.2] 2 cstar_tref ~ silicate + phosphate 1172 1.3707847 4073.2505 25.490993 2000-2009 –> 2010-2019 2.239346 7115.2735
Atlantic (28.15,28.2] 2 cstar_tref ~ silicate + phosphate_star 1172 1.9615685 4913.2488 24.121146 2000-2009 –> 2010-2019 3.627331 9501.2204
Atlantic (28.15,28.2] 2 cstar_tref ~ temp + phosphate 1172 1.3336579 4008.8892 22.729273 2000-2009 –> 2010-2019 2.289798 7278.9736
Atlantic (28.2, Inf] 2 cstar_tref ~ aou + phosphate 2098 0.8645880 5351.3391 9.521783 1982-1999 –> 2000-2009 1.812834 15496.2875
Atlantic (28.2, Inf] 3 cstar_tref ~ aou + silicate + phosphate 2098 0.8639932 5350.4512 9.451563 1982-1999 –> 2000-2009 1.704751 14604.4526
Atlantic (28.2, Inf] 3 cstar_tref ~ aou + silicate + phosphate_star 2098 0.9960966 5947.4552 18.082434 1982-1999 –> 2000-2009 2.129644 17419.1821
Atlantic (28.2, Inf] 2 cstar_tref ~ oxygen + phosphate 2098 0.9543114 5765.6390 14.289005 1982-1999 –> 2000-2009 1.992798 16585.2887
Atlantic (28.2, Inf] 2 cstar_tref ~ oxygen + phosphate_star 2098 0.9543114 5765.6390 14.289005 1982-1999 –> 2000-2009 1.992798 16585.2887
Atlantic (28.2, Inf] 3 cstar_tref ~ oxygen + silicate + phosphate 2098 0.9343785 5679.0681 14.747462 1982-1999 –> 2000-2009 1.963599 16434.1979
Atlantic (28.2, Inf] 3 cstar_tref ~ oxygen + silicate + phosphate_star 2098 0.9343785 5679.0681 14.747462 1982-1999 –> 2000-2009 1.963599 16434.1979
Atlantic (28.2, Inf] 3 cstar_tref ~ sal + aou + phosphate_star 2098 0.9296263 5657.6730 16.943972 1982-1999 –> 2000-2009 2.016167 16815.0514
Atlantic (28.2, Inf] 3 cstar_tref ~ sal + oxygen + phosphate 2098 0.8883097 5466.9136 13.858533 1982-1999 –> 2000-2009 1.915134 16204.7427
Atlantic (28.2, Inf] 3 cstar_tref ~ sal + oxygen + phosphate_star 2098 0.8883097 5466.9136 13.858533 1982-1999 –> 2000-2009 1.915134 16204.7427
Atlantic (28.2, Inf] 2 cstar_tref ~ aou + phosphate 2152 1.1184930 6597.0829 18.135756 2000-2009 –> 2010-2019 1.983081 11948.4219
Atlantic (28.2, Inf] 3 cstar_tref ~ aou + silicate + phosphate 2152 1.1141786 6582.4489 18.349692 2000-2009 –> 2010-2019 1.978172 11932.9001
Atlantic (28.2, Inf] 3 cstar_tref ~ aou + silicate + phosphate_star 2152 1.2923586 7220.9535 26.196527 2000-2009 –> 2010-2019 2.288455 13168.4088
Atlantic (28.2, Inf] 2 cstar_tref ~ oxygen + phosphate 2152 1.2398386 7040.3906 22.189530 2000-2009 –> 2010-2019 2.194150 12806.0297
Atlantic (28.2, Inf] 2 cstar_tref ~ oxygen + phosphate_star 2152 1.2398386 7040.3906 22.189530 2000-2009 –> 2010-2019 2.194150 12806.0297
Atlantic (28.2, Inf] 3 cstar_tref ~ oxygen + silicate + phosphate 2152 1.2084603 6932.0613 22.901356 2000-2009 –> 2010-2019 2.142839 12611.1294
Atlantic (28.2, Inf] 3 cstar_tref ~ oxygen + silicate + phosphate_star 2152 1.2084603 6932.0613 22.901356 2000-2009 –> 2010-2019 2.142839 12611.1294
Atlantic (28.2, Inf] 3 cstar_tref ~ sal + aou + phosphate_star 2152 1.1447596 6698.9891 20.767967 2000-2009 –> 2010-2019 2.074386 12356.6621
Atlantic (28.2, Inf] 3 cstar_tref ~ sal + oxygen + phosphate 2152 1.0979191 6519.1769 18.791518 2000-2009 –> 2010-2019 1.986229 11986.0904
Atlantic (28.2, Inf] 3 cstar_tref ~ sal + oxygen + phosphate_star 2152 1.0979191 6519.1769 18.791518 2000-2009 –> 2010-2019 1.986229 11986.0904
Indo-Pacific (-Inf,26] 4 cstar_tref ~ sal + temp + aou + phosphate_star 4068 6.4428899 26713.6663 30.664457 1982-1999 –> 2000-2009 13.150399 74518.3485
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4068 6.3896745 26648.1875 30.157338 1982-1999 –> 2000-2009 13.029942 74309.9259
Indo-Pacific (-Inf,26] 4 cstar_tref ~ sal + temp + oxygen + phosphate_star 4068 6.4432668 26714.1421 30.718704 1982-1999 –> 2000-2009 13.155426 74528.7958
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 4068 6.3945847 26654.4372 30.238415 1982-1999 –> 2000-2009 13.044120 74336.2396
Indo-Pacific (-Inf,26] 4 cstar_tref ~ sal + temp + phosphate + phosphate_star 4068 6.4432668 26714.1421 30.718704 1982-1999 –> 2000-2009 13.155426 74528.7958
Indo-Pacific (-Inf,26] 3 cstar_tref ~ sal + temp + phosphate_star 4068 6.4436248 26712.5942 30.842472 1982-1999 –> 2000-2009 13.237756 74699.8708
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4068 6.3945847 26654.4372 30.238415 1982-1999 –> 2000-2009 13.044120 74336.2396
Indo-Pacific (-Inf,26] 4 cstar_tref ~ sal + temp + silicate + phosphate_star 4068 6.4287506 26695.7918 31.302451 1982-1999 –> 2000-2009 13.218529 74675.8494
Indo-Pacific (-Inf,26] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 4068 7.3926659 27832.4576 28.698337 1982-1999 –> 2000-2009 15.106969 77649.0008
Indo-Pacific (-Inf,26] 4 cstar_tref ~ temp + silicate + phosphate + phosphate_star 4068 7.3926659 27832.4576 28.698337 1982-1999 –> 2000-2009 15.106969 77649.0008
Indo-Pacific (-Inf,26] 4 cstar_tref ~ sal + temp + aou + phosphate_star 4036 6.2216439 26221.5636 25.833747 2000-2009 –> 2010-2019 12.664534 52935.2299
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4036 6.1462020 26125.0865 25.038689 2000-2009 –> 2010-2019 12.535877 52773.2740
Indo-Pacific (-Inf,26] 4 cstar_tref ~ sal + temp + oxygen + phosphate_star 4036 6.2182336 26217.1379 25.863706 2000-2009 –> 2010-2019 12.661500 52931.2800
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + temp + oxygen + silicate + phosphate_star 4036 6.1481960 26127.7050 25.113367 2000-2009 –> 2010-2019 12.542781 52782.1422
Indo-Pacific (-Inf,26] 4 cstar_tref ~ sal + temp + phosphate + phosphate_star 4036 6.2182336 26217.1379 25.863706 2000-2009 –> 2010-2019 12.661500 52931.2800
Indo-Pacific (-Inf,26] 3 cstar_tref ~ sal + temp + phosphate_star 4036 6.2617745 26271.4621 24.834075 2000-2009 –> 2010-2019 12.705399 52984.0563
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4036 6.1481960 26127.7050 25.113367 2000-2009 –> 2010-2019 12.542781 52782.1422
Indo-Pacific (-Inf,26] 4 cstar_tref ~ sal + temp + silicate + phosphate_star 4036 6.1533062 26132.4114 25.491924 2000-2009 –> 2010-2019 12.582057 52828.2031
Indo-Pacific (-Inf,26] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 4036 7.2224238 27425.5544 24.309044 2000-2009 –> 2010-2019 14.615090 55258.0120
Indo-Pacific (-Inf,26] 4 cstar_tref ~ temp + silicate + phosphate + phosphate_star 4036 7.2224238 27425.5544 24.309044 2000-2009 –> 2010-2019 14.615090 55258.0120
Indo-Pacific (26,26.5] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 4229 6.7509679 28165.5053 56.004693 1982-1999 –> 2000-2009 13.233508 78892.4035
Indo-Pacific (26,26.5] 3 cstar_tref ~ sal + oxygen + phosphate_star 4229 7.3152743 28842.5021 62.076410 1982-1999 –> 2000-2009 14.370359 80872.8289
Indo-Pacific (26,26.5] 4 cstar_tref ~ sal + oxygen + silicate + phosphate_star 4229 6.8340386 28268.9460 55.375743 1982-1999 –> 2000-2009 13.462074 79338.1954
Indo-Pacific (26,26.5] 4 cstar_tref ~ sal + silicate + phosphate + phosphate_star 4229 6.8340386 28268.9460 55.375743 1982-1999 –> 2000-2009 13.462074 79338.1954
Indo-Pacific (26,26.5] 3 cstar_tref ~ sal + silicate + phosphate_star 4229 6.9175803 28369.7126 50.053515 1982-1999 –> 2000-2009 13.865365 80163.6554
Indo-Pacific (26,26.5] 3 cstar_tref ~ temp + aou + phosphate_star 4229 5.4186790 26304.1508 53.395037 1982-1999 –> 2000-2009 10.632185 73668.7779
Indo-Pacific (26,26.5] 3 cstar_tref ~ temp + aou + silicate 4229 7.1484792 28647.4192 27.250623 1982-1999 –> 2000-2009 13.910649 80023.6930
Indo-Pacific (26,26.5] 3 cstar_tref ~ temp + oxygen + phosphate_star 4229 5.4190257 26304.6920 53.509273 1982-1999 –> 2000-2009 10.631787 73667.1159
Indo-Pacific (26,26.5] 3 cstar_tref ~ temp + oxygen + silicate 4229 7.1571040 28657.6178 27.615755 1982-1999 –> 2000-2009 13.929659 80057.5622
Indo-Pacific (26,26.5] 2 cstar_tref ~ temp + phosphate_star 4229 5.4280396 26316.7492 51.683669 1982-1999 –> 2000-2009 11.017555 74753.5037
Indo-Pacific (26,26.5] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 4214 6.7947203 28120.0914 27.054528 2000-2009 –> 2010-2019 13.545688 56285.5968
Indo-Pacific (26,26.5] 3 cstar_tref ~ sal + oxygen + phosphate_star 4214 7.3921771 28828.3731 56.706434 2000-2009 –> 2010-2019 14.707451 57670.8752
Indo-Pacific (26,26.5] 4 cstar_tref ~ sal + oxygen + silicate + phosphate_star 4214 6.8392094 28175.0947 29.698507 2000-2009 –> 2010-2019 13.673248 56444.0408
Indo-Pacific (26,26.5] 4 cstar_tref ~ sal + silicate + phosphate + phosphate_star 4214 6.8392094 28175.0947 29.698507 2000-2009 –> 2010-2019 13.673248 56444.0408
Indo-Pacific (26,26.5] 3 cstar_tref ~ sal + silicate + phosphate_star 4214 6.8571798 28195.2107 33.967019 2000-2009 –> 2010-2019 13.774760 56564.9233
Indo-Pacific (26,26.5] 3 cstar_tref ~ temp + aou + phosphate_star 4214 5.5537623 26418.4142 41.144616 2000-2009 –> 2010-2019 10.972441 52722.5651
Indo-Pacific (26,26.5] 3 cstar_tref ~ temp + aou + silicate 4214 7.2224986 28632.6637 29.412015 2000-2009 –> 2010-2019 14.370978 57280.0829
Indo-Pacific (26,26.5] 3 cstar_tref ~ temp + oxygen + phosphate_star 4214 5.5536243 26418.2048 41.172987 2000-2009 –> 2010-2019 10.972650 52722.8968
Indo-Pacific (26,26.5] 3 cstar_tref ~ temp + oxygen + silicate 4214 7.2276610 28638.6855 29.498648 2000-2009 –> 2010-2019 14.384765 57296.3034
Indo-Pacific (26,26.5] 2 cstar_tref ~ temp + phosphate_star 4214 5.5976571 26482.7641 42.795836 2000-2009 –> 2010-2019 11.025697 52799.5133
Indo-Pacific (26.5,26.75] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 3529 4.8657521 21194.1861 38.819399 1982-1999 –> 2000-2009 9.756065 59012.0931
Indo-Pacific (26.5,26.75] 3 cstar_tref ~ sal + oxygen + phosphate_star 3529 4.9451103 21306.3703 37.975904 1982-1999 –> 2000-2009 9.883167 59244.4590
Indo-Pacific (26.5,26.75] 4 cstar_tref ~ sal + oxygen + silicate + phosphate_star 3529 4.8048445 21105.2792 36.621196 1982-1999 –> 2000-2009 9.675707 58873.0657
Indo-Pacific (26.5,26.75] 4 cstar_tref ~ sal + silicate + phosphate + phosphate_star 3529 4.8048445 21105.2792 36.621196 1982-1999 –> 2000-2009 9.675707 58873.0657
Indo-Pacific (26.5,26.75] 3 cstar_tref ~ temp + aou + phosphate_star 3529 4.4351452 20538.1851 36.125914 1982-1999 –> 2000-2009 8.704023 56644.9451
Indo-Pacific (26.5,26.75] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 3529 4.1687868 20103.0458 32.793647 1982-1999 –> 2000-2009 8.262783 55685.7600
Indo-Pacific (26.5,26.75] 3 cstar_tref ~ temp + oxygen + phosphate_star 3529 4.4243200 20520.9371 35.061218 1982-1999 –> 2000-2009 8.690031 56618.3659
Indo-Pacific (26.5,26.75] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 3529 4.1645262 20095.8287 32.362053 1982-1999 –> 2000-2009 8.261433 55687.4805
Indo-Pacific (26.5,26.75] 2 cstar_tref ~ temp + phosphate_star 3529 5.0639716 21472.0104 53.917048 1982-1999 –> 2000-2009 9.388746 57740.3726
Indo-Pacific (26.5,26.75] 3 cstar_tref ~ temp + silicate + phosphate_star 3529 4.2073011 20165.9533 34.710459 1982-1999 –> 2000-2009 8.323312 55814.1108
Indo-Pacific (26.5,26.75] 3 cstar_tref ~ sal + aou + phosphate_star 3462 5.4247240 21542.9861 24.342689 2000-2009 –> 2010-2019 10.461918 42979.5750
Indo-Pacific (26.5,26.75] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 3462 5.2225762 21282.0388 21.358319 2000-2009 –> 2010-2019 10.088328 42476.2248
Indo-Pacific (26.5,26.75] 3 cstar_tref ~ sal + oxygen + phosphate_star 3462 5.3157512 21402.4795 21.894265 2000-2009 –> 2010-2019 10.260862 42708.8497
Indo-Pacific (26.5,26.75] 4 cstar_tref ~ sal + oxygen + silicate + phosphate_star 3462 5.1459398 21179.6827 20.330712 2000-2009 –> 2010-2019 9.950784 42284.9618
Indo-Pacific (26.5,26.75] 3 cstar_tref ~ sal + silicate + phosphate_star 3462 5.5153247 21657.6718 24.512131 2000-2009 –> 2010-2019 10.556092 43099.2670
Indo-Pacific (26.5,26.75] 3 cstar_tref ~ temp + aou + phosphate_star 3462 4.7794590 20666.1330 19.771771 2000-2009 –> 2010-2019 9.214604 41204.3181
Indo-Pacific (26.5,26.75] 3 cstar_tref ~ temp + oxygen + phosphate_star 3462 4.7672783 20648.4644 19.641583 2000-2009 –> 2010-2019 9.191598 41169.4015
Indo-Pacific (26.5,26.75] 2 cstar_tref ~ temp + phosphate_star 3462 5.6919927 21873.9846 39.173133 2000-2009 –> 2010-2019 10.755964 43345.9950
Indo-Pacific (26.5,26.75] 3 cstar_tref ~ temp + silicate + phosphate 3462 5.6044926 21768.7188 31.462225 2000-2009 –> 2010-2019 11.011486 43705.3262
Indo-Pacific (26.5,26.75] 3 cstar_tref ~ temp + silicate + phosphate_star 3462 4.5787338 20369.0598 28.907657 2000-2009 –> 2010-2019 8.786035 40535.0131
Indo-Pacific (26.75,27] 3 cstar_tref ~ sal + aou + phosphate_star 4017 5.7683009 25488.3535 31.857505 1982-1999 –> 2000-2009 11.147714 72213.4161
Indo-Pacific (26.75,27] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 4017 5.7070008 25404.5188 31.315452 1982-1999 –> 2000-2009 11.037503 71994.0054
Indo-Pacific (26.75,27] 3 cstar_tref ~ sal + oxygen + phosphate_star 4017 5.5029879 25110.0616 28.497004 1982-1999 –> 2000-2009 10.692238 71293.0399
Indo-Pacific (26.75,27] 4 cstar_tref ~ sal + oxygen + silicate + phosphate_star 4017 5.3996200 24959.7159 27.366689 1982-1999 –> 2000-2009 10.491255 70858.6640
Indo-Pacific (26.75,27] 3 cstar_tref ~ sal + phosphate + phosphate_star 4017 5.5029879 25110.0616 28.497004 1982-1999 –> 2000-2009 10.692238 71293.0399
Indo-Pacific (26.75,27] 4 cstar_tref ~ sal + silicate + phosphate + phosphate_star 4017 5.3996200 24959.7159 27.366689 1982-1999 –> 2000-2009 10.491255 70858.6640
Indo-Pacific (26.75,27] 3 cstar_tref ~ temp + aou + phosphate_star 4017 4.8614849 24114.2693 25.034326 1982-1999 –> 2000-2009 9.242701 67747.8100
Indo-Pacific (26.75,27] 3 cstar_tref ~ temp + oxygen + phosphate_star 4017 4.8057173 24021.5761 23.807994 1982-1999 –> 2000-2009 9.151079 67531.3479
Indo-Pacific (26.75,27] 3 cstar_tref ~ temp + phosphate + phosphate_star 4017 4.8057173 24021.5761 23.807994 1982-1999 –> 2000-2009 9.151079 67531.3479
Indo-Pacific (26.75,27] 3 cstar_tref ~ temp + silicate + phosphate_star 4017 6.0562148 25879.6689 35.914239 1982-1999 –> 2000-2009 10.917281 71078.6247
Indo-Pacific (26.75,27] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 4205 6.1237122 27185.6099 25.587650 2000-2009 –> 2010-2019 11.830713 52590.1288
Indo-Pacific (26.75,27] 3 cstar_tref ~ sal + oxygen + phosphate_star 4205 5.8315300 26772.4528 22.523379 2000-2009 –> 2010-2019 11.334518 51882.5145
Indo-Pacific (26.75,27] 4 cstar_tref ~ sal + oxygen + silicate + phosphate_star 4205 5.7543590 26662.4169 23.260730 2000-2009 –> 2010-2019 11.153979 51622.1329
Indo-Pacific (26.75,27] 3 cstar_tref ~ sal + phosphate + phosphate_star 4205 5.8315300 26772.4528 22.523379 2000-2009 –> 2010-2019 11.334518 51882.5145
Indo-Pacific (26.75,27] 3 cstar_tref ~ sal + silicate + phosphate 4205 6.0238256 27045.2996 35.681466 2000-2009 –> 2010-2019 11.873172 52645.7459
Indo-Pacific (26.75,27] 4 cstar_tref ~ sal + silicate + phosphate + phosphate_star 4205 5.7543590 26662.4169 23.260730 2000-2009 –> 2010-2019 11.153979 51622.1329
Indo-Pacific (26.75,27] 3 cstar_tref ~ temp + aou + phosphate_star 4205 5.2337170 25862.8469 19.460267 2000-2009 –> 2010-2019 10.095202 49977.1162
Indo-Pacific (26.75,27] 3 cstar_tref ~ temp + oxygen + phosphate_star 4205 5.1569438 25738.5672 18.103335 2000-2009 –> 2010-2019 9.962661 49760.1432
Indo-Pacific (26.75,27] 2 cstar_tref ~ temp + phosphate 4205 5.8644788 26817.8365 31.713508 2000-2009 –> 2010-2019 11.583151 52234.7687
Indo-Pacific (26.75,27] 3 cstar_tref ~ temp + phosphate + phosphate_star 4205 5.1569438 25738.5672 18.103335 2000-2009 –> 2010-2019 9.962661 49760.1432
Indo-Pacific (27,27.25] 3 cstar_tref ~ sal + aou + phosphate_star 5080 6.0912031 32783.9670 29.870223 1982-1999 –> 2000-2009 11.480349 87642.1408
Indo-Pacific (27,27.25] 3 cstar_tref ~ sal + oxygen + phosphate_star 5080 5.9238060 32500.8437 30.560094 1982-1999 –> 2000-2009 11.211122 87021.8644
Indo-Pacific (27,27.25] 2 cstar_tref ~ sal + phosphate 5080 6.9656220 34144.8425 28.411626 1982-1999 –> 2000-2009 13.653356 92816.6245
Indo-Pacific (27,27.25] 3 cstar_tref ~ sal + phosphate + phosphate_star 5080 5.9238060 32500.8437 30.560094 1982-1999 –> 2000-2009 11.211122 87021.8644
Indo-Pacific (27,27.25] 3 cstar_tref ~ sal + silicate + phosphate_star 5080 7.5903897 35019.5461 46.866258 1982-1999 –> 2000-2009 13.690135 92066.8247
Indo-Pacific (27,27.25] 3 cstar_tref ~ temp + aou + phosphate_star 5080 5.0691973 30917.9495 53.169649 1982-1999 –> 2000-2009 9.302169 81508.1578
Indo-Pacific (27,27.25] 3 cstar_tref ~ temp + oxygen + phosphate_star 5080 5.0411319 30861.5428 52.860320 1982-1999 –> 2000-2009 9.261587 81399.4146
Indo-Pacific (27,27.25] 2 cstar_tref ~ temp + phosphate 5080 6.9537337 34127.4875 27.432742 1982-1999 –> 2000-2009 13.633245 92777.5258
Indo-Pacific (27,27.25] 2 cstar_tref ~ temp + phosphate_star 5080 7.1763435 34447.6418 76.820216 1982-1999 –> 2000-2009 12.038820 87486.2384
Indo-Pacific (27,27.25] 3 cstar_tref ~ temp + silicate + phosphate_star 5080 5.9566829 32557.0753 71.490186 1982-1999 –> 2000-2009 10.516417 84461.5228
Indo-Pacific (27,27.25] 2 cstar_tref ~ oxygen + phosphate_star 4856 7.5920922 33475.9963 26.771087 2000-2009 –> 2010-2019 15.064583 68334.4926
Indo-Pacific (27,27.25] 2 cstar_tref ~ phosphate + phosphate_star 4856 7.5920922 33475.9963 26.771087 2000-2009 –> 2010-2019 15.064583 68334.4926
Indo-Pacific (27,27.25] 3 cstar_tref ~ sal + aou + phosphate_star 4856 6.5770227 32084.0810 26.101046 2000-2009 –> 2010-2019 12.668226 64868.0480
Indo-Pacific (27,27.25] 3 cstar_tref ~ sal + oxygen + phosphate_star 4856 6.3921027 31807.1056 28.220610 2000-2009 –> 2010-2019 12.315909 64307.9493
Indo-Pacific (27,27.25] 2 cstar_tref ~ sal + phosphate 4856 7.0144940 32707.4990 27.753521 2000-2009 –> 2010-2019 13.980116 66852.3415
Indo-Pacific (27,27.25] 3 cstar_tref ~ sal + phosphate + phosphate_star 4856 6.3921027 31807.1056 28.220610 2000-2009 –> 2010-2019 12.315909 64307.9493
Indo-Pacific (27,27.25] 3 cstar_tref ~ temp + aou + phosphate_star 4856 5.9172035 31057.3458 45.800024 2000-2009 –> 2010-2019 10.986401 61975.2953
Indo-Pacific (27,27.25] 3 cstar_tref ~ temp + oxygen + phosphate_star 4856 5.8765262 30990.3510 45.560328 2000-2009 –> 2010-2019 10.917658 61851.8938
Indo-Pacific (27,27.25] 2 cstar_tref ~ temp + phosphate 4856 6.9934646 32678.3388 26.556570 2000-2009 –> 2010-2019 13.947198 66805.8263
Indo-Pacific (27,27.25] 3 cstar_tref ~ temp + silicate + phosphate_star 4856 7.5409022 33412.2909 67.249167 2000-2009 –> 2010-2019 13.497585 65969.3662
Indo-Pacific (27.25,27.5] 2 cstar_tref ~ oxygen + phosphate_star 4470 5.2908243 27587.1186 25.913363 1982-1999 –> 2000-2009 9.933265 74223.8781
Indo-Pacific (27.25,27.5] 3 cstar_tref ~ sal + aou + phosphate_star 4470 4.3921456 25924.8821 26.927543 1982-1999 –> 2000-2009 8.137426 69174.1258
Indo-Pacific (27.25,27.5] 3 cstar_tref ~ sal + oxygen + phosphate_star 4470 4.3740975 25888.0705 28.349046 1982-1999 –> 2000-2009 8.112400 69107.8792
Indo-Pacific (27.25,27.5] 3 cstar_tref ~ sal + phosphate + phosphate_star 4470 4.3740975 25888.0705 28.349046 1982-1999 –> 2000-2009 8.112400 69107.8792
Indo-Pacific (27.25,27.5] 3 cstar_tref ~ sal + silicate + phosphate_star 4470 5.1147838 27286.5986 35.856819 1982-1999 –> 2000-2009 9.081200 71441.3143
Indo-Pacific (27.25,27.5] 3 cstar_tref ~ sal + temp + phosphate 4470 4.9746474 27038.2396 24.882180 1982-1999 –> 2000-2009 9.610546 73654.7378
Indo-Pacific (27.25,27.5] 3 cstar_tref ~ temp + aou + phosphate_star 4470 4.6648918 26463.4884 36.663164 1982-1999 –> 2000-2009 8.609361 70530.6228
Indo-Pacific (27.25,27.5] 3 cstar_tref ~ temp + oxygen + phosphate_star 4470 4.6578100 26449.9062 36.707105 1982-1999 –> 2000-2009 8.600625 70510.4226
Indo-Pacific (27.25,27.5] 3 cstar_tref ~ temp + phosphate + phosphate_star 4470 4.6578100 26449.9062 36.707105 1982-1999 –> 2000-2009 8.600625 70510.4226
Indo-Pacific (27.25,27.5] 3 cstar_tref ~ temp + silicate + phosphate_star 4470 5.2234732 27474.5836 49.887875 1982-1999 –> 2000-2009 9.328952 72173.2040
Indo-Pacific (27.25,27.5] 3 cstar_tref ~ sal + aou + phosphate_star 4504 4.8818259 27074.1562 31.343149 2000-2009 –> 2010-2019 9.273972 52999.0383
Indo-Pacific (27.25,27.5] 3 cstar_tref ~ sal + oxygen + phosphate_star 4504 4.8606423 27034.9829 33.492208 2000-2009 –> 2010-2019 9.234740 52923.0534
Indo-Pacific (27.25,27.5] 3 cstar_tref ~ sal + phosphate + phosphate_star 4504 4.8606423 27034.9829 33.492208 2000-2009 –> 2010-2019 9.234740 52923.0534
Indo-Pacific (27.25,27.5] 3 cstar_tref ~ sal + silicate + phosphate_star 4504 5.9175026 28807.2521 44.539769 2000-2009 –> 2010-2019 11.032286 56093.8508
Indo-Pacific (27.25,27.5] 3 cstar_tref ~ sal + temp + aou 4504 5.7367152 28527.7544 25.504859 2000-2009 –> 2010-2019 11.077631 56201.1158
Indo-Pacific (27.25,27.5] 3 cstar_tref ~ sal + temp + oxygen 4504 5.7181736 28498.5926 25.472837 2000-2009 –> 2010-2019 11.045416 56149.0368
Indo-Pacific (27.25,27.5] 3 cstar_tref ~ sal + temp + phosphate 4504 5.3066825 27825.8519 23.785914 2000-2009 –> 2010-2019 10.281330 54864.0915
Indo-Pacific (27.25,27.5] 3 cstar_tref ~ temp + aou + phosphate_star 4504 5.2415766 27714.6523 42.584680 2000-2009 –> 2010-2019 9.906469 54178.1407
Indo-Pacific (27.25,27.5] 3 cstar_tref ~ temp + oxygen + phosphate_star 4504 5.2304463 27695.5037 42.635810 2000-2009 –> 2010-2019 9.888256 54145.4099
Indo-Pacific (27.25,27.5] 3 cstar_tref ~ temp + phosphate + phosphate_star 4504 5.2304463 27695.5037 42.635810 2000-2009 –> 2010-2019 9.888256 54145.4099
Indo-Pacific (27.5,27.75] 2 cstar_tref ~ oxygen + phosphate_star 4498 3.8600945 24923.5933 14.236840 1982-1999 –> 2000-2009 6.845528 65662.1914
Indo-Pacific (27.5,27.75] 3 cstar_tref ~ sal + aou + phosphate_star 4498 3.4729481 23974.8256 26.741037 1982-1999 –> 2000-2009 6.102738 62659.3309
Indo-Pacific (27.5,27.75] 3 cstar_tref ~ sal + oxygen + phosphate_star 4498 3.4786369 23989.5491 27.177760 1982-1999 –> 2000-2009 6.110164 62684.7618
Indo-Pacific (27.5,27.75] 3 cstar_tref ~ sal + phosphate + phosphate_star 4498 3.4786369 23989.5491 27.177760 1982-1999 –> 2000-2009 6.110164 62684.7618
Indo-Pacific (27.5,27.75] 2 cstar_tref ~ sal + phosphate_star 4498 3.6949373 24530.2156 31.389740 1982-1999 –> 2000-2009 6.665128 65185.8389
Indo-Pacific (27.5,27.75] 3 cstar_tref ~ sal + silicate + phosphate_star 4498 3.5549531 24184.7744 28.587146 1982-1999 –> 2000-2009 6.066395 62122.8616
Indo-Pacific (27.5,27.75] 3 cstar_tref ~ sal + temp + aou 4498 3.4653283 23955.0661 13.445856 1982-1999 –> 2000-2009 6.315683 63945.1174
Indo-Pacific (27.5,27.75] 3 cstar_tref ~ sal + temp + oxygen 4498 3.4646625 23953.3377 13.481193 1982-1999 –> 2000-2009 6.315398 63945.5542
Indo-Pacific (27.5,27.75] 3 cstar_tref ~ sal + temp + phosphate 4498 3.4455761 23903.6427 15.066373 1982-1999 –> 2000-2009 6.330352 64088.2738
Indo-Pacific (27.5,27.75] 3 cstar_tref ~ sal + temp + silicate 4498 3.5824528 24254.0963 16.105669 1982-1999 –> 2000-2009 6.341951 63719.0292
Indo-Pacific (27.5,27.75] 3 cstar_tref ~ sal + aou + phosphate_star 4442 4.1705130 25302.5488 29.528612 2000-2009 –> 2010-2019 7.643461 49277.3744
Indo-Pacific (27.5,27.75] 3 cstar_tref ~ sal + oxygen + phosphate_star 4442 4.1976050 25360.0734 30.626689 2000-2009 –> 2010-2019 7.676242 49349.6224
Indo-Pacific (27.5,27.75] 2 cstar_tref ~ sal + phosphate 4442 4.4800091 25936.5191 19.403209 2000-2009 –> 2010-2019 8.342632 50866.0021
Indo-Pacific (27.5,27.75] 3 cstar_tref ~ sal + phosphate + phosphate_star 4442 4.1976050 25360.0734 30.626689 2000-2009 –> 2010-2019 7.676242 49349.6224
Indo-Pacific (27.5,27.75] 3 cstar_tref ~ sal + silicate + phosphate_star 4442 4.3919294 25762.1145 33.399371 2000-2009 –> 2010-2019 7.946883 49946.8889
Indo-Pacific (27.5,27.75] 3 cstar_tref ~ sal + temp + aou 4442 3.9815894 24890.7046 15.285111 2000-2009 –> 2010-2019 7.446918 48845.7707
Indo-Pacific (27.5,27.75] 3 cstar_tref ~ sal + temp + oxygen 4442 3.9797779 24886.6619 15.373761 2000-2009 –> 2010-2019 7.444440 48839.9996
Indo-Pacific (27.5,27.75] 3 cstar_tref ~ sal + temp + phosphate 4442 3.9801977 24887.5989 19.273494 2000-2009 –> 2010-2019 7.425774 48791.2416
Indo-Pacific (27.5,27.75] 3 cstar_tref ~ sal + temp + silicate 4442 4.2822343 25537.4051 21.421952 2000-2009 –> 2010-2019 7.864687 49791.5014
Indo-Pacific (27.5,27.75] 3 cstar_tref ~ temp + aou + phosphate_star 4442 4.5789278 26132.5439 16.692661 2000-2009 –> 2010-2019 8.435725 51050.4483
Indo-Pacific (27.75,27.85] 2 cstar_tref ~ aou + phosphate_star 1730 2.7833134 8459.3289 18.216894 1982-1999 –> 2000-2009 5.037978 22443.7356
Indo-Pacific (27.75,27.85] 2 cstar_tref ~ sal + aou 1730 2.9058251 8608.3694 28.233312 1982-1999 –> 2000-2009 5.011714 22165.3098
Indo-Pacific (27.75,27.85] 2 cstar_tref ~ sal + oxygen 1730 2.9146786 8618.8954 28.188923 1982-1999 –> 2000-2009 5.036726 22223.6992
Indo-Pacific (27.75,27.85] 2 cstar_tref ~ sal + phosphate 1730 2.8722695 8568.1820 27.191857 1982-1999 –> 2000-2009 5.052721 22343.0033
Indo-Pacific (27.75,27.85] 2 cstar_tref ~ sal + silicate 1730 2.8403034 8529.4589 26.827201 1982-1999 –> 2000-2009 4.947375 22089.9153
Indo-Pacific (27.75,27.85] 2 cstar_tref ~ silicate + phosphate_star 1730 2.7598198 8429.9996 18.374027 1982-1999 –> 2000-2009 4.956597 22251.5330
Indo-Pacific (27.75,27.85] 2 cstar_tref ~ temp + aou 1730 2.8427968 8532.4951 27.115254 1982-1999 –> 2000-2009 4.929520 22032.1826
Indo-Pacific (27.75,27.85] 2 cstar_tref ~ temp + oxygen 1730 2.8433639 8533.1852 27.097327 1982-1999 –> 2000-2009 4.931170 22036.1230
Indo-Pacific (27.75,27.85] 2 cstar_tref ~ temp + phosphate 1730 2.8168275 8500.7422 26.000871 1982-1999 –> 2000-2009 4.987627 22247.7837
Indo-Pacific (27.75,27.85] 2 cstar_tref ~ temp + silicate 1730 2.7930587 8471.4224 25.560890 1982-1999 –> 2000-2009 4.889373 21999.8265
Indo-Pacific (27.75,27.85] 3 cstar_tref ~ sal + aou + phosphate_star 1851 3.0349054 9372.7977 28.870034 2000-2009 –> 2010-2019 5.805542 17818.3324
Indo-Pacific (27.75,27.85] 3 cstar_tref ~ sal + oxygen + phosphate_star 1851 3.0360751 9374.2243 28.721401 2000-2009 –> 2010-2019 5.808585 17822.0975
Indo-Pacific (27.75,27.85] 2 cstar_tref ~ sal + phosphate 1851 3.0628863 9404.7726 33.110125 2000-2009 –> 2010-2019 5.935156 17972.9546
Indo-Pacific (27.75,27.85] 3 cstar_tref ~ sal + phosphate + phosphate_star 1851 3.0360751 9374.2243 28.721401 2000-2009 –> 2010-2019 5.808585 17822.0975
Indo-Pacific (27.75,27.85] 2 cstar_tref ~ sal + silicate 1851 3.0586660 9399.6682 32.840985 2000-2009 –> 2010-2019 5.898969 17929.1272
Indo-Pacific (27.75,27.85] 3 cstar_tref ~ sal + silicate + phosphate_star 1851 3.0305278 9367.4541 28.596015 2000-2009 –> 2010-2019 5.771796 17776.1173
Indo-Pacific (27.75,27.85] 2 cstar_tref ~ temp + aou 1851 3.0847930 9431.1563 31.242089 2000-2009 –> 2010-2019 5.927590 17963.6514
Indo-Pacific (27.75,27.85] 2 cstar_tref ~ temp + oxygen 1851 3.0853105 9431.7773 31.225697 2000-2009 –> 2010-2019 5.928674 17964.9625
Indo-Pacific (27.75,27.85] 2 cstar_tref ~ temp + phosphate 1851 3.0815292 9427.2374 30.681689 2000-2009 –> 2010-2019 5.898357 17927.9796
Indo-Pacific (27.75,27.85] 2 cstar_tref ~ temp + silicate 1851 3.0792238 9424.4667 30.437981 2000-2009 –> 2010-2019 5.872282 17895.8892
Indo-Pacific (27.85,27.95] 3 cstar_tref ~ sal + aou + phosphate_star 2180 2.9503635 10913.7798 29.732656 1982-1999 –> 2000-2009 5.216618 28865.0581
Indo-Pacific (27.85,27.95] 3 cstar_tref ~ sal + silicate + phosphate_star 2180 2.9439356 10904.2703 29.587246 1982-1999 –> 2000-2009 5.183325 28759.9082
Indo-Pacific (27.85,27.95] 3 cstar_tref ~ sal + temp + aou 2180 2.9821661 10960.5256 34.654816 1982-1999 –> 2000-2009 5.163379 28605.0594
Indo-Pacific (27.85,27.95] 3 cstar_tref ~ sal + temp + oxygen 2180 2.9823458 10960.7884 34.634560 1982-1999 –> 2000-2009 5.163456 28604.9441
Indo-Pacific (27.85,27.95] 3 cstar_tref ~ sal + temp + phosphate 2180 2.9542402 10919.5049 33.520031 1982-1999 –> 2000-2009 5.180011 28726.2203
Indo-Pacific (27.85,27.95] 3 cstar_tref ~ sal + temp + silicate 2180 2.9686177 10940.6724 33.449594 1982-1999 –> 2000-2009 5.180063 28695.6014
Indo-Pacific (27.85,27.95] 2 cstar_tref ~ temp + aou 2180 2.9859093 10963.9948 35.204940 1982-1999 –> 2000-2009 5.167263 28607.0479
Indo-Pacific (27.85,27.95] 2 cstar_tref ~ temp + oxygen 2180 2.9859838 10964.1036 35.176852 1982-1999 –> 2000-2009 5.167271 28606.9101
Indo-Pacific (27.85,27.95] 2 cstar_tref ~ temp + phosphate 2180 2.9567577 10921.2188 33.947763 1982-1999 –> 2000-2009 5.182572 28726.0902
Indo-Pacific (27.85,27.95] 2 cstar_tref ~ temp + silicate 2180 2.9687794 10938.9098 33.560272 1982-1999 –> 2000-2009 5.186482 28714.5016
Indo-Pacific (27.85,27.95] 3 cstar_tref ~ sal + aou + phosphate_star 2269 3.0036979 11440.2359 35.787741 2000-2009 –> 2010-2019 5.954061 22354.0157
Indo-Pacific (27.85,27.95] 3 cstar_tref ~ sal + oxygen + phosphate_star 2269 3.0078536 11446.5101 35.548862 2000-2009 –> 2010-2019 5.962644 22366.8267
Indo-Pacific (27.85,27.95] 3 cstar_tref ~ sal + phosphate + phosphate_star 2269 3.0078536 11446.5101 35.548862 2000-2009 –> 2010-2019 5.962644 22366.8267
Indo-Pacific (27.85,27.95] 3 cstar_tref ~ sal + silicate + phosphate_star 2269 3.0122317 11453.1106 35.339140 2000-2009 –> 2010-2019 5.956167 22357.3808
Indo-Pacific (27.85,27.95] 3 cstar_tref ~ sal + temp + aou 2269 2.9963690 11429.1498 40.208596 2000-2009 –> 2010-2019 5.978535 22389.6754
Indo-Pacific (27.85,27.95] 3 cstar_tref ~ sal + temp + oxygen 2269 2.9973557 11430.6440 40.184475 2000-2009 –> 2010-2019 5.979702 22391.4323
Indo-Pacific (27.85,27.95] 3 cstar_tref ~ sal + temp + phosphate 2269 2.9853455 11412.4239 39.167624 2000-2009 –> 2010-2019 5.939586 22331.9287
Indo-Pacific (27.85,27.95] 3 cstar_tref ~ sal + temp + silicate 2269 3.0126933 11453.8059 38.935369 2000-2009 –> 2010-2019 5.981311 22394.4782
Indo-Pacific (27.85,27.95] 2 cstar_tref ~ temp + aou 2269 3.0133310 11452.7664 39.004327 2000-2009 –> 2010-2019 5.999240 22416.7612
Indo-Pacific (27.85,27.95] 2 cstar_tref ~ temp + phosphate 2269 3.0032385 11437.5418 38.009015 2000-2009 –> 2010-2019 5.959996 22358.7607
Indo-Pacific (27.95,28.05] 3 cstar_tref ~ sal + aou + phosphate_star 1874 2.3579336 8543.1623 8.210768 1982-1999 –> 2000-2009 4.455567 23857.1391
Indo-Pacific (27.95,28.05] 3 cstar_tref ~ sal + oxygen + phosphate_star 1874 2.3831694 8583.0621 8.268069 1982-1999 –> 2000-2009 4.512227 24002.4040
Indo-Pacific (27.95,28.05] 3 cstar_tref ~ sal + phosphate + phosphate_star 1874 2.3831694 8583.0621 8.268069 1982-1999 –> 2000-2009 4.512227 24002.4040
Indo-Pacific (27.95,28.05] 2 cstar_tref ~ sal + silicate 1874 2.3894780 8590.9704 9.517444 1982-1999 –> 2000-2009 4.433305 23718.8110
Indo-Pacific (27.95,28.05] 3 cstar_tref ~ sal + silicate + phosphate_star 1874 2.3430278 8519.3939 8.927922 1982-1999 –> 2000-2009 4.385863 23645.7963
Indo-Pacific (27.95,28.05] 2 cstar_tref ~ temp + aou 1874 2.2789511 8413.4665 7.923347 1982-1999 –> 2000-2009 4.309202 23494.0826
Indo-Pacific (27.95,28.05] 3 cstar_tref ~ temp + aou + phosphate_star 1874 2.2786179 8414.9185 7.915353 1982-1999 –> 2000-2009 4.292944 23441.7347
Indo-Pacific (27.95,28.05] 2 cstar_tref ~ temp + oxygen 1874 2.2799806 8415.1594 7.933030 1982-1999 –> 2000-2009 4.310569 23496.9520
Indo-Pacific (27.95,28.05] 2 cstar_tref ~ temp + phosphate 1874 2.2971148 8443.2205 7.873503 1982-1999 –> 2000-2009 4.404428 23787.8221
Indo-Pacific (27.95,28.05] 2 cstar_tref ~ temp + silicate 1874 2.2575211 8378.0557 8.149044 1982-1999 –> 2000-2009 4.285968 23452.3719
Indo-Pacific (27.95,28.05] 3 cstar_tref ~ sal + aou + phosphate_star 1947 2.5390652 9163.7602 9.357185 2000-2009 –> 2010-2019 4.896999 17706.9225
Indo-Pacific (27.95,28.05] 3 cstar_tref ~ sal + oxygen + phosphate_star 1947 2.5579747 9192.6530 9.385134 2000-2009 –> 2010-2019 4.941144 17775.7151
Indo-Pacific (27.95,28.05] 3 cstar_tref ~ sal + phosphate + phosphate_star 1947 2.5579747 9192.6530 9.385134 2000-2009 –> 2010-2019 4.941144 17775.7151
Indo-Pacific (27.95,28.05] 2 cstar_tref ~ sal + silicate 1947 2.6712200 9359.3391 10.956872 2000-2009 –> 2010-2019 5.060698 17950.3095
Indo-Pacific (27.95,28.05] 3 cstar_tref ~ sal + silicate + phosphate_star 1947 2.5480622 9177.5340 9.329082 2000-2009 –> 2010-2019 4.891090 17696.9279
Indo-Pacific (27.95,28.05] 2 cstar_tref ~ temp + aou 1947 2.5254898 9140.8847 9.261142 2000-2009 –> 2010-2019 4.804441 17554.3512
Indo-Pacific (27.95,28.05] 3 cstar_tref ~ temp + aou + phosphate_star 1947 2.5129950 9123.5714 8.023733 2000-2009 –> 2010-2019 4.791613 17538.4899
Indo-Pacific (27.95,28.05] 2 cstar_tref ~ temp + oxygen 1947 2.5268145 9142.9267 9.264138 2000-2009 –> 2010-2019 4.806795 17558.0860
Indo-Pacific (27.95,28.05] 2 cstar_tref ~ temp + phosphate 1947 2.5143818 9123.7197 7.798725 2000-2009 –> 2010-2019 4.811497 17566.9402
Indo-Pacific (27.95,28.05] 2 cstar_tref ~ temp + silicate 1947 2.5076668 9113.3063 7.886716 2000-2009 –> 2010-2019 4.765188 17491.3620
Indo-Pacific (28.05,28.1] 2 cstar_tref ~ sal + aou 1324 2.0871217 5713.7104 9.756770 1982-1999 –> 2000-2009 3.874239 15935.3740
Indo-Pacific (28.05,28.1] 3 cstar_tref ~ sal + aou + phosphate_star 1324 2.0682780 5691.6942 9.298066 1982-1999 –> 2000-2009 3.841853 15876.5033
Indo-Pacific (28.05,28.1] 2 cstar_tref ~ sal + oxygen 1324 2.1146590 5748.4194 9.770060 1982-1999 –> 2000-2009 3.907413 15986.1683
Indo-Pacific (28.05,28.1] 3 cstar_tref ~ sal + oxygen + phosphate_star 1324 2.0857032 5713.9101 9.369106 1982-1999 –> 2000-2009 3.871617 15934.1309
Indo-Pacific (28.05,28.1] 3 cstar_tref ~ sal + phosphate + phosphate_star 1324 2.0857032 5713.9101 9.369106 1982-1999 –> 2000-2009 3.871617 15934.1309
Indo-Pacific (28.05,28.1] 2 cstar_tref ~ sal + silicate 1324 1.9839547 5579.4733 9.253841 1982-1999 –> 2000-2009 3.739676 15710.6018
Indo-Pacific (28.05,28.1] 3 cstar_tref ~ sal + silicate + phosphate_star 1324 1.9825148 5579.5508 9.144331 1982-1999 –> 2000-2009 3.634907 15402.8486
Indo-Pacific (28.05,28.1] 2 cstar_tref ~ temp + aou 1324 2.0377880 5650.3675 9.927150 1982-1999 –> 2000-2009 3.937903 16185.2055
Indo-Pacific (28.05,28.1] 2 cstar_tref ~ temp + oxygen 1324 2.0381058 5650.7804 9.934189 1982-1999 –> 2000-2009 3.938340 16185.9385
Indo-Pacific (28.05,28.1] 2 cstar_tref ~ temp + silicate 1324 1.9835305 5578.9070 8.959606 1982-1999 –> 2000-2009 3.870705 16078.8374
Indo-Pacific (28.05,28.1] 3 cstar_tref ~ sal + aou + phosphate_star 1360 2.1338565 5931.0848 9.707109 2000-2009 –> 2010-2019 4.202134 11622.7790
Indo-Pacific (28.05,28.1] 3 cstar_tref ~ sal + oxygen + phosphate_star 1360 2.1464235 5947.0569 9.644823 2000-2009 –> 2010-2019 4.232127 11660.9670
Indo-Pacific (28.05,28.1] 2 cstar_tref ~ sal + phosphate 1360 2.1632465 5966.2923 10.071374 2000-2009 –> 2010-2019 4.250487 11680.1536
Indo-Pacific (28.05,28.1] 3 cstar_tref ~ sal + phosphate + phosphate_star 1360 2.1464235 5947.0569 9.644823 2000-2009 –> 2010-2019 4.232127 11660.9670
Indo-Pacific (28.05,28.1] 2 cstar_tref ~ sal + silicate 1360 2.0784760 5857.5598 9.822309 2000-2009 –> 2010-2019 4.062431 11437.0331
Indo-Pacific (28.05,28.1] 3 cstar_tref ~ sal + silicate + phosphate_star 1360 2.0615054 5837.2601 9.429180 2000-2009 –> 2010-2019 4.044020 11416.8109
Indo-Pacific (28.05,28.1] 2 cstar_tref ~ temp + aou 1360 2.1658166 5969.5219 10.776863 2000-2009 –> 2010-2019 4.203605 11619.8894
Indo-Pacific (28.05,28.1] 2 cstar_tref ~ temp + oxygen 1360 2.1666773 5970.6027 10.783586 2000-2009 –> 2010-2019 4.204783 11621.3831
Indo-Pacific (28.05,28.1] 2 cstar_tref ~ temp + phosphate 1360 2.1556414 5956.7131 10.142926 2000-2009 –> 2010-2019 4.243729 11671.6492
Indo-Pacific (28.05,28.1] 2 cstar_tref ~ temp + silicate 1360 2.0672644 5842.8480 9.910053 2000-2009 –> 2010-2019 4.050795 11421.7550
Indo-Pacific (28.1, Inf] 3 cstar_tref ~ sal + silicate + phosphate_star 11224 1.8128765 45216.9805 7.901845 1982-1999 –> 2000-2009 3.544603 125267.8789
Indo-Pacific (28.1, Inf] 4 cstar_tref ~ sal + temp + aou + phosphate_star 11224 1.8000356 45059.4112 14.193464 1982-1999 –> 2000-2009 3.510502 124609.9295
Indo-Pacific (28.1, Inf] 4 cstar_tref ~ sal + temp + oxygen + phosphate_star 11224 1.7992024 45049.0180 14.198826 1982-1999 –> 2000-2009 3.508476 124571.1584
Indo-Pacific (28.1, Inf] 4 cstar_tref ~ sal + temp + phosphate + phosphate_star 11224 1.7992024 45049.0180 14.198826 1982-1999 –> 2000-2009 3.508476 124571.1584
Indo-Pacific (28.1, Inf] 4 cstar_tref ~ sal + temp + silicate + phosphate_star 11224 1.3186970 38074.4400 16.001226 1982-1999 –> 2000-2009 2.488354 102168.0855
Indo-Pacific (28.1, Inf] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 11224 1.6447617 43034.3568 17.962279 1982-1999 –> 2000-2009 3.077745 115385.9647
Indo-Pacific (28.1, Inf] 3 cstar_tref ~ temp + oxygen + phosphate_star 11224 1.8052481 45122.3215 13.229311 1982-1999 –> 2000-2009 3.577720 126119.0397
Indo-Pacific (28.1, Inf] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 11224 1.6475027 43071.7351 17.957785 1982-1999 –> 2000-2009 3.082878 115491.1817
Indo-Pacific (28.1, Inf] 3 cstar_tref ~ temp + phosphate + phosphate_star 11224 1.8052481 45122.3215 13.229311 1982-1999 –> 2000-2009 3.577720 126119.0397
Indo-Pacific (28.1, Inf] 4 cstar_tref ~ temp + silicate + phosphate + phosphate_star 11224 1.6475027 43071.7351 17.957785 1982-1999 –> 2000-2009 3.082878 115491.1817
Indo-Pacific (28.1, Inf] 4 cstar_tref ~ sal + temp + aou + phosphate_star 11207 1.9099517 46319.6941 12.285138 2000-2009 –> 2010-2019 3.709987 91379.1053
Indo-Pacific (28.1, Inf] 4 cstar_tref ~ sal + temp + oxygen + phosphate_star 11207 1.9092533 46311.4957 12.288301 2000-2009 –> 2010-2019 3.708456 91360.5137
Indo-Pacific (28.1, Inf] 4 cstar_tref ~ sal + temp + phosphate + phosphate_star 11207 1.9092533 46311.4957 12.288301 2000-2009 –> 2010-2019 3.708456 91360.5137
Indo-Pacific (28.1, Inf] 4 cstar_tref ~ sal + temp + silicate + phosphate_star 11207 1.4656759 40385.3306 13.168099 2000-2009 –> 2010-2019 2.784373 78459.7707
Indo-Pacific (28.1, Inf] 3 cstar_tref ~ temp + aou + phosphate_star 11207 1.9143080 46368.7584 13.200583 2000-2009 –> 2010-2019 3.721533 91515.6516
Indo-Pacific (28.1, Inf] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 11207 1.8503618 45609.2427 16.357517 2000-2009 –> 2010-2019 3.495124 88643.5995
Indo-Pacific (28.1, Inf] 3 cstar_tref ~ temp + oxygen + phosphate_star 11207 1.9145586 46371.6925 13.291798 2000-2009 –> 2010-2019 3.719807 91494.0140
Indo-Pacific (28.1, Inf] 4 cstar_tref ~ temp + oxygen + silicate + phosphate_star 11207 1.8534001 45646.0155 16.358305 2000-2009 –> 2010-2019 3.500903 88717.7506
Indo-Pacific (28.1, Inf] 3 cstar_tref ~ temp + phosphate + phosphate_star 11207 1.9145586 46371.6925 13.291798 2000-2009 –> 2010-2019 3.719807 91494.0140
Indo-Pacific (28.1, Inf] 4 cstar_tref ~ temp + silicate + phosphate + phosphate_star 11207 1.8534001 45646.0155 16.358305 2000-2009 –> 2010-2019 3.500903 88717.7506

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.431632 1107.9005
Atlantic (-Inf,26] 2000-2009 –> 2010-2019 2.395222 754.8143
Atlantic (26,26.5] 1982-1999 –> 2000-2009 9.867481 14138.6935
Atlantic (26,26.5] 2000-2009 –> 2010-2019 9.748212 9702.3883
Atlantic (26.5,26.75] 1982-1999 –> 2000-2009 7.399108 16068.8129
Atlantic (26.5,26.75] 2000-2009 –> 2010-2019 8.045563 11506.2553
Atlantic (26.75,27] 1982-1999 –> 2000-2009 5.760677 23974.5826
Atlantic (26.75,27] 2000-2009 –> 2010-2019 6.521093 18779.4253
Atlantic (27,27.25] 1982-1999 –> 2000-2009 5.076395 19509.1895
Atlantic (27,27.25] 2000-2009 –> 2010-2019 5.489946 15366.0759
Atlantic (27.25,27.5] 1982-1999 –> 2000-2009 8.224204 24858.0064
Atlantic (27.25,27.5] 2000-2009 –> 2010-2019 7.891214 17478.2448
Atlantic (27.5,27.75] 1982-1999 –> 2000-2009 6.840355 31164.9883
Atlantic (27.5,27.75] 2000-2009 –> 2010-2019 7.207101 23281.3589
Atlantic (27.75,27.85] 1982-1999 –> 2000-2009 8.530799 11620.4785
Atlantic (27.75,27.85] 2000-2009 –> 2010-2019 10.100108 9245.5914
Atlantic (27.85,27.95] 1982-1999 –> 2000-2009 9.882091 12442.2744
Atlantic (27.85,27.95] 2000-2009 –> 2010-2019 10.803501 9628.9879
Atlantic (27.95,28.05] 1982-1999 –> 2000-2009 13.080110 17116.6644
Atlantic (27.95,28.05] 2000-2009 –> 2010-2019 11.134138 11718.7934
Atlantic (28.05,28.1] 1982-1999 –> 2000-2009 6.059551 8372.8045
Atlantic (28.05,28.1] 2000-2009 –> 2010-2019 5.691987 6028.9323
Atlantic (28.1,28.15] 1982-1999 –> 2000-2009 4.684336 9003.7327
Atlantic (28.1,28.15] 2000-2009 –> 2010-2019 4.650461 6426.0454
Atlantic (28.15,28.2] 1982-1999 –> 2000-2009 2.440580 11065.2991
Atlantic (28.15,28.2] 2000-2009 –> 2010-2019 2.615747 7789.0458
Atlantic (28.2, Inf] 1982-1999 –> 2000-2009 1.940646 16278.3432
Atlantic (28.2, Inf] 2000-2009 –> 2010-2019 2.097053 12421.2892
Indo-Pacific (-Inf,26] 1982-1999 –> 2000-2009 13.524966 75123.2067
Indo-Pacific (-Inf,26] 2000-2009 –> 2010-2019 13.012661 53346.3632
Indo-Pacific (26,26.5] 1982-1999 –> 2000-2009 12.851522 78077.5931
Indo-Pacific (26,26.5] 2000-2009 –> 2010-2019 13.110093 55623.0838
Indo-Pacific (26.5,26.75] 1982-1999 –> 2000-2009 9.062097 57419.3718
Indo-Pacific (26.5,26.75] 2000-2009 –> 2010-2019 10.027767 42350.8932
Indo-Pacific (26.75,27] 1982-1999 –> 2000-2009 10.301434 70239.9960
Indo-Pacific (26.75,27] 2000-2009 –> 2010-2019 11.028455 51397.7341
Indo-Pacific (27,27.25] 1982-1999 –> 2000-2009 11.599832 87420.2178
Indo-Pacific (27,27.25] 2000-2009 –> 2010-2019 13.075817 65360.7655
Indo-Pacific (27.25,27.5] 1982-1999 –> 2000-2009 8.812680 71043.4486
Indo-Pacific (27.25,27.5] 2000-2009 –> 2010-2019 10.086310 54462.2200
Indo-Pacific (27.5,27.75] 1982-1999 –> 2000-2009 6.320350 63669.7721
Indo-Pacific (27.5,27.75] 2000-2009 –> 2010-2019 7.790300 49610.8472
Indo-Pacific (27.75,27.85] 1982-1999 –> 2000-2009 4.978080 22183.3112
Indo-Pacific (27.75,27.85] 2000-2009 –> 2010-2019 5.865554 17889.3209
Indo-Pacific (27.85,27.95] 1982-1999 –> 2000-2009 5.179044 28691.1341
Indo-Pacific (27.85,27.95] 2000-2009 –> 2010-2019 5.967389 22372.8087
Indo-Pacific (27.95,28.05] 1982-1999 –> 2000-2009 4.390230 23689.9518
Indo-Pacific (27.95,28.05] 2000-2009 –> 2010-2019 4.871061 17661.4819
Indo-Pacific (28.05,28.1] 1982-1999 –> 2000-2009 3.848827 15922.9739
Indo-Pacific (28.05,28.1] 2000-2009 –> 2010-2019 4.172624 11581.3387
Indo-Pacific (28.1, Inf] 1982-1999 –> 2000-2009 3.295935 119979.4618
Indo-Pacific (28.1, Inf] 2000-2009 –> 2010-2019 3.556935 89314.2684

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 - 5).

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
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
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
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
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 oxygen phosphate phosphate_star sal silicate temp
Atlantic - 1982-1999 --> 2000-2009
(-Inf,26] 4 4 2 6 10 5 10
(26,26.5] 5 3 2 7 2 10 6
(26.5,26.75] 4 3 3 10 6 6 4
(26.75,27] 4 3 3 10 6 6 4
(27,27.25] 3 3 3 10 7 7 3
(27.25,27.5] 4 2 4 8 4 0 5
(27.5,27.75] 3 3 5 7 6 2 4
(27.75,27.85] 4 2 6 5 2 1 2
(27.85,27.95] 2 2 6 5 2 1 2
(27.95,28.05] 3 2 6 4 3 2 1
(28.05,28.1] 3 3 4 3 2 4 1
(28.1,28.15] 3 3 3 3 1 5 2
(28.15,28.2] 2 3 4 3 4 4 1
(28.2, Inf] 4 6 5 5 3 4 0
total 48.00 42.00 56.00 86.00 58.00 57.00 45.00
Atlantic - 2000-2009 --> 2010-2019
(-Inf,26] 3 4 2 6 10 5 10
(26,26.5] 4 2 3 5 1 10 8
(26.5,26.75] 4 3 3 10 6 6 4
(26.75,27] 3 3 3 10 7 7 3
(27,27.25] 3 3 3 10 7 7 3
(27.25,27.5] 3 2 6 7 6 4 4
(27.5,27.75] 2 2 5 8 8 4 2
(27.75,27.85] 3 2 6 5 2 2 2
(27.85,27.95] 2 2 6 5 2 1 2
(27.95,28.05] 3 4 7 4 5 0 1
(28.05,28.1] 3 3 5 4 1 3 1
(28.1,28.15] 3 3 3 3 1 5 2
(28.15,28.2] 2 3 5 3 3 4 1
(28.2, Inf] 4 6 5 5 3 4 0
total 42.00 42.00 62.00 85.00 62.00 62.00 43.00
Indo-Pacific - 1982-1999 --> 2000-2009
(-Inf,26] 2 3 3 10 8 6 10
(26,26.5] 3 4 1 8 5 6 5
(26.5,26.75] 3 4 1 10 4 6 6
(26.75,27] 3 3 3 10 6 4 4
(27,27.25] 2 2 3 8 5 2 5
(27.25,27.5] 2 3 3 9 5 2 5
(27.5,27.75] 2 3 2 6 9 2 4
(27.75,27.85] 3 2 2 2 4 3 4
(27.85,27.95] 3 2 2 2 6 3 8
(27.95,28.05] 3 2 2 5 5 3 5
(28.05,28.1] 3 3 1 4 7 3 3
(28.1, Inf] 2 3 3 10 5 5 9
total 31.00 34.00 26.00 84.00 69.00 45.00 68.00
Indo-Pacific - 2000-2009 --> 2010-2019
(-Inf,26] 2 3 3 10 8 6 10
(26,26.5] 3 4 1 8 5 6 5
(26.5,26.75] 3 3 1 9 5 5 5
(26.75,27] 2 3 5 8 6 4 4
(27,27.25] 2 3 4 8 4 1 4
(27.25,27.5] 3 3 3 7 7 1 6
(27.5,27.75] 3 2 3 5 9 2 5
(27.75,27.85] 2 2 3 4 6 3 4
(27.85,27.95] 3 2 3 4 8 2 6
(27.95,28.05] 3 2 2 5 5 3 5
(28.05,28.1] 2 2 3 4 6 3 4
(28.1, Inf] 3 3 3 10 4 4 10
total 31.00 32.00 34.00 82.00 73.00 40.00 68.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
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
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.2      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_0.2.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