eMLR - model fitting

workflowr

• Summary
• Checks
• Past versions

Last updated: 2020-12-09

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: 71c63b0

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 71c63b0. 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:   code/Workflowr_project_managment.R
    Modified:   data/auxillary/params_local.rds

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

---

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

File	Version	Author	Date	Message
Rmd	71c63b0	jens-daniel-mueller	2020-12-09	rerun with variable predictor assignment
html	24a632f	jens-daniel-mueller	2020-12-07	Build site.
html	92dca91	jens-daniel-mueller	2020-12-07	Build site.
html	6a8004b	jens-daniel-mueller	2020-12-07	Build site.
html	70bf1a5	jens-daniel-mueller	2020-12-07	Build site.
html	7555355	jens-daniel-mueller	2020-12-07	Build site.
html	143d6fa	jens-daniel-mueller	2020-12-07	Build site.
html	abc6818	jens-daniel-mueller	2020-12-03	Build site.
Rmd	992ba15	jens-daniel-mueller	2020-12-03	rebuild with variable inventory depth
html	c8c2e7b	jens-daniel-mueller	2020-12-03	Build site.
Rmd	83203db	jens-daniel-mueller	2020-12-03	calculate cant with variable inventory depth
html	090e4d5	jens-daniel-mueller	2020-12-02	Build site.
html	7c25f7a	jens-daniel-mueller	2020-12-02	Build site.
html	ec8dc38	jens-daniel-mueller	2020-12-02	Build site.
html	c987de1	jens-daniel-mueller	2020-12-02	Build site.
html	f8358f8	jens-daniel-mueller	2020-12-02	Build site.
html	b03ddb8	jens-daniel-mueller	2020-12-02	Build site.
Rmd	9183e8f	jens-daniel-mueller	2020-12-02	revised assignment of era to eras
html	22d0127	jens-daniel-mueller	2020-12-01	Build site.
html	0ff728b	jens-daniel-mueller	2020-12-01	Build site.
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

---

1 Required data

Required are:

• cleaned and prepared GLODAPv2_2020 file

GLODAP <-
  read_csv(paste(path_version_data,
                 "GLODAPv2.2020_MLR_fitting_ready.csv",
                 sep = ""))

2 Fitting all models

2.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(
        paste("cstar",
              paste(params_local$MLR_predictors, collapse = " + "),
              sep = " ~ "),
        data = GLODAP_basin_era_slab)
      
      # # old version: 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)
      
      # add currect basin, era, and gamma values
      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")

Past versions of eMLR_predictor_model_quality-1.png

Version	Author	Date
24a632f	jens-daniel-mueller	2020-12-07
6a8004b	jens-daniel-mueller	2020-12-07
70bf1a5	jens-daniel-mueller	2020-12-07
7555355	jens-daniel-mueller	2020-12-07
143d6fa	jens-daniel-mueller	2020-12-07
090e4d5	jens-daniel-mueller	2020-12-02
7c25f7a	jens-daniel-mueller	2020-12-02
b03ddb8	jens-daniel-mueller	2020-12-02
91435ae	jens-daniel-mueller	2020-12-01
196be51	jens-daniel-mueller	2020-11-30

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

Past versions of AIC_RMSE-1.png

Version	Author	Date
24a632f	jens-daniel-mueller	2020-12-07
6a8004b	jens-daniel-mueller	2020-12-07
70bf1a5	jens-daniel-mueller	2020-12-07
7555355	jens-daniel-mueller	2020-12-07
143d6fa	jens-daniel-mueller	2020-12-07
090e4d5	jens-daniel-mueller	2020-12-02
7c25f7a	jens-daniel-mueller	2020-12-02
b03ddb8	jens-daniel-mueller	2020-12-02
91435ae	jens-daniel-mueller	2020-12-01
196be51	jens-daniel-mueller	2020-11-30

2.3 Predictor selection

Within each basin and slab, the 10 linear regression models with lowest summed RMSE across two adjacent (ie compared) eras were selected.

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

# calculate RMSE sum for adjacent eras
lm_all_stats_eras <- 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 sum
lm_best <- lm_all_stats_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, 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	2000-2012	(-Inf,26]	sal aou oxygen silicate phosphate	6.884738	12.968120	1982-1999 –> 2000-2012
Atlantic	2000-2012	(-Inf,26]	sal aou oxygen silicate phosphate_star	6.884738	12.968120	1982-1999 –> 2000-2012
Atlantic	2000-2012	(-Inf,26]	sal aou silicate phosphate phosphate_star	6.884738	12.968120	1982-1999 –> 2000-2012
Atlantic	2000-2012	(-Inf,26]	sal tem aou oxygen phosphate	7.036374	13.098444	1982-1999 –> 2000-2012
Atlantic	2000-2012	(-Inf,26]	sal tem aou oxygen phosphate_star	7.036374	13.098444	1982-1999 –> 2000-2012
Atlantic	2000-2012	(-Inf,26]	sal tem aou silicate phosphate	6.957753	13.006647	1982-1999 –> 2000-2012
Atlantic	2000-2012	(-Inf,26]	sal tem aou silicate phosphate_star	6.966549	13.008554	1982-1999 –> 2000-2012
Atlantic	2000-2012	(-Inf,26]	sal tem oxygen silicate phosphate	6.972745	13.003745	1982-1999 –> 2000-2012
Atlantic	2000-2012	(-Inf,26]	sal tem oxygen silicate phosphate_star	6.972745	13.003745	1982-1999 –> 2000-2012
Atlantic	2000-2012	(-Inf,26]	sal tem silicate phosphate phosphate_star	6.972745	13.003745	1982-1999 –> 2000-2012
Atlantic	2013-2019	(-Inf,26]	sal aou oxygen silicate phosphate	2.950067	9.834805	2000-2012 –> 2013-2019
Atlantic	2013-2019	(-Inf,26]	sal aou oxygen silicate phosphate_star	2.950067	9.834805	2000-2012 –> 2013-2019
Atlantic	2013-2019	(-Inf,26]	sal aou silicate phosphate phosphate_star	2.950067	9.834805	2000-2012 –> 2013-2019
Atlantic	2013-2019	(-Inf,26]	sal tem aou oxygen phosphate	2.912523	9.948897	2000-2012 –> 2013-2019
Atlantic	2013-2019	(-Inf,26]	sal tem aou oxygen phosphate_star	2.912523	9.948897	2000-2012 –> 2013-2019
Atlantic	2013-2019	(-Inf,26]	sal tem aou silicate phosphate	2.936405	9.894159	2000-2012 –> 2013-2019
Atlantic	2013-2019	(-Inf,26]	sal tem aou silicate phosphate_star	2.934164	9.900713	2000-2012 –> 2013-2019
Atlantic	2013-2019	(-Inf,26]	sal tem oxygen silicate phosphate	2.932463	9.905208	2000-2012 –> 2013-2019
Atlantic	2013-2019	(-Inf,26]	sal tem oxygen silicate phosphate_star	2.932463	9.905208	2000-2012 –> 2013-2019
Atlantic	2013-2019	(-Inf,26]	sal tem silicate phosphate phosphate_star	2.932463	9.905208	2000-2012 –> 2013-2019
Atlantic	2000-2012	(26,26.5]	sal tem aou oxygen phosphate	4.993882	9.842365	1982-1999 –> 2000-2012
Atlantic	2000-2012	(26,26.5]	sal tem aou oxygen phosphate_star	4.993882	9.842365	1982-1999 –> 2000-2012
Atlantic	2000-2012	(26,26.5]	sal tem aou phosphate phosphate_star	4.993882	9.842365	1982-1999 –> 2000-2012
Atlantic	2000-2012	(26,26.5]	tem aou oxygen phosphate	4.993882	9.850513	1982-1999 –> 2000-2012
Atlantic	2000-2012	(26,26.5]	tem aou oxygen phosphate phosphate_star	4.993882	9.850513	1982-1999 –> 2000-2012
Atlantic	2000-2012	(26,26.5]	tem aou oxygen phosphate_star	4.993882	9.850513	1982-1999 –> 2000-2012
Atlantic	2000-2012	(26,26.5]	tem aou oxygen silicate phosphate	4.992752	9.817070	1982-1999 –> 2000-2012
Atlantic	2000-2012	(26,26.5]	tem aou oxygen silicate phosphate_star	4.992752	9.817070	1982-1999 –> 2000-2012
Atlantic	2000-2012	(26,26.5]	tem aou phosphate phosphate_star	4.993882	9.850513	1982-1999 –> 2000-2012
Atlantic	2000-2012	(26,26.5]	tem aou silicate phosphate phosphate_star	4.992752	9.817070	1982-1999 –> 2000-2012
Atlantic	2013-2019	(26,26.5]	sal aou oxygen phosphate	2.914985	8.115294	2000-2012 –> 2013-2019
Atlantic	2013-2019	(26,26.5]	sal aou oxygen phosphate phosphate_star	2.914985	8.115294	2000-2012 –> 2013-2019
Atlantic	2013-2019	(26,26.5]	sal aou oxygen phosphate_star	2.914985	8.115294	2000-2012 –> 2013-2019
Atlantic	2013-2019	(26,26.5]	sal aou oxygen silicate phosphate	2.872130	8.071948	2000-2012 –> 2013-2019
Atlantic	2013-2019	(26,26.5]	sal aou oxygen silicate phosphate_star	2.872130	8.071948	2000-2012 –> 2013-2019
Atlantic	2013-2019	(26,26.5]	sal aou phosphate phosphate_star	2.914985	8.115294	2000-2012 –> 2013-2019
Atlantic	2013-2019	(26,26.5]	sal aou silicate phosphate phosphate_star	2.872130	8.071948	2000-2012 –> 2013-2019
Atlantic	2013-2019	(26,26.5]	sal tem aou oxygen phosphate	2.894016	7.887898	2000-2012 –> 2013-2019
Atlantic	2013-2019	(26,26.5]	sal tem aou oxygen phosphate_star	2.894016	7.887898	2000-2012 –> 2013-2019
Atlantic	2013-2019	(26,26.5]	sal tem aou phosphate phosphate_star	2.894016	7.887898	2000-2012 –> 2013-2019
Atlantic	2000-2012	(26.5,26.75]	sal tem aou oxygen phosphate	3.923718	7.708099	1982-1999 –> 2000-2012
Atlantic	2000-2012	(26.5,26.75]	sal tem aou oxygen phosphate_star	3.923718	7.708099	1982-1999 –> 2000-2012
Atlantic	2000-2012	(26.5,26.75]	sal tem aou phosphate phosphate_star	3.923718	7.708099	1982-1999 –> 2000-2012
Atlantic	2000-2012	(26.5,26.75]	tem aou oxygen phosphate	3.956439	7.789039	1982-1999 –> 2000-2012
Atlantic	2000-2012	(26.5,26.75]	tem aou oxygen phosphate phosphate_star	3.956439	7.789039	1982-1999 –> 2000-2012
Atlantic	2000-2012	(26.5,26.75]	tem aou oxygen phosphate_star	3.956439	7.789039	1982-1999 –> 2000-2012
Atlantic	2000-2012	(26.5,26.75]	tem aou oxygen silicate phosphate	3.947135	7.762153	1982-1999 –> 2000-2012
Atlantic	2000-2012	(26.5,26.75]	tem aou oxygen silicate phosphate_star	3.947135	7.762153	1982-1999 –> 2000-2012
Atlantic	2000-2012	(26.5,26.75]	tem aou phosphate phosphate_star	3.956439	7.789039	1982-1999 –> 2000-2012
Atlantic	2000-2012	(26.5,26.75]	tem aou silicate phosphate phosphate_star	3.947135	7.762153	1982-1999 –> 2000-2012
Atlantic	2013-2019	(26.5,26.75]	sal aou oxygen silicate phosphate	3.055061	7.072361	2000-2012 –> 2013-2019
Atlantic	2013-2019	(26.5,26.75]	sal aou oxygen silicate phosphate_star	3.055061	7.072361	2000-2012 –> 2013-2019
Atlantic	2013-2019	(26.5,26.75]	sal tem aou oxygen phosphate	3.145313	7.069030	2000-2012 –> 2013-2019
Atlantic	2013-2019	(26.5,26.75]	sal tem aou oxygen phosphate_star	3.145313	7.069030	2000-2012 –> 2013-2019
Atlantic	2013-2019	(26.5,26.75]	sal tem aou phosphate phosphate_star	3.145313	7.069030	2000-2012 –> 2013-2019
Atlantic	2013-2019	(26.5,26.75]	sal tem aou silicate phosphate	3.072481	7.013966	2000-2012 –> 2013-2019
Atlantic	2013-2019	(26.5,26.75]	sal tem aou silicate phosphate_star	3.086398	7.010048	2000-2012 –> 2013-2019
Atlantic	2013-2019	(26.5,26.75]	sal tem oxygen silicate phosphate	3.088238	7.010238	2000-2012 –> 2013-2019
Atlantic	2013-2019	(26.5,26.75]	sal tem oxygen silicate phosphate_star	3.088238	7.010238	2000-2012 –> 2013-2019
Atlantic	2013-2019	(26.5,26.75]	sal tem silicate phosphate phosphate_star	3.088238	7.010238	2000-2012 –> 2013-2019
Atlantic	2000-2012	(26.75,27]	sal tem aou oxygen phosphate	4.555858	8.931816	1982-1999 –> 2000-2012
Atlantic	2000-2012	(26.75,27]	sal tem aou oxygen phosphate_star	4.555858	8.931816	1982-1999 –> 2000-2012
Atlantic	2000-2012	(26.75,27]	sal tem aou phosphate phosphate_star	4.555858	8.931816	1982-1999 –> 2000-2012
Atlantic	2000-2012	(26.75,27]	sal tem aou phosphate_star	4.557786	8.934167	1982-1999 –> 2000-2012
Atlantic	2000-2012	(26.75,27]	sal tem aou silicate phosphate	4.557144	8.916400	1982-1999 –> 2000-2012
Atlantic	2000-2012	(26.75,27]	sal tem aou silicate phosphate_star	4.540751	8.890241	1982-1999 –> 2000-2012
Atlantic	2000-2012	(26.75,27]	sal tem oxygen silicate phosphate	4.540812	8.890501	1982-1999 –> 2000-2012
Atlantic	2000-2012	(26.75,27]	sal tem oxygen silicate phosphate_star	4.540812	8.890501	1982-1999 –> 2000-2012
Atlantic	2000-2012	(26.75,27]	sal tem phosphate phosphate_star	4.558007	8.934720	1982-1999 –> 2000-2012
Atlantic	2000-2012	(26.75,27]	sal tem silicate phosphate phosphate_star	4.540812	8.890501	1982-1999 –> 2000-2012
Atlantic	2013-2019	(26.75,27]	sal tem aou silicate phosphate	3.397961	7.955105	2000-2012 –> 2013-2019
Atlantic	2013-2019	(26.75,27]	sal tem aou silicate phosphate_star	3.389398	7.930149	2000-2012 –> 2013-2019
Atlantic	2013-2019	(26.75,27]	sal tem oxygen silicate phosphate	3.389205	7.930017	2000-2012 –> 2013-2019
Atlantic	2013-2019	(26.75,27]	sal tem oxygen silicate phosphate_star	3.389205	7.930017	2000-2012 –> 2013-2019
Atlantic	2013-2019	(26.75,27]	sal tem silicate phosphate phosphate_star	3.389205	7.930017	2000-2012 –> 2013-2019
Atlantic	2013-2019	(26.75,27]	tem aou oxygen silicate phosphate	3.439412	8.013144	2000-2012 –> 2013-2019
Atlantic	2013-2019	(26.75,27]	tem aou oxygen silicate phosphate_star	3.439412	8.013144	2000-2012 –> 2013-2019
Atlantic	2013-2019	(26.75,27]	tem aou silicate phosphate phosphate_star	3.439412	8.013144	2000-2012 –> 2013-2019
Atlantic	2013-2019	(26.75,27]	tem aou silicate phosphate_star	3.440094	8.014089	2000-2012 –> 2013-2019
Atlantic	2013-2019	(26.75,27]	tem oxygen silicate phosphate_star	3.440386	8.014411	2000-2012 –> 2013-2019
Atlantic	2000-2012	(27,27.25]	sal tem aou silicate phosphate_star	4.882571	10.104031	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27,27.25]	sal tem oxygen silicate phosphate	4.890111	10.110721	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27,27.25]	sal tem oxygen silicate phosphate_star	4.890111	10.110721	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27,27.25]	sal tem silicate phosphate phosphate_star	4.890111	10.110721	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27,27.25]	tem aou oxygen silicate phosphate	4.849632	10.118722	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27,27.25]	tem aou oxygen silicate phosphate_star	4.849632	10.118722	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27,27.25]	tem aou silicate phosphate phosphate_star	4.849632	10.118722	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27,27.25]	tem aou silicate phosphate_star	4.910129	10.182470	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27,27.25]	tem oxygen silicate phosphate	4.917940	10.189987	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27,27.25]	tem silicate phosphate phosphate_star	4.917940	10.189987	1982-1999 –> 2000-2012
Atlantic	2013-2019	(27,27.25]	sal tem aou silicate phosphate	4.519759	9.483179	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27,27.25]	sal tem aou silicate phosphate_star	4.505911	9.388483	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27,27.25]	sal tem oxygen silicate phosphate	4.505924	9.396036	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27,27.25]	sal tem oxygen silicate phosphate_star	4.505924	9.396036	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27,27.25]	sal tem silicate phosphate phosphate_star	4.505924	9.396036	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27,27.25]	tem aou oxygen silicate phosphate	4.688537	9.538168	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27,27.25]	tem aou oxygen silicate phosphate_star	4.688537	9.538168	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27,27.25]	tem aou silicate phosphate phosphate_star	4.688537	9.538168	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27,27.25]	tem aou silicate phosphate_star	4.698906	9.609035	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27,27.25]	tem oxygen silicate phosphate_star	4.698395	9.616335	2000-2012 –> 2013-2019
Atlantic	2000-2012	(27.25,27.5]	sal tem aou oxygen phosphate	4.219430	8.813273	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.25,27.5]	sal tem aou oxygen phosphate_star	4.219430	8.813273	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.25,27.5]	sal tem aou phosphate phosphate_star	4.219430	8.813273	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.25,27.5]	sal tem aou phosphate_star	4.221209	8.825775	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.25,27.5]	sal tem aou silicate phosphate_star	4.216382	8.817323	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.25,27.5]	sal tem oxygen phosphate_star	4.226925	8.828111	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.25,27.5]	sal tem oxygen silicate phosphate	4.223357	8.821369	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.25,27.5]	sal tem oxygen silicate phosphate_star	4.223357	8.821369	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.25,27.5]	sal tem phosphate phosphate_star	4.226925	8.828111	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.25,27.5]	sal tem silicate phosphate phosphate_star	4.223357	8.821369	1982-1999 –> 2000-2012
Atlantic	2013-2019	(27.25,27.5]	sal aou oxygen silicate phosphate	4.103372	8.427575	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.25,27.5]	sal aou oxygen silicate phosphate_star	4.103372	8.427575	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.25,27.5]	sal aou silicate phosphate phosphate_star	4.103372	8.427575	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.25,27.5]	sal tem aou oxygen phosphate_star	4.231432	8.450862	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.25,27.5]	sal tem aou phosphate phosphate_star	4.231432	8.450862	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.25,27.5]	sal tem aou silicate phosphate	4.082637	8.331970	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.25,27.5]	sal tem aou silicate phosphate_star	4.088200	8.304582	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.25,27.5]	sal tem oxygen silicate phosphate	4.085275	8.308632	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.25,27.5]	sal tem oxygen silicate phosphate_star	4.085275	8.308632	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.25,27.5]	sal tem silicate phosphate phosphate_star	4.085275	8.308632	2000-2012 –> 2013-2019
Atlantic	2000-2012	(27.5,27.75]	sal aou oxygen silicate phosphate	4.431535	8.757825	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.5,27.75]	sal aou oxygen silicate phosphate_star	4.431535	8.757825	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.5,27.75]	sal aou silicate phosphate phosphate_star	4.431535	8.757825	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.5,27.75]	sal tem aou silicate phosphate	4.413151	8.750810	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.5,27.75]	sal tem aou silicate phosphate_star	4.421729	8.794283	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.5,27.75]	sal tem oxygen silicate phosphate	4.412009	8.768759	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.5,27.75]	sal tem oxygen silicate phosphate_star	4.412009	8.768759	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.5,27.75]	sal tem silicate phosphate phosphate_star	4.412009	8.768759	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.5,27.75]	tem aou oxygen silicate phosphate	4.479703	8.807080	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.5,27.75]	tem aou oxygen silicate phosphate_star	4.479703	8.807080	1982-1999 –> 2000-2012
Atlantic	2013-2019	(27.5,27.75]	sal aou oxygen silicate phosphate	4.276944	8.708479	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.5,27.75]	sal aou oxygen silicate phosphate_star	4.276944	8.708479	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.5,27.75]	sal aou silicate phosphate phosphate_star	4.276944	8.708479	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.5,27.75]	sal tem aou silicate phosphate	4.272749	8.685900	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.5,27.75]	sal tem aou silicate phosphate_star	4.269174	8.690904	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.5,27.75]	sal tem oxygen silicate phosphate	4.268578	8.680587	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.5,27.75]	sal tem oxygen silicate phosphate_star	4.268578	8.680587	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.5,27.75]	sal tem silicate phosphate phosphate_star	4.268578	8.680587	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.5,27.75]	tem aou oxygen silicate phosphate	4.273208	8.752911	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.5,27.75]	tem aou silicate phosphate phosphate_star	4.273208	8.752911	2000-2012 –> 2013-2019
Atlantic	2000-2012	(27.75,27.85]	sal aou oxygen silicate phosphate	4.814695	9.117080	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.75,27.85]	sal aou oxygen silicate phosphate_star	4.814695	9.117080	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.75,27.85]	sal aou silicate phosphate	4.815040	9.117776	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.75,27.85]	sal aou silicate phosphate phosphate_star	4.814695	9.117080	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.75,27.85]	sal oxygen silicate phosphate	4.824029	9.146177	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.75,27.85]	sal tem aou silicate phosphate	4.814645	9.115078	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.75,27.85]	sal tem aou silicate phosphate_star	4.820191	9.144181	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.75,27.85]	sal tem oxygen silicate phosphate	4.815453	9.126635	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.75,27.85]	sal tem oxygen silicate phosphate_star	4.815453	9.126635	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.75,27.85]	sal tem silicate phosphate phosphate_star	4.815453	9.126635	1982-1999 –> 2000-2012
Atlantic	2013-2019	(27.75,27.85]	sal aou oxygen silicate phosphate	5.683326	10.498022	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.75,27.85]	sal aou oxygen silicate phosphate_star	5.683326	10.498022	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.75,27.85]	sal aou silicate phosphate phosphate_star	5.683326	10.498022	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.75,27.85]	sal tem aou silicate phosphate_star	5.676858	10.497049	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.75,27.85]	sal tem oxygen silicate phosphate	5.705594	10.521047	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.75,27.85]	sal tem oxygen silicate phosphate_star	5.705594	10.521047	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.75,27.85]	sal tem silicate phosphate phosphate_star	5.705594	10.521047	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.75,27.85]	tem aou oxygen silicate phosphate	5.618862	10.451622	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.75,27.85]	tem aou oxygen silicate phosphate_star	5.618862	10.451622	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.75,27.85]	tem aou silicate phosphate phosphate_star	5.618862	10.451622	2000-2012 –> 2013-2019
Atlantic	2000-2012	(27.85,27.95]	sal aou oxygen silicate phosphate	4.984814	8.694033	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.85,27.95]	sal aou oxygen silicate phosphate_star	4.984814	8.694033	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.85,27.95]	sal aou silicate phosphate phosphate_star	4.984814	8.694033	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.85,27.95]	sal aou silicate phosphate_star	5.021257	8.735509	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.85,27.95]	sal tem aou silicate phosphate	4.985912	8.705685	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.85,27.95]	sal tem aou silicate phosphate_star	4.983893	8.694306	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.85,27.95]	sal tem oxygen silicate phosphate	4.985018	8.701498	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.85,27.95]	sal tem oxygen silicate phosphate_star	4.985018	8.701498	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.85,27.95]	sal tem silicate phosphate	5.097232	8.920300	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.85,27.95]	sal tem silicate phosphate phosphate_star	4.985018	8.701498	1982-1999 –> 2000-2012
Atlantic	2013-2019	(27.85,27.95]	sal aou oxygen silicate phosphate	4.620519	9.605333	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.85,27.95]	sal aou oxygen silicate phosphate_star	4.620519	9.605333	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.85,27.95]	sal aou silicate phosphate phosphate_star	4.620519	9.605333	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.85,27.95]	sal aou silicate phosphate_star	4.797710	9.818967	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.85,27.95]	sal tem aou silicate phosphate	4.661755	9.647668	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.85,27.95]	sal tem aou silicate phosphate_star	4.639728	9.623621	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.85,27.95]	sal tem oxygen silicate phosphate	4.663302	9.648320	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.85,27.95]	sal tem oxygen silicate phosphate_star	4.663302	9.648320	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.85,27.95]	sal tem silicate phosphate	4.669321	9.766553	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.85,27.95]	sal tem silicate phosphate phosphate_star	4.663302	9.648320	2000-2012 –> 2013-2019
Atlantic	2000-2012	(27.95,28.05]	sal aou oxygen silicate phosphate	5.899227	9.447229	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.95,28.05]	sal aou oxygen silicate phosphate_star	5.899227	9.447229	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.95,28.05]	sal aou silicate phosphate phosphate_star	5.899227	9.447229	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.95,28.05]	sal aou silicate phosphate_star	5.899518	9.453076	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.95,28.05]	sal oxygen silicate phosphate	5.946829	9.595335	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.95,28.05]	sal tem aou silicate phosphate	5.927836	9.504712	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.95,28.05]	sal tem aou silicate phosphate_star	5.899333	9.450149	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.95,28.05]	sal tem oxygen silicate phosphate	5.914796	9.480740	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.95,28.05]	sal tem oxygen silicate phosphate_star	5.914796	9.480740	1982-1999 –> 2000-2012
Atlantic	2000-2012	(27.95,28.05]	sal tem silicate phosphate phosphate_star	5.914796	9.480740	1982-1999 –> 2000-2012
Atlantic	2013-2019	(27.95,28.05]	sal aou oxygen silicate phosphate	4.632333	10.531560	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.95,28.05]	sal aou oxygen silicate phosphate_star	4.632333	10.531560	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.95,28.05]	sal aou silicate phosphate phosphate_star	4.632333	10.531560	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.95,28.05]	sal aou silicate phosphate_star	4.681510	10.581028	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.95,28.05]	sal oxygen silicate phosphate	4.810498	10.757327	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.95,28.05]	sal tem aou silicate phosphate	4.728380	10.656217	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.95,28.05]	sal tem aou silicate phosphate_star	4.660985	10.560318	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.95,28.05]	sal tem oxygen silicate phosphate	4.706414	10.621209	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.95,28.05]	sal tem oxygen silicate phosphate_star	4.706414	10.621209	2000-2012 –> 2013-2019
Atlantic	2013-2019	(27.95,28.05]	sal tem silicate phosphate phosphate_star	4.706414	10.621209	2000-2012 –> 2013-2019
Atlantic	2000-2012	(28.05,28.1]	sal tem aou oxygen phosphate	6.223192	9.742593	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.05,28.1]	sal tem aou oxygen phosphate_star	6.223192	9.742593	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.05,28.1]	sal tem aou phosphate phosphate_star	6.223192	9.742593	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.05,28.1]	sal tem aou silicate phosphate_star	6.388273	9.973388	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.05,28.1]	tem aou oxygen phosphate	6.304041	9.973530	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.05,28.1]	tem aou oxygen phosphate phosphate_star	6.304041	9.973530	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.05,28.1]	tem aou oxygen phosphate_star	6.304041	9.973530	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.05,28.1]	tem aou oxygen silicate phosphate	6.228784	9.873586	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.05,28.1]	tem aou oxygen silicate phosphate_star	6.228784	9.873586	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.05,28.1]	tem aou silicate phosphate phosphate_star	6.228784	9.873586	1982-1999 –> 2000-2012
Atlantic	2013-2019	(28.05,28.1]	sal tem aou oxygen phosphate	4.139789	10.362981	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.05,28.1]	sal tem aou oxygen phosphate_star	4.139789	10.362981	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.05,28.1]	sal tem aou phosphate phosphate_star	4.139789	10.362981	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.05,28.1]	tem aou oxygen phosphate	4.204528	10.508569	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.05,28.1]	tem aou oxygen phosphate phosphate_star	4.204528	10.508569	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.05,28.1]	tem aou oxygen phosphate_star	4.204528	10.508569	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.05,28.1]	tem aou oxygen silicate phosphate	4.035564	10.264347	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.05,28.1]	tem aou oxygen silicate phosphate_star	4.035564	10.264347	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.05,28.1]	tem aou phosphate phosphate_star	4.204528	10.508569	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.05,28.1]	tem aou silicate phosphate phosphate_star	4.035564	10.264347	2000-2012 –> 2013-2019
Atlantic	2000-2012	(28.1,28.15]	sal tem aou oxygen phosphate	6.850980	10.312578	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.1,28.15]	sal tem aou oxygen phosphate_star	6.850980	10.312578	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.1,28.15]	sal tem aou phosphate phosphate_star	6.850980	10.312578	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.1,28.15]	sal tem aou phosphate_star	6.859918	10.330397	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.1,28.15]	sal tem aou silicate phosphate_star	6.853668	10.294470	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.1,28.15]	sal tem oxygen silicate phosphate	6.854783	10.294982	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.1,28.15]	sal tem oxygen silicate phosphate_star	6.854783	10.294982	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.1,28.15]	sal tem phosphate phosphate_star	6.859420	10.331939	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.1,28.15]	sal tem silicate phosphate phosphate_star	6.854783	10.294982	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.1,28.15]	sal tem silicate phosphate_star	6.856052	10.303259	1982-1999 –> 2000-2012
Atlantic	2013-2019	(28.1,28.15]	sal tem aou oxygen phosphate	4.136313	10.987293	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.1,28.15]	sal tem aou oxygen phosphate_star	4.136313	10.987293	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.1,28.15]	sal tem aou phosphate phosphate_star	4.136313	10.987293	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.1,28.15]	sal tem aou silicate phosphate_star	4.165875	11.019543	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.1,28.15]	sal tem oxygen silicate phosphate	4.170688	11.025471	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.1,28.15]	sal tem oxygen silicate phosphate_star	4.170688	11.025471	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.1,28.15]	sal tem silicate phosphate phosphate_star	4.170688	11.025471	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.1,28.15]	tem aou oxygen silicate phosphate	4.136640	11.029666	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.1,28.15]	tem aou oxygen silicate phosphate_star	4.136640	11.029666	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.1,28.15]	tem aou silicate phosphate phosphate_star	4.136640	11.029666	2000-2012 –> 2013-2019
Atlantic	2000-2012	(28.15,28.2]	sal aou oxygen silicate phosphate	3.556240	7.039722	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.15,28.2]	sal aou oxygen silicate phosphate_star	3.556240	7.039722	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.15,28.2]	sal aou silicate phosphate phosphate_star	3.556240	7.039722	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.15,28.2]	sal tem aou oxygen phosphate	3.553735	7.020686	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.15,28.2]	sal tem aou oxygen phosphate_star	3.553735	7.020686	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.15,28.2]	sal tem aou phosphate phosphate_star	3.553735	7.020686	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.15,28.2]	tem aou oxygen phosphate_star	3.556786	7.039936	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.15,28.2]	tem aou oxygen silicate phosphate	3.551970	7.034977	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.15,28.2]	tem aou oxygen silicate phosphate_star	3.551970	7.034977	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.15,28.2]	tem aou silicate phosphate phosphate_star	3.551970	7.034977	1982-1999 –> 2000-2012
Atlantic	2013-2019	(28.15,28.2]	sal aou oxygen silicate phosphate	2.832677	6.388917	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.15,28.2]	sal aou oxygen silicate phosphate_star	2.832677	6.388917	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.15,28.2]	sal tem aou oxygen phosphate	2.765004	6.318739	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.15,28.2]	sal tem aou oxygen phosphate_star	2.765004	6.318739	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.15,28.2]	sal tem aou phosphate	2.765311	6.353311	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.15,28.2]	sal tem aou phosphate phosphate_star	2.765004	6.318739	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.15,28.2]	sal tem aou silicate phosphate	2.765305	6.337464	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.15,28.2]	sal tem oxygen silicate phosphate	2.766445	6.377068	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.15,28.2]	sal tem oxygen silicate phosphate_star	2.766445	6.377068	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.15,28.2]	sal tem silicate phosphate phosphate_star	2.766445	6.377068	2000-2012 –> 2013-2019
Atlantic	2000-2012	(28.2, Inf]	sal aou oxygen phosphate	2.978129	5.827349	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.2, Inf]	sal aou oxygen phosphate phosphate_star	2.978129	5.827349	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.2, Inf]	sal aou oxygen phosphate_star	2.978129	5.827349	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.2, Inf]	sal aou oxygen silicate phosphate	2.978066	5.819691	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.2, Inf]	sal aou oxygen silicate phosphate_star	2.978066	5.819691	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.2, Inf]	sal aou phosphate phosphate_star	2.978129	5.827349	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.2, Inf]	sal aou silicate phosphate phosphate_star	2.978066	5.819691	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.2, Inf]	sal tem aou oxygen phosphate	2.977472	5.825378	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.2, Inf]	sal tem aou oxygen phosphate_star	2.977472	5.825378	1982-1999 –> 2000-2012
Atlantic	2000-2012	(28.2, Inf]	sal tem aou phosphate phosphate_star	2.977472	5.825378	1982-1999 –> 2000-2012
Atlantic	2013-2019	(28.2, Inf]	sal aou oxygen phosphate	3.043573	6.021702	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.2, Inf]	sal aou oxygen phosphate phosphate_star	3.043573	6.021702	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.2, Inf]	sal aou oxygen phosphate_star	3.043573	6.021702	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.2, Inf]	sal aou oxygen silicate phosphate	3.042183	6.020250	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.2, Inf]	sal aou oxygen silicate phosphate_star	3.042183	6.020250	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.2, Inf]	sal aou phosphate phosphate_star	3.043573	6.021702	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.2, Inf]	sal aou silicate phosphate phosphate_star	3.042183	6.020250	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.2, Inf]	sal tem aou oxygen phosphate	3.035108	6.012580	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.2, Inf]	sal tem aou oxygen phosphate_star	3.035108	6.012580	2000-2012 –> 2013-2019
Atlantic	2013-2019	(28.2, Inf]	sal tem aou phosphate phosphate_star	3.035108	6.012580	2000-2012 –> 2013-2019
Indo-Pacific	2000-2012	(-Inf,26]	sal aou oxygen silicate phosphate	6.380598	17.782812	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(-Inf,26]	sal aou oxygen silicate phosphate_star	6.380598	17.782812	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(-Inf,26]	sal aou silicate phosphate phosphate_star	6.380598	17.782812	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(-Inf,26]	sal tem aou phosphate phosphate_star	6.763521	18.511314	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(-Inf,26]	sal tem aou silicate phosphate	6.499456	17.934014	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(-Inf,26]	sal tem aou silicate phosphate_star	6.539057	17.990971	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(-Inf,26]	sal tem oxygen silicate phosphate	6.556489	18.014536	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(-Inf,26]	sal tem oxygen silicate phosphate_star	6.556489	18.014536	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(-Inf,26]	sal tem silicate phosphate phosphate_star	6.556489	18.014536	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(-Inf,26]	sal tem silicate phosphate_star	6.783423	18.381306	1982-1999 –> 2000-2012
Indo-Pacific	2013-2019	(-Inf,26]	sal aou oxygen silicate phosphate	5.432557	11.813156	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(-Inf,26]	sal aou oxygen silicate phosphate_star	5.432557	11.813156	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(-Inf,26]	sal aou silicate phosphate phosphate_star	5.432557	11.813156	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(-Inf,26]	sal tem aou oxygen phosphate	5.835285	12.598806	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(-Inf,26]	sal tem aou oxygen phosphate_star	5.835285	12.598806	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(-Inf,26]	sal tem aou silicate phosphate	5.670743	12.170198	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(-Inf,26]	sal tem aou silicate phosphate_star	5.739862	12.278920	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(-Inf,26]	sal tem oxygen silicate phosphate	5.783571	12.340060	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(-Inf,26]	sal tem oxygen silicate phosphate_star	5.783571	12.340060	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(-Inf,26]	sal tem silicate phosphate phosphate_star	5.783571	12.340060	2000-2012 –> 2013-2019
Indo-Pacific	2000-2012	(26,26.5]	sal aou oxygen silicate phosphate	4.711168	12.185750	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(26,26.5]	sal aou oxygen silicate phosphate_star	4.711168	12.185750	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(26,26.5]	sal aou silicate phosphate phosphate_star	4.711168	12.185750	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(26,26.5]	sal tem aou oxygen phosphate	4.921216	12.572455	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(26,26.5]	sal tem aou silicate phosphate	4.723872	12.263838	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(26,26.5]	sal tem aou silicate phosphate_star	4.752327	12.326215	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(26,26.5]	sal tem oxygen silicate phosphate	4.753898	12.329467	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(26,26.5]	sal tem oxygen silicate phosphate_star	4.753898	12.329467	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(26,26.5]	sal tem silicate phosphate phosphate_star	4.753898	12.329467	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(26,26.5]	sal tem silicate phosphate_star	4.778870	12.368894	1982-1999 –> 2000-2012
Indo-Pacific	2013-2019	(26,26.5]	sal aou oxygen silicate phosphate	3.905117	8.616286	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(26,26.5]	sal aou oxygen silicate phosphate_star	3.905117	8.616286	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(26,26.5]	sal aou silicate phosphate phosphate_star	3.905117	8.616286	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(26,26.5]	sal tem aou phosphate phosphate_star	4.245285	9.166501	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(26,26.5]	sal tem aou silicate phosphate	4.004098	8.727969	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(26,26.5]	sal tem aou silicate phosphate_star	4.054300	8.806626	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(26,26.5]	sal tem oxygen silicate phosphate	4.054256	8.808154	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(26,26.5]	sal tem oxygen silicate phosphate_star	4.054256	8.808154	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(26,26.5]	sal tem silicate phosphate phosphate_star	4.054256	8.808154	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(26,26.5]	sal tem silicate phosphate_star	4.054303	8.833174	2000-2012 –> 2013-2019
Indo-Pacific	2000-2012	(26.5,26.75]	sal aou oxygen phosphate	5.295303	10.995513	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(26.5,26.75]	sal aou oxygen phosphate phosphate_star	5.295303	10.995513	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(26.5,26.75]	sal aou oxygen phosphate_star	5.295303	10.995513	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(26.5,26.75]	sal aou oxygen silicate phosphate	5.294420	10.923243	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(26.5,26.75]	sal aou oxygen silicate phosphate_star	5.294420	10.923243	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(26.5,26.75]	sal aou phosphate phosphate_star	5.295303	10.995513	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(26.5,26.75]	sal aou silicate phosphate phosphate_star	5.294420	10.923243	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(26.5,26.75]	sal tem aou oxygen phosphate	5.294693	10.994855	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(26.5,26.75]	sal tem aou oxygen phosphate_star	5.294693	10.994855	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(26.5,26.75]	sal tem aou phosphate phosphate_star	5.294693	10.994855	1982-1999 –> 2000-2012
Indo-Pacific	2013-2019	(26.5,26.75]	sal aou oxygen phosphate	4.567974	9.863277	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(26.5,26.75]	sal aou oxygen phosphate phosphate_star	4.567974	9.863277	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(26.5,26.75]	sal aou oxygen phosphate_star	4.567974	9.863277	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(26.5,26.75]	sal aou oxygen silicate phosphate	4.497617	9.792037	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(26.5,26.75]	sal aou oxygen silicate phosphate_star	4.497617	9.792037	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(26.5,26.75]	sal aou phosphate phosphate_star	4.567974	9.863277	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(26.5,26.75]	sal aou silicate phosphate phosphate_star	4.497617	9.792037	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(26.5,26.75]	sal tem aou oxygen phosphate	4.556095	9.850788	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(26.5,26.75]	sal tem aou oxygen phosphate_star	4.556095	9.850788	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(26.5,26.75]	sal tem aou phosphate phosphate_star	4.556095	9.850788	2000-2012 –> 2013-2019
Indo-Pacific	2000-2012	(26.75,27]	sal tem aou oxygen phosphate	4.923257	11.266423	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(26.75,27]	sal tem aou oxygen phosphate_star	4.923257	11.266423	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(26.75,27]	sal tem aou phosphate phosphate_star	4.923257	11.266423	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(26.75,27]	tem aou oxygen phosphate	4.930573	11.273781	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(26.75,27]	tem aou oxygen phosphate phosphate_star	4.930573	11.273781	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(26.75,27]	tem aou oxygen phosphate_star	4.930573	11.273781	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(26.75,27]	tem aou oxygen silicate phosphate	4.900828	11.220874	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(26.75,27]	tem aou oxygen silicate phosphate_star	4.900828	11.220874	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(26.75,27]	tem aou phosphate phosphate_star	4.930573	11.273781	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(26.75,27]	tem aou silicate phosphate phosphate_star	4.900828	11.220874	1982-1999 –> 2000-2012
Indo-Pacific	2013-2019	(26.75,27]	sal tem aou oxygen phosphate	4.213406	9.136663	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(26.75,27]	sal tem aou oxygen phosphate_star	4.213406	9.136663	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(26.75,27]	sal tem aou phosphate phosphate_star	4.213406	9.136663	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(26.75,27]	tem aou oxygen phosphate	4.224647	9.155220	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(26.75,27]	tem aou oxygen phosphate phosphate_star	4.224647	9.155220	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(26.75,27]	tem aou oxygen phosphate_star	4.224647	9.155220	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(26.75,27]	tem aou oxygen silicate phosphate	4.211624	9.112452	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(26.75,27]	tem aou oxygen silicate phosphate_star	4.211624	9.112452	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(26.75,27]	tem aou phosphate phosphate_star	4.224647	9.155220	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(26.75,27]	tem aou silicate phosphate phosphate_star	4.211624	9.112452	2000-2012 –> 2013-2019
Indo-Pacific	2000-2012	(27,27.25]	sal aou oxygen silicate phosphate	4.676321	10.848347	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27,27.25]	sal aou oxygen silicate phosphate_star	4.676321	10.848347	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27,27.25]	sal aou silicate phosphate phosphate_star	4.676321	10.848347	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27,27.25]	sal tem aou oxygen phosphate	4.665711	10.703915	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27,27.25]	sal tem aou oxygen phosphate_star	4.665711	10.703915	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27,27.25]	sal tem aou phosphate phosphate_star	4.665711	10.703915	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27,27.25]	sal tem aou silicate phosphate	4.685310	10.888076	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27,27.25]	tem aou oxygen silicate phosphate	4.611332	10.567380	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27,27.25]	tem aou oxygen silicate phosphate_star	4.611332	10.567380	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27,27.25]	tem aou silicate phosphate phosphate_star	4.611332	10.567380	1982-1999 –> 2000-2012
Indo-Pacific	2013-2019	(27,27.25]	sal aou oxygen silicate phosphate	4.050835	8.727156	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27,27.25]	sal aou oxygen silicate phosphate_star	4.050835	8.727156	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27,27.25]	sal aou silicate phosphate phosphate_star	4.050835	8.727156	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27,27.25]	sal tem aou oxygen phosphate	4.062458	8.728169	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27,27.25]	sal tem aou oxygen phosphate_star	4.062458	8.728169	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27,27.25]	sal tem aou phosphate phosphate_star	4.062458	8.728169	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27,27.25]	sal tem aou silicate phosphate	4.068251	8.753561	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27,27.25]	tem aou oxygen silicate phosphate	3.969093	8.580425	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27,27.25]	tem aou oxygen silicate phosphate_star	3.969093	8.580425	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27,27.25]	tem aou silicate phosphate phosphate_star	3.969093	8.580425	2000-2012 –> 2013-2019
Indo-Pacific	2000-2012	(27.25,27.5]	sal tem aou oxygen phosphate	4.253828	9.088167	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.25,27.5]	sal tem aou oxygen phosphate_star	4.253828	9.088167	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.25,27.5]	sal tem aou phosphate phosphate_star	4.253828	9.088167	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.25,27.5]	tem aou oxygen phosphate	4.286770	9.139069	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.25,27.5]	tem aou oxygen phosphate phosphate_star	4.286770	9.139069	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.25,27.5]	tem aou oxygen phosphate_star	4.286770	9.139069	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.25,27.5]	tem aou oxygen silicate phosphate	4.040449	8.840374	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.25,27.5]	tem aou oxygen silicate phosphate_star	4.040449	8.840374	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.25,27.5]	tem aou phosphate phosphate_star	4.286770	9.139069	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.25,27.5]	tem aou silicate phosphate phosphate_star	4.040449	8.840374	1982-1999 –> 2000-2012
Indo-Pacific	2013-2019	(27.25,27.5]	sal aou oxygen silicate phosphate	3.224209	7.384177	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.25,27.5]	sal aou oxygen silicate phosphate_star	3.224209	7.384177	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.25,27.5]	sal aou silicate phosphate	3.260740	7.485827	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.25,27.5]	sal aou silicate phosphate phosphate_star	3.224209	7.384177	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.25,27.5]	sal tem aou silicate phosphate	3.228287	7.417095	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.25,27.5]	sal tem aou silicate phosphate_star	3.236659	7.470930	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.25,27.5]	sal tem oxygen silicate phosphate	3.235040	7.460411	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.25,27.5]	sal tem oxygen silicate phosphate_star	3.235040	7.460411	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.25,27.5]	sal tem silicate phosphate phosphate_star	3.235040	7.460411	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.25,27.5]	tem aou oxygen silicate phosphate_star	3.495296	7.535745	2000-2012 –> 2013-2019
Indo-Pacific	2000-2012	(27.5,27.75]	sal aou oxygen silicate phosphate	3.763279	8.474376	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.5,27.75]	sal aou oxygen silicate phosphate_star	3.763279	8.474376	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.5,27.75]	sal aou silicate phosphate phosphate_star	3.763279	8.474376	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.5,27.75]	sal tem aou oxygen phosphate	3.979674	8.458272	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.5,27.75]	sal tem aou oxygen phosphate_star	3.979674	8.458272	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.5,27.75]	sal tem aou phosphate phosphate_star	3.979674	8.458272	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.5,27.75]	tem aou oxygen phosphate_star	4.037971	8.516585	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.5,27.75]	tem aou oxygen silicate phosphate	3.959612	8.420992	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.5,27.75]	tem aou oxygen silicate phosphate_star	3.959612	8.420992	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.5,27.75]	tem aou silicate phosphate phosphate_star	3.959612	8.420992	1982-1999 –> 2000-2012
Indo-Pacific	2013-2019	(27.5,27.75]	sal aou oxygen silicate phosphate	2.958589	6.721868	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.5,27.75]	sal aou oxygen silicate phosphate_star	2.958589	6.721868	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.5,27.75]	sal aou silicate phosphate	3.098905	7.039384	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.5,27.75]	sal aou silicate phosphate phosphate_star	2.958589	6.721868	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.5,27.75]	sal tem aou silicate phosphate	2.958032	6.719449	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.5,27.75]	sal tem aou silicate phosphate_star	2.958393	6.728916	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.5,27.75]	sal tem oxygen silicate phosphate	2.958286	6.727728	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.5,27.75]	sal tem oxygen silicate phosphate_star	2.958286	6.727728	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.5,27.75]	sal tem silicate phosphate phosphate_star	2.958286	6.727728	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.5,27.75]	sal tem silicate phosphate_star	3.053046	6.854119	2000-2012 –> 2013-2019
Indo-Pacific	2000-2012	(27.75,27.85]	sal tem aou oxygen phosphate	3.352937	7.459637	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.75,27.85]	sal tem aou oxygen phosphate_star	3.352937	7.459637	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.75,27.85]	sal tem aou phosphate phosphate_star	3.352937	7.459637	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.75,27.85]	sal tem aou silicate phosphate	3.331757	7.467104	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.75,27.85]	sal tem aou silicate phosphate_star	3.311605	7.465112	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.75,27.85]	tem aou oxygen phosphate_star	3.353727	7.467665	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.75,27.85]	tem aou oxygen silicate phosphate	3.291977	7.405869	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.75,27.85]	tem aou oxygen silicate phosphate_star	3.291977	7.405869	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.75,27.85]	tem aou phosphate phosphate_star	3.353727	7.467665	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.75,27.85]	tem aou silicate phosphate phosphate_star	3.291977	7.405869	1982-1999 –> 2000-2012
Indo-Pacific	2013-2019	(27.75,27.85]	sal tem aou silicate phosphate	2.864350	6.196107	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.75,27.85]	sal tem aou silicate phosphate_star	2.854896	6.166501	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.75,27.85]	sal tem oxygen silicate phosphate	2.854981	6.165294	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.75,27.85]	sal tem oxygen silicate phosphate_star	2.854981	6.165294	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.75,27.85]	sal tem silicate phosphate phosphate_star	2.854981	6.165294	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.75,27.85]	sal tem silicate phosphate_star	2.855286	6.189422	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.75,27.85]	tem aou oxygen silicate phosphate	2.899088	6.191064	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.75,27.85]	tem aou oxygen silicate phosphate_star	2.899088	6.191064	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.75,27.85]	tem aou silicate phosphate phosphate_star	2.899088	6.191064	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.75,27.85]	tem silicate phosphate phosphate_star	2.906081	6.221501	2000-2012 –> 2013-2019
Indo-Pacific	2000-2012	(27.85,27.95]	aou oxygen silicate phosphate	3.011813	7.309681	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.85,27.95]	aou oxygen silicate phosphate phosphate_star	3.011813	7.309681	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.85,27.95]	aou oxygen silicate phosphate_star	3.011813	7.309681	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.85,27.95]	aou silicate phosphate phosphate_star	3.011813	7.309681	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.85,27.95]	sal aou oxygen silicate phosphate	3.011727	7.308427	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.85,27.95]	sal aou oxygen silicate phosphate_star	3.011727	7.308427	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.85,27.95]	sal aou silicate phosphate phosphate_star	3.011727	7.308427	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.85,27.95]	tem aou oxygen silicate phosphate	3.008640	7.306492	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.85,27.95]	tem aou oxygen silicate phosphate_star	3.008640	7.306492	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.85,27.95]	tem aou silicate phosphate phosphate_star	3.008640	7.306492	1982-1999 –> 2000-2012
Indo-Pacific	2013-2019	(27.85,27.95]	sal aou oxygen silicate phosphate	2.875628	5.887356	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.85,27.95]	sal aou oxygen silicate phosphate_star	2.875628	5.887356	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.85,27.95]	sal aou silicate phosphate phosphate_star	2.875628	5.887356	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.85,27.95]	sal tem aou silicate phosphate	2.861003	5.871267	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.85,27.95]	sal tem aou silicate phosphate_star	2.854969	5.873350	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.85,27.95]	sal tem oxygen silicate phosphate	2.854977	5.874309	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.85,27.95]	sal tem oxygen silicate phosphate_star	2.854977	5.874309	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.85,27.95]	sal tem silicate phosphate phosphate_star	2.854977	5.874309	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.85,27.95]	tem aou oxygen silicate phosphate	2.881267	5.889907	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.85,27.95]	tem aou silicate phosphate phosphate_star	2.881267	5.889907	2000-2012 –> 2013-2019
Indo-Pacific	2000-2012	(27.95,28.05]	sal aou oxygen phosphate	2.881016	7.139310	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.95,28.05]	sal aou oxygen phosphate phosphate_star	2.881016	7.139310	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.95,28.05]	sal aou oxygen phosphate_star	2.881016	7.139310	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.95,28.05]	sal aou oxygen silicate phosphate	2.858246	7.114908	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.95,28.05]	sal aou oxygen silicate phosphate_star	2.858246	7.114908	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.95,28.05]	sal aou phosphate phosphate_star	2.881016	7.139310	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.95,28.05]	sal aou silicate phosphate phosphate_star	2.858246	7.114908	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.95,28.05]	sal tem aou oxygen phosphate	2.872498	7.123971	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.95,28.05]	sal tem aou oxygen phosphate_star	2.872498	7.123971	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(27.95,28.05]	sal tem aou phosphate phosphate_star	2.872498	7.123971	1982-1999 –> 2000-2012
Indo-Pacific	2013-2019	(27.95,28.05]	sal aou oxygen silicate phosphate	3.056329	5.914575	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.95,28.05]	sal aou oxygen silicate phosphate_star	3.056329	5.914575	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.95,28.05]	sal aou silicate phosphate phosphate_star	3.056329	5.914575	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.95,28.05]	sal tem aou oxygen phosphate	3.086601	5.959099	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.95,28.05]	sal tem aou oxygen phosphate_star	3.086601	5.959099	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.95,28.05]	sal tem aou phosphate phosphate_star	3.086601	5.959099	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.95,28.05]	sal tem aou silicate phosphate	3.046437	5.932911	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.95,28.05]	tem aou oxygen silicate phosphate	3.089184	5.956841	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.95,28.05]	tem aou oxygen silicate phosphate_star	3.089184	5.956841	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(27.95,28.05]	tem aou silicate phosphate phosphate_star	3.089184	5.956841	2000-2012 –> 2013-2019
Indo-Pacific	2000-2012	(28.05,28.1]	sal aou oxygen phosphate	2.845584	6.905906	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(28.05,28.1]	sal aou oxygen phosphate phosphate_star	2.845584	6.905906	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(28.05,28.1]	sal aou oxygen phosphate_star	2.845584	6.905906	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(28.05,28.1]	sal aou oxygen silicate phosphate	2.827613	6.879271	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(28.05,28.1]	sal aou oxygen silicate phosphate_star	2.827613	6.879271	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(28.05,28.1]	sal aou phosphate phosphate_star	2.845584	6.905906	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(28.05,28.1]	sal aou silicate phosphate phosphate_star	2.827613	6.879271	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(28.05,28.1]	sal tem aou oxygen phosphate	2.843122	6.896059	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(28.05,28.1]	sal tem aou oxygen phosphate_star	2.843122	6.896059	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(28.05,28.1]	sal tem aou phosphate phosphate_star	2.843122	6.896059	1982-1999 –> 2000-2012
Indo-Pacific	2013-2019	(28.05,28.1]	sal aou oxygen phosphate	2.820374	5.665957	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(28.05,28.1]	sal aou oxygen phosphate phosphate_star	2.820374	5.665957	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(28.05,28.1]	sal aou oxygen phosphate_star	2.820374	5.665957	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(28.05,28.1]	sal aou oxygen silicate phosphate	2.818140	5.645753	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(28.05,28.1]	sal aou oxygen silicate phosphate_star	2.818140	5.645753	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(28.05,28.1]	sal aou silicate phosphate phosphate_star	2.818140	5.645753	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(28.05,28.1]	sal tem aou oxygen phosphate	2.809469	5.652591	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(28.05,28.1]	sal tem aou oxygen phosphate_star	2.809469	5.652591	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(28.05,28.1]	sal tem aou phosphate phosphate_star	2.809469	5.652591	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(28.05,28.1]	sal tem aou silicate phosphate	2.827676	5.665405	2000-2012 –> 2013-2019
Indo-Pacific	2000-2012	(28.1, Inf]	aou oxygen silicate phosphate	3.117577	6.933929	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(28.1, Inf]	aou oxygen silicate phosphate phosphate_star	3.117577	6.933929	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(28.1, Inf]	aou oxygen silicate phosphate_star	3.117577	6.933929	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(28.1, Inf]	aou silicate phosphate phosphate_star	3.117577	6.933929	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(28.1, Inf]	sal aou oxygen silicate phosphate	3.116374	6.913547	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(28.1, Inf]	sal aou oxygen silicate phosphate_star	3.116374	6.913547	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(28.1, Inf]	sal aou silicate phosphate phosphate_star	3.116374	6.913547	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(28.1, Inf]	tem aou oxygen silicate phosphate	3.062233	6.830092	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(28.1, Inf]	tem aou oxygen silicate phosphate_star	3.062233	6.830092	1982-1999 –> 2000-2012
Indo-Pacific	2000-2012	(28.1, Inf]	tem aou silicate phosphate phosphate_star	3.062233	6.830092	1982-1999 –> 2000-2012
Indo-Pacific	2013-2019	(28.1, Inf]	aou oxygen silicate phosphate	2.689066	5.806643	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(28.1, Inf]	aou oxygen silicate phosphate phosphate_star	2.689066	5.806643	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(28.1, Inf]	aou oxygen silicate phosphate_star	2.689066	5.806643	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(28.1, Inf]	aou silicate phosphate phosphate_star	2.689066	5.806643	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(28.1, Inf]	sal aou oxygen silicate phosphate	2.600650	5.717024	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(28.1, Inf]	sal aou oxygen silicate phosphate_star	2.600650	5.717024	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(28.1, Inf]	sal aou silicate phosphate phosphate_star	2.600650	5.717024	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(28.1, Inf]	tem aou oxygen silicate phosphate	2.640725	5.702958	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(28.1, Inf]	tem aou oxygen silicate phosphate_star	2.640725	5.702958	2000-2012 –> 2013-2019
Indo-Pacific	2013-2019	(28.1, Inf]	tem aou silicate phosphate phosphate_star	2.640725	5.702958	2000-2012 –> 2013-2019

2.4 Correlations RMSE between eras

2.4.1 All models

max_rmse <-
  max(c(lm_all_stats_eras$rmse,
        lm_all_stats_eras$rmse_sum - lm_all_stats_eras$rmse))

lm_all_stats_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)

Past versions of rmse_correlation_across_eras_all_models-1.png

Version	Author	Date
24a632f	jens-daniel-mueller	2020-12-07
6a8004b	jens-daniel-mueller	2020-12-07
70bf1a5	jens-daniel-mueller	2020-12-07
7555355	jens-daniel-mueller	2020-12-07
143d6fa	jens-daniel-mueller	2020-12-07
090e4d5	jens-daniel-mueller	2020-12-02
7c25f7a	jens-daniel-mueller	2020-12-02
b03ddb8	jens-daniel-mueller	2020-12-02
91435ae	jens-daniel-mueller	2020-12-01
196be51	jens-daniel-mueller	2020-11-30

rm(max_rmse)

2.4.2 Best models

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)

Past versions of rmse_correlation_across_eras_10_models-1.png

Version	Author	Date
24a632f	jens-daniel-mueller	2020-12-07
6a8004b	jens-daniel-mueller	2020-12-07
70bf1a5	jens-daniel-mueller	2020-12-07
7555355	jens-daniel-mueller	2020-12-07
143d6fa	jens-daniel-mueller	2020-12-07
090e4d5	jens-daniel-mueller	2020-12-02
7c25f7a	jens-daniel-mueller	2020-12-02
b03ddb8	jens-daniel-mueller	2020-12-02
91435ae	jens-daniel-mueller	2020-12-01
196be51	jens-daniel-mueller	2020-11-30

rm(max_rmse)

# 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))

3 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))


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

eras_backward <- lm_all_stats %>%
  arrange(era) %>% 
  group_by(basin, gamma_slab, predictors) %>% 
  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 <- lm_best %>% 
  select(-era)

lm_best <- full_join(lm_best, eras_era)
rm(eras_era, eras_forward, eras_backward, lm_all_stats, lm_all_stats_eras)

# 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)
          

          # plotting model diagnostics
          if (params_local$plot_all_figures == "y") {
          
          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(
              "MLR_residuals",
              i_eras,
              i_era,
              i_basin,
              i_gamma_slab,
              i_predictors,
              "predictors.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
          }
          
        }
        
        
      }
    }
    
  }
  
}


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)

3.1 Predictor counts

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

lm_all_stats <- lm_all %>% 
  filter(term != "(Intercept)") %>% 
  group_by(basin, eras, gamma_slab) %>% 
  count(term) %>% 
  ungroup() %>% 
  pivot_wider(values_from = n, names_from = term)


lm_all_stats %>%
  gt(rowname_col = "gamma_slab",
     groupname_col = c("basin", "eras")) %>% 
  summary_rows(
    groups = TRUE,
    fns = list(total = "sum")
  )

	aou	oxygen	phosphate	phosphate_star	sal	silicate	tem
Atlantic - 1982-1999 --> 2000-2012
(-Inf,26]	14	12	12	12	20	16	14
(26,26.5]	20	14	14	14	6	6	20
(26.5,26.75]	20	14	14	14	6	6	20
(26.75,27]	12	8	12	14	20	10	20
(27,27.25]	10	10	12	14	8	20	20
(27.25,27.5]	10	10	10	16	20	8	20
(27.5,27.75]	14	12	12	12	16	20	14
(27.75,27.85]	12	10	14	10	20	20	10
(27.85,27.95]	12	8	12	12	20	20	12
(27.95,28.05]	12	10	12	12	20	20	10
(28.05,28.1]	20	14	12	14	8	8	20
(28.1,28.15]	10	8	10	16	20	10	20
(28.15,28.2]	20	14	12	14	12	12	14
(28.2, Inf]	20	14	14	14	20	6	6
total	206.00	158.00	172.00	188.00	216.00	182.00	220.00
Atlantic - 2000-2012 --> 2013-2019
(-Inf,26]	14	12	12	12	20	16	14
(26,26.5]	20	14	14	14	20	6	6
(26.5,26.75]	14	12	12	12	20	14	16
(26.75,27]	12	10	10	14	10	20	20
(27,27.25]	12	10	10	14	10	20	20
(27.25,27.5]	14	10	12	14	20	16	14
(27.5,27.75]	14	10	14	12	16	20	14
(27.75,27.85]	14	12	12	14	14	20	14
(27.85,27.95]	12	8	12	12	20	20	12
(27.95,28.05]	12	10	12	12	20	20	10
(28.05,28.1]	20	14	14	14	6	6	20
(28.1,28.15]	14	12	12	14	14	14	20
(28.15,28.2]	14	12	14	10	20	12	16
(28.2, Inf]	20	14	14	14	20	6	6
total	206.00	160.00	174.00	182.00	230.00	210.00	202.00
Indo-Pacific - 1982-1999 --> 2000-2012
(-Inf,26]	12	8	12	14	20	18	14
(26,26.5]	12	10	12	12	20	18	14
(26.5,26.75]	20	14	14	14	20	6	6
(26.75,27]	20	14	14	14	6	6	20
(27,27.25]	20	12	14	12	14	14	14
(27.25,27.5]	20	14	14	14	6	6	20
(27.5,27.75]	20	14	12	14	12	12	14
(27.75,27.85]	20	10	12	14	10	10	20
(27.85,27.95]	20	14	14	14	6	20	6
(27.95,28.05]	20	14	14	14	20	6	6
(28.05,28.1]	20	14	14	14	20	6	6
(28.1, Inf]	20	14	14	14	6	20	6
total	224.00	152.00	160.00	164.00	160.00	142.00	146.00
Indo-Pacific - 2000-2012 --> 2013-2019
(-Inf,26]	14	12	12	12	20	16	14
(26,26.5]	12	8	12	14	20	18	14
(26.5,26.75]	20	14	14	14	20	6	6
(26.75,27]	20	14	14	14	6	6	20
(27,27.25]	20	12	14	12	14	14	14
(27.25,27.5]	14	10	12	12	18	20	12
(27.5,27.75]	12	8	12	12	20	20	12
(27.75,27.85]	10	8	12	14	12	20	20
(27.85,27.95]	14	10	14	12	16	20	14
(27.95,28.05]	20	12	14	12	14	14	14
(28.05,28.1]	20	14	14	12	20	8	8
(28.1, Inf]	20	14	14	14	6	20	6
total	196.00	136.00	158.00	154.00	186.00	182.00	154.00

3.2 Prepare coeffcients

Coefficients are prepared for the mapping of cstar and cant.

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


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


# Prepare model coefficients differences for mapping of cant
# this requires to subtract coefficients of adjacent era
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_cant <- lm_all_long %>%
  pivot_wider(values_from = delta_coeff,
              names_from = term,
              names_prefix = "delta_coeff_",
              values_fill = 0)

3.3 Write files

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

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

Session information

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] gt_0.2.2         corrr_0.4.3      broom_0.7.2      kableExtra_1.3.1
 [5] knitr_1.30       olsrr_0.5.3.9000 GGally_2.0.0     lubridate_1.7.9 
 [9] metR_0.9.0       scico_1.2.0      patchwork_1.1.0  collapse_1.4.2  
[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.2   
[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.9              readxl_1.3.1             rstudioapi_0.13         
[37] farver_2.0.3             generics_0.0.2           jsonlite_1.7.1          
[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.0
[67] data.table_1.13.2        modelr_0.1.8             vctrs_0.3.5             
[70] httpuv_1.5.4             cellranger_1.1.0         gtable_0.3.0            
[73] reshape_0.8.8            assertthat_0.2.1         xfun_0.18               
[76] openxlsx_4.2.3           RcppEigen_0.3.3.7.0      later_1.1.0.1           
[79] viridisLite_0.3.0        ellipsis_0.3.1           here_0.1