Last updated: 2021-01-25

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.


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.

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.

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.

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

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

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

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 560dcf7. 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 05ffb0c Donghe-Zhu 2021-01-25 Build site.
html 8b97165 Donghe-Zhu 2021-01-25 Build site.
html c569946 Donghe-Zhu 2021-01-24 Build site.
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
05ffb0c Donghe-Zhu 2021-01-25
8b97165 Donghe-Zhu 2021-01-25
c569946 Donghe-Zhu 2021-01-24
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
05ffb0c Donghe-Zhu 2021-01-25
8b97165 Donghe-Zhu 2021-01-25
c569946 Donghe-Zhu 2021-01-24
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
05ffb0c Donghe-Zhu 2021-01-25
8b97165 Donghe-Zhu 2021-01-25
c569946 Donghe-Zhu 2021-01-24
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 + nitrate Atlantic (-Inf,26] 0.8566630 0.8400495 0.8876974
aou + nitrate Atlantic (26,26.5] 0.9170520 0.9110515 0.9216364
aou + nitrate Atlantic (26.5,26.75] 0.9733858 0.9718730 0.9749248
aou + nitrate Atlantic (26.75,27] 0.9829109 0.9806161 0.9844866
aou + nitrate Atlantic (27,27.25] 0.9612570 0.9568927 0.9645505
aou + nitrate Atlantic (27.25,27.5] 0.8309181 0.8188576 0.8442941
aou + nitrate Atlantic (27.5,27.75] 0.5552282 0.5479099 0.5652643
aou + nitrate Atlantic (27.75,27.85] 0.3547108 0.3377525 0.3831833
aou + nitrate Atlantic (27.85,27.95] 0.5561316 0.4622259 0.6514764
aou + nitrate Atlantic (27.95,28.05] 0.1970693 0.1613187 0.2213080
aou + nitrate Atlantic (28.05,28.1] 0.4788652 0.4240406 0.5424664
aou + nitrate Atlantic (28.1,28.15] 0.8373626 0.8140925 0.8578440
aou + nitrate Atlantic (28.15,28.2] 0.9229742 0.9144140 0.9307013
aou + nitrate Atlantic (28.2, Inf] 0.8770807 0.8702588 0.8838716
aou + nitrate Indo-Pacific (-Inf,26] 0.9322982 0.9304671 0.9343005
aou + nitrate Indo-Pacific (26,26.5] 0.8518576 0.8486442 0.8559370
aou + nitrate Indo-Pacific (26.5,26.75] 0.7885271 0.7837177 0.7977763
aou + nitrate Indo-Pacific (26.75,27] 0.8509457 0.8472746 0.8571177
aou + nitrate Indo-Pacific (27,27.25] 0.9332009 0.9291595 0.9372775
aou + nitrate Indo-Pacific (27.25,27.5] 0.9557019 0.9546184 0.9566444
aou + nitrate Indo-Pacific (27.5,27.75] 0.9756488 0.9755546 0.9756974
aou + nitrate Indo-Pacific (27.75,27.85] 0.9920246 0.9917943 0.9924678
aou + nitrate Indo-Pacific (27.85,27.95] 0.9909981 0.9904821 0.9918483
aou + nitrate Indo-Pacific (27.95,28.05] 0.9863744 0.9861973 0.9865229
aou + nitrate Indo-Pacific (28.05,28.1] 0.9880668 0.9877930 0.9883951
aou + nitrate Indo-Pacific (28.1, Inf] 0.9891921 0.9886841 0.9894492
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
nitrate + phosphate Atlantic (-Inf,26] 0.6913320 0.6358544 0.7735436
nitrate + phosphate Atlantic (26,26.5] 0.8329157 0.8216064 0.8508597
nitrate + phosphate Atlantic (26.5,26.75] 0.9668921 0.9651352 0.9678758
nitrate + phosphate Atlantic (26.75,27] 0.9800420 0.9789957 0.9810574
nitrate + phosphate Atlantic (27,27.25] 0.9863334 0.9859656 0.9867931
nitrate + phosphate Atlantic (27.25,27.5] 0.9912444 0.9910142 0.9914977
nitrate + phosphate Atlantic (27.5,27.75] 0.9945081 0.9943694 0.9946066
nitrate + phosphate Atlantic (27.75,27.85] 0.9968880 0.9968477 0.9969657
nitrate + phosphate Atlantic (27.85,27.95] 0.9987997 0.9986913 0.9989229
nitrate + phosphate Atlantic (27.95,28.05] 0.9994197 0.9993947 0.9994373
nitrate + phosphate Atlantic (28.05,28.1] 0.9997863 0.9997740 0.9997996
nitrate + phosphate Atlantic (28.1,28.15] 0.9998491 0.9998415 0.9998560
nitrate + phosphate Atlantic (28.15,28.2] 0.9983882 0.9983039 0.9984844
nitrate + phosphate Atlantic (28.2, Inf] 0.9866101 0.9856099 0.9877020
nitrate + phosphate Indo-Pacific (-Inf,26] 0.9545388 0.9524337 0.9577585
nitrate + phosphate Indo-Pacific (26,26.5] 0.8699687 0.8678049 0.8742210
nitrate + phosphate Indo-Pacific (26.5,26.75] 0.8110580 0.8050182 0.8217624
nitrate + phosphate Indo-Pacific (26.75,27] 0.8767211 0.8724790 0.8836490
nitrate + phosphate Indo-Pacific (27,27.25] 0.9447745 0.9413047 0.9480273
nitrate + phosphate Indo-Pacific (27.25,27.5] 0.9619597 0.9611442 0.9627674
nitrate + phosphate Indo-Pacific (27.5,27.75] 0.9827379 0.9823686 0.9830756
nitrate + phosphate Indo-Pacific (27.75,27.85] 0.9957061 0.9953913 0.9958890
nitrate + phosphate Indo-Pacific (27.85,27.95] 0.9974660 0.9972014 0.9978710
nitrate + phosphate Indo-Pacific (27.95,28.05] 0.9990242 0.9989610 0.9990911
nitrate + phosphate Indo-Pacific (28.05,28.1] 0.9986147 0.9984581 0.9986949
nitrate + phosphate Indo-Pacific (28.1, Inf] 0.9951368 0.9950077 0.9952227
nitrate + phosphate_star Atlantic (-Inf,26] 0.1391998 0.0453209 0.2437980
nitrate + phosphate_star Atlantic (26,26.5] 0.4147133 0.3635791 0.4528847
nitrate + phosphate_star Atlantic (26.5,26.75] 0.2842813 0.2484607 0.3111819
nitrate + phosphate_star Atlantic (26.75,27] 0.0018870 -0.0148950 0.0150087
nitrate + phosphate_star Atlantic (27,27.25] 0.2106063 0.1764811 0.2474356
nitrate + phosphate_star Atlantic (27.25,27.5] 0.4045003 0.3874761 0.4139486
nitrate + phosphate_star Atlantic (27.5,27.75] 0.5925333 0.5859924 0.5984410
nitrate + phosphate_star Atlantic (27.75,27.85] 0.7466715 0.7265802 0.7636152
nitrate + phosphate_star Atlantic (27.85,27.95] 0.9088651 0.8973759 0.9216017
nitrate + phosphate_star Atlantic (27.95,28.05] 0.9028467 0.9017683 0.9049274
nitrate + phosphate_star Atlantic (28.05,28.1] 0.9350493 0.9314337 0.9379317
nitrate + phosphate_star Atlantic (28.1,28.15] 0.9665714 0.9640628 0.9680389
nitrate + phosphate_star Atlantic (28.15,28.2] 0.9785881 0.9780231 0.9791768
nitrate + phosphate_star Atlantic (28.2, Inf] 0.9587032 0.9572841 0.9595635
nitrate + phosphate_star Indo-Pacific (-Inf,26] 0.7302540 0.7281127 0.7335866
nitrate + phosphate_star Indo-Pacific (26,26.5] 0.5635486 0.5605704 0.5680878
nitrate + phosphate_star Indo-Pacific (26.5,26.75] 0.5741730 0.5667133 0.5828366
nitrate + phosphate_star Indo-Pacific (26.75,27] 0.7061532 0.6930122 0.7180860
nitrate + phosphate_star Indo-Pacific (27,27.25] 0.7905769 0.7779339 0.8048866
nitrate + phosphate_star Indo-Pacific (27.25,27.5] 0.4486140 0.3938967 0.4896526
nitrate + phosphate_star Indo-Pacific (27.5,27.75] -0.2035846 -0.2679761 -0.1476173
nitrate + phosphate_star Indo-Pacific (27.75,27.85] -0.5795626 -0.5841819 -0.5709257
nitrate + phosphate_star Indo-Pacific (27.85,27.95] -0.3993262 -0.4229179 -0.3816230
nitrate + phosphate_star Indo-Pacific (27.95,28.05] -0.3507523 -0.3820853 -0.3241588
nitrate + phosphate_star Indo-Pacific (28.05,28.1] -0.2649867 -0.3323095 -0.2116992
nitrate + phosphate_star Indo-Pacific (28.1, Inf] 0.1672575 0.1388275 0.1852209
nitrate + silicate Atlantic (-Inf,26] 0.5921507 0.5624675 0.6098511
nitrate + silicate Atlantic (26,26.5] 0.5610505 0.5497846 0.5669760
nitrate + silicate Atlantic (26.5,26.75] 0.6989596 0.6753673 0.7239189
nitrate + silicate Atlantic (26.75,27] 0.8608930 0.8535129 0.8659875
nitrate + silicate Atlantic (27,27.25] 0.8800101 0.8780550 0.8820966
nitrate + silicate Atlantic (27.25,27.5] 0.8713724 0.8683804 0.8730735
nitrate + silicate Atlantic (27.5,27.75] 0.8817783 0.8784194 0.8849732
nitrate + silicate Atlantic (27.75,27.85] 0.8613827 0.8604182 0.8632989
nitrate + silicate Atlantic (27.85,27.95] 0.9232382 0.9166427 0.9314203
nitrate + silicate Atlantic (27.95,28.05] 0.9714387 0.9705597 0.9722333
nitrate + silicate Atlantic (28.05,28.1] 0.9835191 0.9825409 0.9843418
nitrate + silicate Atlantic (28.1,28.15] 0.9860056 0.9855119 0.9864498
nitrate + silicate Atlantic (28.15,28.2] 0.9262125 0.9246027 0.9286829
nitrate + silicate Atlantic (28.2, Inf] 0.6475740 0.6443481 0.6502341
nitrate + silicate Indo-Pacific (-Inf,26] 0.8591922 0.8567866 0.8616604
nitrate + silicate Indo-Pacific (26,26.5] 0.7921975 0.7869786 0.8008042
nitrate + silicate Indo-Pacific (26.5,26.75] 0.8328437 0.8284302 0.8383708
nitrate + silicate Indo-Pacific (26.75,27] 0.8862552 0.8827027 0.8901026
nitrate + silicate Indo-Pacific (27,27.25] 0.9464613 0.9436981 0.9487636
nitrate + silicate Indo-Pacific (27.25,27.5] 0.9655265 0.9645432 0.9662119
nitrate + silicate Indo-Pacific (27.5,27.75] 0.9815195 0.9808706 0.9819209
nitrate + silicate Indo-Pacific (27.75,27.85] 0.9910303 0.9904377 0.9915907
nitrate + silicate Indo-Pacific (27.85,27.95] 0.9911812 0.9909684 0.9913855
nitrate + silicate Indo-Pacific (27.95,28.05] 0.9923713 0.9922186 0.9926602
nitrate + silicate Indo-Pacific (28.05,28.1] 0.9915592 0.9911950 0.9918100
nitrate + silicate Indo-Pacific (28.1, Inf] 0.8433853 0.8430264 0.8436477
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 + nitrate Atlantic (-Inf,26] -0.6745164 -0.6845712 -0.6657979
sal + nitrate Atlantic (26,26.5] -0.4262512 -0.4674110 -0.3860388
sal + nitrate Atlantic (26.5,26.75] 0.1902400 0.1732906 0.2220760
sal + nitrate Atlantic (26.75,27] 0.3985253 0.3924883 0.4104851
sal + nitrate Atlantic (27,27.25] 0.2442075 0.1938315 0.2980741
sal + nitrate Atlantic (27.25,27.5] -0.0063061 -0.0563914 0.0541962
sal + nitrate Atlantic (27.5,27.75] 0.2213963 0.1839875 0.2693849
sal + nitrate Atlantic (27.75,27.85] -0.3534892 -0.3961882 -0.3204810
sal + nitrate Atlantic (27.85,27.95] -0.6079697 -0.6345310 -0.5946473
sal + nitrate Atlantic (27.95,28.05] -0.8453962 -0.8520500 -0.8364979
sal + nitrate Atlantic (28.05,28.1] -0.9150235 -0.9174325 -0.9137779
sal + nitrate Atlantic (28.1,28.15] -0.9553978 -0.9576818 -0.9540523
sal + nitrate Atlantic (28.15,28.2] -0.9580031 -0.9600778 -0.9553757
sal + nitrate Atlantic (28.2, Inf] -0.7910473 -0.7984536 -0.7832206
sal + nitrate Indo-Pacific (-Inf,26] -0.3670441 -0.3692076 -0.3650932
sal + nitrate Indo-Pacific (26,26.5] -0.4023262 -0.4056480 -0.3978450
sal + nitrate Indo-Pacific (26.5,26.75] -0.4361355 -0.4388556 -0.4342895
sal + nitrate Indo-Pacific (26.75,27] -0.5862602 -0.6040734 -0.5735445
sal + nitrate Indo-Pacific (27,27.25] -0.7526178 -0.7590138 -0.7467659
sal + nitrate Indo-Pacific (27.25,27.5] -0.4286436 -0.4742571 -0.3778009
sal + nitrate Indo-Pacific (27.5,27.75] 0.4176335 0.3764700 0.4764062
sal + nitrate Indo-Pacific (27.75,27.85] 0.7579713 0.7463217 0.7674463
sal + nitrate Indo-Pacific (27.85,27.95] 0.7703251 0.7698279 0.7708326
sal + nitrate Indo-Pacific (27.95,28.05] 0.8048904 0.7929514 0.8166602
sal + nitrate Indo-Pacific (28.05,28.1] 0.8205091 0.8031262 0.8344112
sal + nitrate Indo-Pacific (28.1, Inf] 0.3807919 0.3660197 0.3933452
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
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
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
05ffb0c Donghe-Zhu 2021-01-25
8b97165 Donghe-Zhu 2021-01-25
c569946 Donghe-Zhu 2021-01-24
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
nitrate Atlantic (-Inf,26] 0.0269936 -0.0131800 0.0546032
nitrate Atlantic (26,26.5] -0.2259605 -0.3718580 -0.0932580
nitrate Atlantic (26.5,26.75] -0.1653403 -0.3181858 0.0387210
nitrate Atlantic (26.75,27] -0.2213816 -0.4183202 0.0159495
nitrate Atlantic (27,27.25] -0.6106341 -0.7145395 -0.4928023
nitrate Atlantic (27.25,27.5] -0.7414915 -0.8045254 -0.6657001
nitrate Atlantic (27.5,27.75] -0.8747258 -0.8988562 -0.8462087
nitrate Atlantic (27.75,27.85] -0.9258115 -0.9365958 -0.9113354
nitrate Atlantic (27.85,27.95] -0.9715810 -0.9767093 -0.9645860
nitrate Atlantic (27.95,28.05] -0.9579371 -0.9601390 -0.9553217
nitrate Atlantic (28.05,28.1] -0.9883567 -0.9897697 -0.9869323
nitrate Atlantic (28.1,28.15] -0.9927733 -0.9940200 -0.9919608
nitrate Atlantic (28.15,28.2] -0.9944769 -0.9961191 -0.9918177
nitrate Atlantic (28.2, Inf] -0.9657490 -0.9718004 -0.9585585
nitrate Indo-Pacific (-Inf,26] -0.1152997 -0.2783666 0.0347572
nitrate Indo-Pacific (26,26.5] -0.4509442 -0.5737159 -0.3051158
nitrate Indo-Pacific (26.5,26.75] -0.6634411 -0.7303425 -0.5723095
nitrate Indo-Pacific (26.75,27] -0.6949987 -0.7358133 -0.6274493
nitrate Indo-Pacific (27,27.25] -0.7471559 -0.8131396 -0.6621821
nitrate Indo-Pacific (27.25,27.5] -0.6685436 -0.8223604 -0.4354675
nitrate Indo-Pacific (27.5,27.75] -0.3070664 -0.7610755 0.3561272
nitrate Indo-Pacific (27.75,27.85] -0.0078834 -0.5645044 0.7098004
nitrate Indo-Pacific (27.85,27.95] 0.1536357 -0.3199871 0.6799997
nitrate Indo-Pacific (27.95,28.05] 0.0415715 -0.3287273 0.4920884
nitrate Indo-Pacific (28.05,28.1] 0.1620539 -0.2011728 0.5790774
nitrate Indo-Pacific (28.1, Inf] 0.5130713 0.4033632 0.6170175
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
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