Last updated: 2020-08-28

Checks: 7 0

Knit directory: Cant_eMLR/

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: 2e6a4ca

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 2e6a4ca. 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:    .Rproj.user/
    Ignored:    data/GLODAPv1_1/
    Ignored:    data/GLODAPv2_2016b_MappedClimatologies/
    Ignored:    data/GLODAPv2_2020/
    Ignored:    data/Gruber_2019/
    Ignored:    data/WOCE/
    Ignored:    data/World_Ocean_Atlas_2013_Clement/
    Ignored:    data/World_Ocean_Atlas_2018/
    Ignored:    data/eMLR/
    Ignored:    data/mapping/
    Ignored:    data/pCO2_atmosphere/
    Ignored:    dump/

Unstaged changes:
    Modified:   output/figure/mapping/GO_new_lon_100.5_model_Cant.png
    Modified:   output/figure/mapping/GO_new_lon_120.5_model_Cant.png
    Modified:   output/figure/mapping/GO_new_lon_140.5_model_Cant.png
    Modified:   output/figure/mapping/GO_new_lon_160.5_model_Cant.png
    Modified:   output/figure/mapping/GO_new_lon_180.5_model_Cant.png
    Modified:   output/figure/mapping/GO_new_lon_20.5_model_Cant.png
    Modified:   output/figure/mapping/GO_new_lon_200.5_model_Cant.png
    Modified:   output/figure/mapping/GO_new_lon_220.5_model_Cant.png
    Modified:   output/figure/mapping/GO_new_lon_240.5_model_Cant.png
    Modified:   output/figure/mapping/GO_new_lon_260.5_model_Cant.png
    Modified:   output/figure/mapping/GO_new_lon_280.5_model_Cant.png
    Modified:   output/figure/mapping/GO_new_lon_300.5_model_Cant.png
    Modified:   output/figure/mapping/GO_new_lon_320.5_model_Cant.png
    Modified:   output/figure/mapping/GO_new_lon_340.5_model_Cant.png
    Modified:   output/figure/mapping/GO_new_lon_40.5_model_Cant.png
    Modified:   output/figure/mapping/GO_new_lon_60.5_model_Cant.png
    Modified:   output/figure/mapping/GO_new_lon_80.5_model_Cant.png
    Modified:   output/figure/mapping/JGOFS_GO_lon_100.5_model_Cant.png
    Modified:   output/figure/mapping/JGOFS_GO_lon_120.5_model_Cant.png
    Modified:   output/figure/mapping/JGOFS_GO_lon_140.5_model_Cant.png
    Modified:   output/figure/mapping/JGOFS_GO_lon_160.5_model_Cant.png
    Modified:   output/figure/mapping/JGOFS_GO_lon_180.5_model_Cant.png
    Modified:   output/figure/mapping/JGOFS_GO_lon_20.5_model_Cant.png
    Modified:   output/figure/mapping/JGOFS_GO_lon_200.5_model_Cant.png
    Modified:   output/figure/mapping/JGOFS_GO_lon_220.5_model_Cant.png
    Modified:   output/figure/mapping/JGOFS_GO_lon_240.5_model_Cant.png
    Modified:   output/figure/mapping/JGOFS_GO_lon_260.5_model_Cant.png
    Modified:   output/figure/mapping/JGOFS_GO_lon_280.5_model_Cant.png
    Modified:   output/figure/mapping/JGOFS_GO_lon_300.5_model_Cant.png
    Modified:   output/figure/mapping/JGOFS_GO_lon_320.5_model_Cant.png
    Modified:   output/figure/mapping/JGOFS_GO_lon_340.5_model_Cant.png
    Modified:   output/figure/mapping/JGOFS_GO_lon_40.5_model_Cant.png
    Modified:   output/figure/mapping/JGOFS_GO_lon_60.5_model_Cant.png
    Modified:   output/figure/mapping/JGOFS_GO_lon_80.5_model_Cant.png

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	2e6a4ca	jens-daniel-mueller	2020-08-28	XXX
html	27404de	jens-daniel-mueller	2020-08-27	Build site.
html	b610bf7	jens-daniel-mueller	2020-08-27	Build site.
Rmd	96ab537	jens-daniel-mueller	2020-08-27	revised
html	b6d0e6a	jens-daniel-mueller	2020-08-27	Build site.
html	f40e48b	jens-daniel-mueller	2020-08-26	Build site.
html	ec20f40	jens-daniel-mueller	2020-08-24	Build site.
Rmd	a804955	jens-daniel-mueller	2020-08-24	split mapping into 2 rmds, po4star selection in parameters, use po4star nitrate
html	5ffe187	jens-daniel-mueller	2020-08-20	Build site.
html	1064ef8	jens-daniel-mueller	2020-08-19	Build site.
Rmd	94f9375	jens-daniel-mueller	2020-08-19	split emlr into data preparation, assumption testing, and model fitting
html	8484d53	jens-daniel-mueller	2020-08-19	Build site.
Rmd	eec0022	jens-daniel-mueller	2020-08-19	split emlr into data preparation, assumption testing, and model fitting

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

1 Required data

Required are:

cleaned and prepared GLODAPv2.2020 file

GLODAP <-
  read_csv(
    here::here(
      "data/GLODAPv2_2020/_summarized_data_files",
      "GLODAP_MLR_fitting_ready.csv"
    )
  )

2 Select phosphate star

Currently, MLR fitting is based on PO4* calculated with:

nitrate

if (parameters$phosphate_star_approach == "nitrate") {
  GLODAP <- GLODAP %>%
    mutate(phosphate_star = phosphate_star_nit) %>%
    select(-c(phosphate_star_oxy, phosphate_star_nit))
} else{
  GLODAP <- GLODAP %>%
    mutate(phosphate_star = phosphate_star_oxy) %>%
    select(-c(phosphate_star_oxy, phosphate_star_nit))
}

3 Model fitting

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.

for (i_basin in unique(GLODAP$basin)) {
  for (i_era in unique(GLODAP$era)) {
    
    # i_basin <- unique(GLODAP$basin)[1]
    # i_era   <- unique(GLODAP$era)[1]
    
    print(i_basin) 
    print(i_era  )
    
    GLODAP_basin_era <- GLODAP %>%
      filter(basin == i_basin, era == i_era)

    for (i_gamma_slab in unique(GLODAP_basin_era$gamma_slab)) {

      # i_gamma_slab <- unique(GLODAP_basin_era$gamma_slab)[1]
      print(i_gamma_slab)
      
      GLODAP_basin_era_slab <- GLODAP_basin_era %>% 
        filter(gamma_slab == i_gamma_slab)
      
      lm_full <- lm(Cstar ~ 
                sal  + 
                tem  + 
                aou  + 
                oxygen  + 
                silicate  + 
                phosphate  + 
                phosphate_star, 
              data = GLODAP_basin_era_slab)
      
      lm_all <- ols_step_all_possible(lm_full)
      
      temp <- as_tibble(lm_all$result)
      
      temp <- temp %>% 
        mutate(basin = i_basin,
               era = i_era,
               gamma_slab = i_gamma_slab)
      
      if (exists("all_lm_stats")) {
        all_lm_stats <- bind_rows(all_lm_stats, temp)
      }
      
      if (!exists("all_lm_stats")) {
        all_lm_stats <- temp
      }
      
      rm(lm_full, lm_all, temp)
      
    }
    
  }
  
}

all_lm_stats %>%
  write_csv(here::here("data/eMLR",
                       "all_lm_stats.csv"))

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

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

all_lm_stats <- read_csv(here::here("data/eMLR",
                                    "all_lm_stats.csv"))

all_lm_stats <- all_lm_stats %>% 
  filter(n >= parameters$predictors_min,
         n <= parameters$predictors_max)

all_lm_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")

3.1 RMSE alternatives

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

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

3.2 Predictor selection

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

all_lm_stats <- all_lm_stats %>%
  select(basin, era, gamma_slab, predictors, rmse)

all_lm_stats_wide <- all_lm_stats %>%
  pivot_wider(names_from = era,
              values_from = rmse,
              names_prefix = "rmse_")

all_lm_stats_wide <- all_lm_stats_wide %>%
  mutate(
    rmse_sum_JGOFS_GO = rmse_JGOFS_WOCE + rmse_GO_SHIP,
    rmse_sum_GO_new = rmse_GO_SHIP + rmse_new_era
  )

all_lm_stats <- all_lm_stats_wide %>%
  pivot_longer(
    c(rmse_sum_JGOFS_GO, rmse_sum_GO_new),
    names_to = "eras",
    values_to = "rmse_sum",
    names_prefix = "rmse_sum_"
  )

lm_best <- all_lm_stats %>%
  group_by(basin, gamma_slab, eras) %>%
  slice_min(order_by = rmse_sum,
            with_ties = FALSE,
            n = parameters$MLR_number) %>%
  ungroup() %>% 
  arrange(basin, gamma_slab, eras, predictors)

kable(lm_best) %>%
  add_header_above() %>%
  kable_styling() %>%
  scroll_box(width = "100%", height = "400px")

basin	gamma_slab	predictors	rmse_JGOFS_WOCE	rmse_GO_SHIP	rmse_new_era	eras	rmse_sum
Atlantic	(-Inf,26]	aou oxygen phosphate phosphate_star	8.223378	6.151569	2.597069	GO_new	8.748638
Atlantic	(-Inf,26]	aou oxygen silicate phosphate phosphate_star	7.551415	6.151535	2.585332	GO_new	8.736868
Atlantic	(-Inf,26]	sal aou oxygen phosphate phosphate_star	7.713991	5.768405	2.539469	GO_new	8.307874
Atlantic	(-Inf,26]	sal tem aou phosphate phosphate_star	7.585762	5.807751	2.536513	GO_new	8.344264
Atlantic	(-Inf,26]	sal tem oxygen phosphate phosphate_star	7.547678	5.837633	2.536054	GO_new	8.373687
Atlantic	(-Inf,26]	tem aou oxygen phosphate phosphate_star	7.150296	6.028173	2.567006	GO_new	8.595179
Atlantic	(-Inf,26]	tem aou phosphate phosphate_star	8.053366	6.103861	2.593184	GO_new	8.697045
Atlantic	(-Inf,26]	tem aou silicate phosphate phosphate_star	7.415798	6.103595	2.582598	GO_new	8.686193
Atlantic	(-Inf,26]	tem oxygen phosphate phosphate_star	7.944696	6.058427	2.586913	GO_new	8.645340
Atlantic	(-Inf,26]	tem oxygen silicate phosphate phosphate_star	7.358462	6.055636	2.578958	GO_new	8.634594
Atlantic	(-Inf,26]	aou oxygen silicate phosphate phosphate_star	7.551415	6.151535	2.585332	JGOFS_GO	13.702951
Atlantic	(-Inf,26]	sal aou oxygen phosphate phosphate_star	7.713991	5.768405	2.539469	JGOFS_GO	13.482396
Atlantic	(-Inf,26]	sal tem aou oxygen phosphate	7.173663	6.639387	2.699587	JGOFS_GO	13.813050
Atlantic	(-Inf,26]	sal tem aou phosphate phosphate_star	7.585762	5.807751	2.536513	JGOFS_GO	13.393513
Atlantic	(-Inf,26]	sal tem oxygen phosphate phosphate_star	7.547678	5.837633	2.536054	JGOFS_GO	13.385311
Atlantic	(-Inf,26]	sal tem oxygen silicate phosphate	7.248967	6.745194	2.784042	JGOFS_GO	13.994160
Atlantic	(-Inf,26]	tem aou oxygen phosphate phosphate_star	7.150296	6.028173	2.567006	JGOFS_GO	13.178469
Atlantic	(-Inf,26]	tem aou silicate phosphate phosphate_star	7.415798	6.103595	2.582598	JGOFS_GO	13.519393
Atlantic	(-Inf,26]	tem oxygen phosphate phosphate_star	7.944696	6.058427	2.586913	JGOFS_GO	14.003123
Atlantic	(-Inf,26]	tem oxygen silicate phosphate phosphate_star	7.358462	6.055636	2.578958	JGOFS_GO	13.414099
Atlantic	(26,26.5]	aou oxygen phosphate phosphate_star	5.479554	7.263998	3.094275	GO_new	10.358274
Atlantic	(26,26.5]	aou oxygen silicate phosphate phosphate_star	5.391757	7.262779	3.059189	GO_new	10.321968
Atlantic	(26,26.5]	sal aou oxygen phosphate phosphate_star	5.456763	7.247045	3.078265	GO_new	10.325310
Atlantic	(26,26.5]	sal tem aou phosphate phosphate_star	5.502307	7.278155	3.085180	GO_new	10.363335
Atlantic	(26,26.5]	sal tem oxygen phosphate phosphate_star	5.558040	7.317006	3.092409	GO_new	10.409414
Atlantic	(26,26.5]	tem aou oxygen phosphate phosphate_star	5.452186	7.236607	3.080579	GO_new	10.317186
Atlantic	(26,26.5]	tem aou oxygen silicate phosphate	5.363554	7.360474	3.065512	GO_new	10.425986
Atlantic	(26,26.5]	tem aou phosphate phosphate_star	5.509470	7.284165	3.103494	GO_new	10.387658
Atlantic	(26,26.5]	tem aou silicate phosphate phosphate_star	5.429777	7.282652	3.070020	GO_new	10.352672
Atlantic	(26,26.5]	tem oxygen silicate phosphate phosphate_star	5.490007	7.317880	3.082355	GO_new	10.400236
Atlantic	(26,26.5]	aou oxygen phosphate phosphate_star	5.479554	7.263998	3.094275	JGOFS_GO	12.743552
Atlantic	(26,26.5]	aou oxygen silicate phosphate phosphate_star	5.391757	7.262779	3.059189	JGOFS_GO	12.654537
Atlantic	(26,26.5]	sal aou oxygen phosphate phosphate_star	5.456763	7.247045	3.078265	JGOFS_GO	12.703808
Atlantic	(26,26.5]	sal tem aou oxygen phosphate	5.477425	7.346736	3.118023	JGOFS_GO	12.824162
Atlantic	(26,26.5]	sal tem aou phosphate phosphate_star	5.502307	7.278155	3.085180	JGOFS_GO	12.780462
Atlantic	(26,26.5]	tem aou oxygen phosphate phosphate_star	5.452186	7.236607	3.080579	JGOFS_GO	12.688793
Atlantic	(26,26.5]	tem aou oxygen silicate phosphate	5.363554	7.360474	3.065512	JGOFS_GO	12.724027
Atlantic	(26,26.5]	tem aou phosphate phosphate_star	5.509470	7.284165	3.103494	JGOFS_GO	12.793635
Atlantic	(26,26.5]	tem aou silicate phosphate phosphate_star	5.429777	7.282652	3.070020	JGOFS_GO	12.712429
Atlantic	(26,26.5]	tem oxygen silicate phosphate phosphate_star	5.490007	7.317880	3.082355	JGOFS_GO	12.807887
Atlantic	(26.5,26.75]	sal aou oxygen silicate phosphate	5.121400	6.278369	3.156785	GO_new	9.435155
Atlantic	(26.5,26.75]	sal tem aou oxygen phosphate	5.318827	6.182233	3.188366	GO_new	9.370599
Atlantic	(26.5,26.75]	sal tem aou phosphate	5.348852	6.294451	3.190597	GO_new	9.485048
Atlantic	(26.5,26.75]	sal tem aou phosphate phosphate_star	5.340770	6.241134	3.189793	GO_new	9.430926
Atlantic	(26.5,26.75]	sal tem aou silicate phosphate	5.080928	6.189493	3.157471	GO_new	9.346963
Atlantic	(26.5,26.75]	sal tem oxygen phosphate	5.332037	6.268202	3.188583	GO_new	9.456786
Atlantic	(26.5,26.75]	sal tem oxygen phosphate phosphate_star	5.324850	6.211578	3.188105	GO_new	9.399684
Atlantic	(26.5,26.75]	sal tem oxygen silicate phosphate	5.073639	6.161326	3.152137	GO_new	9.313462
Atlantic	(26.5,26.75]	tem aou oxygen silicate phosphate	5.153695	6.171505	3.261512	GO_new	9.433017
Atlantic	(26.5,26.75]	tem oxygen silicate phosphate phosphate_star	5.141157	6.276629	3.197287	GO_new	9.473917
Atlantic	(26.5,26.75]	sal aou oxygen silicate phosphate	5.121400	6.278369	3.156785	JGOFS_GO	11.399769
Atlantic	(26.5,26.75]	sal tem aou oxygen phosphate	5.318827	6.182233	3.188366	JGOFS_GO	11.501060
Atlantic	(26.5,26.75]	sal tem aou silicate phosphate	5.080928	6.189493	3.157471	JGOFS_GO	11.270421
Atlantic	(26.5,26.75]	sal tem oxygen phosphate phosphate_star	5.324850	6.211578	3.188105	JGOFS_GO	11.536428
Atlantic	(26.5,26.75]	sal tem oxygen silicate phosphate	5.073639	6.161326	3.152137	JGOFS_GO	11.234964
Atlantic	(26.5,26.75]	tem aou oxygen silicate phosphate	5.153695	6.171505	3.261512	JGOFS_GO	11.325200
Atlantic	(26.5,26.75]	tem aou silicate phosphate	5.173069	6.316120	3.273435	JGOFS_GO	11.489188
Atlantic	(26.5,26.75]	tem aou silicate phosphate phosphate_star	5.154882	6.309885	3.201271	JGOFS_GO	11.464768
Atlantic	(26.5,26.75]	tem oxygen silicate phosphate	5.158997	6.281083	3.267692	JGOFS_GO	11.440080
Atlantic	(26.5,26.75]	tem oxygen silicate phosphate phosphate_star	5.141157	6.276629	3.197287	JGOFS_GO	11.417786
Atlantic	(26.75,27]	aou oxygen silicate phosphate	4.896401	5.405821	3.681006	GO_new	9.086827
Atlantic	(26.75,27]	aou oxygen silicate phosphate phosphate_star	4.835819	5.377290	3.675956	GO_new	9.053246
Atlantic	(26.75,27]	sal aou oxygen silicate phosphate	4.810239	5.344351	3.679790	GO_new	9.024141
Atlantic	(26.75,27]	sal tem aou silicate phosphate	4.708472	5.279945	3.666486	GO_new	8.946431
Atlantic	(26.75,27]	sal tem oxygen silicate phosphate	4.689137	5.255316	3.664798	GO_new	8.920113
Atlantic	(26.75,27]	tem aou oxygen silicate phosphate	4.802227	5.348802	3.665657	GO_new	9.014459
Atlantic	(26.75,27]	tem aou silicate phosphate	4.826200	5.373558	3.667384	GO_new	9.040943
Atlantic	(26.75,27]	tem aou silicate phosphate phosphate_star	4.767595	5.343790	3.662341	GO_new	9.006132
Atlantic	(26.75,27]	tem oxygen silicate phosphate	4.803963	5.356851	3.666006	GO_new	9.022858
Atlantic	(26.75,27]	tem oxygen silicate phosphate phosphate_star	4.739924	5.321954	3.660579	GO_new	8.982533
Atlantic	(26.75,27]	sal aou oxygen silicate phosphate	4.810239	5.344351	3.679790	JGOFS_GO	10.154590
Atlantic	(26.75,27]	sal tem aou oxygen phosphate	4.747838	5.428829	3.705508	JGOFS_GO	10.176667
Atlantic	(26.75,27]	sal tem aou silicate phosphate	4.708472	5.279945	3.666486	JGOFS_GO	9.988417
Atlantic	(26.75,27]	sal tem oxygen phosphate phosphate_star	4.737133	5.414265	3.678545	JGOFS_GO	10.151398
Atlantic	(26.75,27]	sal tem oxygen silicate phosphate	4.689137	5.255316	3.664798	JGOFS_GO	9.944452
Atlantic	(26.75,27]	tem aou oxygen phosphate phosphate_star	4.761392	5.397500	3.708802	JGOFS_GO	10.158892
Atlantic	(26.75,27]	tem aou oxygen silicate phosphate	4.802227	5.348802	3.665657	JGOFS_GO	10.151029
Atlantic	(26.75,27]	tem aou silicate phosphate phosphate_star	4.767595	5.343790	3.662341	JGOFS_GO	10.111385
Atlantic	(26.75,27]	tem oxygen silicate phosphate	4.803963	5.356851	3.666006	JGOFS_GO	10.160814
Atlantic	(26.75,27]	tem oxygen silicate phosphate phosphate_star	4.739924	5.321954	3.660579	JGOFS_GO	10.061879
Atlantic	(27,27.25]	aou oxygen silicate phosphate phosphate_star	5.328621	5.669303	4.933733	GO_new	10.603037
Atlantic	(27,27.25]	sal tem aou phosphate phosphate_star	5.250498	5.696724	5.034078	GO_new	10.730802
Atlantic	(27,27.25]	sal tem aou silicate phosphate	5.078648	5.623468	5.051530	GO_new	10.674998
Atlantic	(27,27.25]	sal tem oxygen phosphate phosphate_star	5.196334	5.658172	5.054511	GO_new	10.712683
Atlantic	(27,27.25]	sal tem oxygen silicate phosphate	5.018688	5.589122	5.077798	GO_new	10.666920
Atlantic	(27,27.25]	tem aou oxygen phosphate phosphate_star	5.147498	5.619345	4.971314	GO_new	10.590659
Atlantic	(27,27.25]	tem aou oxygen silicate phosphate	5.132734	5.610171	5.007417	GO_new	10.617588
Atlantic	(27,27.25]	tem aou silicate phosphate phosphate_star	5.192439	5.582350	4.937753	GO_new	10.520102
Atlantic	(27,27.25]	tem oxygen phosphate phosphate_star	5.234406	5.665572	5.060251	GO_new	10.725823
Atlantic	(27,27.25]	tem oxygen silicate phosphate phosphate_star	5.142108	5.541220	4.952800	GO_new	10.494021
Atlantic	(27,27.25]	sal tem aou silicate phosphate	5.078648	5.623468	5.051530	JGOFS_GO	10.702117
Atlantic	(27,27.25]	sal tem oxygen phosphate phosphate_star	5.196334	5.658172	5.054511	JGOFS_GO	10.854506
Atlantic	(27,27.25]	sal tem oxygen silicate phosphate	5.018688	5.589122	5.077798	JGOFS_GO	10.607810
Atlantic	(27,27.25]	tem aou oxygen phosphate phosphate_star	5.147498	5.619345	4.971314	JGOFS_GO	10.766844
Atlantic	(27,27.25]	tem aou oxygen silicate phosphate	5.132734	5.610171	5.007417	JGOFS_GO	10.742905
Atlantic	(27,27.25]	tem aou silicate phosphate	5.199454	5.645126	5.100418	JGOFS_GO	10.844579
Atlantic	(27,27.25]	tem aou silicate phosphate phosphate_star	5.192439	5.582350	4.937753	JGOFS_GO	10.774788
Atlantic	(27,27.25]	tem oxygen phosphate phosphate_star	5.234406	5.665572	5.060251	JGOFS_GO	10.899978
Atlantic	(27,27.25]	tem oxygen silicate phosphate	5.156878	5.617843	5.133226	JGOFS_GO	10.774722
Atlantic	(27,27.25]	tem oxygen silicate phosphate phosphate_star	5.142108	5.541220	4.952800	JGOFS_GO	10.683329
Atlantic	(27.25,27.5]	aou oxygen phosphate phosphate_star	4.928107	4.547030	3.655045	GO_new	8.202075
Atlantic	(27.25,27.5]	aou oxygen silicate phosphate phosphate_star	4.915736	4.547015	3.648298	GO_new	8.195313
Atlantic	(27.25,27.5]	sal aou oxygen phosphate phosphate_star	4.902251	4.509664	3.633780	GO_new	8.143443
Atlantic	(27.25,27.5]	sal tem aou oxygen phosphate	4.685854	4.506963	3.523918	GO_new	8.030880
Atlantic	(27.25,27.5]	sal tem aou phosphate phosphate_star	4.823760	4.485051	3.663136	GO_new	8.148187
Atlantic	(27.25,27.5]	sal tem oxygen phosphate phosphate_star	4.777170	4.485453	3.700628	GO_new	8.186081
Atlantic	(27.25,27.5]	tem aou oxygen phosphate	4.776702	4.554342	3.528315	GO_new	8.082657
Atlantic	(27.25,27.5]	tem aou oxygen phosphate phosphate_star	4.772737	4.543040	3.421364	GO_new	7.964404
Atlantic	(27.25,27.5]	tem aou oxygen silicate phosphate	4.708352	4.554250	3.467433	GO_new	8.021683
Atlantic	(27.25,27.5]	tem aou silicate phosphate phosphate_star	4.856829	4.544790	3.677480	GO_new	8.222270
Atlantic	(27.25,27.5]	sal tem aou oxygen phosphate	4.685854	4.506963	3.523918	JGOFS_GO	9.192816
Atlantic	(27.25,27.5]	sal tem aou phosphate	4.825263	4.507063	3.906248	JGOFS_GO	9.332326
Atlantic	(27.25,27.5]	sal tem aou phosphate phosphate_star	4.823760	4.485051	3.663136	JGOFS_GO	9.308811
Atlantic	(27.25,27.5]	sal tem aou silicate phosphate	4.824954	4.494151	3.859025	JGOFS_GO	9.319105
Atlantic	(27.25,27.5]	sal tem oxygen phosphate	4.779201	4.514894	3.962830	JGOFS_GO	9.294095
Atlantic	(27.25,27.5]	sal tem oxygen phosphate phosphate_star	4.777170	4.485453	3.700628	JGOFS_GO	9.262623
Atlantic	(27.25,27.5]	sal tem oxygen silicate phosphate	4.779166	4.504624	3.915107	JGOFS_GO	9.283789
Atlantic	(27.25,27.5]	tem aou oxygen phosphate	4.776702	4.554342	3.528315	JGOFS_GO	9.331044
Atlantic	(27.25,27.5]	tem aou oxygen phosphate phosphate_star	4.772737	4.543040	3.421364	JGOFS_GO	9.315776
Atlantic	(27.25,27.5]	tem aou oxygen silicate phosphate	4.708352	4.554250	3.467433	JGOFS_GO	9.262602
Atlantic	(27.5,27.75]	aou oxygen silicate phosphate	4.181349	4.384925	4.678199	GO_new	9.063125
Atlantic	(27.5,27.75]	aou oxygen silicate phosphate phosphate_star	4.165555	4.370717	4.544374	GO_new	8.915091
Atlantic	(27.5,27.75]	sal aou oxygen silicate phosphate	4.170632	4.369707	4.676151	GO_new	9.045858
Atlantic	(27.5,27.75]	sal aou silicate phosphate phosphate_star	4.297026	4.451836	4.564505	GO_new	9.016341
Atlantic	(27.5,27.75]	sal tem aou silicate phosphate	4.157749	4.384594	4.678127	GO_new	9.062721
Atlantic	(27.5,27.75]	tem aou oxygen phosphate phosphate_star	4.221534	4.545312	4.533429	GO_new	9.078741
Atlantic	(27.5,27.75]	tem aou oxygen silicate phosphate	4.148689	4.355735	4.666207	GO_new	9.021941
Atlantic	(27.5,27.75]	tem aou silicate phosphate	4.164381	4.398322	4.681511	GO_new	9.079833
Atlantic	(27.5,27.75]	tem aou silicate phosphate phosphate_star	4.148548	4.383164	4.549982	GO_new	8.933146
Atlantic	(27.5,27.75]	tem oxygen silicate phosphate phosphate_star	4.139702	4.400058	4.556720	GO_new	8.956778
Atlantic	(27.5,27.75]	aou oxygen silicate phosphate	4.181349	4.384925	4.678199	JGOFS_GO	8.566275
Atlantic	(27.5,27.75]	aou oxygen silicate phosphate phosphate_star	4.165555	4.370717	4.544374	JGOFS_GO	8.536272
Atlantic	(27.5,27.75]	sal aou oxygen silicate phosphate	4.170632	4.369707	4.676151	JGOFS_GO	8.540339
Atlantic	(27.5,27.75]	sal tem aou silicate phosphate	4.157749	4.384594	4.678127	JGOFS_GO	8.542343
Atlantic	(27.5,27.75]	sal tem oxygen silicate phosphate	4.148971	4.399401	4.683535	JGOFS_GO	8.548373
Atlantic	(27.5,27.75]	tem aou oxygen silicate phosphate	4.148689	4.355735	4.666207	JGOFS_GO	8.504424
Atlantic	(27.5,27.75]	tem aou silicate phosphate	4.164381	4.398322	4.681511	JGOFS_GO	8.562703
Atlantic	(27.5,27.75]	tem aou silicate phosphate phosphate_star	4.148548	4.383164	4.549982	JGOFS_GO	8.531712
Atlantic	(27.5,27.75]	tem oxygen silicate phosphate	4.155021	4.417619	4.687587	JGOFS_GO	8.572640
Atlantic	(27.5,27.75]	tem oxygen silicate phosphate phosphate_star	4.139702	4.400058	4.556720	JGOFS_GO	8.539760
Atlantic	(27.75,27.85]	aou oxygen silicate phosphate phosphate_star	4.196722	4.307681	7.705766	GO_new	12.013447
Atlantic	(27.75,27.85]	sal aou oxygen silicate phosphate	4.113685	4.285353	7.686739	GO_new	11.972092
Atlantic	(27.75,27.85]	sal aou silicate phosphate	4.120126	4.291980	7.691194	GO_new	11.983174
Atlantic	(27.75,27.85]	sal aou silicate phosphate phosphate_star	4.106682	4.284806	7.682437	GO_new	11.967243
Atlantic	(27.75,27.85]	sal oxygen silicate phosphate	4.121204	4.285377	7.689071	GO_new	11.974448
Atlantic	(27.75,27.85]	sal oxygen silicate phosphate phosphate_star	4.109309	4.277623	7.680367	GO_new	11.957991
Atlantic	(27.75,27.85]	sal tem aou silicate phosphate	4.114571	4.284915	7.687911	GO_new	11.972827
Atlantic	(27.75,27.85]	sal tem oxygen silicate phosphate	4.112834	4.285301	7.687608	GO_new	11.972910
Atlantic	(27.75,27.85]	tem aou silicate phosphate phosphate_star	4.191536	4.306259	7.709321	GO_new	12.015580
Atlantic	(27.75,27.85]	tem oxygen silicate phosphate phosphate_star	4.188746	4.305171	7.708908	GO_new	12.014080
Atlantic	(27.75,27.85]	aou oxygen silicate phosphate phosphate_star	4.196722	4.307681	7.705766	JGOFS_GO	8.504403
Atlantic	(27.75,27.85]	sal aou oxygen silicate phosphate	4.113685	4.285353	7.686739	JGOFS_GO	8.399038
Atlantic	(27.75,27.85]	sal aou silicate phosphate	4.120126	4.291980	7.691194	JGOFS_GO	8.412106
Atlantic	(27.75,27.85]	sal aou silicate phosphate phosphate_star	4.106682	4.284806	7.682437	JGOFS_GO	8.391489
Atlantic	(27.75,27.85]	sal oxygen silicate phosphate	4.121204	4.285377	7.689071	JGOFS_GO	8.406581
Atlantic	(27.75,27.85]	sal oxygen silicate phosphate phosphate_star	4.109309	4.277623	7.680367	JGOFS_GO	8.386932
Atlantic	(27.75,27.85]	sal tem aou silicate phosphate	4.114571	4.284915	7.687911	JGOFS_GO	8.399486
Atlantic	(27.75,27.85]	sal tem oxygen silicate phosphate	4.112834	4.285301	7.687608	JGOFS_GO	8.398135
Atlantic	(27.75,27.85]	tem aou silicate phosphate phosphate_star	4.191536	4.306259	7.709321	JGOFS_GO	8.497795
Atlantic	(27.75,27.85]	tem oxygen silicate phosphate phosphate_star	4.188746	4.305171	7.708908	JGOFS_GO	8.493917
Atlantic	(27.85,27.95]	sal aou oxygen silicate phosphate	3.733203	4.162374	5.126559	GO_new	9.288933
Atlantic	(27.85,27.95]	sal aou silicate phosphate	3.927851	4.640223	5.475140	GO_new	10.115363
Atlantic	(27.85,27.95]	sal aou silicate phosphate phosphate_star	3.895442	4.625502	5.473952	GO_new	10.099454
Atlantic	(27.85,27.95]	sal oxygen silicate phosphate	3.793496	4.463692	5.404331	GO_new	9.868023
Atlantic	(27.85,27.95]	sal oxygen silicate phosphate phosphate_star	3.774106	4.436970	5.397057	GO_new	9.834028
Atlantic	(27.85,27.95]	sal tem aou silicate phosphate	3.755350	4.186637	5.158031	GO_new	9.344668
Atlantic	(27.85,27.95]	sal tem oxygen silicate phosphate	3.744962	4.180818	5.155292	GO_new	9.336111
Atlantic	(27.85,27.95]	sal tem silicate phosphate	3.981432	4.271866	5.181192	GO_new	9.453058
Atlantic	(27.85,27.95]	sal tem silicate phosphate phosphate_star	3.879289	4.271813	5.168314	GO_new	9.440127
Atlantic	(27.85,27.95]	tem aou oxygen silicate phosphate	4.212926	4.831915	5.579193	GO_new	10.411108
Atlantic	(27.85,27.95]	sal aou oxygen silicate phosphate	3.733203	4.162374	5.126559	JGOFS_GO	7.895576
Atlantic	(27.85,27.95]	sal aou silicate phosphate	3.927851	4.640223	5.475140	JGOFS_GO	8.568074
Atlantic	(27.85,27.95]	sal aou silicate phosphate phosphate_star	3.895442	4.625502	5.473952	JGOFS_GO	8.520944
Atlantic	(27.85,27.95]	sal oxygen silicate phosphate	3.793496	4.463692	5.404331	JGOFS_GO	8.257189
Atlantic	(27.85,27.95]	sal oxygen silicate phosphate phosphate_star	3.774106	4.436970	5.397057	JGOFS_GO	8.211076
Atlantic	(27.85,27.95]	sal tem aou silicate phosphate	3.755350	4.186637	5.158031	JGOFS_GO	7.941987
Atlantic	(27.85,27.95]	sal tem oxygen silicate phosphate	3.744962	4.180818	5.155292	JGOFS_GO	7.925780
Atlantic	(27.85,27.95]	sal tem silicate phosphate	3.981432	4.271866	5.181192	JGOFS_GO	8.253298
Atlantic	(27.85,27.95]	sal tem silicate phosphate phosphate_star	3.879289	4.271813	5.168314	JGOFS_GO	8.151102
Atlantic	(27.85,27.95]	tem aou oxygen silicate phosphate	4.212926	4.831915	5.579193	JGOFS_GO	9.044840
Atlantic	(27.95,28.05]	sal aou oxygen silicate phosphate	3.487770	4.308349	5.366509	GO_new	9.674857
Atlantic	(27.95,28.05]	sal aou silicate phosphate	3.757788	4.541320	5.618862	GO_new	10.160182
Atlantic	(27.95,28.05]	sal aou silicate phosphate phosphate_star	3.750648	4.528583	5.554816	GO_new	10.083399
Atlantic	(27.95,28.05]	sal oxygen silicate phosphate	3.570686	4.354920	5.524972	GO_new	9.879893
Atlantic	(27.95,28.05]	sal oxygen silicate phosphate phosphate_star	3.570396	4.330161	5.418256	GO_new	9.748417
Atlantic	(27.95,28.05]	sal tem aou oxygen phosphate	3.755336	4.353971	4.677615	GO_new	9.031586
Atlantic	(27.95,28.05]	sal tem aou silicate phosphate	3.508118	4.354175	5.504209	GO_new	9.858384
Atlantic	(27.95,28.05]	sal tem oxygen silicate phosphate	3.498901	4.330126	5.481394	GO_new	9.811519
Atlantic	(27.95,28.05]	sal tem silicate phosphate phosphate_star	3.755808	4.709496	5.438663	GO_new	10.148158
Atlantic	(27.95,28.05]	tem aou oxygen silicate phosphate	4.117187	4.860501	5.227002	GO_new	10.087504
Atlantic	(27.95,28.05]	sal aou oxygen silicate phosphate	3.487770	4.308349	5.366509	JGOFS_GO	7.796119
Atlantic	(27.95,28.05]	sal aou silicate phosphate	3.757788	4.541320	5.618862	JGOFS_GO	8.299108
Atlantic	(27.95,28.05]	sal aou silicate phosphate phosphate_star	3.750648	4.528583	5.554816	JGOFS_GO	8.279231
Atlantic	(27.95,28.05]	sal oxygen silicate phosphate	3.570686	4.354920	5.524972	JGOFS_GO	7.925607
Atlantic	(27.95,28.05]	sal oxygen silicate phosphate phosphate_star	3.570396	4.330161	5.418256	JGOFS_GO	7.900557
Atlantic	(27.95,28.05]	sal tem aou oxygen phosphate	3.755336	4.353971	4.677615	JGOFS_GO	8.109307
Atlantic	(27.95,28.05]	sal tem aou silicate phosphate	3.508118	4.354175	5.504209	JGOFS_GO	7.862294
Atlantic	(27.95,28.05]	sal tem oxygen phosphate	3.776519	4.617461	5.753398	JGOFS_GO	8.393980
Atlantic	(27.95,28.05]	sal tem oxygen phosphate phosphate_star	3.771038	4.614698	5.706177	JGOFS_GO	8.385736
Atlantic	(27.95,28.05]	sal tem oxygen silicate phosphate	3.498901	4.330126	5.481394	JGOFS_GO	7.829027
Atlantic	(28.05,28.1]	sal tem aou oxygen phosphate	3.539499	4.204534	4.630652	GO_new	8.835185
Atlantic	(28.05,28.1]	sal tem aou phosphate	3.607164	4.569347	5.205477	GO_new	9.774825
Atlantic	(28.05,28.1]	sal tem aou phosphate phosphate_star	3.604413	4.562851	5.185582	GO_new	9.748433
Atlantic	(28.05,28.1]	sal tem aou silicate phosphate	3.586755	4.507157	5.189566	GO_new	9.696724
Atlantic	(28.05,28.1]	sal tem oxygen phosphate	3.557078	4.445505	5.152437	GO_new	9.597943
Atlantic	(28.05,28.1]	sal tem oxygen phosphate phosphate_star	3.554312	4.438866	5.121261	GO_new	9.560127
Atlantic	(28.05,28.1]	sal tem oxygen silicate phosphate	3.531111	4.393917	5.124357	GO_new	9.518275
Atlantic	(28.05,28.1]	tem aou oxygen phosphate	3.737933	4.326767	4.638593	GO_new	8.965360
Atlantic	(28.05,28.1]	tem aou oxygen phosphate phosphate_star	3.737328	4.326389	4.601921	GO_new	8.928310
Atlantic	(28.05,28.1]	tem aou oxygen silicate phosphate	3.643192	4.179970	4.630074	GO_new	8.810044
Atlantic	(28.05,28.1]	sal aou oxygen silicate phosphate	3.554108	4.406193	5.411507	JGOFS_GO	7.960301
Atlantic	(28.05,28.1]	sal oxygen silicate phosphate	3.555467	4.420889	5.419153	JGOFS_GO	7.976356
Atlantic	(28.05,28.1]	sal oxygen silicate phosphate phosphate_star	3.548422	4.401627	5.386608	JGOFS_GO	7.950049
Atlantic	(28.05,28.1]	sal tem aou oxygen phosphate	3.539499	4.204534	4.630652	JGOFS_GO	7.744033
Atlantic	(28.05,28.1]	sal tem oxygen phosphate	3.557078	4.445505	5.152437	JGOFS_GO	8.002583
Atlantic	(28.05,28.1]	sal tem oxygen phosphate phosphate_star	3.554312	4.438866	5.121261	JGOFS_GO	7.993178
Atlantic	(28.05,28.1]	sal tem oxygen silicate phosphate	3.531111	4.393917	5.124357	JGOFS_GO	7.925028
Atlantic	(28.05,28.1]	tem aou oxygen phosphate	3.737933	4.326767	4.638593	JGOFS_GO	8.064700
Atlantic	(28.05,28.1]	tem aou oxygen phosphate phosphate_star	3.737328	4.326389	4.601921	JGOFS_GO	8.063717
Atlantic	(28.05,28.1]	tem aou oxygen silicate phosphate	3.643192	4.179970	4.630074	JGOFS_GO	7.823162
Atlantic	(28.1,28.15]	sal tem aou oxygen phosphate	3.453382	3.943365	4.691597	GO_new	8.634961
Atlantic	(28.1,28.15]	sal tem aou silicate phosphate	3.395913	4.012656	4.782841	GO_new	8.795498
Atlantic	(28.1,28.15]	sal tem oxygen phosphate	3.459589	3.960398	4.753057	GO_new	8.713455
Atlantic	(28.1,28.15]	sal tem oxygen phosphate phosphate_star	3.459522	3.960389	4.749213	GO_new	8.709602
Atlantic	(28.1,28.15]	sal tem oxygen silicate phosphate	3.349985	3.934786	4.738195	GO_new	8.672981
Atlantic	(28.1,28.15]	tem aou oxygen phosphate	3.547872	3.985209	4.716845	GO_new	8.702054
Atlantic	(28.1,28.15]	tem aou oxygen phosphate phosphate_star	3.547716	3.984606	4.715557	GO_new	8.700162
Atlantic	(28.1,28.15]	tem aou oxygen silicate phosphate	3.434468	3.944928	4.715162	GO_new	8.660090
Atlantic	(28.1,28.15]	tem oxygen silicate phosphate	3.434707	3.965040	4.799940	GO_new	8.764981
Atlantic	(28.1,28.15]	tem oxygen silicate phosphate phosphate_star	3.427973	3.964974	4.797528	GO_new	8.762502
Atlantic	(28.1,28.15]	sal tem aou oxygen phosphate	3.453382	3.943365	4.691597	JGOFS_GO	7.396747
Atlantic	(28.1,28.15]	sal tem aou silicate phosphate	3.395913	4.012656	4.782841	JGOFS_GO	7.408570
Atlantic	(28.1,28.15]	sal tem oxygen phosphate	3.459589	3.960398	4.753057	JGOFS_GO	7.419988
Atlantic	(28.1,28.15]	sal tem oxygen phosphate phosphate_star	3.459522	3.960389	4.749213	JGOFS_GO	7.419911
Atlantic	(28.1,28.15]	sal tem oxygen silicate phosphate	3.349985	3.934786	4.738195	JGOFS_GO	7.284771
Atlantic	(28.1,28.15]	tem aou oxygen silicate phosphate	3.434468	3.944928	4.715162	JGOFS_GO	7.379396
Atlantic	(28.1,28.15]	tem aou silicate phosphate	3.459579	4.030995	4.865916	JGOFS_GO	7.490573
Atlantic	(28.1,28.15]	tem aou silicate phosphate phosphate_star	3.454261	4.030922	4.856988	JGOFS_GO	7.485182
Atlantic	(28.1,28.15]	tem oxygen silicate phosphate	3.434707	3.965040	4.799940	JGOFS_GO	7.399747
Atlantic	(28.1,28.15]	tem oxygen silicate phosphate phosphate_star	3.427973	3.964974	4.797528	JGOFS_GO	7.392947
Atlantic	(28.15,28.2]	sal aou oxygen silicate phosphate	3.726384	3.637891	3.023471	GO_new	6.661362
Atlantic	(28.15,28.2]	sal tem aou oxygen phosphate	3.726099	3.629154	3.004142	GO_new	6.633295
Atlantic	(28.15,28.2]	sal tem aou phosphate	3.741007	3.656291	3.004525	GO_new	6.660816
Atlantic	(28.15,28.2]	sal tem aou phosphate phosphate_star	3.732000	3.655905	2.997692	GO_new	6.653596
Atlantic	(28.15,28.2]	sal tem aou silicate phosphate	3.739605	3.638041	2.989969	GO_new	6.628010
Atlantic	(28.15,28.2]	sal tem oxygen phosphate	3.781131	3.701975	3.004907	GO_new	6.706882
Atlantic	(28.15,28.2]	sal tem oxygen phosphate phosphate_star	3.775184	3.701111	2.996741	GO_new	6.697852
Atlantic	(28.15,28.2]	sal tem oxygen silicate phosphate	3.770684	3.671003	2.995966	GO_new	6.666968
Atlantic	(28.15,28.2]	tem aou oxygen phosphate phosphate_star	3.738522	3.641657	3.089357	GO_new	6.731013
Atlantic	(28.15,28.2]	tem aou oxygen silicate phosphate	3.748534	3.636938	3.094567	GO_new	6.731506
Atlantic	(28.15,28.2]	sal aou oxygen phosphate	3.727404	3.639684	3.100446	JGOFS_GO	7.367088
Atlantic	(28.15,28.2]	sal aou oxygen phosphate phosphate_star	3.719971	3.639684	3.099907	JGOFS_GO	7.359655
Atlantic	(28.15,28.2]	sal aou oxygen silicate phosphate	3.726384	3.637891	3.023471	JGOFS_GO	7.364276
Atlantic	(28.15,28.2]	sal tem aou oxygen phosphate	3.726099	3.629154	3.004142	JGOFS_GO	7.355253
Atlantic	(28.15,28.2]	sal tem aou phosphate	3.741007	3.656291	3.004525	JGOFS_GO	7.397297
Atlantic	(28.15,28.2]	sal tem aou phosphate phosphate_star	3.732000	3.655905	2.997692	JGOFS_GO	7.387904
Atlantic	(28.15,28.2]	sal tem aou silicate phosphate	3.739605	3.638041	2.989969	JGOFS_GO	7.377646
Atlantic	(28.15,28.2]	tem aou oxygen phosphate	3.749523	3.641831	3.097179	JGOFS_GO	7.391355
Atlantic	(28.15,28.2]	tem aou oxygen phosphate phosphate_star	3.738522	3.641657	3.089357	JGOFS_GO	7.380178
Atlantic	(28.15,28.2]	tem aou oxygen silicate phosphate	3.748534	3.636938	3.094567	JGOFS_GO	7.385473
Atlantic	(28.2, Inf]	sal aou oxygen phosphate	3.549135	3.174916	3.419848	GO_new	6.594764
Atlantic	(28.2, Inf]	sal aou oxygen phosphate phosphate_star	3.546339	3.160649	3.411306	GO_new	6.571955
Atlantic	(28.2, Inf]	sal aou oxygen silicate phosphate	3.487877	3.158126	3.393186	GO_new	6.551312
Atlantic	(28.2, Inf]	sal oxygen phosphate phosphate_star	3.978629	3.172926	3.414075	GO_new	6.587001
Atlantic	(28.2, Inf]	sal oxygen silicate phosphate	3.765786	3.165659	3.395547	GO_new	6.561207
Atlantic	(28.2, Inf]	sal oxygen silicate phosphate phosphate_star	3.756399	3.137716	3.392982	GO_new	6.530697
Atlantic	(28.2, Inf]	sal tem aou oxygen phosphate	3.522234	3.173511	3.190100	GO_new	6.363611
Atlantic	(28.2, Inf]	sal tem oxygen phosphate	3.651581	3.173715	3.346153	GO_new	6.519868
Atlantic	(28.2, Inf]	sal tem oxygen phosphate phosphate_star	3.649876	3.162480	3.345556	GO_new	6.508036
Atlantic	(28.2, Inf]	sal tem oxygen silicate phosphate	3.544433	3.152200	3.338683	GO_new	6.490883
Atlantic	(28.2, Inf]	sal aou oxygen phosphate	3.549135	3.174916	3.419848	JGOFS_GO	6.724051
Atlantic	(28.2, Inf]	sal aou oxygen phosphate phosphate_star	3.546339	3.160649	3.411306	JGOFS_GO	6.706988
Atlantic	(28.2, Inf]	sal aou oxygen silicate phosphate	3.487877	3.158126	3.393186	JGOFS_GO	6.646004
Atlantic	(28.2, Inf]	sal aou silicate phosphate	3.530184	3.237823	3.463376	JGOFS_GO	6.768007
Atlantic	(28.2, Inf]	sal aou silicate phosphate phosphate_star	3.518405	3.220199	3.456579	JGOFS_GO	6.738604
Atlantic	(28.2, Inf]	sal tem aou oxygen phosphate	3.522234	3.173511	3.190100	JGOFS_GO	6.695745
Atlantic	(28.2, Inf]	sal tem aou phosphate	3.568388	3.192194	3.467972	JGOFS_GO	6.760581
Atlantic	(28.2, Inf]	sal tem aou phosphate phosphate_star	3.565166	3.171852	3.465054	JGOFS_GO	6.737018
Atlantic	(28.2, Inf]	sal tem aou silicate phosphate	3.499887	3.184788	3.462424	JGOFS_GO	6.684675
Atlantic	(28.2, Inf]	sal tem oxygen silicate phosphate	3.544433	3.152200	3.338683	JGOFS_GO	6.696633
Indo-Pacific	(-Inf,26]	sal aou oxygen phosphate phosphate_star	11.656608	6.900704	4.748256	GO_new	11.648960
Indo-Pacific	(-Inf,26]	sal aou oxygen silicate phosphate	11.589172	7.158929	5.649148	GO_new	12.808077
Indo-Pacific	(-Inf,26]	sal tem aou phosphate phosphate_star	11.727854	6.878761	4.968427	GO_new	11.847189
Indo-Pacific	(-Inf,26]	sal tem aou silicate phosphate	11.658578	7.235892	5.852129	GO_new	13.088021
Indo-Pacific	(-Inf,26]	sal tem oxygen phosphate phosphate_star	11.772722	6.912573	5.114563	GO_new	12.027136
Indo-Pacific	(-Inf,26]	sal tem oxygen silicate phosphate	11.696367	7.284738	5.952533	GO_new	13.237271
Indo-Pacific	(-Inf,26]	tem aou oxygen phosphate phosphate_star	12.336939	7.286502	5.639494	GO_new	12.925996
Indo-Pacific	(-Inf,26]	tem aou silicate phosphate phosphate_star	12.281655	7.464928	5.655524	GO_new	13.120452
Indo-Pacific	(-Inf,26]	tem oxygen phosphate phosphate_star	12.406986	7.414955	5.643696	GO_new	13.058651
Indo-Pacific	(-Inf,26]	tem oxygen silicate phosphate phosphate_star	12.209892	7.343438	5.628143	GO_new	12.971580
Indo-Pacific	(-Inf,26]	sal aou oxygen phosphate	11.847558	7.595246	6.038692	JGOFS_GO	19.442804
Indo-Pacific	(-Inf,26]	sal aou oxygen phosphate phosphate_star	11.656608	6.900704	4.748256	JGOFS_GO	18.557312
Indo-Pacific	(-Inf,26]	sal aou oxygen silicate phosphate	11.589172	7.158929	5.649148	JGOFS_GO	18.748101
Indo-Pacific	(-Inf,26]	sal tem aou oxygen phosphate	11.835399	7.592149	6.001991	JGOFS_GO	19.427548
Indo-Pacific	(-Inf,26]	sal tem aou phosphate	11.936886	7.625869	6.165843	JGOFS_GO	19.562755
Indo-Pacific	(-Inf,26]	sal tem aou phosphate phosphate_star	11.727854	6.878761	4.968427	JGOFS_GO	18.606615
Indo-Pacific	(-Inf,26]	sal tem aou silicate phosphate	11.658578	7.235892	5.852129	JGOFS_GO	18.894470
Indo-Pacific	(-Inf,26]	sal tem oxygen phosphate phosphate_star	11.772722	6.912573	5.114563	JGOFS_GO	18.685295
Indo-Pacific	(-Inf,26]	sal tem oxygen silicate phosphate	11.696367	7.284738	5.952533	JGOFS_GO	18.981104
Indo-Pacific	(-Inf,26]	tem oxygen silicate phosphate phosphate_star	12.209892	7.343438	5.628143	JGOFS_GO	19.553330
Indo-Pacific	(26,26.5]	sal aou oxygen phosphate	7.963659	5.384444	4.264629	GO_new	9.649072
Indo-Pacific	(26,26.5]	sal aou oxygen phosphate phosphate_star	7.817898	5.209679	3.993262	GO_new	9.202941
Indo-Pacific	(26,26.5]	sal aou oxygen silicate phosphate	7.855393	5.185799	3.976710	GO_new	9.162509
Indo-Pacific	(26,26.5]	sal tem aou oxygen phosphate	7.962543	5.283401	4.258474	GO_new	9.541875
Indo-Pacific	(26,26.5]	sal tem aou phosphate	8.017736	5.295081	4.301171	GO_new	9.596252
Indo-Pacific	(26,26.5]	sal tem aou phosphate phosphate_star	7.880514	5.108916	4.035301	GO_new	9.144218
Indo-Pacific	(26,26.5]	sal tem aou silicate phosphate	7.929705	5.143958	4.051154	GO_new	9.195112
Indo-Pacific	(26,26.5]	sal tem oxygen phosphate	8.048135	5.284498	4.321445	GO_new	9.605943
Indo-Pacific	(26,26.5]	sal tem oxygen phosphate phosphate_star	7.921449	5.101898	4.066271	GO_new	9.168169
Indo-Pacific	(26,26.5]	sal tem oxygen silicate phosphate	7.966619	5.153265	4.088387	GO_new	9.241652
Indo-Pacific	(26,26.5]	sal aou oxygen phosphate	7.963659	5.384444	4.264629	JGOFS_GO	13.348102
Indo-Pacific	(26,26.5]	sal aou oxygen phosphate phosphate_star	7.817898	5.209679	3.993262	JGOFS_GO	13.027577
Indo-Pacific	(26,26.5]	sal aou oxygen silicate phosphate	7.855393	5.185799	3.976710	JGOFS_GO	13.041193
Indo-Pacific	(26,26.5]	sal tem aou oxygen phosphate	7.962543	5.283401	4.258474	JGOFS_GO	13.245944
Indo-Pacific	(26,26.5]	sal tem aou phosphate	8.017736	5.295081	4.301171	JGOFS_GO	13.312817
Indo-Pacific	(26,26.5]	sal tem aou phosphate phosphate_star	7.880514	5.108916	4.035301	JGOFS_GO	12.989430
Indo-Pacific	(26,26.5]	sal tem aou silicate phosphate	7.929705	5.143958	4.051154	JGOFS_GO	13.073663
Indo-Pacific	(26,26.5]	sal tem oxygen phosphate	8.048135	5.284498	4.321445	JGOFS_GO	13.332633
Indo-Pacific	(26,26.5]	sal tem oxygen phosphate phosphate_star	7.921449	5.101898	4.066271	JGOFS_GO	13.023347
Indo-Pacific	(26,26.5]	sal tem oxygen silicate phosphate	7.966619	5.153265	4.088387	JGOFS_GO	13.119885
Indo-Pacific	(26.5,26.75]	aou oxygen silicate phosphate phosphate_star	6.122193	5.842826	4.471119	GO_new	10.313944
Indo-Pacific	(26.5,26.75]	sal aou oxygen phosphate	5.972964	5.696797	4.631952	GO_new	10.328749
Indo-Pacific	(26.5,26.75]	sal aou oxygen phosphate phosphate_star	5.629568	5.431213	4.210965	GO_new	9.642178
Indo-Pacific	(26.5,26.75]	sal aou oxygen silicate phosphate	5.889764	5.691859	4.563197	GO_new	10.255056
Indo-Pacific	(26.5,26.75]	sal tem aou oxygen phosphate	5.970854	5.696790	4.615782	GO_new	10.312572
Indo-Pacific	(26.5,26.75]	sal tem aou phosphate phosphate_star	5.730080	5.471410	4.292707	GO_new	9.764117
Indo-Pacific	(26.5,26.75]	sal tem aou silicate phosphate	6.032011	5.709349	4.659652	GO_new	10.369001
Indo-Pacific	(26.5,26.75]	sal tem oxygen phosphate phosphate_star	5.794994	5.512058	4.353305	GO_new	9.865362
Indo-Pacific	(26.5,26.75]	tem aou oxygen phosphate phosphate_star	6.099615	5.841145	4.577322	GO_new	10.418466
Indo-Pacific	(26.5,26.75]	tem aou silicate phosphate phosphate_star	6.164107	5.854419	4.582028	GO_new	10.436447
Indo-Pacific	(26.5,26.75]	sal aou oxygen phosphate	5.972964	5.696797	4.631952	JGOFS_GO	11.669761
Indo-Pacific	(26.5,26.75]	sal aou oxygen phosphate phosphate_star	5.629568	5.431213	4.210965	JGOFS_GO	11.060781
Indo-Pacific	(26.5,26.75]	sal aou oxygen silicate phosphate	5.889764	5.691859	4.563197	JGOFS_GO	11.581623
Indo-Pacific	(26.5,26.75]	sal tem aou oxygen phosphate	5.970854	5.696790	4.615782	JGOFS_GO	11.667644
Indo-Pacific	(26.5,26.75]	sal tem aou phosphate	6.059163	5.732700	4.712027	JGOFS_GO	11.791863
Indo-Pacific	(26.5,26.75]	sal tem aou phosphate phosphate_star	5.730080	5.471410	4.292707	JGOFS_GO	11.201489
Indo-Pacific	(26.5,26.75]	sal tem aou silicate phosphate	6.032011	5.709349	4.659652	JGOFS_GO	11.741360
Indo-Pacific	(26.5,26.75]	sal tem oxygen phosphate	6.103238	5.761208	4.756842	JGOFS_GO	11.864446
Indo-Pacific	(26.5,26.75]	sal tem oxygen phosphate phosphate_star	5.794994	5.512058	4.353305	JGOFS_GO	11.307052
Indo-Pacific	(26.5,26.75]	sal tem oxygen silicate phosphate	6.085625	5.726657	4.715326	JGOFS_GO	11.812282
Indo-Pacific	(26.75,27]	aou oxygen phosphate phosphate_star	6.567840	5.032576	4.161568	GO_new	9.194144
Indo-Pacific	(26.75,27]	aou oxygen silicate phosphate phosphate_star	6.488004	4.974323	4.160697	GO_new	9.135019
Indo-Pacific	(26.75,27]	sal aou oxygen phosphate phosphate_star	6.533269	4.972238	4.063475	GO_new	9.035714
Indo-Pacific	(26.75,27]	sal aou oxygen silicate phosphate	6.552579	5.062787	4.421716	GO_new	9.484503
Indo-Pacific	(26.75,27]	sal tem aou oxygen phosphate	6.516743	5.021419	4.250318	GO_new	9.271737
Indo-Pacific	(26.75,27]	sal tem aou phosphate phosphate_star	6.645887	5.160774	4.264093	GO_new	9.424867
Indo-Pacific	(26.75,27]	tem aou oxygen phosphate	6.520101	5.023021	4.264277	GO_new	9.287298
Indo-Pacific	(26.75,27]	tem aou oxygen phosphate phosphate_star	6.482046	4.918858	4.034617	GO_new	8.953476
Indo-Pacific	(26.75,27]	tem aou oxygen silicate phosphate	6.492756	4.983055	4.235837	GO_new	9.218892
Indo-Pacific	(26.75,27]	tem aou silicate phosphate phosphate_star	6.519436	5.038575	4.289332	GO_new	9.327907
Indo-Pacific	(26.75,27]	aou oxygen phosphate phosphate_star	6.567840	5.032576	4.161568	JGOFS_GO	11.600415
Indo-Pacific	(26.75,27]	aou oxygen silicate phosphate phosphate_star	6.488004	4.974323	4.160697	JGOFS_GO	11.462326
Indo-Pacific	(26.75,27]	sal aou oxygen phosphate phosphate_star	6.533269	4.972238	4.063475	JGOFS_GO	11.505507
Indo-Pacific	(26.75,27]	sal aou oxygen silicate phosphate	6.552579	5.062787	4.421716	JGOFS_GO	11.615365
Indo-Pacific	(26.75,27]	sal tem aou oxygen phosphate	6.516743	5.021419	4.250318	JGOFS_GO	11.538162
Indo-Pacific	(26.75,27]	tem aou oxygen phosphate	6.520101	5.023021	4.264277	JGOFS_GO	11.543122
Indo-Pacific	(26.75,27]	tem aou oxygen phosphate phosphate_star	6.482046	4.918858	4.034617	JGOFS_GO	11.400905
Indo-Pacific	(26.75,27]	tem aou oxygen silicate phosphate	6.492756	4.983055	4.235837	JGOFS_GO	11.475811
Indo-Pacific	(26.75,27]	tem aou silicate phosphate phosphate_star	6.519436	5.038575	4.289332	JGOFS_GO	11.558011
Indo-Pacific	(26.75,27]	tem oxygen silicate phosphate phosphate_star	6.555616	5.102646	4.383563	JGOFS_GO	11.658262
Indo-Pacific	(27,27.25]	sal aou oxygen phosphate phosphate_star	6.376268	4.789831	4.178131	GO_new	8.967962
Indo-Pacific	(27,27.25]	sal aou oxygen silicate phosphate	6.377139	4.610352	4.154787	GO_new	8.765140
Indo-Pacific	(27,27.25]	sal aou silicate phosphate	6.502600	4.621347	4.254518	GO_new	8.875865
Indo-Pacific	(27,27.25]	sal aou silicate phosphate phosphate_star	6.447889	4.618712	4.254365	GO_new	8.873077
Indo-Pacific	(27,27.25]	sal tem aou oxygen phosphate	5.996869	4.633139	4.167623	GO_new	8.800762
Indo-Pacific	(27,27.25]	sal tem aou silicate phosphate	6.492264	4.620740	4.179057	GO_new	8.799797
Indo-Pacific	(27,27.25]	sal tem oxygen silicate phosphate	6.598846	4.662267	4.211183	GO_new	8.873450
Indo-Pacific	(27,27.25]	tem aou oxygen phosphate	6.024170	4.748844	4.304139	GO_new	9.052983
Indo-Pacific	(27,27.25]	tem aou oxygen phosphate phosphate_star	6.022258	4.655002	4.175739	GO_new	8.830741
Indo-Pacific	(27,27.25]	tem aou oxygen silicate phosphate	5.981049	4.537003	4.124776	GO_new	8.661779
Indo-Pacific	(27,27.25]	aou oxygen silicate phosphate phosphate_star	6.295462	4.772678	4.374179	JGOFS_GO	11.068141
Indo-Pacific	(27,27.25]	sal aou oxygen phosphate phosphate_star	6.376268	4.789831	4.178131	JGOFS_GO	11.166100
Indo-Pacific	(27,27.25]	sal aou oxygen silicate phosphate	6.377139	4.610352	4.154787	JGOFS_GO	10.987491
Indo-Pacific	(27,27.25]	sal aou silicate phosphate	6.502600	4.621347	4.254518	JGOFS_GO	11.123948
Indo-Pacific	(27,27.25]	sal aou silicate phosphate phosphate_star	6.447889	4.618712	4.254365	JGOFS_GO	11.066601
Indo-Pacific	(27,27.25]	sal tem aou oxygen phosphate	5.996869	4.633139	4.167623	JGOFS_GO	10.630008
Indo-Pacific	(27,27.25]	sal tem aou silicate phosphate	6.492264	4.620740	4.179057	JGOFS_GO	11.113004
Indo-Pacific	(27,27.25]	tem aou oxygen phosphate	6.024170	4.748844	4.304139	JGOFS_GO	10.773014
Indo-Pacific	(27,27.25]	tem aou oxygen phosphate phosphate_star	6.022258	4.655002	4.175739	JGOFS_GO	10.677260
Indo-Pacific	(27,27.25]	tem aou oxygen silicate phosphate	5.981049	4.537003	4.124776	JGOFS_GO	10.518052
Indo-Pacific	(27.25,27.5]	aou oxygen silicate phosphate phosphate_star	5.196364	4.240548	3.681856	GO_new	7.922403
Indo-Pacific	(27.25,27.5]	sal aou oxygen phosphate phosphate_star	5.168633	4.254234	3.572729	GO_new	7.826964
Indo-Pacific	(27.25,27.5]	sal aou oxygen silicate phosphate	5.216429	4.136611	3.512906	GO_new	7.649517
Indo-Pacific	(27.25,27.5]	sal tem aou oxygen phosphate	4.831746	4.162230	3.629966	GO_new	7.792197
Indo-Pacific	(27.25,27.5]	sal tem aou phosphate phosphate_star	5.509279	4.347958	3.586953	GO_new	7.934910
Indo-Pacific	(27.25,27.5]	sal tem aou silicate phosphate	5.541110	4.185755	3.535978	GO_new	7.721733
Indo-Pacific	(27.25,27.5]	sal tem oxygen silicate phosphate	5.764252	4.233532	3.556430	GO_new	7.789962
Indo-Pacific	(27.25,27.5]	tem aou oxygen phosphate	4.835593	4.192637	3.734596	GO_new	7.927233
Indo-Pacific	(27.25,27.5]	tem aou oxygen phosphate phosphate_star	4.832492	4.093420	3.684443	GO_new	7.777862
Indo-Pacific	(27.25,27.5]	tem aou oxygen silicate phosphate	4.804957	3.986913	3.603040	GO_new	7.589953
Indo-Pacific	(27.25,27.5]	aou oxygen phosphate phosphate_star	5.199698	4.327634	3.762548	JGOFS_GO	9.527332
Indo-Pacific	(27.25,27.5]	aou oxygen silicate phosphate	5.256328	4.307021	3.757882	JGOFS_GO	9.563349
Indo-Pacific	(27.25,27.5]	aou oxygen silicate phosphate phosphate_star	5.196364	4.240548	3.681856	JGOFS_GO	9.436912
Indo-Pacific	(27.25,27.5]	sal aou oxygen phosphate	5.224552	4.399037	3.638993	JGOFS_GO	9.623589
Indo-Pacific	(27.25,27.5]	sal aou oxygen phosphate phosphate_star	5.168633	4.254234	3.572729	JGOFS_GO	9.422867
Indo-Pacific	(27.25,27.5]	sal aou oxygen silicate phosphate	5.216429	4.136611	3.512906	JGOFS_GO	9.353040
Indo-Pacific	(27.25,27.5]	sal tem aou oxygen phosphate	4.831746	4.162230	3.629966	JGOFS_GO	8.993977
Indo-Pacific	(27.25,27.5]	tem aou oxygen phosphate	4.835593	4.192637	3.734596	JGOFS_GO	9.028231
Indo-Pacific	(27.25,27.5]	tem aou oxygen phosphate phosphate_star	4.832492	4.093420	3.684443	JGOFS_GO	8.925912
Indo-Pacific	(27.25,27.5]	tem aou oxygen silicate phosphate	4.804957	3.986913	3.603040	JGOFS_GO	8.791870
Indo-Pacific	(27.5,27.75]	sal aou oxygen phosphate phosphate_star	4.656269	3.907055	3.757383	GO_new	7.664438
Indo-Pacific	(27.5,27.75]	sal aou oxygen silicate phosphate	4.655491	3.997240	3.735975	GO_new	7.733215
Indo-Pacific	(27.5,27.75]	sal tem aou oxygen phosphate	4.409996	4.005957	3.525039	GO_new	7.530996
Indo-Pacific	(27.5,27.75]	sal tem aou phosphate	4.826836	4.047735	3.639034	GO_new	7.686769
Indo-Pacific	(27.5,27.75]	sal tem aou phosphate phosphate_star	4.826586	3.894715	3.637182	GO_new	7.531897
Indo-Pacific	(27.5,27.75]	sal tem aou silicate phosphate	4.822537	4.035117	3.627796	GO_new	7.662913
Indo-Pacific	(27.5,27.75]	sal tem oxygen phosphate	4.957364	4.095366	3.574533	GO_new	7.669899
Indo-Pacific	(27.5,27.75]	sal tem oxygen phosphate phosphate_star	4.957076	3.900138	3.561704	GO_new	7.461841
Indo-Pacific	(27.5,27.75]	sal tem oxygen silicate phosphate	4.949934	4.083693	3.565430	GO_new	7.649123
Indo-Pacific	(27.5,27.75]	tem aou oxygen phosphate phosphate_star	4.409247	3.996479	3.684163	GO_new	7.680643
Indo-Pacific	(27.5,27.75]	aou oxygen phosphate phosphate_star	4.663346	3.997608	3.906815	JGOFS_GO	8.660953
Indo-Pacific	(27.5,27.75]	aou oxygen silicate phosphate phosphate_star	4.660946	3.997539	3.901165	JGOFS_GO	8.658485
Indo-Pacific	(27.5,27.75]	sal aou oxygen phosphate	4.656383	4.011924	3.758951	JGOFS_GO	8.668307
Indo-Pacific	(27.5,27.75]	sal aou oxygen phosphate phosphate_star	4.656269	3.907055	3.757383	JGOFS_GO	8.563324
Indo-Pacific	(27.5,27.75]	sal aou oxygen silicate phosphate	4.655491	3.997240	3.735975	JGOFS_GO	8.652731
Indo-Pacific	(27.5,27.75]	sal tem aou oxygen phosphate	4.409996	4.005957	3.525039	JGOFS_GO	8.415953
Indo-Pacific	(27.5,27.75]	sal tem aou phosphate phosphate_star	4.826586	3.894715	3.637182	JGOFS_GO	8.721300
Indo-Pacific	(27.5,27.75]	tem aou oxygen phosphate	4.410507	4.074666	3.760424	JGOFS_GO	8.485173
Indo-Pacific	(27.5,27.75]	tem aou oxygen phosphate phosphate_star	4.409247	3.996479	3.684163	JGOFS_GO	8.405726
Indo-Pacific	(27.5,27.75]	tem aou oxygen silicate phosphate	4.386458	4.074656	3.738199	JGOFS_GO	8.461113
Indo-Pacific	(27.75,27.85]	aou oxygen silicate phosphate	4.249542	3.356084	2.945707	GO_new	6.301792
Indo-Pacific	(27.75,27.85]	aou oxygen silicate phosphate phosphate_star	4.245475	3.332871	2.938537	GO_new	6.271408
Indo-Pacific	(27.75,27.85]	sal aou oxygen silicate phosphate	4.235175	3.347718	2.896917	GO_new	6.244636
Indo-Pacific	(27.75,27.85]	sal tem aou silicate phosphate	4.258533	3.311684	2.876268	GO_new	6.187952
Indo-Pacific	(27.75,27.85]	sal tem oxygen silicate phosphate	4.281560	3.291064	2.865015	GO_new	6.156079
Indo-Pacific	(27.75,27.85]	tem aou oxygen silicate phosphate	4.242882	3.271777	2.903747	GO_new	6.175523
Indo-Pacific	(27.75,27.85]	tem aou silicate phosphate	4.275388	3.316085	2.924736	GO_new	6.240820
Indo-Pacific	(27.75,27.85]	tem aou silicate phosphate phosphate_star	4.273516	3.302046	2.918998	GO_new	6.221044
Indo-Pacific	(27.75,27.85]	tem oxygen silicate phosphate	4.300570	3.296323	2.916579	GO_new	6.212903
Indo-Pacific	(27.75,27.85]	tem oxygen silicate phosphate phosphate_star	4.299485	3.286256	2.911216	GO_new	6.197472
Indo-Pacific	(27.75,27.85]	aou oxygen silicate phosphate phosphate_star	4.245475	3.332871	2.938537	JGOFS_GO	7.578346
Indo-Pacific	(27.75,27.85]	sal aou oxygen silicate phosphate	4.235175	3.347718	2.896917	JGOFS_GO	7.582893
Indo-Pacific	(27.75,27.85]	sal tem aou oxygen phosphate	4.233667	3.327161	3.001071	JGOFS_GO	7.560827
Indo-Pacific	(27.75,27.85]	sal tem aou silicate phosphate	4.258533	3.311684	2.876268	JGOFS_GO	7.570217
Indo-Pacific	(27.75,27.85]	sal tem oxygen silicate phosphate	4.281560	3.291064	2.865015	JGOFS_GO	7.572624
Indo-Pacific	(27.75,27.85]	tem aou oxygen phosphate	4.242883	3.327650	3.002873	JGOFS_GO	7.570533
Indo-Pacific	(27.75,27.85]	tem aou oxygen phosphate phosphate_star	4.234736	3.327543	3.000201	JGOFS_GO	7.562279
Indo-Pacific	(27.75,27.85]	tem aou oxygen silicate phosphate	4.242882	3.271777	2.903747	JGOFS_GO	7.514659
Indo-Pacific	(27.75,27.85]	tem aou silicate phosphate phosphate_star	4.273516	3.302046	2.918998	JGOFS_GO	7.575562
Indo-Pacific	(27.75,27.85]	tem oxygen silicate phosphate phosphate_star	4.299485	3.286256	2.911216	JGOFS_GO	7.585741
Indo-Pacific	(27.85,27.95]	aou oxygen silicate phosphate	4.280926	3.007769	2.912915	GO_new	5.920684
Indo-Pacific	(27.85,27.95]	aou oxygen silicate phosphate phosphate_star	4.256269	2.993750	2.902527	GO_new	5.896277
Indo-Pacific	(27.85,27.95]	sal aou oxygen silicate phosphate	4.272593	3.007732	2.896370	GO_new	5.904102
Indo-Pacific	(27.85,27.95]	sal tem aou silicate phosphate	4.290311	3.000580	2.878296	GO_new	5.878876
Indo-Pacific	(27.85,27.95]	sal tem oxygen silicate phosphate	4.311405	3.004330	2.869502	GO_new	5.873832
Indo-Pacific	(27.85,27.95]	tem aou oxygen silicate phosphate	4.280774	3.000763	2.891621	GO_new	5.892384
Indo-Pacific	(27.85,27.95]	tem aou silicate phosphate	4.295926	3.000801	2.897066	GO_new	5.897867
Indo-Pacific	(27.85,27.95]	tem aou silicate phosphate phosphate_star	4.274004	2.989254	2.887257	GO_new	5.876511
Indo-Pacific	(27.85,27.95]	tem oxygen silicate phosphate	4.315597	3.005155	2.892410	GO_new	5.897565
Indo-Pacific	(27.85,27.95]	tem oxygen silicate phosphate phosphate_star	4.294333	2.994511	2.882504	GO_new	5.877015
Indo-Pacific	(27.85,27.95]	aou oxygen phosphate phosphate_star	4.256599	3.033388	2.960681	JGOFS_GO	7.289987
Indo-Pacific	(27.85,27.95]	aou oxygen silicate phosphate	4.280926	3.007769	2.912915	JGOFS_GO	7.288694
Indo-Pacific	(27.85,27.95]	aou oxygen silicate phosphate phosphate_star	4.256269	2.993750	2.902527	JGOFS_GO	7.250018
Indo-Pacific	(27.85,27.95]	sal aou oxygen phosphate phosphate_star	4.237357	3.016584	2.959140	JGOFS_GO	7.253941
Indo-Pacific	(27.85,27.95]	sal aou oxygen silicate phosphate	4.272593	3.007732	2.896370	JGOFS_GO	7.280326
Indo-Pacific	(27.85,27.95]	sal tem aou phosphate phosphate_star	4.256261	3.030864	2.969139	JGOFS_GO	7.287125
Indo-Pacific	(27.85,27.95]	tem aou oxygen phosphate phosphate_star	4.255961	3.033021	2.960664	JGOFS_GO	7.288981
Indo-Pacific	(27.85,27.95]	tem aou oxygen silicate phosphate	4.280774	3.000763	2.891621	JGOFS_GO	7.281537
Indo-Pacific	(27.85,27.95]	tem aou silicate phosphate phosphate_star	4.274004	2.989254	2.887257	JGOFS_GO	7.263258
Indo-Pacific	(27.85,27.95]	tem oxygen silicate phosphate phosphate_star	4.294333	2.994511	2.882504	JGOFS_GO	7.288844
Indo-Pacific	(27.95,28.05]	aou oxygen phosphate phosphate_star	4.378771	2.877804	3.077571	GO_new	5.955375
Indo-Pacific	(27.95,28.05]	aou oxygen silicate phosphate	4.381257	2.867418	3.113633	GO_new	5.981051
Indo-Pacific	(27.95,28.05]	aou oxygen silicate phosphate phosphate_star	4.377427	2.861323	3.076725	GO_new	5.938048
Indo-Pacific	(27.95,28.05]	sal aou oxygen phosphate phosphate_star	4.359852	2.877276	3.073963	GO_new	5.951239
Indo-Pacific	(27.95,28.05]	sal aou oxygen silicate phosphate	4.359139	2.857154	3.075517	GO_new	5.932671
Indo-Pacific	(27.95,28.05]	sal tem aou oxygen phosphate	4.359973	2.871354	3.107381	GO_new	5.978735
Indo-Pacific	(27.95,28.05]	sal tem aou silicate phosphate	4.377696	2.882277	3.064828	GO_new	5.947104
Indo-Pacific	(27.95,28.05]	tem aou oxygen phosphate phosphate_star	4.378207	2.870660	3.073851	GO_new	5.944512
Indo-Pacific	(27.95,28.05]	tem aou oxygen silicate phosphate	4.381237	2.866160	3.109301	GO_new	5.975461
Indo-Pacific	(27.95,28.05]	tem aou silicate phosphate phosphate_star	4.401881	2.886861	3.065987	GO_new	5.952847
Indo-Pacific	(27.95,28.05]	aou oxygen phosphate phosphate_star	4.378771	2.877804	3.077571	JGOFS_GO	7.256575
Indo-Pacific	(27.95,28.05]	aou oxygen silicate phosphate	4.381257	2.867418	3.113633	JGOFS_GO	7.248675
Indo-Pacific	(27.95,28.05]	aou oxygen silicate phosphate phosphate_star	4.377427	2.861323	3.076725	JGOFS_GO	7.238750
Indo-Pacific	(27.95,28.05]	sal aou oxygen phosphate	4.360181	2.879287	3.107563	JGOFS_GO	7.239469
Indo-Pacific	(27.95,28.05]	sal aou oxygen phosphate phosphate_star	4.359852	2.877276	3.073963	JGOFS_GO	7.237128
Indo-Pacific	(27.95,28.05]	sal aou oxygen silicate phosphate	4.359139	2.857154	3.075517	JGOFS_GO	7.216293
Indo-Pacific	(27.95,28.05]	sal tem aou oxygen phosphate	4.359973	2.871354	3.107381	JGOFS_GO	7.231328
Indo-Pacific	(27.95,28.05]	tem aou oxygen phosphate	4.383037	2.871954	3.114960	JGOFS_GO	7.254991
Indo-Pacific	(27.95,28.05]	tem aou oxygen phosphate phosphate_star	4.378207	2.870660	3.073851	JGOFS_GO	7.248867
Indo-Pacific	(27.95,28.05]	tem aou oxygen silicate phosphate	4.381237	2.866160	3.109301	JGOFS_GO	7.247397
Indo-Pacific	(28.05,28.1]	aou oxygen phosphate phosphate_star	4.232069	2.875188	2.849672	GO_new	5.724859
Indo-Pacific	(28.05,28.1]	aou oxygen silicate phosphate phosphate_star	4.227177	2.874648	2.825582	GO_new	5.700230
Indo-Pacific	(28.05,28.1]	sal aou oxygen phosphate	4.228221	2.863082	2.828633	GO_new	5.691715
Indo-Pacific	(28.05,28.1]	sal aou oxygen phosphate phosphate_star	4.228217	2.860046	2.799367	GO_new	5.659413
Indo-Pacific	(28.05,28.1]	sal aou oxygen silicate phosphate	4.225670	2.848200	2.826261	GO_new	5.674461
Indo-Pacific	(28.05,28.1]	sal aou silicate phosphate	4.362882	2.877440	2.836183	GO_new	5.713623
Indo-Pacific	(28.05,28.1]	sal aou silicate phosphate phosphate_star	4.362458	2.865816	2.804566	GO_new	5.670382
Indo-Pacific	(28.05,28.1]	sal tem aou oxygen phosphate	4.225054	2.860776	2.817179	GO_new	5.677955
Indo-Pacific	(28.05,28.1]	sal tem aou silicate phosphate	4.264970	2.856439	2.836178	GO_new	5.692618
Indo-Pacific	(28.05,28.1]	tem aou oxygen phosphate phosphate_star	4.227073	2.874723	2.843130	GO_new	5.717853
Indo-Pacific	(28.05,28.1]	aou oxygen phosphate phosphate_star	4.232069	2.875188	2.849672	JGOFS_GO	7.107257
Indo-Pacific	(28.05,28.1]	aou oxygen silicate phosphate	4.228213	2.881751	2.870188	JGOFS_GO	7.109964
Indo-Pacific	(28.05,28.1]	aou oxygen silicate phosphate phosphate_star	4.227177	2.874648	2.825582	JGOFS_GO	7.101825
Indo-Pacific	(28.05,28.1]	sal aou oxygen phosphate	4.228221	2.863082	2.828633	JGOFS_GO	7.091304
Indo-Pacific	(28.05,28.1]	sal aou oxygen phosphate phosphate_star	4.228217	2.860046	2.799367	JGOFS_GO	7.088264
Indo-Pacific	(28.05,28.1]	sal aou oxygen silicate phosphate	4.225670	2.848200	2.826261	JGOFS_GO	7.073870
Indo-Pacific	(28.05,28.1]	sal tem aou oxygen phosphate	4.225054	2.860776	2.817179	JGOFS_GO	7.085830
Indo-Pacific	(28.05,28.1]	tem aou oxygen phosphate	4.227670	2.881411	2.893756	JGOFS_GO	7.109081
Indo-Pacific	(28.05,28.1]	tem aou oxygen phosphate phosphate_star	4.227073	2.874723	2.843130	JGOFS_GO	7.101796
Indo-Pacific	(28.05,28.1]	tem aou oxygen silicate phosphate	4.226770	2.880685	2.860748	JGOFS_GO	7.107455
Indo-Pacific	(28.1, Inf]	aou oxygen silicate phosphate	3.640879	3.028695	2.696173	GO_new	5.724868
Indo-Pacific	(28.1, Inf]	aou oxygen silicate phosphate phosphate_star	3.635023	2.938152	2.688069	GO_new	5.626221
Indo-Pacific	(28.1, Inf]	sal aou oxygen phosphate	3.599961	3.230532	2.595467	GO_new	5.825999
Indo-Pacific	(28.1, Inf]	sal aou oxygen phosphate phosphate_star	3.599599	3.136083	2.593341	GO_new	5.729425
Indo-Pacific	(28.1, Inf]	sal aou oxygen silicate phosphate	3.570328	3.028688	2.582084	GO_new	5.610771
Indo-Pacific	(28.1, Inf]	sal aou phosphate phosphate_star	3.713417	3.207210	2.648134	GO_new	5.855343
Indo-Pacific	(28.1, Inf]	sal aou silicate phosphate phosphate_star	3.698036	3.187854	2.643218	GO_new	5.831072
Indo-Pacific	(28.1, Inf]	sal tem aou oxygen phosphate	3.540683	3.225881	2.586405	GO_new	5.812286
Indo-Pacific	(28.1, Inf]	sal tem aou phosphate phosphate_star	3.540570	3.173834	2.644631	GO_new	5.818464
Indo-Pacific	(28.1, Inf]	tem aou oxygen silicate phosphate	3.601769	2.985715	2.648408	GO_new	5.634123
Indo-Pacific	(28.1, Inf]	aou oxygen silicate phosphate	3.640879	3.028695	2.696173	JGOFS_GO	6.669574
Indo-Pacific	(28.1, Inf]	aou oxygen silicate phosphate phosphate_star	3.635023	2.938152	2.688069	JGOFS_GO	6.573176
Indo-Pacific	(28.1, Inf]	sal aou oxygen phosphate phosphate_star	3.599599	3.136083	2.593341	JGOFS_GO	6.735682
Indo-Pacific	(28.1, Inf]	sal aou oxygen silicate phosphate	3.570328	3.028688	2.582084	JGOFS_GO	6.599015
Indo-Pacific	(28.1, Inf]	sal tem aou oxygen phosphate	3.540683	3.225881	2.586405	JGOFS_GO	6.766564
Indo-Pacific	(28.1, Inf]	sal tem aou phosphate phosphate_star	3.540570	3.173834	2.644631	JGOFS_GO	6.714404
Indo-Pacific	(28.1, Inf]	sal tem aou silicate phosphate	3.542161	3.278460	2.644383	JGOFS_GO	6.820621
Indo-Pacific	(28.1, Inf]	sal tem oxygen phosphate phosphate_star	3.594199	3.222924	2.710942	JGOFS_GO	6.817124
Indo-Pacific	(28.1, Inf]	tem aou oxygen silicate phosphate	3.601769	2.985715	2.648408	JGOFS_GO	6.587484
Indo-Pacific	(28.1, Inf]	tem aou silicate phosphate phosphate_star	3.653818	3.124788	2.839805	JGOFS_GO	6.778606

3.3 Correlations RMSE between eras

3.3.1 All models

all_lm_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(all_lm_stats, all_lm_stats_wide)

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

#rm(all_lm_stats, all_lm_stats_wide)

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

lm_best <- lm_best %>% 
  mutate(lm_coeff = str_replace_all(predictors, " ", " + "),
         lm_coeff = paste("Cstar ~", lm_coeff))

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 <- lm_best %>%
      filter(basin == i_basin)

    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_slab <- lm_best_basin %>% 
        filter(gamma_slab == i_gamma_slab)
      
      for (i_eras in unique(lm_best_basin_slab$eras)) {
        
        # i_eras <- unique(lm_best_basin_slab$eras)[1]
        
        lm_best_basin_slab_eras <- lm_best_basin_slab %>%
          filter(eras == i_eras)
      
      for (i_predictors in unique(lm_best_basin_slab_eras$predictors)) {
        
        # i_predictors <- unique(lm_best_basin_slab_eras$predictors)[1]
        
        lm_best_single <- lm_best_basin_slab_eras %>%
          filter(predictors == i_predictors) %>%
          select(lm_coeff) %>%
          pull()
        
        i_rmse_sum <- lm_best_basin_slab_eras %>%
          filter(predictors == i_predictors) %>%
          select(rmse_sum) %>%
          pull()
        
        lm_formula <- as.formula(lm_best_single)
        
        model_single <- lm(lm_formula, 
              data = GLODAP_basin_era_slab)
        
        coefficients <- tidy(model_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)
        
        print(lm_best_single)
      
      if (exists("all_lm")) {
        all_lm <- bind_rows(all_lm, coefficients)
      }
      
      if (!exists("all_lm")) {
        all_lm <- coefficients
      }
      

      }
      
 
      }
    }
    
  }
  
}


rm(lm_formula, model_single, coefficients, i_rmse_sum,
   GLODAP_basin_era, GLODAP_basin_era_slab,
   lm_best, lm_best_basin, lm_best_basin_slab, lm_best_single,
   i_basin, i_era, i_gamma_slab, i_predictors)

all_lm %>%
  write_csv(here::here("data/eMLR",
                       "all_lm.csv"))


# Prepare model coefficients for mapping

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

all_lm <- all_lm %>% 
  mutate(estimate = if_else(is.na(estimate), 0, estimate))

all_lm_wide <- all_lm %>% 
  pivot_wider(names_from = era, values_from = estimate,
              names_prefix = "coeff_")

all_lm_wide <- all_lm_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))

all_lm_long <- all_lm_wide %>% 
  pivot_longer(JGOFS_GO:GO_new, names_to = "eras_fit", values_to = "delta_coeff")

all_lm_long <- all_lm_long %>% 
  filter(eras == eras_fit) %>% 
  select(-eras_fit)

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

all_lm_wide %>%
  write_csv(here::here("data/eMLR",
                       "all_lm_wide.csv"))

rm(all_lm_long, all_lm_wide)

4 Open tasks

combine correlation, rmse and model calculation into one routine, ideally using the purrr map approach
select exactly best 10 models, even if same rmse sum appears more than once (currently, selection is set two ten, even if doublicated rmse occur)
check (phosphate_star) NA’s in coeffcients -> cases of exact colinearity?
test model selection based on AIC

5 Open questions

sessionInfo()

R version 4.0.2 (2020-06-22)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows 10 x64 (build 18363)

Matrix products: default

locale:
[1] LC_COLLATE=English_Germany.1252  LC_CTYPE=English_Germany.1252   
[3] LC_MONETARY=English_Germany.1252 LC_NUMERIC=C                    
[5] LC_TIME=English_Germany.1252    

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
 [1] corrr_0.4.2      broom_0.7.0      kableExtra_1.1.0 knitr_1.29      
 [5] olsrr_0.5.3.9000 GGally_2.0.0     patchwork_1.0.1  lubridate_1.7.9 
 [9] forcats_0.5.0    stringr_1.4.0    dplyr_1.0.0      purrr_0.3.4     
[13] readr_1.3.1      tidyr_1.1.0      tibble_3.0.3     ggplot2_3.3.2   
[17] tidyverse_1.3.0  workflowr_1.6.2 

loaded via a namespace (and not attached):
 [1] fs_1.4.2           webshot_0.5.2      RColorBrewer_1.1-2 httr_1.4.2        
 [5] rprojroot_1.3-2    tools_4.0.2        backports_1.1.8    R6_2.4.1          
 [9] nortest_1.0-4      DBI_1.1.0          colorspace_1.4-1   withr_2.2.0       
[13] tidyselect_1.1.0   gridExtra_2.3      curl_4.3           compiler_4.0.2    
[17] git2r_0.27.1       cli_2.0.2          rvest_0.3.6        xml2_1.3.2        
[21] labeling_0.3       scales_1.1.1       goftest_1.2-2      digest_0.6.25     
[25] foreign_0.8-80     rmarkdown_2.3      rio_0.5.16         pkgconfig_2.0.3   
[29] htmltools_0.5.0    dbplyr_1.4.4       highr_0.8          rlang_0.4.7       
[33] readxl_1.3.1       rstudioapi_0.11    generics_0.0.2     farver_2.0.3      
[37] jsonlite_1.7.0     zip_2.0.4          car_3.0-8          magrittr_1.5      
[41] Rcpp_1.0.5         munsell_0.5.0      fansi_0.4.1        abind_1.4-5       
[45] lifecycle_0.2.0    stringi_1.4.6      whisker_0.4        yaml_2.2.1        
[49] carData_3.0-4      plyr_1.8.6         grid_4.0.2         blob_1.2.1        
[53] promises_1.1.1     crayon_1.3.4       haven_2.3.1        hms_0.5.3         
[57] pillar_1.4.6       reprex_0.3.0       glue_1.4.1         evaluate_0.14     
[61] data.table_1.13.0  modelr_0.1.8       vctrs_0.3.2        httpuv_1.5.4      
[65] cellranger_1.1.0   gtable_0.3.0       reshape_0.8.8      assertthat_0.2.1  
[69] xfun_0.16          openxlsx_4.1.5     later_1.1.0.1      viridisLite_0.3.0 
[73] ellipsis_0.3.1     here_0.1

eMLR - model fitting

Jens Daniel Müller

28 August, 2020