eMLR - assumption testing
Jens Daniel Müller and Donghe Zhu
23 January, 2021
Last updated: 2021-01-23
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 0819197. 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 | 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 |
|---|---|---|
| 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 |
|---|---|---|
| 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 |
|---|---|---|
| 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.8] | 0.9675284 | 0.9669017 | 0.9684557 |
| aou + nitrate | Atlantic | (26.8,27.2] | 0.9607385 | 0.9572266 | 0.9641624 |
| aou + nitrate | Atlantic | (27.2,27.8] | 0.6754123 | 0.6696006 | 0.6790898 |
| aou + nitrate | Atlantic | (27.8,27.9] | 0.4175220 | 0.3733353 | 0.4749548 |
| aou + nitrate | Atlantic | (27.9,28.1] | 0.2985312 | 0.2748615 | 0.3332581 |
| aou + nitrate | Atlantic | (28.1,28.2] | 0.8767484 | 0.8627073 | 0.8911171 |
| 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.8] | 0.8295417 | 0.8245161 | 0.8356921 |
| aou + nitrate | Indo-Pacific | (26.8,27.2] | 0.8882055 | 0.8836577 | 0.8908884 |
| aou + nitrate | Indo-Pacific | (27.2,27.8] | 0.9368402 | 0.9365835 | 0.9369906 |
| aou + nitrate | Indo-Pacific | (27.8,27.9] | 0.9900530 | 0.9896943 | 0.9902486 |
| aou + nitrate | Indo-Pacific | (27.9,28.1] | 0.9868329 | 0.9866892 | 0.9870343 |
| 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.8] | 0.9496684 | 0.9465399 | 0.9513803 |
| aou + phosphate | Atlantic | (26.8,27.2] | 0.9312913 | 0.9280670 | 0.9343924 |
| aou + phosphate | Atlantic | (27.2,27.8] | 0.6039787 | 0.5976782 | 0.6073144 |
| aou + phosphate | Atlantic | (27.8,27.9] | 0.3652683 | 0.3192898 | 0.4258875 |
| aou + phosphate | Atlantic | (27.9,28.1] | 0.2813431 | 0.2573030 | 0.3159401 |
| aou + phosphate | Atlantic | (28.1,28.2] | 0.8804204 | 0.8663831 | 0.8946419 |
| 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.8] | 0.9632100 | 0.9630016 | 0.9634902 |
| aou + phosphate | Indo-Pacific | (26.8,27.2] | 0.9733199 | 0.9722670 | 0.9744779 |
| aou + phosphate | Indo-Pacific | (27.2,27.8] | 0.9716809 | 0.9715684 | 0.9717856 |
| aou + phosphate | Indo-Pacific | (27.8,27.9] | 0.9941190 | 0.9940499 | 0.9942292 |
| aou + phosphate | Indo-Pacific | (27.9,28.1] | 0.9880234 | 0.9876597 | 0.9882100 |
| 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.8] | 0.3953910 | 0.3844395 | 0.4049360 |
| aou + phosphate_star | Atlantic | (26.8,27.2] | -0.0316667 | -0.0439722 | -0.0236420 |
| aou + phosphate_star | Atlantic | (27.2,27.8] | -0.2756144 | -0.2876416 | -0.2675349 |
| aou + phosphate_star | Atlantic | (27.8,27.9] | -0.1633684 | -0.2327585 | -0.0753898 |
| aou + phosphate_star | Atlantic | (27.9,28.1] | -0.1120431 | -0.1329343 | -0.0772664 |
| aou + phosphate_star | Atlantic | (28.1,28.2] | 0.7429988 | 0.7191048 | 0.7661254 |
| 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.8] | 0.4823886 | 0.4715406 | 0.4913736 |
| aou + phosphate_star | Indo-Pacific | (26.8,27.2] | 0.5797755 | 0.5592859 | 0.5994568 |
| aou + phosphate_star | Indo-Pacific | (27.2,27.8] | -0.0108611 | -0.0551713 | 0.0261256 |
| aou + phosphate_star | Indo-Pacific | (27.8,27.9] | -0.5373781 | -0.5446968 | -0.5289580 |
| aou + phosphate_star | Indo-Pacific | (27.9,28.1] | -0.4488572 | -0.4881786 | -0.4171239 |
| 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.8] | 0.5694262 | 0.5371788 | 0.6012561 |
| aou + silicate | Atlantic | (26.8,27.2] | 0.6915728 | 0.6867709 | 0.6982732 |
| aou + silicate | Atlantic | (27.2,27.8] | 0.2845534 | 0.2820066 | 0.2861740 |
| aou + silicate | Atlantic | (27.8,27.9] | 0.0800795 | 0.0367501 | 0.1418396 |
| aou + silicate | Atlantic | (27.9,28.1] | 0.2738368 | 0.2513386 | 0.3068672 |
| aou + silicate | Atlantic | (28.1,28.2] | 0.9011302 | 0.8896725 | 0.9099003 |
| 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.8] | 0.7900078 | 0.7883018 | 0.7912806 |
| aou + silicate | Indo-Pacific | (26.8,27.2] | 0.9020126 | 0.9009748 | 0.9027727 |
| aou + silicate | Indo-Pacific | (27.2,27.8] | 0.9364233 | 0.9359958 | 0.9368634 |
| aou + silicate | Indo-Pacific | (27.8,27.9] | 0.9787601 | 0.9785030 | 0.9792556 |
| aou + silicate | Indo-Pacific | (27.9,28.1] | 0.9817417 | 0.9810080 | 0.9822374 |
| 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.8] | 0.9591475 | 0.9548624 | 0.9621075 |
| nitrate + phosphate | Atlantic | (26.8,27.2] | 0.9836361 | 0.9828995 | 0.9841539 |
| nitrate + phosphate | Atlantic | (27.2,27.8] | 0.9927647 | 0.9925744 | 0.9929690 |
| nitrate + phosphate | Atlantic | (27.8,27.9] | 0.9978195 | 0.9977514 | 0.9979372 |
| nitrate + phosphate | Atlantic | (27.9,28.1] | 0.9994918 | 0.9994712 | 0.9995044 |
| nitrate + phosphate | Atlantic | (28.1,28.2] | 0.9984528 | 0.9983930 | 0.9985359 |
| 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.8] | 0.8566760 | 0.8525113 | 0.8633376 |
| nitrate + phosphate | Indo-Pacific | (26.8,27.2] | 0.9182275 | 0.9143447 | 0.9207776 |
| nitrate + phosphate | Indo-Pacific | (27.2,27.8] | 0.9704398 | 0.9699390 | 0.9708864 |
| nitrate + phosphate | Indo-Pacific | (27.8,27.9] | 0.9962705 | 0.9960133 | 0.9966049 |
| nitrate + phosphate | Indo-Pacific | (27.9,28.1] | 0.9988187 | 0.9987140 | 0.9988962 |
| 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.8] | 0.5073718 | 0.4925110 | 0.5150926 |
| nitrate + phosphate_star | Atlantic | (26.8,27.2] | 0.2100103 | 0.1876533 | 0.2240680 |
| nitrate + phosphate_star | Atlantic | (27.2,27.8] | 0.4954754 | 0.4791266 | 0.5056750 |
| nitrate + phosphate_star | Atlantic | (27.8,27.9] | 0.8240979 | 0.8117461 | 0.8382730 |
| nitrate + phosphate_star | Atlantic | (27.9,28.1] | 0.9129833 | 0.9121098 | 0.9144721 |
| nitrate + phosphate_star | Atlantic | (28.1,28.2] | 0.9713217 | 0.9698060 | 0.9722627 |
| 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.8] | 0.6242012 | 0.6200619 | 0.6269661 |
| nitrate + phosphate_star | Indo-Pacific | (26.8,27.2] | 0.7323070 | 0.7180748 | 0.7436630 |
| nitrate + phosphate_star | Indo-Pacific | (27.2,27.8] | 0.2571960 | 0.2146248 | 0.2913004 |
| nitrate + phosphate_star | Indo-Pacific | (27.8,27.9] | -0.4410618 | -0.4506402 | -0.4304908 |
| nitrate + phosphate_star | Indo-Pacific | (27.9,28.1] | -0.3046793 | -0.3454687 | -0.2707671 |
| 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.8] | 0.6881048 | 0.6663946 | 0.7131034 |
| nitrate + silicate | Atlantic | (26.8,27.2] | 0.8369499 | 0.8359247 | 0.8387284 |
| nitrate + silicate | Atlantic | (27.2,27.8] | 0.8516105 | 0.8465705 | 0.8553136 |
| nitrate + silicate | Atlantic | (27.8,27.9] | 0.8854525 | 0.8820684 | 0.8910475 |
| nitrate + silicate | Atlantic | (27.9,28.1] | 0.9718818 | 0.9716257 | 0.9723141 |
| nitrate + silicate | Atlantic | (28.1,28.2] | 0.9266658 | 0.9242985 | 0.9286199 |
| 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.8] | 0.8298011 | 0.8260938 | 0.8321325 |
| nitrate + silicate | Indo-Pacific | (26.8,27.2] | 0.9239195 | 0.9209835 | 0.9258492 |
| nitrate + silicate | Indo-Pacific | (27.2,27.8] | 0.9707036 | 0.9697212 | 0.9715125 |
| nitrate + silicate | Indo-Pacific | (27.8,27.9] | 0.9897977 | 0.9896575 | 0.9900180 |
| nitrate + silicate | Indo-Pacific | (27.9,28.1] | 0.9901452 | 0.9898784 | 0.9903209 |
| 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.8] | 0.6540787 | 0.6511879 | 0.6587139 |
| phosphate + phosphate_star | Atlantic | (26.8,27.2] | 0.3276651 | 0.3074247 | 0.3404228 |
| phosphate + phosphate_star | Atlantic | (27.2,27.8] | 0.5824100 | 0.5684769 | 0.5920900 |
| phosphate + phosphate_star | Atlantic | (27.8,27.9] | 0.8559218 | 0.8449976 | 0.8679638 |
| phosphate + phosphate_star | Atlantic | (27.9,28.1] | 0.9207654 | 0.9196572 | 0.9223453 |
| phosphate + phosphate_star | Atlantic | (28.1,28.2] | 0.9707592 | 0.9692740 | 0.9718693 |
| 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.8] | 0.6947827 | 0.6868098 | 0.7022858 |
| phosphate + phosphate_star | Indo-Pacific | (26.8,27.2] | 0.7431288 | 0.7294760 | 0.7559827 |
| phosphate + phosphate_star | Indo-Pacific | (27.2,27.8] | 0.2235891 | 0.1804949 | 0.2593321 |
| phosphate + phosphate_star | Indo-Pacific | (27.8,27.9] | -0.4470854 | -0.4555360 | -0.4375290 |
| phosphate + phosphate_star | Indo-Pacific | (27.9,28.1] | -0.3097406 | -0.3493866 | -0.2771317 |
| 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.8] | -0.1692044 | -0.1808224 | -0.1602381 |
| sal + aou | Atlantic | (26.8,27.2] | 0.3572027 | 0.3473294 | 0.3724033 |
| sal + aou | Atlantic | (27.2,27.8] | 0.5488586 | 0.5356102 | 0.5658291 |
| sal + aou | Atlantic | (27.8,27.9] | 0.5169949 | 0.4400365 | 0.5795969 |
| sal + aou | Atlantic | (27.9,28.1] | 0.0928752 | 0.0612795 | 0.1127520 |
| sal + aou | Atlantic | (28.1,28.2] | -0.7154416 | -0.7344539 | -0.6963308 |
| 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.8] | -0.1935401 | -0.1987023 | -0.1890721 |
| sal + aou | Indo-Pacific | (26.8,27.2] | -0.5076043 | -0.5223439 | -0.4922192 |
| sal + aou | Indo-Pacific | (27.2,27.8] | 0.0632330 | 0.0108988 | 0.1251897 |
| sal + aou | Indo-Pacific | (27.8,27.9] | 0.8099125 | 0.8080991 | 0.8113422 |
| sal + aou | Indo-Pacific | (27.9,28.1] | 0.8646913 | 0.8539576 | 0.8752300 |
| 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.8] | -0.2764275 | -0.2846451 | -0.2715476 |
| sal + nitrate | Atlantic | (26.8,27.2] | 0.1952165 | 0.1779190 | 0.2247109 |
| sal + nitrate | Atlantic | (27.2,27.8] | 0.0978980 | 0.0592024 | 0.1461238 |
| sal + nitrate | Atlantic | (27.8,27.9] | -0.4749783 | -0.5067080 | -0.4557831 |
| sal + nitrate | Atlantic | (27.9,28.1] | -0.8676241 | -0.8699896 | -0.8632727 |
| sal + nitrate | Atlantic | (28.1,28.2] | -0.9470035 | -0.9489062 | -0.9459498 |
| 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.8] | -0.4317471 | -0.4362667 | -0.4268347 |
| sal + nitrate | Indo-Pacific | (26.8,27.2] | -0.6558447 | -0.6665197 | -0.6473107 |
| sal + nitrate | Indo-Pacific | (27.2,27.8] | -0.1442509 | -0.1952522 | -0.0798786 |
| sal + nitrate | Indo-Pacific | (27.8,27.9] | 0.7714455 | 0.7652774 | 0.7762223 |
| sal + nitrate | Indo-Pacific | (27.9,28.1] | 0.8088811 | 0.7945059 | 0.8224889 |
| 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.8] | -0.3963564 | -0.4055001 | -0.3897335 |
| sal + phosphate | Atlantic | (26.8,27.2] | 0.0817104 | 0.0612922 | 0.1143760 |
| sal + phosphate | Atlantic | (27.2,27.8] | -0.0068584 | -0.0446348 | 0.0410222 |
| sal + phosphate | Atlantic | (27.8,27.9] | -0.5240815 | -0.5541684 | -0.5066904 |
| sal + phosphate | Atlantic | (27.9,28.1] | -0.8762359 | -0.8789465 | -0.8714096 |
| sal + phosphate | Atlantic | (28.1,28.2] | -0.9352519 | -0.9372711 | -0.9339332 |
| 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.8] | -0.3875433 | -0.3913499 | -0.3833658 |
| sal + phosphate | Indo-Pacific | (26.8,27.2] | -0.6072882 | -0.6210547 | -0.5951134 |
| sal + phosphate | Indo-Pacific | (27.2,27.8] | -0.1089875 | -0.1628505 | -0.0430060 |
| sal + phosphate | Indo-Pacific | (27.8,27.9] | 0.7674341 | 0.7637779 | 0.7707190 |
| sal + phosphate | Indo-Pacific | (27.9,28.1] | 0.8204508 | 0.8075104 | 0.8333058 |
| 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.8] | -0.8741855 | -0.8842990 | -0.8687856 |
| sal + phosphate_star | Atlantic | (26.8,27.2] | -0.7948649 | -0.8077611 | -0.7775857 |
| sal + phosphate_star | Atlantic | (27.2,27.8] | -0.6845440 | -0.7101269 | -0.6589648 |
| sal + phosphate_star | Atlantic | (27.8,27.9] | -0.8578842 | -0.8661396 | -0.8470126 |
| sal + phosphate_star | Atlantic | (27.9,28.1] | -0.9557947 | -0.9587267 | -0.9518122 |
| sal + phosphate_star | Atlantic | (28.1,28.2] | -0.9699695 | -0.9721713 | -0.9685146 |
| 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.8] | -0.8266018 | -0.8346711 | -0.8211797 |
| sal + phosphate_star | Indo-Pacific | (26.8,27.2] | -0.8316525 | -0.8423046 | -0.8255105 |
| sal + phosphate_star | Indo-Pacific | (27.2,27.8] | -0.7734578 | -0.7823279 | -0.7574667 |
| sal + phosphate_star | Indo-Pacific | (27.8,27.9] | -0.7833688 | -0.7892236 | -0.7742338 |
| sal + phosphate_star | Indo-Pacific | (27.9,28.1] | -0.6363064 | -0.6491430 | -0.6265295 |
| 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.8] | -0.0327494 | -0.0336340 | -0.0310878 |
| sal + silicate | Atlantic | (26.8,27.2] | -0.1503560 | -0.1676930 | -0.1217598 |
| sal + silicate | Atlantic | (27.2,27.8] | -0.1885065 | -0.2238839 | -0.1484041 |
| sal + silicate | Atlantic | (27.8,27.9] | -0.7195266 | -0.7454738 | -0.7041747 |
| sal + silicate | Atlantic | (27.9,28.1] | -0.8838302 | -0.8875465 | -0.8771131 |
| sal + silicate | Atlantic | (28.1,28.2] | -0.7858113 | -0.7903575 | -0.7799331 |
| 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.8] | -0.5959117 | -0.6004804 | -0.5928497 |
| sal + silicate | Indo-Pacific | (26.8,27.2] | -0.6956959 | -0.7051174 | -0.6861641 |
| sal + silicate | Indo-Pacific | (27.2,27.8] | -0.1617373 | -0.2127781 | -0.0968063 |
| sal + silicate | Indo-Pacific | (27.8,27.9] | 0.7556772 | 0.7512383 | 0.7598528 |
| sal + silicate | Indo-Pacific | (27.9,28.1] | 0.8399261 | 0.8280678 | 0.8508480 |
| 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.8] | 0.9573163 | 0.9559284 | 0.9584022 |
| sal + temp | Atlantic | (26.8,27.2] | 0.9437031 | 0.9431706 | 0.9440136 |
| sal + temp | Atlantic | (27.2,27.8] | 0.9530933 | 0.9521472 | 0.9539017 |
| sal + temp | Atlantic | (27.8,27.9] | 0.9558762 | 0.9539225 | 0.9572385 |
| sal + temp | Atlantic | (27.9,28.1] | 0.9681062 | 0.9663443 | 0.9690705 |
| sal + temp | Atlantic | (28.1,28.2] | 0.9735150 | 0.9724204 | 0.9747827 |
| 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.8] | 0.9136326 | 0.9123009 | 0.9147414 |
| sal + temp | Indo-Pacific | (26.8,27.2] | 0.9414565 | 0.9398311 | 0.9433111 |
| sal + temp | Indo-Pacific | (27.2,27.8] | 0.8318672 | 0.8237183 | 0.8368833 |
| sal + temp | Indo-Pacific | (27.8,27.9] | 0.8852565 | 0.8826630 | 0.8868695 |
| sal + temp | Indo-Pacific | (27.9,28.1] | 0.8998792 | 0.8992322 | 0.9008397 |
| 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.8] | 0.4985870 | 0.4495776 | 0.5450383 |
| silicate + phosphate | Atlantic | (26.8,27.2] | 0.8490438 | 0.8476862 | 0.8514328 |
| silicate + phosphate | Atlantic | (27.2,27.8] | 0.8868531 | 0.8831422 | 0.8899526 |
| silicate + phosphate | Atlantic | (27.8,27.9] | 0.9108510 | 0.9077812 | 0.9154056 |
| silicate + phosphate | Atlantic | (27.9,28.1] | 0.9767491 | 0.9764550 | 0.9770161 |
| silicate + phosphate | Atlantic | (28.1,28.2] | 0.9438765 | 0.9420624 | 0.9451378 |
| 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.8] | 0.8842304 | 0.8829995 | 0.8856502 |
| silicate + phosphate | Indo-Pacific | (26.8,27.2] | 0.9413279 | 0.9401727 | 0.9421533 |
| silicate + phosphate | Indo-Pacific | (27.2,27.8] | 0.9724652 | 0.9721217 | 0.9728375 |
| silicate + phosphate | Indo-Pacific | (27.8,27.9] | 0.9870199 | 0.9866359 | 0.9875681 |
| silicate + phosphate | Indo-Pacific | (27.9,28.1] | 0.9918442 | 0.9916019 | 0.9920298 |
| 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.8] | 0.1323452 | 0.0790735 | 0.1790904 |
| silicate + phosphate_star | Atlantic | (26.8,27.2] | 0.5521249 | 0.5317020 | 0.5654619 |
| silicate + phosphate_star | Atlantic | (27.2,27.8] | 0.7620599 | 0.7550009 | 0.7668061 |
| silicate + phosphate_star | Atlantic | (27.8,27.9] | 0.9273887 | 0.9180910 | 0.9372661 |
| silicate + phosphate_star | Atlantic | (27.9,28.1] | 0.9055970 | 0.9039666 | 0.9068609 |
| silicate + phosphate_star | Atlantic | (28.1,28.2] | 0.8761764 | 0.8721556 | 0.8793228 |
| 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.8] | 0.7996422 | 0.7969034 | 0.8019918 |
| silicate + phosphate_star | Indo-Pacific | (26.8,27.2] | 0.7826328 | 0.7680422 | 0.7959050 |
| silicate + phosphate_star | Indo-Pacific | (27.2,27.8] | 0.2651950 | 0.2233041 | 0.2994372 |
| silicate + phosphate_star | Indo-Pacific | (27.8,27.9] | -0.4150341 | -0.4231260 | -0.4064925 |
| silicate + phosphate_star | Indo-Pacific | (27.9,28.1] | -0.3200248 | -0.3560192 | -0.2891346 |
| 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.8] | -0.3150902 | -0.3225952 | -0.3091052 |
| temp + aou | Atlantic | (26.8,27.2] | 0.2502941 | 0.2405013 | 0.2598078 |
| temp + aou | Atlantic | (27.2,27.8] | 0.5056865 | 0.4964768 | 0.5216655 |
| temp + aou | Atlantic | (27.8,27.9] | 0.5205917 | 0.4547215 | 0.5681470 |
| temp + aou | Atlantic | (27.9,28.1] | 0.1829542 | 0.1476710 | 0.2020236 |
| temp + aou | Atlantic | (28.1,28.2] | -0.6767043 | -0.6973842 | -0.6523014 |
| 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.8] | -0.2621179 | -0.2704413 | -0.2546859 |
| temp + aou | Indo-Pacific | (26.8,27.2] | -0.4260895 | -0.4433525 | -0.4100984 |
| temp + aou | Indo-Pacific | (27.2,27.8] | 0.0721323 | 0.0435393 | 0.1075568 |
| temp + aou | Indo-Pacific | (27.8,27.9] | 0.7329462 | 0.7308394 | 0.7346405 |
| temp + aou | Indo-Pacific | (27.9,28.1] | 0.8129772 | 0.8039403 | 0.8232848 |
| temp + aou | Indo-Pacific | (28.1, Inf] | 0.8156381 | 0.8060443 | 0.8249069 |
| temp + nitrate | Atlantic | (-Inf,26] | -0.7485724 | -0.7713420 | -0.7270565 |
| temp + nitrate | Atlantic | (26,26.8] | -0.4221636 | -0.4262071 | -0.4175302 |
| temp + nitrate | Atlantic | (26.8,27.2] | 0.0588050 | 0.0468305 | 0.0823800 |
| temp + nitrate | Atlantic | (27.2,27.8] | 0.0770602 | 0.0493553 | 0.1153512 |
| temp + nitrate | Atlantic | (27.8,27.9] | -0.4715488 | -0.4892459 | -0.4623507 |
| temp + nitrate | Atlantic | (27.9,28.1] | -0.8355855 | -0.8397895 | -0.8314256 |
| temp + nitrate | Atlantic | (28.1,28.2] | -0.9291293 | -0.9306516 | -0.9267098 |
| temp + nitrate | Atlantic | (28.2, Inf] | -0.8773801 | -0.8960090 | -0.8542084 |
| temp + nitrate | Indo-Pacific | (-Inf,26] | -0.5290652 | -0.5329882 | -0.5259896 |
| temp + nitrate | Indo-Pacific | (26,26.8] | -0.5105291 | -0.5141808 | -0.5046244 |
| temp + nitrate | Indo-Pacific | (26.8,27.2] | -0.6292579 | -0.6407312 | -0.6207325 |
| temp + nitrate | Indo-Pacific | (27.2,27.8] | -0.2140702 | -0.2412736 | -0.1781833 |
| temp + nitrate | Indo-Pacific | (27.8,27.9] | 0.6677945 | 0.6630710 | 0.6715435 |
| temp + nitrate | Indo-Pacific | (27.9,28.1] | 0.7365298 | 0.7245377 | 0.7490766 |
| temp + nitrate | Indo-Pacific | (28.1, Inf] | 0.7475709 | 0.7333673 | 0.7606749 |
| temp + phosphate | Atlantic | (-Inf,26] | -0.6888977 | -0.7521368 | -0.6452803 |
| temp + phosphate | Atlantic | (26,26.8] | -0.5085501 | -0.5130047 | -0.5043624 |
| temp + phosphate | Atlantic | (26.8,27.2] | -0.0288139 | -0.0421264 | -0.0033763 |
| temp + phosphate | Atlantic | (27.2,27.8] | -0.0135528 | -0.0398423 | 0.0236185 |
| temp + phosphate | Atlantic | (27.8,27.9] | -0.5209357 | -0.5367813 | -0.5125955 |
| temp + phosphate | Atlantic | (27.9,28.1] | -0.8461698 | -0.8506628 | -0.8416567 |
| temp + phosphate | Atlantic | (28.1,28.2] | -0.9248799 | -0.9261595 | -0.9224223 |
| 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.8] | -0.4727513 | -0.4796501 | -0.4692211 |
| temp + phosphate | Indo-Pacific | (26.8,27.2] | -0.5625768 | -0.5767089 | -0.5508198 |
| temp + phosphate | Indo-Pacific | (27.2,27.8] | -0.1490763 | -0.1780413 | -0.1123204 |
| temp + phosphate | Indo-Pacific | (27.8,27.9] | 0.6763486 | 0.6729656 | 0.6796909 |
| temp + phosphate | Indo-Pacific | (27.9,28.1] | 0.7473008 | 0.7367247 | 0.7588515 |
| 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.8] | -0.8731185 | -0.8826163 | -0.8655100 |
| temp + phosphate_star | Atlantic | (26.8,27.2] | -0.8373901 | -0.8466739 | -0.8235777 |
| temp + phosphate_star | Atlantic | (27.2,27.8] | -0.6697622 | -0.6869074 | -0.6539194 |
| temp + phosphate_star | Atlantic | (27.8,27.9] | -0.8638434 | -0.8691584 | -0.8593255 |
| temp + phosphate_star | Atlantic | (27.9,28.1] | -0.9640940 | -0.9651850 | -0.9631882 |
| temp + phosphate_star | Atlantic | (28.1,28.2] | -0.9774569 | -0.9784928 | -0.9766705 |
| 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.8] | -0.9211015 | -0.9253808 | -0.9182789 |
| temp + phosphate_star | Indo-Pacific | (26.8,27.2] | -0.8991644 | -0.9039909 | -0.8963304 |
| temp + phosphate_star | Indo-Pacific | (27.2,27.8] | -0.9680025 | -0.9697839 | -0.9648798 |
| temp + phosphate_star | Indo-Pacific | (27.8,27.9] | -0.9031396 | -0.9087696 | -0.8968506 |
| temp + phosphate_star | Indo-Pacific | (27.9,28.1] | -0.7772454 | -0.7888782 | -0.7683667 |
| 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.8] | -0.1663769 | -0.1805461 | -0.1552982 |
| temp + silicate | Atlantic | (26.8,27.2] | -0.3129500 | -0.3251686 | -0.2886522 |
| temp + silicate | Atlantic | (27.2,27.8] | -0.2054357 | -0.2289148 | -0.1780007 |
| temp + silicate | Atlantic | (27.8,27.9] | -0.7322802 | -0.7448248 | -0.7251895 |
| temp + silicate | Atlantic | (27.9,28.1] | -0.8638862 | -0.8683628 | -0.8576709 |
| temp + silicate | Atlantic | (28.1,28.2] | -0.8054193 | -0.8094234 | -0.8001315 |
| 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.8] | -0.6474695 | -0.6484442 | -0.6460554 |
| temp + silicate | Indo-Pacific | (26.8,27.2] | -0.6803241 | -0.6902847 | -0.6713385 |
| temp + silicate | Indo-Pacific | (27.2,27.8] | -0.2169586 | -0.2438931 | -0.1814548 |
| temp + silicate | Indo-Pacific | (27.8,27.9] | 0.6393020 | 0.6356013 | 0.6420671 |
| temp + silicate | Indo-Pacific | (27.9,28.1] | 0.7499800 | 0.7394144 | 0.7602263 |
| 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 |
|---|---|---|
| 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.8] | 0.1291296 | -0.0284869 | 0.3225885 |
| aou | Atlantic | (26.8,27.2] | -0.2668855 | -0.4220509 | -0.0850814 |
| aou | Atlantic | (27.2,27.8] | -0.1979676 | -0.2744715 | -0.1161445 |
| aou | Atlantic | (27.8,27.9] | -0.1999702 | -0.2122228 | -0.1826573 |
| aou | Atlantic | (27.9,28.1] | -0.2340339 | -0.2807752 | -0.2062425 |
| aou | Atlantic | (28.1,28.2] | -0.8576547 | -0.8607996 | -0.8550606 |
| 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.8] | -0.5985479 | -0.7135437 | -0.4436437 |
| aou | Indo-Pacific | (26.8,27.2] | -0.7831511 | -0.8539798 | -0.6919627 |
| aou | Indo-Pacific | (27.2,27.8] | -0.4196797 | -0.6825199 | -0.0602601 |
| aou | Indo-Pacific | (27.8,27.9] | 0.0214698 | -0.5144762 | 0.6761610 |
| aou | Indo-Pacific | (27.9,28.1] | 0.0024923 | -0.3849993 | 0.4692430 |
| 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.8] | 0.0146606 | -0.1501680 | 0.2166135 |
| nitrate | Atlantic | (26.8,27.2] | -0.4677631 | -0.6021762 | -0.3064964 |
| nitrate | Atlantic | (27.2,27.8] | -0.8062834 | -0.8432769 | -0.7608232 |
| nitrate | Atlantic | (27.8,27.9] | -0.9472903 | -0.9491647 | -0.9437464 |
| nitrate | Atlantic | (27.9,28.1] | -0.9660419 | -0.9679699 | -0.9638472 |
| nitrate | Atlantic | (28.1,28.2] | -0.9920564 | -0.9942188 | -0.9896660 |
| 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.8] | -0.5427093 | -0.6525888 | -0.4015306 |
| nitrate | Indo-Pacific | (26.8,27.2] | -0.7279482 | -0.7841596 | -0.6505054 |
| nitrate | Indo-Pacific | (27.2,27.8] | -0.5582126 | -0.7792233 | -0.2469751 |
| nitrate | Indo-Pacific | (27.8,27.9] | 0.0609063 | -0.4704748 | 0.6994888 |
| nitrate | Indo-Pacific | (27.9,28.1] | 0.0918176 | -0.2797599 | 0.5347383 |
| 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.8] | 0.0689656 | -0.0730130 | 0.2472585 |
| phosphate | Atlantic | (26.8,27.2] | -0.5021491 | -0.6236571 | -0.3560939 |
| phosphate | Atlantic | (27.2,27.8] | -0.8278075 | -0.8581639 | -0.7900335 |
| phosphate | Atlantic | (27.8,27.9] | -0.9573830 | -0.9595211 | -0.9560531 |
| phosphate | Atlantic | (27.9,28.1] | -0.9677872 | -0.9692877 | -0.9648102 |
| phosphate | Atlantic | (28.1,28.2] | -0.9929932 | -0.9937081 | -0.9918837 |
| 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.8] | -0.6775231 | -0.7838261 | -0.5318193 |
| phosphate | Indo-Pacific | (26.8,27.2] | -0.8305531 | -0.8838491 | -0.7599937 |
| phosphate | Indo-Pacific | (27.2,27.8] | -0.5573091 | -0.7818953 | -0.2423390 |
| phosphate | Indo-Pacific | (27.8,27.9] | 0.0616879 | -0.4647111 | 0.6962540 |
| phosphate | Indo-Pacific | (27.9,28.1] | 0.0827574 | -0.2874035 | 0.5235858 |
| 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.8] | 0.0330904 | -0.0022310 | 0.0865298 |
| phosphate_star | Atlantic | (26.8,27.2] | -0.5932136 | -0.6365792 | -0.5443498 |
| phosphate_star | Atlantic | (27.2,27.8] | -0.7051060 | -0.7466127 | -0.6609054 |
| phosphate_star | Atlantic | (27.8,27.9] | -0.8896995 | -0.9250607 | -0.8512384 |
| phosphate_star | Atlantic | (27.9,28.1] | -0.8983374 | -0.9256060 | -0.8701340 |
| phosphate_star | Atlantic | (28.1,28.2] | -0.9693399 | -0.9816035 | -0.9565821 |
| 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.8] | -0.5330718 | -0.5811654 | -0.4558077 |
| phosphate_star | Indo-Pacific | (26.8,27.2] | -0.5910983 | -0.6090677 | -0.5692514 |
| phosphate_star | Indo-Pacific | (27.2,27.8] | -0.6162277 | -0.7438081 | -0.5070049 |
| phosphate_star | Indo-Pacific | (27.8,27.9] | 0.3544758 | -0.0978458 | 0.6722219 |
| phosphate_star | Indo-Pacific | (27.9,28.1] | 0.5507426 | 0.2422048 | 0.7762409 |
| 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.8] | -0.2353248 | -0.2409769 | -0.2306188 |
| sal | Atlantic | (26.8,27.2] | 0.0705823 | -0.0214882 | 0.1623423 |
| sal | Atlantic | (27.2,27.8] | -0.0117092 | -0.1060735 | 0.0833616 |
| sal | Atlantic | (27.8,27.9] | 0.5530231 | 0.4838896 | 0.6324147 |
| sal | Atlantic | (27.9,28.1] | 0.8094648 | 0.7699275 | 0.8459362 |
| sal | Atlantic | (28.1,28.2] | 0.9258609 | 0.9067790 | 0.9479282 |
| 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.8] | 0.2894412 | 0.2524156 | 0.3177122 |
| sal | Indo-Pacific | (26.8,27.2] | 0.3082444 | 0.2511465 | 0.3628761 |
| sal | Indo-Pacific | (27.2,27.8] | 0.2716481 | 0.2138606 | 0.3463489 |
| sal | Indo-Pacific | (27.8,27.9] | -0.2735663 | -0.7410812 | 0.3331410 |
| sal | Indo-Pacific | (27.9,28.1] | -0.2904522 | -0.6378998 | 0.1508245 |
| 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.8] | -0.2548347 | -0.3780540 | -0.0977889 |
| silicate | Atlantic | (26.8,27.2] | -0.7154797 | -0.8094148 | -0.5941910 |
| silicate | Atlantic | (27.2,27.8] | -0.8844360 | -0.8917241 | -0.8738356 |
| silicate | Atlantic | (27.8,27.9] | -0.9040968 | -0.9197067 | -0.8854672 |
| silicate | Atlantic | (27.9,28.1] | -0.9284604 | -0.9303925 | -0.9260989 |
| silicate | Atlantic | (28.1,28.2] | -0.9319784 | -0.9410553 | -0.9201828 |
| 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.8] | -0.6832764 | -0.7695118 | -0.5636973 |
| silicate | Indo-Pacific | (26.8,27.2] | -0.6833458 | -0.7485095 | -0.5998382 |
| silicate | Indo-Pacific | (27.2,27.8] | -0.4884808 | -0.7127864 | -0.1780516 |
| silicate | Indo-Pacific | (27.8,27.9] | 0.0783497 | -0.4505350 | 0.7106207 |
| silicate | Indo-Pacific | (27.9,28.1] | 0.0833938 | -0.2925056 | 0.5335602 |
| 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.8] | -0.3188047 | -0.3575635 | -0.2835453 |
| temp | Atlantic | (26.8,27.2] | 0.1818015 | 0.1281887 | 0.2268658 |
| temp | Atlantic | (27.2,27.8] | -0.0253395 | -0.0919381 | 0.0413942 |
| temp | Atlantic | (27.8,27.9] | 0.5591348 | 0.5042833 | 0.6184287 |
| temp | Atlantic | (27.9,28.1] | 0.7869879 | 0.7477951 | 0.8230029 |
| temp | Atlantic | (28.1,28.2] | 0.9122250 | 0.8918712 | 0.9330506 |
| 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.8] | 0.2451831 | 0.1709819 | 0.2962770 |
| temp | Indo-Pacific | (26.8,27.2] | 0.2784718 | 0.2302637 | 0.3247969 |
| temp | Indo-Pacific | (27.2,27.8] | 0.5251279 | 0.4167914 | 0.6500110 |
| temp | Indo-Pacific | (27.8,27.9] | -0.2768759 | -0.6637713 | 0.2378509 |
| temp | Indo-Pacific | (27.9,28.1] | -0.3828356 | -0.7049990 | 0.0352371 |
| 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