eMLR - model fitting

workflowr

• Summary
• Checks
• Past versions

Last updated: 2021-07-24

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: 414e9ff

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 414e9ff. 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_model_fitting.Rmd) and HTML (docs/eMLR_model_fitting.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	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.
Rmd	89d9fcb	jens-daniel-mueller	2021-07-13	complete revision
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.
Rmd	af271b5	jens-daniel-mueller	2021-06-29	plot spatial residual ditribution
html	147cb07	jens-daniel-mueller	2021-06-29	Build site.
Rmd	80ae560	jens-daniel-mueller	2021-06-29	plot spatial residual ditribution
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.
Rmd	a61014f	jens-daniel-mueller	2021-05-28	test run with cstar target variable
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.
Rmd	e2b3434	jens-daniel-mueller	2021-05-26	generate seperate output
html	b79cb2d	jens-daniel-mueller	2021-05-20	Build site.
Rmd	8adf952	jens-daniel-mueller	2021-05-20	parallelisation tested, not succesful
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.
html	c98be84	jens-daniel-mueller	2021-04-13	Build site.
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.
html	d3c3fde	jens-daniel-mueller	2021-03-26	Build site.
Rmd	88fa519	jens-daniel-mueller	2021-03-26	included model data
html	594d503	jens-daniel-mueller	2021-03-26	Build site.
Rmd	e4ef0f6	jens-daniel-mueller	2021-03-26	included model data
html	06bfaf7	jens-daniel-mueller	2021-03-26	Build site.
Rmd	6d2102a	jens-daniel-mueller	2021-03-26	included model data
html	2eb881a	jens-daniel-mueller	2021-03-26	Build site.
Rmd	a9e976a	jens-daniel-mueller	2021-03-26	included model data
html	03b6009	jens-daniel-mueller	2021-03-23	removed before copying template
html	de6256b	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.
Rmd	6c5618b	jens-daniel-mueller	2021-03-20	rerun with vif threshold
html	108d6b2	jens-daniel-mueller	2021-03-20	Build site.
Rmd	2b99ba0	jens-daniel-mueller	2021-03-20	included vif removel option
html	330dcd0	jens-daniel-mueller	2021-03-20	Build site.
Rmd	4df8990	jens-daniel-mueller	2021-03-20	added vif calculation and plots
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.
Rmd	c9cf1fd	jens-daniel-mueller	2021-03-04	rebuild with NA in Cant replaced by 0
html	27ae473	jens-daniel-mueller	2021-02-24	Build site.
Rmd	7f77d91	jens-daniel-mueller	2021-02-24	removed log10 color scale
html	fec3558	jens-daniel-mueller	2021-02-24	Build site.
Rmd	9ebedac	jens-daniel-mueller	2021-02-24	latitudinal residual plots
html	4bc00ea	jens-daniel-mueller	2021-02-24	Build site.
Rmd	de11bfe	jens-daniel-mueller	2021-02-24	clean up purrr approach and residual plots
html	42eca5d	jens-daniel-mueller	2021-02-24	Build site.
Rmd	06a2f3b	jens-daniel-mueller	2021-02-24	purrr residual plots by basin
html	a1ba577	jens-daniel-mueller	2021-02-24	Build site.
Rmd	9ae7d87	jens-daniel-mueller	2021-02-24	loop residual plots by basin
html	071743d	jens-daniel-mueller	2021-02-24	Build site.
Rmd	c45672c	jens-daniel-mueller	2021-02-24	added residual plots
html	ac1a836	jens-daniel-mueller	2021-02-24	Build site.
Rmd	5f655e0	jens-daniel-mueller	2021-02-24	added plots back to after switching to map aproach
html	b03fbd3	jens-daniel-mueller	2021-02-24	Build site.
Rmd	c69736b	jens-daniel-mueller	2021-02-24	added plots back to after switching to map aproach
html	86406d5	jens-daniel-mueller	2021-02-24	Build site.
Rmd	1b3c171	jens-daniel-mueller	2021-02-24	introduced purrr::map to model fitting, rebuild all
html	3d3b4cc	jens-daniel-mueller	2021-02-23	Build site.
Rmd	cbfc388	jens-daniel-mueller	2021-02-23	introduced purrr::map to model fitting
html	7b672f7	jens-daniel-mueller	2021-01-11	Build site.
html	33ba23c	jens-daniel-mueller	2021-01-07	Build site.
Rmd	0ad30ba	jens-daniel-mueller	2021-01-07	removed GLODAP gamma filter, target variable mapped by eras+era
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.
Rmd	bef9220	jens-daniel-mueller	2020-12-23	rebuild before sensitivity test
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.
Rmd	7e1b1c0	jens-daniel-mueller	2020-12-18	rebuild without na predictors
html	7131186	jens-daniel-mueller	2020-12-17	Build site.
Rmd	737d904	jens-daniel-mueller	2020-12-17	rebuild without na predictors
html	22b07fb	jens-daniel-mueller	2020-12-17	Build site.
html	a84ff3c	jens-daniel-mueller	2020-12-17	Build site.
Rmd	40369db	jens-daniel-mueller	2020-12-17	model selection criterion added
html	5b48ef5	jens-daniel-mueller	2020-12-17	Build site.
Rmd	e6ed2bc	jens-daniel-mueller	2020-12-17	plotted model results
html	f3a708f	jens-daniel-mueller	2020-12-17	Build site.
Rmd	7c8ace9	jens-daniel-mueller	2020-12-17	new MLR fitting routine, rmse corrected
html	e4ca289	jens-daniel-mueller	2020-12-16	Build site.
Rmd	3d5a3e2	jens-daniel-mueller	2020-12-16	new MLR fitting routine
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.
html	e91cebd	jens-daniel-mueller	2020-12-15	Build site.
Rmd	d7992c4	jens-daniel-mueller	2020-12-15	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	ba112d3	jens-daniel-mueller	2020-12-11	Build site.
Rmd	91b2b29	jens-daniel-mueller	2020-12-11	selectable basinmask, try 5
html	b01a367	jens-daniel-mueller	2020-12-09	Build site.
Rmd	71c63b0	jens-daniel-mueller	2020-12-09	rerun with variable predictor assignment
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.
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	0895ad6	jens-daniel-mueller	2020-11-30	rebuild with all plot output
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

---

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 combinations

Find all possible combinations of following considered predictor variables:

• sal, temp, aou, nitrate, silicate, phosphate, phosphate_star

# the following code is a workaround to find all predictor combinations
# using the olsrr package and fit all models for one era, slab, and basin

i_basin <- unique(GLODAP$basin)[1]
i_era   <- unique(GLODAP$era)[1]

# subset one basin and era for fitting
GLODAP_basin_era <- GLODAP %>%
  filter(basin == i_basin, era == i_era,
         data_source == "obs")

i_gamma_slab <- unique(GLODAP_basin_era$gamma_slab)[3]
print(i_gamma_slab)

# subset one gamma slab
GLODAP_basin_era_slab <- GLODAP_basin_era %>%
  filter(gamma_slab == i_gamma_slab)

# fit the full linear model, i.e. all predictor combinations
lm_full <- lm(paste(
  params_local$MLR_target,
  paste(params_local$MLR_predictors, collapse = " + "),
  sep = " ~ "
),
data = GLODAP_basin_era_slab)

# fit linear models for all possible predictor combinations
# unfortunately, this functions does not provide model coefficients (yet)
lm_all <- ols_step_all_possible(lm_full)

# convert to tibble
lm_all <- as_tibble(lm_all)

# format model formula
lm_all <- lm_all %>% 
  select(n, predictors) %>% 
  mutate(model = str_replace_all(predictors, " ", " + "),
         model = paste(params_local$MLR_target, "~", model))

# remove helper objects
rm(i_gamma_slab,
   i_era,
   i_basin,
   GLODAP_basin_era,
   GLODAP_basin_era_slab,
   lm_full)

4 Apply predictor threshold

Select combinations with a total number of predictors in the range:

• Minimum: 2
  
• Maximum: 5

lm_all <- lm_all %>% 
  filter(n >= params_local$MLR_predictors_min,
         n <= params_local$MLR_predictors_max)

This results in a total number of MLR models of:

• 112

5 Fit all models

Individual linear regression models were fitted for the chosen target variable:

• cstar_tref

as a function of each predictor combination. Fitting was performed separately within each basin, era, and slab. Model diagnostics, such as the root mean squared error (RMSE), were calculated for each fitted model.

# prepare nested data frame
GLODAP %>%
  filter_all(any_vars(is.na(.)))

GLODAP_nested <- GLODAP %>% 
  # filter(gamma_slab %in% unique(GLODAP$gamma_slab)[c(5,6)]) %>%
  nest(data = !c(gamma_slab, era, basin, data_source))

# expand with model definitions
GLODAP_nested_lm <- expand_grid(
  GLODAP_nested,
  lm_all#[1:3,]
)

# fit models and extract tidy model output
GLODAP_nested_lm_fit <- GLODAP_nested_lm %>% 
  mutate(
    fit = map2(.x = data, .y = model, 
               ~ lm(as.formula(.y), data = .x)),
    tidied = map(fit, tidy),
    glanced = map(fit, glance),
    augmented = map(fit, augment),
    vif = map(fit, ols_vif_tol)
  )

6 Tidy models

# extract glanced model output (model diagnostics, such as AIC)
GLODAP_glanced <- GLODAP_nested_lm_fit %>%
  select(-c(data, fit, tidied, augmented, vif)) %>%
  unnest(glanced) %>% 
  rename(n_predictors = n)

# extract tidy model output (model coefficients)
GLODAP_tidy <- GLODAP_nested_lm_fit %>%
  select(-c(data, fit, glanced, augmented, vif)) %>%
  unnest(tidied)

# extract augmented model output (fitted values and residuals)
GLODAP_augmented <- GLODAP_nested_lm_fit %>% 
  select(-c(data, fit, tidied, glanced, vif)) %>% 
  unnest(augmented) 

# print(object.size(GLODAP_augmented), units = "MB")

# extract VIC from output
GLODAP_glanced_vif <- GLODAP_nested_lm_fit %>%
  select(-c(data, fit, tidied, augmented, glanced)) %>%
  unnest(vif)

# calculte max vif per model
GLODAP_glanced_vif_max <- GLODAP_glanced_vif %>% 
  group_by(data_source, era, basin, gamma_slab, model, n) %>% 
  summarise(vif_max = max(VIF)) %>% 
  ungroup()

# calculate RMSE from augmented output
GLODAP_glanced_rmse <- GLODAP_augmented %>%
  group_by(data_source, era, basin, gamma_slab, model) %>%
  summarise(rmse = sqrt(c(crossprod(.resid)) / length(.resid))) %>%
  ungroup()

# add RMSE and vif_max to glanced output
GLODAP_glanced <- full_join(GLODAP_glanced, GLODAP_glanced_rmse)
GLODAP_glanced <- full_join(GLODAP_glanced, GLODAP_glanced_vif_max)
rm(GLODAP_glanced_rmse)
rm(GLODAP_glanced_vif_max)


# extract input data
GLODAP_data <- GLODAP_nested_lm_fit %>%
  select(-c(fit, tidied, glanced, augmented, vif)) %>%
  unnest(data)

# append input data with augmented data
GLODAP_augmented <- bind_cols(
  GLODAP_data,
  GLODAP_augmented %>% select(.fitted, .resid)
)

rm(GLODAP, GLODAP_nested, GLODAP_nested_lm, GLODAP_nested_lm_fit, lm_all,
   GLODAP_data)

7 Prepare coeffcients

Coefficients are prepared for the mapping of Cant and the chosen target variable.

7.1 VIF threshold

To avoid multicollinearity among predictors, models were excluded with a VIF above:

• 500

After removing models affected by multicollinearity, the targeted number of MLRs (10) was undercut in following fitting units:

GLODAP_glanced_clean <- GLODAP_glanced %>% 
  filter(vif_max <= params_local$vif_max)

GLODAP_glanced_clean %>% 
  group_by(data_source, era, basin, gamma_slab) %>% 
  count() %>% 
  filter(n < params_local$MLR_number)

# A tibble: 0 x 5
# Groups:   data_source, era, basin, gamma_slab [0]
# … with 5 variables: data_source <chr>, era <chr>, basin <chr>,
#   gamma_slab <chr>, n <int>

7.2 Predictor selection

Within each basin and slab, the following number of best linear regression models was selected:

• 10

The criterion used to select the best models was:

• rmse

The criterion was summed up for two adjacent eras, and the models with lowest summed values were selected.

# calculate RMSE sum for adjacent eras
lm_all <- GLODAP_glanced_clean  %>%
  select(data_source, basin, gamma_slab, model, era, AIC, BIC, rmse) %>% 
  arrange(era) %>% 
  group_by(data_source, basin, gamma_slab, model) %>% 
  mutate(rmse_sum = rmse + lag(rmse),
         aic_sum = AIC + lag(AIC),
         bic_sum = BIC + lag(BIC)
         ) %>% 
  ungroup() %>% 
  select(-c(era)) %>% 
  drop_na()

# subset models with lowest summed criterion
# chose which criterion is applied

if (params_local$MLR_criterion == "aic") {
  lm_best <- lm_all %>%
    group_by(data_source, basin, gamma_slab) %>%
    slice_min(order_by = aic_sum,
              with_ties = FALSE,
              n = params_local$MLR_number) %>%
    ungroup() %>%
    arrange(basin, gamma_slab, model)
} 

if (params_local$MLR_criterion == "bic") {
  lm_best <- lm_all %>%
    group_by(data_source, basin, gamma_slab) %>%
    slice_min(order_by = bic_sum,
              with_ties = FALSE,
              n = params_local$MLR_number) %>%
    ungroup() %>%
    arrange(basin, gamma_slab, model)
}

if (params_local$MLR_criterion == "rmse") {
  lm_best <- lm_all %>%
    group_by(data_source, basin, gamma_slab) %>%
    slice_min(order_by = rmse_sum,
              with_ties = FALSE,
              n = params_local$MLR_number) %>%
    ungroup() %>%
    arrange(basin, gamma_slab, model)
}


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

data_source	basin	gamma_slab	model	AIC	BIC	rmse	rmse_sum	aic_sum	bic_sum
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	2615.9403	2643.2969	8.3004227	15.634728	5604.769	5660.669
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate	2611.3697	2634.8182	8.2714831	15.812556	5622.442	5670.356
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	2611.5563	2638.9129	8.2511283	15.555104	5596.772	5652.672
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	2360.0297	2386.9142	7.3227678	10.906498	4233.816	4287.606
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	2611.4469	2638.8034	8.2499013	15.551830	5596.418	5652.318
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + silicate	2338.5113	2365.3957	7.0972795	10.837308	4241.753	4295.543
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + silicate	2612.8681	2640.2247	8.2658473	15.771996	5621.893	5677.793
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + silicate	2337.5995	2360.6434	7.1085148	10.851124	4239.317	4285.422
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + silicate + phosphate	2337.1624	2364.0469	7.0833788	10.626864	4203.156	4256.946
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + silicate + phosphate	2645.4171	2672.7737	8.6396010	16.288923	5670.917	5726.818
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	2336.6089	2363.4934	7.0776824	10.617154	4201.821	4255.610
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	2644.9557	2672.3123	8.6341868	16.278740	5669.912	5725.812
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	2610.3398	2637.6964	8.2375017	15.533733	5594.630	5650.530
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + nitrate + silicate	2339.3494	2362.3933	7.1266180	10.877864	4242.658	4288.763
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate	2334.0581	2360.9426	7.0514898	10.657239	4212.071	4265.860
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	2339.0287	2365.9132	7.1026195	10.626788	4201.251	4255.040
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + phosphate + phosphate_star	2642.8232	2666.2717	8.6326326	16.307123	5669.188	5717.102
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + silicate + phosphate	2335.5732	2358.6171	7.0876094	10.883320	4247.012	4293.117
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	2336.5975	2363.4820	7.0775648	10.615865	4201.581	4255.370
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	2644.5607	2671.9173	8.6295546	16.268569	5668.885	5724.785
obs	Atlantic	(26,26.5]	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	5961.5560	5990.3167	6.7938266	14.363559	13238.489	13297.018
mod	Atlantic	(26,26.5]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	4560.9722	4593.2753	5.0968871	10.088694	9944.403	10010.221
obs	Atlantic	(26,26.5]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	5962.0462	5995.6005	6.7880796	14.342765	13236.781	13305.064
obs	Atlantic	(26,26.5]	cstar_tref ~ aou + phosphate + phosphate_star	5965.7755	5989.7428	6.8175533	14.389345	13241.283	13290.056
obs	Atlantic	(26,26.5]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	5964.4941	5993.2549	6.8050248	14.363427	13238.267	13296.795
obs	Atlantic	(26,26.5]	cstar_tref ~ sal + aou + phosphate + phosphate_star	5966.9068	5995.6676	6.8142342	14.324487	13227.195	13285.723
obs	Atlantic	(26,26.5]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	5965.1158	5998.6701	6.7997695	14.296083	13223.484	13291.768
mod	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	4475.8047	4508.1078	4.8140908	9.375394	9699.239	9765.057
obs	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	6003.9404	6037.4947	6.9493727	14.538861	13288.373	13356.656
mod	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	4483.0762	4515.3793	4.8376103	9.451879	9726.991	9792.809
obs	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + aou + phosphate	6006.6184	6035.3792	6.9676190	14.568244	13292.145	13350.674
obs	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate	6008.5609	6042.1152	6.9673947	14.560524	13294.006	13362.289
mod	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	4492.8650	4525.1681	4.8694535	9.562667	9766.875	9832.693
mod	Atlantic	(26,26.5]	cstar_tref ~ temp + aou + nitrate + phosphate	4535.4129	4563.1013	5.0170373	9.901773	9878.378	9934.794
mod	Atlantic	(26,26.5]	cstar_tref ~ temp + aou + nitrate + phosphate_star	4531.4850	4559.1733	5.0038464	9.927180	9888.413	9944.828
mod	Atlantic	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	4536.8355	4569.1386	5.0150960	9.862269	9868.106	9933.924
obs	Atlantic	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	6005.4035	6038.9578	6.9550744	14.598014	13304.644	13372.928
mod	Atlantic	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	4533.1061	4565.4091	5.0025757	9.856386	9866.804	9932.622
mod	Atlantic	(26,26.5]	cstar_tref ~ temp + nitrate + phosphate + phosphate_star	4531.2878	4558.9762	5.0031853	9.982627	9908.319	9964.734
mod	Atlantic	(26,26.5]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	4533.1135	4565.4166	5.0026007	9.911579	9886.861	9952.679
mod	Atlantic	(26.5,26.75]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	4245.5894	4278.0777	3.8312669	6.831382	9697.582	9764.963
obs	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	6640.4171	6674.9579	6.0931501	12.768748	15047.432	15117.973
obs	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	6634.6091	6669.1499	6.0759452	12.729266	15033.167	15103.708
obs	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + aou + phosphate	6641.8750	6671.4813	6.1034165	12.803831	15056.278	15116.741
obs	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + aou + phosphate_star	6634.0856	6663.6920	6.0803144	12.774512	15046.140	15106.604
mod	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	4152.2881	4184.7763	3.6048995	6.692604	9666.439	9733.821
obs	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	6643.7885	6678.3293	6.1031597	12.795035	15056.965	15127.506
mod	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	4140.4050	4172.8933	3.5770460	6.696043	9676.337	9743.718
obs	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	6635.0920	6669.6328	6.0773738	12.764808	15046.589	15117.130
obs	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	6634.9048	6669.4456	6.0768200	12.732418	15034.329	15104.869
obs	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + phosphate + phosphate_star	6634.7034	6664.3098	6.0821436	12.780855	15048.463	15108.926
mod	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	4140.8173	4173.3055	3.5780087	6.711843	9687.000	9754.381
obs	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	6635.5090	6670.0498	6.0786079	12.771685	15049.140	15119.681
mod	Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	4142.7557	4175.2440	3.5825388	6.544793	9567.316	9634.698
obs	Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	6640.7638	6675.3046	6.0941788	12.829733	15070.401	15140.941
mod	Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	4113.2922	4145.7805	3.5142976	6.510713	9562.619	9630.001
mod	Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + silicate + phosphate	4187.5177	4215.3648	3.6935764	6.789319	9705.285	9763.040
mod	Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + silicate + phosphate_star	4167.6961	4195.5432	3.6460953	6.772184	9706.534	9764.289
mod	Atlantic	(26.5,26.75]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	4100.7365	4133.2248	3.4856134	6.507376	9568.259	9635.640
mod	Atlantic	(26.5,26.75]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	4166.1687	4194.0158	3.6424620	6.782939	9714.925	9772.681
mod	Atlantic	(26.75,27]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	5344.0834	5379.1176	2.7167318	4.944251	11933.211	12005.348
obs	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	6988.7163	7023.6738	5.9326596	11.972789	16101.422	16173.158
mod	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + aou + silicate + phosphate	5105.2292	5140.2634	2.4377039	4.561644	11553.320	11625.457
mod	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	5164.7350	5199.7692	2.5044160	4.647135	11638.898	11711.036
obs	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	7003.7950	7038.7526	5.9738371	12.076534	16145.644	16217.380
obs	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	6991.9162	7026.8737	5.9413741	11.986743	16107.074	16178.811
obs	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + nitrate + phosphate_star	7015.3895	7045.3531	6.0112063	12.095951	16146.907	16208.395
obs	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	7002.7833	7037.7408	5.9710653	12.002507	16111.418	16183.155
mod	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	5200.5280	5235.5622	2.5454198	4.705720	11698.896	11771.033
obs	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	7003.8404	7038.7979	5.9739614	12.068123	16141.731	16213.468
obs	Atlantic	(26.75,27]	cstar_tref ~ temp + aou + nitrate + phosphate_star	6995.6333	7025.5969	5.9569761	12.010737	16112.713	16174.202
mod	Atlantic	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	5020.9083	5055.9424	2.3462034	4.441362	11428.586	11500.723
mod	Atlantic	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	5003.3205	5038.3547	2.3275553	4.433868	11426.725	11498.863
obs	Atlantic	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	6995.5330	7030.4906	5.9512398	11.982396	16104.034	16175.771
mod	Atlantic	(26.75,27]	cstar_tref ~ temp + aou + silicate + phosphate	5223.8748	5253.9041	2.5748620	4.735752	11721.052	11782.884
mod	Atlantic	(26.75,27]	cstar_tref ~ temp + aou + silicate + phosphate_star	5225.6209	5255.6501	2.5769027	4.745255	11733.009	11794.841
obs	Atlantic	(26.75,27]	cstar_tref ~ temp + nitrate + phosphate + phosphate_star	6996.2565	7026.2201	5.9586793	12.015957	16114.979	16176.468
mod	Atlantic	(26.75,27]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	5049.8312	5084.8653	2.3771953	4.504585	11502.729	11574.866
obs	Atlantic	(26.75,27]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	6996.2421	7031.1996	5.9531758	11.987219	16106.096	16177.833
mod	Atlantic	(26.75,27]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	5245.2916	5275.3209	2.6000045	4.782068	11771.350	11833.182
mod	Atlantic	(27,27.25]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	6108.3909	6144.1141	2.9653016	5.118853	13200.843	13274.293
mod	Atlantic	(27,27.25]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	5942.2484	5977.9716	2.7694911	4.909998	13015.028	13088.478
mod	Atlantic	(27,27.25]	cstar_tref ~ sal + aou + silicate + phosphate	6225.8866	6256.5066	3.1146399	5.280354	13334.575	13397.532
mod	Atlantic	(27,27.25]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	6222.2297	6257.9530	3.1074035	4.878314	12681.249	12754.699
obs	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	9266.6597	9303.5373	6.0935626	12.636384	22093.995	22169.876
obs	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	9248.8084	9285.6860	6.0557523	12.566228	22056.885	22132.766
mod	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + aou + silicate + phosphate	6210.6101	6246.3333	3.0925924	4.895024	12726.757	12800.207
mod	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	6220.2727	6255.9959	3.1049040	5.123785	13103.631	13177.081
obs	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	9275.9753	9312.8529	6.1133875	12.786935	22180.188	22256.069
obs	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	9250.1728	9287.0503	6.0586338	12.572494	22060.269	22136.150
mod	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	6235.8026	6271.5259	3.1247943	5.233173	13259.612	13333.062
obs	Atlantic	(27,27.25]	cstar_tref ~ temp + aou + nitrate + phosphate	9309.9700	9341.5794	6.1905968	12.741338	22140.006	22205.047
obs	Atlantic	(27,27.25]	cstar_tref ~ temp + aou + nitrate + phosphate_star	9312.2055	9343.8149	6.1954240	12.718568	22125.836	22190.877
mod	Atlantic	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	5733.9200	5769.6432	2.5421298	4.370643	12296.586	12370.036
obs	Atlantic	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	9299.2042	9336.0818	6.1631029	12.697940	22121.794	22197.675
mod	Atlantic	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	5722.8240	5758.5472	2.5305578	4.418996	12389.905	12463.356
obs	Atlantic	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	9297.6597	9334.5372	6.1597847	12.658577	22098.756	22174.638
obs	Atlantic	(27,27.25]	cstar_tref ~ temp + nitrate + phosphate + phosphate_star	9306.0092	9337.6186	6.1820533	12.706686	22120.526	22185.567
mod	Atlantic	(27,27.25]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	5743.9502	5779.6734	2.5526359	4.567550	12620.940	12694.390
obs	Atlantic	(27,27.25]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	9292.5957	9329.4733	6.1489181	12.651862	22096.174	22172.055
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	10305.6132	10343.5572	5.2718575	10.795773	24714.787	24792.919
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	8799.3741	8837.1226	3.6179878	6.127618	18731.047	18808.409
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	10284.6928	10322.6368	5.2389400	10.744948	24678.923	24757.055
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + nitrate + silicate	8377.4937	8415.2423	3.1772904	5.308052	17615.267	17692.630
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + phosphate_star	10314.9364	10347.4599	5.2897604	10.868145	24767.265	24834.236
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate	8918.7191	8956.4677	3.7534003	5.952087	18289.547	18366.909
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	8897.7149	8935.4635	3.7292061	5.935051	18282.323	18359.686
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	10299.4706	10337.4147	5.2621711	10.836399	24750.370	24828.501
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	10287.9451	10325.8892	5.2440439	10.750732	24682.744	24760.876
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	8872.5408	8910.2893	3.7004140	5.976928	18390.856	18468.219
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + phosphate + phosphate_star	10315.4215	10347.9450	5.2905289	10.869071	24767.881	24834.851
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	8897.3086	8935.0571	3.7287396	5.943228	18298.499	18375.862
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	10302.1262	10340.0702	5.2663565	10.840926	24753.307	24831.439
mod	Atlantic	(27.25,27.5]	cstar_tref ~ temp + aou + nitrate + silicate	8429.2943	8461.6502	3.2303580	5.362939	17668.686	17734.997
mod	Atlantic	(27.25,27.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	8431.1672	8468.9157	3.2302316	5.306784	17559.674	17637.036
obs	Atlantic	(27.25,27.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	10336.3507	10374.2948	5.3205976	10.853131	24752.698	24830.830
mod	Atlantic	(27.25,27.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	8427.8029	8465.5514	3.2268874	5.250167	17446.114	17523.477
obs	Atlantic	(27.25,27.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	10317.5995	10355.5435	5.2908107	10.808212	24721.343	24799.474
mod	Atlantic	(27.25,27.5]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	8350.7599	8388.5084	3.1512458	5.073889	17152.751	17230.114
obs	Atlantic	(27.25,27.5]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	10319.4598	10357.4038	5.2937584	10.810563	24722.705	24800.837
obs	Atlantic	(27.5,27.75]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	9986.7202	10024.2543	5.7373900	11.183176	21702.308	21778.612
obs	Atlantic	(27.5,27.75]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	9991.7254	10023.8975	5.7501633	11.213065	21717.106	21782.509
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + aou + phosphate + phosphate_star	11269.2187	11303.4745	3.0229522	5.155183	23854.766	23924.834
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	10989.1921	11029.1572	2.8376347	4.784685	23051.790	23133.536
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + aou + phosphate	11367.8174	11402.0733	3.0905564	5.326140	24226.857	24296.925
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + aou + phosphate_star	11617.1124	11651.3682	3.2683065	5.697111	24955.132	25025.200
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	11085.2069	11125.1721	2.8994136	4.948309	23442.397	23524.143
obs	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	9991.2300	10028.7641	5.7456099	11.174334	21695.025	21771.329
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	11304.9352	11344.9003	3.0459019	5.268647	24132.699	24214.445
obs	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	10007.2560	10044.7901	5.7749159	11.213665	21717.984	21794.288
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	11331.4305	11371.3956	3.0640585	5.287665	24161.432	24243.178
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + phosphate + phosphate_star	11586.8744	11621.1303	3.2462131	5.666683	24905.032	24975.100
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	11288.2958	11328.2609	3.0345543	5.255536	24111.474	24193.220
obs	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	9997.9685	10035.5026	5.7579141	11.189382	21703.663	21779.966
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + silicate + phosphate_star	11459.8252	11494.0811	3.1550046	5.381257	24294.699	24364.767
obs	Atlantic	(27.5,27.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	9992.2617	10029.7957	5.7474920	11.189175	21705.017	21781.320
obs	Atlantic	(27.5,27.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	10016.1071	10053.6412	5.7911655	11.268231	21753.218	21829.522
obs	Atlantic	(27.5,27.75]	cstar_tref ~ temp + aou + silicate + phosphate	9995.5224	10027.6944	5.7570987	11.219510	21720.565	21785.968
obs	Atlantic	(27.5,27.75]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	10004.2323	10041.7664	5.7693752	11.225397	21726.877	21803.180
obs	Atlantic	(27.5,27.75]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	10005.9224	10038.0945	5.7761377	11.257756	21744.156	21809.559
mod	Atlantic	(27.75,27.85]	cstar_tref ~ aou + phosphate + phosphate_star	6729.2761	6755.4901	2.6757669	4.324776	13283.451	13336.872
mod	Atlantic	(27.75,27.85]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	6388.4553	6419.9121	2.3670037	3.893629	12681.668	12745.773
obs	Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	8052.0075	8088.1713	5.3899670	10.311012	17844.451	17918.359
obs	Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + silicate + phosphate	8050.0131	8081.0107	5.3899787	10.356664	17870.423	17933.772
obs	Atlantic	(27.75,27.85]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	8059.1990	8095.3628	5.4049537	10.347901	17866.058	17939.966
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + aou + silicate + phosphate	6671.9162	6708.6158	2.6176853	4.308940	13316.351	13391.140
obs	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + aou + silicate + phosphate	8050.7914	8086.9553	5.3874369	10.353805	17872.994	17946.902
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	6814.9842	6851.6838	2.7551152	4.581948	13722.374	13797.163
obs	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	8060.2068	8096.3706	5.4070573	10.393492	17895.498	17969.406
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	6828.5730	6865.2726	2.7685379	4.600096	13744.773	13819.562
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	6801.9100	6838.6096	2.7422622	4.565409	13702.413	13777.202
obs	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	8053.7349	8089.8987	5.3935630	10.367562	17880.921	17954.829
obs	Atlantic	(27.75,27.85]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	8072.9779	8109.1417	5.4337849	10.366973	17873.421	17947.329
obs	Atlantic	(27.75,27.85]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	8067.8590	8104.0229	5.4230562	10.357877	17869.376	17943.284
mod	Atlantic	(27.75,27.85]	cstar_tref ~ temp + aou + silicate + phosphate	6718.8901	6750.3469	2.6639398	4.407076	13464.360	13528.465
mod	Atlantic	(27.75,27.85]	cstar_tref ~ temp + aou + silicate + phosphate_star	6854.2124	6885.6692	2.7960415	4.663898	13835.332	13899.437
obs	Atlantic	(27.75,27.85]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	8068.9365	8105.1003	5.4253128	10.353834	17866.307	17940.215
mod	Atlantic	(27.75,27.85]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	6864.8487	6896.3055	2.8066982	4.678474	13853.115	13917.220
mod	Atlantic	(27.75,27.85]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	6842.8569	6874.3137	2.7847089	4.649354	13818.108	13882.213
obs	Atlantic	(27.75,27.85]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	8067.0593	8098.0569	5.4255700	10.401072	17893.226	17956.576
mod	Atlantic	(27.85,27.95]	cstar_tref ~ aou + phosphate + phosphate_star	6010.9427	6036.2107	3.2361410	5.865201	13516.387	13568.453
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	12225.9065	12264.7758	5.9572241	11.277502	27754.036	27833.702
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	12069.2251	12108.0945	5.7173332	10.894411	27460.385	27540.050
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + aou + silicate + phosphate	5986.1968	6016.5183	3.1989520	5.998742	13691.330	13753.808
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + aou + silicate + phosphate_star	12067.3970	12100.7136	5.7175911	10.945912	27506.001	27574.286
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	12151.9091	12190.7784	5.8426994	11.078242	27599.443	27679.108
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	12154.2360	12187.5525	5.8493350	11.111581	27625.308	27693.593
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate + phosphate	5919.8454	5955.2205	3.1058426	5.613260	13280.444	13353.336
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate + phosphate	12050.9784	12089.8477	5.6900317	10.911771	27485.259	27564.925
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	5954.5040	5989.8791	3.1527114	5.717133	13385.734	13458.626
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	12059.8176	12098.6869	5.7032409	10.931516	27500.378	27580.043
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate	12068.3942	12107.2635	5.7160870	10.955024	27519.181	27598.846
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	5985.5243	6020.8994	3.1952597	5.781194	13443.003	13515.894
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + silicate + phosphate	12066.3944	12099.7110	5.7160875	11.037803	27593.881	27662.165
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	5943.7765	5979.1516	3.1381296	5.692362	13362.497	13435.389
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	12051.6096	12090.4789	5.6909740	10.913357	27486.509	27566.175
mod	Atlantic	(27.85,27.95]	cstar_tref ~ temp + aou + silicate + phosphate	6032.2886	6062.6101	3.2633098	5.807641	13436.806	13499.285
mod	Atlantic	(27.85,27.95]	cstar_tref ~ temp + aou + silicate + phosphate_star	6070.8160	6101.1375	3.3180976	5.921427	13547.359	13609.838
mod	Atlantic	(27.85,27.95]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	6100.8427	6131.1642	3.3614342	5.985816	13602.692	13665.170
mod	Atlantic	(27.85,27.95]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	6059.4490	6089.7705	3.3018383	5.894582	13523.190	13585.668
mod	Atlantic	(27.95,28.05]	cstar_tref ~ aou + phosphate + phosphate_star	7306.0481	7331.5852	4.8008652	8.468245	17075.162	17128.163
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	14869.3465	14909.4022	6.4918894	11.264391	33269.181	33351.472
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	14758.0402	14798.0958	6.3338388	10.874109	32850.290	32932.582
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + aou + silicate + phosphate	7438.1098	7468.7543	5.0634925	8.655292	17134.465	17198.066
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + aou + silicate + phosphate_star	14756.6621	14790.9955	6.3347111	10.877243	32849.983	32920.518
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	14801.8391	14841.8947	6.3955671	11.003804	32985.709	33068.000
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	14799.8398	14834.1732	6.3955681	11.003814	32981.721	33052.256
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + aou + phosphate	7421.8140	7452.4586	5.0298157	8.590039	17086.469	17150.070
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate + phosphate	7348.4048	7384.1567	4.8768681	8.341066	16916.902	16991.103
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate + phosphate	14781.3288	14821.3844	6.3665862	10.937344	32914.845	32997.136
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	7421.3281	7457.0801	5.0246980	8.497706	16998.943	17073.144
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	14757.9674	14798.0231	6.3337368	10.876164	32853.147	32935.438
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate	14818.7708	14858.8264	6.4195908	11.101112	33099.927	33182.218
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + phosphate + phosphate_star	7466.1830	7496.8275	5.1220383	8.700434	17149.115	17212.716
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + silicate + phosphate	14818.8398	14853.1732	6.4225327	11.117845	33116.132	33186.668
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	7384.9304	7420.6824	4.9503610	8.422016	16961.147	17035.348
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	14775.1351	14815.1907	6.3578605	10.918188	32894.565	32976.856
mod	Atlantic	(27.95,28.05]	cstar_tref ~ temp + aou + silicate + phosphate	7347.3637	7378.0082	4.8787835	8.532946	17105.515	17169.116
mod	Atlantic	(27.95,28.05]	cstar_tref ~ temp + aou + silicate + phosphate_star	7423.4651	7454.1096	5.0332175	8.707645	17201.471	17265.072
mod	Atlantic	(27.95,28.05]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	7383.5729	7414.2174	4.9516635	8.624189	17159.720	17223.321
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	12626.6574	12665.6433	6.2657672	10.813722	28499.787	28580.087
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	12631.9143	12670.9002	6.2742710	10.787089	28463.115	28543.416
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + aou + silicate + phosphate_star	12653.0412	12686.4576	6.3118195	10.830989	28489.846	28558.675
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	12624.6058	12663.5917	6.2624517	10.772288	28452.234	28532.535
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	12643.4364	12676.8529	6.2961982	10.813326	28477.797	28546.626
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	12612.7426	12651.7285	6.2433135	10.848669	28553.675	28633.976
mod	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + aou + phosphate	3269.5510	3298.2984	1.5085461	2.870701	8085.409	8145.579
mod	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + aou + phosphate_star	3284.4113	3313.1586	1.5211928	2.897559	8129.124	8189.293
mod	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + aou + silicate + phosphate	3255.0418	3288.5804	1.4946192	2.840781	8040.067	8110.265
mod	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	3274.8989	3308.4374	1.5113860	2.876918	8099.640	8169.838
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	12592.1663	12631.1521	6.2102578	10.695546	28390.288	28470.589
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	12610.2142	12649.2001	6.2392422	10.853163	28561.192	28641.492
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	12581.7799	12620.7658	6.1936387	10.781072	28501.633	28581.934
mod	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + phosphate + phosphate_star	3279.6831	3308.4305	1.5171575	2.886539	8110.251	8170.421
mod	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	3268.1544	3301.6930	1.5056702	2.861177	8072.410	8142.608
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	12603.8270	12642.8129	6.2289691	10.739931	28432.804	28513.105
mod	Atlantic	(28.05,28.1]	cstar_tref ~ temp + aou + phosphate	3289.8131	3313.7692	1.5275316	2.923599	8172.034	8222.175
mod	Atlantic	(28.05,28.1]	cstar_tref ~ temp + aou + silicate + phosphate	3271.7098	3300.4571	1.5103767	2.876135	8094.913	8155.083
mod	Atlantic	(28.05,28.1]	cstar_tref ~ temp + aou + silicate + phosphate_star	3292.7131	3321.4604	1.5283041	2.916429	8161.074	8221.244
mod	Atlantic	(28.05,28.1]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	3286.0401	3314.7874	1.5225854	2.899296	8131.446	8191.616
obs	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	12275.4508	12313.6217	8.4579361	13.049400	26975.707	27054.632
mod	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + aou + phosphate	4370.5221	4401.5151	1.3036774	2.277270	8904.311	8967.660
mod	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + aou + phosphate_star	4378.3569	4409.3498	1.3076300	2.286788	8930.660	8994.009
mod	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + aou + silicate + phosphate	4245.0294	4281.1878	1.2410105	2.171003	8632.006	8705.913
mod	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	4266.0130	4302.1715	1.2511136	2.193876	8697.286	8771.193
obs	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	12279.2079	12317.3788	8.4671519	13.041289	26960.597	27039.522
obs	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	12272.9137	12311.0846	8.4517184	13.043137	26973.120	27052.045
mod	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	4405.8495	4442.0080	1.3205743	2.290010	8927.741	9001.648
obs	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	12269.6578	12307.8286	8.4437459	13.021740	26955.254	27034.179
mod	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + phosphate + phosphate_star	4373.8987	4404.8916	1.3053794	2.282242	8918.579	8981.928
mod	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	4252.9005	4289.0590	1.2447907	2.181714	8663.994	8737.901
obs	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	12279.2282	12317.3991	8.4672018	13.043098	26962.536	27041.461
obs	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + silicate + phosphate_star	12277.2574	12309.9753	8.4672735	13.046933	26962.668	27030.318
obs	Atlantic	(28.1,28.15]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	12266.4182	12304.5891	8.4358209	13.028560	26968.060	27046.985
mod	Atlantic	(28.1,28.15]	cstar_tref ~ temp + aou + silicate + phosphate	4266.9092	4297.9022	1.2525145	2.205978	8732.840	8796.189
mod	Atlantic	(28.1,28.15]	cstar_tref ~ temp + aou + silicate + phosphate_star	4290.8299	4321.8228	1.2641450	2.234248	8812.956	8876.305
obs	Atlantic	(28.1,28.15]	cstar_tref ~ temp + aou + silicate + phosphate_star	12277.2548	12309.9727	8.4672671	13.060545	26977.482	27045.132
obs	Atlantic	(28.1,28.15]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	12264.8469	12303.0177	8.4319796	13.026238	26968.139	27047.064
obs	Atlantic	(28.1,28.15]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	12267.8068	12300.5247	8.4441107	13.040987	26971.942	27039.592
mod	Atlantic	(28.1,28.15]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	4276.2587	4307.2516	1.2570475	2.219689	8773.307	8836.656
mod	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + aou + nitrate	5553.3021	5587.1515	0.9150165	1.743716	12340.952	12410.322
obs	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	5934.6854	5968.1369	7.0210613	11.700040	12198.544	12266.725
mod	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + aou + nitrate + silicate	5477.1931	5516.6840	0.8980205	1.707742	12139.328	12220.260
obs	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + aou + phosphate_star	5935.4248	5964.0975	7.0320102	11.711142	12197.352	12255.793
obs	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + aou + silicate + phosphate	5944.6333	5978.0848	7.0609034	11.751500	12213.725	12281.905
obs	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	5918.8677	5952.3192	6.9581722	11.620309	12175.118	12243.299
obs	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	5937.8076	5971.2591	7.0335415	11.712519	12201.666	12269.846
mod	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + nitrate + phosphate_star	5679.6032	5713.4526	0.9431819	1.803153	12671.125	12740.495
obs	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + nitrate + phosphate_star	5942.9546	5971.6273	7.0621940	11.741220	12204.835	12263.275
mod	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	5577.4744	5616.9654	0.9198994	1.778969	12565.225	12646.157
obs	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	5940.3781	5973.8296	7.0438335	11.707915	12197.508	12265.689
obs	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + phosphate + phosphate_star	5937.0787	5965.7514	7.0386290	11.717763	12199.008	12257.448
obs	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	5924.4102	5957.8617	6.9801439	11.643104	12181.033	12249.214
obs	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + silicate + phosphate_star	5950.2283	5978.9010	7.0914745	11.767175	12210.608	12269.049
mod	Atlantic	(28.15,28.2]	cstar_tref ~ temp + aou + nitrate	5847.5836	5875.7914	0.9824614	1.848043	12872.897	12930.705
mod	Atlantic	(28.15,28.2]	cstar_tref ~ temp + aou + nitrate + phosphate	5437.9687	5471.8180	0.8900322	1.713960	12193.836	12263.206
mod	Atlantic	(28.15,28.2]	cstar_tref ~ temp + aou + nitrate + phosphate_star	5462.6908	5496.5402	0.8953296	1.722366	12239.285	12308.655
mod	Atlantic	(28.15,28.2]	cstar_tref ~ temp + aou + nitrate + silicate	5502.6732	5536.5225	0.9039637	1.714010	12165.015	12234.385
mod	Atlantic	(28.15,28.2]	cstar_tref ~ temp + nitrate + phosphate + phosphate_star	5493.7449	5527.5943	0.9020284	1.733931	12302.625	12371.994
mod	Atlantic	(28.15,28.2]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	5579.4173	5613.2667	0.9207704	1.791943	12642.171	12711.540
mod	Atlantic	(28.2, Inf]	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	4103.4150	4136.2602	0.7726539	1.808555	12881.395	12950.306
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate	3134.3391	3161.7101	0.5872236	1.380743	10303.560	10360.986
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate + phosphate	3092.2922	3125.1374	0.5799295	1.358237	10146.838	10215.750
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate + phosphate_star	3093.2496	3126.0948	0.5800870	1.370015	10237.130	10306.041
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate + silicate	3136.1876	3169.0328	0.5871983	1.380711	10307.357	10376.268
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	4308.6056	4340.9926	4.1585695	8.459956	10187.717	10254.590
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	2986.2779	3024.5973	0.5624237	1.340188	10038.608	10119.005
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + silicate + phosphate	4306.8826	4334.6429	4.1593325	8.460768	10184.017	10241.337
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + aou + nitrate	3114.0749	3146.9202	0.5835252	1.371975	10246.662	10315.573
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	4322.9432	4355.3302	4.1982434	8.481597	10193.492	10260.365
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + aou + nitrate + silicate	3098.1781	3136.4976	0.5805693	1.367691	10222.603	10303.000
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + aou + silicate + phosphate	4297.8571	4330.2441	4.1290730	8.368633	10147.462	10214.335
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	4264.3917	4296.7787	4.0385689	8.328732	10138.178	10205.051
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	4139.6258	4177.9452	0.7801913	1.708438	12258.156	12338.552
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	4274.6850	4307.0720	4.0661928	8.336916	10139.216	10206.089
mod	Atlantic	(28.2, Inf]	cstar_tref ~ temp + aou + nitrate + silicate	4010.0516	4042.8968	0.7524523	1.732299	12452.690	12521.601
obs	Atlantic	(28.2, Inf]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	4296.7377	4329.1248	4.1260132	8.435595	10179.728	10246.601
obs	Atlantic	(28.2, Inf]	cstar_tref ~ temp + aou + silicate + phosphate_star	4297.8466	4325.6069	4.1345168	8.444125	10178.849	10236.169
obs	Atlantic	(28.2, Inf]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	4312.3827	4344.7697	4.1689846	8.461977	10187.513	10254.386
obs	Atlantic	(28.2, Inf]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	4311.7622	4339.5225	4.1727950	8.466701	10185.326	10242.646
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	15500.8535	15540.9215	7.4210982	14.136827	34870.282	34952.186
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	19561.5712	19603.0776	8.1608051	20.553076	45734.304	45818.571
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	15630.3059	15670.3740	7.6365173	13.767737	34469.409	34551.313
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	19344.0714	19385.5778	7.8475076	20.184911	45487.313	45571.580
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	15704.9255	15744.9935	7.7635197	13.959051	34604.799	34686.703
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	19314.3903	19355.8967	7.8056967	20.122900	45446.741	45531.008
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + phosphate	15683.0747	15717.4187	7.7295289	13.957028	34610.922	34681.126
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + phosphate_star	15764.4267	15798.7708	7.8697811	14.166577	34756.723	34826.927
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + phosphate_star	19580.2914	19615.8683	8.1912962	20.782205	45856.709	45928.937
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + silicate + phosphate	15490.1136	15530.1816	7.4035015	13.386426	34186.621	34268.525
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + silicate + phosphate	19322.1424	19363.6488	7.8165954	19.989438	45376.141	45460.407
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	15556.3325	15596.4005	7.5126653	13.534342	34290.440	34372.344
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	19301.4399	19342.9463	7.7875237	19.936768	45342.542	45426.808
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	15722.9324	15763.0004	7.7944824	14.013919	34645.234	34727.139
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	19286.4643	19327.9707	7.7665616	20.068562	45410.607	45494.874
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	15722.4798	15762.5478	7.7937027	14.146801	34768.619	34850.524
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	19508.5154	19550.0218	8.0832461	20.387662	45633.963	45718.230
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + phosphate + phosphate_star	19575.0120	19610.5889	8.1835164	20.772488	45850.407	45922.635
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	15579.6935	15619.7615	7.5515595	13.593255	34333.131	34415.035
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	19296.1659	19337.6723	7.7801349	19.923423	45334.009	45418.276
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	14614.6394	14655.3233	4.6490202	8.760372	29947.829	30029.830
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	14604.8074	14645.4912	4.6397765	8.624475	29768.780	29850.781
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	17357.7438	17398.6453	7.2751533	15.217676	36056.483	36138.629
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	14834.6690	14875.3528	4.8607702	9.071072	30296.474	30378.475
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	17351.4288	17392.3303	7.2661435	15.211475	36052.060	36134.206
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + phosphate_star	17514.9714	17550.0298	7.5060584	15.570070	36292.955	36363.366
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate	14245.5615	14286.2453	4.3143400	8.384839	29524.746	29606.747
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate	17504.5905	17545.4919	7.4878443	15.438167	36208.583	36290.729
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	14446.5502	14487.2340	4.4934929	8.754424	29972.999	30055.000
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	17501.0806	17541.9821	7.4826888	15.446153	36213.912	36296.058
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	14891.4878	14932.1716	4.9170004	9.177069	30416.842	30498.843
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	17351.8964	17392.7978	7.2668101	15.214432	36054.070	36136.216
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + nitrate + phosphate_star	17478.3528	17513.4112	7.4523151	15.532783	36267.247	36337.658
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	17480.2845	17521.1860	7.4522153	15.480049	36236.203	36318.349
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	14490.4003	14531.0841	4.5335571	8.833768	30066.475	30148.476
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	17501.1146	17542.0160	7.4827387	15.447668	36214.931	36297.077
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + silicate + phosphate_star	17501.4414	17536.4997	7.4861560	15.524756	36262.533	36332.944
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	14612.8029	14653.4867	4.6472922	8.818998	30024.804	30106.805
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	14719.7280	14760.4118	4.7489782	9.080612	30335.177	30417.178
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	14766.6160	14807.2998	4.7942677	9.186825	30457.598	30539.599
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	12774.1707	12808.3179	4.4643791	8.205061	25283.345	25351.889
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	12791.8422	12831.6806	4.4803885	8.038370	25074.477	25154.445
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	16374.9584	16415.8571	6.0066736	12.604545	33684.306	33766.290
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	12938.6343	12978.4727	4.6331605	8.393808	25474.104	25554.072
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	16396.2643	16437.1630	6.0318494	12.643883	33716.830	33798.814
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + aou + phosphate	16604.2016	16639.2576	6.2856316	13.030963	34027.194	34097.467
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	12586.6822	12626.5206	4.2752739	7.846465	24886.229	24966.197
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	16588.8498	16629.7485	6.2642572	13.006119	34011.151	34093.135
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	12700.1509	12739.9893	4.3875285	8.135367	25220.049	25300.017
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	16595.6356	16636.5343	6.2726074	13.040533	34038.124	34120.108
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	12946.9651	12986.8035	4.6419852	8.427066	25511.991	25591.959
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	16400.3789	16441.2776	6.0367235	12.653460	33724.665	33806.649
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + nitrate + phosphate_star	16469.5527	16504.6087	6.1216612	12.860083	33887.183	33957.456
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	12763.6103	12803.4487	4.4515893	8.332829	25443.136	25523.104
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	16388.0907	16428.9894	6.0221788	12.749619	33799.188	33881.173
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	12707.9514	12747.7898	4.3953530	8.157537	25245.285	25325.253
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	16594.0206	16634.9193	6.2706191	13.038471	34036.453	34118.437
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + silicate + phosphate_star	12769.8187	12803.9659	4.4599435	8.346582	25453.688	25522.232
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	12873.0927	12912.9311	4.5643156	8.354540	25444.317	25524.285
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	16521.3689	16562.2676	6.1818210	13.022218	34019.584	34101.569
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	14042.9926	14084.2005	3.3740335	6.461367	28173.383	28256.082
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	19938.7390	19981.2963	5.2978304	12.596696	43635.822	43721.455
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	14257.6320	14298.8398	3.5128383	6.847269	28814.877	28897.576
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	20086.6711	20129.2285	5.4206258	12.771251	43832.879	43918.512
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + aou + silicate + phosphate	14842.9625	14884.1704	3.9210765	7.717195	30119.138	30201.838
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + aou + silicate + phosphate	20023.7645	20066.3218	5.3680642	12.772298	43820.490	43906.123
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	14270.2766	14311.4844	3.5211913	6.886297	28878.290	28960.989
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	20113.8617	20156.4191	5.4435039	12.809311	43874.415	43960.047
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + nitrate + phosphate_star	14347.0550	14382.3760	3.5736812	7.201293	29369.507	29440.392
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	14330.2472	14371.4550	3.5610789	7.168901	29324.373	29407.072
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	19927.7817	19970.3390	5.2888463	12.637342	43671.976	43757.609
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + nitrate + phosphate	14535.7948	14571.1158	3.7026433	7.257496	29445.921	29516.806
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + nitrate + phosphate	20018.2087	20054.6864	5.3651084	12.693473	43741.332	43814.732
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	14442.4660	14483.6739	3.6369357	7.167422	29316.460	29399.160
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	19916.7578	19959.3152	5.2798232	12.553325	43589.640	43675.273
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	14664.4445	14705.6523	3.7917794	7.673428	30064.145	30146.844
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	20020.1001	20062.6575	5.3650182	12.680402	43732.899	43818.531
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	14680.9466	14722.1544	3.8035505	7.731935	30146.999	30229.698
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	20022.8224	20065.3797	5.3672810	12.690948	43743.488	43829.121
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	20021.7022	20058.1799	5.3680124	12.746875	43792.565	43865.965
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	22232.3746	22275.9569	4.7294881	11.433282	49636.389	49724.241
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + aou + silicate + phosphate	22232.1497	22269.5060	4.7306115	11.552211	49777.812	49853.114
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	15716.7213	15759.0382	2.9991534	5.707717	31097.222	31181.994
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	22455.6322	22499.2145	4.8728949	11.621651	49914.767	50002.620
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	15950.0660	15992.3829	3.1134673	5.945205	31613.499	31698.271
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + aou + silicate + phosphate	22219.1349	22262.7171	4.7211175	11.515508	49733.841	49821.694
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	22308.9780	22352.5602	4.7782113	11.681652	49954.982	50042.835
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	15967.3787	16009.6956	3.1221203	5.972207	31671.903	31756.674
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	15986.2825	16028.5993	3.1315960	6.170122	32098.122	32182.893
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	22359.5391	22403.1213	4.8106452	11.557767	49816.677	49904.530
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	22294.6121	22338.1944	4.7690359	11.657526	49922.738	50010.591
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ temp + aou + nitrate + phosphate	15836.1289	15872.4005	3.0580964	5.868854	31450.251	31522.912
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ temp + aou + nitrate + phosphate_star	16040.8611	16077.1327	3.1601288	6.091413	31922.102	31994.763
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	15709.8501	15752.1670	2.9958516	5.780162	31265.827	31350.599
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	22479.4511	22523.0334	4.8884491	11.533529	49810.925	49898.778
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	15959.1015	16001.4184	3.1179803	6.037427	31816.598	31901.370
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	22584.2886	22627.8708	4.9575023	11.671825	50001.243	50089.096
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ temp + nitrate + phosphate + phosphate_star	16054.7372	16091.0088	3.1671661	6.115345	31972.540	32045.201
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	15968.2679	16010.5848	3.1225654	6.057455	31859.329	31944.100
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	22575.5937	22619.1759	4.9517384	11.679005	50008.429	50096.282
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	23053.8569	23091.6047	4.3459138	10.347333	50333.113	50408.983
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	18593.6664	18637.2487	2.9065538	5.445736	38086.236	38174.122
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	22658.0791	22702.1181	4.1345426	10.324615	50202.169	50290.684
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	19044.0065	19087.5888	3.0870698	5.794077	39066.892	39154.779
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + aou + silicate + phosphate	22674.1437	22711.8915	4.1439151	10.334023	50216.282	50292.152
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	18898.9549	18942.5372	3.0277352	5.687869	38777.107	38864.993
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	19074.8372	19118.4195	3.0998305	5.659822	38635.035	38722.922
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate	22573.2208	22617.2599	4.0907985	10.161843	49952.429	50040.944
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	22677.9195	22721.9586	4.1448376	10.340372	50229.500	50318.015
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	18694.3984	18737.9806	2.9459925	5.460895	38107.355	38195.242
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	22720.3845	22764.4236	4.1669584	10.328896	50225.789	50314.304
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	22657.7444	22701.7835	4.1343692	10.316175	50190.487	50279.002
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + aou + nitrate + phosphate	19214.7146	19252.0708	3.1592359	5.732491	38815.732	38891.063
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	18839.0056	18882.5879	3.0035467	5.497644	38183.130	38271.016
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	22941.7073	22985.7463	4.2841750	10.341595	50301.837	50390.352
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	19048.6237	19092.2059	3.0889775	5.661612	38649.643	38737.530
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	23057.6179	23101.6570	4.3468733	10.511498	50566.725	50655.241
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + aou + silicate + phosphate	22969.7654	23007.5132	4.3003467	10.382857	50362.996	50438.866
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	19044.6516	19088.2339	3.0873363	5.661929	38651.976	38739.862
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	19051.3756	19088.7318	3.0909421	5.667465	38662.908	38738.239
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	22873.6402	22911.5159	3.9995054	9.375502	51680.436	51756.970
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	22814.5229	22858.7113	3.9696250	9.335389	51605.626	51694.916
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	19130.2944	19174.3282	2.6617692	5.119154	40409.961	40499.011
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	19271.0444	19315.0782	2.7091815	5.274526	40945.146	41034.197
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + phosphate + phosphate_star	22831.8332	22869.7089	3.9790416	9.347215	51625.105	51701.639
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + silicate + phosphate	19586.4440	19624.1872	2.8192217	5.407278	41339.406	41415.735
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	19229.2057	19273.2395	2.6950004	5.250579	40868.316	40957.366
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + nitrate + silicate + phosphate_star	19741.6619	19779.4052	2.8746509	5.452112	41456.991	41533.320
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	19520.0862	19557.8295	2.7958523	5.435778	41455.095	41531.424
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	19831.0827	19875.1165	2.9063478	5.387312	41198.355	41287.405
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	22179.0310	22223.2194	3.6718557	9.054915	51000.017	51089.307
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	19527.7523	19571.7861	2.7978402	5.431015	41441.271	41530.321
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	22226.7262	22270.9146	3.6934070	9.182432	51228.732	51318.022
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	19290.0246	19334.0584	2.7156394	5.212699	40716.621	40805.672
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	22266.5374	22310.7258	3.7114927	9.201992	51271.036	51360.327
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	19504.4053	19548.4391	2.7896583	5.417986	41401.020	41490.071
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	22221.0755	22265.2639	3.6908471	9.182288	51227.167	51316.457
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + silicate + phosphate_star	22350.9282	22388.8039	3.7510442	9.243671	51357.025	51433.559
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	22571.6687	22615.8571	3.8530829	9.235644	51391.795	51481.085
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ temp + aou + silicate + phosphate	22582.1228	22619.9986	3.8589753	9.242323	51401.607	51478.141
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + aou + nitrate + silicate	10891.1987	10925.9947	2.2508761	4.386211	21468.730	21538.290
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	10793.3842	10833.9796	2.2052863	4.321482	21329.229	21410.383
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	10831.9562	10872.5516	2.2227934	4.355969	21406.577	21487.731
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	10765.5101	10806.1055	2.1927207	4.221420	21096.479	21177.633
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + aou + nitrate + silicate	10892.7254	10933.3208	2.2506577	4.354587	21400.366	21481.519
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + aou + phosphate	10821.7182	10856.5142	2.2190426	4.399528	21500.772	21570.332
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + aou + silicate + phosphate	16842.1395	16884.5500	3.4656473	8.312054	34545.472	34629.812
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	10822.7089	10863.3043	2.2185835	4.414777	21538.590	21619.744
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	16849.5857	16891.9963	3.4697316	8.321709	34559.699	34644.039
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	10772.2282	10812.8236	2.1957426	4.237940	21135.374	21216.527
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	16846.4060	16888.8165	3.4679869	8.289604	34519.473	34603.813
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	16849.8379	16892.2484	3.4698700	8.322755	34561.055	34645.395
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	16839.0619	16881.4724	3.4639606	8.275522	34499.806	34584.146
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	11101.0758	11141.6712	2.3488709	4.402267	21490.756	21571.909
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	16849.2972	16891.7077	3.4695732	8.292066	34523.436	34607.776
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + aou + silicate + phosphate	16847.8086	16884.1605	3.4698539	8.318306	34551.633	34623.924
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + aou + silicate + phosphate_star	16853.3887	16889.7405	3.4729179	8.325526	34562.269	34634.560
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	11107.3401	11147.9355	2.3518892	4.415699	21521.564	21602.717
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	16848.9991	16891.4096	3.4694096	8.291563	34522.722	34607.062
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	16847.8359	16884.1877	3.4698689	8.292916	34520.653	34592.944
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + aou + nitrate + silicate	11272.9232	11308.2103	2.0302765	3.977885	22037.972	22108.387
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	11244.5440	11285.7123	2.0186593	3.915118	21874.372	21956.524
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	11175.1782	11216.3464	1.9923819	3.851628	21702.828	21784.980
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	11292.7043	11333.8726	2.0371071	3.959596	21992.819	22074.971
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou + nitrate	11123.8124	11159.0995	1.9738896	3.881302	21781.334	21851.750
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	11120.5392	11161.7075	1.9719244	3.879335	21780.057	21862.208
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	11121.9385	11163.1068	1.9724457	3.879829	21781.381	21863.533
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou + nitrate + silicate	10994.1828	11035.3510	1.9254161	3.759069	21450.365	21532.516
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate + phosphate	20364.5069	20408.5687	3.0759185	7.908978	44639.284	44727.493
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	20380.0686	20424.1304	3.0819046	7.907950	44643.079	44731.287
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	11134.3837	11175.5519	1.9770880	3.890556	21810.247	21892.399
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate	11285.1519	11326.3202	2.0342031	3.973840	22031.054	22113.205
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	20367.8619	20411.9237	3.0772081	7.858905	44556.075	44644.284
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	20380.6160	20424.6778	3.0821154	7.907908	44643.202	44731.411
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + silicate + phosphate_star	20381.1359	20418.9032	3.0830859	7.918670	44658.143	44733.751
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	20349.7486	20393.8104	3.0702521	7.850400	44535.336	44623.545
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	20372.0455	20416.1073	3.0788169	7.856105	44552.786	44640.994
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	20371.0500	20415.1118	3.0784340	7.855702	44551.755	44639.964
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	20372.2545	20410.0218	3.0796668	7.863810	44562.612	44638.220
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	20394.9905	20432.7578	3.0884272	7.937931	44695.287	44770.895
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	26254.4666	26293.7518	3.0860638	7.607907	63145.770	63225.543
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ aou + phosphate + phosphate_star	26269.4095	26302.1471	3.0911404	7.635231	63220.533	63287.011
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	26226.3483	26265.6335	3.0776570	7.619057	63172.009	63251.782
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + aou + nitrate	11268.6645	11303.7954	2.1456262	4.164599	22893.907	22964.525
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + aou + nitrate + silicate	11097.2262	11138.2123	2.0746788	3.989170	22433.595	22515.983
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate	11219.8463	11254.9772	2.1254146	4.155850	22876.077	22946.696
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate + phosphate	11190.3387	11231.3248	2.1124711	4.099259	22729.614	22812.003
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate + phosphate	26208.1397	26253.9724	3.0716293	7.573104	63044.575	63137.643
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	11192.2290	11233.2151	2.1132455	4.102267	22737.657	22820.045
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	26327.2849	26373.1176	3.1073387	7.634824	63236.294	63329.362
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + nitrate + silicate	11233.2749	11268.4059	2.1309553	4.151480	22862.722	22933.340
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	11204.3638	11245.3499	2.1182230	4.107028	22749.195	22831.583
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	26299.6826	26345.5153	3.0990291	7.622686	63198.035	63291.104
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + silicate + phosphate	11237.7331	11272.8640	2.1327979	4.165867	22901.061	22971.680
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	11192.9420	11233.9281	2.1135376	4.102936	22739.404	22821.793
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	26323.2286	26369.0613	3.1061162	7.634129	63233.705	63326.774
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + silicate + phosphate_star	11202.4163	11237.5472	2.1182445	4.110531	22754.821	22825.440
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + silicate + phosphate_star	26328.4713	26367.7565	3.1082994	7.636587	63237.713	63317.487
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	26234.6687	26280.5014	3.0795447	7.622081	63185.482	63278.550
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + aou + silicate + phosphate	26327.1074	26366.3926	3.1078881	7.677772	63351.527	63431.300
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	18825.7163	18869.4376	2.8533547	7.160528	45155.791	45244.504
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	18912.0730	18955.7944	2.8858583	7.132866	45113.572	45202.285
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + aou + phosphate	18921.2513	18958.7268	2.8900926	7.138158	45123.028	45199.067
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + aou + silicate + phosphate	7813.5236	7852.5095	1.8099921	3.471885	15932.046	16010.590
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + aou + silicate + phosphate	18775.4016	18819.1230	2.8345860	7.067099	44945.654	45034.367
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	7817.2175	7856.2034	1.8117178	3.475093	15939.492	16018.035
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	18837.5526	18881.2740	2.8577880	7.122373	45076.803	45165.515
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	7833.4399	7872.4258	1.8193164	3.514008	16034.164	16112.707
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	18829.1998	18872.9212	2.8546588	7.128686	45088.664	45177.377
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	7816.0091	7854.9950	1.8111531	3.473909	15936.717	16015.260
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	18835.0710	18878.7924	2.8568580	7.125580	45083.183	45171.895
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + silicate + phosphate_star	7876.2882	7909.7047	1.8404895	3.599408	16231.468	16298.791
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + silicate + phosphate_star	18837.7070	18875.1825	2.8585957	7.140775	45112.587	45188.627
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	18817.8336	18861.5550	2.8504061	7.133008	45095.616	45184.329
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + aou + silicate + phosphate	7832.8064	7866.2228	1.8199579	3.487028	15962.412	16029.735
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + aou + silicate + phosphate	18861.3825	18898.8579	2.8674865	7.162619	45163.870	45239.909
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + aou + silicate + phosphate_star	7837.7506	7871.1670	1.8222809	3.491230	15972.094	16039.417
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	7860.8936	7894.3101	1.8331940	3.549224	16112.245	16179.569
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	7836.1882	7869.6047	1.8215465	3.489729	15968.599	16035.922
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + silicate + phosphate_star	7895.1554	7923.0025	1.8504249	3.616206	16264.713	16320.816
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	12897.9935	12940.8833	1.6228152	3.378954	28625.733	28712.618
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	16897.5686	16940.8543	2.5543178	6.834445	43819.348	43907.796
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + aou + silicate + phosphate	16896.5403	16933.6423	2.5546642	6.839928	43827.553	43903.366
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + nitrate + silicate + phosphate	12939.3038	12976.0665	1.6332303	3.390908	28671.987	28746.460
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	12928.4411	12971.3308	1.6301301	3.387795	28663.067	28749.952
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	16965.7616	17009.0473	2.5787480	6.870724	43913.434	44001.882
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	16967.9312	17005.0333	2.5802493	6.875872	43921.558	43997.370
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	17045.0841	17088.3698	2.6074595	6.882156	43954.973	44043.421
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + aou + silicate + phosphate	12886.7788	12929.6686	1.6201292	3.393548	28692.147	28779.032
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + aou + silicate + phosphate	16880.0633	16923.3491	2.5480839	6.828398	43802.254	43890.701
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	12893.7498	12936.6396	1.6217983	3.396647	28705.507	28792.392
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	16915.3839	16958.6696	2.5606777	6.832312	43818.562	43907.009
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate	12901.5429	12944.4327	1.6236663	3.368157	28576.524	28663.408
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	12963.2363	13006.1261	1.6385300	3.458801	28975.340	29062.225
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	16978.0843	17021.3700	2.5831875	6.858378	43889.056	43977.504
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	12892.2240	12935.1137	1.6214328	3.395752	28701.616	28788.501
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	16905.1423	16948.4280	2.5570196	6.828752	43808.535	43896.983
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + silicate + phosphate_star	13001.4807	13038.2434	1.6482992	3.499015	29143.083	29217.556
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + silicate + phosphate_star	17010.5386	17047.6406	2.5956409	6.871328	43920.598	43996.411
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	13423.3386	13466.2284	1.7537589	3.519404	29193.881	29280.766
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ aou + nitrate + silicate + phosphate	18017.2270	18057.7206	1.0091829	2.184083	43756.767	43839.289
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	13742.1661	13789.4087	0.7188570	1.891902	39457.981	39554.256
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	15575.4677	15622.7102	0.8313661	2.021984	41533.379	41629.654
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + nitrate + silicate + phosphate	15684.5636	15725.0572	0.8387241	2.035015	41717.875	41800.397
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	13894.5771	13941.8197	0.7275996	1.904773	39667.586	39763.861
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	5362.1783	5397.1486	2.8006766	6.830318	13656.519	13728.550
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	5334.1373	5369.1077	2.7649479	6.822913	13649.099	13721.130
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + temp + aou + silicate + phosphate	5346.7822	5381.7526	2.7810027	6.811930	13642.063	13714.094
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	5378.0529	5413.0232	2.8211077	6.861865	13680.504	13752.535
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate	15400.2674	15447.5100	0.8198933	2.015066	41420.349	41516.624
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	5363.9068	5398.8771	2.8028941	6.840460	13664.032	13736.063
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star	5341.1909	5376.1613	2.7738922	6.860183	13676.631	13748.662
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	16853.3625	16900.6051	0.9200482	2.094227	42584.898	42681.173
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	5358.6127	5393.5831	2.7961081	6.833058	13658.289	13730.320
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	17036.2035	17083.4460	0.9334880	2.111507	42820.910	42917.185
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ temp + aou + silicate + phosphate	5358.0707	5388.0453	2.7979753	6.847730	13665.070	13726.811
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star	5340.4924	5375.4628	2.7730051	6.817790	13645.877	13717.908
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ temp + nitrate + silicate + phosphate	17552.7239	17593.2176	0.9726792	2.170801	43610.931	43693.453
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	16901.9661	16949.2086	0.9236018	2.097155	42624.835	42721.111
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	5376.1649	5411.1352	2.8186700	6.849853	13671.632	13743.663
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	7862.9951	7905.4519	0.8306095	1.729903	18760.049	18846.814
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	8424.9465	8467.4034	0.9072893	1.779101	19064.707	19151.471
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	7862.3826	7904.8394	0.8305296	1.748259	18927.663	19014.428
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + nitrate + phosphate + phosphate_star	7966.5572	8002.9488	0.8445021	1.763399	19040.375	19114.745
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	7959.7049	8002.1618	0.8433282	1.716176	18609.310	18696.075
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + nitrate + silicate + phosphate_star	8056.3806	8092.7722	0.8565061	1.776991	19144.516	19218.885
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	7844.6615	7887.1183	0.8282201	1.727055	18737.487	18824.252
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	1902.7515	1930.4060	2.8302734	6.839196	5173.406	5231.602
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	7868.4746	7910.9315	0.8313250	1.733220	18789.474	18876.238
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	1894.4762	1922.1307	2.7999407	6.807892	5164.850	5223.046
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + aou + phosphate	1909.1532	1932.8570	2.8614056	6.881025	5180.899	5230.781
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	1901.0758	1928.7303	2.8241047	6.839586	5173.627	5231.822
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + aou + phosphate_star	1899.2861	1922.9900	2.8248784	6.844237	5170.957	5220.839
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + aou + silicate + phosphate	1898.4691	1926.1236	2.8145355	6.828715	5170.644	5228.840
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	1898.4659	1926.1204	2.8145240	6.819555	5167.994	5226.190
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	7855.6905	7898.1473	0.8296567	1.728400	18747.671	18834.435
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	1899.2969	1926.9514	2.8175709	6.821869	5168.613	5226.808
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	8002.3229	8044.7798	0.8489947	1.764848	19050.641	19137.405
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + phosphate + phosphate_star	1905.0994	1928.8032	2.8463419	6.861826	5175.651	5225.533
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	1897.8259	1925.4804	2.8121794	6.818065	5167.602	5225.797
obs	Indo-Pacific	(28.4, Inf]	cstar_tref ~ sal + aou + nitrate + phosphate_star	313.2671	326.6717	2.1471818	6.221914	1341.815	1374.377
obs	Indo-Pacific	(28.4, Inf]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	312.1680	327.8067	2.0994995	6.257453	1349.994	1387.984
obs	Indo-Pacific	(28.4, Inf]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	313.7884	329.4272	2.1242980	6.124106	1337.655	1375.644
obs	Indo-Pacific	(28.4, Inf]	cstar_tref ~ sal + aou + silicate + phosphate_star	322.2846	335.6893	2.2921743	6.322792	1346.914	1379.476
obs	Indo-Pacific	(28.4, Inf]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	313.5732	329.2120	2.1209877	6.183236	1343.016	1381.006
obs	Indo-Pacific	(28.4, Inf]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	312.9971	328.6359	2.1121519	6.144075	1339.743	1377.732
obs	Indo-Pacific	(28.4, Inf]	cstar_tref ~ sal + temp + aou + silicate + phosphate	317.6166	333.2554	2.1840520	6.258315	1348.123	1386.112
obs	Indo-Pacific	(28.4, Inf]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	316.3949	332.0337	2.1648022	6.193529	1342.855	1380.844
obs	Indo-Pacific	(28.4, Inf]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	313.3606	328.9993	2.1177224	6.204029	1344.929	1382.919
obs	Indo-Pacific	(28.4, Inf]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	317.4155	333.0542	2.1808707	6.282610	1350.341	1388.331

7.3 Target variable coefficients

A data frame to map the target variable is prepared.

# extend best model selection with coefficients

lm_best <- left_join(
  lm_best,
  GLODAP_tidy %>% select(data_source, basin, gamma_slab, era, model, term, estimate)
)

7.4 dcant coeffcients

A data frame of coefficient offsets is prepared to facilitate the direct mapping of dcant.

# subtract coefficients of adjacent era  
lm_best_dcant <- lm_best %>%
  select(-c(AIC, BIC, rmse, aic_sum, bic_sum, rmse_sum)) %>% 
  arrange(era) %>%
  group_by(data_source, basin, gamma_slab, model, term) %>%
  mutate(delta_coeff = estimate - lag(estimate)) %>%
  ungroup() %>%
  arrange(basin, gamma_slab, model, term) %>%
  drop_na() %>%
  select(-c(era,estimate))

# pivot to wide format
lm_best_dcant <- lm_best_dcant %>%
  pivot_wider(values_from = delta_coeff,
              names_from = term,
              names_prefix = "delta_coeff_",
              values_fill = 0)

# create table of target variable coefficients in wide format
lm_best_target <- lm_best %>% 
  select(-c(AIC, BIC, rmse, aic_sum, bic_sum, rmse_sum)) %>% 
  pivot_wider(names_from = "term",
              names_prefix = "coeff_",
              values_from = "estimate",
              values_fill = 0
              )

GLODAP_augmented_best <- left_join(
  lm_best_target %>% select(data_source, basin, gamma_slab, era, model),
  GLODAP_augmented
)

7.5 Write files

lm_best_target %>%
  write_csv(paste(path_version_data,
                  "lm_best_target.csv",
                  sep = ""))

lm_best_dcant %>%
  write_csv(paste(path_version_data,
                  "lm_best_dcant.csv",
                  sep = ""))

GLODAP_augmented_best %>%
  write_csv(paste(path_version_data,
                  "lm_residuals_best_models.csv",
                  sep = ""))

8 Model diagnotics

8.1 Selection criterion vs predictors

The selection criterion (rmse) was plotted against the number of predictors (limited to 2 - 5).

8.1.1 All models

GLODAP_glanced %>%
  group_split(era) %>%
  # head(1) %>%
  map(
    ~ ggplot(data = .x,
             aes(
               as.factor(n_predictors),
               !!sym(params_local$MLR_criterion),
               col = data_source
             )) +
      geom_hline(yintercept = c(0, 10)) +
      geom_boxplot() +
      facet_grid(gamma_slab ~ basin) +
      scale_color_brewer(palette = "Set1") +
      ylim(c(0, NA)) +
      labs(x = "Number of predictors",
           title = paste("Era:", unique(.x$era)))
  )

[[1]]

Past versions of eMLR_predictor_model_quality_all-1.png

Version	Author	Date
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
7e1f407	jens-daniel-mueller	2021-06-10
2cbe18c	jens-daniel-mueller	2021-06-10
1772903	jens-daniel-mueller	2021-06-07
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
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
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
d3c3fde	jens-daniel-mueller	2021-03-26
594d503	jens-daniel-mueller	2021-03-26
2eb881a	jens-daniel-mueller	2021-03-26
03b6009	jens-daniel-mueller	2021-03-23
de6256b	jens-daniel-mueller	2021-03-23
f155edd	jens-daniel-mueller	2021-03-23
33b385b	jens-daniel-mueller	2021-03-20
108d6b2	jens-daniel-mueller	2021-03-20
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
ac1a836	jens-daniel-mueller	2021-02-24
b03fbd3	jens-daniel-mueller	2021-02-24
3d3b4cc	jens-daniel-mueller	2021-02-23
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
7131186	jens-daniel-mueller	2020-12-17
5b48ef5	jens-daniel-mueller	2020-12-17
f3a708f	jens-daniel-mueller	2020-12-17

[[2]]

Past versions of eMLR_predictor_model_quality_all-2.png

Version	Author	Date
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
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
d3c3fde	jens-daniel-mueller	2021-03-26
594d503	jens-daniel-mueller	2021-03-26
2eb881a	jens-daniel-mueller	2021-03-26

8.1.2 Best models

left_join(lm_best_target %>% select(basin, gamma_slab, era, model),
          GLODAP_glanced) %>%
  group_split(era) %>%
  # head(1) %>%
  imap(
  ~ ggplot(data = .x,
         aes("",
             !!sym(params_local$MLR_criterion),
             col = data_source)) +
  geom_hline(yintercept = c(0, 10)) +
  geom_boxplot() +
  facet_grid(gamma_slab ~ basin) +
  scale_color_brewer(palette = "Set1") +
  ylim(c(0, NA)) +
  labs(x = "Number of predictors pooled",
           title = paste("Era:", unique(.x$era)))
  )

[[1]]

Past versions of eMLR_predictor_model_quality_best-1.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
c6b3da6	jens-daniel-mueller	2021-06-14
7e1f407	jens-daniel-mueller	2021-06-10
2cbe18c	jens-daniel-mueller	2021-06-10
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
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
d3c3fde	jens-daniel-mueller	2021-03-26
594d503	jens-daniel-mueller	2021-03-26
06bfaf7	jens-daniel-mueller	2021-03-26
2eb881a	jens-daniel-mueller	2021-03-26
03b6009	jens-daniel-mueller	2021-03-23
de6256b	jens-daniel-mueller	2021-03-23
f155edd	jens-daniel-mueller	2021-03-23
380d215	jens-daniel-mueller	2021-03-21
33b385b	jens-daniel-mueller	2021-03-20
108d6b2	jens-daniel-mueller	2021-03-20
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
ac1a836	jens-daniel-mueller	2021-02-24
b03fbd3	jens-daniel-mueller	2021-02-24
3d3b4cc	jens-daniel-mueller	2021-02-23
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
7131186	jens-daniel-mueller	2020-12-17
5b48ef5	jens-daniel-mueller	2020-12-17
f3a708f	jens-daniel-mueller	2020-12-17

[[2]]

Past versions of eMLR_predictor_model_quality_best-2.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
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
d3c3fde	jens-daniel-mueller	2021-03-26
594d503	jens-daniel-mueller	2021-03-26

8.2 RMSE correlation between eras

RMSE was plotted to compare the agreement for one model applied to two adjecent eras (ie check whether the same predictor combination performs equal in both eras).

8.2.1 All models

# find max rmse to scale axis
max_rmse <-
  max(c(lm_all$rmse,
        lm_all$rmse_sum - lm_all$rmse))

lm_all %>%
  ggplot(aes(rmse, rmse_sum - rmse, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  coord_equal(xlim = c(0, max_rmse),
              ylim = c(0, max_rmse)) +
  geom_abline(slope = 1,
              col = 'red') +
  facet_grid(data_source ~ basin)

Past versions of rmse_correlation_across_eras_all_models-1.png

Version	Author	Date
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
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
594d503	jens-daniel-mueller	2021-03-26
06bfaf7	jens-daniel-mueller	2021-03-26
2eb881a	jens-daniel-mueller	2021-03-26
03b6009	jens-daniel-mueller	2021-03-23
de6256b	jens-daniel-mueller	2021-03-23
f155edd	jens-daniel-mueller	2021-03-23
380d215	jens-daniel-mueller	2021-03-21
33b385b	jens-daniel-mueller	2021-03-20
108d6b2	jens-daniel-mueller	2021-03-20
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
ac1a836	jens-daniel-mueller	2021-02-24
b03fbd3	jens-daniel-mueller	2021-02-24
3d3b4cc	jens-daniel-mueller	2021-02-23
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
7131186	jens-daniel-mueller	2020-12-17
5b48ef5	jens-daniel-mueller	2020-12-17
e4ca289	jens-daniel-mueller	2020-12-16
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
ba112d3	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
7c25f7a	jens-daniel-mueller	2020-12-02
b03ddb8	jens-daniel-mueller	2020-12-02
91435ae	jens-daniel-mueller	2020-12-01
196be51	jens-daniel-mueller	2020-11-30

rm(max_rmse)

8.2.2 Best models

# find max rmse to scale axis
max_rmse <-
  max(c(lm_best$rmse,
        lm_best$rmse_sum - lm_best$rmse))

lm_best %>%
  ggplot(aes(rmse, rmse_sum - rmse, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  coord_equal(xlim = c(0, max_rmse),
              ylim = c(0, max_rmse)) +
  geom_abline(slope = 1,
              col = 'red') +
  facet_grid(data_source ~ basin)

Past versions of rmse_correlation_across_eras_best_models-1.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
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
594d503	jens-daniel-mueller	2021-03-26
06bfaf7	jens-daniel-mueller	2021-03-26
2eb881a	jens-daniel-mueller	2021-03-26
03b6009	jens-daniel-mueller	2021-03-23
de6256b	jens-daniel-mueller	2021-03-23
f155edd	jens-daniel-mueller	2021-03-23
380d215	jens-daniel-mueller	2021-03-21
33b385b	jens-daniel-mueller	2021-03-20
108d6b2	jens-daniel-mueller	2021-03-20
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
ac1a836	jens-daniel-mueller	2021-02-24
b03fbd3	jens-daniel-mueller	2021-02-24
3d3b4cc	jens-daniel-mueller	2021-02-23
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
7131186	jens-daniel-mueller	2020-12-17
a84ff3c	jens-daniel-mueller	2020-12-17

rm(max_rmse)

8.3 Predictor counts

The number of models where a particular predictor was included were counted for each basin, density slab and compared eras

# calculate cases of predictor used
lm_all_stats <- lm_best_dcant %>% 
  pivot_longer(starts_with("delta_coeff_"),
               names_to = "term",
               names_prefix = "delta_coeff_",
               values_to = "delta_coeff") %>% 
  filter(term != "(Intercept)",
         delta_coeff != 0) %>% 
  group_by(data_source, basin, gamma_slab) %>% 
  count(term) %>% 
  ungroup() %>% 
  pivot_wider(values_from = n, names_from = term)

# print table
lm_all_stats %>%
  gt(rowname_col = "gamma_slab",
     groupname_col = c("basin", "data_source")) %>% 
  summary_rows(
    groups = TRUE,
    fns = list(total = "sum")
  )

	aou	nitrate	phosphate	phosphate_star	sal	silicate	temp
Atlantic - mod
(-Inf,26]	5	5	4	4	10	9	10
(26,26.5]	7	10	7	7	3	4	9
(26.5,26.75]	7	4	7	7	3	10	9
(26.75,27]	7	4	7	7	3	10	9
(27,27.25]	8	5	8	6	6	10	6
(27.25,27.5]	7	7	4	6	6	9	10
(27.5,27.75]	6	1	6	8	10	6	8
(27.75,27.85]	6	2	6	8	4	9	8
(27.85,27.95]	6	2	6	7	5	9	8
(27.95,28.05]	7	NA	8	6	6	7	8
(28.05,28.1]	7	NA	7	6	6	6	10
(28.1,28.15]	6	1	6	7	7	7	10
(28.15,28.2]	6	10	2	5	4	4	10
(28.2, Inf]	9	10	2	4	8	5	4
total	94.00	—	80.00	88.00	81.00	105.00	119.00
Indo-Pacific - mod
(-Inf,26]	7	4	6	7	10	5	9
(26,26.5]	7	7	7	7	6	7	9
(26.5,26.75]	6	5	6	7	8	7	9
(26.75,27]	6	9	6	6	6	5	10
(27,27.25]	6	10	6	7	4	4	10
(27.25,27.5]	6	10	6	6	5	8	7
(27.5,27.75]	5	6	6	7	10	9	4
(27.75,27.85]	8	8	4	6	8	7	7
(27.85,27.95]	7	10	5	4	10	6	6
(27.95,28.05]	5	4	3	4	10	9	10
(28.05,28.1]	4	2	4	8	5	10	10
(28.1,28.2]	4	6	7	5	9	10	7
(28.2,28.3]	5	10	8	4	5	10	5
(28.3,28.4]	5	10	6	8	10	5	4
total	81.00	101.00	80.00	86.00	106.00	102.00	107.00
Atlantic - obs
(-Inf,26]	7	6	6	6	10	4	9
(26,26.5]	10	4	10	6	5	5	4
(26.5,26.75]	7	4	7	6	9	4	10
(26.75,27]	4	8	4	10	6	5	10
(27,27.25]	7	9	6	7	4	4	10
(27.25,27.5]	6	6	6	8	7	5	10
(27.5,27.75]	7	4	8	7	3	10	8
(27.75,27.85]	6	5	8	6	6	10	7
(27.85,27.95]	5	4	7	6	10	10	5
(27.95,28.05]	5	4	7	6	10	10	5
(28.05,28.1]	5	6	5	9	10	8	5
(28.1,28.15]	4	6	3	10	6	8	10
(28.15,28.2]	4	4	4	9	10	5	10
(28.2, Inf]	7	4	7	6	6	9	8
total	84.00	74.00	88.00	102.00	102.00	97.00	111.00
Indo-Pacific - obs
(-Inf,26]	6	4	6	8	10	5	9
(26,26.5]	5	5	4	8	10	5	10
(26.5,26.75]	6	6	6	6	9	5	10
(26.75,27]	6	9	6	6	5	6	10
(27,27.25]	7	6	7	5	7	9	8
(27.25,27.5]	8	5	7	5	6	9	7
(27.5,27.75]	7	4	7	7	7	7	7
(27.75,27.85]	6	5	5	7	4	10	10
(27.85,27.95]	4	5	5	8	5	10	10
(27.95,28.05]	7	3	7	7	5	8	7
(28.05,28.1]	7	4	7	4	8	8	9
(28.1,28.2]	5	4	7	6	10	9	6
(28.2,28.3]	8	4	9	6	5	7	10
(28.3,28.4]	7	3	7	7	10	3	10
(28.4, Inf]	8	6	5	7	10	6	6
total	97.00	73.00	95.00	97.00	111.00	107.00	129.00

8.4 RMSE alternatives

8.4.1 AIC

AIC is an alternative criterion to RMSE to judge model quality, but not (yet) taken into account.

lm_all %>% 
  ggplot(aes(rmse, AIC, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  facet_grid(data_source~basin)

Past versions of aic_vs_rmse_all_models-1.png

Version	Author	Date
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
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
594d503	jens-daniel-mueller	2021-03-26
06bfaf7	jens-daniel-mueller	2021-03-26
2eb881a	jens-daniel-mueller	2021-03-26
03b6009	jens-daniel-mueller	2021-03-23
de6256b	jens-daniel-mueller	2021-03-23
f155edd	jens-daniel-mueller	2021-03-23
380d215	jens-daniel-mueller	2021-03-21
33b385b	jens-daniel-mueller	2021-03-20
108d6b2	jens-daniel-mueller	2021-03-20
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
ac1a836	jens-daniel-mueller	2021-02-24
b03fbd3	jens-daniel-mueller	2021-02-24
3d3b4cc	jens-daniel-mueller	2021-02-23
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
7131186	jens-daniel-mueller	2020-12-17
5b48ef5	jens-daniel-mueller	2020-12-17

lm_best %>% 
  ggplot(aes(rmse, AIC, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  facet_grid(data_source~basin)

Past versions of aic_vs_rmse_selected_models-1.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
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
594d503	jens-daniel-mueller	2021-03-26
06bfaf7	jens-daniel-mueller	2021-03-26
2eb881a	jens-daniel-mueller	2021-03-26
03b6009	jens-daniel-mueller	2021-03-23
de6256b	jens-daniel-mueller	2021-03-23
f155edd	jens-daniel-mueller	2021-03-23
380d215	jens-daniel-mueller	2021-03-21
33b385b	jens-daniel-mueller	2021-03-20
108d6b2	jens-daniel-mueller	2021-03-20
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
ac1a836	jens-daniel-mueller	2021-02-24
b03fbd3	jens-daniel-mueller	2021-02-24
3d3b4cc	jens-daniel-mueller	2021-02-23
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
7131186	jens-daniel-mueller	2020-12-17
5b48ef5	jens-daniel-mueller	2020-12-17

8.4.2 AIC vs BIC

BIC is an alternative criterion to RMSE to judge model quality, but not (yet) taken into account.

lm_all %>% 
  ggplot(aes(AIC, BIC, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  facet_grid(data_source~basin)

Past versions of BIC_vs_AIC_all_models-1.png

Version	Author	Date
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
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
594d503	jens-daniel-mueller	2021-03-26
06bfaf7	jens-daniel-mueller	2021-03-26
2eb881a	jens-daniel-mueller	2021-03-26
03b6009	jens-daniel-mueller	2021-03-23
de6256b	jens-daniel-mueller	2021-03-23
f155edd	jens-daniel-mueller	2021-03-23
380d215	jens-daniel-mueller	2021-03-21
33b385b	jens-daniel-mueller	2021-03-20
108d6b2	jens-daniel-mueller	2021-03-20

8.5 RMSE vs VIF

GLODAP_glanced %>% 
  ggplot(aes(rmse, log10(vif_max),col=data_source)) +
  geom_hline(yintercept = 1) +
  geom_point() +
  scale_color_brewer(palette = "Set1") +
  facet_grid(gamma_slab~basin)

Past versions of rmse_vs_VIF_all_models-1.png

Version	Author	Date
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
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
594d503	jens-daniel-mueller	2021-03-26
2eb881a	jens-daniel-mueller	2021-03-26
03b6009	jens-daniel-mueller	2021-03-23
de6256b	jens-daniel-mueller	2021-03-23
33b385b	jens-daniel-mueller	2021-03-20
108d6b2	jens-daniel-mueller	2021-03-20

8.6 Residual patterns

8.6.1 Fitted vs true

Plotted are fitted vs actual target variable values, here: rparams_local$MLR_target`

i_binwidth <- 2

p_prop_prop(df = GLODAP_augmented,
            var1 = params_local$MLR_target,
            var2 = ".fitted") +
  facet_grid(basin_AIP ~ data_source) +
  labs(title = "All models")

Past versions of fitted_vs_actual_all_property_property-1.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13

GLODAP_augmented %>%
  ggplot(aes(!!sym(params_local$MLR_target) - .fitted)) +
  geom_histogram(binwidth = i_binwidth) +
  labs(title = "All models") +
  scale_x_continuous(limits = c(-50,50)) +
  facet_wrap(~data_source)

Past versions of fitted_vs_actual_all_histogram-1.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13

p_prop_prop(df = GLODAP_augmented_best,
            var1 = params_local$MLR_target,
            var2 = ".fitted") +
  labs(title = "Selected models") +
  facet_grid(basin_AIP ~ data_source)

Past versions of fitted_vs_actual_best_property_property-1.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13

GLODAP_augmented_best %>%
  ggplot(aes(!!sym(params_local$MLR_target) - .fitted)) +
  geom_histogram(binwidth = i_binwidth) +
  labs(title = "Selected models") +
    scale_x_continuous(limits = c(-50,50)) +
  facet_wrap(~data_source)

Past versions of fitted_vs_actual_best_histogram-1.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13

8.6.2 Pooled

In the following, we present residual patterns vs latitude across all domains.

i_ylim <- c(-30,30)

GLODAP_augmented_best %>%
  ggplot(aes(lat, .resid)) +
  geom_bin2d(binwidth = i_binwidth) +
  geom_hline(yintercept = 0, col = "black") +
  scale_fill_viridis_c(trans = "log10") +
  labs(
    title = paste(
      "Target variable:",
      params_local$MLR_target,
      "| Selected models",
      "| All domains"
    )
  ) +
  facet_wrap(~data_source)

Past versions of residual_patterns_full_range-1.png

Version	Author	Date
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
147cb07	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
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
594d503	jens-daniel-mueller	2021-03-26
06bfaf7	jens-daniel-mueller	2021-03-26
2eb881a	jens-daniel-mueller	2021-03-26
03b6009	jens-daniel-mueller	2021-03-23
de6256b	jens-daniel-mueller	2021-03-23
f155edd	jens-daniel-mueller	2021-03-23
380d215	jens-daniel-mueller	2021-03-21
33b385b	jens-daniel-mueller	2021-03-20
108d6b2	jens-daniel-mueller	2021-03-20
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
27ae473	jens-daniel-mueller	2021-02-24
4bc00ea	jens-daniel-mueller	2021-02-24

Due to the few large residuals, we limit the y axis range for the plots below.

GLODAP_augmented_best %>%
  ggplot(aes(lat, .resid)) +
  geom_bin2d(binwidth = i_binwidth) +
  geom_hline(yintercept = 0, col = "black") +
  scale_fill_viridis_c(trans = "log10") +
  coord_cartesian(ylim = i_ylim) +
  labs(
    title = paste(
      "Target variable:",
      params_local$MLR_target,
      "| Selected models",
      "| All domains"
    )
  ) +
  facet_wrap(~data_source)

Past versions of residual_patterns_y_range_limit-1.png

Version	Author	Date
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
147cb07	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
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
d3c3fde	jens-daniel-mueller	2021-03-26
594d503	jens-daniel-mueller	2021-03-26
06bfaf7	jens-daniel-mueller	2021-03-26
2eb881a	jens-daniel-mueller	2021-03-26
03b6009	jens-daniel-mueller	2021-03-23
de6256b	jens-daniel-mueller	2021-03-23
f155edd	jens-daniel-mueller	2021-03-23
380d215	jens-daniel-mueller	2021-03-21
33b385b	jens-daniel-mueller	2021-03-20
108d6b2	jens-daniel-mueller	2021-03-20
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
27ae473	jens-daniel-mueller	2021-02-24
4bc00ea	jens-daniel-mueller	2021-02-24

8.6.3 By model domain

In the following, we present residual patterns vs latitude for separate model domains, ie basins, density slabs and eras.

GLODAP_augmented_best %>% 
  group_split(basin, data_source) %>% 
  # head(1) %>% 
  map(~ ggplot(data = .x, aes(lat, .resid)) +
      geom_bin2d(binwidth = i_binwidth) +
      geom_hline(yintercept = 0, col = "black") +
      scale_fill_viridis_c(trans = "log10") +
      facet_grid(gamma_slab ~ era) +
      coord_cartesian(ylim = i_ylim) +
      labs(
        title = paste(
          "Target:",
          params_local$MLR_target,
          "| best models | basin:",
          unique(.x$basin),
          "| data_source:",
          unique(.x$data_source)
        )
      ))

[[1]]

Past versions of residual_patterns_by_model_domain-1.png

Version	Author	Date
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
147cb07	jens-daniel-mueller	2021-06-29
d054c39	jens-daniel-mueller	2021-06-24
4fcd066	jens-daniel-mueller	2021-06-24
c6b3da6	jens-daniel-mueller	2021-06-14
7e1f407	jens-daniel-mueller	2021-06-10
2cbe18c	jens-daniel-mueller	2021-06-10
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
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
d3c3fde	jens-daniel-mueller	2021-03-26
594d503	jens-daniel-mueller	2021-03-26
06bfaf7	jens-daniel-mueller	2021-03-26
2eb881a	jens-daniel-mueller	2021-03-26
03b6009	jens-daniel-mueller	2021-03-23
de6256b	jens-daniel-mueller	2021-03-23
f155edd	jens-daniel-mueller	2021-03-23
380d215	jens-daniel-mueller	2021-03-21
33b385b	jens-daniel-mueller	2021-03-20
108d6b2	jens-daniel-mueller	2021-03-20
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
27ae473	jens-daniel-mueller	2021-02-24
a1ba577	jens-daniel-mueller	2021-02-24
071743d	jens-daniel-mueller	2021-02-24

[[2]]

Past versions of residual_patterns_by_model_domain-2.png

Version	Author	Date
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
147cb07	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
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
d3c3fde	jens-daniel-mueller	2021-03-26
594d503	jens-daniel-mueller	2021-03-26
06bfaf7	jens-daniel-mueller	2021-03-26
2eb881a	jens-daniel-mueller	2021-03-26
03b6009	jens-daniel-mueller	2021-03-23
de6256b	jens-daniel-mueller	2021-03-23
f155edd	jens-daniel-mueller	2021-03-23
380d215	jens-daniel-mueller	2021-03-21
33b385b	jens-daniel-mueller	2021-03-20
108d6b2	jens-daniel-mueller	2021-03-20
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
27ae473	jens-daniel-mueller	2021-02-24
a1ba577	jens-daniel-mueller	2021-02-24

[[3]]

Past versions of residual_patterns_by_model_domain-3.png

Version	Author	Date
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
147cb07	jens-daniel-mueller	2021-06-29
d054c39	jens-daniel-mueller	2021-06-24
4fcd066	jens-daniel-mueller	2021-06-24
c6b3da6	jens-daniel-mueller	2021-06-14
7e1f407	jens-daniel-mueller	2021-06-10
2cbe18c	jens-daniel-mueller	2021-06-10
1772903	jens-daniel-mueller	2021-06-07
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
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
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
d3c3fde	jens-daniel-mueller	2021-03-26
594d503	jens-daniel-mueller	2021-03-26

[[4]]

Past versions of residual_patterns_by_model_domain-4.png

Version	Author	Date
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
147cb07	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
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
d3c3fde	jens-daniel-mueller	2021-03-26
594d503	jens-daniel-mueller	2021-03-26

8.6.4 Latitudinal mean

GLODAP_augmented_best <- GLODAP_augmented_best %>%
  mutate(lat_grid = as.numeric(as.character(cut(
    lat,
    seq(-90, 90, 10),
    seq(-85, 85, 10)
  ))))

lat_residual <- GLODAP_augmented_best %>%
  group_by(data_source, basin, gamma_slab, era, lat_grid) %>% 
  summarise(.resid_mean = mean(.resid)) %>% 
  ungroup()

lat_residual %>%
  group_split(basin) %>% 
  map(
  ~ggplot(data = .x,
          aes(lat_grid, .resid_mean, col=era)) +
  geom_line() +
  geom_point() +
  geom_hline(yintercept = 0, col = "black") +
  facet_grid(gamma_slab ~ data_source) +
    labs(title = paste("Basin:", unique(.x$basin)))
  )

[[1]]

Past versions of latitudinal_residual-1.png

Version	Author	Date
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
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
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
594d503	jens-daniel-mueller	2021-03-26
06bfaf7	jens-daniel-mueller	2021-03-26
2eb881a	jens-daniel-mueller	2021-03-26
03b6009	jens-daniel-mueller	2021-03-23
de6256b	jens-daniel-mueller	2021-03-23
f155edd	jens-daniel-mueller	2021-03-23
380d215	jens-daniel-mueller	2021-03-21
33b385b	jens-daniel-mueller	2021-03-20
108d6b2	jens-daniel-mueller	2021-03-20
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
fec3558	jens-daniel-mueller	2021-02-24

[[2]]

Past versions of latitudinal_residual-2.png

Version	Author	Date
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
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
594d503	jens-daniel-mueller	2021-03-26

8.6.5 Latitudinal offset

# calculate residual offset for adjacent eras
lat_residual_offset <- lat_residual  %>%
  select(data_source, basin, gamma_slab, era, lat_grid, .resid_mean) %>% 
  arrange(era) %>% 
  group_by(data_source, basin, gamma_slab, lat_grid) %>% 
  mutate(.resid_mean_offset = .resid_mean - lag(.resid_mean)) %>% 
  ungroup() %>% 
  drop_na()

lat_residual_offset %>%
  ggplot(aes(lat_grid, .resid_mean_offset, col=data_source)) +
  geom_line() +
  geom_point() +
  geom_hline(yintercept = 0, col = "black") +
  facet_grid(gamma_slab ~ basin)

Past versions of latitudinal_residual_offset-1.png

Version	Author	Date
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
147cb07	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
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
594d503	jens-daniel-mueller	2021-03-26
06bfaf7	jens-daniel-mueller	2021-03-26
2eb881a	jens-daniel-mueller	2021-03-26
03b6009	jens-daniel-mueller	2021-03-23
de6256b	jens-daniel-mueller	2021-03-23
f155edd	jens-daniel-mueller	2021-03-23
380d215	jens-daniel-mueller	2021-03-21
33b385b	jens-daniel-mueller	2021-03-20
108d6b2	jens-daniel-mueller	2021-03-20
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
fec3558	jens-daniel-mueller	2021-02-24

8.6.6 Spatial mean

GLODAP_augmented_best <- GLODAP_augmented_best %>%
  mutate(lon_grid = as.numeric(as.character(cut(
    lon,
    seq(-360, 500, 10),
    seq(-355, 500, 10)
  ))))

spatial_residual <- GLODAP_augmented_best %>%
  group_by(data_source, gamma_slab, era, lat_grid, lon_grid) %>% 
  summarise(.resid_mean = mean(.resid)) %>% 
  ungroup()

spatial_obs <- GLODAP_augmented_best %>%
  distinct(gamma_slab, era, lat, lon)

spatial_residual %>%
  rename(lat = lat_grid,
         lon = lon_grid) %>%
  group_split(data_source, gamma_slab) %>%
  # head(1) %>%
  map(~ p_map_dcant_slab(
    df = .x,
    var = ".resid_mean",
    col = "divergent",
    title_text = "Residual distribution",
    subtitle_text = paste(
      "Data source:",
      unique(.x$data_source),
      "| density slab:",
      unique(.x$gamma_slab)
      ),
    legend_title = paste(params_local$MLR_target, "\nbias")
  ) +
    geom_raster(data = spatial_obs,
                aes(lon, lat)) +
    facet_grid(era ~ .)
  )

[[1]]

Past versions of spatial_residual-1.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[2]]

Past versions of spatial_residual-2.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[3]]

Past versions of spatial_residual-3.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[4]]

Past versions of spatial_residual-4.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[5]]

Past versions of spatial_residual-5.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[6]]

Past versions of spatial_residual-6.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[7]]

Past versions of spatial_residual-7.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[8]]

Past versions of spatial_residual-8.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[9]]

Past versions of spatial_residual-9.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[10]]

Past versions of spatial_residual-10.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[11]]

Past versions of spatial_residual-11.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[12]]

Past versions of spatial_residual-12.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
910d64a	jens-daniel-mueller	2021-07-02
1cbf907	jens-daniel-mueller	2021-07-02
955f5dc	jens-daniel-mueller	2021-06-29
147cb07	jens-daniel-mueller	2021-06-29

[[13]]

Past versions of spatial_residual-13.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[14]]

Past versions of spatial_residual-14.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[15]]

Past versions of spatial_residual-15.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
910d64a	jens-daniel-mueller	2021-07-02
1cbf907	jens-daniel-mueller	2021-07-02
955f5dc	jens-daniel-mueller	2021-06-29
147cb07	jens-daniel-mueller	2021-06-29

[[16]]

Past versions of spatial_residual-16.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[17]]

Past versions of spatial_residual-17.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[18]]

Past versions of spatial_residual-18.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[19]]

Past versions of spatial_residual-19.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[20]]

Past versions of spatial_residual-20.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[21]]

Past versions of spatial_residual-21.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[22]]

Past versions of spatial_residual-22.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[23]]

Past versions of spatial_residual-23.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[24]]

Past versions of spatial_residual-24.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[25]]

Past versions of spatial_residual-25.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[26]]

Past versions of spatial_residual-26.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[27]]

Past versions of spatial_residual-27.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[28]]

Past versions of spatial_residual-28.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[29]]

Past versions of spatial_residual-29.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[30]]

Past versions of spatial_residual-30.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[31]]

Past versions of spatial_residual-31.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[32]]

Past versions of spatial_residual-32.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
910d64a	jens-daniel-mueller	2021-07-02
1cbf907	jens-daniel-mueller	2021-07-02
955f5dc	jens-daniel-mueller	2021-06-29
147cb07	jens-daniel-mueller	2021-06-29

[[33]]

Past versions of spatial_residual-33.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[34]]

Past versions of spatial_residual-34.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
910d64a	jens-daniel-mueller	2021-07-02
1cbf907	jens-daniel-mueller	2021-07-02
955f5dc	jens-daniel-mueller	2021-06-29
147cb07	jens-daniel-mueller	2021-06-29

[[35]]

Past versions of spatial_residual-35.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
910d64a	jens-daniel-mueller	2021-07-02
1cbf907	jens-daniel-mueller	2021-07-02
955f5dc	jens-daniel-mueller	2021-06-29
147cb07	jens-daniel-mueller	2021-06-29

8.6.7 Spatial offset

# calculate residual offset for adjacent eras
spatial_residual_offset <- spatial_residual  %>%
  select(data_source, gamma_slab, era, lat_grid, lon_grid, .resid_mean) %>% 
  arrange(era) %>% 
  group_by(data_source, gamma_slab, lat_grid, lon_grid) %>% 
  mutate(.resid_mean_offset = .resid_mean - lag(.resid_mean)) %>% 
  ungroup() %>% 
  drop_na()

spatial_obs <- spatial_obs %>% 
  count(lon, lat, gamma_slab) %>% 
  filter(n == 2) %>% 
  select(-n)

spatial_residual_offset %>%
  rename(lat = lat_grid,
         lon = lon_grid) %>%
  group_split(data_source, gamma_slab) %>%
  # head(1) %>%
  map(~ p_map_dcant_slab(
    df = .x,
    var = ".resid_mean_offset",
    col = "divergent",
    title_text = "Residual offset distribution (era 2 - era 1)",
    subtitle_text = paste(
      "Data source:",
      unique(.x$data_source),
      "| density slab:",
      unique(.x$gamma_slab)
      ),
    legend_title = paste(params_local$MLR_target, "\nbias offset")
  ) +
    geom_raster(data = spatial_obs,
                aes(lon, lat))
  )

[[1]]

Past versions of Spatial_residual_offset-1.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[2]]

Past versions of Spatial_residual_offset-2.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[3]]

Past versions of Spatial_residual_offset-3.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[4]]

Past versions of Spatial_residual_offset-4.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[5]]

Past versions of Spatial_residual_offset-5.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[6]]

Past versions of Spatial_residual_offset-6.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[7]]

Past versions of Spatial_residual_offset-7.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[8]]

Past versions of Spatial_residual_offset-8.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[9]]

Past versions of Spatial_residual_offset-9.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[10]]

Past versions of Spatial_residual_offset-10.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[11]]

Past versions of Spatial_residual_offset-11.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[12]]

Past versions of Spatial_residual_offset-12.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
910d64a	jens-daniel-mueller	2021-07-02
1cbf907	jens-daniel-mueller	2021-07-02
955f5dc	jens-daniel-mueller	2021-06-29
147cb07	jens-daniel-mueller	2021-06-29

[[13]]

Past versions of Spatial_residual_offset-13.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[14]]

Past versions of Spatial_residual_offset-14.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
910d64a	jens-daniel-mueller	2021-07-02
1cbf907	jens-daniel-mueller	2021-07-02
955f5dc	jens-daniel-mueller	2021-06-29
147cb07	jens-daniel-mueller	2021-06-29

[[15]]

Past versions of Spatial_residual_offset-15.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[16]]

Past versions of Spatial_residual_offset-16.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[17]]

Past versions of Spatial_residual_offset-17.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
910d64a	jens-daniel-mueller	2021-07-02
1cbf907	jens-daniel-mueller	2021-07-02
955f5dc	jens-daniel-mueller	2021-06-29
147cb07	jens-daniel-mueller	2021-06-29

[[18]]

Past versions of Spatial_residual_offset-18.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[19]]

Past versions of Spatial_residual_offset-19.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[20]]

Past versions of Spatial_residual_offset-20.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[21]]

Past versions of Spatial_residual_offset-21.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[22]]

Past versions of Spatial_residual_offset-22.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[23]]

Past versions of Spatial_residual_offset-23.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[24]]

Past versions of Spatial_residual_offset-24.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[25]]

Past versions of Spatial_residual_offset-25.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[26]]

Past versions of Spatial_residual_offset-26.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[27]]

Past versions of Spatial_residual_offset-27.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
910d64a	jens-daniel-mueller	2021-07-02
1cbf907	jens-daniel-mueller	2021-07-02
955f5dc	jens-daniel-mueller	2021-06-29
147cb07	jens-daniel-mueller	2021-06-29

[[28]]

Past versions of Spatial_residual_offset-28.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[29]]

Past versions of Spatial_residual_offset-29.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[30]]

Past versions of Spatial_residual_offset-30.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
910d64a	jens-daniel-mueller	2021-07-02
1cbf907	jens-daniel-mueller	2021-07-02
955f5dc	jens-daniel-mueller	2021-06-29
147cb07	jens-daniel-mueller	2021-06-29

[[31]]

Past versions of Spatial_residual_offset-31.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[32]]

Past versions of Spatial_residual_offset-32.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
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
147cb07	jens-daniel-mueller	2021-06-29

[[33]]

Past versions of Spatial_residual_offset-33.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
910d64a	jens-daniel-mueller	2021-07-02
1cbf907	jens-daniel-mueller	2021-07-02
955f5dc	jens-daniel-mueller	2021-06-29
147cb07	jens-daniel-mueller	2021-06-29

[[34]]

Past versions of Spatial_residual_offset-34.png

Version	Author	Date
912d90e	jens-daniel-mueller	2021-07-23
2477316	jens-daniel-mueller	2021-07-23
c9ccc00	jens-daniel-mueller	2021-07-22
971ce87	jens-daniel-mueller	2021-07-13
910d64a	jens-daniel-mueller	2021-07-02
1cbf907	jens-daniel-mueller	2021-07-02
955f5dc	jens-daniel-mueller	2021-06-29
147cb07	jens-daniel-mueller	2021-06-29

[[35]]

Past versions of Spatial_residual_offset-35.png

Version	Author	Date
955f5dc	jens-daniel-mueller	2021-06-29
147cb07	jens-daniel-mueller	2021-06-29

9 VIF

GLODAP_glanced %>% 
  ggplot(aes(as.factor(n),
             log10(vif_max),
             col=data_source)) + 
  geom_hline(yintercept = log10(params_local$vif_max)) +
  scale_color_brewer(palette = "Set1") +
  geom_violin() +
  facet_grid(gamma_slab~basin)

Past versions of vif_distribution_n_predictors-1.png

Version	Author	Date
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
594d503	jens-daniel-mueller	2021-03-26
2eb881a	jens-daniel-mueller	2021-03-26
03b6009	jens-daniel-mueller	2021-03-23
de6256b	jens-daniel-mueller	2021-03-23
380d215	jens-daniel-mueller	2021-03-21
33b385b	jens-daniel-mueller	2021-03-20
108d6b2	jens-daniel-mueller	2021-03-20
330dcd0	jens-daniel-mueller	2021-03-20

GLODAP_glanced_vif %>% 
  ggplot(aes(Variables,
             log10(VIF),
             col = data_source)) + 
  geom_hline(yintercept = log10(params_local$vif_max)) +
  geom_violin() +
  scale_color_brewer(palette = "Set1") +
  facet_grid(gamma_slab~basin) +
  coord_flip()

Past versions of vif_distribution_type_predictor-1.png

Version	Author	Date
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
d3c3fde	jens-daniel-mueller	2021-03-26
594d503	jens-daniel-mueller	2021-03-26
2eb881a	jens-daniel-mueller	2021-03-26
03b6009	jens-daniel-mueller	2021-03-23
de6256b	jens-daniel-mueller	2021-03-23
380d215	jens-daniel-mueller	2021-03-21
33b385b	jens-daniel-mueller	2021-03-20
108d6b2	jens-daniel-mueller	2021-03-20
330dcd0	jens-daniel-mueller	2021-03-20

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] tictoc_1.0.1     gt_0.2.2         corrr_0.4.3      broom_0.7.5     
 [5] kableExtra_1.3.1 knitr_1.30       olsrr_0.5.3      GGally_2.0.0    
 [9] lubridate_1.7.9  ggforce_0.3.3    metR_0.9.0       scico_1.2.0     
[13] patchwork_1.1.1  collapse_1.5.0   forcats_0.5.0    stringr_1.4.0   
[17] dplyr_1.0.5      purrr_0.3.4      readr_1.4.0      tidyr_1.1.2     
[21] tibble_3.0.4     ggplot2_3.3.3    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         utf8_1.1.4               R6_2.5.0                
[10] nortest_1.0-4            DBI_1.1.0                colorspace_1.4-1        
[13] withr_2.3.0              gridExtra_2.3            tidyselect_1.1.0        
[16] curl_4.3                 compiler_4.0.3           git2r_0.27.1            
[19] cli_2.1.0                rvest_0.3.6              xml2_1.3.2              
[22] sass_0.2.0               labeling_0.4.2           scales_1.1.1            
[25] checkmate_2.0.0          goftest_1.2-2            digest_0.6.27           
[28] foreign_0.8-80           rmarkdown_2.5            rio_0.5.16              
[31] pkgconfig_2.0.3          htmltools_0.5.0          highr_0.8               
[34] dbplyr_1.4.4             rlang_0.4.10             readxl_1.3.1            
[37] rstudioapi_0.11          farver_2.0.3             generics_0.0.2          
[40] jsonlite_1.7.1           zip_2.1.1                car_3.0-10              
[43] magrittr_1.5             Matrix_1.2-18            Rcpp_1.0.5              
[46] munsell_0.5.0            fansi_0.4.1              abind_1.4-5             
[49] lifecycle_1.0.0          stringi_1.5.3            whisker_0.4             
[52] yaml_2.2.1               carData_3.0-4            MASS_7.3-53             
[55] plyr_1.8.6               grid_4.0.3               blob_1.2.1              
[58] parallel_4.0.3           promises_1.1.1           crayon_1.3.4            
[61] lattice_0.20-41          haven_2.3.1              hms_0.5.3               
[64] pillar_1.4.7             reprex_0.3.0             glue_1.4.2              
[67] evaluate_0.14            RcppArmadillo_0.10.1.2.0 data.table_1.13.2       
[70] modelr_0.1.8             vctrs_0.3.5              tweenr_1.0.2            
[73] httpuv_1.5.4             cellranger_1.1.0         gtable_0.3.0            
[76] polyclip_1.10-0          reshape_0.8.8            assertthat_0.2.1        
[79] openxlsx_4.2.3           xfun_0.18                RcppEigen_0.3.3.7.0     
[82] later_1.1.0.1            viridisLite_0.3.0        ellipsis_0.3.1          
[85] here_0.1