Last updated: 2021-01-24

Checks: 7 0

Knit directory: emlr_mod_v_XXX/

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

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: 846bfc0

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 846bfc0. 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/

Unstaged changes:
    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 a2f0d56 Donghe-Zhu 2021-01-23 Build site.
html 28509fc Donghe-Zhu 2021-01-23 Build site.
html 4c28e4a Donghe-Zhu 2021-01-22 Build site.
html 24cc264 jens-daniel-mueller 2021-01-22 cleaned /docs before creating copies
html 88eb28f Donghe-Zhu 2021-01-21 Build site.
html 2679490 Donghe-Zhu 2021-01-21 Build site.
html 7891955 Donghe-Zhu 2021-01-21 Build site.
html d4cf1cb Donghe-Zhu 2021-01-21 Build site.
html 1f3e5b6 jens-daniel-mueller 2021-01-20 Build site.
html 0e7bdf1 jens-daniel-mueller 2021-01-15 cleaning template repository
html 73cbef3 jens-daniel-mueller 2021-01-15 Build site.
html 4571843 jens-daniel-mueller 2021-01-14 revision and html deleted for template copying
html 23151cd jens-daniel-mueller 2021-01-14 Build site.
html b3564aa jens-daniel-mueller 2021-01-14 Build site.
html 8d032c3 jens-daniel-mueller 2021-01-14 Build site.
html 022871c Donghe-Zhu 2021-01-13 Build site.
Rmd d44f36f Donghe-Zhu 2021-01-13 reorder analysis final
html 17dee1d jens-daniel-mueller 2021-01-13 Build site.
html a076226 Donghe-Zhu 2021-01-11 Build site.
Rmd 52eff18 Donghe-Zhu 2021-01-09 Implemet model_run and subsetting
html 7cdea0c jens-daniel-mueller 2021-01-06 Build site.
html fa85b93 jens-daniel-mueller 2021-01-06 Build site.
html e5cb81a Donghe-Zhu 2021-01-05 Build site.
html a499f10 Donghe-Zhu 2021-01-05 Build site.
Rmd 715bdb4 Donghe-Zhu 2021-01-02 model modification
html fb8a752 Donghe-Zhu 2020-12-23 Build site.
Rmd 82e3c9c Donghe-Zhu 2020-12-23 first build after creating model template
html 8fae0b2 Donghe-Zhu 2020-12-21 Build site.
html c8b76b3 jens-daniel-mueller 2020-12-19 Build site.
Rmd b5fedce jens-daniel-mueller 2020-12-19 first build after creating model template
Rmd 8e8abf5 Jens Müller 2020-12-18 Initial commit

1 Required data

Required are:

  • cleaned and prepared GLODAP-based synthetic cmorized model subsetting 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 coefficients.

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
a2f0d56 Donghe-Zhu 2021-01-23
28509fc Donghe-Zhu 2021-01-23
4c28e4a Donghe-Zhu 2021-01-22
24cc264 jens-daniel-mueller 2021-01-22
7891955 Donghe-Zhu 2021-01-21
d4cf1cb Donghe-Zhu 2021-01-21
1f3e5b6 jens-daniel-mueller 2021-01-20
0e7bdf1 jens-daniel-mueller 2021-01-15
4571843 jens-daniel-mueller 2021-01-14
b3564aa jens-daniel-mueller 2021-01-14
8d032c3 jens-daniel-mueller 2021-01-14
17dee1d jens-daniel-mueller 2021-01-13
7cdea0c jens-daniel-mueller 2021-01-06
fa85b93 jens-daniel-mueller 2021-01-06
e5cb81a Donghe-Zhu 2021-01-05
a499f10 Donghe-Zhu 2021-01-05
fb8a752 Donghe-Zhu 2020-12-23
8fae0b2 Donghe-Zhu 2020-12-21
c8b76b3 jens-daniel-mueller 2020-12-19

Version Author Date
a2f0d56 Donghe-Zhu 2021-01-23
28509fc Donghe-Zhu 2021-01-23
4c28e4a Donghe-Zhu 2021-01-22
24cc264 jens-daniel-mueller 2021-01-22
7891955 Donghe-Zhu 2021-01-21
d4cf1cb Donghe-Zhu 2021-01-21
1f3e5b6 jens-daniel-mueller 2021-01-20
0e7bdf1 jens-daniel-mueller 2021-01-15
4571843 jens-daniel-mueller 2021-01-14
b3564aa jens-daniel-mueller 2021-01-14
8d032c3 jens-daniel-mueller 2021-01-14
17dee1d jens-daniel-mueller 2021-01-13
7cdea0c jens-daniel-mueller 2021-01-06
fa85b93 jens-daniel-mueller 2021-01-06
e5cb81a Donghe-Zhu 2021-01-05
a499f10 Donghe-Zhu 2021-01-05
fb8a752 Donghe-Zhu 2020-12-23
8fae0b2 Donghe-Zhu 2020-12-21
c8b76b3 jens-daniel-mueller 2020-12-19

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

/nfs/kryo/work/jenmueller/emlr_cant/model/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 assessment

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
a2f0d56 Donghe-Zhu 2021-01-23
28509fc Donghe-Zhu 2021-01-23
4c28e4a Donghe-Zhu 2021-01-22
24cc264 jens-daniel-mueller 2021-01-22
7891955 Donghe-Zhu 2021-01-21
d4cf1cb Donghe-Zhu 2021-01-21
1f3e5b6 jens-daniel-mueller 2021-01-20
0e7bdf1 jens-daniel-mueller 2021-01-15
4571843 jens-daniel-mueller 2021-01-14
b3564aa jens-daniel-mueller 2021-01-14
8d032c3 jens-daniel-mueller 2021-01-14
17dee1d jens-daniel-mueller 2021-01-13
e5cb81a Donghe-Zhu 2021-01-05
a499f10 Donghe-Zhu 2021-01-05
8fae0b2 Donghe-Zhu 2020-12-21
c8b76b3 jens-daniel-mueller 2020-12-19 
# 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 + phosphate Atlantic (-Inf,26] 0.9292515 0.9194106 0.9399366
aou + phosphate Atlantic (26,26.5] 0.8912726 0.8893331 0.8940324
aou + phosphate Atlantic (26.5,26.75] 0.9487873 0.9447682 0.9514802
aou + phosphate Atlantic (26.75,27] 0.9633632 0.9601375 0.9656255
aou + phosphate Atlantic (27,27.25] 0.9213875 0.9177527 0.9244745
aou + phosphate Atlantic (27.25,27.5] 0.7645666 0.7509747 0.7810054
aou + phosphate Atlantic (27.5,27.75] 0.4833183 0.4745936 0.4943400
aou + phosphate Atlantic (27.75,27.85] 0.2873566 0.2694339 0.3177601
aou + phosphate Atlantic (27.85,27.95] 0.5235515 0.4253902 0.6237919
aou + phosphate Atlantic (27.95,28.05] 0.1761117 0.1402641 0.2002315
aou + phosphate Atlantic (28.05,28.1] 0.4708000 0.4165720 0.5345300
aou + phosphate Atlantic (28.1,28.15] 0.8367515 0.8136552 0.8569727
aou + phosphate Atlantic (28.15,28.2] 0.9288133 0.9205191 0.9366171
aou + phosphate Atlantic (28.2, Inf] 0.8845363 0.8834007 0.8858820
aou + phosphate Indo-Pacific (-Inf,26] 0.9667103 0.9659406 0.9678654
aou + phosphate Indo-Pacific (26,26.5] 0.9592723 0.9583365 0.9602710
aou + phosphate Indo-Pacific (26.5,26.75] 0.9713336 0.9701722 0.9723817
aou + phosphate Indo-Pacific (26.75,27] 0.9728156 0.9725013 0.9732727
aou + phosphate Indo-Pacific (27,27.25] 0.9857593 0.9851347 0.9863322
aou + phosphate Indo-Pacific (27.25,27.5] 0.9866728 0.9861276 0.9870473
aou + phosphate Indo-Pacific (27.5,27.75] 0.9892080 0.9890966 0.9893364
aou + phosphate Indo-Pacific (27.75,27.85] 0.9944233 0.9941698 0.9945917
aou + phosphate Indo-Pacific (27.85,27.95] 0.9944533 0.9941251 0.9947299
aou + phosphate Indo-Pacific (27.95,28.05] 0.9867609 0.9863242 0.9870202
aou + phosphate Indo-Pacific (28.05,28.1] 0.9899491 0.9895381 0.9901763
aou + phosphate Indo-Pacific (28.1, Inf] 0.9892630 0.9882269 0.9900518
aou + phosphate_star Atlantic (-Inf,26] 0.2394828 0.1624787 0.2829947
aou + phosphate_star Atlantic (26,26.5] 0.3976127 0.3542750 0.4281844
aou + phosphate_star Atlantic (26.5,26.75] 0.1629442 0.1294478 0.1897935
aou + phosphate_star Atlantic (26.75,27] -0.1156820 -0.1279360 -0.0949271
aou + phosphate_star Atlantic (27,27.25] -0.0387544 -0.0627329 -0.0128964
aou + phosphate_star Atlantic (27.25,27.5] -0.1580372 -0.1734956 -0.1500956
aou + phosphate_star Atlantic (27.5,27.75] -0.3202287 -0.3357029 -0.3084762
aou + phosphate_star Atlantic (27.75,27.85] -0.3525517 -0.3923292 -0.2997408
aou + phosphate_star Atlantic (27.85,27.95] 0.1664595 0.0286166 0.3113898
aou + phosphate_star Atlantic (27.95,28.05] -0.2392009 -0.2696142 -0.2176909
aou + phosphate_star Atlantic (28.05,28.1] 0.1405614 0.0684256 0.2212936
aou + phosphate_star Atlantic (28.1,28.15] 0.6718686 0.6329504 0.7028218
aou + phosphate_star Atlantic (28.15,28.2] 0.8289053 0.8148361 0.8421747
aou + phosphate_star Atlantic (28.2, Inf] 0.7314839 0.7245233 0.7439888
aou + phosphate_star Indo-Pacific (-Inf,26] 0.5880253 0.5807937 0.5971065
aou + phosphate_star Indo-Pacific (26,26.5] 0.4176606 0.4135394 0.4201026
aou + phosphate_star Indo-Pacific (26.5,26.75] 0.4555587 0.4324046 0.4696718
aou + phosphate_star Indo-Pacific (26.75,27] 0.5935429 0.5825721 0.6065722
aou + phosphate_star Indo-Pacific (27,27.25] 0.7044325 0.6869792 0.7261065
aou + phosphate_star Indo-Pacific (27.25,27.5] 0.3289182 0.2669157 0.3753428
aou + phosphate_star Indo-Pacific (27.5,27.75] -0.3363708 -0.3976704 -0.2815824
aou + phosphate_star Indo-Pacific (27.75,27.85] -0.6478813 -0.6552659 -0.6377604
aou + phosphate_star Indo-Pacific (27.85,27.95] -0.4977512 -0.5207521 -0.4781334
aou + phosphate_star Indo-Pacific (27.95,28.05] -0.4947872 -0.5250019 -0.4699429
aou + phosphate_star Indo-Pacific (28.05,28.1] -0.4045952 -0.4703552 -0.3538409
aou + phosphate_star Indo-Pacific (28.1, Inf] 0.0431961 0.0165940 0.0576876
aou + silicate Atlantic (-Inf,26] 0.1177800 0.0428611 0.1985943
aou + silicate Atlantic (26,26.5] 0.2719832 0.2386089 0.2982113
aou + silicate Atlantic (26.5,26.75] 0.6146954 0.5880293 0.6407427
aou + silicate Atlantic (26.75,27] 0.8057186 0.7988995 0.8151206
aou + silicate Atlantic (27,27.25] 0.7493190 0.7441654 0.7544708
aou + silicate Atlantic (27.25,27.5] 0.5345390 0.5232626 0.5531612
aou + silicate Atlantic (27.5,27.75] 0.2049043 0.1991621 0.2117942
aou + silicate Atlantic (27.75,27.85] -0.0503226 -0.0713663 -0.0292156
aou + silicate Atlantic (27.85,27.95] 0.2902323 0.1731313 0.4132160
aou + silicate Atlantic (27.95,28.05] 0.1384738 0.1056334 0.1633169
aou + silicate Atlantic (28.05,28.1] 0.5162285 0.4742424 0.5701858
aou + silicate Atlantic (28.1,28.15] 0.8783083 0.8617460 0.8926386
aou + silicate Atlantic (28.15,28.2] 0.9466715 0.9418001 0.9509036
aou + silicate Atlantic (28.2, Inf] 0.6734468 0.6526341 0.6989391
aou + silicate Indo-Pacific (-Inf,26] 0.8361554 0.8329494 0.8385854
aou + silicate Indo-Pacific (26,26.5] 0.7718303 0.7713108 0.7724491
aou + silicate Indo-Pacific (26.5,26.75] 0.7817836 0.7773970 0.7887765
aou + silicate Indo-Pacific (26.75,27] 0.8720541 0.8701964 0.8731531
aou + silicate Indo-Pacific (27,27.25] 0.9442197 0.9437933 0.9447528
aou + silicate Indo-Pacific (27.25,27.5] 0.9578937 0.9574716 0.9583655
aou + silicate Indo-Pacific (27.5,27.75] 0.9683400 0.9676824 0.9688193
aou + silicate Indo-Pacific (27.75,27.85] 0.9844703 0.9836804 0.9850038
aou + silicate Indo-Pacific (27.85,27.95] 0.9846506 0.9840622 0.9850023
aou + silicate Indo-Pacific (27.95,28.05] 0.9835186 0.9829689 0.9838504
aou + silicate Indo-Pacific (28.05,28.1] 0.9886256 0.9875815 0.9892119
aou + silicate Indo-Pacific (28.1, Inf] 0.8363689 0.8338063 0.8383051
phosphate + phosphate_star Atlantic (-Inf,26] 0.5250520 0.4582282 0.5774083
phosphate + phosphate_star Atlantic (26,26.5] 0.7588389 0.7309221 0.7756182
phosphate + phosphate_star Atlantic (26.5,26.75] 0.4609493 0.4412509 0.4814148
phosphate + phosphate_star Atlantic (26.75,27] 0.1488164 0.1327271 0.1605721
phosphate + phosphate_star Atlantic (27,27.25] 0.3463175 0.3162404 0.3789905
phosphate + phosphate_star Atlantic (27.25,27.5] 0.5096226 0.4933008 0.5202769
phosphate + phosphate_star Atlantic (27.5,27.75] 0.6649111 0.6599274 0.6698035
phosphate + phosphate_star Atlantic (27.75,27.85] 0.7928840 0.7748347 0.8076772
phosphate + phosphate_star Atlantic (27.85,27.95] 0.9250857 0.9154900 0.9357191
phosphate + phosphate_star Atlantic (27.95,28.05] 0.9124146 0.9112701 0.9144645
phosphate + phosphate_star Atlantic (28.05,28.1] 0.9385407 0.9346131 0.9415133
phosphate + phosphate_star Atlantic (28.1,28.15] 0.9670449 0.9644274 0.9684368
phosphate + phosphate_star Atlantic (28.15,28.2] 0.9765000 0.9757871 0.9771408
phosphate + phosphate_star Atlantic (28.2, Inf] 0.9609148 0.9585296 0.9643142
phosphate + phosphate_star Indo-Pacific (-Inf,26] 0.7602704 0.7565811 0.7644660
phosphate + phosphate_star Indo-Pacific (26,26.5] 0.6519482 0.6488668 0.6555652
phosphate + phosphate_star Indo-Pacific (26.5,26.75] 0.6469652 0.6316892 0.6547737
phosphate + phosphate_star Indo-Pacific (26.75,27] 0.7511701 0.7421574 0.7595747
phosphate + phosphate_star Indo-Pacific (27,27.25] 0.8086155 0.7966290 0.8239674
phosphate + phosphate_star Indo-Pacific (27.25,27.5] 0.4748480 0.4192468 0.5175535
phosphate + phosphate_star Indo-Pacific (27.5,27.75] -0.1973377 -0.2604249 -0.1413155
phosphate + phosphate_star Indo-Pacific (27.75,27.85] -0.5669797 -0.5729218 -0.5565211
phosphate + phosphate_star Indo-Pacific (27.85,27.95] -0.4078420 -0.4322318 -0.3896299
phosphate + phosphate_star Indo-Pacific (27.95,28.05] -0.3512979 -0.3812845 -0.3259960
phosphate + phosphate_star Indo-Pacific (28.05,28.1] -0.2747125 -0.3413173 -0.2227075
phosphate + phosphate_star Indo-Pacific (28.1, Inf] 0.1835974 0.1518542 0.2049764
sal + aou Atlantic (-Inf,26] -0.4062646 -0.4569132 -0.3689498
sal + aou Atlantic (26,26.5] -0.3776050 -0.4236720 -0.3414765
sal + aou Atlantic (26.5,26.75] 0.2610390 0.2482607 0.2831258
sal + aou Atlantic (26.75,27] 0.4513904 0.4388865 0.4598370
sal + aou Atlantic (27,27.25] 0.3898990 0.3531036 0.4272745
sal + aou Atlantic (27.25,27.5] 0.4297204 0.4207392 0.4447274
sal + aou Atlantic (27.5,27.75] 0.6826882 0.6744922 0.6957271
sal + aou Atlantic (27.75,27.85] 0.6833351 0.6314365 0.7159754
sal + aou Atlantic (27.85,27.95] 0.2311877 0.0960043 0.3666390
sal + aou Atlantic (27.95,28.05] 0.2158888 0.1931916 0.2453183
sal + aou Atlantic (28.05,28.1] -0.1554728 -0.2272560 -0.0888771
sal + aou Atlantic (28.1,28.15] -0.6713629 -0.6998516 -0.6355779
sal + aou Atlantic (28.15,28.2] -0.7965897 -0.8039091 -0.7836805
sal + aou Atlantic (28.2, Inf] -0.5405990 -0.5637477 -0.5150389
sal + aou Indo-Pacific (-Inf,26] -0.2331734 -0.2382909 -0.2297198
sal + aou Indo-Pacific (26,26.5] -0.1600378 -0.1625333 -0.1554548
sal + aou Indo-Pacific (26.5,26.75] -0.1840466 -0.1889324 -0.1753521
sal + aou Indo-Pacific (26.75,27] -0.4161478 -0.4229152 -0.4030782
sal + aou Indo-Pacific (27,27.25] -0.6268290 -0.6455661 -0.6128213
sal + aou Indo-Pacific (27.25,27.5] -0.2994425 -0.3486204 -0.2463665
sal + aou Indo-Pacific (27.5,27.75] 0.5174474 0.4739132 0.5762683
sal + aou Indo-Pacific (27.75,27.85] 0.7988859 0.7942414 0.8038401
sal + aou Indo-Pacific (27.85,27.95] 0.8144371 0.8120544 0.8161352
sal + aou Indo-Pacific (27.95,28.05] 0.8638350 0.8547387 0.8734655
sal + aou Indo-Pacific (28.05,28.1] 0.8822113 0.8696951 0.8918014
sal + aou Indo-Pacific (28.1, Inf] 0.4363003 0.4208270 0.4502529
sal + phosphate Atlantic (-Inf,26] -0.2911109 -0.3888854 -0.2351049
sal + phosphate Atlantic (26,26.5] -0.6415071 -0.6619704 -0.6136488
sal + phosphate Atlantic (26.5,26.75] 0.0130512 0.0005468 0.0334390
sal + phosphate Atlantic (26.75,27] 0.2815793 0.2708340 0.2971710
sal + phosphate Atlantic (27,27.25] 0.1089767 0.0560839 0.1651019
sal + phosphate Atlantic (27.25,27.5] -0.1238997 -0.1728969 -0.0617804
sal + phosphate Atlantic (27.5,27.75] 0.1264645 0.0886744 0.1740041
sal + phosphate Atlantic (27.75,27.85] -0.4203209 -0.4607403 -0.3887917
sal + phosphate Atlantic (27.85,27.95] -0.6392639 -0.6642959 -0.6243918
sal + phosphate Atlantic (27.95,28.05] -0.8570723 -0.8633370 -0.8481127
sal + phosphate Atlantic (28.05,28.1] -0.9197810 -0.9222922 -0.9181476
sal + phosphate Atlantic (28.1,28.15] -0.9557962 -0.9583677 -0.9543802
sal + phosphate Atlantic (28.15,28.2] -0.9460025 -0.9487662 -0.9428381
sal + phosphate Atlantic (28.2, Inf] -0.6991414 -0.7060614 -0.6950667
sal + phosphate Indo-Pacific (-Inf,26] -0.3535559 -0.3558418 -0.3499962
sal + phosphate Indo-Pacific (26,26.5] -0.3782104 -0.3830931 -0.3706121
sal + phosphate Indo-Pacific (26.5,26.75] -0.3629457 -0.3662674 -0.3585739
sal + phosphate Indo-Pacific (26.75,27] -0.5310904 -0.5483642 -0.5154593
sal + phosphate Indo-Pacific (27,27.25] -0.7053193 -0.7164621 -0.6963493
sal + phosphate Indo-Pacific (27.25,27.5] -0.3993491 -0.4490935 -0.3442137
sal + phosphate Indo-Pacific (27.5,27.75] 0.4195491 0.3749865 0.4793682
sal + phosphate Indo-Pacific (27.75,27.85] 0.7489016 0.7407569 0.7566652
sal + phosphate Indo-Pacific (27.85,27.95] 0.7743361 0.7727415 0.7756260
sal + phosphate Indo-Pacific (27.95,28.05] 0.8147132 0.8040411 0.8254605
sal + phosphate Indo-Pacific (28.05,28.1] 0.8380885 0.8226606 0.8512550
sal + phosphate Indo-Pacific (28.1, Inf] 0.4458058 0.4314369 0.4581853
sal + phosphate_star Atlantic (-Inf,26] -0.2430926 -0.3278989 -0.1505424
sal + phosphate_star Atlantic (26,26.5] -0.8757327 -0.8923127 -0.8670023
sal + phosphate_star Atlantic (26.5,26.75] -0.7839335 -0.8038140 -0.7735610
sal + phosphate_star Atlantic (26.75,27] -0.7433979 -0.7666103 -0.7264066
sal + phosphate_star Atlantic (27,27.25] -0.7569740 -0.7746242 -0.7327742
sal + phosphate_star Atlantic (27.25,27.5] -0.8238016 -0.8534970 -0.7910246
sal + phosphate_star Atlantic (27.5,27.75] -0.5298907 -0.5641172 -0.4935860
sal + phosphate_star Atlantic (27.75,27.85] -0.8603089 -0.8704549 -0.8535605
sal + phosphate_star Atlantic (27.85,27.95] -0.8603273 -0.8706289 -0.8440294
sal + phosphate_star Atlantic (27.95,28.05] -0.9485600 -0.9536619 -0.9426898
sal + phosphate_star Atlantic (28.05,28.1] -0.9800714 -0.9805493 -0.9791847
sal + phosphate_star Atlantic (28.1,28.15] -0.9890491 -0.9911631 -0.9867078
sal + phosphate_star Atlantic (28.15,28.2] -0.9788088 -0.9807036 -0.9759401
sal + phosphate_star Atlantic (28.2, Inf] -0.7796586 -0.7866398 -0.7676352
sal + phosphate_star Indo-Pacific (-Inf,26] -0.7110137 -0.7155861 -0.7075781
sal + phosphate_star Indo-Pacific (26,26.5] -0.8690164 -0.8758677 -0.8647775
sal + phosphate_star Indo-Pacific (26.5,26.75] -0.8797529 -0.8912219 -0.8718460
sal + phosphate_star Indo-Pacific (26.75,27] -0.8580166 -0.8763250 -0.8488207
sal + phosphate_star Indo-Pacific (27,27.25] -0.9066968 -0.9100266 -0.9032651
sal + phosphate_star Indo-Pacific (27.25,27.5] -0.8072122 -0.8239302 -0.7797753
sal + phosphate_star Indo-Pacific (27.5,27.75] -0.8208790 -0.8351081 -0.8114991
sal + phosphate_star Indo-Pacific (27.75,27.85] -0.9065033 -0.9181091 -0.8996453
sal + phosphate_star Indo-Pacific (27.85,27.95] -0.7818333 -0.8027658 -0.7683643
sal + phosphate_star Indo-Pacific (27.95,28.05] -0.6961295 -0.7067002 -0.6885578
sal + phosphate_star Indo-Pacific (28.05,28.1] -0.6512503 -0.6744568 -0.6344230
sal + phosphate_star Indo-Pacific (28.1, Inf] -0.0060324 -0.0273991 0.0077986
sal + silicate Atlantic (-Inf,26] -0.7666791 -0.7830523 -0.7520002
sal + silicate Atlantic (26,26.5] 0.1487404 0.1172642 0.1710247
sal + silicate Atlantic (26.5,26.75] 0.3009731 0.2777178 0.3435170
sal + silicate Atlantic (26.75,27] 0.2050417 0.1942691 0.2263134
sal + silicate Atlantic (27,27.25] -0.0317965 -0.0844415 0.0293438
sal + silicate Atlantic (27.25,27.5] -0.2705320 -0.3189942 -0.2052475
sal + silicate Atlantic (27.5,27.75] -0.0844472 -0.1314136 -0.0360830
sal + silicate Atlantic (27.75,27.85] -0.6714757 -0.7067732 -0.6395989
sal + silicate Atlantic (27.85,27.95] -0.7944241 -0.8153488 -0.7679733
sal + silicate Atlantic (27.95,28.05] -0.8944690 -0.8983794 -0.8879716
sal + silicate Atlantic (28.05,28.1] -0.9056015 -0.9097075 -0.8980791
sal + silicate Atlantic (28.1,28.15] -0.9216547 -0.9260961 -0.9178318
sal + silicate Atlantic (28.15,28.2] -0.7961999 -0.8044057 -0.7890895
sal + silicate Atlantic (28.2, Inf] -0.1076821 -0.1117152 -0.1022324
sal + silicate Indo-Pacific (-Inf,26] -0.5419447 -0.5467874 -0.5333976
sal + silicate Indo-Pacific (26,26.5] -0.5840050 -0.5893446 -0.5798730
sal + silicate Indo-Pacific (26.5,26.75] -0.6111466 -0.6174231 -0.6047020
sal + silicate Indo-Pacific (26.75,27] -0.6779054 -0.6898384 -0.6666455
sal + silicate Indo-Pacific (27,27.25] -0.7740539 -0.7842916 -0.7640305
sal + silicate Indo-Pacific (27.25,27.5] -0.4965759 -0.5407315 -0.4446429
sal + silicate Indo-Pacific (27.5,27.75] 0.3570432 0.3095308 0.4206851
sal + silicate Indo-Pacific (27.75,27.85] 0.7269077 0.7170648 0.7361142
sal + silicate Indo-Pacific (27.85,27.95] 0.7702858 0.7681153 0.7738050
sal + silicate Indo-Pacific (27.95,28.05] 0.8291528 0.8172643 0.8405078
sal + silicate Indo-Pacific (28.05,28.1] 0.8549575 0.8415226 0.8651257
sal + silicate Indo-Pacific (28.1, Inf] 0.7352861 0.7218278 0.7465952
sal + temp Atlantic (-Inf,26] 0.6604365 0.6490670 0.6796755
sal + temp Atlantic (26,26.5] 0.9581285 0.9550581 0.9602678
sal + temp Atlantic (26.5,26.75] 0.9696455 0.9677213 0.9730188
sal + temp Atlantic (26.75,27] 0.9722644 0.9697470 0.9739595
sal + temp Atlantic (27,27.25] 0.9657251 0.9653764 0.9660299
sal + temp Atlantic (27.25,27.5] 0.9738815 0.9727263 0.9750673
sal + temp Atlantic (27.5,27.75] 0.9675060 0.9667786 0.9685357
sal + temp Atlantic (27.75,27.85] 0.9868502 0.9864795 0.9870834
sal + temp Atlantic (27.85,27.95] 0.9808757 0.9764118 0.9833204
sal + temp Atlantic (27.95,28.05] 0.9783842 0.9775270 0.9792728
sal + temp Atlantic (28.05,28.1] 0.9878533 0.9872414 0.9884037
sal + temp Atlantic (28.1,28.15] 0.9852256 0.9849785 0.9855916
sal + temp Atlantic (28.15,28.2] 0.9842337 0.9841219 0.9843017
sal + temp Atlantic (28.2, Inf] 0.9165717 0.9088456 0.9215384
sal + temp Indo-Pacific (-Inf,26] 0.8222864 0.8201284 0.8246266
sal + temp Indo-Pacific (26,26.5] 0.9550397 0.9533648 0.9562429
sal + temp Indo-Pacific (26.5,26.75] 0.9809594 0.9807897 0.9810751
sal + temp Indo-Pacific (26.75,27] 0.9807169 0.9797191 0.9815621
sal + temp Indo-Pacific (27,27.25] 0.9774531 0.9768436 0.9782897
sal + temp Indo-Pacific (27.25,27.5] 0.9219953 0.9202496 0.9231127
sal + temp Indo-Pacific (27.5,27.75] 0.8751626 0.8729929 0.8762932
sal + temp Indo-Pacific (27.75,27.85] 0.9725914 0.9709843 0.9743539
sal + temp Indo-Pacific (27.85,27.95] 0.9515315 0.9494284 0.9542711
sal + temp Indo-Pacific (27.95,28.05] 0.9479928 0.9472658 0.9486306
sal + temp Indo-Pacific (28.05,28.1] 0.9772837 0.9765829 0.9783311
sal + temp Indo-Pacific (28.1, Inf] 0.5287939 0.5131957 0.5414809
silicate + phosphate Atlantic (-Inf,26] -0.1108891 -0.2206780 0.0298299
silicate + phosphate Atlantic (26,26.5] 0.0365889 0.0024170 0.0775515
silicate + phosphate Atlantic (26.5,26.75] 0.5612932 0.5119232 0.6052302
silicate + phosphate Atlantic (26.75,27] 0.8472029 0.8421307 0.8528946
silicate + phosphate Atlantic (27,27.25] 0.9079993 0.9050683 0.9103604
silicate + phosphate Atlantic (27.25,27.5] 0.9087311 0.9065501 0.9103433
silicate + phosphate Atlantic (27.5,27.75] 0.9145860 0.9121076 0.9168095
silicate + phosphate Atlantic (27.75,27.85] 0.8956982 0.8945190 0.8974453
silicate + phosphate Atlantic (27.85,27.95] 0.9399793 0.9347996 0.9464344
silicate + phosphate Atlantic (27.95,28.05] 0.9769601 0.9764080 0.9774905
silicate + phosphate Atlantic (28.05,28.1] 0.9857783 0.9850959 0.9864724
silicate + phosphate Atlantic (28.1,28.15] 0.9879030 0.9875228 0.9883290
silicate + phosphate Atlantic (28.15,28.2] 0.9434590 0.9423331 0.9451416
silicate + phosphate Atlantic (28.2, Inf] 0.7602997 0.7531751 0.7665049
silicate + phosphate Indo-Pacific (-Inf,26] 0.8857214 0.8835508 0.8872260
silicate + phosphate Indo-Pacific (26,26.5] 0.8788954 0.8765037 0.8803575
silicate + phosphate Indo-Pacific (26.5,26.75] 0.8668316 0.8639792 0.8687257
silicate + phosphate Indo-Pacific (26.75,27] 0.9228419 0.9222433 0.9232013
silicate + phosphate Indo-Pacific (27,27.25] 0.9601517 0.9591100 0.9607539
silicate + phosphate Indo-Pacific (27.25,27.5] 0.9742334 0.9739650 0.9747515
silicate + phosphate Indo-Pacific (27.5,27.75] 0.9828278 0.9825550 0.9832300
silicate + phosphate Indo-Pacific (27.75,27.85] 0.9900289 0.9890934 0.9907512
silicate + phosphate Indo-Pacific (27.85,27.95] 0.9901705 0.9897852 0.9907545
silicate + phosphate Indo-Pacific (27.95,28.05] 0.9942249 0.9940157 0.9945228
silicate + phosphate Indo-Pacific (28.05,28.1] 0.9949536 0.9947761 0.9950677
silicate + phosphate Indo-Pacific (28.1, Inf] 0.8763306 0.8758965 0.8766181
silicate + phosphate_star Atlantic (-Inf,26] -0.1622457 -0.2465395 -0.0296334
silicate + phosphate_star Atlantic (26,26.5] -0.2772006 -0.3154765 -0.2309406
silicate + phosphate_star Atlantic (26.5,26.75] 0.0171249 -0.0689974 0.0738089
silicate + phosphate_star Atlantic (26.75,27] 0.1777978 0.1590836 0.1961224
silicate + phosphate_star Atlantic (27,27.25] 0.5195742 0.4863158 0.5508853
silicate + phosphate_star Atlantic (27.25,27.5] 0.6637930 0.6531089 0.6741270
silicate + phosphate_star Atlantic (27.5,27.75] 0.8104296 0.8014265 0.8212999
silicate + phosphate_star Atlantic (27.75,27.85] 0.9062917 0.8887847 0.9202127
silicate + phosphate_star Atlantic (27.85,27.95] 0.9626592 0.9587300 0.9681482
silicate + phosphate_star Atlantic (27.95,28.05] 0.9109425 0.9101774 0.9113690
silicate + phosphate_star Atlantic (28.05,28.1] 0.9071439 0.8982928 0.9130375
silicate + phosphate_star Atlantic (28.1,28.15] 0.9313512 0.9262375 0.9340488
silicate + phosphate_star Atlantic (28.15,28.2] 0.8737381 0.8708278 0.8763922
silicate + phosphate_star Atlantic (28.2, Inf] 0.6782217 0.6708205 0.6847330
silicate + phosphate_star Indo-Pacific (-Inf,26] 0.7717461 0.7693337 0.7756372
silicate + phosphate_star Indo-Pacific (26,26.5] 0.7766637 0.7759753 0.7771998
silicate + phosphate_star Indo-Pacific (26.5,26.75] 0.7897679 0.7877944 0.7928896
silicate + phosphate_star Indo-Pacific (26.75,27] 0.8171786 0.8086671 0.8250998
silicate + phosphate_star Indo-Pacific (27,27.25] 0.8209026 0.8031250 0.8395208
silicate + phosphate_star Indo-Pacific (27.25,27.5] 0.5092828 0.4547594 0.5500068
silicate + phosphate_star Indo-Pacific (27.5,27.75] -0.1545258 -0.2204315 -0.0947077
silicate + phosphate_star Indo-Pacific (27.75,27.85] -0.5542522 -0.5617871 -0.5426921
silicate + phosphate_star Indo-Pacific (27.85,27.95] -0.4011510 -0.4223966 -0.3865736
silicate + phosphate_star Indo-Pacific (27.95,28.05] -0.3688805 -0.3968661 -0.3438888
silicate + phosphate_star Indo-Pacific (28.05,28.1] -0.3040564 -0.3642864 -0.2562873
silicate + phosphate_star Indo-Pacific (28.1, Inf] 0.3255520 0.2934792 0.3496408
temp + aou Atlantic (-Inf,26] -0.6706962 -0.7046934 -0.6397576
temp + aou Atlantic (26,26.5] -0.4223480 -0.4511355 -0.3974631
temp + aou Atlantic (26.5,26.75] 0.1502000 0.1340268 0.1704044
temp + aou Atlantic (26.75,27] 0.3928294 0.3769225 0.4057523
temp + aou Atlantic (27,27.25] 0.3632434 0.3302350 0.3960109
temp + aou Atlantic (27.25,27.5] 0.4531016 0.4427140 0.4683424
temp + aou Atlantic (27.5,27.75] 0.6681250 0.6615004 0.6784605
temp + aou Atlantic (27.75,27.85] 0.6887640 0.6425732 0.7191718
temp + aou Atlantic (27.85,27.95] 0.2303808 0.0963479 0.3541226
temp + aou Atlantic (27.95,28.05] 0.3264957 0.3045961 0.3551246
temp + aou Atlantic (28.05,28.1] -0.0851063 -0.1594271 -0.0123631
temp + aou Atlantic (28.1,28.15] -0.6079820 -0.6362557 -0.5681850
temp + aou Atlantic (28.15,28.2] -0.7510685 -0.7625419 -0.7361302
temp + aou Atlantic (28.2, Inf] -0.6364416 -0.6461099 -0.6192575
temp + aou Indo-Pacific (-Inf,26] -0.3770208 -0.3786895 -0.3742847
temp + aou Indo-Pacific (26,26.5] -0.2021793 -0.2054421 -0.1980450
temp + aou Indo-Pacific (26.5,26.75] -0.1623627 -0.1687275 -0.1504496
temp + aou Indo-Pacific (26.75,27] -0.3616401 -0.3685472 -0.3494195
temp + aou Indo-Pacific (27,27.25] -0.5900263 -0.6079131 -0.5786702
temp + aou Indo-Pacific (27.25,27.5] -0.2592334 -0.2958664 -0.2145384
temp + aou Indo-Pacific (27.5,27.75] 0.4256352 0.3856943 0.4762032
temp + aou Indo-Pacific (27.75,27.85] 0.7857201 0.7815670 0.7882150
temp + aou Indo-Pacific (27.85,27.95] 0.7844426 0.7815618 0.7885802
temp + aou Indo-Pacific (27.95,28.05] 0.8194754 0.8117819 0.8295829
temp + aou Indo-Pacific (28.05,28.1] 0.8571042 0.8425968 0.8719163
temp + aou Indo-Pacific (28.1, Inf] 0.8156381 0.8060443 0.8249069
temp + phosphate Atlantic (-Inf,26] -0.6888977 -0.7521368 -0.6452803
temp + phosphate Atlantic (26,26.5] -0.6592327 -0.6685501 -0.6407175
temp + phosphate Atlantic (26.5,26.75] -0.0884151 -0.1000796 -0.0685467
temp + phosphate Atlantic (26.75,27] 0.2144376 0.2082332 0.2266958
temp + phosphate Atlantic (27,27.25] 0.0893949 0.0391779 0.1406970
temp + phosphate Atlantic (27.25,27.5] -0.0847624 -0.1316138 -0.0212447
temp + phosphate Atlantic (27.5,27.75] 0.0941057 0.0647368 0.1309799
temp + phosphate Atlantic (27.75,27.85] -0.4246572 -0.4569037 -0.3969581
temp + phosphate Atlantic (27.85,27.95] -0.6589005 -0.6764810 -0.6493754
temp + phosphate Atlantic (27.95,28.05] -0.8265632 -0.8304338 -0.8210796
temp + phosphate Atlantic (28.05,28.1] -0.8866884 -0.8918404 -0.8822770
temp + phosphate Atlantic (28.1,28.15] -0.9227688 -0.9243513 -0.9203422
temp + phosphate Atlantic (28.15,28.2] -0.9269409 -0.9289907 -0.9236192
temp + phosphate Atlantic (28.2, Inf] -0.8308890 -0.8492777 -0.8128560
temp + phosphate Indo-Pacific (-Inf,26] -0.4993205 -0.5037008 -0.4968231
temp + phosphate Indo-Pacific (26,26.5] -0.4262892 -0.4277956 -0.4246234
temp + phosphate Indo-Pacific (26.5,26.75] -0.3442073 -0.3487591 -0.3415577
temp + phosphate Indo-Pacific (26.75,27] -0.4862052 -0.5022879 -0.4710855
temp + phosphate Indo-Pacific (27,27.25] -0.6782463 -0.6879504 -0.6718818
temp + phosphate Indo-Pacific (27.25,27.5] -0.3827936 -0.4206433 -0.3366120
temp + phosphate Indo-Pacific (27.5,27.75] 0.3063070 0.2660524 0.3577065
temp + phosphate Indo-Pacific (27.75,27.85] 0.7323028 0.7276568 0.7374813
temp + phosphate Indo-Pacific (27.85,27.95] 0.7351505 0.7331591 0.7385016
temp + phosphate Indo-Pacific (27.95,28.05] 0.7492193 0.7402004 0.7597818
temp + phosphate Indo-Pacific (28.05,28.1] 0.8031850 0.7861968 0.8214039
temp + phosphate Indo-Pacific (28.1, Inf] 0.7667281 0.7508399 0.7815563
temp + phosphate_star Atlantic (-Inf,26] -0.6779711 -0.7302908 -0.6230139
temp + phosphate_star Atlantic (26,26.5] -0.8646134 -0.8784317 -0.8574674
temp + phosphate_star Atlantic (26.5,26.75] -0.8045621 -0.8240155 -0.7911259
temp + phosphate_star Atlantic (26.75,27] -0.7823216 -0.8028899 -0.7713480
temp + phosphate_star Atlantic (27,27.25] -0.7555835 -0.7762358 -0.7316394
temp + phosphate_star Atlantic (27.25,27.5] -0.8034838 -0.8353609 -0.7686559
temp + phosphate_star Atlantic (27.5,27.75] -0.5587843 -0.5875648 -0.5311938
temp + phosphate_star Atlantic (27.75,27.85] -0.8690431 -0.8748011 -0.8643248
temp + phosphate_star Atlantic (27.85,27.95] -0.8826153 -0.8900096 -0.8761503
temp + phosphate_star Atlantic (27.95,28.05] -0.9633930 -0.9648657 -0.9607528
temp + phosphate_star Atlantic (28.05,28.1] -0.9733564 -0.9748886 -0.9712625
temp + phosphate_star Atlantic (28.1,28.15] -0.9779750 -0.9786070 -0.9768462
temp + phosphate_star Atlantic (28.15,28.2] -0.9779772 -0.9794525 -0.9754261
temp + phosphate_star Atlantic (28.2, Inf] -0.9086166 -0.9237497 -0.8870265
temp + phosphate_star Indo-Pacific (-Inf,26] -0.8467833 -0.8503010 -0.8442997
temp + phosphate_star Indo-Pacific (26,26.5] -0.9137349 -0.9167031 -0.9111623
temp + phosphate_star Indo-Pacific (26.5,26.75] -0.8851309 -0.8959836 -0.8764429
temp + phosphate_star Indo-Pacific (26.75,27] -0.8613154 -0.8788271 -0.8521545
temp + phosphate_star Indo-Pacific (27,27.25] -0.9222417 -0.9265320 -0.9157792
temp + phosphate_star Indo-Pacific (27.25,27.5] -0.9159857 -0.9279149 -0.8973211
temp + phosphate_star Indo-Pacific (27.5,27.75] -0.9470870 -0.9538829 -0.9408860
temp + phosphate_star Indo-Pacific (27.75,27.85] -0.9341918 -0.9396733 -0.9298138
temp + phosphate_star Indo-Pacific (27.85,27.95] -0.8591334 -0.8734375 -0.8507479
temp + phosphate_star Indo-Pacific (27.95,28.05] -0.8047010 -0.8134910 -0.7989997
temp + phosphate_star Indo-Pacific (28.05,28.1] -0.7151794 -0.7349195 -0.6988595
temp + phosphate_star Indo-Pacific (28.1, Inf] -0.3611979 -0.3685599 -0.3514506
temp + silicate Atlantic (-Inf,26] -0.3651392 -0.3884488 -0.3529087
temp + silicate Atlantic (26,26.5] 0.1118548 0.0987413 0.1315408
temp + silicate Atlantic (26.5,26.75] 0.2050847 0.1819231 0.2505844
temp + silicate Atlantic (26.75,27] 0.1019821 0.0932263 0.1183748
temp + silicate Atlantic (27,27.25] -0.0781814 -0.1264101 -0.0235772
temp + silicate Atlantic (27.25,27.5] -0.2324566 -0.2805431 -0.1650125
temp + silicate Atlantic (27.5,27.75] -0.1267625 -0.1651374 -0.0883813
temp + silicate Atlantic (27.75,27.85] -0.6657131 -0.6950276 -0.6355765
temp + silicate Atlantic (27.85,27.95] -0.8039677 -0.8142140 -0.7907236
temp + silicate Atlantic (27.95,28.05] -0.8594938 -0.8621026 -0.8555638
temp + silicate Atlantic (28.05,28.1] -0.8655830 -0.8728005 -0.8532177
temp + silicate Atlantic (28.1,28.15] -0.8821692 -0.8857513 -0.8761624
temp + silicate Atlantic (28.15,28.2] -0.7813143 -0.7873503 -0.7746203
temp + silicate Atlantic (28.2, Inf] -0.3953665 -0.4063661 -0.3896145
temp + silicate Indo-Pacific (-Inf,26] -0.6017360 -0.6086808 -0.5952674
temp + silicate Indo-Pacific (26,26.5] -0.6075082 -0.6127084 -0.6040884
temp + silicate Indo-Pacific (26.5,26.75] -0.5992953 -0.6022474 -0.5934859
temp + silicate Indo-Pacific (26.75,27] -0.6495632 -0.6600895 -0.6388424
temp + silicate Indo-Pacific (27,27.25] -0.7595886 -0.7687879 -0.7515348
temp + silicate Indo-Pacific (27.25,27.5] -0.4739383 -0.5074598 -0.4301510
temp + silicate Indo-Pacific (27.5,27.75] 0.2406359 0.1968792 0.2943669
temp + silicate Indo-Pacific (27.75,27.85] 0.7102366 0.7055165 0.7156369
temp + silicate Indo-Pacific (27.85,27.95] 0.7217250 0.7200823 0.7237674
temp + silicate Indo-Pacific (27.95,28.05] 0.7645153 0.7539200 0.7756134
temp + silicate Indo-Pacific (28.05,28.1] 0.8228715 0.8086531 0.8371790
temp + silicate Indo-Pacific (28.1, Inf] 0.6401240 0.6178531 0.6589035 
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
a2f0d56 Donghe-Zhu 2021-01-23
28509fc Donghe-Zhu 2021-01-23
4c28e4a Donghe-Zhu 2021-01-22
24cc264 jens-daniel-mueller 2021-01-22
7891955 Donghe-Zhu 2021-01-21
d4cf1cb Donghe-Zhu 2021-01-21
1f3e5b6 jens-daniel-mueller 2021-01-20
0e7bdf1 jens-daniel-mueller 2021-01-15
4571843 jens-daniel-mueller 2021-01-14
b3564aa jens-daniel-mueller 2021-01-14
8d032c3 jens-daniel-mueller 2021-01-14
17dee1d jens-daniel-mueller 2021-01-13
7cdea0c jens-daniel-mueller 2021-01-06
fa85b93 jens-daniel-mueller 2021-01-06
e5cb81a Donghe-Zhu 2021-01-05
a499f10 Donghe-Zhu 2021-01-05
fb8a752 Donghe-Zhu 2020-12-23
8fae0b2 Donghe-Zhu 2020-12-21
c8b76b3 jens-daniel-mueller 2020-12-19 
# 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.3641068 0.3220785 0.4010690
aou Atlantic (26,26.5] 0.0006824 -0.1490483 0.1314943
aou Atlantic (26.5,26.75] -0.0062948 -0.1579981 0.1967323
aou Atlantic (26.75,27] -0.0976286 -0.3026430 0.1521067
aou Atlantic (27,27.25] -0.4094155 -0.5437903 -0.2581406
aou Atlantic (27.25,27.5] -0.3058281 -0.4324190 -0.1659721
aou Atlantic (27.5,27.75] -0.1405842 -0.1973962 -0.0792003
aou Atlantic (27.75,27.85] -0.0363833 -0.0798699 -0.0021351
aou Atlantic (27.85,27.95] -0.4960132 -0.5638601 -0.4402603
aou Atlantic (27.95,28.05] -0.1382259 -0.2071309 -0.0948120
aou Atlantic (28.05,28.1] -0.4321401 -0.4473223 -0.4196087
aou Atlantic (28.1,28.15] -0.8149108 -0.8152321 -0.8146366
aou Atlantic (28.15,28.2] -0.9082657 -0.9106662 -0.9044087
aou Atlantic (28.2, Inf] -0.8005904 -0.8103598 -0.7890016
aou Indo-Pacific (-Inf,26] -0.0963426 -0.2619510 0.0527317
aou Indo-Pacific (26,26.5] -0.4622333 -0.5976599 -0.2874627
aou Indo-Pacific (26.5,26.75] -0.7487964 -0.8166416 -0.6521349
aou Indo-Pacific (26.75,27] -0.8090749 -0.8672779 -0.7325671
aou Indo-Pacific (27,27.25] -0.7913142 -0.8638292 -0.6963558
aou Indo-Pacific (27.25,27.5] -0.6660582 -0.8326780 -0.4150807
aou Indo-Pacific (27.5,27.75] -0.2765598 -0.7500055 0.4099623
aou Indo-Pacific (27.75,27.85] -0.0409437 -0.6024882 0.6915518
aou Indo-Pacific (27.85,27.95] 0.1031628 -0.3818717 0.6528633
aou Indo-Pacific (27.95,28.05] -0.0466444 -0.4305249 0.4255847
aou Indo-Pacific (28.05,28.1] 0.0739391 -0.3094811 0.5178184
aou Indo-Pacific (28.1, Inf] 0.4921644 0.3588916 0.6198919
phosphate Atlantic (-Inf,26] 0.4402535 0.3629020 0.5016541
phosphate Atlantic (26,26.5] 0.0231824 -0.0838596 0.1036402
phosphate Atlantic (26.5,26.75] -0.1495355 -0.2835118 0.0364510
phosphate Atlantic (26.75,27] -0.2481112 -0.4310239 -0.0338118
phosphate Atlantic (27,27.25] -0.6458169 -0.7374433 -0.5386546
phosphate Atlantic (27.25,27.5] -0.7802243 -0.8292070 -0.7188560
phosphate Atlantic (27.5,27.75] -0.8950568 -0.9134885 -0.8730060
phosphate Atlantic (27.75,27.85] -0.9449144 -0.9522141 -0.9344548
phosphate Atlantic (27.85,27.95] -0.9722238 -0.9786287 -0.9633167
phosphate Atlantic (27.95,28.05] -0.9597184 -0.9620782 -0.9558823
phosphate Atlantic (28.05,28.1] -0.9885780 -0.9900367 -0.9867592
phosphate Atlantic (28.1,28.15] -0.9929095 -0.9940464 -0.9921523
phosphate Atlantic (28.15,28.2] -0.9937885 -0.9949452 -0.9924933
phosphate Atlantic (28.2, Inf] -0.9780168 -0.9807725 -0.9762173
phosphate Indo-Pacific (-Inf,26] -0.1432395 -0.3107995 0.0110294
phosphate Indo-Pacific (26,26.5] -0.5866252 -0.7102422 -0.4253541
phosphate Indo-Pacific (26.5,26.75] -0.8143623 -0.8717897 -0.7294770
phosphate Indo-Pacific (26.75,27] -0.8367100 -0.8811836 -0.7761624
phosphate Indo-Pacific (27,27.25] -0.8282547 -0.8861495 -0.7519136
phosphate Indo-Pacific (27.25,27.5] -0.7199561 -0.8639312 -0.5006583
phosphate Indo-Pacific (27.5,27.75] -0.3077663 -0.7568981 0.3480930
phosphate Indo-Pacific (27.75,27.85] 0.0012271 -0.5502468 0.7114011
phosphate Indo-Pacific (27.85,27.95] 0.1437893 -0.3287483 0.6726404
phosphate Indo-Pacific (27.95,28.05] 0.0387412 -0.3299029 0.4870988
phosphate Indo-Pacific (28.05,28.1] 0.1442290 -0.2179524 0.5594772
phosphate Indo-Pacific (28.1, Inf] 0.4638316 0.3455955 0.5767606
phosphate_star Atlantic (-Inf,26] 0.1011514 0.0498067 0.1452297
phosphate_star Atlantic (26,26.5] 0.1075981 0.0932518 0.1347536
phosphate_star Atlantic (26.5,26.75] -0.3313220 -0.3536565 -0.3039910
phosphate_star Atlantic (26.75,27] -0.4306697 -0.5390440 -0.3622338
phosphate_star Atlantic (27,27.25] -0.5758978 -0.6523484 -0.4921346
phosphate_star Atlantic (27.25,27.5] -0.7297849 -0.8069404 -0.6481326
phosphate_star Atlantic (27.5,27.75] -0.8019948 -0.8415493 -0.7617142
phosphate_star Atlantic (27.75,27.85] -0.8866319 -0.9281579 -0.8449938
phosphate_star Atlantic (27.85,27.95] -0.8956414 -0.9282086 -0.8570907
phosphate_star Atlantic (27.95,28.05] -0.8792209 -0.9100347 -0.8489023
phosphate_star Atlantic (28.05,28.1] -0.9357280 -0.9587066 -0.9091095
phosphate_star Atlantic (28.1,28.15] -0.9638590 -0.9773642 -0.9492232
phosphate_star Atlantic (28.15,28.2] -0.9770357 -0.9864388 -0.9669520
phosphate_star Atlantic (28.2, Inf] -0.9694754 -0.9788690 -0.9565256
phosphate_star Indo-Pacific (-Inf,26] -0.0067078 -0.1347902 0.1158792
phosphate_star Indo-Pacific (26,26.5] -0.5570677 -0.6022742 -0.5003744
phosphate_star Indo-Pacific (26.5,26.75] -0.5813075 -0.6059676 -0.5568690
phosphate_star Indo-Pacific (26.75,27] -0.5179954 -0.5369605 -0.4852368
phosphate_star Indo-Pacific (27,27.25] -0.6711656 -0.6945508 -0.6480354
phosphate_star Indo-Pacific (27.25,27.5] -0.5297383 -0.5561789 -0.5107753
phosphate_star Indo-Pacific (27.5,27.75] -0.1399733 -0.4999755 0.1090058
phosphate_star Indo-Pacific (27.75,27.85] 0.3815853 -0.2267951 0.7878502
phosphate_star Indo-Pacific (27.85,27.95] 0.3803007 -0.0181880 0.6902955
phosphate_star Indo-Pacific (27.95,28.05] 0.5629903 0.2472120 0.7891929
phosphate_star Indo-Pacific (28.05,28.1] 0.5188025 0.2076844 0.7609118
phosphate_star Indo-Pacific (28.1, Inf] -0.0214844 -0.1125968 0.0691077
sal Atlantic (-Inf,26] 0.6329389 0.6159084 0.6581242
sal Atlantic (26,26.5] -0.0950453 -0.1155099 -0.0586497
sal Atlantic (26.5,26.75] -0.0610730 -0.1616480 0.0515048
sal Atlantic (26.75,27] -0.1176359 -0.2419380 0.0752622
sal Atlantic (27,27.25] -0.0604729 -0.1930877 0.0621034
sal Atlantic (27.25,27.5] 0.2428429 0.0815750 0.4002206
sal Atlantic (27.5,27.75] -0.0662468 -0.1726693 0.0397511
sal Atlantic (27.75,27.85] 0.5391256 0.4596143 0.6329847
sal Atlantic (27.85,27.95] 0.5714413 0.5123516 0.6365731
sal Atlantic (27.95,28.05] 0.7692805 0.7226039 0.8143549
sal Atlantic (28.05,28.1] 0.8948865 0.8685197 0.9158452
sal Atlantic (28.1,28.15] 0.9402902 0.9260166 0.9592370
sal Atlantic (28.15,28.2] 0.9465060 0.9260655 0.9639575
sal Atlantic (28.2, Inf] 0.6762985 0.6390711 0.7033323
sal Indo-Pacific (-Inf,26] -0.0632689 -0.1296836 0.0033318
sal Indo-Pacific (26,26.5] 0.3003825 0.2763641 0.3304955
sal Indo-Pacific (26.5,26.75] 0.2654368 0.2306789 0.2835971
sal Indo-Pacific (26.75,27] 0.1828136 0.1196013 0.2442309
sal Indo-Pacific (27,27.25] 0.4436219 0.3919447 0.4911330
sal Indo-Pacific (27.25,27.5] 0.2142758 0.1368161 0.2771281
sal Indo-Pacific (27.5,27.75] -0.1593541 -0.4342577 0.2757749
sal Indo-Pacific (27.75,27.85] -0.3330548 -0.8060006 0.3286839
sal Indo-Pacific (27.85,27.95] -0.2070486 -0.6435653 0.3079577
sal Indo-Pacific (27.95,28.05] -0.3531654 -0.6858301 0.0827998
sal Indo-Pacific (28.05,28.1] -0.2501868 -0.5984770 0.1677732
sal Indo-Pacific (28.1, Inf] -0.0200636 -0.2202403 0.1856615
silicate Atlantic (-Inf,26] -0.6235293 -0.6822996 -0.5931277
silicate Atlantic (26,26.5] -0.5093404 -0.5740153 -0.4379810
silicate Atlantic (26.5,26.75] -0.3666501 -0.4750659 -0.2218988
silicate Atlantic (26.75,27] -0.4560659 -0.6294062 -0.2548992
silicate Atlantic (27,27.25] -0.7908733 -0.8530734 -0.7118891
silicate Atlantic (27.25,27.5] -0.8753026 -0.9017861 -0.8384385
silicate Atlantic (27.5,27.75] -0.9171498 -0.9185239 -0.9155889
silicate Atlantic (27.75,27.85] -0.9290193 -0.9410706 -0.9174835
silicate Atlantic (27.85,27.95] -0.9067008 -0.9285901 -0.8796809
silicate Atlantic (27.95,28.05] -0.9174158 -0.9251335 -0.9103238
silicate Atlantic (28.05,28.1] -0.9658083 -0.9681704 -0.9617209
silicate Atlantic (28.1,28.15] -0.9775239 -0.9819129 -0.9730022
silicate Atlantic (28.15,28.2] -0.9265620 -0.9353244 -0.9138045
silicate Atlantic (28.2, Inf] -0.7478464 -0.7758269 -0.7087256
silicate Indo-Pacific (-Inf,26] -0.1261189 -0.2814228 0.0087882
silicate Indo-Pacific (26,26.5] -0.6594363 -0.7532053 -0.5390698
silicate Indo-Pacific (26.5,26.75] -0.7437066 -0.7947383 -0.6664905
silicate Indo-Pacific (26.75,27] -0.6775849 -0.7293239 -0.6097373
silicate Indo-Pacific (27,27.25] -0.6931388 -0.7657052 -0.5987968
silicate Indo-Pacific (27.25,27.5] -0.6167679 -0.7744134 -0.3826825
silicate Indo-Pacific (27.5,27.75] -0.2503160 -0.7028054 0.4027220
silicate Indo-Pacific (27.75,27.85] 0.0275295 -0.5275915 0.7361816
silicate Indo-Pacific (27.85,27.95] 0.1520005 -0.3257949 0.6827522
silicate Indo-Pacific (27.95,28.05] 0.0320706 -0.3444162 0.4920461
silicate Indo-Pacific (28.05,28.1] 0.1409720 -0.2247137 0.5626965
silicate Indo-Pacific (28.1, Inf] 0.3763762 0.2205045 0.5301459
temp Atlantic (-Inf,26] -0.0491292 -0.0822790 -0.0015824
temp Atlantic (26,26.5] -0.2374342 -0.2596759 -0.1960558
temp Atlantic (26.5,26.75] -0.0741267 -0.1555728 0.0127791
temp Atlantic (26.75,27] -0.0490542 -0.1496181 0.1194733
temp Atlantic (27,27.25] -0.0347625 -0.1491391 0.0724962
temp Atlantic (27.25,27.5] 0.2081095 0.0472510 0.3629162
temp Atlantic (27.5,27.75] -0.0282821 -0.1194432 0.0652897
temp Atlantic (27.75,27.85] 0.5506288 0.4761536 0.6358221
temp Atlantic (27.85,27.95] 0.5969218 0.5410243 0.6536703
temp Atlantic (27.95,28.05] 0.7524867 0.7064418 0.7982045
temp Atlantic (28.05,28.1] 0.8605036 0.8267967 0.8856737
temp Atlantic (28.1,28.15] 0.9034743 0.8846602 0.9237863
temp Atlantic (28.15,28.2] 0.9228584 0.9004739 0.9425501
temp Atlantic (28.2, Inf] 0.8171427 0.7975944 0.8462547
temp Indo-Pacific (-Inf,26] -0.3731855 -0.4629043 -0.2790149
temp Indo-Pacific (26,26.5] 0.2568546 0.2169667 0.2973742
temp Indo-Pacific (26.5,26.75] 0.2219722 0.1804235 0.2452792
temp Indo-Pacific (26.75,27] 0.1166388 0.0503751 0.1789004
temp Indo-Pacific (27,27.25] 0.4130346 0.3667944 0.4548654
temp Indo-Pacific (27.25,27.5] 0.3238948 0.2998795 0.3419444
temp Indo-Pacific (27.5,27.75] 0.0779621 -0.1887942 0.4595552
temp Indo-Pacific (27.75,27.85] -0.3414179 -0.8015871 0.3126022
temp Indo-Pacific (27.85,27.95] -0.2656797 -0.6754168 0.2332146
temp Indo-Pacific (27.95,28.05] -0.4304128 -0.7519012 -0.0040522
temp Indo-Pacific (28.05,28.1] -0.2752135 -0.6240377 0.1526569
temp Indo-Pacific (28.1, Inf] 0.2410560 0.0787373 0.4021228 
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.2

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

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

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

other attached packages:
 [1] corrr_0.4.3      broom_0.7.3      kableExtra_1.3.1 knitr_1.30      
 [5] olsrr_0.5.3      GGally_2.0.0     lubridate_1.7.9  metR_0.9.0      
 [9] scico_1.2.0      patchwork_1.1.1  collapse_1.5.0   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.3    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.10            
[34] readxl_1.3.1             rstudioapi_0.13          farver_2.0.3            
[37] generics_0.1.0           jsonlite_1.7.2           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.2 data.table_1.13.6       
[67] modelr_0.1.8             vctrs_0.3.6              httpuv_1.5.4            
[70] cellranger_1.1.0         gtable_0.3.0             reshape_0.8.8           
[73] assertthat_0.2.1         xfun_0.20                openxlsx_4.2.3          
[76] RcppEigen_0.3.3.9.1      later_1.1.0.1            viridisLite_0.3.0       
[79] ellipsis_0.3.1           here_1.0.1