Last updated: 2020-12-23

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: a4531df

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 a4531df. 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_assumption_testing.Rmd) and HTML (docs/eMLR_assumption_testing.html) files. If you’ve configured a remote Git repository (see ?wflow_git_remote), click on the hyperlinks in the table below to view the files as they were in that past version.

File Version Author Date Message
html 0aa2b50 jens-daniel-mueller 2020-12-23 remove html before duplication
html 39113c3 jens-daniel-mueller 2020-12-23 Build site.
html 2886da0 jens-daniel-mueller 2020-12-19 Build site.
html 02f0ee9 jens-daniel-mueller 2020-12-18 cleaned up for copying template
html 965dba3 jens-daniel-mueller 2020-12-18 Build site.
html 5d452fe jens-daniel-mueller 2020-12-18 Build site.
Rmd ca65bf5 jens-daniel-mueller 2020-12-18 rebuild after final cleaning
html 7bcb4eb jens-daniel-mueller 2020-12-18 Build site.
html d397028 jens-daniel-mueller 2020-12-18 Build site.
html 7131186 jens-daniel-mueller 2020-12-17 Build site.
html 22b07fb jens-daniel-mueller 2020-12-17 Build site.
html f3a708f jens-daniel-mueller 2020-12-17 Build site.
html e4ca289 jens-daniel-mueller 2020-12-16 Build site.
html 158fe26 jens-daniel-mueller 2020-12-15 Build site.
html 7a9a4cb jens-daniel-mueller 2020-12-15 Build site.
html 61b263c jens-daniel-mueller 2020-12-15 Build site.
html 4d612dd jens-daniel-mueller 2020-12-15 Build site.
Rmd e7e5ff1 jens-daniel-mueller 2020-12-15 rebuild with eMLR target variable selection
html 953caf3 jens-daniel-mueller 2020-12-15 Build site.
html 42daf5c jens-daniel-mueller 2020-12-14 Build site.
Rmd 923aa7f jens-daniel-mueller 2020-12-14 rebuild with new path and auto folder creation
html 984697e jens-daniel-mueller 2020-12-12 Build site.
html 3ebff89 jens-daniel-mueller 2020-12-12 Build site.
html 5d96d3c jens-daniel-mueller 2020-12-11 Build site.
Rmd 3d33a37 jens-daniel-mueller 2020-12-11 selectable basinmask, try 5
html b01a367 jens-daniel-mueller 2020-12-09 Build site.
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.
Rmd 33b1973 jens-daniel-mueller 2020-12-07 run with WOCE flag 2 only
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 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
html 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 Predictor correlation

The correlation between:

  • pairs of seven potential predictor variables and
  • C* and seven potential predictor variables

were investigated based on:

  • property-property plots and
  • calculated correlation coeffcients.

2.1 Correlation plots

For an overview, a random subset of data from all eras was plotted separately for both basins, with color indicating neutral density slabs (high density = dark-purple color).

for (i_basin in unique(GLODAP$basin)) {
# i_basin <- unique(GLODAP$basin)[1]

print(
GLODAP %>% 
  filter(basin == i_basin) %>% 
  sample_frac(0.05) %>% 
  ggpairs(columns = c(params_local$MLR_target,
                      params_local$MLR_predictors),
          upper = "blank",
          ggplot2::aes(col = gamma_slab, fill = gamma_slab, alpha = 0.01)) +
      scale_fill_viridis_d(direction = -1) +
      scale_color_viridis_d(direction = -1) +
      labs(title = paste("Basin:", i_basin ,"| era: all | subsample size: 5 % of",
                         nrow(GLODAP %>% filter(basin == i_basin))))
)
}

Version Author Date
0aa2b50 jens-daniel-mueller 2020-12-23
2886da0 jens-daniel-mueller 2020-12-19
02f0ee9 jens-daniel-mueller 2020-12-18
7bcb4eb jens-daniel-mueller 2020-12-18
158fe26 jens-daniel-mueller 2020-12-15
7a9a4cb jens-daniel-mueller 2020-12-15
61b263c jens-daniel-mueller 2020-12-15
4d612dd jens-daniel-mueller 2020-12-15
984697e jens-daniel-mueller 2020-12-12
3ebff89 jens-daniel-mueller 2020-12-12
5d96d3c jens-daniel-mueller 2020-12-11

Version Author Date
0aa2b50 jens-daniel-mueller 2020-12-23
2886da0 jens-daniel-mueller 2020-12-19
02f0ee9 jens-daniel-mueller 2020-12-18
7bcb4eb jens-daniel-mueller 2020-12-18
158fe26 jens-daniel-mueller 2020-12-15
7a9a4cb jens-daniel-mueller 2020-12-15
61b263c jens-daniel-mueller 2020-12-15
4d612dd jens-daniel-mueller 2020-12-15
984697e jens-daniel-mueller 2020-12-12
3ebff89 jens-daniel-mueller 2020-12-12
5d96d3c jens-daniel-mueller 2020-12-11

Individual correlation plots for each basin, era and neutral density (gamma) slab are available at:

/nfs/kryo/work/jenmueller/emlr_cant/observations/v_XXX/figures/Observations_correlation/

if (params_local$plot_all_figures == "y") {

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]

    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)[5]

      GLODAP_highlight <- GLODAP_basin_era %>%
        mutate(gamma_highlight = if_else(gamma_slab == i_gamma_slab,
                                         "in", "out")) %>%
        arrange(desc(gamma_highlight))
      
      p <- GLODAP_highlight %>%
        ggpairs(
          columns = c(params_local$MLR_target,
                      params_local$MLR_predictors),
          ggplot2::aes(
            col = gamma_highlight,
            fill = gamma_highlight,
            alpha = 0.01
          )
        ) +
        scale_fill_manual(values = c("red", "grey")) +
        scale_color_manual(values = c("red", "grey")) +
        labs(
          title = paste(
            i_era,
            "|",
            i_basin,
            "| Gamma slab",
            i_gamma_slab,
            "| # obs total",
            nrow(GLODAP_basin_era),
            "| # obs slab",
            nrow(GLODAP_highlight %>%
                   filter(gamma_highlight == "in"))
          )
        )
      
      png(
        filename = paste(
          path_version_figures,
          "Observations_correlation/",
          paste(
            "Predictor_correlation",
            i_era,
            i_basin,
            i_gamma_slab,
            ".png",
            sep = "_"
          ),
          sep = ""),
        width = 12,
        height = 12,
        units = "in",
        res = 300
      )
      
      print(p)
      
      dev.off()
      
    }
  }
}

}

2.2 Correlation assesment

2.2.1 Calculation of correlation coeffcients

Correlation coefficients were calculated individually within each slabs, era and basin.

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]
    
    GLODAP_basin_era <- GLODAP %>%
      filter(basin == i_basin,
             era == i_era) %>%
      select(basin,
             era,
             gamma_slab,
             params_local$MLR_target,
             params_local$MLR_predictors)
    
    for (i_gamma_slab in unique(GLODAP_basin_era$gamma_slab)) {
      # i_gamma_slab <- unique(GLODAP_basin_era$gamma_slab)[5]
      print(i_gamma_slab)
      
      GLODAP_basin_era_slab <- GLODAP_basin_era %>%
        filter(gamma_slab == i_gamma_slab)
      
      # calculate correlation table
      cor_target_predictor_temp <- GLODAP_basin_era_slab %>%
        select(-c(basin, era, gamma_slab)) %>%
        correlate() %>%
        focus(params_local$MLR_target) %>%
        mutate(basin = i_basin,
               era = i_era,
               gamma_slab = i_gamma_slab)
      
      if (exists("cor_target_predictor")) {
        cor_target_predictor <-
          bind_rows(cor_target_predictor, cor_target_predictor_temp)
      }
      
      if (!exists("cor_target_predictor")) {
        cor_target_predictor <- cor_target_predictor_temp
      }
      
      
      cor_predictors_temp <- GLODAP_basin_era_slab %>%
        select(-c(basin, era, gamma_slab)) %>%
        correlate() %>%
        shave %>%
        stretch() %>%
        filter(!is.na(r),
               x != params_local$MLR_target,
               y != params_local$MLR_target) %>%
        mutate(pair = paste(x, y, sep = " + ")) %>%
        select(-c(x, y)) %>%
        mutate(basin = i_basin,
               era = i_era,
               gamma_slab = i_gamma_slab)
      
      if (exists("cor_predictors")) {
        cor_predictors <- bind_rows(cor_predictors, cor_predictors_temp)
      }
      
      if (!exists("cor_predictors")) {
        cor_predictors <- cor_predictors_temp
      }
      
      
      
    }
  }
}

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

2.2.2 Predictor pairs

Below, the range of correlations coefficients for each predictor pair is plotted per basin (facet) and density slab (color). Note that the range indicates the min and max values of in total 3 calculated coefficients (one per era).

# calculate min, max, mean across all eras
cor_predictors_stats <- cor_predictors %>% 
  group_by(pair, basin, gamma_slab) %>% 
  summarise(mean_r = mean(r),
            min_r = min(r),
            max_r = max(r)) %>% 
  ungroup()

# plot figure
cor_predictors_stats %>% 
  mutate(pair = reorder(pair, mean_r)) %>%
  ggplot() +
  geom_vline(xintercept = c(-0.9, 0.9), col = "red") +
  geom_vline(xintercept = 0) +
  geom_linerange(
    aes(y = pair, xmin = min_r, xmax = max_r, col = gamma_slab),
    position = position_dodge(width = 0.6)) +
  facet_wrap(~basin) +
  scale_color_viridis_d(direction = -1) +
  labs(x = "correlation coefficient", y = "") +
  theme(legend.position = "top")

Version Author Date
0aa2b50 jens-daniel-mueller 2020-12-23
2886da0 jens-daniel-mueller 2020-12-19
02f0ee9 jens-daniel-mueller 2020-12-18
158fe26 jens-daniel-mueller 2020-12-15
984697e jens-daniel-mueller 2020-12-12
3ebff89 jens-daniel-mueller 2020-12-12
5d96d3c jens-daniel-mueller 2020-12-11
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
0ff728b jens-daniel-mueller 2020-12-01
91435ae jens-daniel-mueller 2020-12-01
196be51 jens-daniel-mueller 2020-11-30
bc61ce3 Jens Müller 2020-11-30 
# print table
kable(cor_predictors_stats) %>%
  add_header_above() %>%
  kable_styling() %>%
  scroll_box(width = "100%", height = "400px")
pair basin gamma_slab mean_r min_r max_r
aou + nitrate Atlantic (-Inf,26] 0.8882601 0.8639438 0.9185712
aou + nitrate Atlantic (26,26.5] 0.8994326 0.8480338 0.9384440
aou + nitrate Atlantic (26.5,26.75] 0.9486287 0.9364311 0.9728018
aou + nitrate Atlantic (26.75,27] 0.9286820 0.8992595 0.9493538
aou + nitrate Atlantic (27,27.25] 0.8637480 0.7686626 0.9114722
aou + nitrate Atlantic (27.25,27.5] 0.7884053 0.7079764 0.8534782
aou + nitrate Atlantic (27.5,27.75] 0.8712005 0.8376987 0.8958752
aou + nitrate Atlantic (27.75,27.85] 0.9580498 0.9483592 0.9645514
aou + nitrate Atlantic (27.85,27.95] 0.9620424 0.9488690 0.9708406
aou + nitrate Atlantic (27.95,28.05] 0.9630959 0.9448219 0.9858649
aou + nitrate Atlantic (28.05,28.1] 0.9672424 0.9419501 0.9845767
aou + nitrate Atlantic (28.1,28.15] 0.9789566 0.9624539 0.9886620
aou + nitrate Atlantic (28.15,28.2] 0.9819641 0.9590050 0.9950555
aou + nitrate Atlantic (28.2, Inf] 0.8561810 0.7589899 0.9375717
aou + nitrate Indo-Pacific (-Inf,26] 0.9220556 0.8855553 0.9612879
aou + nitrate Indo-Pacific (26,26.5] 0.9058901 0.8969654 0.9172197
aou + nitrate Indo-Pacific (26.5,26.75] 0.9161520 0.8968438 0.9299319
aou + nitrate Indo-Pacific (26.75,27] 0.9396149 0.9348278 0.9431014
aou + nitrate Indo-Pacific (27,27.25] 0.9458293 0.9015387 0.9688879
aou + nitrate Indo-Pacific (27.25,27.5] 0.9353774 0.8812774 0.9712891
aou + nitrate Indo-Pacific (27.5,27.75] 0.9102652 0.8963528 0.9244390
aou + nitrate Indo-Pacific (27.75,27.85] 0.9713027 0.9423407 0.9908923
aou + nitrate Indo-Pacific (27.85,27.95] 0.9748306 0.9569698 0.9910972
aou + nitrate Indo-Pacific (27.95,28.05] 0.9604059 0.9432965 0.9765711
aou + nitrate Indo-Pacific (28.05,28.1] 0.9482872 0.9301783 0.9683823
aou + nitrate Indo-Pacific (28.1, Inf] 0.7810604 0.7215916 0.8611870
aou + silicate Atlantic (-Inf,26] 0.2880178 0.0815925 0.4890832
aou + silicate Atlantic (26,26.5] 0.6885617 0.5569485 0.8265329
aou + silicate Atlantic (26.5,26.75] 0.8941336 0.8634956 0.9228433
aou + silicate Atlantic (26.75,27] 0.8731907 0.8417628 0.9031089
aou + silicate Atlantic (27,27.25] 0.7566277 0.6556533 0.8266663
aou + silicate Atlantic (27.25,27.5] 0.5244898 0.4250217 0.6798948
aou + silicate Atlantic (27.5,27.75] 0.6118017 0.5908351 0.6468691
aou + silicate Atlantic (27.75,27.85] 0.8482958 0.8440302 0.8558311
aou + silicate Atlantic (27.85,27.95] 0.8973798 0.8861594 0.9120410
aou + silicate Atlantic (27.95,28.05] 0.9432373 0.9253467 0.9635956
aou + silicate Atlantic (28.05,28.1] 0.9659902 0.9470616 0.9825390
aou + silicate Atlantic (28.1,28.15] 0.9732607 0.9695009 0.9761102
aou + silicate Atlantic (28.15,28.2] 0.9869606 0.9751151 0.9942627
aou + silicate Atlantic (28.2, Inf] 0.8825833 0.8456145 0.9169515
aou + silicate Indo-Pacific (-Inf,26] 0.7347218 0.6602592 0.8060948
aou + silicate Indo-Pacific (26,26.5] 0.6555734 0.5384547 0.7555346
aou + silicate Indo-Pacific (26.5,26.75] 0.6138018 0.4137384 0.7294312
aou + silicate Indo-Pacific (26.75,27] 0.7486668 0.6950166 0.7800269
aou + silicate Indo-Pacific (27,27.25] 0.8627673 0.8037330 0.9081061
aou + silicate Indo-Pacific (27.25,27.5] 0.8997085 0.8537846 0.9374421
aou + silicate Indo-Pacific (27.5,27.75] 0.9209233 0.8699565 0.9507471
aou + silicate Indo-Pacific (27.75,27.85] 0.9673568 0.9529991 0.9761431
aou + silicate Indo-Pacific (27.85,27.95] 0.9634829 0.9540446 0.9701809
aou + silicate Indo-Pacific (27.95,28.05] 0.9482145 0.9414261 0.9524593
aou + silicate Indo-Pacific (28.05,28.1] 0.9464987 0.9451320 0.9491329
aou + silicate Indo-Pacific (28.1, Inf] 0.6413539 0.5433113 0.6934260
sal + aou Atlantic (-Inf,26] -0.0582286 -0.3694296 0.3402213
sal + aou Atlantic (26,26.5] -0.1350036 -0.2309276 0.0213537
sal + aou Atlantic (26.5,26.75] -0.3233372 -0.4727704 -0.1936560
sal + aou Atlantic (26.75,27] -0.1921542 -0.2257555 -0.1724587
sal + aou Atlantic (27,27.25] -0.1573448 -0.3045530 0.0189123
sal + aou Atlantic (27.25,27.5] -0.1745433 -0.3192984 -0.0824784
sal + aou Atlantic (27.5,27.75] -0.4478260 -0.4883188 -0.3776379
sal + aou Atlantic (27.75,27.85] -0.5338527 -0.5619584 -0.5018946
sal + aou Atlantic (27.85,27.95] -0.3851576 -0.4968361 -0.2625017
sal + aou Atlantic (27.95,28.05] -0.5327608 -0.7918401 -0.2080327
sal + aou Atlantic (28.05,28.1] -0.7288595 -0.8677647 -0.6500999
sal + aou Atlantic (28.1,28.15] -0.8060559 -0.8710340 -0.7594185
sal + aou Atlantic (28.15,28.2] -0.8987392 -0.9757986 -0.7583391
sal + aou Atlantic (28.2, Inf] -0.6318306 -0.8429734 -0.3871254
sal + aou Indo-Pacific (-Inf,26] -0.1024681 -0.1364360 -0.0346833
sal + aou Indo-Pacific (26,26.5] -0.0415126 -0.2469743 0.1386419
sal + aou Indo-Pacific (26.5,26.75] -0.0446696 -0.3659553 0.2606167
sal + aou Indo-Pacific (26.75,27] -0.1499138 -0.3502232 0.0393276
sal + aou Indo-Pacific (27,27.25] -0.0745481 -0.1826713 0.1405027
sal + aou Indo-Pacific (27.25,27.5] 0.1282204 0.0609408 0.2542954
sal + aou Indo-Pacific (27.5,27.75] 0.2796244 0.1327307 0.4868767
sal + aou Indo-Pacific (27.75,27.85] -0.1678690 -0.4059917 0.0695052
sal + aou Indo-Pacific (27.85,27.95] -0.3374984 -0.6468553 -0.1153630
sal + aou Indo-Pacific (27.95,28.05] -0.3231164 -0.3998309 -0.2344439
sal + aou Indo-Pacific (28.05,28.1] -0.3769788 -0.4628420 -0.2746149
sal + aou Indo-Pacific (28.1, Inf] 0.4347919 0.3263842 0.4898025
sal + nitrate Atlantic (-Inf,26] -0.2164247 -0.3742820 0.0480366
sal + nitrate Atlantic (26,26.5] -0.4619060 -0.4992552 -0.4181418
sal + nitrate Atlantic (26.5,26.75] -0.5691143 -0.6403213 -0.4828679
sal + nitrate Atlantic (26.75,27] -0.5147222 -0.5676866 -0.4794929
sal + nitrate Atlantic (27,27.25] -0.6033479 -0.6575971 -0.5577427
sal + nitrate Atlantic (27.25,27.5] -0.7220873 -0.7548862 -0.6684211
sal + nitrate Atlantic (27.5,27.75] -0.7790548 -0.7959070 -0.7559585
sal + nitrate Atlantic (27.75,27.85] -0.7149322 -0.7233348 -0.7106048
sal + nitrate Atlantic (27.85,27.95] -0.5829852 -0.6505871 -0.4734760
sal + nitrate Atlantic (27.95,28.05] -0.6852242 -0.8429564 -0.4226800
sal + nitrate Atlantic (28.05,28.1] -0.8328846 -0.8995744 -0.7632877
sal + nitrate Atlantic (28.1,28.15] -0.8744180 -0.8984825 -0.8419138
sal + nitrate Atlantic (28.15,28.2] -0.9370910 -0.9635752 -0.8894399
sal + nitrate Atlantic (28.2, Inf] -0.9164385 -0.9664373 -0.8528289
sal + nitrate Indo-Pacific (-Inf,26] -0.3533354 -0.3801615 -0.3132694
sal + nitrate Indo-Pacific (26,26.5] -0.3563034 -0.5710368 -0.1890833
sal + nitrate Indo-Pacific (26.5,26.75] -0.3463872 -0.6376549 -0.0844787
sal + nitrate Indo-Pacific (26.75,27] -0.4130527 -0.5889828 -0.2679892
sal + nitrate Indo-Pacific (27,27.25] -0.2935536 -0.3471093 -0.2301284
sal + nitrate Indo-Pacific (27.25,27.5] -0.0446170 -0.1546454 0.0598255
sal + nitrate Indo-Pacific (27.5,27.75] 0.2666655 -0.1448972 0.7280035
sal + nitrate Indo-Pacific (27.75,27.85] -0.2908752 -0.4234049 -0.1769294
sal + nitrate Indo-Pacific (27.85,27.95] -0.4520102 -0.6726286 -0.3005614
sal + nitrate Indo-Pacific (27.95,28.05] -0.5047786 -0.5421106 -0.4562052
sal + nitrate Indo-Pacific (28.05,28.1] -0.5884947 -0.6297631 -0.5256634
sal + nitrate Indo-Pacific (28.1, Inf] -0.0920071 -0.1195623 -0.0596777
sal + silicate Atlantic (-Inf,26] -0.7229515 -0.7963213 -0.6725144
sal + silicate Atlantic (26,26.5] -0.6666292 -0.6961385 -0.6125911
sal + silicate Atlantic (26.5,26.75] -0.6234854 -0.7013916 -0.5055632
sal + silicate Atlantic (26.75,27] -0.5271098 -0.5482018 -0.5082031
sal + silicate Atlantic (27,27.25] -0.6502900 -0.6804575 -0.6223765
sal + silicate Atlantic (27.25,27.5] -0.7930190 -0.8100851 -0.7824100
sal + silicate Atlantic (27.5,27.75] -0.8131489 -0.8323011 -0.7794771
sal + silicate Atlantic (27.75,27.85] -0.7939878 -0.8147856 -0.7780823
sal + silicate Atlantic (27.85,27.95] -0.6904396 -0.7405664 -0.6125805
sal + silicate Atlantic (27.95,28.05] -0.7312631 -0.8881178 -0.4822718
sal + silicate Atlantic (28.05,28.1] -0.8568889 -0.9273929 -0.8074986
sal + silicate Atlantic (28.1,28.15] -0.8947062 -0.9387988 -0.8695910
sal + silicate Atlantic (28.15,28.2] -0.9247385 -0.9729301 -0.8350616
sal + silicate Atlantic (28.2, Inf] -0.8796120 -0.9797038 -0.7821316
sal + silicate Indo-Pacific (-Inf,26] -0.5262763 -0.5734616 -0.4400617
sal + silicate Indo-Pacific (26,26.5] -0.5585775 -0.6796337 -0.4557575
sal + silicate Indo-Pacific (26.5,26.75] -0.6683266 -0.8095328 -0.5574581
sal + silicate Indo-Pacific (26.75,27] -0.6615296 -0.7640584 -0.4927603
sal + silicate Indo-Pacific (27,27.25] -0.4669166 -0.5431525 -0.3513330
sal + silicate Indo-Pacific (27.25,27.5] -0.1588531 -0.1968028 -0.1146590
sal + silicate Indo-Pacific (27.5,27.75] 0.1108211 -0.0798353 0.3808631
sal + silicate Indo-Pacific (27.75,27.85] -0.2197543 -0.3581084 -0.0755570
sal + silicate Indo-Pacific (27.85,27.95] -0.3710245 -0.5933264 -0.2065742
sal + silicate Indo-Pacific (27.95,28.05] -0.4001010 -0.4457158 -0.3433892
sal + silicate Indo-Pacific (28.05,28.1] -0.4874820 -0.5271515 -0.4130308
sal + silicate Indo-Pacific (28.1, Inf] -0.0105220 -0.1576516 0.1555784
sal + temp Atlantic (-Inf,26] 0.9216414 0.8422006 0.9786562
sal + temp Atlantic (26,26.5] 0.9600484 0.9470189 0.9714949
sal + temp Atlantic (26.5,26.75] 0.9864911 0.9811261 0.9901624
sal + temp Atlantic (26.75,27] 0.9786516 0.9737876 0.9842533
sal + temp Atlantic (27,27.25] 0.9765338 0.9733752 0.9795262
sal + temp Atlantic (27.25,27.5] 0.9721161 0.9697519 0.9733243
sal + temp Atlantic (27.5,27.75] 0.9579710 0.9539439 0.9630964
sal + temp Atlantic (27.75,27.85] 0.9768976 0.9758545 0.9776904
sal + temp Atlantic (27.85,27.95] 0.9421017 0.9311908 0.9477183
sal + temp Atlantic (27.95,28.05] 0.8833438 0.8067063 0.9356883
sal + temp Atlantic (28.05,28.1] 0.9502657 0.9161398 0.9677069
sal + temp Atlantic (28.1,28.15] 0.9487678 0.9248997 0.9669567
sal + temp Atlantic (28.15,28.2] 0.8834110 0.8634964 0.9035039
sal + temp Atlantic (28.2, Inf] 0.2479905 0.1081040 0.4471581
sal + temp Indo-Pacific (-Inf,26] 0.7696141 0.7410403 0.8176761
sal + temp Indo-Pacific (26,26.5] 0.9287046 0.8849628 0.9536762
sal + temp Indo-Pacific (26.5,26.75] 0.9785285 0.9692090 0.9832380
sal + temp Indo-Pacific (26.75,27] 0.9770281 0.9721445 0.9798567
sal + temp Indo-Pacific (27,27.25] 0.9503584 0.9326264 0.9715997
sal + temp Indo-Pacific (27.25,27.5] 0.8810632 0.8282727 0.9587561
sal + temp Indo-Pacific (27.5,27.75] 0.8636152 0.7962705 0.9159071
sal + temp Indo-Pacific (27.75,27.85] 0.8135013 0.6546288 0.9398146
sal + temp Indo-Pacific (27.85,27.95] 0.7215049 0.4948250 0.8739283
sal + temp Indo-Pacific (27.95,28.05] 0.7435356 0.6602919 0.8266420
sal + temp Indo-Pacific (28.05,28.1] 0.7320792 0.6427383 0.8009832
sal + temp Indo-Pacific (28.1, Inf] 0.7274017 0.6248405 0.8291044
silicate + nitrate Atlantic (-Inf,26] 0.4974689 0.3482023 0.6168467
silicate + nitrate Atlantic (26,26.5] 0.8568300 0.8072608 0.9335862
silicate + nitrate Atlantic (26.5,26.75] 0.9679547 0.9606611 0.9766322
silicate + nitrate Atlantic (26.75,27] 0.9466879 0.9296813 0.9564043
silicate + nitrate Atlantic (27,27.25] 0.9266034 0.9104018 0.9439514
silicate + nitrate Atlantic (27.25,27.5] 0.8654226 0.8273825 0.9022854
silicate + nitrate Atlantic (27.5,27.75] 0.8683352 0.8529238 0.8803526
silicate + nitrate Atlantic (27.75,27.85] 0.9434500 0.9358477 0.9570492
silicate + nitrate Atlantic (27.85,27.95] 0.9678043 0.9636127 0.9742486
silicate + nitrate Atlantic (27.95,28.05] 0.9788492 0.9758933 0.9845790
silicate + nitrate Atlantic (28.05,28.1] 0.9880763 0.9857313 0.9903917
silicate + nitrate Atlantic (28.1,28.15] 0.9904300 0.9893853 0.9913419
silicate + nitrate Atlantic (28.15,28.2] 0.9929119 0.9863975 0.9971804
silicate + nitrate Atlantic (28.2, Inf] 0.9751292 0.9534844 0.9949288
silicate + nitrate Indo-Pacific (-Inf,26] 0.8805434 0.8631933 0.8923245
silicate + nitrate Indo-Pacific (26,26.5] 0.8191346 0.7423559 0.8788096
silicate + nitrate Indo-Pacific (26.5,26.75] 0.7822461 0.6314419 0.8683984
silicate + nitrate Indo-Pacific (26.75,27] 0.8456070 0.8019618 0.8703006
silicate + nitrate Indo-Pacific (27,27.25] 0.8916237 0.8670918 0.9249746
silicate + nitrate Indo-Pacific (27.25,27.5] 0.8755929 0.8521580 0.9189094
silicate + nitrate Indo-Pacific (27.5,27.75] 0.8423682 0.8158117 0.8606030
silicate + nitrate Indo-Pacific (27.75,27.85] 0.9463954 0.9209059 0.9598690
silicate + nitrate Indo-Pacific (27.85,27.95] 0.9586162 0.9424249 0.9687097
silicate + nitrate Indo-Pacific (27.95,28.05] 0.9497408 0.9345232 0.9667133
silicate + nitrate Indo-Pacific (28.05,28.1] 0.9551038 0.9486495 0.9609530
silicate + nitrate Indo-Pacific (28.1, Inf] 0.8041022 0.7748250 0.8522371
temp + aou Atlantic (-Inf,26] -0.0851759 -0.3981193 0.3090765
temp + aou Atlantic (26,26.5] -0.2067585 -0.3782715 -0.0333001
temp + aou Atlantic (26.5,26.75] -0.3439030 -0.4869077 -0.2278605
temp + aou Atlantic (26.75,27] -0.1902410 -0.2333278 -0.1663666
temp + aou Atlantic (27,27.25] -0.0945087 -0.2251903 0.0698207
temp + aou Atlantic (27.25,27.5] -0.1282161 -0.2779215 -0.0258653
temp + aou Atlantic (27.5,27.75] -0.3965518 -0.4385052 -0.3236885
temp + aou Atlantic (27.75,27.85] -0.5115646 -0.5307704 -0.4810074
temp + aou Atlantic (27.85,27.95] -0.3200279 -0.4537339 -0.1401074
temp + aou Atlantic (27.95,28.05] -0.5080419 -0.7118815 -0.3201445
temp + aou Atlantic (28.05,28.1] -0.5811548 -0.8030842 -0.4451464
temp + aou Atlantic (28.1,28.15] -0.6238241 -0.7216720 -0.5358002
temp + aou Atlantic (28.15,28.2] -0.6640638 -0.8202708 -0.4465653
temp + aou Atlantic (28.2, Inf] 0.4122638 0.3210534 0.4936508
temp + aou Indo-Pacific (-Inf,26] -0.2751309 -0.4113818 -0.1829156
temp + aou Indo-Pacific (26,26.5] -0.1146005 -0.2965725 -0.0148728
temp + aou Indo-Pacific (26.5,26.75] -0.0348987 -0.3357194 0.2515573
temp + aou Indo-Pacific (26.75,27] -0.1679451 -0.3472432 -0.0187972
temp + aou Indo-Pacific (27,27.25] -0.1526322 -0.2800714 0.0525694
temp + aou Indo-Pacific (27.25,27.5] 0.1001891 0.0002744 0.2363866
temp + aou Indo-Pacific (27.5,27.75] 0.3122538 0.1872311 0.5049224
temp + aou Indo-Pacific (27.75,27.85] -0.0086509 -0.1992955 0.1632033
temp + aou Indo-Pacific (27.85,27.95] -0.0303555 -0.2115065 0.0697391
temp + aou Indo-Pacific (27.95,28.05] 0.1290613 0.0869409 0.1677201
temp + aou Indo-Pacific (28.05,28.1] 0.2155092 0.0630482 0.3241956
temp + aou Indo-Pacific (28.1, Inf] 0.8481228 0.8021672 0.8741366
temp + nitrate Atlantic (-Inf,26] -0.1809098 -0.3100413 0.0622657
temp + nitrate Atlantic (26,26.5] -0.5477495 -0.6252393 -0.4962847
temp + nitrate Atlantic (26.5,26.75] -0.5973295 -0.6596867 -0.5285584
temp + nitrate Atlantic (26.75,27] -0.5228838 -0.5689871 -0.4862217
temp + nitrate Atlantic (27,27.25] -0.5580248 -0.6008066 -0.5128627
temp + nitrate Atlantic (27.25,27.5] -0.6855015 -0.7215558 -0.6342618
temp + nitrate Atlantic (27.5,27.75] -0.7346149 -0.7503145 -0.7187927
temp + nitrate Atlantic (27.75,27.85] -0.7020555 -0.7062854 -0.6983008
temp + nitrate Atlantic (27.85,27.95] -0.5229448 -0.6045090 -0.3600604
temp + nitrate Atlantic (27.95,28.05] -0.6378063 -0.7563390 -0.4554136
temp + nitrate Atlantic (28.05,28.1] -0.7070521 -0.8346832 -0.5635383
temp + nitrate Atlantic (28.1,28.15] -0.7129853 -0.7665529 -0.6108489
temp + nitrate Atlantic (28.15,28.2] -0.7118151 -0.7864857 -0.6243332
temp + nitrate Atlantic (28.2, Inf] 0.0681218 0.0092170 0.1142905
temp + nitrate Indo-Pacific (-Inf,26] -0.5291880 -0.6779130 -0.3799719
temp + nitrate Indo-Pacific (26,26.5] -0.4385356 -0.6243743 -0.3258660
temp + nitrate Indo-Pacific (26.5,26.75] -0.3481661 -0.6216810 -0.1084174
temp + nitrate Indo-Pacific (26.75,27] -0.4461526 -0.6039536 -0.3303970
temp + nitrate Indo-Pacific (27,27.25] -0.3913654 -0.4276749 -0.3230531
temp + nitrate Indo-Pacific (27.25,27.5] -0.1030011 -0.1859559 0.0119417
temp + nitrate Indo-Pacific (27.5,27.75] 0.2945481 -0.1168033 0.7296973
temp + nitrate Indo-Pacific (27.75,27.85] -0.1279082 -0.2183487 -0.0750038
temp + nitrate Indo-Pacific (27.85,27.95] -0.1450673 -0.2432755 -0.0783159
temp + nitrate Indo-Pacific (27.95,28.05] -0.0664517 -0.1486362 0.0046347
temp + nitrate Indo-Pacific (28.05,28.1] -0.0135400 -0.1764795 0.0744593
temp + nitrate Indo-Pacific (28.1, Inf] 0.4072338 0.2704676 0.5369527
temp + silicate Atlantic (-Inf,26] -0.6260580 -0.7368512 -0.5701898
temp + silicate Atlantic (26,26.5] -0.6784155 -0.6953405 -0.6685095
temp + silicate Atlantic (26.5,26.75] -0.6467624 -0.7140280 -0.5514632
temp + silicate Atlantic (26.75,27] -0.5447929 -0.5585652 -0.5359397
temp + silicate Atlantic (27,27.25] -0.6325754 -0.6560349 -0.6136864
temp + silicate Atlantic (27.25,27.5] -0.8293745 -0.8366160 -0.8182826
temp + silicate Atlantic (27.5,27.75] -0.8468240 -0.8709653 -0.8196601
temp + silicate Atlantic (27.75,27.85] -0.8078386 -0.8203330 -0.7908282
temp + silicate Atlantic (27.85,27.95] -0.6547381 -0.7267480 -0.5239783
temp + silicate Atlantic (27.95,28.05] -0.7263671 -0.8499067 -0.5716594
temp + silicate Atlantic (28.05,28.1] -0.7415263 -0.8780690 -0.6242234
temp + silicate Atlantic (28.1,28.15] -0.7456072 -0.8217856 -0.6641956
temp + silicate Atlantic (28.15,28.2] -0.6983869 -0.8073258 -0.5440499
temp + silicate Atlantic (28.2, Inf] 0.0909309 0.0584582 0.1328502
temp + silicate Indo-Pacific (-Inf,26] -0.5975045 -0.7261138 -0.4491996
temp + silicate Indo-Pacific (26,26.5] -0.6082548 -0.7122705 -0.4698433
temp + silicate Indo-Pacific (26.5,26.75] -0.6671407 -0.8019071 -0.5607259
temp + silicate Indo-Pacific (26.75,27] -0.6734037 -0.7577238 -0.5469896
temp + silicate Indo-Pacific (27,27.25] -0.5354191 -0.6210765 -0.4415797
temp + silicate Indo-Pacific (27.25,27.5] -0.2364291 -0.2696803 -0.1745648
temp + silicate Indo-Pacific (27.5,27.75] 0.0463766 -0.1387695 0.3060006
temp + silicate Indo-Pacific (27.75,27.85] -0.1288690 -0.2543698 -0.0170969
temp + silicate Indo-Pacific (27.85,27.95] -0.1636326 -0.3029008 -0.0895579
temp + silicate Indo-Pacific (27.95,28.05] -0.0686761 -0.1189552 -0.0261005
temp + silicate Indo-Pacific (28.05,28.1] 0.0721979 -0.0806949 0.1549224
temp + silicate Indo-Pacific (28.1, Inf] 0.2740956 0.0745190 0.3871172 
rm(cor_predictors, cor_predictors_stats)

2.2.3 C* vs individual predictors

Below, the range of correlations coefficients for C* with each predictor is plotted per basin (facet) and density slab (color). Note that the range indicates the min and max values of in total 3 calculated coefficients (one per era).

cor_target_predictor <- cor_target_predictor %>% 
  rename(predictor = term)

# calculate min, max, mean across all eras
cor_target_predictor_stats <- cor_target_predictor %>%
  select(-era) %>% 
  group_by(predictor, basin, gamma_slab) %>% 
  summarise_all(list(mean_r = mean, min_r = min, max_r = max)) %>% 
  ungroup()

# plot figure
cor_target_predictor_stats %>% 
  mutate(predictor = reorder(predictor, mean_r)) %>%
  ggplot() +
  geom_vline(xintercept = c(-0.9, 0.9), col = "red") +
  geom_vline(xintercept = 0) +
  geom_linerange(
    aes(y = predictor, xmin = min_r, xmax = max_r, col = gamma_slab),
    position = position_dodge(width = 0.6)) +
  facet_wrap(~basin) +
  scale_color_viridis_d(direction = -1) +
  labs(x = "correlation coefficient", y = "C* correlation with...") +
  theme(legend.position = "top")

Version Author Date
0aa2b50 jens-daniel-mueller 2020-12-23
2886da0 jens-daniel-mueller 2020-12-19
02f0ee9 jens-daniel-mueller 2020-12-18
7bcb4eb jens-daniel-mueller 2020-12-18
158fe26 jens-daniel-mueller 2020-12-15
7a9a4cb jens-daniel-mueller 2020-12-15
61b263c jens-daniel-mueller 2020-12-15
984697e jens-daniel-mueller 2020-12-12
3ebff89 jens-daniel-mueller 2020-12-12
5d96d3c jens-daniel-mueller 2020-12-11
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
0ff728b jens-daniel-mueller 2020-12-01
91435ae jens-daniel-mueller 2020-12-01
196be51 jens-daniel-mueller 2020-11-30
bc61ce3 Jens Müller 2020-11-30 
# print table
kable(cor_target_predictor_stats) %>%
  add_header_above() %>%
  kable_styling() %>%
  scroll_box(width = "100%", height = "400px")
predictor basin gamma_slab mean_r min_r max_r
aou Atlantic (-Inf,26] 0.3830181 0.1440842 0.6380857
aou Atlantic (26,26.5] -0.1245743 -0.3991308 0.0969375
aou Atlantic (26.5,26.75] -0.3948662 -0.4249988 -0.3773173
aou Atlantic (26.75,27] -0.2535462 -0.3437865 -0.1113619
aou Atlantic (27,27.25] -0.3074978 -0.3885808 -0.1925615
aou Atlantic (27.25,27.5] -0.4089852 -0.4851971 -0.3033817
aou Atlantic (27.5,27.75] -0.7634588 -0.7886942 -0.7244788
aou Atlantic (27.75,27.85] -0.9249006 -0.9315170 -0.9160734
aou Atlantic (27.85,27.95] -0.9342809 -0.9460473 -0.9280482
aou Atlantic (27.95,28.05] -0.8825705 -0.8890758 -0.8768636
aou Atlantic (28.05,28.1] -0.8444629 -0.8563246 -0.8254472
aou Atlantic (28.1,28.15] -0.8816436 -0.8880498 -0.8717577
aou Atlantic (28.15,28.2] -0.9218585 -0.9325341 -0.9042445
aou Atlantic (28.2, Inf] -0.1114958 -0.1728027 -0.0451601
aou Indo-Pacific (-Inf,26] -0.1627743 -0.3449584 -0.0017882
aou Indo-Pacific (26,26.5] -0.7168869 -0.7262033 -0.7101609
aou Indo-Pacific (26.5,26.75] -0.7319914 -0.7543229 -0.7115013
aou Indo-Pacific (26.75,27] -0.7116657 -0.8105528 -0.6131552
aou Indo-Pacific (27,27.25] 0.0174095 -0.0424298 0.0752846
aou Indo-Pacific (27.25,27.5] 0.5698873 0.5010153 0.6431101
aou Indo-Pacific (27.5,27.75] 0.7087383 0.6053759 0.7812682
aou Indo-Pacific (27.75,27.85] 0.7579569 0.6333551 0.8314406
aou Indo-Pacific (27.85,27.95] 0.7025995 0.5646737 0.7978928
aou Indo-Pacific (27.95,28.05] 0.4637161 0.2827374 0.6419535
aou Indo-Pacific (28.05,28.1] 0.2483325 0.1664652 0.3752323
aou Indo-Pacific (28.1, Inf] 0.3738607 0.2159182 0.6453168
nitrate Atlantic (-Inf,26] 0.1050073 -0.1017984 0.3701643
nitrate Atlantic (26,26.5] -0.4401480 -0.5929758 -0.3488691
nitrate Atlantic (26.5,26.75] -0.6193335 -0.6865261 -0.5395627
nitrate Atlantic (26.75,27] -0.5509072 -0.6573778 -0.4038497
nitrate Atlantic (27,27.25] -0.7103834 -0.7486744 -0.6700314
nitrate Atlantic (27.25,27.5] -0.8743415 -0.8948801 -0.8553441
nitrate Atlantic (27.5,27.75] -0.9706184 -0.9753393 -0.9618843
nitrate Atlantic (27.75,27.85] -0.9652244 -0.9730676 -0.9562539
nitrate Atlantic (27.85,27.95] -0.9385916 -0.9488667 -0.9304706
nitrate Atlantic (27.95,28.05] -0.9079981 -0.9323673 -0.8875973
nitrate Atlantic (28.05,28.1] -0.8876957 -0.9223522 -0.8662311
nitrate Atlantic (28.1,28.15] -0.9182939 -0.9404562 -0.9068616
nitrate Atlantic (28.15,28.2] -0.9574305 -0.9727023 -0.9448310
nitrate Atlantic (28.2, Inf] -0.7063784 -0.7516513 -0.6625676
nitrate Indo-Pacific (-Inf,26] -0.2047389 -0.3961106 -0.0445297
nitrate Indo-Pacific (26,26.5] -0.8096827 -0.8363954 -0.7649020
nitrate Indo-Pacific (26.5,26.75] -0.8436600 -0.8774341 -0.8181932
nitrate Indo-Pacific (26.75,27] -0.8235685 -0.8839218 -0.7345178
nitrate Indo-Pacific (27,27.25] -0.0425817 -0.0621112 -0.0232184
nitrate Indo-Pacific (27.25,27.5] 0.5107187 0.4474103 0.5841877
nitrate Indo-Pacific (27.5,27.75] 0.6716743 0.6125077 0.7285477
nitrate Indo-Pacific (27.75,27.85] 0.7402431 0.6587653 0.7934720
nitrate Indo-Pacific (27.85,27.95] 0.6707149 0.5646581 0.7607693
nitrate Indo-Pacific (27.95,28.05] 0.3947247 0.2152161 0.6022420
nitrate Indo-Pacific (28.05,28.1] 0.1599242 0.0665017 0.3399900
nitrate Indo-Pacific (28.1, Inf] 0.1445219 -0.0121114 0.4262703
sal Atlantic (-Inf,26] 0.2512223 0.1793040 0.3639105
sal Atlantic (26,26.5] 0.7445338 0.6336010 0.8292115
sal Atlantic (26.5,26.75] 0.8806279 0.8657431 0.9018502
sal Atlantic (26.75,27] 0.9202670 0.9018156 0.9343887
sal Atlantic (27,27.25] 0.9363668 0.9192656 0.9508385
sal Atlantic (27.25,27.5] 0.9287015 0.9047237 0.9535289
sal Atlantic (27.5,27.75] 0.8370623 0.8204456 0.8585839
sal Atlantic (27.75,27.85] 0.7105591 0.6933716 0.7221110
sal Atlantic (27.85,27.95] 0.4598247 0.3613666 0.5310005
sal Atlantic (27.95,28.05] 0.5066564 0.2265576 0.6841424
sal Atlantic (28.05,28.1] 0.6707803 0.6164069 0.7419700
sal Atlantic (28.1,28.15] 0.7542408 0.7112357 0.8088706
sal Atlantic (28.15,28.2] 0.8671350 0.8607639 0.8785950
sal Atlantic (28.2, Inf] 0.7020920 0.6453872 0.7364102
sal Indo-Pacific (-Inf,26] 0.0330003 -0.0341625 0.0717814
sal Indo-Pacific (26,26.5] 0.4436728 0.3271161 0.6572028
sal Indo-Pacific (26.5,26.75] 0.4884841 0.2880822 0.7408681
sal Indo-Pacific (26.75,27] 0.4569162 0.3511283 0.6682632
sal Indo-Pacific (27,27.25] -0.2323326 -0.4662042 -0.0965192
sal Indo-Pacific (27.25,27.5] -0.3504481 -0.5338893 -0.1766816
sal Indo-Pacific (27.5,27.75] -0.2902758 -0.3770012 -0.2283024
sal Indo-Pacific (27.75,27.85] -0.4178095 -0.5370460 -0.3212982
sal Indo-Pacific (27.85,27.95] -0.3996362 -0.5567612 -0.2418701
sal Indo-Pacific (27.95,28.05] -0.1101145 -0.1885055 0.0006139
sal Indo-Pacific (28.05,28.1] 0.0232885 -0.1033014 0.1632659
sal Indo-Pacific (28.1, Inf] 0.2893544 0.2504861 0.3431585
silicate Atlantic (-Inf,26] -0.2456739 -0.4269937 -0.1006930
silicate Atlantic (26,26.5] -0.6214511 -0.6892365 -0.5766186
silicate Atlantic (26.5,26.75] -0.6765171 -0.7149006 -0.6148891
silicate Atlantic (26.75,27] -0.5781199 -0.6466591 -0.4652937
silicate Atlantic (27,27.25] -0.7590046 -0.8043119 -0.7119611
silicate Atlantic (27.25,27.5] -0.8797179 -0.9028570 -0.8564637
silicate Atlantic (27.5,27.75] -0.8729908 -0.8909023 -0.8528956
silicate Atlantic (27.75,27.85] -0.8878451 -0.9201603 -0.8491464
silicate Atlantic (27.85,27.95] -0.8742003 -0.9029090 -0.8528623
silicate Atlantic (27.95,28.05] -0.8603167 -0.9050954 -0.8182550
silicate Atlantic (28.05,28.1] -0.8579624 -0.9023960 -0.8268018
silicate Atlantic (28.1,28.15] -0.9005775 -0.9266522 -0.8775091
silicate Atlantic (28.15,28.2] -0.9450009 -0.9514886 -0.9377657
silicate Atlantic (28.2, Inf] -0.6164544 -0.7600736 -0.5216039
silicate Indo-Pacific (-Inf,26] -0.2165038 -0.3871915 -0.0766448
silicate Indo-Pacific (26,26.5] -0.7321365 -0.7626486 -0.6953515
silicate Indo-Pacific (26.5,26.75] -0.7347999 -0.7969260 -0.6161383
silicate Indo-Pacific (26.75,27] -0.6551105 -0.7057277 -0.6246816
silicate Indo-Pacific (27,27.25] 0.2352792 0.2015347 0.2926645
silicate Indo-Pacific (27.25,27.5] 0.7645062 0.6981399 0.8011277
silicate Indo-Pacific (27.5,27.75] 0.8401187 0.7771956 0.8733525
silicate Indo-Pacific (27.75,27.85] 0.7932830 0.6950045 0.8567088
silicate Indo-Pacific (27.85,27.95] 0.6930884 0.5593271 0.7923990
silicate Indo-Pacific (27.95,28.05] 0.4315117 0.2433007 0.6287303
silicate Indo-Pacific (28.05,28.1] 0.1820392 0.1053170 0.3193432
silicate Indo-Pacific (28.1, Inf] 0.0397268 -0.0900561 0.2574856
temp Atlantic (-Inf,26] 0.0415334 -0.1571372 0.2013624
temp Atlantic (26,26.5] 0.6938081 0.5541504 0.7815567
temp Atlantic (26.5,26.75] 0.8689923 0.8401891 0.9067421
temp Atlantic (26.75,27] 0.9037761 0.8821850 0.9189779
temp Atlantic (27,27.25] 0.9133703 0.8971567 0.9260691
temp Atlantic (27.25,27.5] 0.9026622 0.8842077 0.9267237
temp Atlantic (27.5,27.75] 0.7858647 0.7747505 0.8052038
temp Atlantic (27.75,27.85] 0.6884914 0.6595455 0.7044793
temp Atlantic (27.85,27.95] 0.4278572 0.2865497 0.5442972
temp Atlantic (27.95,28.05] 0.4994601 0.3044851 0.6370820
temp Atlantic (28.05,28.1] 0.5517579 0.4410296 0.6470826
temp Atlantic (28.1,28.15] 0.5812781 0.5262238 0.6896696
temp Atlantic (28.15,28.2] 0.6021209 0.5747708 0.6248101
temp Atlantic (28.2, Inf] 0.5844192 0.5185888 0.6683384
temp Indo-Pacific (-Inf,26] -0.2581827 -0.3236845 -0.1595184
temp Indo-Pacific (26,26.5] 0.4549465 0.3218139 0.6371093
temp Indo-Pacific (26.5,26.75] 0.4810164 0.2932067 0.7209958
temp Indo-Pacific (26.75,27] 0.5034505 0.3891555 0.6920926
temp Indo-Pacific (27,27.25] -0.1204423 -0.3273641 0.0384837
temp Indo-Pacific (27.25,27.5] -0.3223411 -0.4669013 -0.1784236
temp Indo-Pacific (27.5,27.75] -0.3382037 -0.4377703 -0.2707134
temp Indo-Pacific (27.75,27.85] -0.3581289 -0.4758661 -0.2418659
temp Indo-Pacific (27.85,27.95] -0.2861855 -0.3936964 -0.1333148
temp Indo-Pacific (27.95,28.05] 0.0155185 -0.0678413 0.0589147
temp Indo-Pacific (28.05,28.1] 0.1462424 0.0764143 0.1857852
temp Indo-Pacific (28.1, Inf] 0.4207951 0.2699418 0.6957901 
rm(cor_target_predictor, cor_target_predictor_stats)

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

Matrix products: default
BLAS:   /usr/local/R-4.0.3/lib64/R/lib/libRblas.so
LAPACK: /usr/local/R-4.0.3/lib64/R/lib/libRlapack.so

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
 [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8    
 [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8   
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C            
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       

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

other attached packages:
 [1] corrr_0.4.3      broom_0.7.2      kableExtra_1.3.1 knitr_1.30      
 [5] olsrr_0.5.3.9000 GGally_2.0.0     lubridate_1.7.9  metR_0.9.0      
 [9] scico_1.2.0      patchwork_1.1.0  collapse_1.4.2   forcats_0.5.0   
[13] stringr_1.4.0    dplyr_1.0.2      purrr_0.3.4      readr_1.4.0     
[17] tidyr_1.1.2      tibble_3.0.4     ggplot2_3.3.2    tidyverse_1.3.0 
[21] workflowr_1.6.2 

loaded via a namespace (and not attached):
 [1] fs_1.5.0                 webshot_0.5.2            RColorBrewer_1.1-2      
 [4] httr_1.4.2               rprojroot_2.0.2          tools_4.0.3             
 [7] backports_1.1.10         R6_2.5.0                 nortest_1.0-4           
[10] DBI_1.1.0                colorspace_2.0-0         withr_2.3.0             
[13] gridExtra_2.3            tidyselect_1.1.0         curl_4.3                
[16] compiler_4.0.3           git2r_0.27.1             cli_2.2.0               
[19] rvest_0.3.6              xml2_1.3.2               labeling_0.4.2          
[22] scales_1.1.1             checkmate_2.0.0          goftest_1.2-2           
[25] digest_0.6.27            foreign_0.8-80           rmarkdown_2.5           
[28] rio_0.5.16               pkgconfig_2.0.3          htmltools_0.5.0         
[31] highr_0.8                dbplyr_1.4.4             rlang_0.4.9             
[34] readxl_1.3.1             rstudioapi_0.13          farver_2.0.3            
[37] generics_0.0.2           jsonlite_1.7.1           zip_2.1.1               
[40] car_3.0-10               magrittr_2.0.1           Matrix_1.2-18           
[43] Rcpp_1.0.5               munsell_0.5.0            fansi_0.4.1             
[46] abind_1.4-5              lifecycle_0.2.0          stringi_1.5.3           
[49] whisker_0.4              yaml_2.2.1               carData_3.0-4           
[52] plyr_1.8.6               grid_4.0.3               blob_1.2.1              
[55] parallel_4.0.3           promises_1.1.1           crayon_1.3.4            
[58] lattice_0.20-41          haven_2.3.1              hms_0.5.3               
[61] pillar_1.4.7             reprex_0.3.0             glue_1.4.2              
[64] evaluate_0.14            RcppArmadillo_0.10.1.2.0 data.table_1.13.2       
[67] modelr_0.1.8             vctrs_0.3.5              httpuv_1.5.4            
[70] cellranger_1.1.0         gtable_0.3.0             reshape_0.8.8           
[73] assertthat_0.2.1         xfun_0.18                openxlsx_4.2.3          
[76] RcppEigen_0.3.3.7.0      later_1.1.0.1            viridisLite_0.3.0       
[79] ellipsis_0.3.1           here_0.1