eMLR - assumption testing

workflowr

• Summary
• Checks
• Past versions

Last updated: 2021-08-09

Checks: 7 0

Knit directory: emlr_obs_v_XXX/

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

---

R Markdown file: up-to-date

Great! Since the R Markdown file has been committed to the Git repository, you know the exact version of the code that produced these results.

Environment: empty

Great job! The global environment was empty. Objects defined in the global environment can affect the analysis in your R Markdown file in unknown ways. For reproduciblity it’s best to always run the code in an empty environment.

Seed: set.seed(20200707)

The command set.seed(20200707) was run prior to running the code in the R Markdown file. Setting a seed ensures that any results that rely on randomness, e.g. subsampling or permutations, are reproducible.

Session information: recorded

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

Cache: none

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

File paths: relative

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

Repository version: dd6419c

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 dd6419c. See the Past versions tab to see a history of the changes made to the R Markdown and HTML files.

Note that you need to be careful to ensure that all relevant files for the analysis have been committed to Git prior to generating the results (you can use wflow_publish or wflow_git_commit). workflowr only checks the R Markdown file, but you know if there are other scripts or data files that it depends on. Below is the status of the Git repository when the results were generated:

Ignored files:
    Ignored:    .Rhistory
    Ignored:    .Rproj.user/

Unstaged changes:
    Modified:   code/Workflowr_project_managment.R
    Modified:   data/auxillary/params_local.rds

Note that any generated files, e.g. HTML, png, CSS, etc., are not included in this status report because it is ok for generated content to have uncommitted changes.

---

These are the previous versions of the repository in which changes were made to the R Markdown (analysis/eMLR_assumption_testing.Rmd) and HTML (docs/eMLR_assumption_testing.html) files. If you’ve configured a remote Git repository (see ?wflow_git_remote), click on the hyperlinks in the table below to view the files as they were in that past version.

File	Version	Author	Date	Message
html	cd8e0d5	jens-daniel-mueller	2021-08-06	Build site.
html	15773a0	jens-daniel-mueller	2021-08-06	included calculation of revelle factor
html	da61d1a	jens-daniel-mueller	2021-08-06	Build site.
html	340d731	jens-daniel-mueller	2021-08-06	Build site.
html	71546e4	jens-daniel-mueller	2021-08-06	test with stricter CANYON-B filtering
html	29444a1	jens-daniel-mueller	2021-08-05	Build site.
html	42e80c0	jens-daniel-mueller	2021-08-04	Build site.
html	48f6eed	jens-daniel-mueller	2021-08-04	Build site.
html	1c597ab	jens-daniel-mueller	2021-08-04	Build site.
html	81a46a4	jens-daniel-mueller	2021-08-03	Build site.
html	b88c61b	jens-daniel-mueller	2021-08-03	Build site.
html	0f0d5e5	jens-daniel-mueller	2021-08-03	Build site.
html	a53656d	jens-daniel-mueller	2021-08-03	Build site.
html	88f7356	jens-daniel-mueller	2021-08-02	Build site.
html	d759279	jens-daniel-mueller	2021-08-02	Build site.
html	127b801	jens-daniel-mueller	2021-07-24	Build site.
html	1f65ef1	jens-daniel-mueller	2021-07-23	Build site.
html	912d90e	jens-daniel-mueller	2021-07-23	Build site.
html	2477316	jens-daniel-mueller	2021-07-23	rebuild: surface dcant mapping seperate
html	c9ccc00	jens-daniel-mueller	2021-07-22	Build site.
html	426b2df	jens-daniel-mueller	2021-07-21	Build site.
html	971ce87	jens-daniel-mueller	2021-07-13	Build site.
html	c18a8b1	jens-daniel-mueller	2021-07-09	Build site.
Rmd	795a3b0	jens-daniel-mueller	2021-07-09	complete revision
html	910d64a	jens-daniel-mueller	2021-07-02	Build site.
html	1cbf907	jens-daniel-mueller	2021-07-02	Build site.
html	955f5dc	jens-daniel-mueller	2021-06-29	Build site.
html	9480f94	jens-daniel-mueller	2021-06-29	Build site.
html	d054c39	jens-daniel-mueller	2021-06-24	Build site.
html	4fcd066	jens-daniel-mueller	2021-06-24	Build site.
html	3b6387f	jens-daniel-mueller	2021-06-18	Build site.
html	3052a6c	jens-daniel-mueller	2021-06-15	Build site.
html	3cb64bd	jens-daniel-mueller	2021-06-15	Build site.
html	c6b3da6	jens-daniel-mueller	2021-06-14	Build site.
html	48c73fc	jens-daniel-mueller	2021-06-14	Build site.
html	439ee80	jens-daniel-mueller	2021-06-11	Build site.
html	33ffcab	jens-daniel-mueller	2021-06-10	Build site.
html	7e1f407	jens-daniel-mueller	2021-06-10	Build site.
html	2cbe18c	jens-daniel-mueller	2021-06-10	added zonal mean section control plots
html	69c79d0	jens-daniel-mueller	2021-06-08	Build site.
html	1772903	jens-daniel-mueller	2021-06-07	Build site.
html	594ed9a	jens-daniel-mueller	2021-06-04	Build site.
html	db7df0e	jens-daniel-mueller	2021-06-04	rebuild without overlapping eras
html	2edc791	jens-daniel-mueller	2021-06-04	Build site.
html	207339d	jens-daniel-mueller	2021-06-03	Build site.
html	315710b	jens-daniel-mueller	2021-06-03	include anomalous changes
html	be90356	jens-daniel-mueller	2021-06-02	Build site.
html	d37a85d	jens-daniel-mueller	2021-05-31	Build site.
html	4b7a5ee	jens-daniel-mueller	2021-05-28	Build site.
html	12b455a	jens-daniel-mueller	2021-05-27	Build site.
Rmd	ee13efb	jens-daniel-mueller	2021-05-27	optional source of local params fully implemented
html	8c736a6	jens-daniel-mueller	2021-05-27	Build site.
Rmd	dc8e4e1	jens-daniel-mueller	2021-05-27	optional source of local params implemented
html	25bd183	jens-daniel-mueller	2021-05-26	Build site.
html	b79cb2d	jens-daniel-mueller	2021-05-20	Build site.
html	62bd574	jens-daniel-mueller	2021-05-20	Build site.
html	7c56c39	jens-daniel-mueller	2021-05-19	Build site.
html	0de759e	jens-daniel-mueller	2021-05-13	Build site.
html	52e7583	jens-daniel-mueller	2021-05-12	Build site.
html	969e631	jens-daniel-mueller	2021-05-12	Build site.
html	d2a83bc	jens-daniel-mueller	2021-04-16	Build site.
html	c0a47df	jens-daniel-mueller	2021-04-16	Build site.
html	50290e8	jens-daniel-mueller	2021-04-16	overlapping eras
html	a00ec94	jens-daniel-mueller	2021-04-16	Build site.
html	b6fe355	jens-daniel-mueller	2021-04-16	Build site.
html	81b7c6d	jens-daniel-mueller	2021-04-16	Build site.
html	ddec5b7	jens-daniel-mueller	2021-04-15	Build site.
html	29edae5	jens-daniel-mueller	2021-04-14	Build site.
html	099d566	jens-daniel-mueller	2021-04-14	Build site.
html	bb44686	jens-daniel-mueller	2021-04-14	Build site.
html	bf40480	jens-daniel-mueller	2021-04-13	Build site.
html	9f31fe3	jens-daniel-mueller	2021-04-13	Build site.
Rmd	80e1ee3	jens-daniel-mueller	2021-04-13	rerun post-2000 with model data
html	c98be84	jens-daniel-mueller	2021-04-13	Build site.
html	ea17c90	jens-daniel-mueller	2021-04-09	Build site.
Rmd	4d61e3d	jens-daniel-mueller	2021-04-09	included model data
html	338dd3c	jens-daniel-mueller	2021-04-09	Build site.
html	a79ca2c	jens-daniel-mueller	2021-04-09	included model data
html	7831fd3	jens-daniel-mueller	2021-04-09	Build site.
html	20b26d4	jens-daniel-mueller	2021-04-08	Build site.
html	c0895f8	jens-daniel-mueller	2021-04-07	Build site.
html	156d5b7	jens-daniel-mueller	2021-04-07	Build site.
html	eb827c9	jens-daniel-mueller	2021-04-07	Build site.
Rmd	49be8ed	jens-daniel-mueller	2021-03-26	included model data
html	ab5c8cf	jens-daniel-mueller	2021-03-24	Build site.
Rmd	da3c481	jens-daniel-mueller	2021-03-24	included model data
html	19faa17	jens-daniel-mueller	2021-03-24	Build site.
Rmd	ae3b4d7	jens-daniel-mueller	2021-03-24	included model data
html	03b6009	jens-daniel-mueller	2021-03-23	removed before copying template
html	98d4a6d	jens-daniel-mueller	2021-03-23	Build site.
html	f155edd	jens-daniel-mueller	2021-03-23	Build site.
html	380d215	jens-daniel-mueller	2021-03-21	Build site.
html	33b385b	jens-daniel-mueller	2021-03-20	Build site.
html	330dcd0	jens-daniel-mueller	2021-03-20	Build site.
html	83a13de	jens-daniel-mueller	2021-03-20	Build site.
html	cf98c6d	jens-daniel-mueller	2021-03-16	Build site.
html	a1d52ff	jens-daniel-mueller	2021-03-15	Build site.
html	0bade3b	jens-daniel-mueller	2021-03-15	Build site.
html	27c1f4b	jens-daniel-mueller	2021-03-14	Build site.
html	af75ebf	jens-daniel-mueller	2021-03-14	Build site.
html	5017709	jens-daniel-mueller	2021-03-11	Build site.
html	585b07f	jens-daniel-mueller	2021-03-11	Build site.
html	6482ed7	jens-daniel-mueller	2021-03-11	Build site.
html	85a5ed2	jens-daniel-mueller	2021-03-10	Build site.
html	00688a1	jens-daniel-mueller	2021-03-05	Build site.
html	6c0bec6	jens-daniel-mueller	2021-03-05	Build site.
html	3c2ec33	jens-daniel-mueller	2021-03-05	Build site.
html	af70b94	jens-daniel-mueller	2021-03-04	Build site.
html	86406d5	jens-daniel-mueller	2021-02-24	Build site.
html	3d3b4cc	jens-daniel-mueller	2021-02-23	Build site.
html	7b672f7	jens-daniel-mueller	2021-01-11	Build site.
html	33ba23c	jens-daniel-mueller	2021-01-07	Build site.
html	318609d	jens-daniel-mueller	2020-12-23	adapted more variable predictor selection
html	9d0b2d0	jens-daniel-mueller	2020-12-23	Build site.
html	0aa2b50	jens-daniel-mueller	2020-12-23	remove html before duplication
html	39113c3	jens-daniel-mueller	2020-12-23	Build site.
html	2886da0	jens-daniel-mueller	2020-12-19	Build site.
html	02f0ee9	jens-daniel-mueller	2020-12-18	cleaned up for copying template
html	965dba3	jens-daniel-mueller	2020-12-18	Build site.
html	5d452fe	jens-daniel-mueller	2020-12-18	Build site.
Rmd	ca65bf5	jens-daniel-mueller	2020-12-18	rebuild after final cleaning
html	7bcb4eb	jens-daniel-mueller	2020-12-18	Build site.
html	d397028	jens-daniel-mueller	2020-12-18	Build site.
html	7131186	jens-daniel-mueller	2020-12-17	Build site.
html	22b07fb	jens-daniel-mueller	2020-12-17	Build site.
html	f3a708f	jens-daniel-mueller	2020-12-17	Build site.
html	e4ca289	jens-daniel-mueller	2020-12-16	Build site.
html	158fe26	jens-daniel-mueller	2020-12-15	Build site.
html	7a9a4cb	jens-daniel-mueller	2020-12-15	Build site.
html	61b263c	jens-daniel-mueller	2020-12-15	Build site.
html	4d612dd	jens-daniel-mueller	2020-12-15	Build site.
Rmd	e7e5ff1	jens-daniel-mueller	2020-12-15	rebuild with eMLR target variable selection
html	953caf3	jens-daniel-mueller	2020-12-15	Build site.
html	42daf5c	jens-daniel-mueller	2020-12-14	Build site.
Rmd	923aa7f	jens-daniel-mueller	2020-12-14	rebuild with new path and auto folder creation
html	984697e	jens-daniel-mueller	2020-12-12	Build site.
html	3ebff89	jens-daniel-mueller	2020-12-12	Build site.
html	5d96d3c	jens-daniel-mueller	2020-12-11	Build site.
Rmd	3d33a37	jens-daniel-mueller	2020-12-11	selectable basinmask, try 5
html	b01a367	jens-daniel-mueller	2020-12-09	Build site.
html	24a632f	jens-daniel-mueller	2020-12-07	Build site.
html	92dca91	jens-daniel-mueller	2020-12-07	Build site.
html	6a8004b	jens-daniel-mueller	2020-12-07	Build site.
html	70bf1a5	jens-daniel-mueller	2020-12-07	Build site.
html	7555355	jens-daniel-mueller	2020-12-07	Build site.
html	143d6fa	jens-daniel-mueller	2020-12-07	Build site.
Rmd	33b1973	jens-daniel-mueller	2020-12-07	run with WOCE flag 2 only
html	abc6818	jens-daniel-mueller	2020-12-03	Build site.
Rmd	992ba15	jens-daniel-mueller	2020-12-03	rebuild with variable inventory depth
html	c8c2e7b	jens-daniel-mueller	2020-12-03	Build site.
Rmd	83203db	jens-daniel-mueller	2020-12-03	calculate cant with variable inventory depth
html	090e4d5	jens-daniel-mueller	2020-12-02	Build site.
html	7c25f7a	jens-daniel-mueller	2020-12-02	Build site.
html	ec8dc38	jens-daniel-mueller	2020-12-02	Build site.
html	c987de1	jens-daniel-mueller	2020-12-02	Build site.
html	f8358f8	jens-daniel-mueller	2020-12-02	Build site.
html	b03ddb8	jens-daniel-mueller	2020-12-02	Build site.
Rmd	9183e8f	jens-daniel-mueller	2020-12-02	revised assignment of era to eras
html	22d0127	jens-daniel-mueller	2020-12-01	Build site.
html	0ff728b	jens-daniel-mueller	2020-12-01	Build site.
html	91435ae	jens-daniel-mueller	2020-12-01	Build site.
Rmd	17d09be	jens-daniel-mueller	2020-12-01	auto eras naming
html	cf19652	jens-daniel-mueller	2020-11-30	Build site.
Rmd	2842970	jens-daniel-mueller	2020-11-30	cleaned for eMLR part only
html	196be51	jens-daniel-mueller	2020-11-30	Build site.
Rmd	7a4b015	jens-daniel-mueller	2020-11-30	first rebuild on ETH server
Rmd	bc61ce3	Jens Müller	2020-11-30	Initial commit
html	bc61ce3	Jens Müller	2020-11-30	Initial commit

---

1 Version ID

The results displayed on this site correspond to the Version_ID: v_XXX

2 Required data

Required are:

• cleaned and prepared GLODAPv2.2020 file

GLODAP <-
  read_csv(paste(path_version_data,
                 "GLODAPv2.2020_MLR_fitting_ready.csv",
                 sep = ""))

3 Predictor correlation

The correlation between:

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

were investigated based on:

• property-property plots and
• calculated correlation coeffcients.

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

plot_fun <- function(df){
  ggpairs(
    data = df,
    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:",
      unique(df$basin),
      "| era: all | random subsample:",
      nrow(df),
      "| data source:",
      unique(df$data_source)
    ))
}

GLODAP %>%
  slice_sample(n = 5e3) %>% 
  group_split(basin, data_source) %>%
  # head(1) %>% 
  map(plot_fun)

[[1]]

Past versions of predictor_correlation_per_basin-1.png

Version	Author	Date
cd8e0d5	jens-daniel-mueller	2021-08-06
15773a0	jens-daniel-mueller	2021-08-06
da61d1a	jens-daniel-mueller	2021-08-06
340d731	jens-daniel-mueller	2021-08-06
71546e4	jens-daniel-mueller	2021-08-06
29444a1	jens-daniel-mueller	2021-08-05
42e80c0	jens-daniel-mueller	2021-08-04
48f6eed	jens-daniel-mueller	2021-08-04
81a46a4	jens-daniel-mueller	2021-08-03
b88c61b	jens-daniel-mueller	2021-08-03
a53656d	jens-daniel-mueller	2021-08-03
88f7356	jens-daniel-mueller	2021-08-02
d759279	jens-daniel-mueller	2021-08-02
127b801	jens-daniel-mueller	2021-07-24
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
c18a8b1	jens-daniel-mueller	2021-07-09
910d64a	jens-daniel-mueller	2021-07-02
1cbf907	jens-daniel-mueller	2021-07-02
955f5dc	jens-daniel-mueller	2021-06-29
9480f94	jens-daniel-mueller	2021-06-29
d054c39	jens-daniel-mueller	2021-06-24
4fcd066	jens-daniel-mueller	2021-06-24
3b6387f	jens-daniel-mueller	2021-06-18
c6b3da6	jens-daniel-mueller	2021-06-14
439ee80	jens-daniel-mueller	2021-06-11
33ffcab	jens-daniel-mueller	2021-06-10
7e1f407	jens-daniel-mueller	2021-06-10
2cbe18c	jens-daniel-mueller	2021-06-10
69c79d0	jens-daniel-mueller	2021-06-08
1772903	jens-daniel-mueller	2021-06-07
594ed9a	jens-daniel-mueller	2021-06-04
db7df0e	jens-daniel-mueller	2021-06-04
2edc791	jens-daniel-mueller	2021-06-04
207339d	jens-daniel-mueller	2021-06-03
315710b	jens-daniel-mueller	2021-06-03
be90356	jens-daniel-mueller	2021-06-02
d37a85d	jens-daniel-mueller	2021-05-31
4b7a5ee	jens-daniel-mueller	2021-05-28
25bd183	jens-daniel-mueller	2021-05-26
62bd574	jens-daniel-mueller	2021-05-20
7c56c39	jens-daniel-mueller	2021-05-19
52e7583	jens-daniel-mueller	2021-05-12
969e631	jens-daniel-mueller	2021-05-12
d2a83bc	jens-daniel-mueller	2021-04-16
c0a47df	jens-daniel-mueller	2021-04-16
50290e8	jens-daniel-mueller	2021-04-16
a00ec94	jens-daniel-mueller	2021-04-16
b6fe355	jens-daniel-mueller	2021-04-16
81b7c6d	jens-daniel-mueller	2021-04-16
ddec5b7	jens-daniel-mueller	2021-04-15
29edae5	jens-daniel-mueller	2021-04-14
9f31fe3	jens-daniel-mueller	2021-04-13
ea17c90	jens-daniel-mueller	2021-04-09
338dd3c	jens-daniel-mueller	2021-04-09
a79ca2c	jens-daniel-mueller	2021-04-09
eb827c9	jens-daniel-mueller	2021-04-07
ab5c8cf	jens-daniel-mueller	2021-03-24
03b6009	jens-daniel-mueller	2021-03-23
98d4a6d	jens-daniel-mueller	2021-03-23
83a13de	jens-daniel-mueller	2021-03-20
cf98c6d	jens-daniel-mueller	2021-03-16
a1d52ff	jens-daniel-mueller	2021-03-15
0bade3b	jens-daniel-mueller	2021-03-15
27c1f4b	jens-daniel-mueller	2021-03-14
af75ebf	jens-daniel-mueller	2021-03-14
5017709	jens-daniel-mueller	2021-03-11
585b07f	jens-daniel-mueller	2021-03-11
85a5ed2	jens-daniel-mueller	2021-03-10
6c0bec6	jens-daniel-mueller	2021-03-05
af70b94	jens-daniel-mueller	2021-03-04
7b672f7	jens-daniel-mueller	2021-01-11
33ba23c	jens-daniel-mueller	2021-01-07
318609d	jens-daniel-mueller	2020-12-23
9d0b2d0	jens-daniel-mueller	2020-12-23
0aa2b50	jens-daniel-mueller	2020-12-23
2886da0	jens-daniel-mueller	2020-12-19
02f0ee9	jens-daniel-mueller	2020-12-18
7bcb4eb	jens-daniel-mueller	2020-12-18
158fe26	jens-daniel-mueller	2020-12-15
7a9a4cb	jens-daniel-mueller	2020-12-15
61b263c	jens-daniel-mueller	2020-12-15
4d612dd	jens-daniel-mueller	2020-12-15
984697e	jens-daniel-mueller	2020-12-12
3ebff89	jens-daniel-mueller	2020-12-12
5d96d3c	jens-daniel-mueller	2020-12-11

[[2]]

Past versions of predictor_correlation_per_basin-2.png

Version	Author	Date
cd8e0d5	jens-daniel-mueller	2021-08-06
15773a0	jens-daniel-mueller	2021-08-06
da61d1a	jens-daniel-mueller	2021-08-06
340d731	jens-daniel-mueller	2021-08-06
71546e4	jens-daniel-mueller	2021-08-06
29444a1	jens-daniel-mueller	2021-08-05
42e80c0	jens-daniel-mueller	2021-08-04
48f6eed	jens-daniel-mueller	2021-08-04
81a46a4	jens-daniel-mueller	2021-08-03
b88c61b	jens-daniel-mueller	2021-08-03
a53656d	jens-daniel-mueller	2021-08-03
88f7356	jens-daniel-mueller	2021-08-02
d759279	jens-daniel-mueller	2021-08-02
127b801	jens-daniel-mueller	2021-07-24
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
c18a8b1	jens-daniel-mueller	2021-07-09
910d64a	jens-daniel-mueller	2021-07-02
1cbf907	jens-daniel-mueller	2021-07-02
955f5dc	jens-daniel-mueller	2021-06-29
9480f94	jens-daniel-mueller	2021-06-29
d054c39	jens-daniel-mueller	2021-06-24
4fcd066	jens-daniel-mueller	2021-06-24
3b6387f	jens-daniel-mueller	2021-06-18
c6b3da6	jens-daniel-mueller	2021-06-14
439ee80	jens-daniel-mueller	2021-06-11
33ffcab	jens-daniel-mueller	2021-06-10
7e1f407	jens-daniel-mueller	2021-06-10
2cbe18c	jens-daniel-mueller	2021-06-10
69c79d0	jens-daniel-mueller	2021-06-08
1772903	jens-daniel-mueller	2021-06-07
594ed9a	jens-daniel-mueller	2021-06-04
db7df0e	jens-daniel-mueller	2021-06-04
2edc791	jens-daniel-mueller	2021-06-04
207339d	jens-daniel-mueller	2021-06-03
315710b	jens-daniel-mueller	2021-06-03
d37a85d	jens-daniel-mueller	2021-05-31
4b7a5ee	jens-daniel-mueller	2021-05-28
25bd183	jens-daniel-mueller	2021-05-26
62bd574	jens-daniel-mueller	2021-05-20
7c56c39	jens-daniel-mueller	2021-05-19
52e7583	jens-daniel-mueller	2021-05-12
969e631	jens-daniel-mueller	2021-05-12
d2a83bc	jens-daniel-mueller	2021-04-16
c0a47df	jens-daniel-mueller	2021-04-16
50290e8	jens-daniel-mueller	2021-04-16
a00ec94	jens-daniel-mueller	2021-04-16
b6fe355	jens-daniel-mueller	2021-04-16
81b7c6d	jens-daniel-mueller	2021-04-16
ddec5b7	jens-daniel-mueller	2021-04-15
29edae5	jens-daniel-mueller	2021-04-14
9f31fe3	jens-daniel-mueller	2021-04-13
ea17c90	jens-daniel-mueller	2021-04-09
338dd3c	jens-daniel-mueller	2021-04-09
a79ca2c	jens-daniel-mueller	2021-04-09
eb827c9	jens-daniel-mueller	2021-04-07
ab5c8cf	jens-daniel-mueller	2021-03-24
03b6009	jens-daniel-mueller	2021-03-23
98d4a6d	jens-daniel-mueller	2021-03-23
83a13de	jens-daniel-mueller	2021-03-20
cf98c6d	jens-daniel-mueller	2021-03-16
a1d52ff	jens-daniel-mueller	2021-03-15
0bade3b	jens-daniel-mueller	2021-03-15
27c1f4b	jens-daniel-mueller	2021-03-14
af75ebf	jens-daniel-mueller	2021-03-14
5017709	jens-daniel-mueller	2021-03-11
585b07f	jens-daniel-mueller	2021-03-11
85a5ed2	jens-daniel-mueller	2021-03-10
6c0bec6	jens-daniel-mueller	2021-03-05
af70b94	jens-daniel-mueller	2021-03-04
7b672f7	jens-daniel-mueller	2021-01-11
33ba23c	jens-daniel-mueller	2021-01-07
318609d	jens-daniel-mueller	2020-12-23
9d0b2d0	jens-daniel-mueller	2020-12-23
0aa2b50	jens-daniel-mueller	2020-12-23
2886da0	jens-daniel-mueller	2020-12-19
02f0ee9	jens-daniel-mueller	2020-12-18
7bcb4eb	jens-daniel-mueller	2020-12-18
158fe26	jens-daniel-mueller	2020-12-15
7a9a4cb	jens-daniel-mueller	2020-12-15
61b263c	jens-daniel-mueller	2020-12-15
4d612dd	jens-daniel-mueller	2020-12-15
984697e	jens-daniel-mueller	2020-12-12
3ebff89	jens-daniel-mueller	2020-12-12
5d96d3c	jens-daniel-mueller	2020-12-11

GLODAP <- GLODAP %>%
  filter(data_source == "obs") %>%
  select(-data_source)

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

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

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

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

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

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

}

3.2 Correlation assesment

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

3.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")

Past versions of correlation_predictor_pairs-1.png

Version	Author	Date
cd8e0d5	jens-daniel-mueller	2021-08-06
15773a0	jens-daniel-mueller	2021-08-06
da61d1a	jens-daniel-mueller	2021-08-06
340d731	jens-daniel-mueller	2021-08-06
71546e4	jens-daniel-mueller	2021-08-06
29444a1	jens-daniel-mueller	2021-08-05
42e80c0	jens-daniel-mueller	2021-08-04
48f6eed	jens-daniel-mueller	2021-08-04
81a46a4	jens-daniel-mueller	2021-08-03
b88c61b	jens-daniel-mueller	2021-08-03
a53656d	jens-daniel-mueller	2021-08-03
88f7356	jens-daniel-mueller	2021-08-02
d759279	jens-daniel-mueller	2021-08-02
127b801	jens-daniel-mueller	2021-07-24
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
910d64a	jens-daniel-mueller	2021-07-02
1cbf907	jens-daniel-mueller	2021-07-02
955f5dc	jens-daniel-mueller	2021-06-29
9480f94	jens-daniel-mueller	2021-06-29
d054c39	jens-daniel-mueller	2021-06-24
4fcd066	jens-daniel-mueller	2021-06-24
7e1f407	jens-daniel-mueller	2021-06-10
2cbe18c	jens-daniel-mueller	2021-06-10
594ed9a	jens-daniel-mueller	2021-06-04
db7df0e	jens-daniel-mueller	2021-06-04
207339d	jens-daniel-mueller	2021-06-03
315710b	jens-daniel-mueller	2021-06-03
d37a85d	jens-daniel-mueller	2021-05-31
25bd183	jens-daniel-mueller	2021-05-26
62bd574	jens-daniel-mueller	2021-05-20
7c56c39	jens-daniel-mueller	2021-05-19
52e7583	jens-daniel-mueller	2021-05-12
969e631	jens-daniel-mueller	2021-05-12
d2a83bc	jens-daniel-mueller	2021-04-16
c0a47df	jens-daniel-mueller	2021-04-16
50290e8	jens-daniel-mueller	2021-04-16
b6fe355	jens-daniel-mueller	2021-04-16
81b7c6d	jens-daniel-mueller	2021-04-16
ddec5b7	jens-daniel-mueller	2021-04-15
29edae5	jens-daniel-mueller	2021-04-14
099d566	jens-daniel-mueller	2021-04-14
bb44686	jens-daniel-mueller	2021-04-14
bf40480	jens-daniel-mueller	2021-04-13
9f31fe3	jens-daniel-mueller	2021-04-13
338dd3c	jens-daniel-mueller	2021-04-09
a79ca2c	jens-daniel-mueller	2021-04-09
eb827c9	jens-daniel-mueller	2021-04-07
19faa17	jens-daniel-mueller	2021-03-24
03b6009	jens-daniel-mueller	2021-03-23
98d4a6d	jens-daniel-mueller	2021-03-23
f155edd	jens-daniel-mueller	2021-03-23
83a13de	jens-daniel-mueller	2021-03-20
cf98c6d	jens-daniel-mueller	2021-03-16
a1d52ff	jens-daniel-mueller	2021-03-15
0bade3b	jens-daniel-mueller	2021-03-15
27c1f4b	jens-daniel-mueller	2021-03-14
af75ebf	jens-daniel-mueller	2021-03-14
5017709	jens-daniel-mueller	2021-03-11
585b07f	jens-daniel-mueller	2021-03-11
85a5ed2	jens-daniel-mueller	2021-03-10
6c0bec6	jens-daniel-mueller	2021-03-05
7b672f7	jens-daniel-mueller	2021-01-11
33ba23c	jens-daniel-mueller	2021-01-07
318609d	jens-daniel-mueller	2020-12-23
9d0b2d0	jens-daniel-mueller	2020-12-23
0aa2b50	jens-daniel-mueller	2020-12-23
2886da0	jens-daniel-mueller	2020-12-19
02f0ee9	jens-daniel-mueller	2020-12-18
158fe26	jens-daniel-mueller	2020-12-15
984697e	jens-daniel-mueller	2020-12-12
3ebff89	jens-daniel-mueller	2020-12-12
5d96d3c	jens-daniel-mueller	2020-12-11
24a632f	jens-daniel-mueller	2020-12-07
6a8004b	jens-daniel-mueller	2020-12-07
70bf1a5	jens-daniel-mueller	2020-12-07
7555355	jens-daniel-mueller	2020-12-07
143d6fa	jens-daniel-mueller	2020-12-07
0ff728b	jens-daniel-mueller	2020-12-01
91435ae	jens-daniel-mueller	2020-12-01
196be51	jens-daniel-mueller	2020-11-30
bc61ce3	Jens Müller	2020-11-30

# print table
kable(cor_predictors_stats) %>%
  add_header_above() %>%
  kable_styling() %>%
  scroll_box(width = "100%", height = "400px")

pair	basin	gamma_slab	mean_r	min_r	max_r
aou + nitrate	global	(-Inf,26]	0.9178018	0.9012309	0.9343727
aou + nitrate	global	(26,26.5]	0.9081675	0.9044419	0.9118931
aou + nitrate	global	(26.5,26.75]	0.9110977	0.9070899	0.9151055
aou + nitrate	global	(26.75,27]	0.9184310	0.9176688	0.9191932
aou + nitrate	global	(27,27.25]	0.9007433	0.8991015	0.9023850
aou + nitrate	global	(27.25,27.5]	0.8802511	0.8743037	0.8861984
aou + nitrate	global	(27.5,27.75]	0.9021361	0.8975054	0.9067669
aou + nitrate	global	(27.75,27.85]	0.9485019	0.9476013	0.9494024
aou + nitrate	global	(27.85,27.95]	0.9640856	0.9635069	0.9646644
aou + nitrate	global	(27.95,28.05]	0.9698045	0.9685493	0.9710598
aou + nitrate	global	(28.05,28.1]	0.9801019	0.9797650	0.9804389
aou + nitrate	global	(28.1,28.2]	0.9721410	0.9717692	0.9725128
aou + nitrate	global	(28.2,28.3]	0.9478889	0.9457441	0.9500337
aou + nitrate	global	(28.3,28.4]	0.3898885	0.2248323	0.5549447
aou + nitrate	global	(28.4, Inf]	0.8149304	0.7571257	0.8727351
aou + oxygen	global	(-Inf,26]	-0.9251749	-0.9444603	-0.9058896
aou + oxygen	global	(26,26.5]	-0.9635362	-0.9659030	-0.9611693
aou + oxygen	global	(26.5,26.75]	-0.9672285	-0.9682174	-0.9662397
aou + oxygen	global	(26.75,27]	-0.9719759	-0.9733693	-0.9705826
aou + oxygen	global	(27,27.25]	-0.9817002	-0.9817390	-0.9816614
aou + oxygen	global	(27.25,27.5]	-0.9847989	-0.9853036	-0.9842942
aou + oxygen	global	(27.5,27.75]	-0.9861305	-0.9867943	-0.9854668
aou + oxygen	global	(27.75,27.85]	-0.9918427	-0.9918446	-0.9918408
aou + oxygen	global	(27.85,27.95]	-0.9941138	-0.9943016	-0.9939261
aou + oxygen	global	(27.95,28.05]	-0.9950001	-0.9950957	-0.9949046
aou + oxygen	global	(28.05,28.1]	-0.9974717	-0.9975956	-0.9973477
aou + oxygen	global	(28.1,28.2]	-0.9950791	-0.9952005	-0.9949576
aou + oxygen	global	(28.2,28.3]	-0.9884561	-0.9898019	-0.9871102
aou + oxygen	global	(28.3,28.4]	-0.9973535	-0.9983196	-0.9963874
aou + oxygen	global	(28.4, Inf]	-0.9985233	-0.9985955	-0.9984512
aou + silicate	global	(-Inf,26]	0.7225558	0.6738058	0.7713057
aou + silicate	global	(26,26.5]	0.7333268	0.7255349	0.7411187
aou + silicate	global	(26.5,26.75]	0.6790251	0.6711666	0.6868837
aou + silicate	global	(26.75,27]	0.7366545	0.7274843	0.7458246
aou + silicate	global	(27,27.25]	0.8422811	0.8236554	0.8609068
aou + silicate	global	(27.25,27.5]	0.8956376	0.8809732	0.9103020
aou + silicate	global	(27.5,27.75]	0.9356841	0.9265712	0.9447971
aou + silicate	global	(27.75,27.85]	0.9749299	0.9727873	0.9770725
aou + silicate	global	(27.85,27.95]	0.9786639	0.9783521	0.9789757
aou + silicate	global	(27.95,28.05]	0.9838061	0.9836708	0.9839414
aou + silicate	global	(28.05,28.1]	0.9919029	0.9912089	0.9925970
aou + silicate	global	(28.1,28.2]	0.9812852	0.9791533	0.9834172
aou + silicate	global	(28.2,28.3]	0.9413194	0.9343511	0.9482877
aou + silicate	global	(28.3,28.4]	0.5575278	0.3853448	0.7297107
aou + silicate	global	(28.4, Inf]	0.9136217	0.8819225	0.9453208
nitrate + silicate	global	(-Inf,26]	0.8590376	0.8532402	0.8648351
nitrate + silicate	global	(26,26.5]	0.8618053	0.8546268	0.8689837
nitrate + silicate	global	(26.5,26.75]	0.8325953	0.8266031	0.8385874
nitrate + silicate	global	(26.75,27]	0.8230943	0.8086901	0.8374984
nitrate + silicate	global	(27,27.25]	0.8306752	0.8139144	0.8474360
nitrate + silicate	global	(27.25,27.5]	0.8401735	0.8305940	0.8497530
nitrate + silicate	global	(27.5,27.75]	0.8897701	0.8845433	0.8949969
nitrate + silicate	global	(27.75,27.85]	0.9504074	0.9478235	0.9529912
nitrate + silicate	global	(27.85,27.95]	0.9686553	0.9665884	0.9707223
nitrate + silicate	global	(27.95,28.05]	0.9751250	0.9744653	0.9757848
nitrate + silicate	global	(28.05,28.1]	0.9789056	0.9779279	0.9798832
nitrate + silicate	global	(28.1,28.2]	0.9827300	0.9825179	0.9829420
nitrate + silicate	global	(28.2,28.3]	0.9655737	0.9621075	0.9690399
nitrate + silicate	global	(28.3,28.4]	0.6741324	0.5932117	0.7550530
nitrate + silicate	global	(28.4, Inf]	0.9013465	0.8789942	0.9236987
oxygen + nitrate	global	(-Inf,26]	-0.7319623	-0.7944662	-0.6694585
oxygen + nitrate	global	(26,26.5]	-0.7771262	-0.7875149	-0.7667375
oxygen + nitrate	global	(26.5,26.75]	-0.7887882	-0.7919455	-0.7856310
oxygen + nitrate	global	(26.75,27]	-0.8099448	-0.8128625	-0.8070271
oxygen + nitrate	global	(27,27.25]	-0.8143641	-0.8174105	-0.8113177
oxygen + nitrate	global	(27.25,27.5]	-0.7991794	-0.8092315	-0.7891273
oxygen + nitrate	global	(27.5,27.75]	-0.8314068	-0.8392088	-0.8236047
oxygen + nitrate	global	(27.75,27.85]	-0.9074591	-0.9088904	-0.9060278
oxygen + nitrate	global	(27.85,27.95]	-0.9342536	-0.9354429	-0.9330644
oxygen + nitrate	global	(27.95,28.05]	-0.9441030	-0.9460672	-0.9421389
oxygen + nitrate	global	(28.05,28.1]	-0.9659589	-0.9666669	-0.9652509
oxygen + nitrate	global	(28.1,28.2]	-0.9491085	-0.9511024	-0.9471147
oxygen + nitrate	global	(28.2,28.3]	-0.9086696	-0.9108294	-0.9065099
oxygen + nitrate	global	(28.3,28.4]	-0.3531524	-0.5028046	-0.2035002
oxygen + nitrate	global	(28.4, Inf]	-0.8073301	-0.8771217	-0.7375385
oxygen + silicate	global	(-Inf,26]	-0.5001776	-0.6052520	-0.3951032
oxygen + silicate	global	(26,26.5]	-0.5764384	-0.5949574	-0.5579193
oxygen + silicate	global	(26.5,26.75]	-0.5057479	-0.5159968	-0.4954990
oxygen + silicate	global	(26.75,27]	-0.5908345	-0.5953242	-0.5863448
oxygen + silicate	global	(27,27.25]	-0.7565669	-0.7789494	-0.7341844
oxygen + silicate	global	(27.25,27.5]	-0.8332802	-0.8531969	-0.8133636
oxygen + silicate	global	(27.5,27.75]	-0.8861153	-0.9003300	-0.8719006
oxygen + silicate	global	(27.75,27.85]	-0.9493474	-0.9521364	-0.9465583
oxygen + silicate	global	(27.85,27.95]	-0.9573471	-0.9581120	-0.9565823
oxygen + silicate	global	(27.95,28.05]	-0.9660398	-0.9669724	-0.9651073
oxygen + silicate	global	(28.05,28.1]	-0.9843405	-0.9847743	-0.9839066
oxygen + silicate	global	(28.1,28.2]	-0.9628581	-0.9665027	-0.9592135
oxygen + silicate	global	(28.2,28.3]	-0.9073461	-0.9213661	-0.8933262
oxygen + silicate	global	(28.3,28.4]	-0.5208724	-0.6902963	-0.3514484
oxygen + silicate	global	(28.4, Inf]	-0.9030201	-0.9407477	-0.8652926
sal + aou	global	(-Inf,26]	-0.0885511	-0.1026593	-0.0744430
sal + aou	global	(26,26.5]	-0.3731743	-0.3944948	-0.3518539
sal + aou	global	(26.5,26.75]	-0.3883769	-0.3898626	-0.3868912
sal + aou	global	(26.75,27]	-0.3239785	-0.3628419	-0.2851151
sal + aou	global	(27,27.25]	-0.2321376	-0.2579643	-0.2063109
sal + aou	global	(27.25,27.5]	-0.2068642	-0.2714360	-0.1422923
sal + aou	global	(27.5,27.75]	-0.4180241	-0.4697746	-0.3662737
sal + aou	global	(27.75,27.85]	-0.6005078	-0.6005457	-0.6004699
sal + aou	global	(27.85,27.95]	-0.7069768	-0.7204966	-0.6934570
sal + aou	global	(27.95,28.05]	-0.7913593	-0.8062037	-0.7765150
sal + aou	global	(28.05,28.1]	-0.8880307	-0.8948649	-0.8811965
sal + aou	global	(28.1,28.2]	-0.8790276	-0.8821164	-0.8759388
sal + aou	global	(28.2,28.3]	-0.7576780	-0.7803583	-0.7349977
sal + aou	global	(28.3,28.4]	0.5724448	0.2396820	0.9052076
sal + aou	global	(28.4, Inf]	0.7472819	0.7144056	0.7801583
sal + nitrate	global	(-Inf,26]	-0.3436612	-0.4020244	-0.2852981
sal + nitrate	global	(26,26.5]	-0.6582791	-0.6791520	-0.6374061
sal + nitrate	global	(26.5,26.75]	-0.6765526	-0.6766854	-0.6764198
sal + nitrate	global	(26.75,27]	-0.6249800	-0.6572995	-0.5926605
sal + nitrate	global	(27,27.25]	-0.5728465	-0.5997932	-0.5458997
sal + nitrate	global	(27.25,27.5]	-0.5769425	-0.6181047	-0.5357803
sal + nitrate	global	(27.5,27.75]	-0.7083688	-0.7382589	-0.6784786
sal + nitrate	global	(27.75,27.85]	-0.7744305	-0.7808572	-0.7680037
sal + nitrate	global	(27.85,27.95]	-0.8397656	-0.8519455	-0.8275858
sal + nitrate	global	(27.95,28.05]	-0.8996204	-0.9097324	-0.8895084
sal + nitrate	global	(28.05,28.1]	-0.9492920	-0.9547598	-0.9438241
sal + nitrate	global	(28.1,28.2]	-0.9473236	-0.9524792	-0.9421679
sal + nitrate	global	(28.2,28.3]	-0.9101925	-0.9190295	-0.9013554
sal + nitrate	global	(28.3,28.4]	-0.2454841	-0.5754172	0.0844491
sal + nitrate	global	(28.4, Inf]	0.5221268	0.2617044	0.7825492
sal + oxygen	global	(-Inf,26]	-0.2083222	-0.2436785	-0.1729658
sal + oxygen	global	(26,26.5]	0.1221505	0.1077527	0.1365484
sal + oxygen	global	(26.5,26.75]	0.1477210	0.1455552	0.1498868
sal + oxygen	global	(26.75,27]	0.1001178	0.0645826	0.1356530
sal + oxygen	global	(27,27.25]	0.0490111	0.0227298	0.0752924
sal + oxygen	global	(27.25,27.5]	0.0419404	-0.0268554	0.1107362
sal + oxygen	global	(27.5,27.75]	0.2721320	0.2142202	0.3300439
sal + oxygen	global	(27.75,27.85]	0.4978849	0.4976732	0.4980965
sal + oxygen	global	(27.85,27.95]	0.6314786	0.6155079	0.6474492
sal + oxygen	global	(27.95,28.05]	0.7307152	0.7134392	0.7479912
sal + oxygen	global	(28.05,28.1]	0.8548633	0.8461974	0.8635291
sal + oxygen	global	(28.1,28.2]	0.8315531	0.8286934	0.8344129
sal + oxygen	global	(28.2,28.3]	0.6867516	0.6702182	0.7032851
sal + oxygen	global	(28.3,28.4]	-0.6038070	-0.9071371	-0.3004769
sal + oxygen	global	(28.4, Inf]	-0.7623235	-0.7886076	-0.7360394
sal + silicate	global	(-Inf,26]	-0.4937315	-0.5642770	-0.4231860
sal + silicate	global	(26,26.5]	-0.6850370	-0.7130180	-0.6570560
sal + silicate	global	(26.5,26.75]	-0.7629453	-0.7790600	-0.7468306
sal + silicate	global	(26.75,27]	-0.6912532	-0.7153167	-0.6671896
sal + silicate	global	(27,27.25]	-0.5149951	-0.5192517	-0.5107386
sal + silicate	global	(27.25,27.5]	-0.4183956	-0.4518198	-0.3849715
sal + silicate	global	(27.5,27.75]	-0.5371207	-0.5702974	-0.5039439
sal + silicate	global	(27.75,27.85]	-0.6775556	-0.6793904	-0.6757208
sal + silicate	global	(27.85,27.95]	-0.7761513	-0.7841121	-0.7681905
sal + silicate	global	(27.95,28.05]	-0.8410509	-0.8486516	-0.8334503
sal + silicate	global	(28.05,28.1]	-0.9057595	-0.9104113	-0.9011076
sal + silicate	global	(28.1,28.2]	-0.9216825	-0.9263553	-0.9170097
sal + silicate	global	(28.2,28.3]	-0.8382229	-0.8470638	-0.8293819
sal + silicate	global	(28.3,28.4]	0.0594360	-0.1333245	0.2521965
sal + silicate	global	(28.4, Inf]	0.6083910	0.4060495	0.8107325
sal + temp	global	(-Inf,26]	0.7291594	0.6977445	0.7605742
sal + temp	global	(26,26.5]	0.9677063	0.9672392	0.9681735
sal + temp	global	(26.5,26.75]	0.9870407	0.9870112	0.9870701
sal + temp	global	(26.75,27]	0.9809128	0.9807554	0.9810701
sal + temp	global	(27,27.25]	0.9744051	0.9737302	0.9750801
sal + temp	global	(27.25,27.5]	0.9605281	0.9600890	0.9609671
sal + temp	global	(27.5,27.75]	0.9454098	0.9401576	0.9506620
sal + temp	global	(27.75,27.85]	0.9790320	0.9764564	0.9816077
sal + temp	global	(27.85,27.95]	0.9547016	0.9522253	0.9571778
sal + temp	global	(27.95,28.05]	0.9309894	0.9266075	0.9353713
sal + temp	global	(28.05,28.1]	0.9401634	0.9392934	0.9410334
sal + temp	global	(28.1,28.2]	0.8097474	0.7741767	0.8453182
sal + temp	global	(28.2,28.3]	0.3498770	0.3199282	0.3798259
sal + temp	global	(28.3,28.4]	0.7009309	0.5147202	0.8871415
sal + temp	global	(28.4, Inf]	0.7448607	0.6854389	0.8042826
temp + aou	global	(-Inf,26]	-0.1512994	-0.1578146	-0.1447843
temp + aou	global	(26,26.5]	-0.3975093	-0.4178681	-0.3771506
temp + aou	global	(26.5,26.75]	-0.3698098	-0.3744261	-0.3651934
temp + aou	global	(26.75,27]	-0.3262933	-0.3660724	-0.2865143
temp + aou	global	(27,27.25]	-0.2316771	-0.2512453	-0.2121089
temp + aou	global	(27.25,27.5]	-0.1923700	-0.2548301	-0.1299099
temp + aou	global	(27.5,27.75]	-0.3694283	-0.4304853	-0.3083713
temp + aou	global	(27.75,27.85]	-0.5634464	-0.5726864	-0.5542064
temp + aou	global	(27.85,27.95]	-0.5836562	-0.6062653	-0.5610470
temp + aou	global	(27.95,28.05]	-0.5928343	-0.6214102	-0.5642584
temp + aou	global	(28.05,28.1]	-0.7097145	-0.7243766	-0.6950524
temp + aou	global	(28.1,28.2]	-0.4940891	-0.5414950	-0.4466832
temp + aou	global	(28.2,28.3]	0.0342311	-0.0058651	0.0743272
temp + aou	global	(28.3,28.4]	0.9128200	0.8848488	0.9407912
temp + aou	global	(28.4, Inf]	0.9653064	0.9641727	0.9664401
temp + nitrate	global	(-Inf,26]	-0.4308160	-0.4595457	-0.4020864
temp + nitrate	global	(26,26.5]	-0.6947572	-0.7137891	-0.6757252
temp + nitrate	global	(26.5,26.75]	-0.6737968	-0.6743516	-0.6732420
temp + nitrate	global	(26.75,27]	-0.6360868	-0.6687650	-0.6034086
temp + nitrate	global	(27,27.25]	-0.5725912	-0.5945977	-0.5505846
temp + nitrate	global	(27.25,27.5]	-0.5565825	-0.5931701	-0.5199950
temp + nitrate	global	(27.5,27.75]	-0.6578894	-0.6983933	-0.6173854
temp + nitrate	global	(27.75,27.85]	-0.7472239	-0.7501015	-0.7443464
temp + nitrate	global	(27.85,27.95]	-0.7441325	-0.7629535	-0.7253115
temp + nitrate	global	(27.95,28.05]	-0.7452313	-0.7667607	-0.7237020
temp + nitrate	global	(28.05,28.1]	-0.8160674	-0.8293678	-0.8027671
temp + nitrate	global	(28.1,28.2]	-0.6425329	-0.7002300	-0.5848358
temp + nitrate	global	(28.2,28.3]	-0.1263867	-0.1675726	-0.0852008
temp + nitrate	global	(28.3,28.4]	0.1865812	0.1044667	0.2686957
temp + nitrate	global	(28.4, Inf]	0.8357710	0.7920920	0.8794500
temp + oxygen	global	(-Inf,26]	-0.2296056	-0.2853727	-0.1738386
temp + oxygen	global	(26,26.5]	0.1382558	0.1250243	0.1514872
temp + oxygen	global	(26.5,26.75]	0.1225386	0.1214785	0.1235986
temp + oxygen	global	(26.75,27]	0.0959133	0.0598670	0.1319597
temp + oxygen	global	(27,27.25]	0.0425266	0.0226259	0.0624273
temp + oxygen	global	(27.25,27.5]	0.0195518	-0.0469951	0.0860987
temp + oxygen	global	(27.5,27.75]	0.2106935	0.1425604	0.2788266
temp + oxygen	global	(27.75,27.85]	0.4537812	0.4438099	0.4637524
temp + oxygen	global	(27.85,27.95]	0.4924697	0.4667348	0.5182047
temp + oxygen	global	(27.95,28.05]	0.5098346	0.4784055	0.5412636
temp + oxygen	global	(28.05,28.1]	0.6590148	0.6417557	0.6762739
temp + oxygen	global	(28.1,28.2]	0.4090282	0.3575879	0.4604686
temp + oxygen	global	(28.2,28.3]	-0.1764396	-0.2085699	-0.1443092
temp + oxygen	global	(28.3,28.4]	-0.9271344	-0.9457274	-0.9085415
temp + oxygen	global	(28.4, Inf]	-0.9709472	-0.9747891	-0.9671054
temp + silicate	global	(-Inf,26]	-0.5272130	-0.5787868	-0.4756393
temp + silicate	global	(26,26.5]	-0.7235496	-0.7501254	-0.6969738
temp + silicate	global	(26.5,26.75]	-0.7848380	-0.8001591	-0.7695168
temp + silicate	global	(26.75,27]	-0.7248557	-0.7517341	-0.6979773
temp + silicate	global	(27,27.25]	-0.5431633	-0.5449297	-0.5413968
temp + silicate	global	(27.25,27.5]	-0.4417687	-0.4734581	-0.4100794
temp + silicate	global	(27.5,27.75]	-0.5466940	-0.5813571	-0.5120310
temp + silicate	global	(27.75,27.85]	-0.6619160	-0.6666263	-0.6572058
temp + silicate	global	(27.85,27.95]	-0.6854326	-0.6997786	-0.6710866
temp + silicate	global	(27.95,28.05]	-0.6809790	-0.7001608	-0.6617971
temp + silicate	global	(28.05,28.1]	-0.7422450	-0.7523022	-0.7321879
temp + silicate	global	(28.1,28.2]	-0.5841102	-0.6387170	-0.5295033
temp + silicate	global	(28.2,28.3]	-0.0636541	-0.1268127	-0.0004955
temp + silicate	global	(28.3,28.4]	0.4927263	0.3843513	0.6011013
temp + silicate	global	(28.4, Inf]	0.9164134	0.9091916	0.9236352

rm(cor_predictors, cor_predictors_stats)

3.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")

Past versions of correlation_predictors-1.png

Version	Author	Date
cd8e0d5	jens-daniel-mueller	2021-08-06
15773a0	jens-daniel-mueller	2021-08-06
da61d1a	jens-daniel-mueller	2021-08-06
340d731	jens-daniel-mueller	2021-08-06
71546e4	jens-daniel-mueller	2021-08-06
29444a1	jens-daniel-mueller	2021-08-05
42e80c0	jens-daniel-mueller	2021-08-04
48f6eed	jens-daniel-mueller	2021-08-04
81a46a4	jens-daniel-mueller	2021-08-03
b88c61b	jens-daniel-mueller	2021-08-03
a53656d	jens-daniel-mueller	2021-08-03
88f7356	jens-daniel-mueller	2021-08-02
d759279	jens-daniel-mueller	2021-08-02
127b801	jens-daniel-mueller	2021-07-24
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
c18a8b1	jens-daniel-mueller	2021-07-09
910d64a	jens-daniel-mueller	2021-07-02
1cbf907	jens-daniel-mueller	2021-07-02
955f5dc	jens-daniel-mueller	2021-06-29
9480f94	jens-daniel-mueller	2021-06-29
d054c39	jens-daniel-mueller	2021-06-24
4fcd066	jens-daniel-mueller	2021-06-24
3b6387f	jens-daniel-mueller	2021-06-18
c6b3da6	jens-daniel-mueller	2021-06-14
439ee80	jens-daniel-mueller	2021-06-11
33ffcab	jens-daniel-mueller	2021-06-10
7e1f407	jens-daniel-mueller	2021-06-10
2cbe18c	jens-daniel-mueller	2021-06-10
69c79d0	jens-daniel-mueller	2021-06-08
1772903	jens-daniel-mueller	2021-06-07
594ed9a	jens-daniel-mueller	2021-06-04
db7df0e	jens-daniel-mueller	2021-06-04
2edc791	jens-daniel-mueller	2021-06-04
207339d	jens-daniel-mueller	2021-06-03
315710b	jens-daniel-mueller	2021-06-03
d37a85d	jens-daniel-mueller	2021-05-31
4b7a5ee	jens-daniel-mueller	2021-05-28
25bd183	jens-daniel-mueller	2021-05-26
62bd574	jens-daniel-mueller	2021-05-20
7c56c39	jens-daniel-mueller	2021-05-19
52e7583	jens-daniel-mueller	2021-05-12
969e631	jens-daniel-mueller	2021-05-12
d2a83bc	jens-daniel-mueller	2021-04-16
c0a47df	jens-daniel-mueller	2021-04-16
50290e8	jens-daniel-mueller	2021-04-16
a00ec94	jens-daniel-mueller	2021-04-16
b6fe355	jens-daniel-mueller	2021-04-16
81b7c6d	jens-daniel-mueller	2021-04-16
ddec5b7	jens-daniel-mueller	2021-04-15
29edae5	jens-daniel-mueller	2021-04-14
099d566	jens-daniel-mueller	2021-04-14
bb44686	jens-daniel-mueller	2021-04-14
bf40480	jens-daniel-mueller	2021-04-13
9f31fe3	jens-daniel-mueller	2021-04-13
338dd3c	jens-daniel-mueller	2021-04-09
a79ca2c	jens-daniel-mueller	2021-04-09
eb827c9	jens-daniel-mueller	2021-04-07
19faa17	jens-daniel-mueller	2021-03-24
03b6009	jens-daniel-mueller	2021-03-23
98d4a6d	jens-daniel-mueller	2021-03-23
f155edd	jens-daniel-mueller	2021-03-23
83a13de	jens-daniel-mueller	2021-03-20
cf98c6d	jens-daniel-mueller	2021-03-16
a1d52ff	jens-daniel-mueller	2021-03-15
0bade3b	jens-daniel-mueller	2021-03-15
27c1f4b	jens-daniel-mueller	2021-03-14
af75ebf	jens-daniel-mueller	2021-03-14
5017709	jens-daniel-mueller	2021-03-11
585b07f	jens-daniel-mueller	2021-03-11
85a5ed2	jens-daniel-mueller	2021-03-10
6c0bec6	jens-daniel-mueller	2021-03-05
3c2ec33	jens-daniel-mueller	2021-03-05
af70b94	jens-daniel-mueller	2021-03-04
7b672f7	jens-daniel-mueller	2021-01-11
33ba23c	jens-daniel-mueller	2021-01-07
318609d	jens-daniel-mueller	2020-12-23
9d0b2d0	jens-daniel-mueller	2020-12-23
0aa2b50	jens-daniel-mueller	2020-12-23
2886da0	jens-daniel-mueller	2020-12-19
02f0ee9	jens-daniel-mueller	2020-12-18
7bcb4eb	jens-daniel-mueller	2020-12-18
158fe26	jens-daniel-mueller	2020-12-15
7a9a4cb	jens-daniel-mueller	2020-12-15
61b263c	jens-daniel-mueller	2020-12-15
984697e	jens-daniel-mueller	2020-12-12
3ebff89	jens-daniel-mueller	2020-12-12
5d96d3c	jens-daniel-mueller	2020-12-11
24a632f	jens-daniel-mueller	2020-12-07
6a8004b	jens-daniel-mueller	2020-12-07
70bf1a5	jens-daniel-mueller	2020-12-07
7555355	jens-daniel-mueller	2020-12-07
143d6fa	jens-daniel-mueller	2020-12-07
090e4d5	jens-daniel-mueller	2020-12-02
0ff728b	jens-daniel-mueller	2020-12-01
91435ae	jens-daniel-mueller	2020-12-01
196be51	jens-daniel-mueller	2020-11-30
bc61ce3	Jens Müller	2020-11-30

# print table
kable(cor_target_predictor_stats) %>%
  add_header_above() %>%
  kable_styling() %>%
  scroll_box(width = "100%", height = "400px")

predictor	basin	gamma_slab	mean_r	min_r	max_r
aou	global	(-Inf,26]	0.0488900	0.0433129	0.0544670
aou	global	(26,26.5]	-0.1025362	-0.1877981	-0.0172742
aou	global	(26.5,26.75]	-0.1196878	-0.1704218	-0.0689538
aou	global	(26.75,27]	-0.0623009	-0.0847374	-0.0398643
aou	global	(27,27.25]	0.3374152	0.3143583	0.3604721
aou	global	(27.25,27.5]	0.4201674	0.3974180	0.4429168
aou	global	(27.5,27.75]	0.2779122	0.2212127	0.3346116
aou	global	(27.75,27.85]	0.0564590	0.0252250	0.0876931
aou	global	(27.85,27.95]	-0.0496954	-0.0668815	-0.0325092
aou	global	(27.95,28.05]	-0.1177178	-0.1389601	-0.0964755
aou	global	(28.05,28.1]	-0.3579795	-0.4191434	-0.2968156
aou	global	(28.1,28.2]	-0.6139942	-0.6626396	-0.5653487
aou	global	(28.2,28.3]	-0.8651535	-0.8777957	-0.8525112
aou	global	(28.3,28.4]	0.1262044	-0.1125515	0.3649603
aou	global	(28.4, Inf]	-0.0698988	-0.2905262	0.1507286
nitrate	global	(-Inf,26]	-0.0766444	-0.0874777	-0.0658110
nitrate	global	(26,26.5]	-0.4275834	-0.5042853	-0.3508815
nitrate	global	(26.5,26.75]	-0.4202878	-0.4598440	-0.3807315
nitrate	global	(26.75,27]	-0.3384728	-0.3562878	-0.3206577
nitrate	global	(27,27.25]	-0.0140052	-0.0296997	0.0016894
nitrate	global	(27.25,27.5]	0.0136596	-0.0096866	0.0370057
nitrate	global	(27.5,27.75]	-0.1072729	-0.1639475	-0.0505984
nitrate	global	(27.75,27.85]	-0.2198442	-0.2534860	-0.1862024
nitrate	global	(27.85,27.95]	-0.2600396	-0.2837974	-0.2362818
nitrate	global	(27.95,28.05]	-0.3014217	-0.3306051	-0.2722383
nitrate	global	(28.05,28.1]	-0.4936976	-0.5551002	-0.4322950
nitrate	global	(28.1,28.2]	-0.7394726	-0.7825317	-0.6964134
nitrate	global	(28.2,28.3]	-0.9537417	-0.9543911	-0.9530923
nitrate	global	(28.3,28.4]	-0.6642601	-0.7261054	-0.6024149
nitrate	global	(28.4, Inf]	-0.4204623	-0.5022068	-0.3387179
oxygen	global	(-Inf,26]	-0.0097804	-0.0145106	-0.0050503
oxygen	global	(26,26.5]	-0.0346052	-0.1094063	0.0401960
oxygen	global	(26.5,26.75]	0.0050674	-0.0386730	0.0488078
oxygen	global	(26.75,27]	-0.0298964	-0.0473346	-0.0124583
oxygen	global	(27,27.25]	-0.3997853	-0.4275945	-0.3719760
oxygen	global	(27.25,27.5]	-0.4803891	-0.5003317	-0.4604465
oxygen	global	(27.5,27.75]	-0.3515601	-0.4030024	-0.3001179
oxygen	global	(27.75,27.85]	-0.1218975	-0.1508324	-0.0929626
oxygen	global	(27.85,27.95]	0.0003451	-0.0147427	0.0154329
oxygen	global	(27.95,28.05]	0.0717788	0.0528788	0.0906789
oxygen	global	(28.05,28.1]	0.3296346	0.2690382	0.3902310
oxygen	global	(28.1,28.2]	0.5832975	0.5271512	0.6394439
oxygen	global	(28.2,28.3]	0.8270232	0.8144463	0.8396001
oxygen	global	(28.3,28.4]	-0.1611396	-0.3835087	0.0612295
oxygen	global	(28.4, Inf]	0.0609673	-0.1692803	0.2912149
sal	global	(-Inf,26]	0.3225546	0.3154836	0.3296256
sal	global	(26,26.5]	0.5044746	0.4317806	0.5771686
sal	global	(26.5,26.75]	0.4297000	0.3882114	0.4711885
sal	global	(26.75,27]	0.4158748	0.3988367	0.4329128
sal	global	(27,27.25]	0.3340385	0.3225387	0.3455383
sal	global	(27.25,27.5]	0.3629750	0.3607383	0.3652116
sal	global	(27.5,27.75]	0.3909749	0.3709527	0.4109971
sal	global	(27.75,27.85]	0.3932655	0.3753569	0.4111742
sal	global	(27.85,27.95]	0.3792057	0.3544555	0.4039559
sal	global	(27.95,28.05]	0.3953032	0.3691564	0.4214500
sal	global	(28.05,28.1]	0.5096476	0.4562836	0.5630116
sal	global	(28.1,28.2]	0.6870618	0.6591751	0.7149485
sal	global	(28.2,28.3]	0.8854291	0.8759054	0.8949529
sal	global	(28.3,28.4]	0.4743767	0.3739824	0.5747710
sal	global	(28.4, Inf]	0.1741940	-0.1401216	0.4885096
silicate	global	(-Inf,26]	-0.2012661	-0.2153992	-0.1871329
silicate	global	(26,26.5]	-0.4374258	-0.4919690	-0.3828826
silicate	global	(26.5,26.75]	-0.3798776	-0.4074487	-0.3523066
silicate	global	(26.75,27]	-0.0814610	-0.1051551	-0.0577669
silicate	global	(27,27.25]	0.4024840	0.3996346	0.4053335
silicate	global	(27.25,27.5]	0.4780560	0.4504111	0.5057009
silicate	global	(27.5,27.75]	0.3030964	0.2431790	0.3630138
silicate	global	(27.75,27.85]	0.0676646	0.0250422	0.1102869
silicate	global	(27.85,27.95]	-0.0387983	-0.0710705	-0.0065261
silicate	global	(27.95,28.05]	-0.1173879	-0.1470870	-0.0876889
silicate	global	(28.05,28.1]	-0.3441757	-0.4112729	-0.2770785
silicate	global	(28.1,28.2]	-0.6495028	-0.6957970	-0.6032085
silicate	global	(28.2,28.3]	-0.9057747	-0.9064284	-0.9051210
silicate	global	(28.3,28.4]	-0.3687800	-0.4741142	-0.2634458
silicate	global	(28.4, Inf]	-0.2240889	-0.2874281	-0.1607497
temp	global	(-Inf,26]	-0.1446571	-0.1618702	-0.1274441
temp	global	(26,26.5]	0.4927533	0.4155923	0.5699143
temp	global	(26.5,26.75]	0.4440702	0.4078070	0.4803334
temp	global	(26.75,27]	0.4095483	0.3982906	0.4208060
temp	global	(27,27.25]	0.3005917	0.2890447	0.3121388
temp	global	(27.25,27.5]	0.3145379	0.3070597	0.3220161
temp	global	(27.5,27.75]	0.3359671	0.3114170	0.3605171
temp	global	(27.75,27.85]	0.3857119	0.3623430	0.4090809
temp	global	(27.85,27.95]	0.4000249	0.3773910	0.4226588
temp	global	(27.95,28.05]	0.4534182	0.4266189	0.4802176
temp	global	(28.05,28.1]	0.5610956	0.5294205	0.5927707
temp	global	(28.1,28.2]	0.5629728	0.4891487	0.6367970
temp	global	(28.2,28.3]	0.1118979	0.1006318	0.1231641
temp	global	(28.3,28.4]	0.2886986	0.1181721	0.4592251
temp	global	(28.4, Inf]	-0.0849493	-0.2898858	0.1199872

rm(cor_target_predictor, cor_target_predictor_stats)

Session information

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.5      kableExtra_1.3.1 knitr_1.30      
 [5] olsrr_0.5.3      GGally_2.0.0     lubridate_1.7.9  ggforce_0.3.3   
 [9] metR_0.9.0       scico_1.2.0      patchwork_1.1.1  collapse_1.5.0  
[13] forcats_0.5.0    stringr_1.4.0    dplyr_1.0.5      purrr_0.3.4     
[17] readr_1.4.0      tidyr_1.1.2      tibble_3.0.4     ggplot2_3.3.3   
[21] tidyverse_1.3.0  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_1.4-1         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.1.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.11          farver_2.0.3            
[37] generics_0.0.2           jsonlite_1.7.1           zip_2.1.1               
[40] car_3.0-10               magrittr_1.5             Matrix_1.2-18           
[43] Rcpp_1.0.5               munsell_0.5.0            fansi_0.4.1             
[46] abind_1.4-5              lifecycle_1.0.0          stringi_1.5.3           
[49] whisker_0.4              yaml_2.2.1               carData_3.0-4           
[52] MASS_7.3-53              plyr_1.8.6               grid_4.0.3              
[55] blob_1.2.1               parallel_4.0.3           promises_1.1.1          
[58] crayon_1.3.4             lattice_0.20-41          haven_2.3.1             
[61] hms_0.5.3                pillar_1.4.7             reprex_0.3.0            
[64] glue_1.4.2               evaluate_0.14            RcppArmadillo_0.10.1.2.0
[67] data.table_1.13.2        modelr_0.1.8             vctrs_0.3.5             
[70] tweenr_1.0.2             httpuv_1.5.4             cellranger_1.1.0        
[73] gtable_0.3.0             polyclip_1.10-0          reshape_0.8.8           
[76] assertthat_0.2.1         openxlsx_4.2.3           xfun_0.18               
[79] RcppEigen_0.3.3.7.0      later_1.1.0.1            viridisLite_0.3.0       
[82] ellipsis_0.3.1           here_0.1