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.


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.

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.

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.

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

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

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

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

1 Required data

Required are:

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

2 Predictor combinations

Find all possible combinations of following considered predictor variables:

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

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

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

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

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

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

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

# convert to tibble
lm_all <- as_tibble(lm_all)

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

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

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

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

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

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

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

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

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

3 Apply predictor threshold

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

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

This results in a total number of MLR models of:

  • 64

4 Fit all models

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

  • cstar_tref

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

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

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

5 Prepare coeffcients

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

5.1 Formatting

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

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

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

5.2 Predictor selection

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

  • 5

The criterion used to select the best models was:

  • rmse

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

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

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

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

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

5.3 RMSE tables

5.3.1 per model

lm_best %>%
  kable() %>%
  add_header_above() %>%
  kable_styling() %>%
  scroll_box(width = "100%", height = "400px")
basin gamma_slab n_predictors model nr_obs rmse aic resid_max eras rmse_sum aic_sum
Atlantic (-Inf,26] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 158 1.5921923 609.3599345 6.1184078 1982-1999 –> 2000-2012 3.2779521 1430.88036
Atlantic (-Inf,26] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 158 1.0492038 475.5625582 4.3578441 1982-1999 –> 2000-2012 2.3855522 1198.45637
Atlantic (-Inf,26] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 158 1.0480832 477.2248679 4.2990206 1982-1999 –> 2000-2012 2.3837757 1201.91446
Atlantic (-Inf,26] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 158 1.4654213 581.1416848 7.4170479 1982-1999 –> 2000-2012 3.1276914 1394.82476
Atlantic (-Inf,26] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 158 1.4545749 580.7941074 7.0740442 1982-1999 –> 2000-2012 3.1160634 1396.28150
Atlantic (-Inf,26] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 72 1.6825475 293.2516466 3.6000125 2000-2012 –> 2013-2019 3.2747398 902.61158
Atlantic (-Inf,26] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 72 1.2725075 251.0296204 3.6583670 2000-2012 –> 2013-2019 2.3217113 726.59218
Atlantic (-Inf,26] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 72 1.1936365 243.8158001 3.2379740 2000-2012 –> 2013-2019 2.2417197 721.04067
Atlantic (-Inf,26] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 72 1.5761474 281.8447781 5.7132647 2000-2012 –> 2013-2019 3.0415687 862.98646
Atlantic (-Inf,26] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 72 1.5044066 277.1365424 5.0034141 2000-2012 –> 2013-2019 2.9589816 857.93065
Atlantic (26,26.5] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1040 3.5443539 5597.3324051 13.4662233 1982-1999 –> 2000-2012 6.9995324 13651.58188
Atlantic (26,26.5] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 1040 4.2012197 5948.9719164 10.6709624 1982-1999 –> 2000-2012 8.3550943 14558.14368
Atlantic (26,26.5] 4 cstar_tref ~ temp + aou + nitrate + silicate 1040 4.0312783 5863.0858536 14.5182831 1982-1999 –> 2000-2012 8.0820619 14396.26092
Atlantic (26,26.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1040 4.0161192 5857.2495844 16.2265490 1982-1999 –> 2000-2012 8.0668993 14392.42196
Atlantic (26,26.5] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1040 4.2762271 5987.7800473 13.2008700 1982-1999 –> 2000-2012 8.5234582 14666.16238
Atlantic (26,26.5] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 562 3.6054876 3050.3686130 11.1243922 2000-2012 –> 2013-2019 7.1498416 8647.70102
Atlantic (26,26.5] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 562 4.0918975 3190.6127891 12.0407091 2000-2012 –> 2013-2019 8.2931172 9139.58471
Atlantic (26,26.5] 4 cstar_tref ~ temp + aou + nitrate + silicate 562 3.8522413 3122.7752745 13.6121936 2000-2012 –> 2013-2019 7.8835195 8985.86113
Atlantic (26,26.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 562 3.7832546 3104.4640637 14.8635127 2000-2012 –> 2013-2019 7.7993739 8961.71365
Atlantic (26,26.5] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 562 4.2226843 3227.9763371 10.8682033 2000-2012 –> 2013-2019 8.4989114 9215.75638
Atlantic (26.5,26.75] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 1360 3.3727467 7178.2914763 10.7243667 1982-1999 –> 2000-2012 6.4872186 16577.35466
Atlantic (26.5,26.75] 4 cstar_tref ~ sal + aou + silicate + phosphate 1360 3.3303244 7143.8624556 10.1582800 1982-1999 –> 2000-2012 6.4360999 16532.65245
Atlantic (26.5,26.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1360 3.2843587 7108.0590248 10.5247223 1982-1999 –> 2000-2012 6.3355501 16433.70428
Atlantic (26.5,26.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1360 3.3435115 7156.6115389 11.0105664 1982-1999 –> 2000-2012 6.4965794 16602.92486
Atlantic (26.5,26.75] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1360 3.2415884 7072.4053982 10.0491820 1982-1999 –> 2000-2012 6.2239552 16314.22857
Atlantic (26.5,26.75] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 646 3.4450544 3443.3946871 11.3466443 2000-2012 –> 2013-2019 6.8178011 10621.68616
Atlantic (26.5,26.75] 4 cstar_tref ~ sal + aou + silicate + phosphate 646 3.3611477 3411.5375638 10.4660887 2000-2012 –> 2013-2019 6.6914721 10555.40002
Atlantic (26.5,26.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 646 3.3363780 3403.9810106 10.7144118 2000-2012 –> 2013-2019 6.6207366 10512.04004
Atlantic (26.5,26.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 646 3.4054258 3430.4466004 9.5018855 2000-2012 –> 2013-2019 6.7489373 10587.05814
Atlantic (26.5,26.75] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 646 3.3285183 3400.9337733 11.1260032 2000-2012 –> 2013-2019 6.5701066 10473.33917
Atlantic (26.75,27] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 2134 2.1608919 9356.6134164 15.2322922 1982-1999 –> 2000-2012 3.8655083 20016.34118
Atlantic (26.75,27] 4 cstar_tref ~ sal + aou + silicate + phosphate 2134 2.1205644 9276.2096645 17.9450752 1982-1999 –> 2000-2012 3.9632502 20360.71696
Atlantic (26.75,27] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2134 2.0150942 9060.4718086 14.1929877 1982-1999 –> 2000-2012 3.6685116 19555.87536
Atlantic (26.75,27] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2134 2.2560885 9542.6134942 15.8146255 1982-1999 –> 2000-2012 4.0882000 20597.73284
Atlantic (26.75,27] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 2134 2.1557967 9348.5380810 19.2601021 1982-1999 –> 2000-2012 3.6479228 19284.13041
Atlantic (26.75,27] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 1121 2.7579341 5467.7285635 27.8630855 2000-2012 –> 2013-2019 4.9188259 14824.34198
Atlantic (26.75,27] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1121 2.6413897 5372.9263797 21.1206682 2000-2012 –> 2013-2019 5.0063466 15116.67896
Atlantic (26.75,27] 4 cstar_tref ~ sal + aou + silicate + phosphate 1121 2.6063943 5341.0239223 24.8434412 2000-2012 –> 2013-2019 4.7269587 14617.23359
Atlantic (26.75,27] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1121 2.5950428 5333.2381803 25.4741707 2000-2012 –> 2013-2019 4.6101370 14393.70999
Atlantic (26.75,27] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1121 2.3110503 5073.3883143 22.0043753 2000-2012 –> 2013-2019 4.4668470 14421.92640
Atlantic (27,27.25] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 1717 2.3872244 7872.6659685 12.3327089 1982-1999 –> 2000-2012 4.2396475 17015.61877
Atlantic (27,27.25] 4 cstar_tref ~ sal + aou + silicate + phosphate 1717 2.0250912 7307.7159173 9.3807341 1982-1999 –> 2000-2012 3.7447144 16116.95820
Atlantic (27,27.25] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1717 2.0072484 7279.3253051 9.7017253 1982-1999 –> 2000-2012 3.6183821 15798.22804
Atlantic (27,27.25] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1717 2.1502220 7515.6061171 12.8861292 1982-1999 –> 2000-2012 3.8839356 16363.45622
Atlantic (27,27.25] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1717 2.2212354 7627.1850767 15.8645705 1982-1999 –> 2000-2012 4.2651964 17213.53794
Atlantic (27,27.25] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 982 2.3112528 4446.2142611 10.6624255 2000-2012 –> 2013-2019 4.4358419 11920.63745
Atlantic (27,27.25] 4 cstar_tref ~ sal + aou + silicate + phosphate 982 2.4264509 4539.7431240 13.9625358 2000-2012 –> 2013-2019 4.4515421 11847.45904
Atlantic (27,27.25] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 982 2.4198641 4536.4044470 14.2305742 2000-2012 –> 2013-2019 4.4271125 11815.72975
Atlantic (27,27.25] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 982 2.4537027 4561.6780993 14.9835525 2000-2012 –> 2013-2019 4.6724715 12183.04779
Atlantic (27,27.25] 4 cstar_tref ~ sal + silicate + phosphate + phosphate_star 982 2.4701167 4574.7725503 15.1052250 2000-2012 –> 2013-2019 4.7371689 12270.06903
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1621 2.2249112 7206.8812527 13.1379972 1982-1999 –> 2000-2012 4.4362899 17172.95232
Atlantic (27.25,27.5] 3 cstar_tref ~ temp + aou + nitrate 1621 2.4518761 7517.7976652 15.0667238 1982-1999 –> 2000-2012 4.2959606 16662.88251
Atlantic (27.25,27.5] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 1621 2.4080396 7461.3103609 17.0152980 1982-1999 –> 2000-2012 4.2365015 16570.12710
Atlantic (27.25,27.5] 4 cstar_tref ~ temp + aou + nitrate + silicate 1621 2.4361730 7498.9674242 14.2153496 1982-1999 –> 2000-2012 4.1836374 16403.98301
Atlantic (27.25,27.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1621 2.3121875 7331.6238130 16.5099918 1982-1999 –> 2000-2012 4.0506846 16215.49782
Atlantic (27.25,27.5] 3 cstar_tref ~ sal + aou + nitrate 904 2.5315398 4254.7614075 17.3212618 2000-2012 –> 2013-2019 4.8109801 11536.14096
Atlantic (27.25,27.5] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 904 2.5244498 4251.6906907 18.1263712 2000-2012 –> 2013-2019 4.8027876 11533.50176
Atlantic (27.25,27.5] 4 cstar_tref ~ sal + aou + nitrate + silicate 904 2.5311276 4256.4669738 17.4690433 2000-2012 –> 2013-2019 4.7699467 11481.55093
Atlantic (27.25,27.5] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 904 2.5215176 4251.5894594 18.0824185 2000-2012 –> 2013-2019 4.7464288 11458.47071
Atlantic (27.25,27.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 904 2.8515321 4473.9652800 17.6060112 2000-2012 –> 2013-2019 5.1637196 11805.58909
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 2079 2.5437146 9795.9611173 14.8081133 1982-1999 –> 2000-2012 5.0655363 22893.72285
Atlantic (27.5,27.75] 3 cstar_tref ~ sal + aou + phosphate 2079 2.4557608 9645.6458851 17.2009018 1982-1999 –> 2000-2012 4.8210934 22381.80850
Atlantic (27.5,27.75] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 2079 2.4488381 9635.9080260 16.8614853 1982-1999 –> 2000-2012 4.7627250 22251.32328
Atlantic (27.5,27.75] 4 cstar_tref ~ sal + aou + silicate + phosphate 2079 2.4557583 9647.6415993 17.2170993 1982-1999 –> 2000-2012 4.7833793 22296.09080
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2079 2.4487769 9637.8041398 16.7778467 1982-1999 –> 2000-2012 4.7338073 22185.16839
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1115 2.8834250 5539.7557331 17.3414728 2000-2012 –> 2013-2019 5.4271397 15335.71685
Atlantic (27.5,27.75] 3 cstar_tref ~ sal + aou + phosphate 1115 2.8730540 5527.7204339 19.2847791 2000-2012 –> 2013-2019 5.3288148 15173.36632
Atlantic (27.5,27.75] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 1115 2.7822286 5458.0854966 18.3159335 2000-2012 –> 2013-2019 5.2310667 15093.99352
Atlantic (27.5,27.75] 4 cstar_tref ~ sal + aou + silicate + phosphate 1115 2.8564479 5516.7937746 20.6399300 2000-2012 –> 2013-2019 5.3122062 15164.43537
Atlantic (27.5,27.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1115 2.7700300 5450.2866318 19.4859769 2000-2012 –> 2013-2019 5.2188069 15088.09077
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 722 1.6904824 2821.0673683 12.5742519 1982-1999 –> 2000-2012 3.5382366 6730.12190
Atlantic (27.75,27.85] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 722 1.6354065 2771.2383808 13.3517370 1982-1999 –> 2000-2012 3.3492292 6534.12439
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 722 1.6267617 2765.5851457 14.5734851 1982-1999 –> 2000-2012 3.2965550 6480.60424
Atlantic (27.75,27.85] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 722 1.7037266 2830.3363492 11.9220717 1982-1999 –> 2000-2012 3.4207775 6596.82798
Atlantic (27.75,27.85] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 722 1.6657015 2799.7430000 14.2357594 1982-1999 –> 2000-2012 3.3491948 6530.41802
Atlantic (27.75,27.85] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 364 1.8007426 1473.1962392 11.3084900 2000-2012 –> 2013-2019 3.4912705 4292.30242
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 364 1.8006139 1475.1441738 11.2744869 2000-2012 –> 2013-2019 3.4910963 4296.21154
Atlantic (27.75,27.85] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 364 1.7438081 1449.8071331 12.2711012 2000-2012 –> 2013-2019 3.3792145 4221.04551
Atlantic (27.75,27.85] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 364 1.7268302 1444.6845377 13.7518672 2000-2012 –> 2013-2019 3.3535919 4210.26968
Atlantic (27.75,27.85] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 364 1.7767316 1465.4238298 10.2149472 2000-2012 –> 2013-2019 3.4424332 4265.16683
Atlantic (27.85,27.95] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 694 2.8412740 3430.9092372 20.0203261 1982-1999 –> 2000-2012 5.7998788 8259.94399
Atlantic (27.85,27.95] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 694 2.7709000 3398.0976801 20.1418519 1982-1999 –> 2000-2012 5.7244719 8225.85675
Atlantic (27.85,27.95] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 694 2.9086574 3465.4426035 21.2167232 1982-1999 –> 2000-2012 5.8434073 8280.90157
Atlantic (27.85,27.95] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 694 2.7871773 3404.2274548 20.7413865 1982-1999 –> 2000-2012 5.8806820 8319.04731
Atlantic (27.85,27.95] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 694 2.6878487 3355.8593868 20.1028452 1982-1999 –> 2000-2012 5.6412180 8183.48649
Atlantic (27.85,27.95] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 415 3.0542817 2116.4508686 18.0269148 2000-2012 –> 2013-2019 5.8955557 5547.36011
Atlantic (27.85,27.95] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 415 3.0317519 2112.3057046 17.9822303 2000-2012 –> 2013-2019 5.8026519 5510.40338
Atlantic (27.85,27.95] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 415 3.1182188 2135.6463927 18.1103786 2000-2012 –> 2013-2019 6.0268762 5601.08900
Atlantic (27.85,27.95] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 415 3.1908755 2152.7640965 18.7561375 2000-2012 –> 2013-2019 5.9780528 5556.99155
Atlantic (27.85,27.95] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 415 3.0207158 2109.2788576 17.8578045 2000-2012 –> 2013-2019 5.7085645 5465.13824
Atlantic (27.95,28.05] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 723 2.5928883 3441.4940712 13.9368647 1982-1999 –> 2000-2012 4.8700740 7905.86546
Atlantic (27.95,28.05] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 723 2.5849087 3439.0371527 13.5046374 1982-1999 –> 2000-2012 4.8323074 7879.25904
Atlantic (27.95,28.05] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 723 3.0943625 3699.1605551 18.3736306 1982-1999 –> 2000-2012 6.8421502 9155.00761
Atlantic (27.95,28.05] 3 cstar_tref ~ temp + nitrate + phosphate_star 723 3.0944464 3695.1997685 16.3877482 1982-1999 –> 2000-2012 6.7831416 9115.48344
Atlantic (27.95,28.05] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 723 2.8215315 3563.6917628 20.6751203 1982-1999 –> 2000-2012 6.0658508 8731.03808
Atlantic (27.95,28.05] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 386 1.9424018 1619.9708427 14.5879054 2000-2012 –> 2013-2019 4.5352901 5061.46491
Atlantic (27.95,28.05] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 386 1.9183328 1612.3449546 15.4301667 2000-2012 –> 2013-2019 4.5032416 5051.38211
Atlantic (27.95,28.05] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 386 2.8196164 1909.6763894 13.5039640 2000-2012 –> 2013-2019 5.9139789 5608.83694
Atlantic (27.95,28.05] 3 cstar_tref ~ temp + nitrate + phosphate_star 386 2.6073429 1845.2525848 16.7126542 2000-2012 –> 2013-2019 5.7017893 5540.45235
Atlantic (27.95,28.05] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 386 2.4606514 1802.5495198 13.8311274 2000-2012 –> 2013-2019 5.2821830 5366.24128
Atlantic (28.05,28.1] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 436 0.2179047 -79.3496622 1.3640876 1982-1999 –> 2000-2012 0.3337049 -947.26878
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 436 0.2109412 -105.6708211 1.4763185 1982-1999 –> 2000-2012 0.3259563 -979.71285
Atlantic (28.05,28.1] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 436 0.2026977 -140.4320462 1.1989259 1982-1999 –> 2000-2012 0.3184343 -1007.00668
Atlantic (28.05,28.1] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 436 0.2176361 -78.4253437 1.3751833 1982-1999 –> 2000-2012 0.3331923 -946.86274
Atlantic (28.05,28.1] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 436 0.2117907 -104.1663554 1.3461511 1982-1999 –> 2000-2012 0.3300322 -947.17451
Atlantic (28.05,28.1] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 229 0.2500674 27.0745728 0.9325065 2000-2012 –> 2013-2019 0.4621040 -76.08002
Atlantic (28.05,28.1] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 229 0.2472173 23.8245532 1.0717180 2000-2012 –> 2013-2019 0.4581585 -81.84627
Atlantic (28.05,28.1] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 229 0.2412331 12.6016075 1.0576928 2000-2012 –> 2013-2019 0.4439308 -127.83044
Atlantic (28.05,28.1] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 229 0.2589000 44.9724007 1.3338851 2000-2012 –> 2013-2019 0.4765361 -33.45294
Atlantic (28.05,28.1] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 229 0.2498995 26.7669599 1.1561502 2000-2012 –> 2013-2019 0.4616902 -77.39940
Atlantic (28.1,28.15] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 541 0.1480631 -519.4550915 1.2063423 1982-1999 –> 2000-2012 0.3108018 -1121.50164
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 541 0.1394443 -582.3456119 1.0824605 1982-1999 –> 2000-2012 0.2743304 -1472.97554
Atlantic (28.1,28.15] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 541 0.1576027 -451.8958141 1.4966383 1982-1999 –> 2000-2012 0.2858731 -1422.37448
Atlantic (28.1,28.15] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 541 0.1549911 -467.9756571 1.5279244 1982-1999 –> 2000-2012 0.2748469 -1541.48871
Atlantic (28.1,28.15] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 541 0.1566336 -456.5697948 1.2740811 1982-1999 –> 2000-2012 0.3008900 -1243.23369
Atlantic (28.1,28.15] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 288 0.1409161 -299.4157065 0.6463000 2000-2012 –> 2013-2019 0.2889791 -818.87080
Atlantic (28.1,28.15] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 288 0.1282569 -351.6339749 0.4661753 2000-2012 –> 2013-2019 0.2677013 -933.97959
Atlantic (28.1,28.15] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 288 0.1624875 -215.3720860 1.0609876 2000-2012 –> 2013-2019 0.3174787 -683.34774
Atlantic (28.1,28.15] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 288 0.1592978 -228.7916378 0.4962206 2000-2012 –> 2013-2019 0.3237730 -634.50572
Atlantic (28.1,28.15] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 288 0.1506654 -258.8829283 0.6799218 2000-2012 –> 2013-2019 0.3072990 -715.45272
Atlantic (28.15,28.2] 3 cstar_tref ~ sal + aou + nitrate 1155 0.1574471 -982.6698392 0.9830553 1982-1999 –> 2000-2012 0.3636269 -1527.46540
Atlantic (28.15,28.2] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 1155 0.1550542 -1016.0477142 0.9483180 1982-1999 –> 2000-2012 0.3595782 -1586.79067
Atlantic (28.15,28.2] 4 cstar_tref ~ sal + aou + nitrate + silicate 1155 0.1567453 -990.9892289 0.9661276 1982-1999 –> 2000-2012 0.3491822 -1772.87154
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1155 0.1338912 -1353.0327195 0.6859787 1982-1999 –> 2000-2012 0.3082551 -2474.74523
Atlantic (28.15,28.2] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1155 0.1703992 -796.0542748 0.7124229 1982-1999 –> 2000-2012 0.3927223 -1075.57265
Atlantic (28.15,28.2] 3 cstar_tref ~ sal + aou + nitrate 636 0.2156021 -136.7663848 1.9725799 2000-2012 –> 2013-2019 0.3730492 -1119.43622
Atlantic (28.15,28.2] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 636 0.2000269 -230.1440140 1.4705162 2000-2012 –> 2013-2019 0.3550811 -1246.19173
Atlantic (28.15,28.2] 4 cstar_tref ~ sal + aou + nitrate + silicate 636 0.1903257 -293.3813460 1.3866890 2000-2012 –> 2013-2019 0.3470711 -1284.37057
Atlantic (28.15,28.2] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 636 0.1879213 -307.5531446 1.4575343 2000-2012 –> 2013-2019 0.3218125 -1660.58586
Atlantic (28.15,28.2] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 636 0.2394319 0.5830451 1.9671533 2000-2012 –> 2013-2019 0.4098310 -795.47123
Atlantic (28.2, Inf] 3 cstar_tref ~ sal + aou + nitrate 955 0.1973403 -379.4239149 1.9207165 1982-1999 –> 2000-2012 0.4817501 31.01892
Atlantic (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 955 0.1541931 -848.6675542 0.7135736 1982-1999 –> 2000-2012 0.3663878 -1162.05962
Atlantic (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + silicate 955 0.1744896 -612.4775343 0.9944731 1982-1999 –> 2000-2012 0.4081182 -687.41598
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 955 0.1427395 -994.0896712 0.8634543 1982-1999 –> 2000-2012 0.3506627 -1355.87253
Atlantic (28.2, Inf] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 955 0.1862717 -485.6754946 1.3306437 1982-1999 –> 2000-2012 0.5404322 472.27744
Atlantic (28.2, Inf] 3 cstar_tref ~ sal + aou + nitrate 552 0.2864800 196.4127031 2.5941915 2000-2012 –> 2013-2019 0.4838204 -183.01121
Atlantic (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 552 0.2005924 -195.0462133 1.3544749 2000-2012 –> 2013-2019 0.3547854 -1043.71377
Atlantic (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + silicate 552 0.2343110 -23.5130276 1.6064160 2000-2012 –> 2013-2019 0.4088006 -635.99056
Atlantic (28.2, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 552 0.1875683 -267.1597642 0.8704042 2000-2012 –> 2013-2019 0.3303078 -1261.24944
Atlantic (28.2, Inf] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 552 0.2982717 244.9436937 2.0512928 2000-2012 –> 2013-2019 0.4845435 -240.73180
Indo-Pacific (-Inf,26] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 5288 6.2360939 34376.5179250 29.9685673 1982-1999 –> 2000-2012 12.7664702 81783.00523
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 5288 6.1928478 34304.9199551 29.5273536 1982-1999 –> 2000-2012 12.6725136 81601.29062
Indo-Pacific (-Inf,26] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5288 7.1987438 35896.7294426 31.3718077 1982-1999 –> 2000-2012 14.7274968 85351.27007
Indo-Pacific (-Inf,26] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 5288 7.2933052 36034.7494052 29.1665119 1982-1999 –> 2000-2012 14.8782495 85596.23226
Indo-Pacific (-Inf,26] 4 cstar_tref ~ temp + silicate + phosphate + phosphate_star 5288 7.3907730 36173.1510635 28.4925990 1982-1999 –> 2000-2012 15.1057280 85977.06154
Indo-Pacific (-Inf,26] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 2815 5.9530030 18043.9572913 24.8937268 2000-2012 –> 2013-2019 12.1890969 52420.47522
Indo-Pacific (-Inf,26] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2815 5.8787943 17975.3339916 24.1274152 2000-2012 –> 2013-2019 12.0716422 52280.25395
Indo-Pacific (-Inf,26] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2815 6.9518302 18919.2219176 25.9076786 2000-2012 –> 2013-2019 14.1505740 54815.95136
Indo-Pacific (-Inf,26] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 2815 7.0008395 18958.7732280 25.4625103 2000-2012 –> 2013-2019 14.2941447 54993.52263
Indo-Pacific (-Inf,26] 4 cstar_tref ~ temp + silicate + phosphate + phosphate_star 2815 7.1346964 19063.4033807 24.4397815 2000-2012 –> 2013-2019 14.5254694 55236.55444
Indo-Pacific (26,26.5] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 5428 4.8772756 32618.2708412 41.7349156 1982-1999 –> 2000-2012 9.5093971 77476.34842
Indo-Pacific (26,26.5] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 5428 4.5387228 31839.2778896 43.7366220 1982-1999 –> 2000-2012 8.7811733 75364.37494
Indo-Pacific (26,26.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5428 4.7468521 32326.0177685 49.2620011 1982-1999 –> 2000-2012 9.1449428 76398.47442
Indo-Pacific (26,26.5] 4 cstar_tref ~ temp + aou + silicate + phosphate 5428 4.8836514 32632.4530772 40.4064523 1982-1999 –> 2000-2012 9.4934799 77417.23967
Indo-Pacific (26,26.5] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 5428 4.8704336 32605.0310057 41.2998501 1982-1999 –> 2000-2012 9.4796435 77389.77888
Indo-Pacific (26,26.5] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 2876 4.8221258 17222.8663824 40.3488812 2000-2012 –> 2013-2019 9.6994014 49841.13722
Indo-Pacific (26,26.5] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2876 4.4815941 16803.6126080 46.1900332 2000-2012 –> 2013-2019 9.0203169 48642.89050
Indo-Pacific (26,26.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2876 4.5772124 16925.0450303 45.7195416 2000-2012 –> 2013-2019 9.3240644 49251.06280
Indo-Pacific (26,26.5] 4 cstar_tref ~ temp + aou + silicate + phosphate 2876 4.8221384 17222.8814144 36.9760613 2000-2012 –> 2013-2019 9.7057898 49855.33449
Indo-Pacific (26,26.5] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 2876 4.8127496 17213.6712211 38.8375973 2000-2012 –> 2013-2019 9.6831832 49818.70223
Indo-Pacific (26.5,26.75] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 4269 4.0546161 24078.8678518 22.3581891 1982-1999 –> 2000-2012 8.0016719 57191.73449
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 4269 3.9438120 23844.2949199 20.4830181 1982-1999 –> 2000-2012 7.6993981 56369.61516
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 4269 4.0277915 24024.1942408 33.2108036 1982-1999 –> 2000-2012 8.0184863 57269.42280
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 4269 3.9458514 23848.7088413 23.7921959 1982-1999 –> 2000-2012 7.8106187 56713.78794
Indo-Pacific (26.5,26.75] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 4269 4.0285176 24023.7333576 32.9352042 1982-1999 –> 2000-2012 8.0934769 57485.56713
Indo-Pacific (26.5,26.75] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 2324 4.4254119 13520.4911820 24.9112112 2000-2012 –> 2013-2019 8.4800280 37599.35903
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2324 4.3199526 13410.3863975 24.6916724 2000-2012 –> 2013-2019 8.2637646 37254.68132
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2324 4.3368286 13428.5085435 32.2595908 2000-2012 –> 2013-2019 8.3646201 37452.70278
Indo-Pacific (26.5,26.75] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 2324 4.3137865 13403.7473329 23.1641121 2000-2012 –> 2013-2019 8.2596379 37252.45617
Indo-Pacific (26.5,26.75] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 2324 4.3478487 13438.3043881 34.5430589 2000-2012 –> 2013-2019 8.3763663 37462.03775
Indo-Pacific (26.75,27] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 4557 4.5583246 26769.7349207 19.0779671 1982-1999 –> 2000-2012 8.5369273 62782.45338
Indo-Pacific (26.75,27] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 4557 4.3904749 26427.7981926 21.4038278 1982-1999 –> 2000-2012 8.0962279 61527.03030
Indo-Pacific (26.75,27] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 4557 4.3862266 26420.9752077 21.0605286 1982-1999 –> 2000-2012 8.0872756 61505.87489
Indo-Pacific (26.75,27] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 4557 4.1598161 25935.9472159 23.1614170 1982-1999 –> 2000-2012 7.7437292 60605.32053
Indo-Pacific (26.75,27] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 4557 4.1362281 25886.1196857 21.9180570 1982-1999 –> 2000-2012 7.7019358 60492.01106
Indo-Pacific (26.75,27] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 2460 4.9804504 14892.3376484 18.4452760 2000-2012 –> 2013-2019 9.5387750 41662.07257
Indo-Pacific (26.75,27] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 2460 4.7927943 14703.3765563 19.7027300 2000-2012 –> 2013-2019 9.1832692 41131.17475
Indo-Pacific (26.75,27] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2460 4.7910605 14703.5963954 19.4544536 2000-2012 –> 2013-2019 9.1772872 41124.57160
Indo-Pacific (26.75,27] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 2460 4.6128201 14515.0674866 16.6134065 2000-2012 –> 2013-2019 8.7726362 40451.01470
Indo-Pacific (26.75,27] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 2460 4.5994232 14502.7576801 16.6240534 2000-2012 –> 2013-2019 8.7356513 40388.87737
Indo-Pacific (27,27.25] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 4267 4.2336742 24436.3826880 43.0543060 1982-1999 –> 2000-2012 7.6345327 54932.35764
Indo-Pacific (27,27.25] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 4267 4.3158154 24600.3724211 51.0072982 1982-1999 –> 2000-2012 7.7016451 55045.26427
Indo-Pacific (27,27.25] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 4267 4.2143064 24399.2526775 44.7114638 1982-1999 –> 2000-2012 7.5507703 54676.73929
Indo-Pacific (27,27.25] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 4267 3.7387791 23375.5116467 42.2525000 1982-1999 –> 2000-2012 6.7614943 52511.94640
Indo-Pacific (27,27.25] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 4267 3.4041802 22577.4068316 30.9255112 1982-1999 –> 2000-2012 6.2015505 50822.23662
Indo-Pacific (27,27.25] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 2344 5.0787920 14282.3280763 32.6832487 2000-2012 –> 2013-2019 9.3124661 38718.71076
Indo-Pacific (27,27.25] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2344 5.0701758 14276.3681434 33.6277311 2000-2012 –> 2013-2019 9.2844822 38675.62082
Indo-Pacific (27,27.25] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2344 5.2740156 14461.1529911 34.1511699 2000-2012 –> 2013-2019 9.5983174 39080.28973
Indo-Pacific (27,27.25] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 2344 4.5558378 13772.9113115 35.4736376 2000-2012 –> 2013-2019 8.2946169 37148.42296
Indo-Pacific (27,27.25] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 2344 3.9533412 13109.9262055 24.4569335 2000-2012 –> 2013-2019 7.3575214 35687.33304
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 3520 3.0170713 17777.5047944 30.6141183 1982-1999 –> 2000-2012 5.8968935 41156.38222
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 3520 3.1346725 18046.7012210 29.1767786 1982-1999 –> 2000-2012 6.0487445 41537.11519
Indo-Pacific (27.25,27.5] 4 cstar_tref ~ sal + silicate + phosphate + phosphate_star 3520 3.1493898 18077.6766494 29.8455275 1982-1999 –> 2000-2012 6.0757661 41605.84076
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 3520 3.0094328 17759.6586184 29.1482995 1982-1999 –> 2000-2012 5.9385710 41298.72235
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 3520 2.8774108 17443.8388755 30.1555714 1982-1999 –> 2000-2012 5.6901545 40600.37457
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1879 3.4636269 10014.9955686 28.8274843 2000-2012 –> 2013-2019 6.4806982 27792.50036
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1879 3.6867091 10249.5622039 27.5141850 2000-2012 –> 2013-2019 6.8213816 28296.26342
Indo-Pacific (27.25,27.5] 4 cstar_tref ~ sal + silicate + phosphate + phosphate_star 1879 3.6875866 10248.4565706 27.7368313 2000-2012 –> 2013-2019 6.8369764 28326.13322
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1879 3.5497395 10107.2843807 33.8094442 2000-2012 –> 2013-2019 6.5591723 27866.94300
Indo-Pacific (27.25,27.5] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1879 3.3526141 9892.5752799 28.9794819 2000-2012 –> 2013-2019 6.2300248 27336.41416
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 3170 1.8542474 12924.8866415 13.2838580 1982-1999 –> 2000-2012 3.4799789 30094.38329
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 3170 1.8834035 13023.8007420 9.5549910 1982-1999 –> 2000-2012 3.5338959 30329.43049
Indo-Pacific (27.5,27.75] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 3170 1.9027986 13086.7555124 12.4473296 1982-1999 –> 2000-2012 3.5535107 30391.58410
Indo-Pacific (27.5,27.75] 4 cstar_tref ~ sal + silicate + phosphate + phosphate_star 3170 1.9247016 13159.3181099 14.1925496 1982-1999 –> 2000-2012 3.5760715 30467.73487
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 3170 1.9773131 13332.2950506 12.8773083 1982-1999 –> 2000-2012 3.6743555 30888.41122
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1683 2.1594104 7381.4124538 13.7085041 2000-2012 –> 2013-2019 4.0136578 20306.29910
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1683 2.1483716 7364.1614621 12.8217994 2000-2012 –> 2013-2019 4.0317751 20387.96220
Indo-Pacific (27.5,27.75] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 1683 2.2370458 7498.3027878 13.7718692 2000-2012 –> 2013-2019 4.1398444 20585.05830
Indo-Pacific (27.5,27.75] 4 cstar_tref ~ sal + silicate + phosphate + phosphate_star 1683 2.2955622 7585.2182760 14.9967179 2000-2012 –> 2013-2019 4.2202638 20744.53639
Indo-Pacific (27.5,27.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1683 2.2907164 7580.1053726 13.2792998 2000-2012 –> 2013-2019 4.2680295 20912.40042
Indo-Pacific (27.75,27.85] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 1216 1.4393241 4348.5287948 8.0041666 1982-1999 –> 2000-2012 2.7883962 10161.69710
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1216 1.4268879 4329.4241893 7.6562558 1982-1999 –> 2000-2012 2.7551730 10092.16905
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1216 1.4429788 4356.6962523 8.1238690 1982-1999 –> 2000-2012 2.7877140 10160.99422
Indo-Pacific (27.75,27.85] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 1216 1.4426877 4354.2054884 6.9725516 1982-1999 –> 2000-2012 2.7872451 10156.05712
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1216 1.3124934 4126.1891604 6.8457618 1982-1999 –> 2000-2012 2.5618830 9682.20926
Indo-Pacific (27.75,27.85] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 691 1.6407154 2657.2460163 13.7031008 2000-2012 –> 2013-2019 3.0800396 7005.77481
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 691 1.6140386 2636.5909926 13.0991822 2000-2012 –> 2013-2019 3.0409264 6966.01518
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 691 1.6406123 2659.1591558 13.9456772 2000-2012 –> 2013-2019 3.0835911 7015.85541
Indo-Pacific (27.75,27.85] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 691 1.6225333 2641.8453942 11.4420005 2000-2012 –> 2013-2019 3.0652209 6996.05088
Indo-Pacific (27.75,27.85] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 691 1.5393786 2571.1385959 10.1972442 2000-2012 –> 2013-2019 2.8518720 6697.32776
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1797 1.4125952 6355.1355196 11.3038103 1982-1999 –> 2000-2012 2.6595697 14669.08780
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1797 1.3767808 6262.8394080 10.8206192 1982-1999 –> 2000-2012 2.6534586 14695.95693
Indo-Pacific (27.85,27.95] 4 cstar_tref ~ temp + aou + nitrate + silicate 1797 1.3683349 6238.7241576 9.2082326 1982-1999 –> 2000-2012 2.6475550 14679.91157
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1797 1.3568501 6210.4314395 9.3793503 1982-1999 –> 2000-2012 2.5976136 14499.10798
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1797 1.3914252 6300.8660064 9.9533251 1982-1999 –> 2000-2012 2.6565451 14687.94739
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 987 1.6494324 3802.8359632 12.1825689 2000-2012 –> 2013-2019 3.0620277 10157.97148
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 987 1.5919784 3732.8502661 10.5660006 2000-2012 –> 2013-2019 2.9687592 9995.68967
Indo-Pacific (27.85,27.95] 4 cstar_tref ~ temp + aou + nitrate + silicate 987 1.6038431 3745.5076106 8.9600325 2000-2012 –> 2013-2019 2.9721781 9984.23177
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 987 1.5960835 3737.9339228 9.2774322 2000-2012 –> 2013-2019 2.9529336 9948.36536
Indo-Pacific (27.85,27.95] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 987 1.6427665 3794.8421339 10.0005178 2000-2012 –> 2013-2019 3.0341917 10095.70814
Indo-Pacific (27.95,28.05] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 1162 0.6613106 2348.5655420 3.2626451 1982-1999 –> 2000-2012 1.2578203 5498.69125
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1162 0.5801738 2046.3633101 3.0668291 1982-1999 –> 2000-2012 1.1222674 4865.79437
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1162 0.6596354 2344.6710734 5.9778740 1982-1999 –> 2000-2012 1.2522936 5474.26713
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1162 0.5749226 2025.2329788 2.9259808 1982-1999 –> 2000-2012 1.1074497 4782.73873
Indo-Pacific (27.95,28.05] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 1162 0.6598787 2343.5280283 5.7094931 1982-1999 –> 2000-2012 1.2528152 5472.75732
Indo-Pacific (27.95,28.05] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 648 0.6657190 1323.6179563 3.0901877 2000-2012 –> 2013-2019 1.3270296 3672.18350
Indo-Pacific (27.95,28.05] 4 cstar_tref ~ sal + aou + silicate + phosphate 648 0.6534666 1299.5431965 3.1086636 2000-2012 –> 2013-2019 1.2877207 3551.02612
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 648 0.6051146 1201.9150760 2.6426945 2000-2012 –> 2013-2019 1.1852884 3248.27839
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 648 0.6715663 1336.9516306 3.3640626 2000-2012 –> 2013-2019 1.3312017 3681.62270
Indo-Pacific (27.95,28.05] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 648 0.5922331 1174.0283013 2.7169734 2000-2012 –> 2013-2019 1.1671557 3199.26128
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1015 0.4795310 1402.5034846 3.1382674 1982-1999 –> 2000-2012 0.9490005 3430.04886
Indo-Pacific (28.05,28.1] 4 cstar_tref ~ sal + aou + silicate + phosphate 1015 0.4621065 1325.3667754 2.5240688 1982-1999 –> 2000-2012 0.9190226 3268.57196
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1015 0.4567175 1303.5539755 2.4275171 1982-1999 –> 2000-2012 0.9128334 3243.43396
Indo-Pacific (28.05,28.1] 4 cstar_tref ~ sal + nitrate + silicate + phosphate_star 1015 0.4832957 1416.3783333 2.5369659 1982-1999 –> 2000-2012 0.9535442 3446.96091
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1015 0.4546336 1294.2705038 2.4086220 1982-1999 –> 2000-2012 0.9095606 3226.22099
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 565 0.4790009 785.6609647 2.5225697 2000-2012 –> 2013-2019 0.9585320 2188.16445
Indo-Pacific (28.05,28.1] 4 cstar_tref ~ sal + aou + silicate + phosphate 565 0.4782438 781.8733163 2.7103826 2000-2012 –> 2013-2019 0.9403503 2107.24009
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 565 0.4757060 777.8611535 2.6220126 2000-2012 –> 2013-2019 0.9324235 2081.41513
Indo-Pacific (28.05,28.1] 4 cstar_tref ~ sal + nitrate + silicate + phosphate_star 565 0.4904495 810.3512236 2.7325905 2000-2012 –> 2013-2019 0.9737452 2226.72956
Indo-Pacific (28.05,28.1] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 565 0.4756293 777.6789761 2.6025948 2000-2012 –> 2013-2019 0.9302630 2071.94948
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 9318 0.7861340 21972.9954474 3.7460043 1982-1999 –> 2000-2012 1.5770507 52723.91650
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 9318 0.7764916 21742.9989329 4.8420893 1982-1999 –> 2000-2012 1.5383986 51524.14490
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 9318 0.6753077 19141.0900046 2.6646187 1982-1999 –> 2000-2012 1.3451269 45579.18684
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 9318 0.7744993 21695.1218258 4.9051251 1982-1999 –> 2000-2012 1.5417070 51656.19560
Indo-Pacific (28.1, Inf] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 9318 0.7901431 22065.7914418 3.2330378 1982-1999 –> 2000-2012 1.5802157 52787.00533
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 5093 0.7791371 11925.2057059 4.2437545 2000-2012 –> 2013-2019 1.5652711 33898.20115
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 5093 0.7763424 11888.6031698 4.2785787 2000-2012 –> 2013-2019 1.5528339 33631.60210
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5093 0.6800660 10539.9378299 3.0769840 2000-2012 –> 2013-2019 1.3553737 29681.02783
Indo-Pacific (28.1, Inf] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 5093 0.7692867 11795.6056405 4.3360410 2000-2012 –> 2013-2019 1.5437859 33490.72747
Indo-Pacific (28.1, Inf] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 5093 0.7933509 12107.3541598 2.7897695 2000-2012 –> 2013-2019 1.5834939 34173.14560
SO (26,26.5] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 103 2.7820684 517.0794443 11.7335300 1982-1999 –> 2000-2012 6.3296968 1223.86565
SO (26,26.5] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 103 1.0023937 304.7938596 4.5747847 1982-1999 –> 2000-2012 2.3237890 754.78155
SO (26,26.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 103 0.9768158 301.4691489 3.5663094 1982-1999 –> 2000-2012 2.1945946 732.38866
SO (26,26.5] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 103 2.8670780 521.2797745 11.5979108 1982-1999 –> 2000-2012 6.3550799 1221.69282
SO (26,26.5] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 103 2.8011962 518.4909224 11.6880130 1982-1999 –> 2000-2012 6.2394349 1217.19657
SO (26,26.5] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 47 3.2156777 255.1758056 8.1981778 2000-2012 –> 2013-2019 6.0144042 771.48503
SO (26,26.5] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 47 3.0953135 253.5898064 8.2440313 2000-2012 –> 2013-2019 5.8773819 770.66925
SO (26,26.5] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 47 1.1019644 154.5070979 2.3246978 2000-2012 –> 2013-2019 2.1043582 459.30096
SO (26,26.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 47 1.0381247 150.8973184 2.2629736 2000-2012 –> 2013-2019 2.0149405 452.36647
SO (26,26.5] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 47 3.2135703 257.1141819 7.9236145 2000-2012 –> 2013-2019 6.0147664 775.60510
SO (26.5,26.75] 4 cstar_tref ~ aou + nitrate + silicate + phosphate_star 312 3.0634697 1596.0156948 14.8316212 1982-1999 –> 2000-2012 6.6583175 3874.70545
SO (26.5,26.75] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 312 3.1102237 1605.4670996 13.3865164 1982-1999 –> 2000-2012 6.5963633 3858.36326
SO (26.5,26.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 312 3.0429219 1593.8161999 14.2419897 1982-1999 –> 2000-2012 6.5266828 3848.13901
SO (26.5,26.75] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 312 1.6818786 1221.8423319 7.4624349 1982-1999 –> 2000-2012 3.5967902 2971.47474
SO (26.5,26.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 312 1.6697360 1219.3209394 7.1133689 1982-1999 –> 2000-2012 3.5488640 2955.10787
SO (26.5,26.75] 4 cstar_tref ~ aou + nitrate + silicate + phosphate_star 133 2.8586117 668.8270445 10.7899231 2000-2012 –> 2013-2019 5.9220813 2264.84274
SO (26.5,26.75] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 133 2.9251771 674.9500825 10.0430783 2000-2012 –> 2013-2019 6.0354007 2280.41718
SO (26.5,26.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 133 2.8346133 668.5845216 10.4915004 2000-2012 –> 2013-2019 5.8775352 2262.40072
SO (26.5,26.75] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 133 1.6548865 523.4304747 4.9197792 2000-2012 –> 2013-2019 3.3367651 1745.27281
SO (26.5,26.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 133 1.6547816 525.4136173 4.9257339 2000-2012 –> 2013-2019 3.3245176 1744.73456
SO (26.75,27] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 1061 2.4677656 4939.8180549 13.5857513 1982-1999 –> 2000-2012 5.1625431 12254.88820
SO (26.75,27] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1061 2.4672141 4941.3437400 13.4599518 1982-1999 –> 2000-2012 5.1618710 12258.27828
SO (26.75,27] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 1061 2.3658852 4850.3526345 11.7790894 1982-1999 –> 2000-2012 4.8536773 11923.26504
SO (26.75,27] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1061 2.3191394 4810.0058453 10.6756905 1982-1999 –> 2000-2012 4.7542973 11820.12478
SO (26.75,27] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 1061 2.9578329 5324.2050702 10.9471230 1982-1999 –> 2000-2012 6.0140839 13020.67143
SO (26.75,27] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 550 2.5663729 2609.5753272 17.6995339 2000-2012 –> 2013-2019 5.0341385 7549.39338
SO (26.75,27] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 550 2.5605544 2609.0785756 17.0705169 2000-2012 –> 2013-2019 5.0277685 7550.42232
SO (26.75,27] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 550 2.3957372 2533.8924621 15.4998931 2000-2012 –> 2013-2019 4.7616224 7384.24510
SO (26.75,27] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 550 2.2748295 2478.9280053 12.6508529 2000-2012 –> 2013-2019 4.5939689 7288.93385
SO (26.75,27] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 550 3.0086847 2784.4856997 12.4224158 2000-2012 –> 2013-2019 5.9665176 8108.69077
SO (27,27.25] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 2513 1.9002603 10370.2312203 19.2728362 1982-1999 –> 2000-2012 3.7237639 24269.67624
SO (27,27.25] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2513 1.8530732 10245.8502101 18.9619811 1982-1999 –> 2000-2012 3.6704902 24124.30576
SO (27,27.25] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 2513 1.7159114 9857.3454571 16.6672621 1982-1999 –> 2000-2012 3.3256604 22899.48106
SO (27,27.25] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 2513 1.6615541 9697.5535747 16.3323879 1982-1999 –> 2000-2012 3.2687331 22730.70242
SO (27,27.25] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 2513 1.9463574 10490.6980852 19.9217964 1982-1999 –> 2000-2012 3.8581729 24715.33421
SO (27,27.25] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 1286 2.1008436 5570.8058048 13.9602984 2000-2012 –> 2013-2019 4.0011039 15941.03703
SO (27,27.25] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1286 2.0262434 5479.8139370 13.8567921 2000-2012 –> 2013-2019 3.8793165 15725.66415
SO (27,27.25] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 1286 1.9277877 5349.7014449 11.5137595 2000-2012 –> 2013-2019 3.6436991 15207.04690
SO (27,27.25] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1286 1.8504129 5246.3413614 11.4263441 2000-2012 –> 2013-2019 3.5119670 14943.89494
SO (27,27.25] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 1286 2.1184198 5592.2342269 14.6526181 2000-2012 –> 2013-2019 4.0647771 16082.93231
SO (27.25,27.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2731 1.6825482 10606.1722440 8.7665619 1982-1999 –> 2000-2012 3.1211793 24182.42563
SO (27.25,27.5] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 2731 1.6893569 10628.2307746 10.5492429 1982-1999 –> 2000-2012 3.1276338 24202.61058
SO (27.25,27.5] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 2731 1.6893649 10626.2566473 10.5431309 1982-1999 –> 2000-2012 3.1306658 24214.61540
SO (27.25,27.5] 4 cstar_tref ~ temp + nitrate + silicate + phosphate_star 2731 1.6834191 10606.9989239 8.1155062 1982-1999 –> 2000-2012 3.1248601 24196.09752
SO (27.25,27.5] 4 cstar_tref ~ temp + silicate + phosphate + phosphate_star 2731 1.6897063 10627.3602422 10.5063777 1982-1999 –> 2000-2012 3.1326599 24224.43821
SO (27.25,27.5] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1453 1.8604040 5941.4617956 9.5764858 2000-2012 –> 2013-2019 3.5429522 16547.63404
SO (27.25,27.5] 4 cstar_tref ~ temp + aou + silicate + phosphate 1453 1.8580136 5935.7255654 9.5664542 2000-2012 –> 2013-2019 3.5508334 16573.14093
SO (27.25,27.5] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1453 1.8579855 5937.6815988 9.5635309 2000-2012 –> 2013-2019 3.5473425 16565.91237
SO (27.25,27.5] 4 cstar_tref ~ temp + aou + silicate + phosphate_star 1453 1.8608548 5940.1658157 9.5050972 2000-2012 –> 2013-2019 3.5502197 16566.42246
SO (27.25,27.5] 4 cstar_tref ~ temp + silicate + phosphate + phosphate_star 1453 1.8592539 5937.6646737 9.5274690 2000-2012 –> 2013-2019 3.5489602 16565.02492
SO (27.5,27.75] 4 cstar_tref ~ aou + nitrate + silicate + phosphate_star 3391 1.4983094 12377.4572107 8.4631334 1982-1999 –> 2000-2012 2.7575978 27383.27911
SO (27.5,27.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 3391 1.4809706 12300.5166921 8.2778025 1982-1999 –> 2000-2012 2.7401872 27307.82001
SO (27.5,27.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 3391 1.5320138 12530.3276953 9.1542418 1982-1999 –> 2000-2012 2.8147095 27705.56033
SO (27.5,27.75] 4 cstar_tref ~ temp + aou + nitrate + silicate 3391 1.5314112 12525.6595815 8.9583230 1982-1999 –> 2000-2012 2.7897750 27524.80452
SO (27.5,27.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 3391 1.4865171 12325.8693341 8.1453800 1982-1999 –> 2000-2012 2.7423144 27308.45606
SO (27.5,27.75] 4 cstar_tref ~ aou + nitrate + silicate + phosphate_star 1821 1.7424835 7202.2181577 10.3871064 2000-2012 –> 2013-2019 3.2407928 19579.67537
SO (27.5,27.75] 4 cstar_tref ~ sal + aou + nitrate + silicate 1821 1.7135136 7141.1585933 9.9489584 2000-2012 –> 2013-2019 3.2205821 19558.14808
SO (27.5,27.75] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1821 1.6982942 7110.6659963 9.9919662 2000-2012 –> 2013-2019 3.1792648 19411.18269
SO (27.5,27.75] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1821 1.7554240 7231.1654646 10.3901005 2000-2012 –> 2013-2019 3.2874378 19761.49316
SO (27.5,27.75] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1821 1.6931774 7099.6764854 9.9569404 2000-2012 –> 2013-2019 3.1796946 19425.54582
SO (27.75,27.85] 4 cstar_tref ~ aou + nitrate + silicate + phosphate_star 1415 1.4036591 4987.1994177 11.2501595 1982-1999 –> 2000-2012 2.6328581 10866.53706
SO (27.75,27.85] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1415 1.3799763 4941.0437367 10.6092986 1982-1999 –> 2000-2012 2.6003361 10796.32800
SO (27.75,27.85] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1415 1.4740467 5127.6684405 13.4027185 1982-1999 –> 2000-2012 2.7633662 11181.39033
SO (27.75,27.85] 4 cstar_tref ~ temp + aou + nitrate + silicate 1415 1.5056329 5185.6697430 13.4580336 1982-1999 –> 2000-2012 2.7650010 11152.54027
SO (27.75,27.85] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1415 1.3062806 4785.7263346 8.6276320 1982-1999 –> 2000-2012 2.5318731 10656.45667
SO (27.75,27.85] 4 cstar_tref ~ aou + nitrate + silicate + phosphate_star 722 1.5404965 2684.9065454 8.3689795 2000-2012 –> 2013-2019 2.9441556 7672.10596
SO (27.75,27.85] 4 cstar_tref ~ sal + aou + nitrate + silicate 722 1.5764019 2718.1765705 8.5554720 2000-2012 –> 2013-2019 3.0324771 7809.12992
SO (27.75,27.85] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 722 1.5043391 2652.6099150 8.2718276 2000-2012 –> 2013-2019 2.8843154 7593.65365
SO (27.75,27.85] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 722 1.5913038 2733.7627824 8.6142494 2000-2012 –> 2013-2019 3.0653505 7861.43122
SO (27.75,27.85] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 722 1.4200481 2569.3447020 8.5861197 2000-2012 –> 2013-2019 2.7263287 7355.07104
SO (27.85,27.95] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1399 1.3661540 4857.1645441 8.5867618 1982-1999 –> 2000-2012 2.7006066 11199.01016
SO (27.85,27.95] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 1399 1.5214807 5156.4658834 8.1358516 1982-1999 –> 2000-2012 2.9663087 11790.82404
SO (27.85,27.95] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1399 1.4396129 5003.7092447 9.8030845 1982-1999 –> 2000-2012 2.8546010 11562.72636
SO (27.85,27.95] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 1399 1.4579776 5037.1766853 8.4189791 1982-1999 –> 2000-2012 2.9684893 11836.29687
SO (27.85,27.95] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1399 1.4352405 4995.1981575 9.3064311 1982-1999 –> 2000-2012 2.9394913 11780.92549
SO (27.85,27.95] 4 cstar_tref ~ aou + nitrate + silicate + phosphate_star 709 1.5653688 2659.4910609 8.8715048 2000-2012 –> 2013-2019 3.0679953 7781.06745
SO (27.85,27.95] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 709 1.5223861 2622.0102956 7.8651764 2000-2012 –> 2013-2019 2.8885400 7479.17484
SO (27.85,27.95] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 709 1.6286558 2717.6914766 10.1355952 2000-2012 –> 2013-2019 3.0682688 7721.40072
SO (27.85,27.95] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 709 1.5516698 2649.0271053 9.0966674 2000-2012 –> 2013-2019 3.0542668 7772.54860
SO (27.85,27.95] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 709 1.5409988 2639.2416998 8.7822837 2000-2012 –> 2013-2019 2.9762393 7634.43986
SO (27.95,28.05] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 1786 1.6882407 6951.0583598 9.4364874 1982-1999 –> 2000-2012 3.1710443 15980.28933
SO (27.95,28.05] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1786 1.5483054 6643.9876245 7.5549297 1982-1999 –> 2000-2012 2.8861036 15163.34679
SO (27.95,28.05] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 1786 1.7323419 7043.1701928 10.1736604 1982-1999 –> 2000-2012 3.2195838 16087.26729
SO (27.95,28.05] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1786 1.7234592 7026.8075133 9.6018392 1982-1999 –> 2000-2012 3.1951756 16020.70736
SO (27.95,28.05] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1786 1.8015674 7185.1315361 10.7855346 1982-1999 –> 2000-2012 3.4847877 16846.93732
SO (27.95,28.05] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 986 1.9327722 4109.6067463 9.1048682 2000-2012 –> 2013-2019 3.6210129 11060.66511
SO (27.95,28.05] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 986 1.6648161 3817.3040710 7.0125027 2000-2012 –> 2013-2019 3.2131214 10461.29170
SO (27.95,28.05] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 986 1.8748363 4049.5909257 9.4099333 2000-2012 –> 2013-2019 3.6071782 11092.76112
SO (27.95,28.05] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 986 1.8071065 3979.0323733 8.8048453 2000-2012 –> 2013-2019 3.5305658 11005.83989
SO (27.95,28.05] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 986 1.8908856 4068.4000974 8.6762070 2000-2012 –> 2013-2019 3.6771146 11222.98969
SO (28.05,28.1] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 1071 1.2084510 3456.9312497 6.8893057 1982-1999 –> 2000-2012 2.3093792 8056.63029
SO (28.05,28.1] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1071 1.1155347 3287.5593365 8.1328311 1982-1999 –> 2000-2012 2.1524995 7708.01563
SO (28.05,28.1] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 1071 1.1929018 3429.1911457 10.1579725 1982-1999 –> 2000-2012 2.2787248 7987.05738
SO (28.05,28.1] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1071 1.1927699 3430.9543570 10.1729214 1982-1999 –> 2000-2012 2.2759473 7983.43376
SO (28.05,28.1] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1071 1.2196875 3478.7562091 8.1179298 1982-1999 –> 2000-2012 2.3991905 8289.20418
SO (28.05,28.1] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 606 1.4409135 2174.4695378 7.1949964 2000-2012 –> 2013-2019 2.6493644 5631.40079
SO (28.05,28.1] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 606 1.3483022 2095.9550940 8.8409161 2000-2012 –> 2013-2019 2.4638369 5383.51443
SO (28.05,28.1] 4 cstar_tref ~ sal + aou + phosphate + phosphate_star 606 1.4121339 2150.0170644 11.3700248 2000-2012 –> 2013-2019 2.6050356 5579.20821
SO (28.05,28.1] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 606 1.4121202 2152.0053508 11.3570568 2000-2012 –> 2013-2019 2.6048901 5582.95971
SO (28.05,28.1] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 606 1.3968083 2138.7915666 9.0012617 2000-2012 –> 2013-2019 2.6164958 5617.54778
SO (28.1, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 5309 1.0966504 16059.9105945 6.0094365 1982-1999 –> 2000-2012 2.0858171 36797.53418
SO (28.1, Inf] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 5309 1.1018498 16110.1338606 6.1005720 1982-1999 –> 2000-2012 2.0809310 36696.92150
SO (28.1, Inf] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 5309 1.1018770 16108.3952774 6.0859562 1982-1999 –> 2000-2012 2.0912319 36846.81840
SO (28.1, Inf] 4 cstar_tref ~ sal + nitrate + silicate + phosphate_star 5309 1.0966518 16057.9241425 6.0074321 1982-1999 –> 2000-2012 2.0897732 36852.27237
SO (28.1, Inf] 4 cstar_tref ~ sal + silicate + phosphate + phosphate_star 5309 1.1077529 16164.8671348 5.8381245 1982-1999 –> 2000-2012 2.1175228 37203.89999
SO (28.1, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 2711 1.1833450 8620.2523271 6.1959706 2000-2012 –> 2013-2019 2.2799954 24680.16292
SO (28.1, Inf] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2711 1.1918719 8659.1818134 6.0524338 2000-2012 –> 2013-2019 2.2937217 24769.31567
SO (28.1, Inf] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 2711 1.1946469 8669.7908607 6.0803943 2000-2012 –> 2013-2019 2.2965238 24778.18614
SO (28.1, Inf] 4 cstar_tref ~ sal + nitrate + silicate + phosphate_star 2711 1.1852384 8626.9207957 6.1789071 2000-2012 –> 2013-2019 2.2818902 24684.84494
SO (28.1, Inf] 4 cstar_tref ~ sal + silicate + phosphate + phosphate_star 2711 1.1919888 8657.7136946 6.0393123 2000-2012 –> 2013-2019 2.2997417 24822.58083
SO (28.1,28.15] 4 cstar_tref ~ sal + aou + nitrate + silicate 395 0.5288565 629.7012772 1.8285447 1982-1999 –> 2000-2012 1.0213642 1382.24384
SO (28.1,28.15] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 395 0.4829853 560.0239392 1.5318770 1982-1999 –> 2000-2012 0.9517188 1263.01268
SO (28.1,28.15] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 395 0.4692853 537.2913431 1.7413747 1982-1999 –> 2000-2012 0.9103638 1176.91460
SO (28.1,28.15] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 395 0.5103244 601.5215519 1.7907538 1982-1999 –> 2000-2012 0.9951592 1337.70278
SO (28.1,28.15] 4 cstar_tref ~ sal + nitrate + silicate + phosphate_star 395 0.5279427 628.3351075 1.8926109 1982-1999 –> 2000-2012 1.0267214 1394.06143
SO (28.1,28.15] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 209 0.5031920 320.0407883 2.0966657 2000-2012 –> 2013-2019 0.9861773 880.06473
SO (28.1,28.15] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 209 0.5007130 317.9764705 2.3778112 2000-2012 –> 2013-2019 0.9699983 855.26781
SO (28.1,28.15] 4 cstar_tref ~ sal + aou + silicate + phosphate_star 209 0.5573440 360.7649404 2.1406620 2000-2012 –> 2013-2019 1.0676684 962.28649
SO (28.1,28.15] 4 cstar_tref ~ sal + nitrate + silicate + phosphate_star 209 0.5768634 375.1536760 2.1192378 2000-2012 –> 2013-2019 1.1048061 1003.48878
SO (28.1,28.15] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 209 0.5580602 363.3017219 2.6731375 2000-2012 –> 2013-2019 1.0983891 1011.95719
SO (28.15,28.2] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 385 0.3060639 194.9324461 1.0624911 1982-1999 –> 2000-2012 0.6130546 442.64918
SO (28.15,28.2] 4 cstar_tref ~ temp + aou + nitrate + phosphate_star 385 0.3077637 197.1970155 1.1575377 1982-1999 –> 2000-2012 0.6198704 459.14375
SO (28.15,28.2] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 385 0.3073698 198.2108824 1.1858788 1982-1999 –> 2000-2012 0.5773798 319.87950
SO (28.15,28.2] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 385 0.3122313 208.2941302 1.4335968 1982-1999 –> 2000-2012 0.6434714 528.66859
SO (28.15,28.2] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 385 0.3121632 210.1262775 1.4181146 1982-1999 –> 2000-2012 0.6425998 530.11557
SO (28.15,28.2] 4 cstar_tref ~ aou + silicate + phosphate + phosphate_star 183 0.3459452 142.8317039 0.8221319 2000-2012 –> 2013-2019 0.6645772 366.75122
SO (28.15,28.2] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 183 0.3423788 141.0389371 0.8031784 2000-2012 –> 2013-2019 0.6484427 335.97138
SO (28.15,28.2] 5 cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 183 0.3470133 145.9599817 0.8256096 2000-2012 –> 2013-2019 0.6543831 344.17086
SO (28.15,28.2] 4 cstar_tref ~ temp + aou + phosphate + phosphate_star 183 0.3366745 132.8896984 0.8481580 2000-2012 –> 2013-2019 0.6489057 341.18383
SO (28.15,28.2] 5 cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 183 0.3337812 131.7308238 0.8025278 2000-2012 –> 2013-2019 0.6459444 341.85710
SO (28.2, Inf] 3 cstar_tref ~ sal + aou + nitrate 1738 0.3752615 1535.2908035 2.0011297 1982-1999 –> 2000-2012 0.6790183 2669.60844
SO (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 1738 0.3665194 1455.3553954 1.8883521 1982-1999 –> 2000-2012 0.6702525 2591.28860
SO (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + silicate 1738 0.3672295 1462.0838414 1.9673041 1982-1999 –> 2000-2012 0.6709769 2598.24970
SO (28.2, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1738 0.3664835 1457.0151862 1.8977416 1982-1999 –> 2000-2012 0.6701959 2594.61023
SO (28.2, Inf] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1738 0.3918565 1689.7063043 2.1137168 1982-1999 –> 2000-2012 0.7169227 3163.23502
SO (28.2, Inf] 3 cstar_tref ~ sal + aou + nitrate 1005 0.4164558 1101.3567107 2.5953775 2000-2012 –> 2013-2019 0.7917173 2636.64751
SO (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + phosphate_star 1005 0.3988626 1016.5985475 2.4052531 2000-2012 –> 2013-2019 0.7653820 2471.95394
SO (28.2, Inf] 4 cstar_tref ~ sal + aou + nitrate + silicate 1005 0.3952080 998.0966826 2.5161713 2000-2012 –> 2013-2019 0.7624375 2460.18052
SO (28.2, Inf] 5 cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1005 0.3950006 999.0416755 2.5557103 2000-2012 –> 2013-2019 0.7614841 2456.05686
SO (28.2, Inf] 5 cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1005 0.4243307 1143.0096194 2.6684674 2000-2012 –> 2013-2019 0.8161872 2832.71592

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.8582070 1324.47149
Atlantic (-Inf,26] 2000-2012 –> 2013-2019 2.7677442 814.23231
Atlantic (26,26.5] 1982-1999 –> 2000-2012 8.0054093 14332.91416
Atlantic (26,26.5] 2000-2012 –> 2013-2019 7.9249527 8990.12338
Atlantic (26.5,26.75] 1982-1999 –> 2000-2012 6.3958806 16492.17296
Atlantic (26.5,26.75] 2000-2012 –> 2013-2019 6.6898108 10549.90471
Atlantic (26.75,27] 1982-1999 –> 2000-2012 3.8466786 19962.95935
Atlantic (26.75,27] 2000-2012 –> 2013-2019 4.7458231 14674.77818
Atlantic (27,27.25] 1982-1999 –> 2000-2012 3.9503752 16501.55983
Atlantic (27,27.25] 2000-2012 –> 2013-2019 4.5448274 12007.38861
Atlantic (27.25,27.5] 1982-1999 –> 2000-2012 4.2406148 16605.08855
Atlantic (27.25,27.5] 2000-2012 –> 2013-2019 4.8587725 11563.05069
Atlantic (27.5,27.75] 1982-1999 –> 2000-2012 4.8333082 22401.62276
Atlantic (27.5,27.75] 2000-2012 –> 2013-2019 5.3036069 15171.12057
Atlantic (27.75,27.85] 1982-1999 –> 2000-2012 3.3907986 6574.41931
Atlantic (27.75,27.85] 2000-2012 –> 2013-2019 3.4315213 4256.99920
Atlantic (27.85,27.95] 1982-1999 –> 2000-2012 5.7779316 8253.84722
Atlantic (27.85,27.95] 2000-2012 –> 2013-2019 5.8823402 5536.19646
Atlantic (27.95,28.05] 1982-1999 –> 2000-2012 5.8787048 8557.33073
Atlantic (27.95,28.05] 2000-2012 –> 2013-2019 5.1872966 5325.67552
Atlantic (28.05,28.1] 1982-1999 –> 2000-2012 0.3282640 -965.60511
Atlantic (28.05,28.1] 2000-2012 –> 2013-2019 0.4604839 -79.32181
Atlantic (28.1,28.15] 1982-1999 –> 2000-2012 0.2893484 -1360.31481
Atlantic (28.1,28.15] 2000-2012 –> 2013-2019 0.3010462 -757.23131
Atlantic (28.15,28.2] 1982-1999 –> 2000-2012 0.3546729 -1687.48910
Atlantic (28.15,28.2] 2000-2012 –> 2013-2019 0.3613690 -1221.21112
Atlantic (28.2, Inf] 1982-1999 –> 2000-2012 0.4294702 -540.41035
Atlantic (28.2, Inf] 2000-2012 –> 2013-2019 0.4124515 -672.93936
Indo-Pacific (-Inf,26] 1982-1999 –> 2000-2012 14.0300916 84061.77194
Indo-Pacific (-Inf,26] 2000-2012 –> 2013-2019 13.4461854 53949.35152
Indo-Pacific (26,26.5] 1982-1999 –> 2000-2012 9.2817273 76809.24327
Indo-Pacific (26,26.5] 2000-2012 –> 2013-2019 9.4865512 49481.82545
Indo-Pacific (26.5,26.75] 1982-1999 –> 2000-2012 7.9247304 57006.02550
Indo-Pacific (26.5,26.75] 2000-2012 –> 2013-2019 8.3488834 37404.24741
Indo-Pacific (26.75,27] 1982-1999 –> 2000-2012 8.0332192 61382.53803
Indo-Pacific (26.75,27] 2000-2012 –> 2013-2019 9.0815238 40951.54220
Indo-Pacific (27,27.25] 1982-1999 –> 2000-2012 7.1699986 53597.70884
Indo-Pacific (27,27.25] 2000-2012 –> 2013-2019 8.7694808 37862.07546
Indo-Pacific (27.25,27.5] 1982-1999 –> 2000-2012 5.9300259 41239.68701
Indo-Pacific (27.25,27.5] 2000-2012 –> 2013-2019 6.5856507 27923.65083
Indo-Pacific (27.5,27.75] 1982-1999 –> 2000-2012 3.5635625 30434.30879
Indo-Pacific (27.5,27.75] 2000-2012 –> 2013-2019 4.1347141 20587.25128
Indo-Pacific (27.75,27.85] 1982-1999 –> 2000-2012 2.7360823 10050.62535
Indo-Pacific (27.75,27.85] 2000-2012 –> 2013-2019 3.0243300 6936.20481
Indo-Pacific (27.85,27.95] 1982-1999 –> 2000-2012 2.6429484 14646.40233
Indo-Pacific (27.85,27.95] 2000-2012 –> 2013-2019 2.9980181 10036.39329
Indo-Pacific (27.95,28.05] 1982-1999 –> 2000-2012 1.1985292 5218.84976
Indo-Pacific (27.95,28.05] 2000-2012 –> 2013-2019 1.2596792 3470.47440
Indo-Pacific (28.05,28.1] 1982-1999 –> 2000-2012 0.9287923 3323.04733
Indo-Pacific (28.05,28.1] 2000-2012 –> 2013-2019 0.9470628 2135.09974
Indo-Pacific (28.1, Inf] 1982-1999 –> 2000-2012 1.5164998 50854.08984
Indo-Pacific (28.1, Inf] 2000-2012 –> 2013-2019 1.5201517 32974.94083
SO (26,26.5] 1982-1999 –> 2000-2012 4.6885190 1029.98505
SO (26,26.5] 2000-2012 –> 2013-2019 4.4051702 645.88536
SO (26.5,26.75] 1982-1999 –> 2000-2012 5.3854035 3501.55807
SO (26.5,26.75] 2000-2012 –> 2013-2019 4.8992600 2059.53360
SO (26.75,27] 1982-1999 –> 2000-2012 5.1892945 12255.44554
SO (26.75,27] 2000-2012 –> 2013-2019 5.0768032 7576.33708
SO (27,27.25] 1982-1999 –> 2000-2012 3.5693641 23747.89994
SO (27,27.25] 2000-2012 –> 2013-2019 3.8201728 15580.11506
SO (27.25,27.5] 1982-1999 –> 2000-2012 3.1273998 24204.03747
SO (27.25,27.5] 2000-2012 –> 2013-2019 3.5480616 16563.62695
SO (27.5,27.75] 1982-1999 –> 2000-2012 2.7689168 27445.98401
SO (27.5,27.75] 2000-2012 –> 2013-2019 3.2215544 19547.20902
SO (27.75,27.85] 1982-1999 –> 2000-2012 2.6586869 10930.65047
SO (27.75,27.85] 2000-2012 –> 2013-2019 2.9305255 7658.27836
SO (27.85,27.95] 1982-1999 –> 2000-2012 2.8858994 11633.95658
SO (27.85,27.95] 2000-2012 –> 2013-2019 3.0110620 7677.72629
SO (27.95,28.05] 1982-1999 –> 2000-2012 3.1913390 16019.70962
SO (27.95,28.05] 2000-2012 –> 2013-2019 3.5297986 10968.70950
SO (28.05,28.1] 1982-1999 –> 2000-2012 2.2831483 8004.86825
SO (28.05,28.1] 2000-2012 –> 2013-2019 2.5879246 5558.92618
SO (28.1, Inf] 1982-1999 –> 2000-2012 2.0930552 36879.48929
SO (28.1, Inf] 2000-2012 –> 2013-2019 2.2903746 24747.01810
SO (28.1,28.15] 1982-1999 –> 2000-2012 0.9810655 1310.78707
SO (28.1,28.15] 2000-2012 –> 2013-2019 1.0454078 942.61300
SO (28.15,28.2] 1982-1999 –> 2000-2012 0.6192752 456.09132
SO (28.15,28.2] 2000-2012 –> 2013-2019 0.6524506 345.98688
SO (28.2, Inf] 1982-1999 –> 2000-2012 0.6814732 2723.39840
SO (28.2, Inf] 2000-2012 –> 2013-2019 0.7794416 2571.51095

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

6.3 Predictor counts

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

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

# print table
lm_all_stats %>%
  gt(rowname_col = "gamma_slab",
     groupname_col = c("basin", "eras")) %>% 
  summary_rows(
    groups = TRUE,
    fns = list(total = "sum")
  )
aou nitrate phosphate phosphate_star sal silicate temp
Atlantic - 1982-1999 --> 2000-2012
(-Inf,26] 5 1 4 5 3 3 2
(26,26.5] 5 3 2 4 1 4 4
(26.5,26.75] 5 1 4 4 2 5 2
(26.75,27] 5 1 4 4 2 5 2
(27,27.25] 5 1 4 4 3 4 2
(27.25,27.5] 5 5 0 3 1 3 4
(27.5,27.75] 5 1 4 3 5 3 0
(27.75,27.85] 5 1 4 5 3 3 2
(27.85,27.95] 5 1 4 5 1 4 3
(27.95,28.05] 3 4 1 5 0 3 5
(28.05,28.1] 4 2 3 5 1 5 3
(28.1,28.15] 5 2 3 5 2 3 3
(28.15,28.2] 5 5 0 3 4 3 1
(28.2, Inf] 5 5 0 3 4 3 1
total 67.00 33.00 37.00 58.00 32.00 51.00 34.00
Atlantic - 2000-2012 --> 2013-2019
(-Inf,26] 5 1 4 5 3 3 2
(26,26.5] 5 3 2 4 1 4 4
(26.5,26.75] 5 1 4 4 2 5 2
(26.75,27] 5 1 4 4 3 5 1
(27,27.25] 4 1 3 4 5 5 0
(27.25,27.5] 5 5 0 3 4 3 1
(27.5,27.75] 5 1 4 3 5 3 0
(27.75,27.85] 5 2 3 5 4 3 1
(27.85,27.95] 5 1 4 5 1 4 3
(27.95,28.05] 3 4 1 5 0 3 5
(28.05,28.1] 4 2 3 5 2 4 3
(28.1,28.15] 5 1 4 5 2 3 3
(28.15,28.2] 5 5 0 3 4 3 1
(28.2, Inf] 5 5 0 3 4 3 1
total 66.00 33.00 36.00 58.00 40.00 51.00 27.00
Indo-Pacific - 1982-1999 --> 2000-2012
(-Inf,26] 4 1 4 5 2 4 3
(26,26.5] 5 1 4 4 1 5 3
(26.5,26.75] 4 2 3 5 1 5 3
(26.75,27] 5 0 5 5 2 3 2
(27,27.25] 5 0 5 5 2 3 2
(27.25,27.5] 4 2 3 5 3 5 2
(27.5,27.75] 4 2 2 5 4 5 1
(27.75,27.85] 5 1 4 5 1 4 3
(27.85,27.95] 5 3 2 4 2 5 3
(27.95,28.05] 4 2 3 5 1 5 3
(28.05,28.1] 4 2 3 4 4 5 1
(28.1, Inf] 4 3 2 5 2 5 3
total 53.00 19.00 40.00 57.00 25.00 54.00 29.00
Indo-Pacific - 2000-2012 --> 2013-2019
(-Inf,26] 4 1 4 5 2 4 3
(26,26.5] 5 1 4 4 1 5 3
(26.5,26.75] 4 2 3 5 1 5 3
(26.75,27] 5 0 5 5 2 3 2
(27,27.25] 5 1 4 5 1 4 3
(27.25,27.5] 4 2 3 5 3 5 2
(27.5,27.75] 4 2 2 5 4 5 1
(27.75,27.85] 5 1 4 5 1 4 3
(27.85,27.95] 5 3 2 4 2 5 3
(27.95,28.05] 5 1 4 4 2 5 2
(28.05,28.1] 4 2 3 4 4 5 1
(28.1, Inf] 4 3 2 5 2 5 3
total 54.00 19.00 40.00 56.00 25.00 55.00 29.00
SO - 1982-1999 --> 2000-2012
(26,26.5] 5 2 3 5 1 3 4
(26.5,26.75] 5 5 0 5 2 3 2
(26.75,27] 5 2 2 5 0 3 5
(27,27.25] 5 2 2 5 0 3 5
(27.25,27.5] 3 2 2 5 0 5 5
(27.5,27.75] 5 4 1 4 2 5 2
(27.75,27.85] 5 4 1 4 2 5 2
(27.85,27.95] 5 1 4 5 3 3 2
(27.95,28.05] 5 2 3 5 4 3 1
(28.05,28.1] 5 2 3 5 4 3 1
(28.1, Inf] 3 2 2 5 5 5 0
(28.1,28.15] 4 3 1 4 5 5 0
(28.15,28.2] 5 2 3 5 1 3 4
(28.2, Inf] 5 4 1 3 5 3 0
total 65.00 37.00 28.00 65.00 34.00 52.00 33.00
SO - 2000-2012 --> 2013-2019
(26,26.5] 5 2 3 5 2 3 3
(26.5,26.75] 5 5 0 5 2 3 2
(26.75,27] 5 2 2 5 0 3 5
(27,27.25] 5 2 2 5 0 3 5
(27.25,27.5] 4 1 3 4 0 5 5
(27.5,27.75] 5 4 1 4 3 5 1
(27.75,27.85] 5 4 1 4 3 5 1
(27.85,27.95] 5 3 2 5 2 5 2
(27.95,28.05] 5 3 2 5 4 3 1
(28.05,28.1] 5 2 3 5 4 3 1
(28.1, Inf] 3 2 2 5 5 5 0
(28.1,28.15] 4 3 1 5 4 5 1
(28.15,28.2] 5 1 4 5 1 4 3
(28.2, Inf] 5 4 1 3 5 3 0
total 66.00 38.00 27.00 65.00 35.00 55.00 30.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
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
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