Last updated: 2021-03-11

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

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

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


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

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

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

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

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

1 Required data

Required are:

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

2 Predictor combinations

Find all possible combinations of following considered predictor variables:

  • sal, temp, aou, nitrate, silicate, phosphate, phosphate_star
# the following code is a workaround to find all predictor combinations
# using the olsrr package and fit all models for one era, slab, and basin

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

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

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

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

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

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

# convert to tibble
lm_all <- as_tibble(lm_all)

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

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

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

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

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

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

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

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

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

3 Apply predictor threshold

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

  • Minimum: 2
  • Maximum: 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)[1]
        
        # extract one model definition
        i_lm <- lm_all %>%
          filter(predictors == i_predictors) %>%
          select(lm_coeff) %>%
          pull()
        
        # extract number of predictors
        i_n_predictors <- lm_all %>%
          filter(predictors == i_predictors) %>%
          select(n) %>%
          pull()
        
        if (i_nr_obs > i_n_predictors) {
          # fit model
          if (params_local$MLR_type == "rlm") {
            i_lm_fit <- MASS::rlm(as.formula(i_lm),
                                  data = GLODAP_basin_era_slab)
          }
          
          if (params_local$MLR_type == "lm") {
            i_lm_fit <- lm(as.formula(i_lm),
                           data = GLODAP_basin_era_slab)
          }
          
          # find max predictor correlation
          i_cor_max <- GLODAP_basin_era_slab %>%
            select(!!!syms(str_split(i_predictors, " ",
                                     simplify = TRUE))) %>%
            correlate(quiet = TRUE) %>%
            select(-term) %>%
            abs() %>%
            max(na.rm = TRUE)
          
          # calculate root mean squared error
          i_rmse <- sqrt(c(crossprod(i_lm_fit$residuals)) /
                           length(i_lm_fit$residuals))
          
          # calculate maximum residual
          i_resid_max <- max(abs(i_lm_fit$residuals))
          
          # calculate Akaike information criterion aic
          i_aic <- AIC(i_lm_fit)
          
          # collect model coefficients and diagnostics
          coefficients <- tidy(i_lm_fit)
          
          coefficients <- coefficients %>%
            mutate(
              basin = i_basin,
              era = i_era,
              gamma_slab = i_gamma_slab,
              model = i_lm,
              nr_obs = i_nr_obs,
              rmse = i_rmse,
              aic = i_aic,
              resid_max = i_resid_max,
              n_predictors = i_n_predictors,
              na_predictor = anyNA(coefficients$estimate),
              cor_max = i_cor_max
            )
          
          if (exists("lm_all_fitted")) {
            lm_all_fitted <- bind_rows(lm_all_fitted, coefficients)
          }
          
          if (!exists("lm_all_fitted")) {
            lm_all_fitted <- coefficients
          }
        }
      }
    }
  }
}

rm(
  i_lm_fit,
  coefficients,
  i_rmse,
  GLODAP_basin_era,
  GLODAP_basin_era_slab,
  i_lm,
  i_basin,
  i_era,
  i_gamma_slab,
  i_nr_obs,
  i_predictors,
  i_aic,
  i_n_predictors,
  i_resid_max
)

5 Prepare coeffcients

Coefficients are prepared for the mapping of Cant and the chosen target variable.

5.1 Formatting

# select relevant columns
lm_all_fitted <- lm_all_fitted %>%
  select(
    basin,
    gamma_slab,
    era,
    model,
    nr_obs,
    n_predictors,
    term,
    estimate,
    rmse,
    aic,
    resid_max,
    na_predictor,
    cor_max
  )

# set coefficient to zero if not fitted (=NA)
lm_all_fitted <- lm_all_fitted %>%
  mutate(estimate = if_else(is.na(estimate), 0, estimate))

# Prepare model coefficients for mapping of target variable
lm_all_fitted_wide <- lm_all_fitted %>%
  pivot_wider(
    values_from = estimate,
    names_from = term,
    names_prefix = "coeff_",
    values_fill = 0
  )

5.2 Predictor selection

Within each basin and slab, the following number of best linear regression models was selected:

  • 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)
# calculate RMSE sum for adjacent eras
lm_all_fitted_wide_eras <- lm_all_fitted_wide %>%
  select(basin, gamma_slab, n_predictors, model, era, nr_obs, rmse, aic, resid_max) %>%
  arrange(era) %>%
  group_by(basin, gamma_slab, model) %>%
  mutate(
    eras = paste(lag(era), era, sep = " --> "),
    rmse_sum = rmse + lag(rmse),
    aic_sum = aic + lag(aic)
  ) %>%
  ungroup() %>%
  select(-c(era)) %>%
  drop_na()

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

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

5.3 RMSE tables

5.3.1 per model

lm_best %>%
  kable() %>%
  add_header_above() %>%
  kable_styling() %>%
  scroll_box(width = "100%", height = "400px")
basin gamma_slab n_predictors model nr_obs rmse aic resid_max eras rmse_sum aic_sum
Atlantic (-Inf,26] 5 cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 158 1.0382074 474.2332 3.944222 1987-1999 –> 2000-2012 2.3374558 975.6433
Atlantic (-Inf,26] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 158 1.0492038 475.5626 4.357844 1987-1999 –> 2000-2012 2.3484591 974.9742
Atlantic (-Inf,26] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 158 1.0480832 477.2249 4.299021 1987-1999 –> 2000-2012 2.3469374 978.5470
Atlantic (-Inf,26] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate 158 1.0440638 476.0107 4.148260 1987-1999 –> 2000-2012 2.3453010 977.8644
Atlantic (-Inf,26] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 158 1.0458289 476.5445 4.199297 1987-1999 –> 2000-2012 2.3474545 978.4847
Atlantic (-Inf,26] 4 cstar_tref ~ sal + temp + aou + phosphate 158 1.0547002 477.2137 4.547311 1987-1999 –> 2000-2012 2.3559722 977.0751
Atlantic (-Inf,26] 5 cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 158 1.0447877 476.2297 3.991338 1987-1999 –> 2000-2012 2.3437009 977.5650
Atlantic (-Inf,26] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 158 1.0538228 478.9507 4.496846 1987-1999 –> 2000-2012 2.3544566 980.6698
Atlantic (-Inf,26] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 158 1.0553928 479.4211 4.540758 1987-1999 –> 2000-2012 2.3563544 981.2133
Atlantic (-Inf,26] 5 cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 158 1.0478179 477.1449 4.250178 1987-1999 –> 2000-2012 2.3504756 979.3150
Atlantic (-Inf,26] 5 cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 72 1.1326926 236.2692 2.855476 2000-2012 –> 2013-2019 2.1709001 710.5024
Atlantic (-Inf,26] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 72 1.1936365 243.8158 3.237974 2000-2012 –> 2013-2019 2.2417197 721.0407
Atlantic (-Inf,26] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate 72 1.1361477 236.7078 2.802670 2000-2012 –> 2013-2019 2.1802115 712.7185
Atlantic (-Inf,26] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 72 1.1377961 236.9166 2.786570 2000-2012 –> 2013-2019 2.1836250 713.4610
Atlantic (-Inf,26] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 72 1.2083604 245.5812 2.691226 2000-2012 –> 2013-2019 2.2954998 734.3676
Atlantic (-Inf,26] 5 cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 72 1.2718853 252.9592 3.547879 2000-2012 –> 2013-2019 2.3166730 729.1889
Atlantic (-Inf,26] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 72 1.1990801 244.4710 3.379652 2000-2012 –> 2013-2019 2.2529029 723.4217
Atlantic (-Inf,26] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 72 1.2012446 244.7307 3.425080 2000-2012 –> 2013-2019 2.2566374 724.1518
Atlantic (-Inf,26] 5 cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 72 1.1396017 237.1449 2.772460 2000-2012 –> 2013-2019 2.1874195 714.2898
Atlantic (-Inf,26] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 72 1.2024324 244.8730 3.445295 2000-2012 –> 2013-2019 2.2595421 724.8078
Atlantic (26,26.5] 5 cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 1046 3.3042550 5482.8009 13.811611 1987-1999 –> 2000-2012 6.5814104 11203.7533
Atlantic (26,26.5] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1046 3.5518213 5633.9468 13.477555 1987-1999 –> 2000-2012 7.0260971 11482.8176
Atlantic (26,26.5] 4 cstar_tref ~ sal + temp + aou + nitrate 1046 3.4447540 5567.9148 15.626500 1987-1999 –> 2000-2012 6.8877638 11394.9881
Atlantic (26,26.5] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate 1046 3.3726755 5525.6771 14.465145 1987-1999 –> 2000-2012 6.7472012 11310.7504
Atlantic (26,26.5] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 1046 3.3758183 5527.6256 14.515788 1987-1999 –> 2000-2012 6.7543773 11315.3149
Atlantic (26,26.5] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 1046 3.3139957 5488.9589 14.956164 1987-1999 –> 2000-2012 6.6586043 11254.5301
Atlantic (26,26.5] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 1046 3.6117860 5668.9708 14.137864 1987-1999 –> 2000-2012 7.1751096 11573.2654
Atlantic (26,26.5] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1046 3.6195773 5673.4788 14.246703 1987-1999 –> 2000-2012 7.1938295 11584.4798
Atlantic (26,26.5] 5 cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 1046 3.3921889 5537.7460 14.699641 1987-1999 –> 2000-2012 6.7939822 11340.4443
Atlantic (26,26.5] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 1046 3.6210378 5674.3227 14.283746 1987-1999 –> 2000-2012 7.1996234 11587.9772
Atlantic (26,26.5] 5 cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 567 3.2403908 2956.3132 10.534106 2000-2012 –> 2013-2019 6.5446458 8439.1141
Atlantic (26,26.5] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 567 3.6070302 3077.8676 10.902628 2000-2012 –> 2013-2019 7.1588515 8711.8144
Atlantic (26,26.5] 4 cstar_tref ~ sal + temp + aou + nitrate 567 3.3369968 2987.6271 10.231047 2000-2012 –> 2013-2019 6.7817508 8555.5419
Atlantic (26,26.5] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate 567 3.2965285 2975.7908 9.955222 2000-2012 –> 2013-2019 6.6692041 8501.4679
Atlantic (26,26.5] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 567 3.2984682 2976.4579 9.956192 2000-2012 –> 2013-2019 6.6742865 8504.0835
Atlantic (26,26.5] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 567 3.1723161 2932.2362 9.432214 2000-2012 –> 2013-2019 6.4863118 8421.1950
Atlantic (26,26.5] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 567 3.6470365 3090.3758 10.325180 2000-2012 –> 2013-2019 7.2588226 8759.3466
Atlantic (26,26.5] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 567 3.6519756 3091.9105 10.331554 2000-2012 –> 2013-2019 7.2715530 8765.3893
Atlantic (26,26.5] 5 cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 567 3.3102918 2980.5155 9.977380 2000-2012 –> 2013-2019 6.7024807 8518.2615
Atlantic (26,26.5] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 567 3.6516653 3091.8142 10.332889 2000-2012 –> 2013-2019 7.2727031 8766.1369
Atlantic (26.5,26.75] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 1397 3.2423787 7265.1167 9.560288 1987-1999 –> 2000-2012 6.3261737 14435.9365
Atlantic (26.5,26.75] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 1397 3.3429532 7348.4661 10.697715 1987-1999 –> 2000-2012 6.4791235 14564.6440
Atlantic (26.5,26.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate 1397 3.3033338 7317.1548 9.343758 1987-1999 –> 2000-2012 6.4531807 14547.5690
Atlantic (26.5,26.75] 4 cstar_tref ~ sal + aou + silicate + phosphate 1397 3.3316510 7339.0038 10.073920 1987-1999 –> 2000-2012 6.4918567 14576.6505
Atlantic (26.5,26.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1397 3.2700730 7288.8799 10.507891 1987-1999 –> 2000-2012 6.3573995 14462.9181
Atlantic (26.5,26.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 1397 3.2846060 7301.2695 10.500101 1987-1999 –> 2000-2012 6.3930466 14494.4734
Atlantic (26.5,26.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1397 3.3014232 7315.5384 10.543653 1987-1999 –> 2000-2012 6.4318399 14528.5527
Atlantic (26.5,26.75] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 1397 3.3059150 7319.3372 10.550809 1987-1999 –> 2000-2012 6.4430536 14538.3832
Atlantic (26.5,26.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate 1397 3.2989930 7313.4809 9.947012 1987-1999 –> 2000-2012 6.4429440 14538.6266
Atlantic (26.5,26.75] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1397 3.2249399 7250.0489 10.037456 1987-1999 –> 2000-2012 6.2363263 14354.0544
Atlantic (26.5,26.75] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 656 3.3433340 3459.1900 9.454244 2000-2012 –> 2013-2019 6.5857127 10724.3067
Atlantic (26.5,26.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate 656 3.3807232 3473.7810 10.074129 2000-2012 –> 2013-2019 6.6840570 10790.9358
Atlantic (26.5,26.75] 4 cstar_tref ~ sal + aou + silicate + phosphate 656 3.3850399 3473.4551 10.375376 2000-2012 –> 2013-2019 6.7166909 10812.4589
Atlantic (26.5,26.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 656 3.3455096 3460.0435 10.735554 2000-2012 –> 2013-2019 6.6155826 10748.9234
Atlantic (26.5,26.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 656 3.3565540 3464.3677 10.767824 2000-2012 –> 2013-2019 6.6411600 10765.6372
Atlantic (26.5,26.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 656 3.3723676 3470.5343 10.887326 2000-2012 –> 2013-2019 6.6737908 10786.0727
Atlantic (26.5,26.75] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 656 3.3762721 3472.0524 10.919163 2000-2012 –> 2013-2019 6.6821871 10791.3896
Atlantic (26.5,26.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate 656 3.3892046 3477.0684 9.580435 2000-2012 –> 2013-2019 6.6881976 10790.5493
Atlantic (26.5,26.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 656 3.4093443 3484.8416 9.603997 2000-2012 –> 2013-2019 6.7344151 10820.3216
Atlantic (26.5,26.75] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 656 3.3238240 3451.5114 11.113656 2000-2012 –> 2013-2019 6.5487639 10701.5603
Atlantic (26.75,27] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 2398 2.0422179 10243.7478 13.072427 1987-1999 –> 2000-2012 3.7023665 18916.3499
Atlantic (26.75,27] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 2398 2.0935628 10360.8368 14.739529 1987-1999 –> 2000-2012 3.7859338 19117.8672
Atlantic (26.75,27] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2398 2.0001820 10143.9994 14.081577 1987-1999 –> 2000-2012 3.6755570 18857.6495
Atlantic (26.75,27] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 2398 2.0669164 10301.4026 14.196373 1987-1999 –> 2000-2012 3.8121755 19198.7861
Atlantic (26.75,27] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2398 2.1361037 10459.3140 22.290539 1987-1999 –> 2000-2012 3.9210989 19457.9146
Atlantic (26.75,27] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 2398 2.1629225 10519.1528 25.264079 1987-1999 –> 2000-2012 3.9620563 19553.2254
Atlantic (26.75,27] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate 2398 2.1370107 10461.3499 13.373762 1987-1999 –> 2000-2012 3.9024816 19410.5024
Atlantic (26.75,27] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2398 2.1989701 10598.4250 16.309198 1987-1999 –> 2000-2012 4.0099938 19662.1122
Atlantic (26.75,27] 4 cstar_tref ~ temp + aou + silicate + phosphate 2398 2.2004206 10599.5875 15.280820 1987-1999 –> 2000-2012 4.0054217 19646.2983
Atlantic (26.75,27] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 2398 2.0860637 10345.6268 19.975942 1987-1999 –> 2000-2012 3.6269684 18683.1087
Atlantic (26.75,27] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 1263 2.5145378 5927.4356 22.300060 2000-2012 –> 2013-2019 4.5567556 16171.1833
Atlantic (26.75,27] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 1263 2.6364762 6045.0524 27.471250 2000-2012 –> 2013-2019 4.7300390 16405.8891
Atlantic (26.75,27] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate 1263 2.5534340 5966.2100 22.124060 2000-2012 –> 2013-2019 4.7278634 16510.8101
Atlantic (26.75,27] 4 cstar_tref ~ sal + aou + silicate + phosphate 1263 2.5846959 5994.9482 24.744636 2000-2012 –> 2013-2019 4.7623441 16544.6427
Atlantic (26.75,27] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1263 2.5330762 5945.9902 25.891764 2000-2012 –> 2013-2019 4.5332582 16089.9896
Atlantic (26.75,27] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 1263 2.5653762 5977.9963 25.592327 2000-2012 –> 2013-2019 4.6322925 16279.3988
Atlantic (26.75,27] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1263 2.5985873 6010.4878 25.748151 2000-2012 –> 2013-2019 4.7346911 16469.8018
Atlantic (26.75,27] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 1263 2.6096389 6021.2079 25.770703 2000-2012 –> 2013-2019 4.7725614 16540.3607
Atlantic (26.75,27] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate 1263 2.6260789 6037.0711 22.609023 2000-2012 –> 2013-2019 4.7630896 16498.4210
Atlantic (26.75,27] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1263 2.2552518 5652.5376 22.583567 2000-2012 –> 2013-2019 4.3413155 15998.1644
Atlantic (27,27.25] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 2031 1.9602837 8511.8166 15.249340 1987-1999 –> 2000-2012 3.6882340 16102.2969
Atlantic (27,27.25] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate 2031 1.8626152 8304.2173 11.021443 1987-1999 –> 2000-2012 3.6130762 15944.5808
Atlantic (27,27.25] 4 cstar_tref ~ sal + aou + silicate + phosphate 2031 1.9148200 8414.4995 9.235730 1987-1999 –> 2000-2012 3.6840903 16094.0545
Atlantic (27,27.25] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2031 1.8904789 8364.5325 9.689213 1987-1999 –> 2000-2012 3.5155361 15718.4044
Atlantic (27,27.25] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 2031 1.8979717 8380.6003 9.355072 1987-1999 –> 2000-2012 3.5191781 15725.3287
Atlantic (27,27.25] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2031 1.9187787 8424.8886 8.642107 1987-1999 –> 2000-2012 3.5586300 15813.6875
Atlantic (27,27.25] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 2031 1.9234272 8434.7175 8.493771 1987-1999 –> 2000-2012 3.5707920 15841.1343
Atlantic (27,27.25] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate 2031 2.0463518 8686.3583 13.496221 1987-1999 –> 2000-2012 3.7569690 16237.9836
Atlantic (27,27.25] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2031 2.0594140 8712.2044 13.384239 1987-1999 –> 2000-2012 3.7830359 16293.0183
Atlantic (27,27.25] 5 cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 2031 2.1067869 8804.5844 12.394761 1987-1999 –> 2000-2012 3.8609144 16453.0121
Atlantic (27,27.25] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 1142 2.4331867 5285.7925 20.133740 2000-2012 –> 2013-2019 4.3934703 13797.6091
Atlantic (27,27.25] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate 1142 2.1593133 5013.0566 10.465472 2000-2012 –> 2013-2019 4.0219286 13317.2739
Atlantic (27,27.25] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1142 2.2408412 5097.7041 11.555324 2000-2012 –> 2013-2019 4.3317934 13871.6430
Atlantic (27,27.25] 4 cstar_tref ~ sal + aou + silicate + phosphate 1142 2.3003827 5155.6000 13.777336 2000-2012 –> 2013-2019 4.2152027 13570.0995
Atlantic (27,27.25] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1142 2.2946948 5151.9457 14.080782 2000-2012 –> 2013-2019 4.1851737 13516.4782
Atlantic (27,27.25] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 1142 2.3424794 5197.0192 15.033313 2000-2012 –> 2013-2019 4.4776528 14053.9689
Atlantic (27,27.25] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 1142 2.2882847 5145.5566 14.231726 2000-2012 –> 2013-2019 4.1862565 13526.1569
Atlantic (27,27.25] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1142 2.3093184 5166.4549 14.219293 2000-2012 –> 2013-2019 4.2280971 13591.3435
Atlantic (27,27.25] 5 cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 1142 2.3856862 5240.7635 15.221464 2000-2012 –> 2013-2019 4.4201181 13903.3917
Atlantic (27,27.25] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 1142 2.3127467 5169.8432 14.213647 2000-2012 –> 2013-2019 4.2361739 13604.5606
Atlantic (27.25,27.5] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 2025 2.2571186 9057.7617 21.333316 1987-1999 –> 2000-2012 3.9161317 17172.8410
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate 2025 2.4835165 9444.8868 16.071294 1987-1999 –> 2000-2012 4.9362847 19205.2704
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 2025 2.4514644 9392.2776 12.813379 1987-1999 –> 2000-2012 4.9420057 19216.9712
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2025 2.6905754 9769.2091 21.422922 1987-1999 –> 2000-2012 4.9620143 19206.4025
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 2025 2.1625100 8884.3428 12.574713 1987-1999 –> 2000-2012 3.8895412 17168.5036
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 2025 2.7127786 9802.4934 22.028851 1987-1999 –> 2000-2012 5.0507483 19361.1692
Atlantic (27.25,27.5] 4 cstar_tref ~ temp + aou + nitrate + silicate 2025 2.4063497 9315.0504 11.614926 1987-1999 –> 2000-2012 4.1655142 17674.7861
Atlantic (27.25,27.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate 2025 2.2664591 9074.4871 13.372191 1987-1999 –> 2000-2012 4.0158947 17412.8863
Atlantic (27.25,27.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2025 2.2547324 9053.4778 15.509020 1987-1999 –> 2000-2012 3.9954999 17370.9755
Atlantic (27.25,27.5] 5 cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 2025 2.3326197 9191.0186 10.601623 1987-1999 –> 2000-2012 4.2011875 17806.6367
Atlantic (27.25,27.5] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 1109 2.7884162 5435.7065 20.513112 2000-2012 –> 2013-2019 5.0455349 14493.4682
Atlantic (27.25,27.5] 4 cstar_tref ~ aou + nitrate + silicate + phosphate_star 1109 3.0051254 5599.7139 21.412798 2000-2012 –> 2013-2019 5.7555576 15456.0350
Atlantic (27.25,27.5] 4 cstar_tref ~ sal + aou + nitrate + silicate 1109 3.0025580 5597.8181 16.245244 2000-2012 –> 2013-2019 5.5311321 15113.5242
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate 1109 2.9863807 5587.8356 14.517111 2000-2012 –> 2013-2019 5.4698972 15032.7224
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1109 2.7952620 5441.1451 16.992150 2000-2012 –> 2013-2019 5.2467264 14833.4228
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 1109 2.8205406 5461.1132 14.090309 2000-2012 –> 2013-2019 4.9830506 14345.4560
Atlantic (27.25,27.5] 4 cstar_tref ~ temp + aou + nitrate + silicate 1109 3.1081847 5674.5038 12.222568 2000-2012 –> 2013-2019 5.5145344 14989.5543
Atlantic (27.25,27.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate 1109 2.8123363 5454.6522 14.111336 2000-2012 –> 2013-2019 5.0787954 14529.1392
Atlantic (27.25,27.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1109 2.7923328 5438.8197 16.610907 2000-2012 –> 2013-2019 5.0470652 14492.2975
Atlantic (27.25,27.5] 5 cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 1109 2.8746758 5503.2804 14.300893 2000-2012 –> 2013-2019 5.2072955 14694.2990
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 2804 2.5430119 13205.6294 18.070467 1987-1999 –> 2000-2012 4.8148742 25362.4068
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate 2804 2.3843613 12844.3739 16.074344 1987-1999 –> 2000-2012 4.7117043 25131.9701
Atlantic (27.5,27.75] 4 cstar_tref ~ sal + aou + silicate + phosphate 2804 2.4538915 13003.5694 13.901382 1987-1999 –> 2000-2012 4.8369150 25417.3567
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2804 2.3944064 12867.9502 14.685866 1987-1999 –> 2000-2012 4.5455109 24728.5863
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate 2804 2.5420080 13203.4151 18.096781 1987-1999 –> 2000-2012 4.8090346 25348.6393
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 2804 2.3920617 12862.4560 14.587684 1987-1999 –> 2000-2012 4.5190958 24662.0790
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2804 2.4132697 12911.9574 14.216479 1987-1999 –> 2000-2012 4.5476976 24730.3953
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 2804 2.4148672 12915.6685 14.181558 1987-1999 –> 2000-2012 4.5508765 24738.1219
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 2804 2.4122515 12909.5908 14.214259 1987-1999 –> 2000-2012 4.5475429 24730.2217
Atlantic (27.5,27.75] 4 cstar_tref ~ sal + temp + silicate + phosphate_star 2804 2.4316201 12952.4392 14.640607 1987-1999 –> 2000-2012 4.5888204 24826.4182
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 1515 2.8849775 7523.7206 21.626774 2000-2012 –> 2013-2019 5.4279894 20729.3500
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate 1515 2.7078202 7331.6999 19.383046 2000-2012 –> 2013-2019 5.0921815 20176.0738
Atlantic (27.5,27.75] 4 cstar_tref ~ sal + aou + silicate + phosphate 1515 2.8745087 7510.7055 15.545879 2000-2012 –> 2013-2019 5.3284002 20514.2749
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1515 2.8711798 7509.1946 15.461950 2000-2012 –> 2013-2019 5.2655862 20377.1448
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate 1515 2.8717546 7509.8010 21.184103 2000-2012 –> 2013-2019 5.4137626 20713.2162
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 1515 2.8690598 7506.9564 15.594453 2000-2012 –> 2013-2019 5.2611214 20369.4124
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1515 2.8670393 7504.8218 16.960962 2000-2012 –> 2013-2019 5.2803090 20416.7793
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 1515 2.8833164 7521.9755 16.841839 2000-2012 –> 2013-2019 5.2981836 20437.6440
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 1515 2.8660787 7503.8064 16.855872 2000-2012 –> 2013-2019 5.2783302 20413.3972
Atlantic (27.5,27.75] 4 cstar_tref ~ sal + temp + silicate + phosphate_star 1515 2.9784908 7618.3764 19.270332 2000-2012 –> 2013-2019 5.4101109 20570.8156
Atlantic (27.75,27.85] 4 cstar_tref ~ sal + aou + nitrate + phosphate 1036 1.6166443 3947.3312 14.036706 1987-1999 –> 2000-2012 3.2976999 7653.8462
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 1036 1.5100084 3807.9434 14.539112 1987-1999 –> 2000-2012 3.0456767 7344.0487
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate 1036 1.6095022 3940.1571 13.113728 1987-1999 –> 2000-2012 3.2773491 7633.6369
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1036 1.5641536 3880.9392 12.798383 1987-1999 –> 2000-2012 3.1675100 7499.2570
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate 1036 1.5197119 3821.2156 14.808950 1987-1999 –> 2000-2012 3.0697649 7375.0917
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 1036 1.6216292 3955.7103 16.087961 1987-1999 –> 2000-2012 3.2479679 7601.1547
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 1036 1.6994499 4052.8317 11.282662 1987-1999 –> 2000-2012 3.2700765 7631.8394
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 1036 1.5917451 3917.1705 12.897616 1987-1999 –> 2000-2012 3.2284271 7574.6983
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 1036 1.5430795 3852.8330 15.324821 1987-1999 –> 2000-2012 3.1220204 7441.9036
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 1036 1.6298056 3966.1312 13.022129 1987-1999 –> 2000-2012 3.3095272 7673.1332
Atlantic (27.75,27.85] 4 cstar_tref ~ sal + aou + nitrate + phosphate 513 1.7818523 2060.5034 13.626298 2000-2012 –> 2013-2019 3.3984967 6007.8346
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 513 1.7094323 2019.9324 13.736344 2000-2012 –> 2013-2019 3.2194407 5827.8758
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate 513 1.7344893 2034.8625 11.391059 2000-2012 –> 2013-2019 3.3439915 5975.0196
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 513 1.7475521 2042.5606 12.039326 2000-2012 –> 2013-2019 3.3117057 5923.4997
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate 513 1.7242945 2028.8142 13.967079 2000-2012 –> 2013-2019 3.2440064 5850.0298
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 513 1.7754726 2058.8233 12.075398 2000-2012 –> 2013-2019 3.3672178 5975.9938
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 513 1.8242119 2086.6088 12.213063 2000-2012 –> 2013-2019 3.4614785 6062.2038
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 513 1.7556052 2047.2777 14.468303 2000-2012 –> 2013-2019 3.2986847 5900.1107
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 513 1.8185913 2083.4428 12.157352 2000-2012 –> 2013-2019 3.4483969 6049.5740
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 513 1.8275241 2088.4700 12.223997 2000-2012 –> 2013-2019 3.4690021 6069.3878
Atlantic (27.85,27.95] 4 cstar_tref ~ aou + nitrate + phosphate + phosphate_star 983 3.4411154 5231.2074 39.186076 1987-1999 –> 2000-2012 6.3536260 10005.1514
Atlantic (27.85,27.95] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 983 3.4269579 5225.1022 39.102145 1987-1999 –> 2000-2012 6.2954509 9971.8985
Atlantic (27.85,27.95] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 983 3.4332878 5226.7302 38.381356 1987-1999 –> 2000-2012 6.3174152 9981.9303
Atlantic (27.85,27.95] 5 cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 983 3.4355650 5230.0338 35.139570 1987-1999 –> 2000-2012 6.1339844 9859.8452
Atlantic (27.85,27.95] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 983 3.4396213 5230.3536 36.781212 1987-1999 –> 2000-2012 6.1567023 9871.3563
Atlantic (27.85,27.95] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 983 3.4272905 5225.2930 34.214054 1987-1999 –> 2000-2012 6.1443635 9868.2900
Atlantic (27.85,27.95] 5 cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 983 2.9304503 4917.3900 22.456892 1987-1999 –> 2000-2012 5.6457831 9559.1608
Atlantic (27.85,27.95] 5 cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 983 3.0821699 5016.6293 23.370605 1987-1999 –> 2000-2012 5.9862528 9787.0269
Atlantic (27.85,27.95] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 983 3.0946626 5022.5819 23.332862 1987-1999 –> 2000-2012 6.0211490 9805.6882
Atlantic (27.85,27.95] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 983 3.0757549 5012.5332 23.130965 1987-1999 –> 2000-2012 5.9457038 9760.3007
Atlantic (27.85,27.95] 5 cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 569 3.2198068 2959.4398 28.682477 2000-2012 –> 2013-2019 6.6553718 8189.4735
Atlantic (27.85,27.95] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 569 3.4103854 3022.8793 30.447244 2000-2012 –> 2013-2019 6.8500067 8253.2329
Atlantic (27.85,27.95] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 569 3.3215129 2994.8304 29.207068 2000-2012 –> 2013-2019 6.7488034 8220.1234
Atlantic (27.85,27.95] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate 569 3.2167027 2958.3421 28.783818 2000-2012 –> 2013-2019 6.9174456 8334.5523
Atlantic (27.85,27.95] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 569 3.2269879 2961.9750 29.065190 2000-2012 –> 2013-2019 7.0153766 8384.2037
Atlantic (27.85,27.95] 5 cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 569 3.3261262 2996.4099 29.121490 2000-2012 –> 2013-2019 6.2565764 7913.7999
Atlantic (27.85,27.95] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 569 3.3417360 3001.7382 29.366995 2000-2012 –> 2013-2019 7.0328338 8372.8178
Atlantic (27.85,27.95] 5 cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 569 3.6712299 3108.7515 22.690294 2000-2012 –> 2013-2019 6.7533998 8125.3808
Atlantic (27.85,27.95] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 569 3.6730463 3107.3144 22.693963 2000-2012 –> 2013-2019 6.7677089 8129.8962
Atlantic (27.85,27.95] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 569 3.6600265 3105.2733 23.151675 2000-2012 –> 2013-2019 6.7357814 8117.8065
Atlantic (27.95,28.05] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 1204 4.0178610 6779.7292 25.762442 1987-1999 –> 2000-2012 7.6113769 13366.2268
Atlantic (27.95,28.05] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate 1204 4.1624092 6864.8384 29.399782 1987-1999 –> 2000-2012 7.6408571 13372.0567
Atlantic (27.95,28.05] 4 cstar_tref ~ sal + nitrate + silicate + phosphate 1204 4.2522303 6914.2482 32.185459 1987-1999 –> 2000-2012 7.7900598 13460.7011
Atlantic (27.95,28.05] 5 cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 1204 4.1957717 6884.0621 30.223807 1987-1999 –> 2000-2012 7.6850894 13398.8807
Atlantic (27.95,28.05] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 1204 3.9558395 6742.2684 22.787482 1987-1999 –> 2000-2012 7.7117210 13436.4174
Atlantic (27.95,28.05] 5 cstar_tref ~ sal + temp + nitrate + silicate + phosphate 1204 4.2268913 6901.8560 32.270727 1987-1999 –> 2000-2012 7.7644409 13450.1162
Atlantic (27.95,28.05] 5 cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 1204 3.8735532 6691.6507 19.419274 1987-1999 –> 2000-2012 7.7277431 13448.7408
Atlantic (27.95,28.05] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate 1204 3.9682210 6749.7934 25.265565 1987-1999 –> 2000-2012 7.5342656 13317.5971
Atlantic (27.95,28.05] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1204 3.8590720 6682.6316 25.099675 1987-1999 –> 2000-2012 7.3109407 13171.1646
Atlantic (27.95,28.05] 5 cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 1204 3.9591041 6744.2548 25.125796 1987-1999 –> 2000-2012 7.5205468 13308.9127
Atlantic (27.95,28.05] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 666 3.5624057 3596.2470 28.359612 2000-2012 –> 2013-2019 7.5802666 10375.9761
Atlantic (27.95,28.05] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate 666 3.4782139 3564.3893 31.711472 2000-2012 –> 2013-2019 7.6406231 10429.2278
Atlantic (27.95,28.05] 5 cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 666 3.4861342 3567.4190 32.154120 2000-2012 –> 2013-2019 7.6819059 10451.4811
Atlantic (27.95,28.05] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 666 3.5357472 3586.2418 25.590197 2000-2012 –> 2013-2019 7.4915867 10328.5101
Atlantic (27.95,28.05] 5 cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 666 3.6379372 3624.1933 25.188475 2000-2012 –> 2013-2019 7.6972651 10428.6472
Atlantic (27.95,28.05] 5 cstar_tref ~ sal + temp + nitrate + silicate + phosphate 666 3.4602252 3557.4826 32.737359 2000-2012 –> 2013-2019 7.6871165 10459.3386
Atlantic (27.95,28.05] 5 cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 666 3.7631817 3669.2788 21.509432 2000-2012 –> 2013-2019 7.6367350 10360.9296
Atlantic (27.95,28.05] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate 666 3.5417278 3588.4929 27.962736 2000-2012 –> 2013-2019 7.5099488 10338.2863
Atlantic (27.95,28.05] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 666 3.5060655 3575.0128 26.986021 2000-2012 –> 2013-2019 7.3651375 10257.6443
Atlantic (27.95,28.05] 5 cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 666 3.5394291 3587.6281 27.891571 2000-2012 –> 2013-2019 7.4985332 10331.8828
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 804 1.0231680 2332.4823 7.031639 1987-1999 –> 2000-2012 1.9182588 4410.6774
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 804 1.0324411 2346.9902 6.945256 1987-1999 –> 2000-2012 1.9298612 4429.2864
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate 804 1.0272609 2338.9019 6.919786 1987-1999 –> 2000-2012 1.9279534 4426.9417
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 804 1.0374154 2354.7189 6.772617 1987-1999 –> 2000-2012 1.9480301 4460.0471
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 804 0.9934328 2285.0583 8.087713 1987-1999 –> 2000-2012 1.8635240 4318.5531
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 804 1.0379653 2355.5710 6.801673 1987-1999 –> 2000-2012 1.9431652 4451.4881
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 804 1.0296888 2342.6978 6.877082 1987-1999 –> 2000-2012 1.9344359 4437.8253
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 804 1.0412160 2360.5991 6.747646 1987-1999 –> 2000-2012 1.9518821 4466.0165
Atlantic (28.05,28.1] 5 cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 804 1.0479786 2371.0092 11.443033 1987-1999 –> 2000-2012 1.8779453 4330.0028
Atlantic (28.05,28.1] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 804 1.0436142 2364.2985 12.188672 1987-1999 –> 2000-2012 1.8317195 4241.6245
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 434 1.0460351 1284.7046 7.627744 2000-2012 –> 2013-2019 2.0692031 3617.1869
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 434 1.0586072 1295.0748 7.429722 2000-2012 –> 2013-2019 2.0910483 3642.0650
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate 434 1.0500656 1288.0427 7.506976 2000-2012 –> 2013-2019 2.0773265 3626.9446
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 434 1.0605845 1296.6946 7.328068 2000-2012 –> 2013-2019 2.0979999 3651.4134
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 434 1.0174317 1260.6390 9.169936 2000-2012 –> 2013-2019 2.0108645 3545.6973
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 434 1.0645856 1299.9629 7.272475 2000-2012 –> 2013-2019 2.1025508 3655.5339
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 434 1.0701724 1304.5062 7.163185 2000-2012 –> 2013-2019 2.1143246 3669.6335
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 434 1.0517810 1289.4595 7.474466 2000-2012 –> 2013-2019 2.0814698 3632.1573
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 434 1.0832185 1315.0237 7.100796 2000-2012 –> 2013-2019 2.1382142 3696.7638
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 434 1.0672164 1302.1053 7.231208 2000-2012 –> 2013-2019 2.1084324 3662.7045
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 936 0.7619305 2161.2561 6.055788 1987-1999 –> 2000-2012 1.4129829 3990.5005
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 936 0.7598609 2156.1646 5.927820 1987-1999 –> 2000-2012 1.4065541 3973.0879
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate 936 0.7657655 2170.6548 5.940365 1987-1999 –> 2000-2012 1.4217844 4013.8368
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 936 0.7694118 2179.5475 5.861659 1987-1999 –> 2000-2012 1.4297727 4034.8280
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 936 0.7641038 2166.5883 5.828897 1987-1999 –> 2000-2012 1.4161699 3998.6861
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 936 0.7691115 2178.8167 5.746491 1987-1999 –> 2000-2012 1.4283996 4031.1155
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 936 0.7682084 2176.6172 5.890183 1987-1999 –> 2000-2012 1.4283267 4031.2241
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 936 0.7666596 2172.8394 5.792427 1987-1999 –> 2000-2012 1.4228638 4016.5390
Atlantic (28.1,28.15] 5 cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 936 0.8033422 2260.3327 10.678091 1987-1999 –> 2000-2012 1.4323220 4026.3207
Atlantic (28.1,28.15] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 936 0.7856962 2218.7544 10.025704 1987-1999 –> 2000-2012 1.4043107 3954.2671
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 497 0.8713460 1287.5350 6.170119 2000-2012 –> 2013-2019 1.6332764 3448.7912
Atlantic (28.1,28.15] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 497 0.8714406 1285.6430 6.181914 2000-2012 –> 2013-2019 1.6447306 3472.6026
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 497 0.8712585 1287.4353 6.182326 2000-2012 –> 2013-2019 1.6311195 3443.5998
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate 497 0.8772737 1294.2742 6.030801 2000-2012 –> 2013-2019 1.6430391 3464.9290
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 497 0.8847238 1302.6800 5.915366 2000-2012 –> 2013-2019 1.6541356 3482.2275
Atlantic (28.1,28.15] 4 cstar_tref ~ sal + temp + aou + phosphate 497 0.8778716 1292.9515 6.059913 2000-2012 –> 2013-2019 1.6516901 3481.1901
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 497 0.8492568 1262.0116 6.971050 2000-2012 –> 2013-2019 1.6225461 3450.9696
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 497 0.8778689 1294.9484 6.059947 2000-2012 –> 2013-2019 1.6419727 3461.5366
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 497 0.8803306 1297.7319 5.989721 2000-2012 –> 2013-2019 1.6485390 3474.3491
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 497 0.8812870 1298.8112 6.030679 2000-2012 –> 2013-2019 1.6479466 3471.6506
Atlantic (28.15,28.2] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 1540 0.5531775 2560.7353 2.247783 1987-1999 –> 2000-2012 1.1303358 5359.9720
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 1540 0.5735292 2672.0157 1.995751 1987-1999 –> 2000-2012 1.1134623 5257.6373
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate 1540 0.5708735 2657.7209 2.069091 1987-1999 –> 2000-2012 1.1297679 5353.9296
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1540 0.5471615 2527.0561 2.040783 1987-1999 –> 2000-2012 1.0845927 5097.7965
Atlantic (28.15,28.2] 4 cstar_tref ~ sal + temp + aou + nitrate 1540 0.5816010 2713.0610 2.059205 1987-1999 –> 2000-2012 1.1321035 5358.7965
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate 1540 0.5815984 2715.0475 2.055947 1987-1999 –> 2000-2012 1.1270613 5333.3163
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 1540 0.5814763 2714.4005 2.081606 1987-1999 –> 2000-2012 1.1254939 5324.1688
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 1540 0.5693363 2649.4158 2.047829 1987-1999 –> 2000-2012 1.1198001 5296.9262
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 1540 0.5879976 2748.7509 2.066613 1987-1999 –> 2000-2012 1.1421749 5417.8031
Atlantic (28.15,28.2] 5 cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 1540 0.3447362 1104.2051 1.681506 1987-1999 –> 2000-2012 0.7824170 3017.1298
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 819 0.7377046 1839.9223 4.453066 2000-2012 –> 2013-2019 1.3112338 4511.9380
Atlantic (28.15,28.2] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 819 0.7431790 1850.0329 4.262643 2000-2012 –> 2013-2019 1.3371546 4627.9389
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate 819 0.7644207 1898.1939 4.832071 2000-2012 –> 2013-2019 1.3352942 4555.9148
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 819 0.7372762 1838.9708 4.731734 2000-2012 –> 2013-2019 1.2844377 4366.0269
Atlantic (28.15,28.2] 4 cstar_tref ~ sal + temp + aou + nitrate 819 0.7638188 1894.9037 5.281723 2000-2012 –> 2013-2019 1.3454198 4607.9647
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate 819 0.7440089 1853.8609 4.415910 2000-2012 –> 2013-2019 1.3256073 4568.9084
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 819 0.7409497 1847.1121 4.370070 2000-2012 –> 2013-2019 1.3224260 4561.5126
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 819 0.7570242 1882.2675 4.592474 2000-2012 –> 2013-2019 1.3263604 4531.6833
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 819 0.7504528 1867.9868 4.442750 2000-2012 –> 2013-2019 1.3384504 4616.7377
Atlantic (28.15,28.2] 5 cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 819 0.6693303 1680.6011 8.052277 2000-2012 –> 2013-2019 1.0140665 2784.8062
Atlantic (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + phosphate 2693 0.3492775 1988.9312 1.996120 1987-1999 –> 2000-2012 0.6605172 3334.1811
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 2693 0.3469830 1955.4330 2.125336 1987-1999 –> 2000-2012 0.6580261 3299.3355
Atlantic (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 2693 0.3540505 2062.0339 1.937681 1987-1999 –> 2000-2012 0.6652696 3406.9342
Atlantic (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + silicate 2693 0.3588037 2133.8610 1.963165 1987-1999 –> 2000-2012 0.6730852 3530.6199
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate 2693 0.3242828 1591.0173 2.115462 1987-1999 –> 2000-2012 0.6331479 2897.7064
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 2693 0.3463114 1944.9979 1.900529 1987-1999 –> 2000-2012 0.6552064 3252.1991
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate 2693 0.3471287 1957.6942 2.120895 1987-1999 –> 2000-2012 0.6580626 3299.7368
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 2693 0.3473235 1960.7157 2.119236 1987-1999 –> 2000-2012 0.6584686 3306.3555
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 2693 0.3584097 2129.9437 2.006883 1987-1999 –> 2000-2012 0.6719880 3516.8388
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 2693 0.3457583 1936.3893 2.124533 1987-1999 –> 2000-2012 0.6529867 3214.9394
Atlantic (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + phosphate 1557 0.3868691 1473.3053 2.525540 2000-2012 –> 2013-2019 0.7361466 3462.2365
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 1557 0.3797467 1417.4419 2.705317 2000-2012 –> 2013-2019 0.7267298 3372.8749
Atlantic (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 1557 0.3991609 1570.7062 2.411345 2000-2012 –> 2013-2019 0.7532114 3632.7401
Atlantic (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + silicate 1557 0.4110188 1661.8663 2.465092 2000-2012 –> 2013-2019 0.7698225 3795.7273
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate 1557 0.3502557 1165.7023 2.704232 2000-2012 –> 2013-2019 0.6745385 2756.7196
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1557 0.3901368 1501.4981 2.457959 2000-2012 –> 2013-2019 0.7364483 3446.4960
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate 1557 0.3800118 1419.6148 2.702494 2000-2012 –> 2013-2019 0.7271406 3377.3089
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 1557 0.3803071 1422.0335 2.702200 2000-2012 –> 2013-2019 0.7276306 3382.7492
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 1557 0.4078439 1639.7187 2.585072 2000-2012 –> 2013-2019 0.7662535 3769.6624
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 1557 0.3788592 1410.1550 2.700717 2000-2012 –> 2013-2019 0.7246175 3346.5443
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 5288 5.9081446 33807.1795 33.261331 1987-1999 –> 2000-2012 12.1438898 67736.7113
Indo-Pacific (-Inf,26] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 5288 6.2360939 34376.5179 29.968567 1987-1999 –> 2000-2012 12.7645112 68782.8038
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 5288 6.1928478 34304.9200 29.527354 1987-1999 –> 2000-2012 12.6659614 68624.4065
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate 5288 6.0145280 33995.9192 33.811166 1987-1999 –> 2000-2012 12.3488462 68089.1616
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 5288 6.0656694 34085.4664 34.021650 1987-1999 –> 2000-2012 12.4469069 68255.7399
Indo-Pacific (-Inf,26] 4 cstar_tref ~ sal + temp + aou + phosphate 5288 6.3644017 34591.9111 30.317607 1987-1999 –> 2000-2012 13.0100704 69184.0015
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 5288 6.2323348 34372.1408 29.869481 1987-1999 –> 2000-2012 12.7593392 68778.1674
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 5288 6.3160008 34513.1739 29.849776 1987-1999 –> 2000-2012 12.8971080 69005.3651
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5288 6.3788089 34617.8251 29.965279 1987-1999 –> 2000-2012 13.0127648 69193.5027
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 5288 6.0863610 34121.4827 34.020695 1987-1999 –> 2000-2012 12.4880597 68325.1716
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 2816 5.6382388 17746.4099 27.667683 2000-2012 –> 2013-2019 11.5463834 51553.5894
Indo-Pacific (-Inf,26] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 2816 5.9520247 18049.4373 24.894024 2000-2012 –> 2013-2019 12.1881186 52425.9552
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2816 5.8778725 17980.8313 24.128023 2000-2012 –> 2013-2019 12.0707203 52285.7513
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate 2816 5.7297950 17837.1302 28.456699 2000-2012 –> 2013-2019 11.7443230 51833.0493
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 2816 5.7800044 17886.2676 29.028201 2000-2012 –> 2013-2019 11.8456737 51971.7340
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 2816 5.9502121 18049.7219 24.771202 2000-2012 –> 2013-2019 12.1825468 52421.8627
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 2816 5.9876308 18085.0286 24.907000 2000-2012 –> 2013-2019 12.3036316 52598.2025
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 2816 5.7921918 17898.1305 28.600115 2000-2012 –> 2013-2019 11.8785528 52019.6131
Indo-Pacific (-Inf,26] 4 cstar_tref ~ sal + temp + nitrate + phosphate_star 2816 5.9297425 18028.3137 27.542353 2000-2012 –> 2013-2019 12.3132039 52651.8497
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 2816 5.9257831 18026.5518 27.359693 2000-2012 –> 2013-2019 12.3060994 52646.8759
Indo-Pacific (26,26.5] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 5525 4.7575940 32928.4208 48.010441 1987-1999 –> 2000-2012 9.1806115 64911.5999
Indo-Pacific (26,26.5] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 5525 4.5360650 32401.5327 43.772780 1987-1999 –> 2000-2012 8.7333493 63808.3783
Indo-Pacific (26,26.5] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate 5525 4.8200707 33072.5849 52.205037 1987-1999 –> 2000-2012 9.1224501 64751.5152
Indo-Pacific (26,26.5] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 5525 4.8628318 33170.1823 52.418978 1987-1999 –> 2000-2012 9.1928551 64919.5772
Indo-Pacific (26,26.5] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 5525 4.5877220 32526.6595 39.934460 1987-1999 –> 2000-2012 8.8395546 64075.5670
Indo-Pacific (26,26.5] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5525 4.6579284 32694.4781 38.557719 1987-1999 –> 2000-2012 8.9676104 64392.0669
Indo-Pacific (26,26.5] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5525 4.6658266 32713.1992 39.207721 1987-1999 –> 2000-2012 8.9864870 64438.7785
Indo-Pacific (26,26.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate 5525 4.7292783 32862.4581 47.412368 1987-1999 –> 2000-2012 9.0854923 64678.1998
Indo-Pacific (26,26.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5525 4.7555301 32923.6260 49.041150 1987-1999 –> 2000-2012 9.1099591 64734.8587
Indo-Pacific (26,26.5] 5 cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 5525 4.7637711 32942.7584 51.203206 1987-1999 –> 2000-2012 9.1218326 64763.1651
Indo-Pacific (26,26.5] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 2918 4.5881541 17185.9417 50.798645 2000-2012 –> 2013-2019 9.3457481 50114.3625
Indo-Pacific (26,26.5] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2918 4.4724496 17036.8813 46.310491 2000-2012 –> 2013-2019 9.0085146 49438.4140
Indo-Pacific (26,26.5] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate 2918 4.7212174 17352.7865 54.105676 2000-2012 –> 2013-2019 9.5412882 50425.3714
Indo-Pacific (26,26.5] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 2918 4.5594237 17149.2824 39.840874 2000-2012 –> 2013-2019 9.1471456 49675.9420
Indo-Pacific (26,26.5] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2918 4.6381187 17249.1516 36.657339 2000-2012 –> 2013-2019 9.2960471 49943.6297
Indo-Pacific (26,26.5] 5 cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 2918 4.6824845 17304.7104 36.962646 2000-2012 –> 2013-2019 9.5754407 50543.1344
Indo-Pacific (26,26.5] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 2918 4.6426236 17254.8173 36.476302 2000-2012 –> 2013-2019 9.3084502 49968.0165
Indo-Pacific (26,26.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate 2918 4.5489045 17135.8025 48.107779 2000-2012 –> 2013-2019 9.2781828 49998.2607
Indo-Pacific (26,26.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2918 4.5715543 17164.7889 45.715284 2000-2012 –> 2013-2019 9.3270844 50088.4149
Indo-Pacific (26,26.5] 5 cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 2918 4.5779873 17172.9954 45.249093 2000-2012 –> 2013-2019 9.3417584 50115.7538
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 4544 3.8913517 25257.6932 27.853242 1987-1999 –> 2000-2012 7.8243975 50773.9666
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 4544 3.9266514 25339.7617 26.472095 1987-1999 –> 2000-2012 7.7936031 50701.0322
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 4544 3.9265162 25339.4488 20.226531 1987-1999 –> 2000-2012 7.7560685 50611.8329
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate 4544 3.9690317 25437.3229 27.777704 1987-1999 –> 2000-2012 7.8854836 50914.9268
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 4544 4.0209885 25555.5176 28.598406 1987-1999 –> 2000-2012 7.9881840 51150.8611
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 4544 4.0407862 25600.1536 20.538812 1987-1999 –> 2000-2012 7.9731465 51114.8329
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 4544 4.0193476 25551.8084 28.313853 1987-1999 –> 2000-2012 7.9881765 51150.9165
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 4544 3.9296720 25346.7499 27.085865 1987-1999 –> 2000-2012 7.8831072 50910.3151
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate 4544 3.9535838 25401.8824 28.364134 1987-1999 –> 2000-2012 7.9461242 51055.4701
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 4544 3.9477481 25388.4579 23.151958 1987-1999 –> 2000-2012 7.8833195 50910.6027
Indo-Pacific (26.5,26.75] 4 cstar_tref ~ aou + nitrate + phosphate + phosphate_star 2447 4.3286710 14127.2704 31.590767 2000-2012 –> 2013-2019 8.3395286 39657.8620
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 2447 4.2048977 13987.2925 31.207461 2000-2012 –> 2013-2019 8.0962494 39244.9857
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 2447 4.2585434 14049.3348 30.985390 2000-2012 –> 2013-2019 8.1851948 39389.0965
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2447 4.2938981 14089.7973 24.434470 2000-2012 –> 2013-2019 8.2204143 39429.2462
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate 2447 4.3013059 14098.2332 32.116974 2000-2012 –> 2013-2019 8.2703377 39535.5561
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 2447 4.2783866 14072.0860 34.710261 2000-2012 –> 2013-2019 8.2080586 39418.8360
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 2447 4.3165997 14115.6036 30.418179 2000-2012 –> 2013-2019 8.3097998 39608.0975
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate 2447 4.2621159 14053.4387 31.613974 2000-2012 –> 2013-2019 8.2156997 39455.3210
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 2447 4.2923135 14087.9909 22.815590 2000-2012 –> 2013-2019 8.2400615 39476.4489
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 2447 4.3205048 14120.0290 31.479242 2000-2012 –> 2013-2019 8.3423820 39677.5550
Indo-Pacific (26.75,27] 4 cstar_tref ~ aou + nitrate + phosphate + phosphate_star 5354 4.4461238 31182.6795 21.544599 1987-1999 –> 2000-2012 8.4330785 61500.7464
Indo-Pacific (26.75,27] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 5354 4.3860284 31038.9594 20.143089 1987-1999 –> 2000-2012 8.3266577 61232.6165
Indo-Pacific (26.75,27] 5 cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 5354 4.2395987 30675.3639 23.246791 1987-1999 –> 2000-2012 7.9521937 60224.2864
Indo-Pacific (26.75,27] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 5354 4.3077986 30844.2463 21.980567 1987-1999 –> 2000-2012 8.0607655 60508.1504
Indo-Pacific (26.75,27] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 5354 4.2906105 30803.4359 21.052395 1987-1999 –> 2000-2012 8.0273091 60422.3531
Indo-Pacific (26.75,27] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate 5354 4.4804315 31266.9886 23.528787 1987-1999 –> 2000-2012 8.4430071 61520.7158
Indo-Pacific (26.75,27] 5 cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 5354 4.1643941 30483.7138 22.734855 1987-1999 –> 2000-2012 7.8118146 59841.0756
Indo-Pacific (26.75,27] 5 cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 5354 4.1734420 30506.9538 22.914051 1987-1999 –> 2000-2012 7.8834244 60048.2622
Indo-Pacific (26.75,27] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 5354 4.2057274 30587.4712 23.104196 1987-1999 –> 2000-2012 7.9322585 60174.9183
Indo-Pacific (26.75,27] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 5354 4.1702723 30498.8181 21.453131 1987-1999 –> 2000-2012 7.8654423 59996.8565
Indo-Pacific (26.75,27] 4 cstar_tref ~ aou + nitrate + phosphate + phosphate_star 2868 4.7973577 17145.4539 21.705143 2000-2012 –> 2013-2019 9.2434815 48328.1334
Indo-Pacific (26.75,27] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 2868 4.7432090 17082.3425 20.342378 2000-2012 –> 2013-2019 9.1292375 48121.3019
Indo-Pacific (26.75,27] 5 cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 2868 4.5834202 16885.7787 21.941741 2000-2012 –> 2013-2019 8.8230189 47561.1426
Indo-Pacific (26.75,27] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 2868 4.6612052 16980.3073 20.095371 2000-2012 –> 2013-2019 8.9690038 47824.5536
Indo-Pacific (26.75,27] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2868 4.6499950 16968.4956 19.282213 2000-2012 –> 2013-2019 8.9406055 47771.9315
Indo-Pacific (26.75,27] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate 2868 4.8368720 17194.5060 22.935453 2000-2012 –> 2013-2019 9.3173034 48461.4946
Indo-Pacific (26.75,27] 5 cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 2868 4.5360260 16826.1577 17.606576 2000-2012 –> 2013-2019 8.7004201 47309.8716
Indo-Pacific (26.75,27] 5 cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 2868 4.5371637 16827.5961 18.186152 2000-2012 –> 2013-2019 8.7106057 47334.5499
Indo-Pacific (26.75,27] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 2868 4.5774010 16876.2410 16.522089 2000-2012 –> 2013-2019 8.7831284 47463.7121
Indo-Pacific (26.75,27] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 2868 4.5537010 16848.4650 17.066955 2000-2012 –> 2013-2019 8.7239733 47347.2831
Indo-Pacific (27,27.25] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 6466 3.9057692 35982.9778 48.713871 1987-1999 –> 2000-2012 7.1752995 69192.6850
Indo-Pacific (27,27.25] 5 cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 6466 4.3554983 37392.3626 65.794302 1987-1999 –> 2000-2012 7.9401910 71775.4065
Indo-Pacific (27,27.25] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 6466 4.3929768 37503.1646 41.945012 1987-1999 –> 2000-2012 8.0699453 72210.2610
Indo-Pacific (27,27.25] 5 cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 6466 4.2956539 37213.4453 48.286060 1987-1999 –> 2000-2012 7.9243301 71751.9773
Indo-Pacific (27,27.25] 5 cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 6466 4.2296282 37013.1327 54.406272 1987-1999 –> 2000-2012 7.8059700 71366.4393
Indo-Pacific (27,27.25] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate 6466 4.0268865 36377.9050 48.755703 1987-1999 –> 2000-2012 7.3968635 69973.4217
Indo-Pacific (27,27.25] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 6466 4.1535561 36778.4265 49.359080 1987-1999 –> 2000-2012 7.6352783 70789.8628
Indo-Pacific (27,27.25] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 6466 4.3599345 37403.5276 46.003290 1987-1999 –> 2000-2012 8.0771692 72247.4890
Indo-Pacific (27,27.25] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 6466 3.9500266 36128.6900 31.512269 1987-1999 –> 2000-2012 7.3263654 69748.2534
Indo-Pacific (27,27.25] 5 cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 6466 4.1253304 36690.2467 48.506869 1987-1999 –> 2000-2012 7.5924634 70648.1454
Indo-Pacific (27,27.25] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 3470 4.7284767 20643.4390 37.054814 2000-2012 –> 2013-2019 8.6342459 56626.4168
Indo-Pacific (27,27.25] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate 3470 4.9830367 21007.3475 27.071879 2000-2012 –> 2013-2019 9.3334653 58384.6488
Indo-Pacific (27,27.25] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 3470 5.0789520 21139.6617 33.597697 2000-2012 –> 2013-2019 9.4719288 58642.8263
Indo-Pacific (27,27.25] 5 cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 3470 5.0651721 21120.8069 42.030410 2000-2012 –> 2013-2019 9.3608260 58334.2522
Indo-Pacific (27,27.25] 5 cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 3470 5.1142281 21187.6972 42.108737 2000-2012 –> 2013-2019 9.3438563 58200.8298
Indo-Pacific (27,27.25] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate 3470 4.8388654 20803.5946 35.817177 2000-2012 –> 2013-2019 8.8657519 57181.4995
Indo-Pacific (27,27.25] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 3470 4.9716940 20991.5322 35.708885 2000-2012 –> 2013-2019 9.1252501 57769.9588
Indo-Pacific (27,27.25] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 3470 5.1322261 21210.0776 39.784583 2000-2012 –> 2013-2019 9.4921606 58613.6052
Indo-Pacific (27,27.25] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 3470 4.4723380 20256.9380 25.616236 2000-2012 –> 2013-2019 8.4223646 56385.6280
Indo-Pacific (27,27.25] 5 cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 3470 4.9160028 20913.3540 34.970533 2000-2012 –> 2013-2019 9.0413332 57603.6007
Indo-Pacific (27.25,27.5] 4 cstar_tref ~ aou + nitrate + silicate + phosphate 5847 3.6227030 31657.8230 34.491357 1987-1999 –> 2000-2012 6.7579985 61144.3054
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 5847 3.5830901 31531.2492 32.454766 1987-1999 –> 2000-2012 6.6949058 60933.2405
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate 5847 3.4547220 31104.6102 30.937998 1987-1999 –> 2000-2012 6.4725755 60153.8209
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 5847 3.4261572 31007.5182 32.305003 1987-1999 –> 2000-2012 6.4043943 59904.6857
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 5847 3.4222130 30994.0484 31.732058 1987-1999 –> 2000-2012 6.3952098 59870.9527
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 5847 3.6349919 31699.4240 32.897992 1987-1999 –> 2000-2012 6.7367134 61064.0314
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5847 3.6434273 31726.5299 33.075508 1987-1999 –> 2000-2012 6.7582801 61139.7456
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate 5847 3.5979957 31579.7952 32.877171 1987-1999 –> 2000-2012 6.7170617 61008.5633
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 5847 3.3461863 30731.3292 33.390589 1987-1999 –> 2000-2012 6.3351483 59669.8566
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 5847 3.6204935 31652.6884 33.074778 1987-1999 –> 2000-2012 6.7589833 61152.8877
Indo-Pacific (27.25,27.5] 4 cstar_tref ~ aou + nitrate + silicate + phosphate 3127 4.2427270 17924.3614 33.367400 2000-2012 –> 2013-2019 7.8654301 49582.1844
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 3127 4.1827766 17837.3609 30.427009 2000-2012 –> 2013-2019 7.7658666 49368.6101
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate 3127 4.0139968 17579.7723 29.726264 2000-2012 –> 2013-2019 7.4687188 48684.3824
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 3127 3.9776245 17522.8442 31.036835 2000-2012 –> 2013-2019 7.4037817 48530.3624
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 3127 3.9789503 17524.9283 30.123196 2000-2012 –> 2013-2019 7.4011633 48518.9767
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 3127 4.2298938 17907.4158 32.750414 2000-2012 –> 2013-2019 7.8648856 49606.8398
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + temp + nitrate + silicate + phosphate 3127 4.0652536 17659.1273 34.878738 2000-2012 –> 2013-2019 7.6535356 49207.3087
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate 3127 4.2035697 17868.3735 30.955479 2000-2012 –> 2013-2019 7.8015654 49448.1686
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 3127 3.8437243 17308.6883 33.232474 2000-2012 –> 2013-2019 7.1899105 48040.0175
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 3127 4.2264297 17902.2919 31.188098 2000-2012 –> 2013-2019 7.8469231 49554.9803
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 5836 3.2683124 30398.6938 24.119856 1987-1999 –> 2000-2012 5.6276908 57326.2596
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 5836 3.2280103 30253.8699 9.417658 1987-1999 –> 2000-2012 5.6206063 57346.6598
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 5836 3.2980912 30504.5604 26.307453 1987-1999 –> 2000-2012 5.6641455 57465.5180
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 5836 3.1915247 30121.1923 10.602973 1987-1999 –> 2000-2012 5.7370142 57946.0367
Indo-Pacific (27.5,27.75] 4 cstar_tref ~ sal + temp + aou + silicate 5836 3.2235153 30235.6053 10.906227 1987-1999 –> 2000-2012 5.7799780 58109.2855
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 5836 3.1735450 30055.2512 8.689085 1987-1999 –> 2000-2012 5.5355040 56995.7355
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5836 3.1804225 30080.5185 8.700790 1987-1999 –> 2000-2012 5.5552316 57085.1239
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + temp + nitrate + silicate + phosphate 5836 3.1294481 29891.9296 12.455483 1987-1999 –> 2000-2012 5.5232934 56990.8888
Indo-Pacific (27.5,27.75] 4 cstar_tref ~ sal + temp + silicate + phosphate 5836 3.1824417 30085.9266 11.982022 1987-1999 –> 2000-2012 5.5856721 57229.1278
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5836 3.1698959 30041.8225 8.686571 1987-1999 –> 2000-2012 5.5383912 57014.9659
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 3104 3.9417327 17337.7899 25.386634 2000-2012 –> 2013-2019 7.2100451 47736.4837
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 3104 3.7804110 17078.3721 11.020697 2000-2012 –> 2013-2019 7.0084214 47332.2420
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 3104 4.0075075 17440.5265 28.440018 2000-2012 –> 2013-2019 7.3055987 47945.0868
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 3104 3.6736333 16900.5032 10.209719 2000-2012 –> 2013-2019 6.8651581 47021.6955
Indo-Pacific (27.5,27.75] 4 cstar_tref ~ sal + temp + aou + silicate 3104 3.7197021 16975.8697 10.325293 2000-2012 –> 2013-2019 6.9432174 47211.4751
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 3104 3.7040372 16951.6707 10.351471 2000-2012 –> 2013-2019 6.8775822 47006.9218
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 3104 3.7079082 16958.1551 10.354380 2000-2012 –> 2013-2019 6.8883307 47038.6736
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + temp + nitrate + silicate + phosphate 3104 3.7023122 16948.7788 17.138991 2000-2012 –> 2013-2019 6.8317603 46840.7084
Indo-Pacific (27.5,27.75] 4 cstar_tref ~ sal + temp + silicate + phosphate 3104 3.7970417 17103.6222 16.783216 2000-2012 –> 2013-2019 6.9794834 47189.5487
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 3104 3.6928602 16932.9096 10.332716 2000-2012 –> 2013-2019 6.8627561 46974.7321
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 2317 2.7212475 11228.4142 21.917997 1987-1999 –> 2000-2012 4.6423537 20631.9916
Indo-Pacific (27.75,27.85] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 2317 2.7304690 11242.0909 21.046900 1987-1999 –> 2000-2012 4.6699469 20686.8024
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2317 2.7290365 11241.6591 21.325228 1987-1999 –> 2000-2012 4.6624643 20674.2110
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 2317 2.7364990 11254.3133 24.038402 1987-1999 –> 2000-2012 4.6548228 20651.3223
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 2317 2.7349709 11251.7249 22.598195 1987-1999 –> 2000-2012 4.6726675 20694.2721
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate 2317 2.7312774 11245.4627 21.597290 1987-1999 –> 2000-2012 4.6666295 20682.5230
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2317 2.7271269 11238.4154 21.453348 1987-1999 –> 2000-2012 4.6670710 20686.2160
Indo-Pacific (27.75,27.85] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 2317 2.7408402 11259.6589 23.487260 1987-1999 –> 2000-2012 4.6793675 20702.1486
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 2317 2.6406178 11089.0349 20.968772 1987-1999 –> 2000-2012 4.5455194 20454.2227
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 2317 2.7269514 11238.1171 21.350081 1987-1999 –> 2000-2012 4.6695784 20692.1811
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 1264 3.1715991 6518.9850 31.495482 2000-2012 –> 2013-2019 5.8928466 17747.3991
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate 1264 3.1499466 6501.6671 31.458612 2000-2012 –> 2013-2019 5.8945181 17769.6304
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1264 3.1358059 6490.2930 31.989646 2000-2012 –> 2013-2019 5.8629052 17728.6613
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 1264 3.1390093 6492.8741 31.279189 2000-2012 –> 2013-2019 5.8653527 17729.9579
Indo-Pacific (27.75,27.85] 4 cstar_tref ~ sal + nitrate + silicate + phosphate_star 1264 3.1488202 6498.7630 29.448400 2000-2012 –> 2013-2019 5.8762396 17735.6753
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 1264 3.0327595 6405.8242 34.844332 2000-2012 –> 2013-2019 5.7692585 17660.1374
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 1264 3.1487768 6500.7281 29.551108 2000-2012 –> 2013-2019 5.8739654 17735.8488
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 1264 3.1085536 6468.2269 32.890347 2000-2012 –> 2013-2019 5.8435245 17719.9518
Indo-Pacific (27.75,27.85] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 1264 3.1128492 6469.7178 31.766441 2000-2012 –> 2013-2019 5.8536894 17729.3767
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1264 3.0665742 6433.8549 28.861392 2000-2012 –> 2013-2019 5.7071919 17522.8898
Indo-Pacific (27.85,27.95] 4 cstar_tref ~ aou + nitrate + phosphate + phosphate_star 2907 2.7391027 14120.0717 28.616999 1987-1999 –> 2000-2012 4.9389199 26921.2420
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 2907 2.7369084 14117.4123 28.359649 1987-1999 –> 2000-2012 4.9362933 26919.4439
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 2907 2.7387015 14121.2200 28.728209 1987-1999 –> 2000-2012 4.9384200 26924.1304
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate 2907 2.7464377 14137.6201 28.678001 1987-1999 –> 2000-2012 4.9507150 26952.5257
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 2907 2.7459052 14136.4927 28.587098 1987-1999 –> 2000-2012 4.9518292 26955.7252
Indo-Pacific (27.85,27.95] 4 cstar_tref ~ temp + aou + nitrate + phosphate 2907 2.7478439 14138.5961 28.465495 1987-1999 –> 2000-2012 4.9526352 26952.8529
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 2907 2.7026569 14044.1928 28.417644 1987-1999 –> 2000-2012 4.8938578 26824.6247
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate 2907 2.7448988 14134.3614 28.185033 1987-1999 –> 2000-2012 4.9495465 26950.2407
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2907 2.7457460 14136.1557 28.208152 1987-1999 –> 2000-2012 4.9526693 26958.0122
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 2907 2.7445406 14133.6026 28.109511 1987-1999 –> 2000-2012 4.9509877 26954.2090
Indo-Pacific (27.85,27.95] 4 cstar_tref ~ aou + nitrate + phosphate + phosphate_star 1542 3.3028283 8072.7054 32.389813 2000-2012 –> 2013-2019 6.0419310 22192.7771
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 1542 3.2823201 8055.4962 31.823285 2000-2012 –> 2013-2019 6.0192285 22172.9085
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 1542 3.2955181 8067.8720 32.895254 2000-2012 –> 2013-2019 6.0342196 22189.0919
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 1542 3.2978683 8070.0705 32.004378 2000-2012 –> 2013-2019 6.0457481 22210.7427
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 1542 3.3028408 8074.7171 32.542471 2000-2012 –> 2013-2019 6.0487460 22211.2098
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 1542 3.1767686 7954.6926 31.850384 2000-2012 –> 2013-2019 5.8794254 21998.8853
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate 1542 3.2971709 8069.4182 31.533669 2000-2012 –> 2013-2019 6.0420696 22203.7796
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1542 3.2946985 8067.1049 31.554855 2000-2012 –> 2013-2019 6.0404446 22203.2606
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1542 3.2637446 8037.9935 37.426960 2000-2012 –> 2013-2019 6.0494318 22258.1136
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 1542 3.2939987 8066.4498 31.381712 2000-2012 –> 2013-2019 6.0385393 22200.0524
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 2467 2.1861625 10874.1595 7.408724 1987-1999 –> 2000-2012 4.1608748 21582.3057
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2467 2.2115502 10931.1275 8.468880 1987-1999 –> 2000-2012 4.1169999 21457.3941
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate 2467 2.1856536 10873.0109 7.438055 1987-1999 –> 2000-2012 4.1607681 21582.1945
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 2467 2.1868435 10875.6962 8.583791 1987-1999 –> 2000-2012 4.0828446 21376.6405
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 2467 2.2194485 10948.7174 8.518833 1987-1999 –> 2000-2012 4.1317894 21493.3734
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2467 2.2200573 10950.0706 8.511640 1987-1999 –> 2000-2012 4.1339878 21498.9594
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 2467 2.2177921 10945.0336 8.776951 1987-1999 –> 2000-2012 4.1300492 21489.4667
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 2467 2.2201767 10950.3359 8.520844 1987-1999 –> 2000-2012 4.1341116 21499.2364
Indo-Pacific (27.95,28.05] 4 cstar_tref ~ sal + temp + silicate + phosphate_star 2467 2.2207729 10949.6606 8.612052 1987-1999 –> 2000-2012 4.1347078 21496.5612
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 2467 2.1886299 10879.7253 7.878228 1987-1999 –> 2000-2012 4.0668310 21332.6199
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 1354 2.3937485 6220.1919 8.306747 2000-2012 –> 2013-2019 4.5799110 17094.3514
Indo-Pacific (27.95,28.05] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 1354 2.4060015 6232.0181 8.104133 2000-2012 –> 2013-2019 4.6173063 17160.5981
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1354 2.4007452 6228.0956 8.259744 2000-2012 –> 2013-2019 4.5943605 17119.0470
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate 1354 2.3949270 6221.5248 8.259385 2000-2012 –> 2013-2019 4.5805806 17094.5357
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 1354 2.3941382 6220.6327 8.284045 2000-2012 –> 2013-2019 4.5799766 17094.0607
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 1354 2.4833521 6319.7076 8.863057 2000-2012 –> 2013-2019 4.6701956 17195.4038
Indo-Pacific (27.95,28.05] 4 cstar_tref ~ sal + temp + nitrate + phosphate 1354 2.3989137 6224.0289 7.921483 2000-2012 –> 2013-2019 4.6236462 17182.4790
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 1354 2.3949943 6221.6009 8.251733 2000-2012 –> 2013-2019 4.5805940 17094.4901
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ sal + temp + nitrate + silicate + phosphate 1354 2.3891259 6214.9574 8.045546 2000-2012 –> 2013-2019 4.5812165 17102.4783
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1354 2.4266912 6257.2052 8.523746 2000-2012 –> 2013-2019 4.6153211 17136.9304
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 1718 1.8463208 6996.4105 10.225349 1987-1999 –> 2000-2012 3.4685949 14054.4557
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1718 1.8450332 6994.0133 9.467590 1987-1999 –> 2000-2012 3.4644540 14045.5418
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 1718 1.8511670 7005.4174 9.474435 1987-1999 –> 2000-2012 3.4730789 14062.6359
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1718 1.8524675 7007.8305 9.472903 1987-1999 –> 2000-2012 3.4760130 14068.7758
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 1718 1.8567455 7015.7563 9.427843 1987-1999 –> 2000-2012 3.4879817 14094.1964
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 1718 1.8521019 7007.1522 9.467448 1987-1999 –> 2000-2012 3.4752800 14067.2599
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate 1718 1.8576083 7017.3526 10.254032 1987-1999 –> 2000-2012 3.4892097 14096.6216
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1718 1.8612616 7024.1034 10.244079 1987-1999 –> 2000-2012 3.4988977 14117.0393
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1718 1.8032691 6915.3424 9.697645 1987-1999 –> 2000-2012 3.4566627 14043.7291
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 1718 1.8584907 7018.9844 10.244915 1987-1999 –> 2000-2012 3.4914553 14101.3451
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 966 2.0709715 4161.9197 10.526175 2000-2012 –> 2013-2019 3.9172924 11158.3302
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 966 2.0390556 4131.9137 9.747515 2000-2012 –> 2013-2019 3.8840888 11125.9270
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + temp + aou + nitrate + silicate 966 2.0456515 4138.1531 9.661353 2000-2012 –> 2013-2019 3.9261877 11197.6557
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 966 2.0453970 4137.9128 9.764284 2000-2012 –> 2013-2019 3.8965640 11143.3302
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 966 2.0456543 4138.1559 9.760669 2000-2012 –> 2013-2019 3.8981218 11145.9864
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + temp + nitrate + silicate + phosphate 966 2.0483772 4140.7257 9.785200 2000-2012 –> 2013-2019 3.9169687 11178.3340
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 966 2.0462976 4138.7633 9.594209 2000-2012 –> 2013-2019 3.9030431 11154.5196
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 966 2.0451030 4137.6351 9.755580 2000-2012 –> 2013-2019 3.8972049 11144.7873
Indo-Pacific (28.05,28.1] 4 cstar_tref ~ sal + temp + silicate + phosphate_star 966 2.0655686 4154.8728 9.967760 2000-2012 –> 2013-2019 3.9297216 11182.3096
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 966 2.0001288 4094.6740 9.111786 2000-2012 –> 2013-2019 3.8033978 11010.0164
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 14627 1.3445893 50185.4032 13.266137 1987-1999 –> 2000-2012 2.4794754 95679.5203
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate 14627 1.2494600 48038.8292 11.951281 1987-1999 –> 2000-2012 2.3340743 92199.6248
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 14627 1.3133686 49498.1312 9.366747 1987-1999 –> 2000-2012 2.4815263 95842.5854
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 14627 1.3376755 50034.5935 16.199862 1987-1999 –> 2000-2012 2.5307848 97001.0068
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 14627 1.3331929 49936.3988 16.466261 1987-1999 –> 2000-2012 2.5162084 96652.7875
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 14627 1.3376755 50034.5945 16.199920 1987-1999 –> 2000-2012 2.5307756 97000.7777
Indo-Pacific (28.1, Inf] 4 cstar_tref ~ sal + temp + silicate + phosphate_star 14627 1.3376891 50032.8912 16.220724 1987-1999 –> 2000-2012 2.5313778 97011.5894
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate 14627 1.3529440 50366.6138 13.068005 1987-1999 –> 2000-2012 2.4941824 96024.9866
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 14627 1.3692488 50717.0570 13.004762 1987-1999 –> 2000-2012 2.5253684 96756.6773
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 14627 1.3591254 50499.9672 12.907049 1987-1999 –> 2000-2012 2.5070021 96329.0199
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 7804 1.5608153 29109.6035 11.237085 2000-2012 –> 2013-2019 2.9054045 79295.0066
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate 7804 1.4123708 27549.7620 9.529089 2000-2012 –> 2013-2019 2.6618308 75588.5913
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 7804 1.4923380 28409.3615 12.644005 2000-2012 –> 2013-2019 2.8302151 78448.3652
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 7804 1.4569278 28034.5504 8.140875 2000-2012 –> 2013-2019 2.7702964 77532.6816
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate 7804 1.4925685 28411.7727 12.628703 2000-2012 –> 2013-2019 2.8303164 78447.9491
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + temp + aou + silicate + phosphate_star 7804 1.4928857 28415.0892 12.611716 2000-2012 –> 2013-2019 2.8305612 78449.6827
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 7804 1.4959833 28447.4412 13.045112 2000-2012 –> 2013-2019 2.8291763 78383.8400
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 7804 1.4928764 28414.9925 12.611567 2000-2012 –> 2013-2019 2.8305520 78449.5870
Indo-Pacific (28.1, Inf] 4 cstar_tref ~ sal + temp + silicate + phosphate_star 7804 1.4961291 28446.9616 12.985788 2000-2012 –> 2013-2019 2.8338182 78479.8528
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate 7804 1.5700546 29201.7234 11.045333 2000-2012 –> 2013-2019 2.9229986 79568.3372

5.3.2 per fitting unit

lm_best %>%
  select(-n_predictors) %>%
  group_by(basin, gamma_slab, eras) %>% 
  summarise(rmse_sum_mean = mean(rmse_sum),
            ais_sum_mean = mean(aic_sum)) %>% 
  ungroup() %>% 
  kable() %>%
  add_header_above() %>%
  kable_styling() %>%
  scroll_box(width = "100%", height = "400px")
basin gamma_slab eras rmse_sum_mean ais_sum_mean
Atlantic (-Inf,26] 1987-1999 –> 2000-2012 2.3486567 978.1352
Atlantic (-Inf,26] 2000-2012 –> 2013-2019 2.2345131 720.7950
Atlantic (26,26.5] 1987-1999 –> 2000-2012 6.9017999 11404.8321
Atlantic (26,26.5] 2000-2012 –> 2013-2019 6.8820610 8594.2351
Atlantic (26.5,26.75] 1987-1999 –> 2000-2012 6.4054944 14504.1808
Atlantic (26.5,26.75] 2000-2012 –> 2013-2019 6.6570558 10773.2155
Atlantic (26.75,27] 1987-1999 –> 2000-2012 3.8404053 19250.3814
Atlantic (26.75,27] 2000-2012 –> 2013-2019 4.6554210 16350.8662
Atlantic (27,27.25] 1987-1999 –> 2000-2012 3.6550456 16022.3501
Atlantic (27,27.25] 2000-2012 –> 2013-2019 4.2695867 13675.2525
Atlantic (27.25,27.5] 1987-1999 –> 2000-2012 4.4074822 18159.6443
Atlantic (27.25,27.5] 2000-2012 –> 2013-2019 5.2879589 14797.9919
Atlantic (27.5,27.75] 1987-1999 –> 2000-2012 4.6472072 24967.6195
Atlantic (27.5,27.75] 2000-2012 –> 2013-2019 5.3055975 20471.8108
Atlantic (27.75,27.85] 1987-1999 –> 2000-2012 3.2036020 7542.8610
Atlantic (27.75,27.85] 2000-2012 –> 2013-2019 3.3562421 5964.1530
Atlantic (27.85,27.95] 1987-1999 –> 2000-2012 6.1000431 9847.0648
Atlantic (27.85,27.95] 2000-2012 –> 2013-2019 6.7733304 8204.1287
Atlantic (27.95,28.05] 1987-1999 –> 2000-2012 7.6297041 13373.0814
Atlantic (27.95,28.05] 2000-2012 –> 2013-2019 7.5789118 10376.1924
Atlantic (28.05,28.1] 1987-1999 –> 2000-2012 1.9126776 4397.2463
Atlantic (28.05,28.1] 2000-2012 –> 2013-2019 2.0891434 3640.0100
Atlantic (28.1,28.15] 1987-1999 –> 2000-2012 1.4203487 4007.0406
Atlantic (28.1,28.15] 2000-2012 –> 2013-2019 1.6418996 3465.1846
Atlantic (28.15,28.2] 1987-1999 –> 2000-2012 1.0887209 5081.7476
Atlantic (28.15,28.2] 2000-2012 –> 2013-2019 1.2940451 4373.3431
Atlantic (28.2, Inf] 1987-1999 –> 2000-2012 0.6586758 3305.8847
Atlantic (28.2, Inf] 2000-2012 –> 2013-2019 0.7342539 3434.3059
Indo-Pacific (-Inf,26] 1987-1999 –> 2000-2012 12.6537458 68597.5031
Indo-Pacific (-Inf,26] 2000-2012 –> 2013-2019 12.0379254 52240.8483
Indo-Pacific (26,26.5] 1987-1999 –> 2000-2012 9.0340202 64547.3707
Indo-Pacific (26,26.5] 2000-2012 –> 2013-2019 9.3169660 50031.1300
Indo-Pacific (26.5,26.75] 1987-1999 –> 2000-2012 7.8921611 50929.4757
Indo-Pacific (26.5,26.75] 2000-2012 –> 2013-2019 8.2427726 39489.3005
Indo-Pacific (26.75,27] 1987-1999 –> 2000-2012 8.0735951 60546.9981
Indo-Pacific (26.75,27] 2000-2012 –> 2013-2019 8.9340778 47752.3974
Indo-Pacific (27,27.25] 1987-1999 –> 2000-2012 7.6943876 70970.3942
Indo-Pacific (27,27.25] 2000-2012 –> 2013-2019 9.1091183 57774.3266
Indo-Pacific (27.25,27.5] 1987-1999 –> 2000-2012 6.6031271 60604.2090
Indo-Pacific (27.25,27.5] 2000-2012 –> 2013-2019 7.6261781 49054.1831
Indo-Pacific (27.5,27.75] 1987-1999 –> 2000-2012 5.6167527 57350.9601
Indo-Pacific (27.5,27.75] 2000-2012 –> 2013-2019 6.9772353 47229.7568
Indo-Pacific (27.75,27.85] 1987-1999 –> 2000-2012 4.6530421 20655.5891
Indo-Pacific (27.75,27.85] 2000-2012 –> 2013-2019 5.8439492 17707.9529
Indo-Pacific (27.85,27.95] 1987-1999 –> 2000-2012 4.9415874 26931.3007
Indo-Pacific (27.85,27.95] 2000-2012 –> 2013-2019 6.0239784 22184.0822
Indo-Pacific (27.95,28.05] 1987-1999 –> 2000-2012 4.1252964 21480.8752
Indo-Pacific (27.95,28.05] 2000-2012 –> 2013-2019 4.6023108 17127.4375
Indo-Pacific (28.05,28.1] 1987-1999 –> 2000-2012 3.4781628 14075.1600
Indo-Pacific (28.05,28.1] 2000-2012 –> 2013-2019 3.8972591 11144.1196
Indo-Pacific (28.1, Inf] 1987-1999 –> 2000-2012 2.4930776 96049.8576
Indo-Pacific (28.1, Inf] 2000-2012 –> 2013-2019 2.8245169 78264.3893

5.4 Target variable coefficients

A data frame to map the target variable is prepared.

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

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

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

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

lm_best_target <- left_join(lm_best_target, lm_all_fitted_wide)

rm(eras_era, eras_forward, eras_backward,
   lm_all_fitted)

5.5 Plot selected model residuals

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

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

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

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

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

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

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

5.6 Cant coeffcients

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

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

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

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

5.7 Write files

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

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

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

6 Model diagnotics

6.1 Selection criterion vs predictors

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

6.3 Predictor counts

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

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

# print table
lm_all_stats %>%
  gt(rowname_col = "gamma_slab",
     groupname_col = c("basin", "eras")) %>% 
  summary_rows(
    groups = TRUE,
    fns = list(total = "sum")
  )
aou nitrate phosphate phosphate_star sal silicate temp
Atlantic - 1987-1999 --> 2000-2012
(-Inf,26] 9 4 8 7 10 3 7
(26,26.5] 8 6 6 6 10 5 8
(26.5,26.75] 9 3 9 6 6 10 5
(26.75,27] 9 3 8 7 4 10 7
(27,27.25] 8 5 8 6 6 10 6
(27.25,27.5] 9 8 6 5 5 10 6
(27.5,27.75] 7 4 7 6 10 8 6
(27.75,27.85] 8 7 8 6 10 4 6
(27.85,27.95] 10 4 10 10 4 4 4
(27.95,28.05] 6 10 8 6 5 8 6
(28.05,28.1] 8 5 9 8 8 4 8
(28.1,28.15] 8 5 8 8 8 5 8
(28.15,28.2] 9 10 6 6 8 4 6
(28.2, Inf] 9 10 5 5 10 4 4
total 117.00 84.00 106.00 92.00 104.00 89.00 87.00
Atlantic - 2000-2012 --> 2013-2019
(-Inf,26] 8 5 7 7 10 5 8
(26,26.5] 8 6 6 6 10 5 8
(26.5,26.75] 9 4 8 6 6 10 6
(26.75,27] 9 3 9 6 6 10 5
(27,27.25] 8 4 6 7 9 10 4
(27.25,27.5] 9 10 4 5 4 10 5
(27.5,27.75] 7 4 7 6 10 8 6
(27.75,27.85] 7 6 8 6 10 6 6
(27.85,27.95] 10 4 9 8 7 3 7
(27.95,28.05] 6 10 8 7 5 7 7
(28.05,28.1] 7 5 7 8 10 5 8
(28.1,28.15] 8 4 9 7 10 3 7
(28.15,28.2] 9 10 5 6 9 3 6
(28.2, Inf] 9 10 5 5 10 4 4
total 114.00 85.00 98.00 90.00 116.00 89.00 87.00
Indo-Pacific - 1987-1999 --> 2000-2012
(-Inf,26] 9 4 8 7 10 3 7
(26,26.5] 8 6 7 7 6 8 8
(26.5,26.75] 8 7 8 7 7 6 7
(26.75,27] 10 5 10 9 5 3 5
(27,27.25] 9 6 9 9 3 6 7
(27.25,27.5] 7 7 9 6 5 10 5
(27.5,27.75] 6 4 6 5 10 10 7
(27.75,27.85] 9 5 9 9 2 7 7
(27.85,27.95] 8 10 9 7 3 4 7
(27.95,28.05] 7 4 6 7 9 8 8
(28.05,28.1] 7 5 7 8 5 10 8
(28.1, Inf] 5 7 6 8 6 10 7
total 93.00 70.00 94.00 89.00 71.00 85.00 83.00
Indo-Pacific - 2000-2012 --> 2013-2019
(-Inf,26] 7 6 7 8 10 3 7
(26,26.5] 7 6 7 7 6 9 8
(26.5,26.75] 8 8 9 8 4 6 6
(26.75,27] 10 5 10 9 5 3 5
(27,27.25] 9 6 9 8 3 7 7
(27.25,27.5] 7 8 9 5 5 10 5
(27.5,27.75] 6 4 6 5 10 10 7
(27.75,27.85] 7 7 7 9 6 7 5
(27.85,27.95] 8 9 8 8 3 6 7
(27.95,28.05] 7 9 6 6 9 4 7
(28.05,28.1] 6 4 6 7 8 10 8
(28.1, Inf] 6 5 7 7 8 10 6
total 88.00 77.00 91.00 87.00 77.00 85.00 78.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
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
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] gt_0.2.2         corrr_0.4.3      broom_0.7.5      kableExtra_1.3.1
 [5] knitr_1.30       olsrr_0.5.3      GGally_2.0.0     lubridate_1.7.9 
 [9] metR_0.9.0       scico_1.2.0      patchwork_1.1.1  collapse_1.5.0  
[13] forcats_0.5.0    stringr_1.4.0    dplyr_1.0.2      purrr_0.3.4     
[17] readr_1.4.0      tidyr_1.1.2      tibble_3.0.4     ggplot2_3.3.3   
[21] tidyverse_1.3.0  workflowr_1.6.2 

loaded via a namespace (and not attached):
 [1] fs_1.5.0                 webshot_0.5.2            RColorBrewer_1.1-2      
 [4] httr_1.4.2               rprojroot_2.0.2          tools_4.0.3             
 [7] backports_1.1.10         R6_2.5.0                 nortest_1.0-4           
[10] DBI_1.1.0                colorspace_2.0-0         withr_2.3.0             
[13] gridExtra_2.3            tidyselect_1.1.0         curl_4.3                
[16] compiler_4.0.3           git2r_0.27.1             cli_2.2.0               
[19] rvest_0.3.6              xml2_1.3.2               sass_0.2.0              
[22] labeling_0.4.2           scales_1.1.1             checkmate_2.0.0         
[25] goftest_1.2-2            digest_0.6.27            foreign_0.8-80          
[28] rmarkdown_2.5            rio_0.5.16               pkgconfig_2.0.3         
[31] htmltools_0.5.0          highr_0.8                dbplyr_1.4.4            
[34] rlang_0.4.10             readxl_1.3.1             rstudioapi_0.13         
[37] farver_2.0.3             generics_0.1.0           jsonlite_1.7.2          
[40] zip_2.1.1                car_3.0-10               magrittr_2.0.1          
[43] Matrix_1.2-18            Rcpp_1.0.5               munsell_0.5.0           
[46] fansi_0.4.1              abind_1.4-5              lifecycle_0.2.0         
[49] stringi_1.5.3            whisker_0.4              yaml_2.2.1              
[52] carData_3.0-4            plyr_1.8.6               grid_4.0.3              
[55] blob_1.2.1               parallel_4.0.3           promises_1.1.1          
[58] crayon_1.3.4             lattice_0.20-41          haven_2.3.1             
[61] hms_0.5.3                pillar_1.4.7             reprex_0.3.0            
[64] glue_1.4.2               evaluate_0.14            RcppArmadillo_0.10.1.2.2
[67] data.table_1.13.6        modelr_0.1.8             vctrs_0.3.6             
[70] httpuv_1.5.4             cellranger_1.1.0         gtable_0.3.0            
[73] reshape_0.8.8            assertthat_0.2.1         xfun_0.20               
[76] openxlsx_4.2.3           RcppEigen_0.3.3.9.1      later_1.1.0.1           
[79] viridisLite_0.3.0        ellipsis_0.3.1           here_1.0.1