Last updated: 2020-12-01

Checks: 7 0

Knit directory: emlr_obs_v_XXX/

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

R Markdown file: up-to-date

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

Environment: empty

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

Seed: set.seed(20200707)

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

Session information: recorded

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

Cache: none

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

File paths: relative

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

Repository version: 6287a18

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 6287a18. 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:
    Modified:   analysis/_site.yml
    Deleted:    code/README.md
    Modified:   code/Workflowr_project_managment.R
    Deleted:    code/biogeochemical_functions.R
    Deleted:    code/eda.R
    Deleted:    code/mapping_functions.R
    Deleted:    code/plotting_functions.R
    Deleted:    code/python_scripts/Gamma_GLODAP_python.py
    Deleted:    code/test_scripts/analysis_previous_studies.Rmd
    Deleted:    code/test_scripts/analysis_this_study_3d.Rmd
    Deleted:    code/test_scripts/model_nested_df.R
    Deleted:    code/test_scripts/raster_interpolation.R
    Deleted:    code/test_scripts/read_WOCE.Rmd
    Deleted:    code/test_scripts/read_World_Ocean_Atlas_2013_Clement.Rmd
    Deleted:    code/test_scripts/read_dclement_nc.Rmd
    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 91435ae jens-daniel-mueller 2020-12-01 Build site.
Rmd 17d09be jens-daniel-mueller 2020-12-01 auto eras naming
html cf19652 jens-daniel-mueller 2020-11-30 Build site.
Rmd 0895ad6 jens-daniel-mueller 2020-11-30 rebuild with all plot output
Rmd 2842970 jens-daniel-mueller 2020-11-30 cleaned for eMLR part only
html 196be51 jens-daniel-mueller 2020-11-30 Build site.
Rmd 7a4b015 jens-daniel-mueller 2020-11-30 first rebuild on ETH server
Rmd bc61ce3 Jens Müller 2020-11-30 Initial commit 
path_functions      <- "/nfs/kryo/work/updata/emlr_cant/utilities/functions/"
path_files          <- "/nfs/kryo/work/updata/emlr_cant/utilities/files/"
path_preprocessing    <-
  "/nfs/kryo/work/updata/emlr_cant/observations/preprocessing/"

path_version_data     <-
  paste(
    "/nfs/kryo/work/updata/emlr_cant/observations/",
    params_local$Version_ID,
    "/data/",
    sep = ""
  )

path_version_figures  <-
  paste(
    "/nfs/kryo/work/updata/emlr_cant/observations/",
    params_local$Version_ID,
    "/figures/",
    sep = ""
  )

1 Libraries

Loading libraries specific to the the analysis performed in this section.

library(lubridate)
library(GGally)
library(olsrr)
library(knitr)
library(kableExtra)
library(broom)
library(corrr)
library(GGally)

2 Required data

Required are:

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

3 Fitting all models

3.1 RMSE calculation

Individual linear regression models were fitted for C* as a function of each predictor combination. Fitting was performed separately within each basin, era, and slab. The root mean squared error (RMSE) was calculated for each fitted model.

# loops across all basins and eras
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  )
    
    # subset one basin and era for fitting
    GLODAP_basin_era <- GLODAP %>%
      filter(basin == i_basin, era == i_era)
    
    # loop across all gamma slabs in the current basin and 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)
      
      # 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(
        params_local$MLR_full_model,
        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)
      
      # extract diagnostics of each linear model
      lm_all <- as_tibble(lm_all$result)
      
      lm_all <- lm_all %>% 
        mutate(basin = i_basin,
               era = i_era,
               gamma_slab = i_gamma_slab)
      
      # expanding table with model diagnostics for all eras, basins, slabs
      if (exists("lm_all_stats")) {
        lm_all_stats <- bind_rows(lm_all_stats, lm_all)
      }
      
      if (!exists("lm_all_stats")) {
        lm_all_stats <- lm_all
      }
      
      rm(lm_full, lm_all)
      
    }
    
  }
  
}

# write model diagnostics to file

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

rm(i_gamma_slab, i_era, i_basin,
   GLODAP_basin_era, GLODAP_basin_era_slab)

RMSE was plotted against the number of predictors (limited to 2 - 5).

lm_all_stats <- read_csv(paste(path_version_data,
                               "lm_all_stats.csv",
                               sep = ""))


lm_all_stats <- lm_all_stats %>% 
  filter(n >= params_local$predictors_min,
         n <= params_local$predictors_max)

lm_all_stats %>% 
  ggplot(aes(n, rmse, col = basin)) +
  geom_hline(yintercept = 10) + 
  geom_point(shape = 21) +
  facet_grid(gamma_slab~era) +
  scale_color_brewer(palette = "Set1")

Version Author Date
91435ae jens-daniel-mueller 2020-12-01
196be51 jens-daniel-mueller 2020-11-30

3.2 RMSE alternatives

AIC is an alternative criterion to RMSE to judge model quality, but not (yet) taken into account.

lm_all_stats %>% 
  ggplot(aes(rmse, aic, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  facet_grid(era~basin)

Version Author Date
91435ae jens-daniel-mueller 2020-12-01
196be51 jens-daniel-mueller 2020-11-30

3.3 Predictor selection

Within each basin and slab, the 10 linear regression models with lowest summed RMSE across all eras were selected.

# select relevant columns
lm_all_stats <- lm_all_stats %>%
  select(basin, era, gamma_slab, predictors, rmse)

# pivot table
lm_all_stats_wide <- lm_all_stats %>%
  pivot_wider(names_from = era,
              values_from = rmse,
              names_prefix = "rmse_")

# calculate coefficient difference
# it would be good to recode this for a variable number of eras
# currently it is hard coded by era names
lm_all_stats_wide <- lm_all_stats_wide %>%
  mutate(
    rmse_sum_JGOFS_GO = rmse_JGOFS_WOCE + rmse_GO_SHIP,
    rmse_sum_GO_new = rmse_GO_SHIP + rmse_new_era
  )

# pivot back to long format
lm_all_stats <- lm_all_stats_wide %>%
  pivot_longer(
    c(rmse_sum_JGOFS_GO, rmse_sum_GO_new),
    names_to = "eras",
    values_to = "rmse_sum",
    names_prefix = "rmse_sum_"
  )

# subset models with lowest RMSE
lm_best <- lm_all_stats %>%
  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, predictors)

# plot table for website
kable(lm_best) %>%
  add_header_above() %>%
  kable_styling() %>%
  scroll_box(width = "100%", height = "400px")
# select relevant columns
lm_all_stats <- lm_all_stats %>%
  select(basin, era, gamma_slab, predictors, rmse)

# select relevant columns
lm_all_stats <- lm_all_stats %>%
  arrange(era) %>% 
  group_by(basin, gamma_slab, predictors) %>% 
  mutate(rmse_sum = rmse + lag(rmse),
         eras = paste(lag(era), era, sep = " --> ")) %>% 
  ungroup() %>% 
  arrange(basin, gamma_slab, predictors, era) %>% 
  drop_na()

# subset models with lowest RMSE
lm_best <- lm_all_stats %>%
  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, predictors, era)

# plot table for website
kable(lm_best) %>%
  add_header_above() %>%
  kable_styling() %>%
  scroll_box(width = "100%", height = "400px")
basin era gamma_slab predictors rmse rmse_sum eras
Atlantic (1999,2012] (-Inf,26] sal aou oxygen silicate phosphate 6.883345 12.955019 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (-Inf,26] sal aou oxygen silicate phosphate_star 6.883345 12.955019 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (-Inf,26] sal aou silicate phosphate phosphate_star 6.883345 12.955019 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (-Inf,26] sal tem aou oxygen phosphate 7.039840 13.098128 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (-Inf,26] sal tem aou oxygen phosphate_star 7.039840 13.098128 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (-Inf,26] sal tem aou silicate phosphate 6.955823 12.993571 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (-Inf,26] sal tem aou silicate phosphate_star 6.964611 12.995587 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (-Inf,26] sal tem oxygen silicate phosphate 6.970726 12.991033 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (-Inf,26] sal tem oxygen silicate phosphate_star 6.970726 12.991033 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (-Inf,26] sal tem silicate phosphate phosphate_star 6.970726 12.991033 (1981,1999] –> (1999,2012]
Atlantic (2012,Inf] (-Inf,26] sal aou oxygen silicate phosphate 2.940339 9.823684 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (-Inf,26] sal aou oxygen silicate phosphate_star 2.940339 9.823684 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (-Inf,26] sal aou silicate phosphate phosphate_star 2.940339 9.823684 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (-Inf,26] sal tem aou oxygen phosphate 2.901429 9.941270 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (-Inf,26] sal tem aou oxygen phosphate_star 2.901429 9.941270 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (-Inf,26] sal tem aou silicate phosphate 2.926536 9.882359 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (-Inf,26] sal tem aou silicate phosphate_star 2.924241 9.888852 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (-Inf,26] sal tem oxygen silicate phosphate 2.922521 9.893247 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (-Inf,26] sal tem oxygen silicate phosphate_star 2.922521 9.893247 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (-Inf,26] sal tem silicate phosphate phosphate_star 2.922521 9.893247 (1999,2012] –> (2012,Inf]
Atlantic (1999,2012] (26,26.5] sal tem aou oxygen phosphate 5.005534 9.832733 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (26,26.5] sal tem aou oxygen phosphate_star 5.005534 9.832733 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (26,26.5] sal tem aou phosphate phosphate_star 5.005534 9.832733 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (26,26.5] tem aou oxygen phosphate 5.005534 9.841368 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (26,26.5] tem aou oxygen phosphate phosphate_star 5.005534 9.841368 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (26,26.5] tem aou oxygen phosphate_star 5.005534 9.841368 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (26,26.5] tem aou oxygen silicate phosphate 5.004072 9.807359 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (26,26.5] tem aou oxygen silicate phosphate_star 5.004072 9.807359 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (26,26.5] tem aou phosphate phosphate_star 5.005534 9.841368 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (26,26.5] tem aou silicate phosphate phosphate_star 5.004072 9.807359 (1981,1999] –> (1999,2012]
Atlantic (2012,Inf] (26,26.5] sal aou oxygen phosphate 2.906765 8.118811 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (26,26.5] sal aou oxygen phosphate phosphate_star 2.906765 8.118811 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (26,26.5] sal aou oxygen phosphate_star 2.906765 8.118811 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (26,26.5] sal aou oxygen silicate phosphate 2.864980 8.076720 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (26,26.5] sal aou oxygen silicate phosphate_star 2.864980 8.076720 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (26,26.5] sal aou phosphate phosphate_star 2.906765 8.118811 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (26,26.5] sal aou silicate phosphate phosphate_star 2.864980 8.076720 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (26,26.5] sal tem aou oxygen phosphate 2.887114 7.892648 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (26,26.5] sal tem aou oxygen phosphate_star 2.887114 7.892648 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (26,26.5] sal tem aou phosphate phosphate_star 2.887114 7.892648 (1999,2012] –> (2012,Inf]
Atlantic (1999,2012] (26.5,26.75] sal tem aou oxygen phosphate 3.929664 7.689631 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (26.5,26.75] sal tem aou oxygen phosphate_star 3.929664 7.689631 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (26.5,26.75] sal tem aou phosphate phosphate_star 3.929664 7.689631 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (26.5,26.75] tem aou oxygen phosphate 3.962835 7.773281 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (26.5,26.75] tem aou oxygen phosphate phosphate_star 3.962835 7.773281 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (26.5,26.75] tem aou oxygen phosphate_star 3.962835 7.773281 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (26.5,26.75] tem aou oxygen silicate phosphate 3.953607 7.746285 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (26.5,26.75] tem aou oxygen silicate phosphate_star 3.953607 7.746285 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (26.5,26.75] tem aou phosphate phosphate_star 3.962835 7.773281 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (26.5,26.75] tem aou silicate phosphate phosphate_star 3.953607 7.746285 (1981,1999] –> (1999,2012]
Atlantic (2012,Inf] (26.5,26.75] sal aou oxygen silicate phosphate 3.056592 7.079905 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (26.5,26.75] sal aou oxygen silicate phosphate_star 3.056592 7.079905 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (26.5,26.75] sal tem aou oxygen phosphate 3.143854 7.073518 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (26.5,26.75] sal tem aou oxygen phosphate_star 3.143854 7.073518 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (26.5,26.75] sal tem aou phosphate phosphate_star 3.143854 7.073518 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (26.5,26.75] sal tem aou silicate phosphate 3.074346 7.021933 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (26.5,26.75] sal tem aou silicate phosphate_star 3.088551 7.018308 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (26.5,26.75] sal tem oxygen silicate phosphate 3.090373 7.018476 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (26.5,26.75] sal tem oxygen silicate phosphate_star 3.090373 7.018476 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (26.5,26.75] sal tem silicate phosphate phosphate_star 3.090373 7.018476 (1999,2012] –> (2012,Inf]
Atlantic (1999,2012] (26.75,27] sal tem aou oxygen phosphate 4.551459 8.904452 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (26.75,27] sal tem aou oxygen phosphate_star 4.551459 8.904452 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (26.75,27] sal tem aou phosphate phosphate_star 4.551459 8.904452 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (26.75,27] sal tem aou phosphate_star 4.553368 8.906802 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (26.75,27] sal tem aou silicate phosphate 4.552190 8.888517 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (26.75,27] sal tem aou silicate phosphate_star 4.535845 8.862458 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (26.75,27] sal tem oxygen phosphate 4.553586 8.907357 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (26.75,27] sal tem oxygen silicate phosphate 4.535899 8.862723 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (26.75,27] sal tem oxygen silicate phosphate_star 4.535899 8.862723 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (26.75,27] sal tem silicate phosphate phosphate_star 4.535899 8.862723 (1981,1999] –> (1999,2012]
Atlantic (2012,Inf] (26.75,27] sal tem aou silicate phosphate 3.393347 7.945537 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (26.75,27] sal tem aou silicate phosphate_star 3.384849 7.920693 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (26.75,27] sal tem oxygen silicate phosphate 3.384667 7.920566 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (26.75,27] sal tem oxygen silicate phosphate_star 3.384667 7.920566 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (26.75,27] sal tem silicate phosphate phosphate_star 3.384667 7.920566 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (26.75,27] tem aou oxygen silicate phosphate 3.432645 8.002118 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (26.75,27] tem aou oxygen silicate phosphate_star 3.432645 8.002118 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (26.75,27] tem aou silicate phosphate phosphate_star 3.432645 8.002118 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (26.75,27] tem aou silicate phosphate_star 3.433313 8.003036 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (26.75,27] tem oxygen silicate phosphate 3.433597 8.003349 (1999,2012] –> (2012,Inf]
Atlantic (1999,2012] (27,27.25] sal tem aou silicate phosphate_star 4.873592 10.072126 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27,27.25] sal tem oxygen silicate phosphate 4.881072 10.078736 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27,27.25] sal tem oxygen silicate phosphate_star 4.881072 10.078736 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27,27.25] sal tem silicate phosphate phosphate_star 4.881072 10.078736 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27,27.25] tem aou oxygen silicate phosphate 4.841409 10.089543 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27,27.25] tem aou oxygen silicate phosphate_star 4.841409 10.089543 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27,27.25] tem aou silicate phosphate phosphate_star 4.841409 10.089543 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27,27.25] tem aou silicate phosphate_star 4.900926 10.152563 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27,27.25] tem oxygen silicate phosphate_star 4.908671 10.160012 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27,27.25] tem silicate phosphate phosphate_star 4.908671 10.160012 (1981,1999] –> (1999,2012]
Atlantic (2012,Inf] (27,27.25] sal tem aou silicate phosphate 4.515975 9.470028 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27,27.25] sal tem aou silicate phosphate_star 4.502162 9.375754 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27,27.25] sal tem oxygen silicate phosphate 4.502181 9.383254 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27,27.25] sal tem oxygen silicate phosphate_star 4.502181 9.383254 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27,27.25] sal tem silicate phosphate phosphate_star 4.502181 9.383254 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27,27.25] tem aou oxygen silicate phosphate 4.682454 9.523863 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27,27.25] tem aou oxygen silicate phosphate_star 4.682454 9.523863 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27,27.25] tem aou silicate phosphate phosphate_star 4.682454 9.523863 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27,27.25] tem aou silicate phosphate_star 4.692716 9.593642 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27,27.25] tem oxygen silicate phosphate_star 4.692199 9.600870 (1999,2012] –> (2012,Inf]
Atlantic (1999,2012] (27.25,27.5] sal tem aou oxygen phosphate 4.211690 8.795424 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.25,27.5] sal tem aou oxygen phosphate_star 4.211690 8.795424 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.25,27.5] sal tem aou phosphate phosphate_star 4.211690 8.795424 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.25,27.5] sal tem aou phosphate_star 4.213212 8.806263 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.25,27.5] sal tem aou silicate phosphate_star 4.208525 8.798014 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.25,27.5] sal tem oxygen phosphate 4.218653 8.808569 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.25,27.5] sal tem oxygen phosphate phosphate_star 4.218653 8.808569 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.25,27.5] sal tem oxygen silicate phosphate 4.215196 8.802014 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.25,27.5] sal tem oxygen silicate phosphate_star 4.215196 8.802014 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.25,27.5] sal tem silicate phosphate phosphate_star 4.215196 8.802014 (1981,1999] –> (1999,2012]
Atlantic (2012,Inf] (27.25,27.5] sal aou oxygen silicate phosphate 4.094948 8.409799 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.25,27.5] sal aou oxygen silicate phosphate_star 4.094948 8.409799 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.25,27.5] sal aou silicate phosphate phosphate_star 4.094948 8.409799 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.25,27.5] sal tem aou oxygen phosphate_star 4.222095 8.433785 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.25,27.5] sal tem aou phosphate phosphate_star 4.222095 8.433785 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.25,27.5] sal tem aou silicate phosphate 4.074540 8.315126 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.25,27.5] sal tem aou silicate phosphate_star 4.080304 8.288829 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.25,27.5] sal tem oxygen silicate phosphate 4.077325 8.292522 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.25,27.5] sal tem oxygen silicate phosphate_star 4.077325 8.292522 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.25,27.5] sal tem silicate phosphate phosphate_star 4.077325 8.292522 (1999,2012] –> (2012,Inf]
Atlantic (1999,2012] (27.5,27.75] sal aou oxygen silicate phosphate 4.423711 8.743900 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.5,27.75] sal aou oxygen silicate phosphate_star 4.423711 8.743900 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.5,27.75] sal aou silicate phosphate phosphate_star 4.423711 8.743900 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.5,27.75] sal tem aou silicate phosphate 4.405532 8.736716 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.5,27.75] sal tem aou silicate phosphate_star 4.414472 8.779972 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.5,27.75] sal tem oxygen silicate phosphate 4.404547 8.754463 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.5,27.75] sal tem oxygen silicate phosphate_star 4.404547 8.754463 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.5,27.75] sal tem silicate phosphate phosphate_star 4.404547 8.754463 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.5,27.75] tem aou oxygen silicate phosphate 4.472363 8.793834 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.5,27.75] tem aou oxygen silicate phosphate_star 4.472363 8.793834 (1981,1999] –> (1999,2012]
Atlantic (2012,Inf] (27.5,27.75] sal aou oxygen silicate phosphate 4.271770 8.695481 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.5,27.75] sal aou oxygen silicate phosphate_star 4.271770 8.695481 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.5,27.75] sal aou silicate phosphate phosphate_star 4.271770 8.695481 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.5,27.75] sal tem aou silicate phosphate 4.267580 8.673112 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.5,27.75] sal tem aou silicate phosphate_star 4.263957 8.678429 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.5,27.75] sal tem oxygen silicate phosphate 4.263383 8.667930 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.5,27.75] sal tem oxygen silicate phosphate_star 4.263383 8.667930 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.5,27.75] sal tem silicate phosphate phosphate_star 4.263383 8.667930 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.5,27.75] tem aou oxygen silicate phosphate 4.268160 8.740523 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.5,27.75] tem aou silicate phosphate phosphate_star 4.268160 8.740523 (1999,2012] –> (2012,Inf]
Atlantic (1999,2012] (27.75,27.85] sal aou oxygen silicate phosphate 4.805949 9.103865 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.75,27.85] sal aou oxygen silicate phosphate_star 4.805949 9.103865 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.75,27.85] sal aou silicate phosphate 4.806319 9.104645 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.75,27.85] sal aou silicate phosphate phosphate_star 4.805949 9.103865 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.75,27.85] sal oxygen silicate phosphate 4.815138 9.132461 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.75,27.85] sal tem aou silicate phosphate 4.805899 9.101797 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.75,27.85] sal tem aou silicate phosphate_star 4.811399 9.130734 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.75,27.85] sal tem oxygen silicate phosphate 4.806687 9.113266 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.75,27.85] sal tem oxygen silicate phosphate_star 4.806687 9.113266 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.75,27.85] sal tem silicate phosphate phosphate_star 4.806687 9.113266 (1981,1999] –> (1999,2012]
Atlantic (2012,Inf] (27.75,27.85] sal aou oxygen silicate phosphate 5.680123 10.486072 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.75,27.85] sal aou oxygen silicate phosphate_star 5.680123 10.486072 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.75,27.85] sal aou silicate phosphate phosphate_star 5.680123 10.486072 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.75,27.85] sal tem aou silicate phosphate_star 5.673680 10.485079 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.75,27.85] sal tem oxygen silicate phosphate 5.702448 10.509135 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.75,27.85] sal tem oxygen silicate phosphate_star 5.702448 10.509135 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.75,27.85] sal tem silicate phosphate phosphate_star 5.702448 10.509135 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.75,27.85] tem aou oxygen silicate phosphate 5.615715 10.439808 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.75,27.85] tem aou oxygen silicate phosphate_star 5.615715 10.439808 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.75,27.85] tem aou silicate phosphate phosphate_star 5.615715 10.439808 (1999,2012] –> (2012,Inf]
Atlantic (1999,2012] (27.85,27.95] sal aou oxygen silicate phosphate 4.978568 8.684765 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.85,27.95] sal aou oxygen silicate phosphate_star 4.978568 8.684765 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.85,27.95] sal aou silicate phosphate phosphate_star 4.978568 8.684765 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.85,27.95] sal aou silicate phosphate_star 5.015519 8.726941 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.85,27.95] sal tem aou silicate phosphate 4.979564 8.696440 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.85,27.95] sal tem aou silicate phosphate_star 4.977605 8.685033 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.85,27.95] sal tem oxygen silicate phosphate 4.978684 8.692249 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.85,27.95] sal tem oxygen silicate phosphate_star 4.978684 8.692249 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.85,27.95] sal tem silicate phosphate 5.090386 8.910207 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.85,27.95] sal tem silicate phosphate phosphate_star 4.978684 8.692249 (1981,1999] –> (1999,2012]
Atlantic (2012,Inf] (27.85,27.95] sal aou oxygen silicate phosphate 4.616917 9.595485 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.85,27.95] sal aou oxygen silicate phosphate_star 4.616917 9.595485 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.85,27.95] sal aou silicate phosphate phosphate_star 4.616917 9.595485 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.85,27.95] sal aou silicate phosphate_star 4.794481 9.810001 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.85,27.95] sal tem aou silicate phosphate 4.658115 9.637678 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.85,27.95] sal tem aou silicate phosphate_star 4.636104 9.613709 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.85,27.95] sal tem oxygen silicate phosphate 4.659656 9.638341 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.85,27.95] sal tem oxygen silicate phosphate_star 4.659656 9.638341 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.85,27.95] sal tem silicate phosphate 4.665649 9.756035 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.85,27.95] sal tem silicate phosphate phosphate_star 4.659656 9.638341 (1999,2012] –> (2012,Inf]
Atlantic (1999,2012] (27.95,28.05] sal aou oxygen silicate phosphate 5.898927 9.446795 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.95,28.05] sal aou oxygen silicate phosphate_star 5.898927 9.446795 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.95,28.05] sal aou silicate phosphate phosphate_star 5.898927 9.446795 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.95,28.05] sal aou silicate phosphate_star 5.899193 9.452696 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.95,28.05] sal oxygen silicate phosphate_star 5.946224 9.594932 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.95,28.05] sal tem aou silicate phosphate 5.927518 9.504339 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.95,28.05] sal tem aou silicate phosphate_star 5.898986 9.449695 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.95,28.05] sal tem oxygen silicate phosphate 5.914457 9.480325 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.95,28.05] sal tem oxygen silicate phosphate_star 5.914457 9.480325 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (27.95,28.05] sal tem silicate phosphate phosphate_star 5.914457 9.480325 (1981,1999] –> (1999,2012]
Atlantic (2012,Inf] (27.95,28.05] sal aou oxygen silicate phosphate 4.630927 10.529854 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.95,28.05] sal aou oxygen silicate phosphate_star 4.630927 10.529854 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.95,28.05] sal aou silicate phosphate phosphate_star 4.630927 10.529854 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.95,28.05] sal aou silicate phosphate_star 4.680224 10.579417 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.95,28.05] sal oxygen silicate phosphate_star 4.809349 10.755573 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.95,28.05] sal tem aou silicate phosphate 4.727004 10.654522 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.95,28.05] sal tem aou silicate phosphate_star 4.659611 10.558597 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.95,28.05] sal tem oxygen silicate phosphate 4.705040 10.619497 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.95,28.05] sal tem oxygen silicate phosphate_star 4.705040 10.619497 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (27.95,28.05] sal tem silicate phosphate phosphate_star 4.705040 10.619497 (1999,2012] –> (2012,Inf]
Atlantic (1999,2012] (28.05,28.1] sal tem aou oxygen phosphate 6.223403 9.743365 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.05,28.1] sal tem aou oxygen phosphate_star 6.223403 9.743365 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.05,28.1] sal tem aou phosphate phosphate_star 6.223403 9.743365 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.05,28.1] sal tem aou silicate phosphate_star 6.388721 9.974259 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.05,28.1] tem aou oxygen phosphate 6.304099 9.974618 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.05,28.1] tem aou oxygen phosphate phosphate_star 6.304099 9.974618 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.05,28.1] tem aou oxygen phosphate_star 6.304099 9.974618 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.05,28.1] tem aou oxygen silicate phosphate 6.228880 9.874605 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.05,28.1] tem aou oxygen silicate phosphate_star 6.228880 9.874605 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.05,28.1] tem aou silicate phosphate phosphate_star 6.228880 9.874605 (1981,1999] –> (1999,2012]
Atlantic (2012,Inf] (28.05,28.1] sal tem aou oxygen phosphate 4.139601 10.363004 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.05,28.1] sal tem aou oxygen phosphate_star 4.139601 10.363004 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.05,28.1] sal tem aou phosphate phosphate_star 4.139601 10.363004 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.05,28.1] tem aou oxygen phosphate 4.204511 10.508610 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.05,28.1] tem aou oxygen phosphate phosphate_star 4.204511 10.508610 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.05,28.1] tem aou oxygen phosphate_star 4.204511 10.508610 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.05,28.1] tem aou oxygen silicate phosphate 4.035435 10.264315 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.05,28.1] tem aou oxygen silicate phosphate_star 4.035435 10.264315 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.05,28.1] tem aou phosphate phosphate_star 4.204511 10.508610 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.05,28.1] tem aou silicate phosphate phosphate_star 4.035435 10.264315 (1999,2012] –> (2012,Inf]
Atlantic (1999,2012] (28.1,28.15] sal tem aou oxygen phosphate 6.851585 10.313701 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.1,28.15] sal tem aou oxygen phosphate_star 6.851585 10.313701 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.1,28.15] sal tem aou phosphate phosphate_star 6.851585 10.313701 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.1,28.15] sal tem aou phosphate_star 6.860555 10.331526 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.1,28.15] sal tem aou silicate phosphate_star 6.854317 10.295464 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.1,28.15] sal tem oxygen phosphate 6.860056 10.333059 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.1,28.15] sal tem oxygen silicate phosphate 6.855432 10.295980 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.1,28.15] sal tem oxygen silicate phosphate_star 6.855432 10.295980 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.1,28.15] sal tem silicate phosphate phosphate_star 6.855432 10.295980 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.1,28.15] sal tem silicate phosphate_star 6.856695 10.304339 (1981,1999] –> (1999,2012]
Atlantic (2012,Inf] (28.1,28.15] sal tem aou oxygen phosphate 4.136313 10.987898 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.1,28.15] sal tem aou oxygen phosphate_star 4.136313 10.987898 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.1,28.15] sal tem aou phosphate phosphate_star 4.136313 10.987898 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.1,28.15] sal tem aou silicate phosphate_star 4.165873 11.020190 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.1,28.15] sal tem oxygen silicate phosphate 4.170693 11.026125 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.1,28.15] sal tem oxygen silicate phosphate_star 4.170693 11.026125 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.1,28.15] sal tem silicate phosphate phosphate_star 4.170693 11.026125 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.1,28.15] tem aou oxygen silicate phosphate 4.136659 11.030173 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.1,28.15] tem aou oxygen silicate phosphate_star 4.136659 11.030173 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.1,28.15] tem aou silicate phosphate phosphate_star 4.136659 11.030173 (1999,2012] –> (2012,Inf]
Atlantic (1999,2012] (28.15,28.2] sal aou oxygen silicate phosphate 3.557614 7.040826 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.15,28.2] sal aou oxygen silicate phosphate_star 3.557614 7.040826 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.15,28.2] sal aou silicate phosphate phosphate_star 3.557614 7.040826 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.15,28.2] sal tem aou oxygen phosphate 3.555102 7.021909 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.15,28.2] sal tem aou oxygen phosphate_star 3.555102 7.021909 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.15,28.2] sal tem aou phosphate phosphate_star 3.555102 7.021909 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.15,28.2] tem aou oxygen phosphate 3.558133 7.041274 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.15,28.2] tem aou oxygen silicate phosphate 3.553318 7.036315 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.15,28.2] tem aou oxygen silicate phosphate_star 3.553318 7.036315 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.15,28.2] tem aou silicate phosphate phosphate_star 3.553318 7.036315 (1981,1999] –> (1999,2012]
Atlantic (2012,Inf] (28.15,28.2] sal aou oxygen silicate phosphate 2.832647 6.390261 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.15,28.2] sal aou oxygen silicate phosphate_star 2.832647 6.390261 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.15,28.2] sal tem aou oxygen phosphate 2.764973 6.320075 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.15,28.2] sal tem aou oxygen phosphate_star 2.764973 6.320075 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.15,28.2] sal tem aou phosphate 2.765279 6.354653 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.15,28.2] sal tem aou phosphate phosphate_star 2.764973 6.320075 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.15,28.2] sal tem aou silicate phosphate 2.765274 6.338808 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.15,28.2] sal tem oxygen silicate phosphate 2.766409 6.378395 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.15,28.2] sal tem oxygen silicate phosphate_star 2.766409 6.378395 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.15,28.2] sal tem silicate phosphate phosphate_star 2.766409 6.378395 (1999,2012] –> (2012,Inf]
Atlantic (1999,2012] (28.2, Inf] sal aou oxygen phosphate 2.978002 5.827047 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.2, Inf] sal aou oxygen phosphate phosphate_star 2.978002 5.827047 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.2, Inf] sal aou oxygen phosphate_star 2.978002 5.827047 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.2, Inf] sal aou oxygen silicate phosphate 2.977938 5.819401 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.2, Inf] sal aou oxygen silicate phosphate_star 2.977938 5.819401 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.2, Inf] sal aou phosphate phosphate_star 2.978002 5.827047 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.2, Inf] sal aou silicate phosphate phosphate_star 2.977938 5.819401 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.2, Inf] sal tem aou oxygen phosphate 2.977333 5.825104 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.2, Inf] sal tem aou oxygen phosphate_star 2.977333 5.825104 (1981,1999] –> (1999,2012]
Atlantic (1999,2012] (28.2, Inf] sal tem aou phosphate phosphate_star 2.977333 5.825104 (1981,1999] –> (1999,2012]
Atlantic (2012,Inf] (28.2, Inf] sal aou oxygen phosphate 3.043592 6.021594 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.2, Inf] sal aou oxygen phosphate phosphate_star 3.043592 6.021594 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.2, Inf] sal aou oxygen phosphate_star 3.043592 6.021594 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.2, Inf] sal aou oxygen silicate phosphate 3.042204 6.020142 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.2, Inf] sal aou oxygen silicate phosphate_star 3.042204 6.020142 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.2, Inf] sal aou phosphate phosphate_star 3.043592 6.021594 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.2, Inf] sal aou silicate phosphate phosphate_star 3.042204 6.020142 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.2, Inf] sal tem aou oxygen phosphate 3.035138 6.012471 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.2, Inf] sal tem aou oxygen phosphate_star 3.035138 6.012471 (1999,2012] –> (2012,Inf]
Atlantic (2012,Inf] (28.2, Inf] sal tem aou phosphate phosphate_star 3.035138 6.012471 (1999,2012] –> (2012,Inf]
Indo-Pacific (1999,2012] (-Inf,26] sal aou oxygen silicate phosphate 6.382776 17.781950 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (-Inf,26] sal aou oxygen silicate phosphate_star 6.382776 17.781950 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (-Inf,26] sal aou silicate phosphate phosphate_star 6.382776 17.781950 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (-Inf,26] sal tem aou oxygen phosphate_star 6.766883 18.511261 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (-Inf,26] sal tem aou silicate phosphate 6.500945 17.932462 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (-Inf,26] sal tem aou silicate phosphate_star 6.540434 17.989323 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (-Inf,26] sal tem oxygen silicate phosphate 6.557887 18.012918 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (-Inf,26] sal tem oxygen silicate phosphate_star 6.557887 18.012918 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (-Inf,26] sal tem silicate phosphate phosphate_star 6.557887 18.012918 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (-Inf,26] sal tem silicate phosphate_star 6.786607 18.381810 (1981,1999] –> (1999,2012]
Indo-Pacific (2012,Inf] (-Inf,26] sal aou oxygen silicate phosphate 5.431352 11.814129 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (-Inf,26] sal aou oxygen silicate phosphate_star 5.431352 11.814129 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (-Inf,26] sal aou silicate phosphate phosphate_star 5.431352 11.814129 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (-Inf,26] sal tem aou oxygen phosphate 5.833742 12.600625 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (-Inf,26] sal tem aou oxygen phosphate_star 5.833742 12.600625 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (-Inf,26] sal tem aou silicate phosphate 5.669234 12.170179 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (-Inf,26] sal tem aou silicate phosphate_star 5.738305 12.278740 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (-Inf,26] sal tem oxygen silicate phosphate 5.781974 12.339862 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (-Inf,26] sal tem oxygen silicate phosphate_star 5.781974 12.339862 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (-Inf,26] sal tem silicate phosphate phosphate_star 5.781974 12.339862 (1999,2012] –> (2012,Inf]
Indo-Pacific (1999,2012] (26,26.5] sal aou oxygen silicate phosphate 4.711026 12.184997 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (26,26.5] sal aou oxygen silicate phosphate_star 4.711026 12.184997 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (26,26.5] sal aou silicate phosphate phosphate_star 4.711026 12.184997 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (26,26.5] sal tem aou oxygen phosphate_star 4.920095 12.569779 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (26,26.5] sal tem aou silicate phosphate 4.722661 12.262082 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (26,26.5] sal tem aou silicate phosphate_star 4.750734 12.324104 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (26,26.5] sal tem oxygen silicate phosphate 4.752288 12.327336 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (26,26.5] sal tem oxygen silicate phosphate_star 4.752288 12.327336 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (26,26.5] sal tem silicate phosphate phosphate_star 4.752288 12.327336 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (26,26.5] sal tem silicate phosphate_star 4.777258 12.366672 (1981,1999] –> (1999,2012]
Indo-Pacific (2012,Inf] (26,26.5] sal aou oxygen silicate phosphate 3.902907 8.613934 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (26,26.5] sal aou oxygen silicate phosphate_star 3.902907 8.613934 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (26,26.5] sal aou silicate phosphate phosphate_star 3.902907 8.613934 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (26,26.5] sal tem aou oxygen phosphate_star 4.242711 9.162805 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (26,26.5] sal tem aou silicate phosphate 4.001587 8.724248 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (26,26.5] sal tem aou silicate phosphate_star 4.051659 8.802393 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (26,26.5] sal tem oxygen silicate phosphate 4.051596 8.803884 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (26,26.5] sal tem oxygen silicate phosphate_star 4.051596 8.803884 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (26,26.5] sal tem silicate phosphate phosphate_star 4.051596 8.803884 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (26,26.5] sal tem silicate phosphate_star 4.051660 8.828918 (1999,2012] –> (2012,Inf]
Indo-Pacific (1999,2012] (26.5,26.75] sal aou oxygen phosphate 5.294290 10.995301 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (26.5,26.75] sal aou oxygen phosphate phosphate_star 5.294290 10.995301 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (26.5,26.75] sal aou oxygen phosphate_star 5.294290 10.995301 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (26.5,26.75] sal aou oxygen silicate phosphate 5.293367 10.923436 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (26.5,26.75] sal aou oxygen silicate phosphate_star 5.293367 10.923436 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (26.5,26.75] sal aou phosphate phosphate_star 5.294290 10.995301 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (26.5,26.75] sal aou silicate phosphate phosphate_star 5.293367 10.923436 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (26.5,26.75] sal tem aou oxygen phosphate 5.293685 10.994656 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (26.5,26.75] sal tem aou oxygen phosphate_star 5.293685 10.994656 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (26.5,26.75] sal tem aou phosphate phosphate_star 5.293685 10.994656 (1981,1999] –> (1999,2012]
Indo-Pacific (2012,Inf] (26.5,26.75] sal aou oxygen phosphate 4.567155 9.861444 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (26.5,26.75] sal aou oxygen phosphate phosphate_star 4.567155 9.861444 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (26.5,26.75] sal aou oxygen phosphate_star 4.567155 9.861444 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (26.5,26.75] sal aou oxygen silicate phosphate 4.496899 9.790266 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (26.5,26.75] sal aou oxygen silicate phosphate_star 4.496899 9.790266 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (26.5,26.75] sal aou phosphate phosphate_star 4.567155 9.861444 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (26.5,26.75] sal aou silicate phosphate phosphate_star 4.496899 9.790266 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (26.5,26.75] sal tem aou oxygen phosphate 4.555157 9.848842 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (26.5,26.75] sal tem aou oxygen phosphate_star 4.555157 9.848842 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (26.5,26.75] sal tem aou phosphate phosphate_star 4.555157 9.848842 (1999,2012] –> (2012,Inf]
Indo-Pacific (1999,2012] (26.75,27] sal tem aou oxygen phosphate 4.917465 11.260619 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (26.75,27] sal tem aou oxygen phosphate_star 4.917465 11.260619 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (26.75,27] sal tem aou phosphate phosphate_star 4.917465 11.260619 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (26.75,27] tem aou oxygen phosphate 4.924678 11.267885 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (26.75,27] tem aou oxygen phosphate phosphate_star 4.924678 11.267885 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (26.75,27] tem aou oxygen phosphate_star 4.924678 11.267885 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (26.75,27] tem aou oxygen silicate phosphate 4.894896 11.214805 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (26.75,27] tem aou oxygen silicate phosphate_star 4.894896 11.214805 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (26.75,27] tem aou phosphate phosphate_star 4.924678 11.267885 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (26.75,27] tem aou silicate phosphate phosphate_star 4.894896 11.214805 (1981,1999] –> (1999,2012]
Indo-Pacific (2012,Inf] (26.75,27] sal tem aou oxygen phosphate 4.212845 9.130311 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (26.75,27] sal tem aou oxygen phosphate_star 4.212845 9.130311 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (26.75,27] sal tem aou phosphate phosphate_star 4.212845 9.130311 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (26.75,27] tem aou oxygen phosphate 4.224145 9.148823 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (26.75,27] tem aou oxygen phosphate phosphate_star 4.224145 9.148823 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (26.75,27] tem aou oxygen phosphate_star 4.224145 9.148823 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (26.75,27] tem aou oxygen silicate phosphate 4.210817 9.105713 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (26.75,27] tem aou oxygen silicate phosphate_star 4.210817 9.105713 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (26.75,27] tem aou phosphate phosphate_star 4.224145 9.148823 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (26.75,27] tem aou silicate phosphate phosphate_star 4.210817 9.105713 (1999,2012] –> (2012,Inf]
Indo-Pacific (1999,2012] (27,27.25] sal aou oxygen silicate phosphate 4.671700 10.840355 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27,27.25] sal aou oxygen silicate phosphate_star 4.671700 10.840355 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27,27.25] sal aou silicate phosphate phosphate_star 4.671700 10.840355 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27,27.25] sal tem aou oxygen phosphate 4.661474 10.695686 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27,27.25] sal tem aou oxygen phosphate_star 4.661474 10.695686 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27,27.25] sal tem aou phosphate phosphate_star 4.661474 10.695686 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27,27.25] sal tem aou silicate phosphate 4.680606 10.880070 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27,27.25] tem aou oxygen silicate phosphate 4.606714 10.559137 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27,27.25] tem aou oxygen silicate phosphate_star 4.606714 10.559137 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27,27.25] tem aou silicate phosphate phosphate_star 4.606714 10.559137 (1981,1999] –> (1999,2012]
Indo-Pacific (2012,Inf] (27,27.25] sal aou oxygen silicate phosphate 4.050682 8.722383 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27,27.25] sal aou oxygen silicate phosphate_star 4.050682 8.722383 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27,27.25] sal aou silicate phosphate phosphate_star 4.050682 8.722383 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27,27.25] sal tem aou oxygen phosphate 4.063004 8.724478 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27,27.25] sal tem aou oxygen phosphate_star 4.063004 8.724478 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27,27.25] sal tem aou phosphate phosphate_star 4.063004 8.724478 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27,27.25] sal tem aou silicate phosphate 4.068014 8.748620 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27,27.25] tem aou oxygen silicate phosphate 3.969332 8.576046 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27,27.25] tem aou oxygen silicate phosphate_star 3.969332 8.576046 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27,27.25] tem aou silicate phosphate phosphate_star 3.969332 8.576046 (1999,2012] –> (2012,Inf]
Indo-Pacific (1999,2012] (27.25,27.5] sal tem aou oxygen phosphate 4.251711 9.083022 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.25,27.5] sal tem aou oxygen phosphate_star 4.251711 9.083022 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.25,27.5] sal tem aou phosphate phosphate_star 4.251711 9.083022 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.25,27.5] tem aou oxygen phosphate 4.284580 9.133860 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.25,27.5] tem aou oxygen phosphate phosphate_star 4.284580 9.133860 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.25,27.5] tem aou oxygen phosphate_star 4.284580 9.133860 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.25,27.5] tem aou oxygen silicate phosphate 4.037922 8.834995 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.25,27.5] tem aou oxygen silicate phosphate_star 4.037922 8.834995 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.25,27.5] tem aou phosphate phosphate_star 4.284580 9.133860 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.25,27.5] tem aou silicate phosphate phosphate_star 4.037922 8.834995 (1981,1999] –> (1999,2012]
Indo-Pacific (2012,Inf] (27.25,27.5] sal aou oxygen silicate phosphate 3.224223 7.381656 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.25,27.5] sal aou oxygen silicate phosphate_star 3.224223 7.381656 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.25,27.5] sal aou silicate phosphate 3.260701 7.483265 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.25,27.5] sal aou silicate phosphate phosphate_star 3.224223 7.381656 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.25,27.5] sal tem aou silicate phosphate 3.228306 7.414579 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.25,27.5] sal tem aou silicate phosphate_star 3.236686 7.468433 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.25,27.5] sal tem oxygen silicate phosphate 3.235064 7.457907 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.25,27.5] sal tem oxygen silicate phosphate_star 3.235064 7.457907 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.25,27.5] sal tem silicate phosphate phosphate_star 3.235064 7.457907 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.25,27.5] tem aou oxygen silicate phosphate_star 3.495450 7.533372 (1999,2012] –> (2012,Inf]
Indo-Pacific (1999,2012] (27.5,27.75] sal aou oxygen silicate phosphate 3.762133 8.471531 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.5,27.75] sal aou oxygen silicate phosphate_star 3.762133 8.471531 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.5,27.75] sal aou silicate phosphate phosphate_star 3.762133 8.471531 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.5,27.75] sal tem aou oxygen phosphate 3.978677 8.455887 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.5,27.75] sal tem aou oxygen phosphate_star 3.978677 8.455887 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.5,27.75] sal tem aou phosphate phosphate_star 3.978677 8.455887 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.5,27.75] tem aou oxygen phosphate 4.036864 8.514089 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.5,27.75] tem aou oxygen silicate phosphate 3.958379 8.418432 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.5,27.75] tem aou oxygen silicate phosphate_star 3.958379 8.418432 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.5,27.75] tem aou silicate phosphate phosphate_star 3.958379 8.418432 (1981,1999] –> (1999,2012]
Indo-Pacific (2012,Inf] (27.5,27.75] sal aou oxygen silicate phosphate 2.958522 6.720656 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.5,27.75] sal aou oxygen silicate phosphate_star 2.958522 6.720656 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.5,27.75] sal aou silicate phosphate 3.098919 7.038367 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.5,27.75] sal aou silicate phosphate phosphate_star 2.958522 6.720656 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.5,27.75] sal tem aou silicate phosphate 2.957954 6.718272 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.5,27.75] sal tem aou silicate phosphate_star 2.958306 6.727781 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.5,27.75] sal tem oxygen silicate phosphate 2.958200 6.726590 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.5,27.75] sal tem oxygen silicate phosphate_star 2.958200 6.726590 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.5,27.75] sal tem silicate phosphate phosphate_star 2.958200 6.726590 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.5,27.75] sal tem silicate phosphate_star 3.053052 6.853125 (1999,2012] –> (2012,Inf]
Indo-Pacific (1999,2012] (27.75,27.85] sal tem aou oxygen phosphate 3.352787 7.458802 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.75,27.85] sal tem aou oxygen phosphate_star 3.352787 7.458802 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.75,27.85] sal tem aou phosphate phosphate_star 3.352787 7.458802 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.75,27.85] sal tem aou silicate phosphate 3.331480 7.466189 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.75,27.85] sal tem aou silicate phosphate_star 3.311349 7.464236 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.75,27.85] tem aou oxygen phosphate 3.353579 7.466826 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.75,27.85] tem aou oxygen phosphate_star 3.353579 7.466826 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.75,27.85] tem aou oxygen silicate phosphate 3.291787 7.404989 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.75,27.85] tem aou oxygen silicate phosphate_star 3.291787 7.404989 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.75,27.85] tem aou silicate phosphate phosphate_star 3.291787 7.404989 (1981,1999] –> (1999,2012]
Indo-Pacific (2012,Inf] (27.75,27.85] sal tem aou silicate phosphate 2.864436 6.195915 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.75,27.85] sal tem aou silicate phosphate_star 2.854956 6.166305 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.75,27.85] sal tem oxygen silicate phosphate 2.855040 6.165101 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.75,27.85] sal tem oxygen silicate phosphate_star 2.855040 6.165101 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.75,27.85] sal tem silicate phosphate phosphate_star 2.855040 6.165101 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.75,27.85] sal tem silicate phosphate_star 2.855337 6.189176 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.75,27.85] tem aou oxygen silicate phosphate 2.899065 6.190852 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.75,27.85] tem aou oxygen silicate phosphate_star 2.899065 6.190852 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.75,27.85] tem aou silicate phosphate phosphate_star 2.899065 6.190852 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.75,27.85] tem oxygen silicate phosphate_star 2.906132 6.221301 (1999,2012] –> (2012,Inf]
Indo-Pacific (1999,2012] (27.85,27.95] aou oxygen silicate phosphate 3.011375 7.308924 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.85,27.95] aou oxygen silicate phosphate phosphate_star 3.011375 7.308924 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.85,27.95] aou oxygen silicate phosphate_star 3.011375 7.308924 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.85,27.95] aou silicate phosphate phosphate_star 3.011375 7.308924 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.85,27.95] sal aou oxygen silicate phosphate 3.011290 7.307671 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.85,27.95] sal aou oxygen silicate phosphate_star 3.011290 7.307671 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.85,27.95] sal aou silicate phosphate phosphate_star 3.011290 7.307671 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.85,27.95] tem aou oxygen silicate phosphate 3.008174 7.305708 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.85,27.95] tem aou oxygen silicate phosphate_star 3.008174 7.305708 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.85,27.95] tem aou silicate phosphate phosphate_star 3.008174 7.305708 (1981,1999] –> (1999,2012]
Indo-Pacific (2012,Inf] (27.85,27.95] sal aou oxygen silicate phosphate 2.875658 5.886949 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.85,27.95] sal aou oxygen silicate phosphate_star 2.875658 5.886949 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.85,27.95] sal aou silicate phosphate phosphate_star 2.875658 5.886949 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.85,27.95] sal tem aou silicate phosphate 2.860993 5.870774 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.85,27.95] sal tem aou silicate phosphate_star 2.854937 5.872812 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.85,27.95] sal tem oxygen silicate phosphate 2.854944 5.873768 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.85,27.95] sal tem oxygen silicate phosphate_star 2.854944 5.873768 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.85,27.95] sal tem silicate phosphate phosphate_star 2.854944 5.873768 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.85,27.95] tem aou oxygen silicate phosphate 2.881215 5.889390 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.85,27.95] tem aou oxygen silicate phosphate_star 2.881215 5.889390 (1999,2012] –> (2012,Inf]
Indo-Pacific (1999,2012] (27.95,28.05] sal aou oxygen phosphate 2.880616 7.138803 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.95,28.05] sal aou oxygen phosphate phosphate_star 2.880616 7.138803 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.95,28.05] sal aou oxygen phosphate_star 2.880616 7.138803 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.95,28.05] sal aou oxygen silicate phosphate 2.857866 7.114416 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.95,28.05] sal aou oxygen silicate phosphate_star 2.857866 7.114416 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.95,28.05] sal aou phosphate phosphate_star 2.880616 7.138803 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.95,28.05] sal aou silicate phosphate phosphate_star 2.857866 7.114416 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.95,28.05] sal tem aou oxygen phosphate 2.872134 7.123483 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.95,28.05] sal tem aou oxygen phosphate_star 2.872134 7.123483 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (27.95,28.05] sal tem aou phosphate phosphate_star 2.872134 7.123483 (1981,1999] –> (1999,2012]
Indo-Pacific (2012,Inf] (27.95,28.05] sal aou oxygen silicate phosphate 3.056238 5.914104 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.95,28.05] sal aou oxygen silicate phosphate_star 3.056238 5.914104 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.95,28.05] sal aou silicate phosphate phosphate_star 3.056238 5.914104 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.95,28.05] sal tem aou oxygen phosphate 3.086481 5.958615 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.95,28.05] sal tem aou oxygen phosphate_star 3.086481 5.958615 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.95,28.05] sal tem aou phosphate phosphate_star 3.086481 5.958615 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.95,28.05] sal tem aou silicate phosphate 3.046334 5.932391 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.95,28.05] tem aou oxygen silicate phosphate 3.089078 5.956376 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.95,28.05] tem aou oxygen silicate phosphate_star 3.089078 5.956376 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (27.95,28.05] tem aou silicate phosphate phosphate_star 3.089078 5.956376 (1999,2012] –> (2012,Inf]
Indo-Pacific (1999,2012] (28.05,28.1] sal aou oxygen phosphate 2.845344 6.905545 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (28.05,28.1] sal aou oxygen phosphate phosphate_star 2.845344 6.905545 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (28.05,28.1] sal aou oxygen phosphate_star 2.845344 6.905545 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (28.05,28.1] sal aou oxygen silicate phosphate 2.827384 6.878933 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (28.05,28.1] sal aou oxygen silicate phosphate_star 2.827384 6.878933 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (28.05,28.1] sal aou phosphate phosphate_star 2.845344 6.905545 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (28.05,28.1] sal aou silicate phosphate phosphate_star 2.827384 6.878933 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (28.05,28.1] sal tem aou oxygen phosphate 2.842890 6.895716 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (28.05,28.1] sal tem aou oxygen phosphate_star 2.842890 6.895716 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (28.05,28.1] sal tem aou phosphate phosphate_star 2.842890 6.895716 (1981,1999] –> (1999,2012]
Indo-Pacific (2012,Inf] (28.05,28.1] sal aou oxygen phosphate 2.820324 5.665668 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (28.05,28.1] sal aou oxygen phosphate phosphate_star 2.820324 5.665668 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (28.05,28.1] sal aou oxygen phosphate_star 2.820324 5.665668 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (28.05,28.1] sal aou oxygen silicate phosphate 2.818086 5.645470 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (28.05,28.1] sal aou oxygen silicate phosphate_star 2.818086 5.645470 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (28.05,28.1] sal aou silicate phosphate phosphate_star 2.818086 5.645470 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (28.05,28.1] sal tem aou oxygen phosphate 2.809427 5.652317 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (28.05,28.1] sal tem aou oxygen phosphate_star 2.809427 5.652317 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (28.05,28.1] sal tem aou phosphate phosphate_star 2.809427 5.652317 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (28.05,28.1] sal tem aou silicate phosphate 2.827612 5.665100 (1999,2012] –> (2012,Inf]
Indo-Pacific (1999,2012] (28.1, Inf] aou oxygen silicate phosphate 3.116853 6.933099 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (28.1, Inf] aou oxygen silicate phosphate phosphate_star 3.116853 6.933099 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (28.1, Inf] aou oxygen silicate phosphate_star 3.116853 6.933099 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (28.1, Inf] aou silicate phosphate phosphate_star 3.116853 6.933099 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (28.1, Inf] sal aou oxygen silicate phosphate 3.115648 6.912688 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (28.1, Inf] sal aou oxygen silicate phosphate_star 3.115648 6.912688 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (28.1, Inf] sal aou silicate phosphate phosphate_star 3.115648 6.912688 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (28.1, Inf] tem aou oxygen silicate phosphate 3.061562 6.829327 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (28.1, Inf] tem aou oxygen silicate phosphate_star 3.061562 6.829327 (1981,1999] –> (1999,2012]
Indo-Pacific (1999,2012] (28.1, Inf] tem aou silicate phosphate phosphate_star 3.061562 6.829327 (1981,1999] –> (1999,2012]
Indo-Pacific (2012,Inf] (28.1, Inf] aou oxygen silicate phosphate 2.689332 5.806185 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (28.1, Inf] aou oxygen silicate phosphate phosphate_star 2.689332 5.806185 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (28.1, Inf] aou oxygen silicate phosphate_star 2.689332 5.806185 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (28.1, Inf] aou silicate phosphate phosphate_star 2.689332 5.806185 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (28.1, Inf] sal aou oxygen silicate phosphate 2.600926 5.716575 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (28.1, Inf] sal aou oxygen silicate phosphate_star 2.600926 5.716575 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (28.1, Inf] sal aou silicate phosphate phosphate_star 2.600926 5.716575 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (28.1, Inf] tem aou oxygen silicate phosphate 2.640982 5.702544 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (28.1, Inf] tem aou oxygen silicate phosphate_star 2.640982 5.702544 (1999,2012] –> (2012,Inf]
Indo-Pacific (2012,Inf] (28.1, Inf] tem aou silicate phosphate phosphate_star 2.640982 5.702544 (1999,2012] –> (2012,Inf]

3.4 Correlations RMSE between eras

3.4.1 All models

lm_all_stats %>%
  ggplot(aes(rmse_GO_SHIP, rmse_sum - rmse_GO_SHIP, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  facet_grid(eras ~ basin)

rm(lm_all_stats, lm_all_stats_wide)

3.4.2 10 best models

lm_best %>%
  ggplot(aes(rmse_GO_SHIP, rmse_sum - rmse_GO_SHIP, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  facet_grid(eras ~ basin)
# this code tests the rmse output of the ols_step_all_possible() function used
lm_1 <- GLODAP %>% 
  filter(basin == "Atlantic",
         gamma_slab == "(-Inf,26]",
         era == "new_era") %>% 
  lm(cstar ~ sal + tem + aou + oxygen + phosphate_star, data = .)

lm_2 <- GLODAP %>% 
  filter(basin == "Atlantic",
         gamma_slab == "(-Inf,26]",
         era == "new_era") %>% 
  lm(cstar ~ sal + tem + aou + oxygen + phosphate, data = .)

lm_1
lm_2

sqrt(mean(lm_1$residuals^2))
sqrt(mean(lm_2$residuals^2))

4 Fitting best models

After selecting 10 linear regression models with lowest summed RMSE across all eras, models are fitted again and model coefficients are saved to file.

# format model formula
lm_best <- lm_best %>% 
  mutate(lm_coeff = str_replace_all(predictors, " ", " + "),
         lm_coeff = paste("cstar ~", lm_coeff))


# define combined eras that will be analysed against each other
eras_era <-
  tibble(
    eras = c("JGOFS_GO", "GO_new", "JGOFS_GO", "GO_new"),
    era = c("JGOFS_WOCE", "GO_SHIP", "GO_SHIP", "new_era")
  )

lm_best <- full_join(lm_best, eras_era)
rm(eras_era)
# format model formula
lm_best <- lm_best %>% 
  mutate(lm_coeff = str_replace_all(predictors, " ", " + "),
         lm_coeff = paste("cstar ~", lm_coeff))


era <- sort(unique(GLODAP$era))
eras <- sort(unique(lm_best$eras))
# define combined eras that will be analysed against each other
eras_era <-
  tibble(
    eras = c(eras[1], eras[2], eras[1], eras[2]),
    era = c(era[1], era[2], era[2], era[3])
  )

lm_best <- lm_best %>% 
  select(-era)

lm_best <- full_join(lm_best, eras_era)
rm(eras_era, era, eras)
# remove model objects not required for prediction
strip_lm = function(cm) {
  cm$y = c()
  cm$model = c()
  
  cm$residuals = c()
  cm$fitted.values = c()
  cm$effects = c()
  cm$qr$qr = c()
  cm$linear.predictors = c()
  cm$weights = c()
  cm$prior.weights = c()
  cm$data = c()
  
  cm$family$variance = c()
  cm$family$dev.resids = c()
  cm$family$aic = c()
  cm$family$validmu = c()
  cm$family$simulate = c()
  attr(cm$terms, ".Environment") = c()
  attr(cm$formula, ".Environment") = c()
  
  cm
}


# select required columns
lm_best <- lm_best %>% 
  select(basin, gamma_slab, era, eras, predictors, lm_coeff, rmse_sum)

# run similar loop as for rmse calculation
# to get model coefficients for 10 best 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)
    
    lm_best_basin_era <- lm_best %>%
      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_basin_era %>%
        filter(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]
        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_eras$predictors)[2]
          print(i_predictors)
          
          lm_best_single <- lm_best_basin_era_slab_eras %>%
            filter(predictors == i_predictors) %>%
            select(lm_coeff) %>%
            pull()
          
          i_rmse_sum <- lm_best_basin_era_slab_eras %>%
            filter(predictors == i_predictors) %>%
            select(rmse_sum) %>%
            pull()
          
          lm_best_single_formula <- as.formula(lm_best_single)
          
          lm_fit_single <- lm(lm_best_single_formula,
                              data = GLODAP_basin_era_slab)
          

          if (params_local$plot_all_figures == "y") {
          # plotting model diagnostics
          p_model <- ggnostic(
            lm_fit_single,
            columnsY = c("cstar", ".fitted", ".resid"),
            title = paste("eras:", i_eras,
                          "| era:", i_era,
                          "| basin:", i_basin,
                          "| gamma slab:", i_gamma_slab,
                          "| predictors:", i_predictors)
          )

          ggsave(
            plot = p_model,
            path = paste(path_version_figures, "eMLR_diagnostics/", sep = ""),
            filename = paste(
              i_eras,
              i_era,
              i_basin,
              i_gamma_slab,
              i_predictors,
              "eMLR_cstar_residuals.png",
              sep = "_"
          ),
          width = 14,
          height = 8
          )
          
          rm(p_model)
          
          }

          # collect model coefficients in data frame
          coefficients <- tidy(lm_fit_single)
          coefficients <- coefficients %>%
            mutate(
              basin = i_basin,
              era = i_era,
              eras = i_eras,
              gamma_slab = i_gamma_slab,
              model = lm_best_single,
              rmse = i_rmse_sum
            )
          
          if (exists("lm_all")) {
            lm_all <- bind_rows(lm_all, coefficients)
          }
          
          if (!exists("lm_all")) {
            lm_all <- coefficients
          }
          
          # collect full models in data frame
          # models <- as_tibble(
          #   bind_cols(
          #     basin = i_basin,
          #     era = i_era,
          #     eras = i_eras,
          #     gamma_slab = i_gamma_slab,
          #     model = lm_best_single,
          #     rmse = i_rmse_sum))
          # 
          # lm_fit_single <- strip_lm(lm_fit_single)
          # 
          # models <- models %>% 
          #   mutate(lm = list(lm_fit_single))
          # 
          # if (exists("lm_all_models")) {
          #   lm_all_models <- bind_rows(lm_all_models, models)
          # }
          # 
          # if (!exists("lm_all_models")) {
          #   lm_all_models <- models
          # }
          
          
        }
        
        
      }
    }
    
  }
  
}


rm(lm_fit_single, coefficients, i_rmse_sum,
   GLODAP_basin_era, GLODAP_basin_era_slab,
   lm_best, lm_best_basin_era, lm_best_basin_era_slab, lm_best_basin_era_slab_eras,
   lm_best_single, lm_best_single_formula,
   i_basin, i_era, i_gamma_slab, i_predictors)


lm_all <- lm_all %>% 
  select(basin, gamma_slab, era, eras, model, term, estimate)

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

# write lm models to file
lm_all %>% write_csv(paste(path_version_data,
                           "lm_all.csv",
                           sep = ""))
# write lm models to file
# lm_all_models %>% write_rds(paste(path_version_data,
#                            "lm_all_models.rds",
#                            sep = ""))

# Prepare model coefficients for mapping of cstar
lm_all_wide_cstar <- lm_all %>% 
  pivot_wider(values_from = estimate,
              names_from = term,
              names_prefix = "coeff_",
              values_fill = 0)

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

rm(lm_all_wide_cstar)

# Prepare model coefficients differences for mapping of Cant

# select relevant columns
lm_all_long <- lm_all %>% 
  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))


# lm_all_wide <- lm_all_wide %>% 
#   mutate(JGOFS_GO = coeff_GO_SHIP - coeff_JGOFS_WOCE,
#          GO_new = coeff_new_era - coeff_GO_SHIP) %>% 
#   select(-c(coeff_JGOFS_WOCE,
#             coeff_GO_SHIP,
#             coeff_new_era))
# 
# lm_all_long <- lm_all_wide %>% 
#   pivot_longer(JGOFS_GO:GO_new, names_to = "eras_fit", values_to = "delta_coeff")
# 
# lm_all_long <- lm_all_long %>% 
#   filter(eras == eras_fit) %>% 
#   select(-eras_fit)
# 

lm_all_wide <- lm_all_long %>%
  pivot_wider(values_from = delta_coeff,
              names_from = term,
              names_prefix = "delta_coeff_",
              values_fill = 0)

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

rm(lm_all_long, lm_all_wide)

sessionInfo()
R version 4.0.3 (2020-10-10)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: openSUSE Leap 15.1

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] corrr_0.4.3      broom_0.7.2      kableExtra_1.3.1 knitr_1.30      
 [5] olsrr_0.5.3.9000 GGally_2.0.0     lubridate_1.7.9  metR_0.9.0      
 [9] scico_1.2.0      patchwork_1.1.0  collapse_1.4.2   forcats_0.5.0   
[13] stringr_1.4.0    dplyr_1.0.2      purrr_0.3.4      readr_1.4.0     
[17] tidyr_1.1.2      tibble_3.0.4     ggplot2_3.3.2    tidyverse_1.3.0 
[21] 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               labeling_0.4.2          
[22] scales_1.1.1             checkmate_2.0.0          goftest_1.2-2           
[25] digest_0.6.27            foreign_0.8-80           rmarkdown_2.5           
[28] rio_0.5.16               pkgconfig_2.0.3          htmltools_0.5.0         
[31] highr_0.8                dbplyr_1.4.4             rlang_0.4.9             
[34] readxl_1.3.1             rstudioapi_0.13          farver_2.0.3            
[37] generics_0.0.2           jsonlite_1.7.1           zip_2.1.1               
[40] car_3.0-10               magrittr_2.0.1           Matrix_1.2-18           
[43] Rcpp_1.0.5               munsell_0.5.0            fansi_0.4.1             
[46] abind_1.4-5              lifecycle_0.2.0          stringi_1.5.3           
[49] whisker_0.4              yaml_2.2.1               carData_3.0-4           
[52] plyr_1.8.6               grid_4.0.3               blob_1.2.1              
[55] parallel_4.0.3           promises_1.1.1           crayon_1.3.4            
[58] lattice_0.20-41          haven_2.3.1              hms_0.5.3               
[61] pillar_1.4.7             reprex_0.3.0             glue_1.4.2              
[64] evaluate_0.14            RcppArmadillo_0.10.1.2.0 data.table_1.13.2       
[67] modelr_0.1.8             vctrs_0.3.5              httpuv_1.5.4            
[70] cellranger_1.1.0         gtable_0.3.0             reshape_0.8.8           
[73] assertthat_0.2.1         xfun_0.18                openxlsx_4.2.3          
[76] RcppEigen_0.3.3.7.0      later_1.1.0.1            viridisLite_0.3.0       
[79] ellipsis_0.3.1           here_0.1