Last updated: 2021-03-10

Checks: 7 0

Knit directory: emlr_mod_v_XXX/

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

R Markdown file: up-to-date

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

Environment: empty

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

Seed: set.seed(20200707)

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

Session information: recorded

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

Cache: none

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

File paths: relative

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

Repository version: d12b395

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

1 Required data

Required are:

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

2 Predictor combinations

Find all possible combinations of following considered predictor variables:

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

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

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

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

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

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

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

# convert to tibble
lm_all <- as_tibble(lm_all)

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

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

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

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

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

lm_rm <- lm_rm_ph %>%
  select(-lm_coeff_filter) %>%
  unique()

lm_all <- anti_join(lm_all, lm_rm)

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

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

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

3 Apply predictor threshold

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

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

This results in a total number of MLR models of:

  • 120

4 Fit all models

Individual linear regression models were fitted for the chosen target variable:

  • cstar_tref

as a function of each predictor combination. Fitting was performed separately within each basin, era, and slab. Model diagnostics, such as the root mean squared error (RMSE), were calculated for each fitted model.

# loop across all basins, era, gamma slabs, and MLRs
# fit all MLR models
for (i_basin in unique(GLODAP$basin)) {
  for (i_era in unique(GLODAP$era)) {
    #i_basin <- unique(GLODAP$basin)[1]
    #i_era   <- unique(GLODAP$era)[1]
    print(i_basin)
    print(i_era)
    
    GLODAP_basin_era <- GLODAP %>%
      filter(basin == i_basin, era == i_era)
    
    for (i_gamma_slab in unique(GLODAP_basin_era$gamma_slab)) {
      #i_gamma_slab <- unique(GLODAP_basin_era$gamma_slab)[1]
      print(i_gamma_slab)
      
      GLODAP_basin_era_slab <- GLODAP_basin_era %>%
        filter(gamma_slab == i_gamma_slab)
      
      # number of observations used for each fitting model
      i_nr_obs = nrow(GLODAP_basin_era_slab)
      
      for (i_predictors in unique(lm_all$predictors)) {
        #i_predictors <- unique(lm_all$predictors)[1]
        
        # extract one model definition
        i_lm <- lm_all %>%
          filter(predictors == i_predictors) %>%
          select(lm_coeff) %>%
          pull()
        
        # extract number of predictors
        i_n_predictors <- lm_all %>%
          filter(predictors == i_predictors) %>%
          select(n) %>%
          pull()
        
        if (i_nr_obs > i_n_predictors) {
          # fit model
          if (params_local$MLR_type == "rlm") {
            i_lm_fit <- MASS::rlm(as.formula(i_lm),
                                  data = GLODAP_basin_era_slab)
          }
          
          if (params_local$MLR_type == "lm") {
            i_lm_fit <- lm(as.formula(i_lm),
                           data = GLODAP_basin_era_slab)
          }
          
          # find max predictor correlation
          i_cor_max <- GLODAP_basin_era_slab %>%
            select(!!!syms(str_split(i_predictors, " ",
                                     simplify = TRUE))) %>%
            correlate(quiet = TRUE) %>%
            select(-term) %>%
            abs() %>%
            max(na.rm = TRUE)
          
          # calculate root mean squared error
          i_rmse <- sqrt(c(crossprod(i_lm_fit$residuals)) /
                           length(i_lm_fit$residuals))
          
          # calculate maximum residual
          i_resid_max <- max(abs(i_lm_fit$residuals))
          
          # calculate Akaike information criterion aic
          i_aic <- AIC(i_lm_fit)
          
          # collect model coefficients and diagnostics
          coefficients <- tidy(i_lm_fit)
          
          coefficients <- coefficients %>%
            mutate(
              basin = i_basin,
              era = i_era,
              gamma_slab = i_gamma_slab,
              model = i_lm,
              nr_obs = i_nr_obs,
              rmse = i_rmse,
              aic = i_aic,
              resid_max = i_resid_max,
              n_predictors = i_n_predictors,
              na_predictor = anyNA(coefficients$estimate),
              cor_max = i_cor_max
            )
          
          if (exists("lm_all_fitted")) {
            lm_all_fitted <- bind_rows(lm_all_fitted, coefficients)
          }
          
          if (!exists("lm_all_fitted")) {
            lm_all_fitted <- coefficients
          }
        }
      }
    }
  }
}

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

5 Prepare coeffcients

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

5.1 Formatting

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

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

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

5.2 Predictor selection

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

  • 5

The criterion used to select the best models was:

  • rmse

The criterion was summed up for two adjacent eras, and the models with lowest summed values were selected.

Please note, that currently the lm() function produces NAs for some predictors. It is not yet entirely clear when this happens, but presumably it is caused by some form of collinearity between predictors, such that including another predictor does not help to explain the target variable any better. The issues also expresses as exactly identical rmse values of different models. As an interim solution, models with fitted NA predictors were not included.

# remove models with predictors fitted as NA
lm_all_fitted_wide <- lm_all_fitted_wide %>%
  filter(na_predictor == FALSE)
# calculate RMSE sum for adjacent eras
lm_all_fitted_wide_eras <- lm_all_fitted_wide %>%
  select(basin, gamma_slab, n_predictors, model, era, nr_obs, rmse, aic, resid_max) %>%
  arrange(era) %>%
  group_by(basin, gamma_slab, model) %>%
  mutate(
    eras = paste(lag(era), era, sep = " --> "),
    rmse_sum = rmse + lag(rmse),
    aic_sum = aic + lag(aic)
  ) %>%
  ungroup() %>%
  select(-c(era)) %>%
  drop_na()

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

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

5.3 RMSE tables

5.3.1 per model

lm_best %>%
  kable() %>%
  add_header_above() %>%
  kable_styling() %>%
  scroll_box(width = "100%", height = "400px")
basin gamma_slab n_predictors model nr_obs rmse aic resid_max eras rmse_sum aic_sum
Atlantic (-Inf,26] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 158 1.0102185 467.597241 3.3397502 1982-1999 –> 2000-2012 2.3458438 1194.2659
Atlantic (-Inf,26] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 158 1.0146566 468.982441 3.3276683 1982-1999 –> 2000-2012 2.3517436 1196.1061
Atlantic (-Inf,26] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 158 1.0337825 474.883476 3.6842107 1982-1999 –> 2000-2012 2.3695815 1201.6062
Atlantic (-Inf,26] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 158 1.0093505 469.325595 3.3926925 1982-1999 –> 2000-2012 2.3449212 1197.9773
Atlantic (-Inf,26] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 158 1.0170935 469.740462 3.4137854 1982-1999 –> 2000-2012 2.3550073 1197.1213
Atlantic (-Inf,26] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 72 1.0834728 231.871842 2.6507960 2000-2012 –> 2013-2019 2.0936914 699.4691
Atlantic (-Inf,26] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 72 1.0780946 231.155259 2.5660869 2000-2012 –> 2013-2019 2.0927512 700.1377
Atlantic (-Inf,26] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 72 1.1315155 238.119480 2.9046838 2000-2012 –> 2013-2019 2.1652979 713.0030
Atlantic (-Inf,26] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 72 1.0760286 232.879045 2.6532792 2000-2012 –> 2013-2019 2.0853791 702.2046
Atlantic (-Inf,26] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 72 1.0766499 230.962161 2.5864506 2000-2012 –> 2013-2019 2.0937433 700.7026
Atlantic (26,26.5] 5 cstar_tref ~ sal + oxygen + aou + nitrate + silicate 1040 3.2353488 5407.596546 14.5591803 1982-1999 –> 2000-2012 6.4866697 13277.9482
Atlantic (26,26.5] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 1040 3.2323733 5407.682748 14.8452854 1982-1999 –> 2000-2012 6.4813028 13277.8094
Atlantic (26,26.5] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 1040 3.2672168 5429.984224 12.8619385 1982-1999 –> 2000-2012 6.4813008 13267.5026
Atlantic (26,26.5] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 1040 3.1872705 5378.455168 13.6598787 1982-1999 –> 2000-2012 6.3907896 13206.0170
Atlantic (26,26.5] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 1040 3.1869539 5380.248550 13.7664353 1982-1999 –> 2000-2012 6.3765469 13196.6361
Atlantic (26,26.5] 5 cstar_tref ~ sal + oxygen + aou + nitrate + silicate 562 3.0957021 2879.023463 9.4138668 2000-2012 –> 2013-2019 6.3310508 8286.6200
Atlantic (26,26.5] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 562 3.0352995 2858.875478 9.1979710 2000-2012 –> 2013-2019 6.2676728 8266.5582
Atlantic (26,26.5] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 562 3.1001098 2882.622708 8.8926345 2000-2012 –> 2013-2019 6.4061789 8337.1955
Atlantic (26,26.5] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 562 3.0399853 2860.609307 10.8756225 2000-2012 –> 2013-2019 6.2272558 8239.0645
Atlantic (26,26.5] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 562 3.0022601 2848.573600 10.3922038 2000-2012 –> 2013-2019 6.1892141 8228.8221
Atlantic (26.5,26.75] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 1360 3.2193481 7055.679376 9.2926631 1982-1999 –> 2000-2012 6.2440270 16351.2609
Atlantic (26.5,26.75] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 1360 3.1419857 6991.518410 8.7432236 1982-1999 –> 2000-2012 6.0838053 16187.0486
Atlantic (26.5,26.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1360 3.2276124 7062.652885 10.1005384 1982-1999 –> 2000-2012 6.2017135 16296.2793
Atlantic (26.5,26.75] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 1360 3.1421405 6989.652422 8.8940091 1982-1999 –> 2000-2012 6.0852724 16184.8210
Atlantic (26.5,26.75] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 1360 3.2415884 7072.405398 10.0491820 1982-1999 –> 2000-2012 6.2239552 16314.2286
Atlantic (26.5,26.75] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 646 3.3037649 3393.289581 9.7111224 2000-2012 –> 2013-2019 6.5231130 10448.9690
Atlantic (26.5,26.75] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 646 3.2593836 3377.815822 9.1602872 2000-2012 –> 2013-2019 6.4013694 10369.3342
Atlantic (26.5,26.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 646 3.2958735 3390.199829 10.5343935 2000-2012 –> 2013-2019 6.5234860 10452.8527
Atlantic (26.5,26.75] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 646 3.2640580 3377.667400 9.1208579 2000-2012 –> 2013-2019 6.4061986 10367.3198
Atlantic (26.5,26.75] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 646 3.3285183 3400.933773 11.1260032 2000-2012 –> 2013-2019 6.5701066 10473.3392
Atlantic (26.75,27] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 2134 1.9902912 9009.612817 13.1902563 1982-1999 –> 2000-2012 3.6193377 19426.0277
Atlantic (26.75,27] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 2134 1.9598014 8945.724193 12.8450595 1982-1999 –> 2000-2012 3.4255096 18787.8892
Atlantic (26.75,27] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 2134 1.9756584 8978.118057 13.5842460 1982-1999 –> 2000-2012 3.4606942 18889.7325
Atlantic (26.75,27] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 2134 2.0962710 9231.032860 17.5357342 1982-1999 –> 2000-2012 3.5641818 19079.3875
Atlantic (26.75,27] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 2134 2.1557967 9348.538081 19.2601021 1982-1999 –> 2000-2012 3.6479228 19284.1304
Atlantic (26.75,27] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 1121 2.5305889 5278.849696 22.1725059 2000-2012 –> 2013-2019 4.5208801 14288.4625
Atlantic (26.75,27] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 1121 2.2297870 4997.133558 19.1589435 2000-2012 –> 2013-2019 4.1895884 13942.8578
Atlantic (26.75,27] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1121 2.2698743 5035.082409 21.4501194 2000-2012 –> 2013-2019 4.2455326 14013.2005
Atlantic (26.75,27] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 1121 2.2490469 5014.415793 19.0619091 2000-2012 –> 2013-2019 4.3453179 14245.4487
Atlantic (26.75,27] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 1121 2.3110503 5073.388314 22.0043753 2000-2012 –> 2013-2019 4.4668470 14421.9264
Atlantic (27,27.25] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 1717 1.8684401 7035.241779 13.3529386 1982-1999 –> 2000-2012 3.3709959 15243.1534
Atlantic (27,27.25] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 1717 1.8968952 7087.145042 12.2722372 1982-1999 –> 2000-2012 3.3784479 15231.9081
Atlantic (27,27.25] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 1717 1.8340106 6973.373554 16.0425252 1982-1999 –> 2000-2012 3.3034158 15083.2036
Atlantic (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1717 1.9602218 7199.914902 13.1630682 1982-1999 –> 2000-2012 3.5598163 15688.5723
Atlantic (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 1717 1.9121553 7114.660302 11.6847551 1982-1999 –> 2000-2012 3.4199590 15338.2129
Atlantic (27,27.25] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 982 2.2533101 4398.349480 9.8669171 2000-2012 –> 2013-2019 4.1217502 11433.5913
Atlantic (27,27.25] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 982 2.2697140 4412.595423 9.7561790 2000-2012 –> 2013-2019 4.1666092 11499.7405
Atlantic (27,27.25] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 982 2.1901333 4344.497459 9.6942991 2000-2012 –> 2013-2019 4.0241439 11317.8710
Atlantic (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 982 2.3482249 4479.382903 14.3974130 2000-2012 –> 2013-2019 4.3084468 11679.2978
Atlantic (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 982 2.2871678 4427.640559 9.7187497 2000-2012 –> 2013-2019 4.1993231 11542.3009
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 1621 2.0062322 6871.468535 13.6015772 1982-1999 –> 2000-2012 3.6434264 15485.2959
Atlantic (27.25,27.5] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 1621 2.0044104 6870.523344 13.6394556 1982-1999 –> 2000-2012 3.6274420 15447.2711
Atlantic (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 1621 2.0007200 6864.548797 11.1046357 1982-1999 –> 2000-2012 3.6266986 15449.4563
Atlantic (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 1621 2.0072942 6875.184353 10.8817923 1982-1999 –> 2000-2012 3.6512482 15509.5451
Atlantic (27.25,27.5] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 1621 1.9987512 6863.356999 11.1114802 1982-1999 –> 2000-2012 3.6107239 15411.3516
Atlantic (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 904 2.4549943 4205.249886 14.1424201 2000-2012 –> 2013-2019 4.4557143 11069.7987
Atlantic (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 904 2.4552021 4205.402918 13.6176582 2000-2012 –> 2013-2019 4.4624964 11080.5873
Atlantic (27.25,27.5] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 904 2.4545073 4206.891136 14.0278754 2000-2012 –> 2013-2019 4.4532585 11070.2481
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + temp + oxygen + nitrate + phosphate_star 904 2.4610657 4207.715646 14.4768702 2000-2012 –> 2013-2019 4.4767074 11094.3541
Atlantic (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 904 2.4606156 4209.384994 14.3676036 2000-2012 –> 2013-2019 4.4744249 11095.0749
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + temp + oxygen + aou + nitrate 2079 2.1493013 9095.410309 13.2353104 1982-1999 –> 2000-2012 4.0590191 20641.2306
Atlantic (27.5,27.75] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 2079 2.1479725 9094.838855 13.4747665 1982-1999 –> 2000-2012 4.0573653 20641.7090
Atlantic (27.5,27.75] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 2079 2.1421629 9083.577489 13.9206011 1982-1999 –> 2000-2012 4.0510182 20628.8763
Atlantic (27.5,27.75] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 2079 2.1398842 9081.151964 14.2861377 1982-1999 –> 2000-2012 4.0485002 20627.7510
Atlantic (27.5,27.75] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 2079 2.3001465 9379.447344 16.4944804 1982-1999 –> 2000-2012 4.3698246 21376.2706
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + temp + oxygen + aou + nitrate 1115 2.5428241 5259.436901 15.6212383 2000-2012 –> 2013-2019 4.6921255 14354.8472
Atlantic (27.5,27.75] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 1115 2.5196651 5241.033868 16.4283287 2000-2012 –> 2013-2019 4.6676376 14335.8727
Atlantic (27.5,27.75] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 1115 2.5158530 5237.657523 16.3872161 2000-2012 –> 2013-2019 4.6580160 14321.2350
Atlantic (27.5,27.75] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 1115 2.4881845 5214.996829 17.5841716 2000-2012 –> 2013-2019 4.6280686 14296.1488
Atlantic (27.5,27.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1115 2.5622095 5278.372978 18.9011113 2000-2012 –> 2013-2019 4.8689954 14669.8052
Atlantic (27.75,27.85] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 722 1.6263944 2767.259069 14.3703602 1982-1999 –> 2000-2012 3.2794688 6464.9975
Atlantic (27.75,27.85] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 722 1.5992190 2742.927442 12.1298694 1982-1999 –> 2000-2012 3.2742684 6465.9682
Atlantic (27.75,27.85] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 722 1.5915372 2737.974535 13.3977571 1982-1999 –> 2000-2012 3.2283323 6418.7508
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 722 1.5992193 2740.927704 12.1339144 1982-1999 –> 2000-2012 3.2778246 6466.0316
Atlantic (27.75,27.85] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 722 1.5916813 2736.105232 13.2839690 1982-1999 –> 2000-2012 3.2452848 6434.4569
Atlantic (27.75,27.85] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 364 1.6780436 1425.820870 10.6323642 2000-2012 –> 2013-2019 3.2772626 4168.7483
Atlantic (27.75,27.85] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 364 1.7133169 1440.965149 11.0768950 2000-2012 –> 2013-2019 3.3409938 4209.3625
Atlantic (27.75,27.85] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 364 1.6644845 1421.914514 10.8788216 2000-2012 –> 2013-2019 3.2560217 4159.8890
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 364 1.6780668 1423.830913 10.6183742 2000-2012 –> 2013-2019 3.2772860 4164.7586
Atlantic (27.75,27.85] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 364 1.6649520 1420.118948 10.6920877 2000-2012 –> 2013-2019 3.2566332 4156.2242
Atlantic (27.85,27.95] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 694 2.7700878 3399.690774 20.0287808 1982-1999 –> 2000-2012 5.6973757 8212.2514
Atlantic (27.85,27.95] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 694 2.6216352 3325.238651 18.8864708 1982-1999 –> 2000-2012 5.5488708 8139.7650
Atlantic (27.85,27.95] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 694 2.6858576 3356.830803 20.1208556 1982-1999 –> 2000-2012 5.6386261 8186.0664
Atlantic (27.85,27.95] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 694 2.6286282 3326.936087 18.9155308 1982-1999 –> 2000-2012 5.5560063 8139.5560
Atlantic (27.85,27.95] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 694 2.6878487 3355.859387 20.1028452 1982-1999 –> 2000-2012 5.6412180 8183.4865
Atlantic (27.85,27.95] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 415 2.9980518 2105.028014 17.1978515 2000-2012 –> 2013-2019 5.7681397 5504.7188
Atlantic (27.85,27.95] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 415 2.9854323 2103.526959 17.0507619 2000-2012 –> 2013-2019 5.6070675 5428.7656
Atlantic (27.85,27.95] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 415 3.0207079 2111.276685 17.8599673 2000-2012 –> 2013-2019 5.7065655 5468.1075
Atlantic (27.85,27.95] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 415 2.9858817 2101.651898 17.0435155 2000-2012 –> 2013-2019 5.6145099 5428.5880
Atlantic (27.85,27.95] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 415 3.0207158 2109.278858 17.8578045 2000-2012 –> 2013-2019 5.7085645 5465.1382
Atlantic (27.95,28.05] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 723 2.5486408 3420.605178 13.5656845 1982-1999 –> 2000-2012 4.7958426 7862.6531
Atlantic (27.95,28.05] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 723 2.3956083 3331.064462 11.5680140 1982-1999 –> 2000-2012 4.6539653 7782.9466
Atlantic (27.95,28.05] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 723 2.3733651 3319.575685 12.0861046 1982-1999 –> 2000-2012 4.5747223 7722.6887
Atlantic (27.95,28.05] 5 cstar_tref ~ temp + oxygen + aou + nitrate + phosphate_star 723 2.4844583 3381.724146 12.6800624 1982-1999 –> 2000-2012 4.7616370 7848.0894
Atlantic (27.95,28.05] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 723 2.4691283 3374.774164 13.1508492 1982-1999 –> 2000-2012 4.7055423 7807.2653
Atlantic (27.95,28.05] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 386 1.9182619 1614.316404 15.5159025 2000-2012 –> 2013-2019 4.4669027 5034.9216
Atlantic (27.95,28.05] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 386 1.8983477 1606.260107 14.3241618 2000-2012 –> 2013-2019 4.2939560 4937.3246
Atlantic (27.95,28.05] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 386 1.8741838 1598.370322 14.5567935 2000-2012 –> 2013-2019 4.2475490 4917.9460
Atlantic (27.95,28.05] 5 cstar_tref ~ temp + oxygen + aou + nitrate + phosphate_star 386 1.9421109 1621.855231 14.3765572 2000-2012 –> 2013-2019 4.4265693 5003.5794
Atlantic (27.95,28.05] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 386 1.9131818 1612.269225 14.9940476 2000-2012 –> 2013-2019 4.3823101 4987.0434
Atlantic (28.05,28.1] 5 cstar_tref ~ oxygen + aou + nitrate + silicate + phosphate_star 436 0.2006132 -149.446120 1.1479509 1982-1999 –> 2000-2012 0.3133988 -1046.8590
Atlantic (28.05,28.1] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 436 0.1883627 -202.389941 1.2633372 1982-1999 –> 2000-2012 0.2993422 -1117.0789
Atlantic (28.05,28.1] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 436 0.1727137 -276.022281 1.4177362 1982-1999 –> 2000-2012 0.2836425 -1189.2566
Atlantic (28.05,28.1] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 436 0.1931541 -180.486583 1.6293925 1982-1999 –> 2000-2012 0.3078192 -1056.1670
Atlantic (28.05,28.1] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 436 0.1987386 -155.632485 1.1631731 1982-1999 –> 2000-2012 0.3104055 -1062.9485
Atlantic (28.05,28.1] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 229 0.2249240 -17.458790 0.8208726 2000-2012 –> 2013-2019 0.4132868 -219.8487
Atlantic (28.05,28.1] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 229 0.2357929 4.154713 1.0833354 2000-2012 –> 2013-2019 0.4370684 -140.4168
Atlantic (28.05,28.1] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 229 0.2037662 -60.704379 0.8292810 2000-2012 –> 2013-2019 0.3764799 -336.7267
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 229 0.2362004 2.945556 1.0819569 2000-2012 –> 2013-2019 0.4378010 -142.2187
Atlantic (28.05,28.1] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 229 0.2182829 -31.185279 0.9320609 2000-2012 –> 2013-2019 0.4114370 -211.6719
Atlantic (28.1,28.15] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 541 0.1345088 -619.336477 1.2143151 1982-1999 –> 2000-2012 0.2509305 -1735.8502
Atlantic (28.1,28.15] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 541 0.1382898 -589.341173 1.2495530 1982-1999 –> 2000-2012 0.2572541 -1672.4118
Atlantic (28.1,28.15] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 541 0.1303683 -653.166270 1.1005939 1982-1999 –> 2000-2012 0.2531673 -1687.1256
Atlantic (28.1,28.15] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 541 0.1281185 -670.001838 1.1338428 1982-1999 –> 2000-2012 0.2410768 -1831.2647
Atlantic (28.1,28.15] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 541 0.1387514 -585.735481 1.2504771 1982-1999 –> 2000-2012 0.2582470 -1661.9077
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + oxygen + aou + nitrate + phosphate_star 288 0.1316401 -336.637307 0.5482708 2000-2012 –> 2013-2019 0.2670058 -951.1023
Atlantic (28.1,28.15] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 288 0.1271996 -354.401902 0.4635671 2000-2012 –> 2013-2019 0.2617084 -973.7384
Atlantic (28.1,28.15] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 288 0.1287069 -347.616789 0.5801724 2000-2012 –> 2013-2019 0.2590752 -1000.7831
Atlantic (28.1,28.15] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 288 0.1219786 -376.543574 0.4857358 2000-2012 –> 2013-2019 0.2500970 -1046.5454
Atlantic (28.1,28.15] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 288 0.1247025 -365.822431 0.4871925 2000-2012 –> 2013-2019 0.2601790 -977.4022
Atlantic (28.15,28.2] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 1155 0.1331367 -1364.086548 0.6733656 1982-1999 –> 2000-2012 0.3074876 -2484.0577
Atlantic (28.15,28.2] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 1155 0.1330873 -1364.944479 0.6730089 1982-1999 –> 2000-2012 0.3074425 -2484.8298
Atlantic (28.15,28.2] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 1155 0.1313965 -1392.480390 0.6733169 1982-1999 –> 2000-2012 0.3048434 -2528.4671
Atlantic (28.15,28.2] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 1155 0.1334071 -1359.399767 0.6730578 1982-1999 –> 2000-2012 0.3080253 -2474.0611
Atlantic (28.15,28.2] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 1155 0.1319719 -1384.386811 0.6732519 1982-1999 –> 2000-2012 0.3073385 -2484.2238
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 636 0.1879213 -307.553145 1.4575343 2000-2012 –> 2013-2019 0.3218125 -1660.5859
Atlantic (28.15,28.2] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 636 0.1863196 -316.441232 1.4592968 2000-2012 –> 2013-2019 0.3194564 -1680.5278
Atlantic (28.15,28.2] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 636 0.1862961 -316.601851 1.4570795 2000-2012 –> 2013-2019 0.3193834 -1681.5463
Atlantic (28.15,28.2] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 636 0.1861145 -315.842051 1.4217230 2000-2012 –> 2013-2019 0.3175110 -1708.3224
Atlantic (28.15,28.2] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 636 0.1862056 -317.219763 1.4466140 2000-2012 –> 2013-2019 0.3196127 -1676.6195
Atlantic (28.2, Inf] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 955 0.1348665 -1100.455035 0.8777874 1982-1999 –> 2000-2012 0.2978964 -2062.9752
Atlantic (28.2, Inf] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 955 0.1353162 -1094.096142 0.9017636 1982-1999 –> 2000-2012 0.3019585 -2002.3092
Atlantic (28.2, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 955 0.1405463 -1021.664584 0.6550805 1982-1999 –> 2000-2012 0.2890069 -2216.1595
Atlantic (28.2, Inf] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 955 0.1334735 -1118.285100 0.7137703 1982-1999 –> 2000-2012 0.2747034 -2434.5100
Atlantic (28.2, Inf] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 955 0.1340784 -1111.648172 0.8191361 1982-1999 –> 2000-2012 0.2935903 -2128.2266
Atlantic (28.2, Inf] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 552 0.1517723 -498.944890 0.9275123 2000-2012 –> 2013-2019 0.2866388 -1599.3999
Atlantic (28.2, Inf] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 552 0.1550890 -475.078380 0.9388129 2000-2012 –> 2013-2019 0.2904053 -1569.1745
Atlantic (28.2, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 552 0.1569488 -461.918306 1.0878117 2000-2012 –> 2013-2019 0.2974951 -1483.5829
Atlantic (28.2, Inf] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 552 0.1443427 -552.355889 0.8257393 2000-2012 –> 2013-2019 0.2778162 -1670.6410
Atlantic (28.2, Inf] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 552 0.1503309 -509.479552 0.9141270 2000-2012 –> 2013-2019 0.2844093 -1621.1277
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + oxygen + aou + nitrate + phosphate_star 5288 5.9081446 33807.179508 33.2613315 1982-1999 –> 2000-2012 12.1263367 80511.1594
Indo-Pacific (-Inf,26] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 5288 5.8704994 33741.576188 32.8231038 1982-1999 –> 2000-2012 12.0412064 80337.3068
Indo-Pacific (-Inf,26] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 5288 6.0193373 34006.372533 33.5251746 1982-1999 –> 2000-2012 12.3199691 80901.7736
Indo-Pacific (-Inf,26] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 5288 5.9081398 33809.170907 33.2662064 1982-1999 –> 2000-2012 12.1260301 80514.4528
Indo-Pacific (-Inf,26] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 5288 5.8703920 33743.382691 32.8451163 1982-1999 –> 2000-2012 12.0400346 80338.6323
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + oxygen + aou + nitrate + phosphate_star 2815 5.6389160 17740.789029 27.6780849 2000-2012 –> 2013-2019 11.5470606 51547.9685
Indo-Pacific (-Inf,26] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 2815 5.5785102 17682.153450 25.7404923 2000-2012 –> 2013-2019 11.4490096 51423.7296
Indo-Pacific (-Inf,26] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 2815 5.7096684 17812.990108 26.8866125 2000-2012 –> 2013-2019 11.7290058 51819.3626
Indo-Pacific (-Inf,26] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 2815 5.6384269 17742.300738 27.7405532 2000-2012 –> 2013-2019 11.5465668 51551.4716
Indo-Pacific (-Inf,26] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 2815 5.5774824 17683.116042 25.8372410 2000-2012 –> 2013-2019 11.4478744 51426.4987
Indo-Pacific (26,26.5] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 5428 4.3894781 31478.304168 50.5213809 1982-1999 –> 2000-2012 8.4316284 74270.6513
Indo-Pacific (26,26.5] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 5428 4.4763284 31691.003687 45.9884146 1982-1999 –> 2000-2012 8.5701782 74676.4264
Indo-Pacific (26,26.5] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 5428 4.3686353 31428.633324 49.4844000 1982-1999 –> 2000-2012 8.3888228 74140.2108
Indo-Pacific (26,26.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 5428 4.5340876 31830.185580 43.0233728 1982-1999 –> 2000-2012 8.7741402 75348.6932
Indo-Pacific (26,26.5] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 5428 4.4903869 31725.045180 45.6927227 1982-1999 –> 2000-2012 8.5991679 74765.7752
Indo-Pacific (26,26.5] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 2876 4.2202535 16460.011635 56.9391393 2000-2012 –> 2013-2019 8.6097316 47938.3158
Indo-Pacific (26,26.5] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 2876 4.3452833 16627.945616 48.1945701 2000-2012 –> 2013-2019 8.8216116 48318.9493
Indo-Pacific (26,26.5] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 2876 4.2122415 16451.081280 55.7738871 2000-2012 –> 2013-2019 8.5808769 47879.7146
Indo-Pacific (26,26.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 2876 4.4815444 16805.548776 46.3106843 2000-2012 –> 2013-2019 9.0156321 48635.7344
Indo-Pacific (26,26.5] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 2876 4.3577945 16644.483316 47.6441715 2000-2012 –> 2013-2019 8.8481814 48369.5285
Indo-Pacific (26.5,26.75] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 4269 3.7987203 23526.260526 26.3099686 1982-1999 –> 2000-2012 7.4930968 55858.7566
Indo-Pacific (26.5,26.75] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 4269 3.9099557 23772.682655 26.8010367 1982-1999 –> 2000-2012 7.6834415 56356.3761
Indo-Pacific (26.5,26.75] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 4269 3.7864982 23500.745889 26.3956705 1982-1999 –> 2000-2012 7.4772907 55823.7346
Indo-Pacific (26.5,26.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 4269 3.8963414 23742.901689 26.7023232 1982-1999 –> 2000-2012 7.6360524 56220.0000
Indo-Pacific (26.5,26.75] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 4269 3.8401319 23618.833639 25.9985422 1982-1999 –> 2000-2012 7.6582559 56341.9540
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ oxygen + aou + nitrate + silicate + phosphate_star 2324 4.2164585 13297.677654 31.7627484 2000-2012 –> 2013-2019 8.1023350 37015.6169
Indo-Pacific (26.5,26.75] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 2324 4.1409652 13215.703841 31.0550390 2000-2012 –> 2013-2019 7.9396855 36741.9644
Indo-Pacific (26.5,26.75] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 2324 4.1247434 13199.460017 29.7101718 2000-2012 –> 2013-2019 7.9112416 36700.2059
Indo-Pacific (26.5,26.75] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 2324 4.2544954 13341.419522 23.3613336 2000-2012 –> 2013-2019 8.1508368 37084.3212
Indo-Pacific (26.5,26.75] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 2324 4.1797917 13259.081345 29.6661189 2000-2012 –> 2013-2019 8.0199236 36877.9150
Indo-Pacific (26.75,27] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 4557 4.1353200 25886.118633 22.9077644 1982-1999 –> 2000-2012 7.6662364 60367.9356
Indo-Pacific (26.75,27] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 4557 4.0876903 25782.536319 20.9673455 1982-1999 –> 2000-2012 7.5839847 60139.6536
Indo-Pacific (26.75,27] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 4557 4.1556325 25928.776454 23.0206730 1982-1999 –> 2000-2012 7.6982624 60451.1786
Indo-Pacific (26.75,27] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 4557 4.1322056 25879.252031 21.7852477 1982-1999 –> 2000-2012 7.6563812 60336.4988
Indo-Pacific (26.75,27] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 4557 4.0936601 25793.836953 21.1472190 1982-1999 –> 2000-2012 7.6348763 60313.1069
Indo-Pacific (26.75,27] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 2460 4.5716148 14474.920680 18.5454907 2000-2012 –> 2013-2019 8.7069348 40361.0393
Indo-Pacific (26.75,27] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 2460 4.5357389 14438.158604 17.7419483 2000-2012 –> 2013-2019 8.6234292 40220.6949
Indo-Pacific (26.75,27] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 2460 4.5918239 14496.621995 16.9140190 2000-2012 –> 2013-2019 8.7240295 40375.8740
Indo-Pacific (26.75,27] 5 cstar_tref ~ temp + oxygen + aou + nitrate + phosphate_star 2460 4.5812181 14483.245002 17.9419312 2000-2012 –> 2013-2019 8.7219819 40379.3536
Indo-Pacific (26.75,27] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 2460 4.5453248 14446.545621 17.4037469 2000-2012 –> 2013-2019 8.6389849 40240.3826
Indo-Pacific (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 4267 3.7264465 23351.315142 45.1418759 1982-1999 –> 2000-2012 6.6897195 52262.6718
Indo-Pacific (27,27.25] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 4267 3.3354520 22407.347822 35.3437306 1982-1999 –> 2000-2012 6.0441809 50284.7794
Indo-Pacific (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 4267 3.4036330 22578.034898 30.6158677 1982-1999 –> 2000-2012 6.1906243 50781.9909
Indo-Pacific (27,27.25] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 4267 3.3354527 22405.349588 35.3558695 1982-1999 –> 2000-2012 6.0610193 50352.2554
Indo-Pacific (27,27.25] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 4267 3.4041802 22577.406832 30.9255112 1982-1999 –> 2000-2012 6.2015505 50822.2366
Indo-Pacific (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 2344 4.5216543 13741.603590 32.3418461 2000-2012 –> 2013-2019 8.2481008 37092.9187
Indo-Pacific (27,27.25] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 2344 3.9285194 13084.398956 25.0306980 2000-2012 –> 2013-2019 7.2639714 35491.7468
Indo-Pacific (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 2344 3.9472749 13104.727079 23.9342640 2000-2012 –> 2013-2019 7.3509079 35682.7620
Indo-Pacific (27,27.25] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 2344 3.9337598 13088.648245 25.8018276 2000-2012 –> 2013-2019 7.2692124 35493.9978
Indo-Pacific (27,27.25] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 2344 3.9533412 13109.926205 24.4569335 2000-2012 –> 2013-2019 7.3575214 35687.3330
Indo-Pacific (27.25,27.5] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 3520 2.9102241 17525.667191 29.0753107 1982-1999 –> 2000-2012 5.7375263 40732.9062
Indo-Pacific (27.25,27.5] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 3520 2.7975878 17249.780387 29.7923176 1982-1999 –> 2000-2012 5.5567744 40228.9521
Indo-Pacific (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 3520 2.8617081 17407.314891 30.1089756 1982-1999 –> 2000-2012 5.6400680 40449.8156
Indo-Pacific (27.25,27.5] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 3520 2.8157828 17293.418946 29.8473473 1982-1999 –> 2000-2012 5.6130436 40399.8810
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 3520 2.8774108 17443.838876 30.1555714 1982-1999 –> 2000-2012 5.6901545 40600.3746
Indo-Pacific (27.25,27.5] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 1879 3.3805403 9925.748570 28.8306410 2000-2012 –> 2013-2019 6.2907644 27451.4158
Indo-Pacific (27.25,27.5] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 1879 3.2338324 9761.015166 28.3297859 2000-2012 –> 2013-2019 6.0314202 27010.7956
Indo-Pacific (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1879 3.3074238 9843.576259 29.1273995 2000-2012 –> 2013-2019 6.1691319 27250.8912
Indo-Pacific (27.25,27.5] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 1879 3.2836899 9816.511804 31.7122089 2000-2012 –> 2013-2019 6.0994726 27109.9307
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 1879 3.3526141 9892.575280 28.9794819 2000-2012 –> 2013-2019 6.2300248 27336.4142
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 3170 1.8542474 12924.886641 13.2838580 1982-1999 –> 2000-2012 3.4799789 30094.3833
Indo-Pacific (27.5,27.75] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 3170 1.8424389 12886.382103 10.4441578 1982-1999 –> 2000-2012 3.4663426 30047.7479
Indo-Pacific (27.5,27.75] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 3170 1.8427988 12887.620256 10.5806054 1982-1999 –> 2000-2012 3.4662898 30046.6975
Indo-Pacific (27.5,27.75] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 3170 1.8420371 12886.999210 10.2134889 1982-1999 –> 2000-2012 3.4619347 30028.1202
Indo-Pacific (27.5,27.75] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 3170 1.8431057 12888.676087 10.6548695 1982-1999 –> 2000-2012 3.4669502 30049.7137
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1683 2.1594104 7381.412454 13.7085041 2000-2012 –> 2013-2019 4.0136578 20306.2991
Indo-Pacific (27.5,27.75] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 1683 2.0871546 7268.855677 11.9976954 2000-2012 –> 2013-2019 3.9295935 20155.2378
Indo-Pacific (27.5,27.75] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 1683 2.0884149 7270.887475 11.9868971 2000-2012 –> 2013-2019 3.9312136 20158.5077
Indo-Pacific (27.5,27.75] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 1683 2.0863413 7269.543688 12.0287717 2000-2012 –> 2013-2019 3.9283783 20156.5429
Indo-Pacific (27.5,27.75] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 1683 2.0891042 7271.998315 11.9822574 2000-2012 –> 2013-2019 3.9322099 20160.6744
Indo-Pacific (27.75,27.85] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 1216 1.1844481 3878.540094 5.9028516 1982-1999 –> 2000-2012 2.3573798 9223.3700
Indo-Pacific (27.75,27.85] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 1216 1.0964322 3692.752023 6.0472671 1982-1999 –> 2000-2012 2.1895942 8801.8038
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 1216 1.2257896 3959.977795 5.7133287 1982-1999 –> 2000-2012 2.4231295 9372.3398
Indo-Pacific (27.75,27.85] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1216 1.1484687 3803.519113 5.5653347 1982-1999 –> 2000-2012 2.2742935 9009.9652
Indo-Pacific (27.75,27.85] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 1216 1.2976712 4100.567778 6.7996925 1982-1999 –> 2000-2012 2.5425351 9646.3369
Indo-Pacific (27.75,27.85] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 691 1.2773911 2315.313992 8.9754138 2000-2012 –> 2013-2019 2.4618392 6193.8541
Indo-Pacific (27.75,27.85] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 691 1.2492490 2286.526882 7.7397274 2000-2012 –> 2013-2019 2.3456812 5979.2789
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 691 1.3286664 2367.703966 8.2042737 2000-2012 –> 2013-2019 2.5544561 6327.6818
Indo-Pacific (27.75,27.85] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 691 1.3067937 2346.763836 7.0585944 2000-2012 –> 2013-2019 2.4552624 6150.2829
Indo-Pacific (27.75,27.85] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 691 1.5287827 2563.593078 10.4519985 2000-2012 –> 2013-2019 2.8264539 6664.1609
Indo-Pacific (27.85,27.95] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 1797 1.2569456 5937.557787 10.2423669 1982-1999 –> 2000-2012 2.3129625 13412.1248
Indo-Pacific (27.85,27.95] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 1797 1.2314246 5863.834129 6.7020452 1982-1999 –> 2000-2012 2.3640740 13693.0211
Indo-Pacific (27.85,27.95] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 1797 1.1203810 5526.190843 7.5717876 1982-1999 –> 2000-2012 2.0757088 12495.5204
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 1797 1.2870458 6020.609331 8.7933519 1982-1999 –> 2000-2012 2.3725611 13632.6372
Indo-Pacific (27.85,27.95] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1797 1.1268453 5544.867661 6.8163957 1982-1999 –> 2000-2012 2.0904038 12555.6225
Indo-Pacific (27.85,27.95] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 987 1.4148279 3501.978286 10.5168893 2000-2012 –> 2013-2019 2.6717735 9439.5361
Indo-Pacific (27.85,27.95] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 987 1.4141378 3501.015168 6.9780896 2000-2012 –> 2013-2019 2.6455624 9364.8493
Indo-Pacific (27.85,27.95] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 987 1.3054295 3345.119105 8.3031159 2000-2012 –> 2013-2019 2.4258106 8871.3099
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ sal + temp + oxygen + aou + silicate 987 1.4185912 3505.221899 7.0604103 2000-2012 –> 2013-2019 2.6559788 9384.4175
Indo-Pacific (27.85,27.95] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 987 1.3133146 3355.006629 7.3269651 2000-2012 –> 2013-2019 2.4401599 8899.8743
Indo-Pacific (27.95,28.05] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 1162 0.5791361 2044.203012 3.6621709 1982-1999 –> 2000-2012 1.1194923 4854.4693
Indo-Pacific (27.95,28.05] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 1162 0.5701419 2009.827139 3.1159933 1982-1999 –> 2000-2012 1.0958787 4726.6997
Indo-Pacific (27.95,28.05] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1162 0.5749179 2027.213849 2.9262058 1982-1999 –> 2000-2012 1.1061869 4778.4935
Indo-Pacific (27.95,28.05] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 1162 0.5712227 2012.228499 3.2380940 1982-1999 –> 2000-2012 1.1008975 4753.0557
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 1162 0.5749226 2025.232979 2.9259808 1982-1999 –> 2000-2012 1.1074497 4782.7387
Indo-Pacific (27.95,28.05] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 648 0.6048430 1203.333200 2.5954975 2000-2012 –> 2013-2019 1.1839791 3247.5362
Indo-Pacific (27.95,28.05] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 648 0.5866412 1165.733340 2.7215486 2000-2012 –> 2013-2019 1.1567831 3175.5605
Indo-Pacific (27.95,28.05] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 648 0.5887663 1168.419405 2.7201223 2000-2012 –> 2013-2019 1.1636841 3195.6333
Indo-Pacific (27.95,28.05] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 648 0.5916062 1174.655608 2.6354676 2000-2012 –> 2013-2019 1.1628288 3186.8841
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 648 0.5922331 1174.028301 2.7169734 2000-2012 –> 2013-2019 1.1671557 3199.2613
Indo-Pacific (28.05,28.1] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 1015 0.4475786 1264.521751 3.1558350 1982-1999 –> 2000-2012 0.8951399 3148.8818
Indo-Pacific (28.05,28.1] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 1015 0.4482077 1267.373242 3.1714147 1982-1999 –> 2000-2012 0.8962295 3154.8573
Indo-Pacific (28.05,28.1] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 1015 0.4343993 1205.849169 2.4951723 1982-1999 –> 2000-2012 0.8698140 3008.6203
Indo-Pacific (28.05,28.1] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 1015 0.4481986 1267.332079 3.1780378 1982-1999 –> 2000-2012 0.8962258 3154.8527
Indo-Pacific (28.05,28.1] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 1015 0.4378458 1219.891662 2.7719565 1982-1999 –> 2000-2012 0.8752364 3034.4176
Indo-Pacific (28.05,28.1] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 565 0.4493063 715.343489 2.4314701 2000-2012 –> 2013-2019 0.8968848 1979.8652
Indo-Pacific (28.05,28.1] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 565 0.4499026 716.842192 2.4409361 2000-2012 –> 2013-2019 0.8981103 1984.2154
Indo-Pacific (28.05,28.1] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 565 0.4351399 681.141559 2.1322710 2000-2012 –> 2013-2019 0.8695392 1886.9907
Indo-Pacific (28.05,28.1] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 565 0.4498014 716.587923 2.4435171 2000-2012 –> 2013-2019 0.8980000 1983.9200
Indo-Pacific (28.05,28.1] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 565 0.4395227 690.466051 2.2371776 2000-2012 –> 2013-2019 0.8773685 1910.3577
Indo-Pacific (28.1, Inf] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 9318 0.6578925 18656.188640 2.8836217 1982-1999 –> 2000-2012 1.3009615 44038.5898
Indo-Pacific (28.1, Inf] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 9318 0.6595803 18703.939183 2.9337738 1982-1999 –> 2000-2012 1.3041520 44146.9097
Indo-Pacific (28.1, Inf] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 9318 0.6427485 18224.193202 2.2172137 1982-1999 –> 2000-2012 1.2769741 43249.2293
Indo-Pacific (28.1, Inf] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 9318 0.6591098 18690.639245 2.9222470 1982-1999 –> 2000-2012 1.3033084 44118.5877
Indo-Pacific (28.1, Inf] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 9318 0.6502613 18438.759039 2.6265795 1982-1999 –> 2000-2012 1.2867521 43554.3169
Indo-Pacific (28.1, Inf] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 5093 0.6618557 10265.467303 3.2550523 2000-2012 –> 2013-2019 1.3197482 28921.6559
Indo-Pacific (28.1, Inf] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 5093 0.6637355 10294.355861 3.3154593 2000-2012 –> 2013-2019 1.3233158 28998.2950
Indo-Pacific (28.1, Inf] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 5093 0.6410498 9942.121503 2.2137429 2000-2012 –> 2013-2019 1.2837983 28166.3147
Indo-Pacific (28.1, Inf] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 5093 0.6631944 10286.049205 3.2998121 2000-2012 –> 2013-2019 1.3223042 28976.6884
Indo-Pacific (28.1, Inf] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 5093 0.6520006 10112.655841 2.8653159 2000-2012 –> 2013-2019 1.3022619 28551.4149
SO (26,26.5] 5 cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 103 0.7055448 234.451622 1.8539995 1982-1999 –> 2000-2012 1.6941305 611.5759
SO (26,26.5] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 103 0.7003038 234.915698 2.1163299 1982-1999 –> 2000-2012 1.6845459 612.9039
SO (26,26.5] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 103 0.7045483 236.160473 1.9249255 1982-1999 –> 2000-2012 1.6930212 615.2554
SO (26,26.5] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 103 0.7002576 236.902081 2.1117984 1982-1999 –> 2000-2012 1.6844189 616.8691
SO (26,26.5] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 103 0.7142789 238.986090 2.1855914 1982-1999 –> 2000-2012 1.7045995 618.5628
SO (26,26.5] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 47 0.7766527 125.620592 1.6121850 2000-2012 –> 2013-2019 1.4769565 360.5363
SO (26,26.5] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 47 0.8127027 129.885565 1.5126948 2000-2012 –> 2013-2019 1.5172510 366.0460
SO (26,26.5] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 47 0.7556873 125.048221 1.5320040 2000-2012 –> 2013-2019 1.4559448 361.9503
SO (26,26.5] 5 cstar_tref ~ sal + temp + oxygen + nitrate + phosphate_star 47 0.8132086 127.944071 1.5202774 2000-2012 –> 2013-2019 1.5275474 364.9474
SO (26,26.5] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 47 0.7585524 123.403948 1.5618342 2000-2012 –> 2013-2019 1.4728313 362.3900
SO (26.5,26.75] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 312 1.4083536 1115.088585 7.4275774 1982-1999 –> 2000-2012 3.0615996 2745.2991
SO (26.5,26.75] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 312 1.2717338 1051.415482 5.9028790 1982-1999 –> 2000-2012 2.7349528 2579.0604
SO (26.5,26.75] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 312 1.2281813 1031.671130 6.3708166 1982-1999 –> 2000-2012 2.6553109 2540.3382
SO (26.5,26.75] 5 cstar_tref ~ temp + oxygen + aou + nitrate + phosphate_star 312 1.3885330 1104.244260 6.3863351 1982-1999 –> 2000-2012 2.9751399 2697.8947
SO (26.5,26.75] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 312 1.3843720 1104.371537 6.2836537 1982-1999 –> 2000-2012 2.9683453 2698.6265
SO (26.5,26.75] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 133 1.2241241 447.229644 6.1629627 2000-2012 –> 2013-2019 2.6324777 1562.3182
SO (26.5,26.75] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 133 1.2387812 450.395698 4.5972243 2000-2012 –> 2013-2019 2.5105150 1501.8112
SO (26.5,26.75] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 133 1.1506016 432.753437 3.8170539 2000-2012 –> 2013-2019 2.3787828 1464.4246
SO (26.5,26.75] 5 cstar_tref ~ temp + oxygen + aou + nitrate + phosphate_star 133 1.3368487 468.661481 5.1221114 2000-2012 –> 2013-2019 2.7253817 1572.9057
SO (26.5,26.75] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 133 1.3230060 467.892748 4.8994666 2000-2012 –> 2013-2019 2.7073780 1572.2643
SO (26.75,27] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 1061 2.0677789 4568.567596 15.0650630 1982-1999 –> 2000-2012 4.3600517 11397.4691
SO (26.75,27] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 1061 1.9578212 4454.615506 12.5612675 1982-1999 –> 2000-2012 4.1070569 11090.2895
SO (26.75,27] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1061 2.2880755 4783.390414 9.3213970 1982-1999 –> 2000-2012 4.7000755 11766.5567
SO (26.75,27] 5 cstar_tref ~ temp + oxygen + aou + nitrate + phosphate_star 1061 2.0678244 4566.614247 15.0903489 1982-1999 –> 2000-2012 4.3601033 11393.5238
SO (26.75,27] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 1061 1.9930102 4490.416732 13.9007340 1982-1999 –> 2000-2012 4.1824842 11180.2943
SO (26.75,27] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 550 2.1237258 2405.321609 19.3107447 2000-2012 –> 2013-2019 4.1915047 6973.8892
SO (26.75,27] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 550 1.9412886 2308.519586 15.1544756 2000-2012 –> 2013-2019 3.8991098 6763.1351
SO (26.75,27] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 550 2.2521270 2469.895006 11.5622287 2000-2012 –> 2013-2019 4.5402025 7253.2854
SO (26.75,27] 5 cstar_tref ~ temp + oxygen + aou + nitrate + phosphate_star 550 2.1304645 2406.806408 19.0820438 2000-2012 –> 2013-2019 4.1982888 6973.4207
SO (26.75,27] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 550 1.9601216 2317.139572 16.1108570 2000-2012 –> 2013-2019 3.9531319 6807.5563
SO (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 2513 1.7142081 9856.353884 16.6315870 1982-1999 –> 2000-2012 3.3021064 22808.5153
SO (27,27.25] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 2513 1.6110858 9546.526494 12.7970347 1982-1999 –> 2000-2012 3.1569488 22316.2121
SO (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 2513 1.6118710 9546.975590 12.9579746 1982-1999 –> 2000-2012 3.1585654 22318.3580
SO (27,27.25] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 2513 1.6593137 9692.772001 16.3534096 1982-1999 –> 2000-2012 3.2638856 22716.7578
SO (27,27.25] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 2513 1.6615541 9697.553575 16.3323879 1982-1999 –> 2000-2012 3.2687331 22730.7024
SO (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 1286 1.9250058 5349.987299 11.2278081 2000-2012 –> 2013-2019 3.6392139 15206.3412
SO (27,27.25] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 1286 1.7517142 5109.359907 7.4475239 2000-2012 –> 2013-2019 3.3628000 14655.8864
SO (27,27.25] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1286 1.7553735 5112.727102 7.7602872 2000-2012 –> 2013-2019 3.3672445 14659.7027
SO (27,27.25] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 1286 1.8477238 5244.600909 11.4645683 2000-2012 –> 2013-2019 3.5070375 14937.3729
SO (27,27.25] 5 cstar_tref ~ temp + oxygen + aou + silicate + phosphate_star 1286 1.8504129 5246.341361 11.4263441 2000-2012 –> 2013-2019 3.5119670 14943.8949
SO (27.25,27.5] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 2731 1.6619933 10541.034809 7.7450245 1982-1999 –> 2000-2012 3.0787039 24002.4725
SO (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 2731 1.6617833 10540.344591 7.3898985 1982-1999 –> 2000-2012 3.0863767 24043.9971
SO (27.25,27.5] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 2731 1.6597560 10535.677129 7.5354082 1982-1999 –> 2000-2012 3.0703213 23966.0415
SO (27.25,27.5] 5 cstar_tref ~ sal + temp + oxygen + nitrate + silicate 2731 1.6623528 10540.215983 7.4913929 1982-1999 –> 2000-2012 3.0903673 24060.1171
SO (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 2731 1.6603021 10535.474074 7.6355850 1982-1999 –> 2000-2012 3.0747466 23984.7321
SO (27.25,27.5] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 1453 1.8185695 5877.369215 8.5799014 2000-2012 –> 2013-2019 3.4805628 16418.4040
SO (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 1453 1.8209112 5881.108844 8.3718926 2000-2012 –> 2013-2019 3.4826945 16421.4534
SO (27.25,27.5] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 1453 1.8181682 5878.727950 8.6011021 2000-2012 –> 2013-2019 3.4779242 16414.4051
SO (27.25,27.5] 5 cstar_tref ~ sal + temp + oxygen + nitrate + silicate 1453 1.8228962 5882.274894 8.3460548 2000-2012 –> 2013-2019 3.4852489 16422.4909
SO (27.25,27.5] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 1453 1.8200698 5879.765771 8.5793746 2000-2012 –> 2013-2019 3.4803720 16415.2398
SO (27.5,27.75] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 3391 1.4729465 12265.671022 7.9500627 1982-1999 –> 2000-2012 2.7307449 27264.7308
SO (27.5,27.75] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 3391 1.4672435 12241.361201 7.7881314 1982-1999 –> 2000-2012 2.7196997 27203.7313
SO (27.5,27.75] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 3391 1.4674316 12240.230996 7.8118195 1982-1999 –> 2000-2012 2.7199581 27201.1108
SO (27.5,27.75] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 3391 1.4728162 12265.070928 7.7521860 1982-1999 –> 2000-2012 2.7254069 27226.4176
SO (27.5,27.75] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 3391 1.4728687 12263.313048 7.7650825 1982-1999 –> 2000-2012 2.7255384 27223.2322
SO (27.5,27.75] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 1821 1.6691334 7049.587393 9.6582328 2000-2012 –> 2013-2019 3.1420799 19315.2584
SO (27.5,27.75] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 1821 1.6678684 7048.826151 9.6098696 2000-2012 –> 2013-2019 3.1351119 19290.1874
SO (27.5,27.75] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 1821 1.6680863 7047.301950 9.5902240 2000-2012 –> 2013-2019 3.1355180 19287.5329
SO (27.5,27.75] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 1821 1.6781054 7069.111499 9.6845747 2000-2012 –> 2013-2019 3.1509215 19334.1824
SO (27.5,27.75] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 1821 1.6789480 7068.939886 9.6484521 2000-2012 –> 2013-2019 3.1518168 19332.2529
SO (27.75,27.85] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 1415 1.2886716 4749.317831 8.6923687 1982-1999 –> 2000-2012 2.4939862 10561.8196
SO (27.75,27.85] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 1415 1.2748827 4720.873270 9.3744209 1982-1999 –> 2000-2012 2.4774237 10527.0584
SO (27.75,27.85] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 1415 1.2820314 4734.697745 8.8776002 1982-1999 –> 2000-2012 2.4846208 10539.0283
SO (27.75,27.85] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 1415 1.2769205 4723.393277 9.4127066 1982-1999 –> 2000-2012 2.4979546 10582.6716
SO (27.75,27.85] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 1415 1.2832711 4735.433018 8.9326433 1982-1999 –> 2000-2012 2.5044896 10593.2564
SO (27.75,27.85] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 722 1.3897211 2540.172099 8.3374398 2000-2012 –> 2013-2019 2.6783927 7289.4899
SO (27.75,27.85] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 722 1.3708261 2522.404455 7.9316494 2000-2012 –> 2013-2019 2.6457088 7243.2777
SO (27.75,27.85] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 722 1.3811068 2531.193485 8.2222107 2000-2012 –> 2013-2019 2.6631382 7265.8912
SO (27.75,27.85] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 722 1.3768681 2526.754899 8.0018925 2000-2012 –> 2013-2019 2.6537886 7250.1482
SO (27.75,27.85] 5 cstar_tref ~ temp + oxygen + nitrate + silicate + phosphate_star 722 1.3849071 2533.161377 8.2492850 2000-2012 –> 2013-2019 2.6681782 7268.5944
SO (27.85,27.95] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 1399 1.3260264 4775.748645 8.5642886 1982-1999 –> 2000-2012 2.6449160 11076.1192
SO (27.85,27.95] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 1399 1.2870785 4692.335133 8.5164180 1982-1999 –> 2000-2012 2.5480100 10826.1600
SO (27.85,27.95] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 1399 1.2744015 4666.639762 8.1267517 1982-1999 –> 2000-2012 2.5313459 10790.7277
SO (27.85,27.95] 5 cstar_tref ~ sal + temp + oxygen + aou + silicate 1399 1.2912832 4699.460807 8.0313297 1982-1999 –> 2000-2012 2.5525237 10832.1940
SO (27.85,27.95] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1399 1.2836451 4684.861093 8.5215752 1982-1999 –> 2000-2012 2.5448741 10819.5605
SO (27.85,27.95] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 709 1.4718952 2576.183923 7.4972945 2000-2012 –> 2013-2019 2.7979216 7351.9326
SO (27.85,27.95] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 709 1.4977451 2600.871146 8.8681972 2000-2012 –> 2013-2019 2.7848237 7293.2063
SO (27.85,27.95] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 709 1.4376205 2544.773643 7.5356670 2000-2012 –> 2013-2019 2.7120220 7211.4134
SO (27.85,27.95] 5 cstar_tref ~ sal + temp + oxygen + aou + silicate 709 1.4996937 2600.714760 8.5827164 2000-2012 –> 2013-2019 2.7909769 7300.1756
SO (27.85,27.95] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 709 1.4880087 2591.623013 9.0787780 2000-2012 –> 2013-2019 2.7716538 7276.4841
SO (27.95,28.05] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 1786 1.4829516 6491.939519 5.7404493 1982-1999 –> 2000-2012 2.7849116 14878.2330
SO (27.95,28.05] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 1786 1.4719049 6465.231467 5.5467278 1982-1999 –> 2000-2012 2.7642890 14814.8065
SO (27.95,28.05] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 1786 1.5258734 6593.857755 8.3812420 1982-1999 –> 2000-2012 2.8092581 14908.6751
SO (27.95,28.05] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 1786 1.4364245 6380.073446 5.8268292 1982-1999 –> 2000-2012 2.6759000 14523.7332
SO (27.95,28.05] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 1786 1.4936018 6517.500963 5.9646450 1982-1999 –> 2000-2012 2.8074255 14948.9135
SO (27.95,28.05] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 986 1.6047359 3746.822359 5.6205654 2000-2012 –> 2013-2019 3.0876876 10238.7619
SO (27.95,28.05] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 986 1.5938488 3733.397940 5.6119621 2000-2012 –> 2013-2019 3.0657536 10198.6294
SO (27.95,28.05] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 986 1.6261350 3772.945116 8.0450596 2000-2012 –> 2013-2019 3.1520085 10366.8029
SO (27.95,28.05] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 986 1.5531954 3684.446680 6.0836268 2000-2012 –> 2013-2019 2.9896199 10064.5201
SO (27.95,28.05] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 986 1.6150504 3759.456807 5.6719144 2000-2012 –> 2013-2019 3.1086522 10276.9578
SO (28.05,28.1] 5 cstar_tref ~ sal + temp + oxygen + aou + nitrate 1071 1.0450829 3147.820495 8.4427617 1982-1999 –> 2000-2012 1.9983465 7313.4345
SO (28.05,28.1] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 1071 1.0447507 3149.139516 8.5646872 1982-1999 –> 2000-2012 1.9976873 7315.7149
SO (28.05,28.1] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 1071 1.0389462 3137.205569 8.6555860 1982-1999 –> 2000-2012 1.9898140 7297.2000
SO (28.05,28.1] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 1071 1.0293063 3119.238088 7.8770763 1982-1999 –> 2000-2012 1.9778296 7273.7576
SO (28.05,28.1] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 1071 1.0473850 3154.533537 9.0409967 1982-1999 –> 2000-2012 2.0055492 7337.6746
SO (28.05,28.1] 5 cstar_tref ~ sal + temp + oxygen + aou + nitrate 606 1.2264599 1981.161326 9.4690189 2000-2012 –> 2013-2019 2.2715428 5128.9818
SO (28.05,28.1] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 606 1.2258531 1982.561539 9.7027730 2000-2012 –> 2013-2019 2.2706038 5131.7011
SO (28.05,28.1] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 606 1.2236721 1980.403317 9.7409069 2000-2012 –> 2013-2019 2.2626183 5117.6089
SO (28.05,28.1] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 606 1.2220173 1980.763201 9.3383523 2000-2012 –> 2013-2019 2.2513236 5100.0013
SO (28.05,28.1] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 606 1.2282847 1984.963330 10.1968303 2000-2012 –> 2013-2019 2.2756697 5139.4969
SO (28.1, Inf] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 5309 1.0911712 16008.726846 6.2189622 1982-1999 –> 2000-2012 2.0576363 36406.6328
SO (28.1, Inf] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 5309 1.0913161 16010.136574 6.2151315 1982-1999 –> 2000-2012 2.0580529 36412.1785
SO (28.1, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 5309 1.0896709 15994.117646 6.2333995 1982-1999 –> 2000-2012 2.0589958 36435.5098
SO (28.1, Inf] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 5309 1.0896439 15995.855042 6.2403672 1982-1999 –> 2000-2012 2.0549755 36378.4873
SO (28.1, Inf] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 5309 1.0912186 16009.188354 6.2165223 1982-1999 –> 2000-2012 2.0578536 36409.6807
SO (28.1, Inf] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 2711 1.1818780 8615.526404 6.2714947 2000-2012 –> 2013-2019 2.2730492 24624.2532
SO (28.1, Inf] 6 cstar_tref ~ sal + temp + aou + nitrate + silicate + phosphate_star 2711 1.1819169 8615.705040 6.2707898 2000-2012 –> 2013-2019 2.2732330 24625.8416
SO (28.1, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 2711 1.1818836 8615.551963 6.2901289 2000-2012 –> 2013-2019 2.2715545 24609.6696
SO (28.1, Inf] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 2711 1.1814651 8615.632015 6.2643722 2000-2012 –> 2013-2019 2.2711091 24611.4871
SO (28.1, Inf] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 2711 1.1818567 8615.428461 6.2712108 2000-2012 –> 2013-2019 2.2730753 24624.6168
SO (28.1,28.15] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 395 0.4678912 536.941101 1.9309444 1982-1999 –> 2000-2012 0.8900713 1132.9342
SO (28.1,28.15] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 395 0.4342795 478.048612 1.8794440 1982-1999 –> 2000-2012 0.8276507 1000.3947
SO (28.1,28.15] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 395 0.4483007 503.151413 1.3556998 1982-1999 –> 2000-2012 0.8827308 1128.9490
SO (28.1,28.15] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 395 0.4324991 476.803113 1.7823800 1982-1999 –> 2000-2012 0.8239699 996.1034
SO (28.1,28.15] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 395 0.4412700 490.663828 1.3656130 1982-1999 –> 2000-2012 0.8604524 1079.2317
SO (28.1,28.15] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 209 0.4664365 290.335578 2.2055021 2000-2012 –> 2013-2019 0.9007161 768.3842
SO (28.1,28.15] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 209 0.4655560 289.545781 1.9653861 2000-2012 –> 2013-2019 0.9138567 792.6972
SO (28.1,28.15] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 209 0.4616073 287.985301 2.1074848 2000-2012 –> 2013-2019 0.8941064 764.7884
SO (28.1,28.15] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 209 0.4863454 305.806779 1.9335059 2000-2012 –> 2013-2019 0.9551804 842.3397
SO (28.1,28.15] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 209 0.4627191 286.990868 2.0315187 2000-2012 –> 2013-2019 0.9039892 777.6547
SO (28.15,28.2] 6 cstar_tref ~ sal + oxygen + aou + nitrate + silicate + phosphate_star 385 0.2987639 178.344514 1.1612381 1982-1999 –> 2000-2012 0.5576342 260.6397
SO (28.15,28.2] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 385 0.2801783 128.889346 0.9612678 1982-1999 –> 2000-2012 0.5513348 256.7189
SO (28.15,28.2] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 385 0.2767731 121.473654 1.0089753 1982-1999 –> 2000-2012 0.5244894 162.5184
SO (28.15,28.2] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 385 0.2997862 180.974638 1.1660905 1982-1999 –> 2000-2012 0.5593828 266.0212
SO (28.15,28.2] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 385 0.2855045 143.389637 1.2148661 1982-1999 –> 2000-2012 0.5435167 222.4242
SO (28.15,28.2] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 183 0.3211214 119.578890 0.8292375 2000-2012 –> 2013-2019 0.6012997 248.4682
SO (28.15,28.2] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 183 0.3202129 120.542025 0.8052040 2000-2012 –> 2013-2019 0.5969861 242.0157
SO (28.15,28.2] 5 cstar_tref ~ sal + temp + oxygen + aou + phosphate_star 183 0.3300672 127.635476 0.8682035 2000-2012 –> 2013-2019 0.6258955 296.3768
SO (28.15,28.2] 6 cstar_tref ~ sal + temp + oxygen + aou + silicate + phosphate_star 183 0.3212016 119.670270 0.7924252 2000-2012 –> 2013-2019 0.6120461 277.3287
SO (28.15,28.2] 6 cstar_tref ~ temp + oxygen + aou + nitrate + silicate + phosphate_star 183 0.3277625 127.070986 0.8280325 2000-2012 –> 2013-2019 0.6132671 270.4606
SO (28.2, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 1738 0.3490755 1289.855319 1.8275044 1982-1999 –> 2000-2012 0.6397759 2212.9619
SO (28.2, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 1738 0.3508194 1307.177393 1.8206271 1982-1999 –> 2000-2012 0.6420929 2240.0211
SO (28.2, Inf] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 1738 0.3484018 1285.139614 1.8445807 1982-1999 –> 2000-2012 0.6387632 2204.4781
SO (28.2, Inf] 5 cstar_tref ~ sal + temp + oxygen + nitrate + phosphate_star 1738 0.3649501 1442.440521 1.8935493 1982-1999 –> 2000-2012 0.6672204 2556.5054
SO (28.2, Inf] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 1738 0.3649328 1444.276128 1.8966743 1982-1999 –> 2000-2012 0.6671980 2560.2572
SO (28.2, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + phosphate_star 1005 0.3782848 914.129374 2.5487447 2000-2012 –> 2013-2019 0.7273603 2203.9847
SO (28.2, Inf] 6 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate 1005 0.3808359 927.639329 2.5626963 2000-2012 –> 2013-2019 0.7316554 2234.8167
SO (28.2, Inf] 7 cstar_tref ~ sal + temp + oxygen + aou + nitrate + silicate + phosphate_star 1005 0.3772774 910.769725 2.5394974 2000-2012 –> 2013-2019 0.7256792 2195.9093
SO (28.2, Inf] 5 cstar_tref ~ sal + temp + oxygen + nitrate + phosphate_star 1005 0.3929413 988.535520 2.5681945 2000-2012 –> 2013-2019 0.7578914 2430.9760
SO (28.2, Inf] 6 cstar_tref ~ sal + temp + oxygen + nitrate + silicate + phosphate_star 1005 0.3928015 989.820073 2.5650272 2000-2012 –> 2013-2019 0.7577343 2434.0962

5.3.2 per fitting unit

lm_best %>%
  select(-n_predictors) %>%
  group_by(basin, gamma_slab, eras) %>% 
  summarise(rmse_sum_mean = mean(rmse_sum),
            ais_sum_mean = mean(aic_sum)) %>% 
  ungroup() %>% 
  kable() %>%
  add_header_above() %>%
  kable_styling() %>%
  scroll_box(width = "100%", height = "400px")
basin gamma_slab eras rmse_sum_mean ais_sum_mean
Atlantic (-Inf,26] 1982-1999 –> 2000-2012 2.3534195 1197.4154
Atlantic (-Inf,26] 2000-2012 –> 2013-2019 2.1061726 703.1034
Atlantic (26,26.5] 1982-1999 –> 2000-2012 6.4433220 13245.1826
Atlantic (26,26.5] 2000-2012 –> 2013-2019 6.2842745 8271.6521
Atlantic (26.5,26.75] 1982-1999 –> 2000-2012 6.1677547 16266.7277
Atlantic (26.5,26.75] 2000-2012 –> 2013-2019 6.4848547 10422.3630
Atlantic (26.75,27] 1982-1999 –> 2000-2012 3.5435292 19093.4335
Atlantic (26.75,27] 2000-2012 –> 2013-2019 4.3536332 14182.3792
Atlantic (27,27.25] 1982-1999 –> 2000-2012 3.4065270 15317.0101
Atlantic (27,27.25] 2000-2012 –> 2013-2019 4.1640546 11494.5603
Atlantic (27.25,27.5] 1982-1999 –> 2000-2012 3.6319078 15460.5840
Atlantic (27.25,27.5] 2000-2012 –> 2013-2019 4.4645203 11082.0126
Atlantic (27.5,27.75] 1982-1999 –> 2000-2012 4.1171455 20783.1675
Atlantic (27.5,27.75] 2000-2012 –> 2013-2019 4.7029686 14395.5818
Atlantic (27.75,27.85] 1982-1999 –> 2000-2012 3.2610358 6450.0410
Atlantic (27.75,27.85] 2000-2012 –> 2013-2019 3.2816395 4171.7965
Atlantic (27.85,27.95] 1982-1999 –> 2000-2012 5.6164194 8172.2251
Atlantic (27.85,27.95] 2000-2012 –> 2013-2019 5.6809694 5459.0636
Atlantic (27.95,28.05] 1982-1999 –> 2000-2012 4.6983419 7804.7286
Atlantic (27.95,28.05] 2000-2012 –> 2013-2019 4.3634574 4976.1630
Atlantic (28.05,28.1] 1982-1999 –> 2000-2012 0.3029216 -1094.4620
Atlantic (28.05,28.1] 2000-2012 –> 2013-2019 0.4152146 -210.1766
Atlantic (28.1,28.15] 1982-1999 –> 2000-2012 0.2521351 -1717.7120
Atlantic (28.1,28.15] 2000-2012 –> 2013-2019 0.2596131 -989.9143
Atlantic (28.15,28.2] 1982-1999 –> 2000-2012 0.3070275 -2491.1279
Atlantic (28.15,28.2] 2000-2012 –> 2013-2019 0.3195552 -1681.5204
Atlantic (28.2, Inf] 1982-1999 –> 2000-2012 0.2914311 -2168.8361
Atlantic (28.2, Inf] 2000-2012 –> 2013-2019 0.2873529 -1588.7852
Indo-Pacific (-Inf,26] 1982-1999 –> 2000-2012 12.1307154 80520.6650
Indo-Pacific (-Inf,26] 2000-2012 –> 2013-2019 11.5439034 51553.8062
Indo-Pacific (26,26.5] 1982-1999 –> 2000-2012 8.5527875 74640.3514
Indo-Pacific (26,26.5] 2000-2012 –> 2013-2019 8.7752067 48228.4485
Indo-Pacific (26.5,26.75] 1982-1999 –> 2000-2012 7.5896275 56120.1642
Indo-Pacific (26.5,26.75] 2000-2012 –> 2013-2019 8.0248045 36884.0047
Indo-Pacific (26.75,27] 1982-1999 –> 2000-2012 7.6479482 60321.6747
Indo-Pacific (26.75,27] 2000-2012 –> 2013-2019 8.6830721 40315.4689
Indo-Pacific (27,27.25] 1982-1999 –> 2000-2012 6.2374189 50900.7868
Indo-Pacific (27,27.25] 2000-2012 –> 2013-2019 7.4979428 35889.7517
Indo-Pacific (27.25,27.5] 1982-1999 –> 2000-2012 5.6475134 40482.3859
Indo-Pacific (27.25,27.5] 2000-2012 –> 2013-2019 6.1641628 27231.8895
Indo-Pacific (27.5,27.75] 1982-1999 –> 2000-2012 3.4682992 30053.3325
Indo-Pacific (27.5,27.75] 2000-2012 –> 2013-2019 3.9470106 20187.4524
Indo-Pacific (27.75,27.85] 1982-1999 –> 2000-2012 2.3573864 9210.7631
Indo-Pacific (27.75,27.85] 2000-2012 –> 2013-2019 2.5287386 6263.0517
Indo-Pacific (27.85,27.95] 1982-1999 –> 2000-2012 2.2431420 13157.7852
Indo-Pacific (27.85,27.95] 2000-2012 –> 2013-2019 2.5678570 9191.9974
Indo-Pacific (27.95,28.05] 1982-1999 –> 2000-2012 1.1059810 4779.0914
Indo-Pacific (27.95,28.05] 2000-2012 –> 2013-2019 1.1668862 3200.9751
Indo-Pacific (28.05,28.1] 1982-1999 –> 2000-2012 0.8865291 3100.3259
Indo-Pacific (28.05,28.1] 2000-2012 –> 2013-2019 0.8879806 1949.0698
Indo-Pacific (28.1, Inf] 1982-1999 –> 2000-2012 1.2944296 43821.5267
Indo-Pacific (28.1, Inf] 2000-2012 –> 2013-2019 1.3102857 28722.8738
SO (26,26.5] 1982-1999 –> 2000-2012 1.6921432 615.0334
SO (26,26.5] 2000-2012 –> 2013-2019 1.4901062 363.1740
SO (26.5,26.75] 1982-1999 –> 2000-2012 2.8790697 2652.2438
SO (26.5,26.75] 2000-2012 –> 2013-2019 2.5909070 1534.7448
SO (26.75,27] 1982-1999 –> 2000-2012 4.3419543 11365.6267
SO (26.75,27] 2000-2012 –> 2013-2019 4.1564475 6954.2573
SO (27,27.25] 1982-1999 –> 2000-2012 3.2300479 22578.1091
SO (27,27.25] 2000-2012 –> 2013-2019 3.4776526 14880.6396
SO (27.25,27.5] 1982-1999 –> 2000-2012 3.0801032 24011.4720
SO (27.25,27.5] 2000-2012 –> 2013-2019 3.4813605 16418.3987
SO (27.5,27.75] 1982-1999 –> 2000-2012 2.7242696 27223.8445
SO (27.5,27.75] 2000-2012 –> 2013-2019 3.1430896 19311.8828
SO (27.75,27.85] 1982-1999 –> 2000-2012 2.4916950 10560.7669
SO (27.75,27.85] 2000-2012 –> 2013-2019 2.6618413 7263.4803
SO (27.85,27.95] 1982-1999 –> 2000-2012 2.5643340 10868.9523
SO (27.85,27.95] 2000-2012 –> 2013-2019 2.7714796 7286.6424
SO (27.95,28.05] 1982-1999 –> 2000-2012 2.7683568 14814.8723
SO (27.95,28.05] 2000-2012 –> 2013-2019 3.0807444 10229.1344
SO (28.05,28.1] 1982-1999 –> 2000-2012 1.9938453 7307.5563
SO (28.05,28.1] 2000-2012 –> 2013-2019 2.2663516 5123.5580
SO (28.1, Inf] 1982-1999 –> 2000-2012 2.0575028 36408.4978
SO (28.1, Inf] 2000-2012 –> 2013-2019 2.2724042 24619.1737
SO (28.1,28.15] 1982-1999 –> 2000-2012 0.8569750 1067.5226
SO (28.1,28.15] 2000-2012 –> 2013-2019 0.9135697 789.1728
SO (28.15,28.2] 1982-1999 –> 2000-2012 0.5472716 233.6645
SO (28.15,28.2] 2000-2012 –> 2013-2019 0.6098989 266.9300
SO (28.2, Inf] 1982-1999 –> 2000-2012 0.6510101 2354.8447
SO (28.2, Inf] 2000-2012 –> 2013-2019 0.7400641 2299.9566

5.4 Target variable coefficients

A data frame to map the target variable is prepared.

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

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

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

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

lm_best_target <- left_join(lm_best_target, lm_all_fitted_wide)

rm(eras_era, eras_forward, eras_backward,
   lm_all_fitted)

5.5 Plot selected model residuals

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

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

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

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

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

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

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

5.6 Cant coeffcients

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

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

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

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

5.7 Write files

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

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

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

6 Model diagnotics

6.1 Selection criterion vs predictors

The selection criterion (rmse) was plotted against the number of predictors (limited to 2 - 9).

6.1.1 All models

lm_all_fitted_wide %>%
  ggplot(aes(as.factor(n_predictors),
             !!sym(params_local$MLR_criterion),
             col = basin)) +
  geom_hline(yintercept = 10) +
  geom_boxplot() +
  facet_grid(gamma_slab~era) +
  scale_color_brewer(palette = "Set1") +
  labs(x="Number of predictors")

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

6.3 Predictor counts

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

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

# print table
lm_all_stats %>%
  gt(rowname_col = "gamma_slab",
     groupname_col = c("basin", "eras")) %>% 
  summary_rows(
    groups = TRUE,
    fns = list(total = "sum")
  )
aou nitrate oxygen phosphate_star sal silicate temp
Atlantic - 1982-1999 --> 2000-2012
(-Inf,26] 4 5 4 5 5 4 4
(26,26.5] 5 5 5 3 5 4 3
(26.5,26.75] 5 3 5 5 3 5 4
(26.75,27] 5 3 5 5 3 5 4
(27,27.25] 4 4 4 5 5 5 4
(27.25,27.5] 5 5 3 4 5 3 5
(27.5,27.75] 4 5 5 3 5 3 5
(27.75,27.85] 5 3 5 5 5 3 4
(27.85,27.95] 5 3 5 5 3 5 4
(27.95,28.05] 5 5 4 5 3 3 5
(28.05,28.1] 5 4 5 5 3 5 3
(28.1,28.15] 4 5 4 5 5 4 4
(28.15,28.2] 4 5 4 5 4 5 4
(28.2, Inf] 4 5 4 5 5 4 4
total 64.00 60.00 62.00 65.00 59.00 58.00 57.00
Atlantic - 2000-2012 --> 2013-2019
(-Inf,26] 4 5 4 5 5 4 4
(26,26.5] 5 5 4 3 5 5 3
(26.5,26.75] 5 3 5 5 3 5 4
(26.75,27] 5 3 5 5 3 5 4
(27,27.25] 4 4 4 5 5 5 4
(27.25,27.5] 3 5 5 4 5 3 5
(27.5,27.75] 5 4 5 3 5 3 5
(27.75,27.85] 5 3 5 4 5 3 5
(27.85,27.95] 5 3 5 5 3 5 4
(27.95,28.05] 5 5 4 5 3 3 5
(28.05,28.1] 5 3 5 5 5 3 4
(28.1,28.15] 5 4 5 5 5 3 3
(28.15,28.2] 4 5 3 5 5 5 3
(28.2, Inf] 4 5 4 5 5 4 4
total 64.00 57.00 63.00 64.00 62.00 56.00 57.00
Indo-Pacific - 1982-1999 --> 2000-2012
(-Inf,26] 5 5 4 5 5 3 3
(26,26.5] 4 4 4 5 5 5 4
(26.5,26.75] 5 4 5 5 4 4 4
(26.75,27] 5 3 5 5 4 3 5
(27,27.25] 5 3 5 5 3 4 5
(27.25,27.5] 5 3 5 5 3 5 4
(27.5,27.75] 4 5 3 5 5 5 3
(27.75,27.85] 5 3 5 5 4 3 5
(27.85,27.95] 5 3 5 4 5 3 5
(27.95,28.05] 5 3 5 5 3 5 4
(28.05,28.1] 4 5 4 5 4 5 4
(28.1, Inf] 4 5 4 5 4 5 4
total 56.00 46.00 54.00 59.00 49.00 50.00 50.00
Indo-Pacific - 2000-2012 --> 2013-2019
(-Inf,26] 5 5 4 5 5 3 3
(26,26.5] 4 4 4 5 5 5 4
(26.5,26.75] 5 4 5 5 3 5 3
(26.75,27] 5 4 5 5 3 3 5
(27,27.25] 5 3 5 5 3 4 5
(27.25,27.5] 5 3 5 5 3 5 4
(27.5,27.75] 4 5 3 5 5 5 3
(27.75,27.85] 5 3 5 5 4 3 5
(27.85,27.95] 5 3 5 3 5 4 5
(27.95,28.05] 5 3 5 5 3 5 4
(28.05,28.1] 4 5 4 5 4 5 4
(28.1, Inf] 4 5 4 5 4 5 4
total 56.00 47.00 54.00 58.00 47.00 52.00 49.00
SO - 1982-1999 --> 2000-2012
(26,26.5] 4 5 3 5 5 3 5
(26.5,26.75] 5 5 4 5 3 3 5
(26.75,27] 5 4 5 5 3 3 5
(27,27.25] 5 3 5 5 3 4 5
(27.25,27.5] 3 5 4 3 5 5 5
(27.5,27.75] 3 5 5 5 3 5 4
(27.75,27.85] 3 5 5 5 3 5 4
(27.85,27.95] 5 3 4 3 5 5 5
(27.95,28.05] 4 5 4 4 5 5 4
(28.05,28.1] 5 4 5 3 5 3 5
(28.1, Inf] 4 5 4 4 5 5 4
(28.1,28.15] 5 4 5 4 5 4 4
(28.15,28.2] 4 5 5 5 4 4 4
(28.2, Inf] 3 5 5 4 5 3 5
total 58.00 63.00 63.00 60.00 59.00 57.00 64.00
SO - 2000-2012 --> 2013-2019
(26,26.5] 3 5 4 5 5 3 5
(26.5,26.75] 5 5 4 5 3 3 5
(26.75,27] 5 4 5 5 3 3 5
(27,27.25] 5 3 5 5 3 4 5
(27.25,27.5] 3 5 4 3 5 5 5
(27.5,27.75] 3 5 5 5 3 5 4
(27.75,27.85] 3 5 5 5 3 5 4
(27.85,27.95] 5 3 4 3 5 5 5
(27.95,28.05] 4 5 4 4 5 5 4
(28.05,28.1] 5 4 5 3 5 3 5
(28.1, Inf] 4 5 4 4 5 5 4
(28.1,28.15] 5 3 5 4 5 3 5
(28.15,28.2] 5 3 5 5 4 3 5
(28.2, Inf] 3 5 5 4 5 3 5
total 58.00 60.00 64.00 60.00 59.00 55.00 66.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
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
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