eMLR - model fitting

workflowr

• Summary
• Checks
• Past versions

Last updated: 2021-05-12

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: 5f542d0

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 5f542d0. 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

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	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 Required data

Required are:

• cleaned and prepared GLODAPv2_2020 file

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

2 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)

3 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

4 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)]) %>%
  group_by(gamma_slab, era, basin, data_source) %>% 
  nest()

# expand with model definitions

GLODAP_nested_lm <- expand_grid(
  GLODAP_nested,
  lm_all
)

# 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)
  )


# test parallel processing
# 
# library(parallel)
# library(multidplyr)
# 
# n_cores <- detectCores()
# system.time(
# cluster <- new_cluster(5)
# )
# 
# cluster
# 
# system.time(
#   GLODAP_nested_lm <-
#     expand_grid(GLODAP,
#                 lm_all) %>%
#     group_by(gamma_slab, era, basin, data_source, model, predictors, n) %>%
#     partition(cluster)
# )
# 
# GLODAP_nested_lm 
# 
# cluster %>%
#   cluster_library(c("tidyverse",
#                     "broom",
#                     "olsrr",
#                     "purrr"))
# 
# cluster_call(cluster, search())
# 
# system.time(
# 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)
#   )
# )
# 
# system.time(
# GLODAP_nested_lm_fit <- GLODAP_nested_lm_fit %>% 
#   collect()
# )
# 
# 
# 
# library(furrr)
# 
# plan(sequential)
# plan(multisession, workers = 5)
# 
# # This does run in parallel!
# system.time(
# GLODAP_nested_lm_fit <- GLODAP_nested_lm %>% 
#   mutate(
#     fit = future_map2(.x = data, .y = model, 
#                ~ lm(as.formula(.y), data = .x))
#     # tidied = future_map(fit, tidy),
#     # glanced = future_map(fit, glance),
#     # augmented = future_map(fit, augment),
#     # vif = future_map(fit, ols_vif_tol)
#   )
# )



# print(object.size(GLODAP_nested), units = "MB")
# print(object.size(GLODAP_nested_lm), units = "MB")
# print(object.size(GLODAP_nested_lm_fit), units = "MB")

5 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)

6 Prepare coeffcients

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

6.1 VIF threshold

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

• 200

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>

6.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_eras <- 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(eras = paste(lag(era), era, sep = " --> "),
         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_eras <- lm_all_eras %>%
    group_by(data_source, basin, gamma_slab, eras) %>%
    slice_min(order_by = aic_sum,
              with_ties = FALSE,
              n = params_local$MLR_number) %>%
    ungroup() %>%
    arrange(basin, gamma_slab, eras, model)
} 

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

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


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

data_source	basin	gamma_slab	model	AIC	BIC	rmse	eras	rmse_sum	aic_sum	bic_sum
mod	Atlantic	(-Inf,26]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	258.7583	276.0196	0.9877932	2000-2009 –> 2010-2019	4.670380	1035.073	1072.926
obs	Atlantic	(-Inf,26]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	895.4113	915.6438	6.6501395	2000-2009 –> 2010-2019	14.216337	1370.724	1406.389
mod	Atlantic	(-Inf,26]	cstar_tref ~ nitrate + silicate + phosphate + phosphate_star	259.3824	274.1779	1.0028034	2000-2009 –> 2010-2019	4.687030	1033.822	1066.267
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	898.5943	918.8267	6.7301934	2000-2009 –> 2010-2019	13.830070	1365.382	1401.048
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	222.2509	239.5123	0.8008399	2000-2009 –> 2010-2019	4.562658	1004.526	1042.379
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	231.9165	249.1778	0.8465846	2000-2009 –> 2010-2019	4.513834	1007.063	1044.916
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	232.4643	247.2598	0.8590722	2000-2009 –> 2010-2019	4.751332	1022.284	1054.730
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate	899.9583	917.3004	6.8158486	2000-2009 –> 2010-2019	13.980517	1365.964	1396.534
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	898.8621	919.0945	6.7369719	2000-2009 –> 2010-2019	13.897097	1366.782	1402.448
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	898.8896	919.1221	6.7376702	2000-2009 –> 2010-2019	13.896362	1366.783	1402.449
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + silicate	892.9116	913.1440	6.5879377	2000-2009 –> 2010-2019	13.686193	1359.669	1395.334
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + aou + silicate + phosphate	221.9033	239.1647	0.7992418	2000-2009 –> 2010-2019	4.458538	996.442	1034.295
obs	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	899.1354	919.3678	6.7438983	2000-2009 –> 2010-2019	13.940994	1367.746	1403.411
mod	Atlantic	(-Inf,26]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	222.5504	239.8118	0.8022196	2000-2009 –> 2010-2019	4.693630	1014.309	1052.162
mod	Atlantic	(-Inf,26]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	256.1562	273.4175	0.9731313	2000-2009 –> 2010-2019	4.643551	1031.545	1069.398
obs	Atlantic	(-Inf,26]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	894.1126	914.3451	6.6177504	2000-2009 –> 2010-2019	14.158448	1368.973	1404.638
obs	Atlantic	(-Inf,26]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	893.8311	914.0636	6.6107515	2000-2009 –> 2010-2019	14.144870	1368.574	1404.239
mod	Atlantic	(-Inf,26]	cstar_tref ~ temp + aou + silicate + phosphate	262.7603	277.5558	1.0224613	2000-2009 –> 2010-2019	4.694249	1036.253	1068.698
mod	Atlantic	(-Inf,26]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	259.4053	276.6667	0.9914734	2000-2009 –> 2010-2019	4.629811	1032.336	1070.188
obs	Atlantic	(-Inf,26]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	892.4417	912.6741	6.5763095	2000-2009 –> 2010-2019	14.087951	1366.784	1402.450
mod	Atlantic	(26,26.5]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	2393.0650	2421.5925	3.7270573	2000-2009 –> 2010-2019	8.212230	4719.480	4775.877
mod	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	2353.5687	2382.0961	3.5616392	2000-2009 –> 2010-2019	8.035177	4677.926	4734.324
obs	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	4545.8321	4577.6092	6.3949374	2000-2009 –> 2010-2019	12.079851	9168.694	9232.622
mod	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	2353.8085	2382.3359	3.5626210	2000-2009 –> 2010-2019	8.094049	4688.349	4744.746
obs	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	4554.4813	4586.2584	6.4350269	2000-2009 –> 2010-2019	12.210675	9200.461	9264.390
mod	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate	2307.1838	2335.7113	3.3767205	2000-2009 –> 2010-2019	7.888743	4638.326	4694.723
obs	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate	4645.0682	4676.8453	6.8703083	2000-2009 –> 2010-2019	12.685393	9300.983	9364.912
mod	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	2307.8698	2336.3972	3.3793840	2000-2009 –> 2010-2019	7.997402	4657.402	4713.799
mod	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	2354.6677	2383.1951	3.5661412	2000-2009 –> 2010-2019	8.164602	4700.839	4757.236
obs	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	4582.0081	4613.7852	6.5642964	2000-2009 –> 2010-2019	12.448609	9255.202	9319.130
mod	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	2296.8428	2325.3702	3.3368216	2000-2009 –> 2010-2019	8.241893	4694.136	4750.533
mod	Atlantic	(26,26.5]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	2308.2210	2336.7485	3.3807485	2000-2009 –> 2010-2019	8.033178	4663.633	4720.030
obs	Atlantic	(26,26.5]	cstar_tref ~ temp + aou + nitrate + phosphate	4554.7457	4581.9833	6.4455644	2000-2009 –> 2010-2019	12.139660	9177.964	9232.759
obs	Atlantic	(26,26.5]	cstar_tref ~ temp + aou + nitrate + phosphate_star	4565.3273	4592.5648	6.4950337	2000-2009 –> 2010-2019	12.275020	9210.404	9265.199
mod	Atlantic	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	2389.1576	2417.6851	3.7103556	2000-2009 –> 2010-2019	8.254206	4725.867	4782.264
obs	Atlantic	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	4556.0934	4587.8705	6.4425269	2000-2009 –> 2010-2019	12.106244	9173.501	9237.429
mod	Atlantic	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	2383.2053	2411.7327	3.6850567	2000-2009 –> 2010-2019	8.275241	4727.949	4784.347
obs	Atlantic	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	4566.2228	4597.9999	6.4898522	2000-2009 –> 2010-2019	12.232117	9203.739	9267.668
obs	Atlantic	(26,26.5]	cstar_tref ~ temp + nitrate + phosphate + phosphate_star	4586.2586	4613.4961	6.5940094	2000-2009 –> 2010-2019	12.500020	9262.826	9317.622
obs	Atlantic	(26,26.5]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	4585.0326	4616.8097	6.5786576	2000-2009 –> 2010-2019	12.425589	9248.922	9312.850
obs	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	6962.0882	6996.8569	6.3940801	2000-2009 –> 2010-2019	11.502330	13154.974	13224.187
obs	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	6956.5076	6991.2764	6.3772866	2000-2009 –> 2010-2019	11.492945	13152.329	13221.543
mod	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	3139.0141	3169.1580	4.1885377	2000-2009 –> 2010-2019	7.163124	6224.083	6285.144
obs	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	6986.5426	7021.3114	6.4681932	2000-2009 –> 2010-2019	11.574249	13178.558	13247.771
mod	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	3146.2528	3176.3967	4.2162933	2000-2009 –> 2010-2019	7.192985	6232.188	6293.249
obs	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	6987.3952	7022.1639	6.4707925	2000-2009 –> 2010-2019	11.586464	13183.222	13252.435
mod	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	3108.1984	3138.3423	4.0724110	2000-2009 –> 2010-2019	7.445130	6347.020	6408.081
mod	Atlantic	(26.5,26.75]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	3147.5229	3177.6669	4.2211823	2000-2009 –> 2010-2019	7.208415	6237.785	6298.846
obs	Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + phosphate	6978.4380	7008.2398	6.4496122	2000-2009 –> 2010-2019	11.560960	13170.552	13229.878
obs	Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + phosphate_star	6973.9638	7003.7656	6.4360274	2000-2009 –> 2010-2019	11.557438	13170.062	13229.388
mod	Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	3110.0631	3140.2071	4.0793456	2000-2009 –> 2010-2019	7.180054	6245.960	6307.021
obs	Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	6976.2480	7011.0168	6.4368896	2000-2009 –> 2010-2019	11.544354	13168.822	13238.035
mod	Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	3114.4615	3144.6054	4.0957492	2000-2009 –> 2010-2019	7.177869	6242.998	6304.059
obs	Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	6975.1453	7009.9141	6.4335455	2000-2009 –> 2010-2019	11.553567	13172.694	13241.907
mod	Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + silicate + phosphate	3148.1213	3173.9589	4.2312014	2000-2009 –> 2010-2019	7.352431	6290.092	6342.430
mod	Atlantic	(26.5,26.75]	cstar_tref ~ temp + aou + silicate + phosphate_star	3154.2141	3180.0517	4.2547887	2000-2009 –> 2010-2019	7.373760	6295.299	6347.637
obs	Atlantic	(26.5,26.75]	cstar_tref ~ temp + nitrate + phosphate + phosphate_star	6973.7407	7003.5425	6.4353509	2000-2009 –> 2010-2019	11.558146	13170.387	13229.713
mod	Atlantic	(26.5,26.75]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	3118.7194	3148.8634	4.1116922	2000-2009 –> 2010-2019	7.196800	6248.442	6309.503
obs	Atlantic	(26.5,26.75]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	6975.5096	7010.2784	6.4346502	2000-2009 –> 2010-2019	11.555775	13173.495	13242.708
mod	Atlantic	(26.5,26.75]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	3155.3084	3181.1461	4.2590392	2000-2009 –> 2010-2019	7.385732	6299.420	6351.758
mod	Atlantic	(26.75,27]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	5576.3392	5606.4227	2.9535639	2000-2009 –> 2010-2019	5.571076	10612.592	10672.444
mod	Atlantic	(26.75,27]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	5328.6148	5363.7122	2.6398614	2000-2009 –> 2010-2019	5.277637	10383.140	10452.967
mod	Atlantic	(26.75,27]	cstar_tref ~ sal + aou + silicate + phosphate	5486.6991	5516.7826	2.8368853	2000-2009 –> 2010-2019	5.551820	10600.060	10659.911
mod	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + aou + silicate + phosphate	5452.9193	5488.0168	2.7916107	2000-2009 –> 2010-2019	5.309322	10409.149	10478.975
obs	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + aou + silicate + phosphate	7791.2248	7827.1411	5.4300665	2000-2009 –> 2010-2019	10.965432	14810.380	14881.437
mod	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	5479.8064	5514.9039	2.8255647	2000-2009 –> 2010-2019	5.398042	10481.441	10551.268
obs	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	7781.4234	7817.3397	5.4088194	2000-2009 –> 2010-2019	10.884194	14776.191	14847.249
obs	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	7779.4877	7815.4040	5.4046330	2000-2009 –> 2010-2019	10.877015	14773.032	14844.090
mod	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	5490.4175	5525.5149	2.8390781	2000-2009 –> 2010-2019	5.436482	10512.399	10582.226
obs	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	7781.0955	7817.0118	5.4081099	2000-2009 –> 2010-2019	10.881213	14774.935	14845.992
obs	Atlantic	(26.75,27]	cstar_tref ~ sal + temp + silicate + phosphate_star	7793.1468	7823.9322	5.4385919	2000-2009 –> 2010-2019	10.946080	14799.002	14859.909
mod	Atlantic	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	5167.6502	5202.7476	2.4555487	2000-2009 –> 2010-2019	4.777258	9952.871	10022.698
mod	Atlantic	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	5177.2078	5212.3052	2.4661241	2000-2009 –> 2010-2019	4.794813	9968.763	10038.589
obs	Atlantic	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	7779.3597	7815.2760	5.4043563	2000-2009 –> 2010-2019	10.925004	14792.556	14863.614
mod	Atlantic	(26.75,27]	cstar_tref ~ temp + aou + silicate + phosphate	5627.8859	5657.9694	3.0228195	2000-2009 –> 2010-2019	5.666936	10685.478	10745.329
obs	Atlantic	(26.75,27]	cstar_tref ~ temp + aou + silicate + phosphate_star	7779.4535	7810.2389	5.4088845	2000-2009 –> 2010-2019	10.964576	14804.812	14865.718
mod	Atlantic	(26.75,27]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	5201.7861	5236.8835	2.4935293	2000-2009 –> 2010-2019	4.849948	10018.318	10088.144
obs	Atlantic	(26.75,27]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	7779.4913	7815.4076	5.4046409	2000-2009 –> 2010-2019	10.933396	14795.972	14867.030
obs	Atlantic	(26.75,27]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	7777.5624	7808.3478	5.4047944	2000-2009 –> 2010-2019	10.948435	14798.061	14858.967
obs	Atlantic	(26.75,27]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	7779.1318	7809.9172	5.4081884	2000-2009 –> 2010-2019	10.962144	14803.791	14864.697
mod	Atlantic	(27,27.25]	cstar_tref ~ sal + aou + nitrate + silicate	6872.4998	6904.4111	2.3531692	2000-2009 –> 2010-2019	4.960853	13717.150	13780.683
mod	Atlantic	(27,27.25]	cstar_tref ~ sal + aou + silicate + phosphate	6825.5799	6857.4912	2.3168442	2000-2009 –> 2010-2019	4.794204	13522.884	13586.417
obs	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	13310.6882	13350.1184	6.0532696	2000-2009 –> 2010-2019	11.653858	26180.975	26259.771
mod	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + aou + nitrate + silicate	6620.3312	6657.5610	2.1629858	2000-2009 –> 2010-2019	4.476941	13123.526	13197.648
mod	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + aou + silicate + phosphate	6827.3717	6864.6014	2.3166842	2000-2009 –> 2010-2019	4.762213	13489.508	13563.630
mod	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	6923.6872	6960.9170	2.3918612	2000-2009 –> 2010-2019	4.863055	13615.828	13689.950
obs	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	13346.9909	13386.4211	6.1067124	2000-2009 –> 2010-2019	11.605389	26142.132	26220.928
obs	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	13323.4866	13362.9168	6.0720571	2000-2009 –> 2010-2019	11.690360	26206.696	26285.492
mod	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	7100.7036	7137.9334	2.5364469	2000-2009 –> 2010-2019	5.220235	14030.030	14104.152
obs	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	13372.4096	13411.8398	6.1444130	2000-2009 –> 2010-2019	11.679137	26194.289	26273.084
mod	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	6947.8610	6985.0907	2.4111094	2000-2009 –> 2010-2019	4.896615	13656.598	13730.720
obs	Atlantic	(27,27.25]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	13347.9424	13387.3726	6.1081195	2000-2009 –> 2010-2019	11.613293	26147.915	26226.711
mod	Atlantic	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate	6835.1515	6867.0628	2.3242087	2000-2009 –> 2010-2019	4.702640	13415.333	13478.866
mod	Atlantic	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	6835.7561	6872.9859	2.3231336	2000-2009 –> 2010-2019	4.680256	13392.072	13466.194
obs	Atlantic	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	13287.4172	13326.8474	6.0192577	2000-2009 –> 2010-2019	11.663472	26189.456	26268.251
mod	Atlantic	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	6837.1333	6874.3631	2.3241946	2000-2009 –> 2010-2019	4.667118	13376.085	13450.207
obs	Atlantic	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	13244.5771	13284.0073	5.9571432	2000-2009 –> 2010-2019	11.516173	26084.387	26163.183
obs	Atlantic	(27,27.25]	cstar_tref ~ temp + aou + silicate + phosphate_star	13349.0125	13382.8098	6.1126617	2000-2009 –> 2010-2019	11.690080	26200.336	26267.875
obs	Atlantic	(27,27.25]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	13255.2644	13294.6946	5.9725785	2000-2009 –> 2010-2019	11.545979	26105.639	26184.434
obs	Atlantic	(27,27.25]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	13349.9397	13383.7370	6.1140342	2000-2009 –> 2010-2019	11.699279	26207.001	26274.540
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + aou + nitrate + silicate	14747.7245	14783.2188	3.5610743	2000-2009 –> 2010-2019	6.940302	29788.262	29859.487
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + aou + silicate + phosphate	14537.3457	14572.8400	3.4269546	2000-2009 –> 2010-2019	7.086255	30031.746	30102.971
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	16538.4436	16579.5026	5.7633715	2000-2009 –> 2010-2019	10.250372	31153.059	31234.886
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	16541.1749	16582.2339	5.7663926	2000-2009 –> 2010-2019	10.219733	31118.140	31199.967
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + nitrate + silicate	14706.2446	14747.6546	3.5329316	2000-2009 –> 2010-2019	6.827928	29604.903	29687.999
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + phosphate	16574.8896	16610.0831	5.8060419	2000-2009 –> 2010-2019	10.306799	31202.810	31272.948
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + phosphate_star	16576.1514	16611.3449	5.8074477	2000-2009 –> 2010-2019	10.280962	31173.715	31243.853
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate	14538.2838	14579.6938	3.4262906	2000-2009 –> 2010-2019	7.074272	30017.025	30100.121
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate	16565.8961	16606.9551	5.7938083	2000-2009 –> 2010-2019	10.293017	31194.097	31275.924
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	14390.1822	14431.5922	3.3349323	2000-2009 –> 2010-2019	6.903831	29743.998	29827.093
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	16569.4901	16610.5491	5.7978049	2000-2009 –> 2010-2019	10.271314	31169.049	31250.876
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	16534.6755	16575.7345	5.7592063	2000-2009 –> 2010-2019	10.220844	31120.948	31202.775
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + nitrate + phosphate_star	14751.5010	14786.9953	3.5635293	2000-2009 –> 2010-2019	7.109980	30067.351	30138.575
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	14359.5396	14400.9495	3.3163363	2000-2009 –> 2010-2019	6.819022	29606.611	29689.706
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + phosphate + phosphate_star	14772.6995	14808.1938	3.5773410	2000-2009 –> 2010-2019	7.209778	30225.100	30296.325
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + phosphate + phosphate_star	16572.2366	16607.4300	5.8030872	2000-2009 –> 2010-2019	10.281498	31175.271	31245.409
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	14392.4220	14433.8320	3.3362956	2000-2009 –> 2010-2019	6.905336	29746.463	29829.558
obs	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	16564.9750	16606.0340	5.7927844	2000-2009 –> 2010-2019	10.270976	31169.764	31251.592
mod	Atlantic	(27.25,27.5]	cstar_tref ~ sal + temp + silicate + phosphate_star	14512.2613	14547.7556	3.4113038	2000-2009 –> 2010-2019	7.049432	29973.586	30044.811
obs	Atlantic	(27.25,27.5]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	16560.8362	16601.8952	5.7881863	2000-2009 –> 2010-2019	10.400923	31313.635	31395.462
obs	Atlantic	(27.5,27.75]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	18264.0968	18306.1133	5.1292521	2000-2009 –> 2010-2019	10.007905	35892.425	35976.325
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + aou + phosphate	24906.2982	24939.1585	2.5542861	2000-2009 –> 2010-2019	4.995759	48840.737	48906.356
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + aou + silicate + phosphate	24033.4900	24072.9224	2.3512861	2000-2009 –> 2010-2019	4.741816	47751.684	47830.426
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + aou + phosphate	24904.9358	24944.3681	2.5534732	2000-2009 –> 2010-2019	4.977937	48769.015	48847.757
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + aou + phosphate_star	24730.7163	24770.1487	2.5117072	2000-2009 –> 2010-2019	4.910291	48483.721	48562.463
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	24025.6303	24071.6347	2.3490926	2000-2009 –> 2010-2019	4.728100	47695.815	47787.681
obs	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	18255.3615	18297.3781	5.1217600	2000-2009 –> 2010-2019	9.832473	35678.274	35762.175
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	23773.6074	23819.6118	2.2937140	2000-2009 –> 2010-2019	4.644567	47320.577	47412.443
obs	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	18268.7406	18310.7571	5.1332394	2000-2009 –> 2010-2019	9.857773	35708.832	35792.733
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	23911.2024	23957.2068	2.3237848	2000-2009 –> 2010-2019	4.682422	47492.385	47584.250
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + phosphate + phosphate_star	24732.7376	24772.1699	2.5121878	2000-2009 –> 2010-2019	4.910632	48485.142	48563.884
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	23782.5140	23828.5184	2.2956487	2000-2009 –> 2010-2019	4.646832	47330.937	47422.803
obs	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	18255.8114	18297.8279	5.1221456	2000-2009 –> 2010-2019	9.834982	35681.367	35765.268
mod	Atlantic	(27.5,27.75]	cstar_tref ~ sal + temp + silicate + phosphate_star	24489.1964	24528.6288	2.4549346	2000-2009 –> 2010-2019	4.961089	48696.268	48775.010
obs	Atlantic	(27.5,27.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	18226.5288	18268.5453	5.0971083	2000-2009 –> 2010-2019	9.860255	35714.352	35798.253
obs	Atlantic	(27.5,27.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	18249.1893	18291.2058	5.1164728	2000-2009 –> 2010-2019	9.900830	35763.067	35846.968
obs	Atlantic	(27.5,27.75]	cstar_tref ~ temp + aou + silicate + phosphate	18253.6085	18289.6227	5.1219717	2000-2009 –> 2010-2019	9.886071	35740.604	35812.518
obs	Atlantic	(27.5,27.75]	cstar_tref ~ temp + aou + silicate + phosphate_star	18266.8713	18302.8855	5.1333518	2000-2009 –> 2010-2019	9.918516	35779.738	35851.653
obs	Atlantic	(27.5,27.75]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	18228.6680	18270.6846	5.0989333	2000-2009 –> 2010-2019	9.867198	35722.789	35806.690
obs	Atlantic	(27.5,27.75]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	18253.9910	18290.0052	5.1222996	2000-2009 –> 2010-2019	9.891843	35747.684	35819.598
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + nitrate + silicate	6918.0340	6951.5121	1.4114314	2000-2009 –> 2010-2019	2.725508	13755.855	13822.989
obs	Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	11889.4401	11928.2542	5.5904005	2000-2009 –> 2010-2019	10.291673	23533.185	23611.064
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + silicate + phosphate	6639.8777	6673.3558	1.3146542	2000-2009 –> 2010-2019	2.583873	13337.588	13404.722
obs	Atlantic	(27.75,27.85]	cstar_tref ~ sal + aou + silicate + phosphate	11888.9518	11922.2210	5.5926354	2000-2009 –> 2010-2019	10.342787	23571.243	23637.996
obs	Atlantic	(27.75,27.85]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	11879.2671	11918.0811	5.5753833	2000-2009 –> 2010-2019	10.288008	23532.467	23610.346
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + nitrate + silicate + phosphate_star	7346.1454	7379.6235	1.5744846	2000-2009 –> 2010-2019	3.095319	14773.435	14840.569
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	7333.9042	7367.3823	1.5695706	2000-2009 –> 2010-2019	3.092267	14766.130	14833.264
obs	Atlantic	(27.75,27.85]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	11878.5298	11911.7990	5.5772451	2000-2009 –> 2010-2019	10.344444	23574.864	23641.617
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + aou + nitrate + silicate	6881.0103	6920.0680	1.3974361	2000-2009 –> 2010-2019	2.705238	13701.523	13779.846
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + aou + silicate + phosphate	6525.9888	6565.0466	1.2763188	2000-2009 –> 2010-2019	2.511251	13115.226	13193.550
obs	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + aou + silicate + phosphate	11882.7751	11921.5891	5.5805571	2000-2009 –> 2010-2019	10.329740	23566.267	23644.146
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	6668.3189	6707.3767	1.3235609	2000-2009 –> 2010-2019	2.609258	13420.066	13498.389
obs	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	11868.7486	11907.5626	5.5598986	2000-2009 –> 2010-2019	10.315871	23557.840	23635.719
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	6676.6906	6715.7483	1.3263934	2000-2009 –> 2010-2019	2.611871	13427.749	13506.073
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	6660.3063	6699.3641	1.3208555	2000-2009 –> 2010-2019	2.603950	13403.877	13482.200
obs	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	11879.1769	11917.9910	5.5752504	2000-2009 –> 2010-2019	10.325983	23563.948	23641.827
mod	Atlantic	(27.75,27.85]	cstar_tref ~ sal + temp + silicate + phosphate_star	6797.2208	6830.6988	1.3685519	2000-2009 –> 2010-2019	2.681833	13632.600	13699.734
obs	Atlantic	(27.75,27.85]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	11907.5260	11946.3401	5.6171983	2000-2009 –> 2010-2019	10.324021	23555.897	23633.776
obs	Atlantic	(27.75,27.85]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	11892.9805	11931.7945	5.5956362	2000-2009 –> 2010-2019	10.305920	23544.232	23622.111
obs	Atlantic	(27.75,27.85]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	11903.7523	11942.5664	5.6115963	2000-2009 –> 2010-2019	10.318440	23552.140	23630.019
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	24685.6561	24729.7511	5.2013191	2000-2009 –> 2010-2019	10.422374	50439.983	50528.461
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	23995.8642	24039.9592	4.7737807	2000-2009 –> 2010-2019	9.704907	49271.425	49359.903
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + aou + silicate + phosphate	9356.4160	9389.6501	2.9045500	2000-2009 –> 2010-2019	5.557458	19291.618	19358.670
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + aou + silicate + phosphate_star	24005.9304	24043.7262	4.7809487	2000-2009 –> 2010-2019	9.712098	49279.529	49355.367
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	24432.7685	24476.8635	5.0403040	2000-2009 –> 2010-2019	10.140560	49990.929	50079.407
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	24481.4171	24519.2128	5.0721481	2000-2009 –> 2010-2019	10.176791	50044.783	50120.621
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + aou + nitrate + silicate	9466.8446	9505.6178	2.9895291	2000-2009 –> 2010-2019	5.803419	19648.175	19726.402
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + aou + phosphate	9487.9112	9521.1453	3.0079254	2000-2009 –> 2010-2019	5.688534	19466.158	19533.209
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + aou + phosphate_star	9568.4355	9601.6697	3.0730381	2000-2009 –> 2010-2019	5.806364	19627.388	19694.440
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate + phosphate	9212.0510	9250.8242	2.7936574	2000-2009 –> 2010-2019	5.401902	19078.892	19157.119
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate + phosphate	23901.2314	23945.3264	4.7179353	2000-2009 –> 2010-2019	9.627180	49139.524	49228.002
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	9263.2407	9302.0139	2.8319511	2000-2009 –> 2010-2019	5.477825	19189.439	19267.666
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	23850.4686	23894.5636	4.6882485	2000-2009 –> 2010-2019	9.581534	49061.474	49149.952
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate	23957.7934	24001.8884	4.7512350	2000-2009 –> 2010-2019	9.734863	49322.104	49410.583
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	9278.0511	9316.8243	2.8431280	2000-2009 –> 2010-2019	5.508399	19234.520	19312.747
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + phosphate + phosphate_star	9548.7709	9582.0050	3.0570081	2000-2009 –> 2010-2019	5.775869	19585.735	19652.786
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + silicate + phosphate	23975.2644	24013.0601	4.7627525	2000-2009 –> 2010-2019	9.757315	49355.941	49431.780
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	9256.7633	9295.5365	2.8270767	2000-2009 –> 2010-2019	5.465814	19171.770	19249.997
obs	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	23893.3012	23937.3963	4.7132853	2000-2009 –> 2010-2019	9.618415	49124.566	49213.044
mod	Atlantic	(27.85,27.95]	cstar_tref ~ sal + temp + silicate + phosphate_star	9335.6563	9368.8904	2.8885576	2000-2009 –> 2010-2019	5.638082	19419.095	19486.147
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	26679.7154	26724.1817	5.6159085	2000-2009 –> 2010-2019	11.693422	54899.639	54988.791
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	26202.1030	26246.5692	5.3083505	2000-2009 –> 2010-2019	11.194001	54141.341	54230.493
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + aou + silicate + phosphate_star	26270.0929	26308.2068	5.3523445	2000-2009 –> 2010-2019	11.255449	54233.241	54309.656
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	26346.8862	26391.3525	5.3997608	2000-2009 –> 2010-2019	11.329967	54352.114	54441.266
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + nitrate + silicate + phosphate_star	26614.8569	26652.9708	5.5744342	2000-2009 –> 2010-2019	11.716908	54925.811	55002.227
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	26436.0023	26474.1162	5.4580930	2000-2009 –> 2010-2019	11.404404	54462.961	54539.377
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + aou + phosphate	11720.4879	11754.6734	3.4502881	2000-2009 –> 2010-2019	7.547934	24866.336	24935.020
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + aou + phosphate_star	11753.4094	11787.5948	3.4761651	2000-2009 –> 2010-2019	7.582732	24909.352	24978.036
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate + phosphate	11575.3932	11615.2763	3.3370015	2000-2009 –> 2010-2019	7.411516	24696.963	24777.094
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate + phosphate	26362.2824	26406.7487	5.4095734	2000-2009 –> 2010-2019	11.344579	54374.589	54463.741
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	11586.8024	11626.6854	3.3456537	2000-2009 –> 2010-2019	7.422285	24710.783	24790.914
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	26266.6058	26311.0721	5.3488823	2000-2009 –> 2010-2019	11.251113	54230.459	54319.611
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + nitrate + phosphate_star	11843.5720	11877.7574	3.5480327	2000-2009 –> 2010-2019	7.685616	25034.444	25103.128
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	11578.6039	11618.4870	3.3394341	2000-2009 –> 2010-2019	7.379056	24660.250	24740.381
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	26615.1469	26659.6132	5.5733102	2000-2009 –> 2010-2019	11.706815	54915.314	55004.466
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + phosphate + phosphate_star	11747.5525	11781.7379	3.4715473	2000-2009 –> 2010-2019	7.577166	24902.424	24971.108
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + phosphate_star	11880.9823	11909.4702	3.5799111	2000-2009 –> 2010-2019	7.717513	25069.875	25127.112
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	11586.5282	11626.4112	3.3454455	2000-2009 –> 2010-2019	7.422074	24710.505	24790.636
obs	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	26326.9588	26371.4251	5.3870866	2000-2009 –> 2010-2019	11.306983	54316.961	54406.113
mod	Atlantic	(27.95,28.05]	cstar_tref ~ sal + temp + silicate + phosphate_star	11585.9265	11620.1119	3.3465074	2000-2009 –> 2010-2019	7.423139	24707.908	24776.592
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	21500.5795	21543.5228	5.6443931	2000-2009 –> 2010-2019	12.041563	42551.857	42637.323
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + aou + silicate + phosphate_star	21551.0423	21587.8509	5.6879672	2000-2009 –> 2010-2019	12.093619	42608.832	42682.089
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	21596.5368	21633.3454	5.7260258	2000-2009 –> 2010-2019	12.146418	42669.093	42742.349
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	21571.8186	21614.7619	5.7036438	2000-2009 –> 2010-2019	12.109157	42631.470	42716.935
mod	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + aou + nitrate + silicate	4364.2435	4400.2045	1.3635520	2000-2009 –> 2010-2019	2.703987	8805.301	8877.415
mod	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + aou + phosphate	4177.1197	4207.9434	1.2668259	2000-2009 –> 2010-2019	2.507369	8415.906	8477.718
mod	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + aou + phosphate_star	4228.9181	4259.7418	1.2931771	2000-2009 –> 2010-2019	2.553554	8508.721	8570.533
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + aou + phosphate_star	21587.6861	21624.4947	5.7186018	2000-2009 –> 2010-2019	12.124763	42645.987	42719.243
mod	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + aou + silicate + phosphate	4172.9474	4208.9084	1.2637219	2000-2009 –> 2010-2019	2.497733	8400.081	8472.195
mod	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	4219.2387	4255.1996	1.2871880	2000-2009 –> 2010-2019	2.543590	8492.873	8564.987
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	21503.1852	21546.1285	5.6465494	2000-2009 –> 2010-2019	12.017517	42528.097	42613.563
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	21562.4571	21605.4004	5.6958223	2000-2009 –> 2010-2019	12.108309	42629.097	42714.563
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + nitrate + phosphate_star	21569.8980	21606.7065	5.7037102	2000-2009 –> 2010-2019	12.117629	42635.974	42709.230
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	21506.0355	21548.9788	5.6489091	2000-2009 –> 2010-2019	12.051189	42562.443	42647.909
obs	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	21571.0033	21613.9466	5.7029622	2000-2009 –> 2010-2019	12.104990	42627.157	42712.623
mod	Atlantic	(28.05,28.1]	cstar_tref ~ sal + temp + silicate + phosphate_star	4364.0634	4394.8871	1.3645387	2000-2009 –> 2010-2019	2.803700	8986.900	9048.712
mod	Atlantic	(28.05,28.1]	cstar_tref ~ temp + aou + phosphate	4222.9847	4248.6711	1.2911570	2000-2009 –> 2010-2019	2.565672	8529.635	8581.145
mod	Atlantic	(28.05,28.1]	cstar_tref ~ temp + aou + phosphate_star	4276.2970	4301.9834	1.3188076	2000-2009 –> 2010-2019	2.616952	8630.452	8681.962
mod	Atlantic	(28.05,28.1]	cstar_tref ~ temp + aou + silicate + phosphate	4224.9722	4255.7959	1.2911506	2000-2009 –> 2010-2019	2.555423	8512.756	8574.568
mod	Atlantic	(28.05,28.1]	cstar_tref ~ temp + aou + silicate + phosphate_star	4277.4411	4308.2648	1.3183591	2000-2009 –> 2010-2019	2.609350	8619.306	8681.118
mod	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + aou + nitrate	6324.1899	6357.6095	1.2321413	2000-2009 –> 2010-2019	2.586005	13141.143	13208.095
obs	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	16952.3511	16993.5194	5.9334088	2000-2009 –> 2010-2019	12.860154	34039.884	34121.945
mod	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + aou + nitrate + silicate	6159.7128	6198.7023	1.1803665	2000-2009 –> 2010-2019	2.531782	12971.514	13049.626
mod	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + aou + phosphate	5839.4547	5872.8743	1.0873650	2000-2009 –> 2010-2019	2.419482	12592.412	12659.365
mod	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + aou + phosphate_star	5898.0032	5931.4227	1.1039062	2000-2009 –> 2010-2019	2.445302	12678.392	12745.344
obs	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + aou + phosphate_star	16956.8689	16992.1560	5.9407183	2000-2009 –> 2010-2019	12.867925	34042.741	34113.079
mod	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + aou + silicate + phosphate	5800.8382	5839.8277	1.0760359	2000-2009 –> 2010-2019	2.314995	12269.285	12347.396
mod	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	5870.9603	5909.9498	1.0956698	2000-2009 –> 2010-2019	2.349369	12386.148	12464.259
obs	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	16945.4502	16986.6185	5.9256794	2000-2009 –> 2010-2019	12.835128	34020.257	34102.318
obs	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	16938.9110	16980.0793	5.9183645	2000-2009 –> 2010-2019	12.828025	34013.874	34095.935
obs	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + phosphate + phosphate_star	16950.4565	16985.7436	5.9335269	2000-2009 –> 2010-2019	12.859204	34035.204	34105.542
obs	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	16946.3499	16987.5182	5.9266866	2000-2009 –> 2010-2019	12.840125	34024.095	34106.156
obs	Atlantic	(28.1,28.15]	cstar_tref ~ sal + temp + silicate + phosphate_star	16944.4996	16979.7866	5.9268541	2000-2009 –> 2010-2019	12.844501	34023.342	34093.681
obs	Atlantic	(28.1,28.15]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	16949.3375	16990.5058	5.9300322	2000-2009 –> 2010-2019	12.860913	34039.908	34121.969
mod	Atlantic	(28.1,28.15]	cstar_tref ~ temp + aou + phosphate	5944.8169	5972.6665	1.1178892	2000-2009 –> 2010-2019	2.481809	12789.015	12844.808
mod	Atlantic	(28.1,28.15]	cstar_tref ~ temp + aou + phosphate_star	6020.6281	6048.4777	1.1399579	2000-2009 –> 2010-2019	2.517756	12904.838	12960.632
mod	Atlantic	(28.1,28.15]	cstar_tref ~ temp + aou + silicate + phosphate	5925.1399	5958.5594	1.1116580	2000-2009 –> 2010-2019	2.421276	12610.781	12677.734
mod	Atlantic	(28.1,28.15]	cstar_tref ~ temp + aou + silicate + phosphate_star	6009.7786	6043.1982	1.1361870	2000-2009 –> 2010-2019	2.465857	12755.469	12822.422
obs	Atlantic	(28.1,28.15]	cstar_tref ~ temp + aou + silicate + phosphate_star	16954.7039	16989.9910	5.9382894	2000-2009 –> 2010-2019	12.869734	34043.689	34114.027
obs	Atlantic	(28.1,28.15]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	16948.9992	16984.2863	5.9318938	2000-2009 –> 2010-2019	12.865177	34039.334	34109.672
mod	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + aou + nitrate	9014.7749	9051.6010	0.9020114	2000-2009 –> 2010-2019	1.909723	18907.548	18981.274
obs	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	6374.3065	6409.7359	3.7004915	2000-2009 –> 2010-2019	7.460486	12094.487	12164.545
mod	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + aou + phosphate	9481.2281	9518.0542	0.9656507	2000-2009 –> 2010-2019	2.033756	19777.128	19850.854
obs	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + aou + phosphate	6348.9897	6379.3577	3.6636759	2000-2009 –> 2010-2019	7.393197	12050.244	12110.293
obs	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + aou + silicate + phosphate	6333.1946	6368.6240	3.6358255	2000-2009 –> 2010-2019	7.342726	12023.794	12093.853
mod	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	9497.8924	9540.8563	0.9677226	2000-2009 –> 2010-2019	2.035440	19793.273	19879.286
obs	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	6348.9640	6384.3933	3.6604949	2000-2009 –> 2010-2019	7.393638	12054.238	12124.296
mod	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + nitrate + phosphate_star	9163.2830	9200.1091	0.9218039	2000-2009 –> 2010-2019	1.952209	19210.300	19284.025
mod	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	9150.3557	9193.3196	0.9197950	2000-2009 –> 2010-2019	1.945180	19165.542	19251.555
obs	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + phosphate + phosphate_star	6368.1714	6398.5394	3.6939356	2000-2009 –> 2010-2019	7.449460	12083.878	12143.928
mod	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	9450.7601	9493.7239	0.9610791	2000-2009 –> 2010-2019	2.020940	19694.992	19781.005
obs	Atlantic	(28.15,28.2]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	6343.2455	6378.6748	3.6515296	2000-2009 –> 2010-2019	7.374290	12042.726	12112.784
mod	Atlantic	(28.15,28.2]	cstar_tref ~ temp + aou + nitrate	9249.7633	9280.4518	0.9338020	2000-2009 –> 2010-2019	1.979182	19394.711	19456.149
mod	Atlantic	(28.15,28.2]	cstar_tref ~ temp + aou + nitrate + silicate	8883.5491	8920.3752	0.8848762	2000-2009 –> 2010-2019	1.894461	18789.187	18862.912
obs	Atlantic	(28.15,28.2]	cstar_tref ~ temp + aou + phosphate	6362.9371	6388.2437	3.6888158	2000-2009 –> 2010-2019	7.422502	12064.513	12114.554
obs	Atlantic	(28.15,28.2]	cstar_tref ~ temp + aou + silicate + phosphate	6346.3167	6376.6847	3.6594790	2000-2009 –> 2010-2019	7.370243	12037.083	12097.133
obs	Atlantic	(28.15,28.2]	cstar_tref ~ temp + aou + silicate + phosphate_star	6367.7317	6398.0997	3.6932391	2000-2009 –> 2010-2019	7.436070	12076.396	12136.446
mod	Atlantic	(28.15,28.2]	cstar_tref ~ temp + nitrate + phosphate_star	9328.6573	9359.3457	0.9446319	2000-2009 –> 2010-2019	1.982072	19420.800	19482.238
mod	Atlantic	(28.15,28.2]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	9201.7034	9238.5296	0.9269948	2000-2009 –> 2010-2019	1.952622	19216.527	19290.252
obs	Atlantic	(28.15,28.2]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	6360.5235	6390.8915	3.6818410	2000-2009 –> 2010-2019	7.412686	12062.516	12122.566
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate	13408.6945	13442.6061	0.6763017	2000-2009 –> 2010-2019	1.356597	24448.387	24515.209
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + aou + nitrate + silicate	13262.7537	13303.4476	0.6686700	2000-2009 –> 2010-2019	1.326939	23953.264	24033.450
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + nitrate + phosphate_star	16366.7176	16400.6292	0.8485688	2000-2009 –> 2010-2019	1.749457	30403.122	30469.944
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	14117.9746	14159.3282	3.2401240	2000-2009 –> 2010-2019	6.490172	26668.032	26749.906
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + nitrate + silicate + phosphate_star	16336.7926	16377.4866	0.8464932	2000-2009 –> 2010-2019	1.734302	30219.153	30299.339
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + aou + nitrate	13345.0866	13385.7805	0.6729065	2000-2009 –> 2010-2019	1.352322	24372.963	24453.149
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	14074.5551	14115.9087	3.2143470	2000-2009 –> 2010-2019	6.484972	26655.071	26736.944
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	14118.6651	14160.0186	3.2405356	2000-2009 –> 2010-2019	6.484545	26659.748	26741.622
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + aou + nitrate + silicate	13210.7065	13258.1827	0.6659034	2000-2009 –> 2010-2019	1.318961	23818.404	23911.955
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + aou + phosphate	14074.3475	14109.7934	3.2154070	2000-2009 –> 2010-2019	6.493675	26664.128	26734.306
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + aou + phosphate_star	14117.0780	14152.5239	3.2407818	2000-2009 –> 2010-2019	6.485759	26657.600	26727.778
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + aou + silicate + phosphate	14045.2872	14086.6408	3.1970872	2000-2009 –> 2010-2019	6.474746	26636.171	26718.045
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	14031.1048	14072.4583	3.1887569	2000-2009 –> 2010-2019	6.430239	26568.426	26650.300
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	14091.6063	14132.9599	3.2244453	2000-2009 –> 2010-2019	6.471209	26636.785	26718.659
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + nitrate + phosphate_star	15809.0041	15849.6981	0.8129057	2000-2009 –> 2010-2019	1.562715	27888.792	27968.978
mod	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	15777.9970	15825.4732	0.8108501	2000-2009 –> 2010-2019	1.559200	27839.009	27932.560
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + phosphate + phosphate_star	14089.6077	14125.0536	3.2244461	2000-2009 –> 2010-2019	6.473999	26636.930	26707.108
obs	Atlantic	(28.2, Inf]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	14019.3735	14060.7271	3.1818827	2000-2009 –> 2010-2019	6.430692	26567.591	26649.465
mod	Atlantic	(28.2, Inf]	cstar_tref ~ temp + aou + nitrate + silicate	15327.0790	15367.7729	0.7834023	2000-2009 –> 2010-2019	1.586070	28133.732	28213.918
mod	Atlantic	(28.2, Inf]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	16779.3874	16820.0814	0.8757266	2000-2009 –> 2010-2019	1.723521	30169.661	30249.847
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	26326.0798	26370.6607	5.1218050	2000-2009 –> 2010-2019	10.436245	48101.576	48189.324
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	34295.4773	34341.8055	5.3630452	2000-2009 –> 2010-2019	11.303151	63365.119	63456.390
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + aou + phosphate + phosphate_star	27255.2758	27293.4880	5.7060916	2000-2009 –> 2010-2019	11.469459	49599.999	49675.212
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	26934.1887	26978.7696	5.4961794	2000-2009 –> 2010-2019	11.029732	48994.279	49082.027
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	35493.9695	35540.2976	5.9766220	2000-2009 –> 2010-2019	12.232596	65033.941	65125.213
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate	27384.0748	27422.2870	5.7919914	2000-2009 –> 2010-2019	11.956505	50202.770	50277.983
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate	34613.5291	34653.2390	5.5204493	2000-2009 –> 2010-2019	11.624128	63927.772	64006.004
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	26762.8991	26807.4800	5.3880418	2000-2009 –> 2010-2019	11.035077	48965.987	49053.736
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	34519.1274	34565.4556	5.4725581	2000-2009 –> 2010-2019	11.469412	63675.083	63766.354
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	26854.1127	26898.6936	5.4453587	2000-2009 –> 2010-2019	11.174306	49158.641	49246.390
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	34536.4077	34582.7359	5.4811120	2000-2009 –> 2010-2019	11.488220	63707.873	63799.144
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + silicate	26903.5569	26948.1378	5.4766830	2000-2009 –> 2010-2019	11.285900	49306.096	49393.845
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + nitrate + silicate	34613.6310	34659.9591	5.5195023	2000-2009 –> 2010-2019	11.622259	63928.501	64019.773
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + silicate + phosphate	27655.0072	27699.5881	5.9755422	2000-2009 –> 2010-2019	12.071187	50396.565	50484.313
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + silicate + phosphate	35736.3558	35782.6839	6.1090006	2000-2009 –> 2010-2019	12.493976	65461.616	65552.887
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	27825.2072	27869.7881	6.0947004	2000-2009 –> 2010-2019	12.359469	50759.539	50847.288
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	35803.0724	35849.4006	6.1459495	2000-2009 –> 2010-2019	12.571900	65586.404	65677.676
mod	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	27001.2907	27045.8716	5.5391311	2000-2009 –> 2010-2019	11.376616	49438.023	49525.772
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	34668.5663	34714.8945	5.5469761	2000-2009 –> 2010-2019	11.607017	63919.674	64010.945
obs	Indo-Pacific	(-Inf,26]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	35831.4918	35877.8199	6.1617565	2000-2009 –> 2010-2019	12.606763	65641.704	65732.976
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	19576.4414	19619.2500	4.5010887	2000-2009 –> 2010-2019	8.798270	37594.068	37679.216
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	20085.8884	20128.7429	4.7628912	2000-2009 –> 2010-2019	9.413708	38385.589	38470.702
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	19841.0924	19883.9010	4.6826616	2000-2009 –> 2010-2019	9.170631	38130.573	38215.721
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	20118.5768	20161.4313	4.7860607	2000-2009 –> 2010-2019	9.454062	38441.093	38526.206
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + phosphate	20494.5158	20531.2483	5.0622700	2000-2009 –> 2010-2019	10.013876	39179.888	39252.841
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate	18744.2435	18787.0521	3.9747512	2000-2009 –> 2010-2019	7.865010	36139.302	36224.450
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate	20255.2860	20298.1405	4.8841876	2000-2009 –> 2010-2019	9.649559	38705.508	38790.621
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	18997.8547	19040.6633	4.1282758	2000-2009 –> 2010-2019	8.197299	36674.067	36759.215
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	20303.8085	20346.6631	4.9194977	2000-2009 –> 2010-2019	9.721892	38801.905	38887.018
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	19893.0902	19935.8988	4.7191883	2000-2009 –> 2010-2019	9.258392	38253.607	38338.755
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	20125.3353	20168.1898	4.7908651	2000-2009 –> 2010-2019	9.464608	38455.455	38540.568
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + nitrate + phosphate_star	20251.9268	20288.6592	4.8832022	2000-2009 –> 2010-2019	9.746666	38826.192	38899.145
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	20170.6008	20213.4553	4.8231679	2000-2009 –> 2010-2019	9.635126	38681.005	38766.117
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	19045.7872	19088.5958	4.1579514	2000-2009 –> 2010-2019	8.265390	36780.805	36865.953
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	20302.8570	20345.7115	4.9188028	2000-2009 –> 2010-2019	9.723176	38803.502	38888.615
obs	Indo-Pacific	(26,26.5]	cstar_tref ~ sal + temp + silicate + phosphate_star	20305.4887	20342.2211	4.9221862	2000-2009 –> 2010-2019	9.750948	38835.456	38908.410
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	19329.6466	19372.4553	4.3381163	2000-2009 –> 2010-2019	8.611298	37312.225	37397.373
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	19466.3728	19509.1814	4.4276612	2000-2009 –> 2010-2019	8.743764	37511.494	37596.642
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ temp + aou + silicate + phosphate	19711.8300	19748.5231	4.5944524	2000-2009 –> 2010-2019	9.131940	38067.981	38140.965
mod	Indo-Pacific	(26,26.5]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	19516.6637	19559.4723	4.4610608	2000-2009 –> 2010-2019	8.804606	37601.449	37686.597
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	19327.0614	19369.6728	4.7094048	2000-2009 –> 2010-2019	9.139825	35866.911	35951.188
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	21694.4907	21737.6063	5.3753202	2000-2009 –> 2010-2019	10.431126	41130.742	41216.343
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	21767.6535	21810.7691	5.4318617	2000-2009 –> 2010-2019	10.545849	41277.019	41362.620
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + aou + phosphate	22186.9765	22223.9327	5.7692368	2000-2009 –> 2010-2019	11.173248	42046.828	42120.200
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	18841.1899	18883.8012	4.3705157	2000-2009 –> 2010-2019	8.518413	35006.503	35090.780
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	22182.9785	22226.0942	5.7642899	2000-2009 –> 2010-2019	11.160216	42035.263	42120.864
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	19034.5467	19077.1581	4.5023558	2000-2009 –> 2010-2019	8.766338	35356.751	35441.028
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	21776.2538	21819.3694	5.4385471	2000-2009 –> 2010-2019	10.562217	41297.707	41383.308
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + nitrate + phosphate_star	21847.2510	21884.2073	5.4956236	2000-2009 –> 2010-2019	10.729878	41503.153	41576.525
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	19257.1980	19299.8094	4.6591044	2000-2009 –> 2010-2019	9.153406	35878.418	35962.695
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	21803.3199	21846.4355	5.4596406	2000-2009 –> 2010-2019	10.631038	41384.020	41469.621
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	19048.1657	19090.7771	4.5117905	2000-2009 –> 2010-2019	8.786451	35384.586	35468.863
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	22229.3384	22272.4540	5.8026365	2000-2009 –> 2010-2019	11.235598	42125.331	42210.932
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ sal + temp + silicate + phosphate_star	19255.2157	19291.7397	4.6591170	2000-2009 –> 2010-2019	9.154677	35876.027	35948.265
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + phosphate	22189.0206	22225.9768	5.7709236	2000-2009 –> 2010-2019	11.239215	42124.432	42197.804
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	19189.2277	19231.8390	4.6106825	2000-2009 –> 2010-2019	9.045294	35734.452	35818.729
obs	Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	22189.5519	22232.6675	5.7697116	2000-2009 –> 2010-2019	11.226615	42113.642	42199.243
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	19358.5877	19401.1990	4.7322807	2000-2009 –> 2010-2019	9.276688	36042.827	36127.104
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + aou + silicate + phosphate	19390.3469	19426.8709	4.7569000	2000-2009 –> 2010-2019	9.364835	36151.493	36223.731
mod	Indo-Pacific	(26.5,26.75]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	19381.9271	19424.5384	4.7492875	2000-2009 –> 2010-2019	9.312518	36089.661	36173.938
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ aou + nitrate + silicate + phosphate	33121.5483	33161.0787	5.2828271	2000-2009 –> 2010-2019	10.683907	63232.417	63310.863
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	33072.8145	33118.9333	5.2579262	2000-2009 –> 2010-2019	10.615170	63106.701	63198.221
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	24214.1902	24258.8598	3.8695323	2000-2009 –> 2010-2019	7.485539	45960.016	46048.775
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	32988.1175	33034.2363	5.2166169	2000-2009 –> 2010-2019	10.430724	62759.527	62851.047
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	24484.6792	24529.3488	3.9913020	2000-2009 –> 2010-2019	7.689270	46410.616	46499.376
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	33248.0333	33294.1521	5.3444269	2000-2009 –> 2010-2019	10.671347	63226.904	63318.424
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	24498.4197	24543.0893	3.9975891	2000-2009 –> 2010-2019	7.701554	46437.380	46526.139
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + nitrate + phosphate_star	24605.0919	24643.3801	4.0476626	2000-2009 –> 2010-2019	7.808444	46664.383	46740.462
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	24473.0355	24517.7051	3.9859822	2000-2009 –> 2010-2019	7.727466	46492.985	46581.744
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	32942.1088	32988.2276	5.1943133	2000-2009 –> 2010-2019	10.391223	62681.500	62773.020
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + nitrate + phosphate	25075.5067	25113.7949	4.2717536	2000-2009 –> 2010-2019	8.208455	47502.174	47578.254
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + nitrate + phosphate	33032.5499	33072.0803	5.2392230	2000-2009 –> 2010-2019	10.503332	62894.460	62972.906
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	24729.4548	24774.1244	4.1047956	2000-2009 –> 2010-2019	7.957511	46984.827	47073.587
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	32821.2850	32867.4038	5.1361946	2000-2009 –> 2010-2019	10.308978	62515.576	62607.096
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	24967.7179	25012.3875	4.2183683	2000-2009 –> 2010-2019	8.174127	47435.194	47523.953
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	32998.8137	33044.9325	5.2218158	2000-2009 –> 2010-2019	10.469368	62832.192	62923.712
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	24980.5937	25025.2633	4.2245945	2000-2009 –> 2010-2019	8.188932	47465.477	47554.237
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	32996.3125	33042.4313	5.2205996	2000-2009 –> 2010-2019	10.470538	62834.097	62925.617
mod	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	25065.9718	25104.2600	4.2670906	2000-2009 –> 2010-2019	8.296263	47679.219	47755.298
obs	Indo-Pacific	(26.75,27]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	32996.8248	33036.3552	5.2218212	2000-2009 –> 2010-2019	10.472493	62833.964	62912.410
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	40552.8711	40600.6181	4.8204152	2000-2009 –> 2010-2019	9.612730	79748.534	79843.803
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	40748.9845	40796.7315	4.8906900	2000-2009 –> 2010-2019	9.796344	80251.372	80346.642
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + aou + silicate + phosphate	40551.1189	40592.0449	4.8205034	2000-2009 –> 2010-2019	9.612893	79744.986	79826.646
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + aou + silicate + phosphate_star	40776.9682	40817.8942	4.9015241	2000-2009 –> 2010-2019	9.818009	80306.296	80387.956
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	40672.7357	40720.4827	4.8632463	2000-2009 –> 2010-2019	9.727825	80064.791	80160.061
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	40697.1499	40738.0759	4.8727358	2000-2009 –> 2010-2019	9.751485	80125.373	80207.032
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	28609.0067	28655.5345	2.9829741	2000-2009 –> 2010-2019	5.909444	54652.149	54744.600
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	41172.7001	41220.4471	5.0460407	2000-2009 –> 2010-2019	10.010355	80831.068	80926.337
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	29032.0206	29078.5483	3.0959024	2000-2009 –> 2010-2019	6.111191	55387.366	55479.817
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	40710.8632	40758.6102	4.8769499	2000-2009 –> 2010-2019	9.721167	80047.879	80143.148
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	29056.1896	29102.7173	3.1024822	2000-2009 –> 2010-2019	6.125766	55439.184	55531.635
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	29174.7521	29221.2798	3.1349629	2000-2009 –> 2010-2019	6.199525	55699.354	55791.805
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	40895.3789	40943.1259	4.9438155	2000-2009 –> 2010-2019	9.811196	80294.989	80390.258
obs	Indo-Pacific	(27,27.25]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	40673.5745	40721.3215	4.8635473	2000-2009 –> 2010-2019	9.697596	79983.015	80078.284
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ temp + aou + nitrate + phosphate	28918.3062	28958.1871	3.0656700	2000-2009 –> 2010-2019	6.079818	55267.699	55346.943
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ temp + aou + nitrate + phosphate_star	29283.4431	29323.3240	3.1655940	2000-2009 –> 2010-2019	6.255134	55890.807	55970.050
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	28787.3003	28833.8281	3.0300605	2000-2009 –> 2010-2019	5.968352	54872.540	54964.991
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	29227.3698	29273.8976	3.1494865	2000-2009 –> 2010-2019	6.194395	55684.788	55777.239
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ temp + nitrate + phosphate + phosphate_star	29299.5135	29339.3944	3.1700660	2000-2009 –> 2010-2019	6.265718	55927.514	56006.757
mod	Indo-Pacific	(27,27.25]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	29238.9616	29285.4893	3.1526951	2000-2009 –> 2010-2019	6.202570	55713.398	55805.849
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	38346.4038	38394.7267	3.2758620	2000-2009 –> 2010-2019	6.260960	72094.185	72190.190
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	42992.7627	43041.2743	4.1566134	2000-2009 –> 2010-2019	8.458990	82874.288	82970.701
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	43480.0769	43528.5885	4.2928173	2000-2009 –> 2010-2019	8.767103	83904.311	84000.724
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + aou + silicate + phosphate	42992.4492	43034.0306	4.1570772	2000-2009 –> 2010-2019	8.469444	82904.103	82986.743
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + aou + silicate + phosphate_star	43491.0968	43532.6781	4.2965169	2000-2009 –> 2010-2019	8.794391	83986.168	84068.808
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	38876.5658	38924.8887	3.3960638	2000-2009 –> 2010-2019	6.546779	73349.309	73445.315
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	43260.5386	43309.0502	4.2309127	2000-2009 –> 2010-2019	8.616782	83408.305	83504.718
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	43293.2388	43334.8202	4.2406376	2000-2009 –> 2010-2019	8.666334	83564.222	83646.862
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	38832.4995	38880.8224	3.3859069	2000-2009 –> 2010-2019	6.460299	72976.008	73072.013
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate	42925.6199	42974.1315	4.1381889	2000-2009 –> 2010-2019	8.407784	82701.198	82797.611
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	43085.8475	43134.3591	4.1822923	2000-2009 –> 2010-2019	8.494452	82998.836	83095.249
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	38847.9980	38896.3209	3.3894757	2000-2009 –> 2010-2019	6.469830	73017.516	73113.521
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	38238.9792	38287.3021	3.2520293	2000-2009 –> 2010-2019	6.149032	71584.570	71680.576
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	43106.9992	43155.5108	4.1881494	2000-2009 –> 2010-2019	8.510053	83051.252	83147.664
obs	Indo-Pacific	(27.25,27.5]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	43051.6805	43100.1921	4.1728484	2000-2009 –> 2010-2019	8.476987	82938.879	83035.292
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + aou + nitrate + phosphate_star	39041.2197	39082.6393	3.4347523	2000-2009 –> 2010-2019	6.561221	73408.242	73490.533
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	38599.7370	38648.0599	3.3327591	2000-2009 –> 2010-2019	6.357117	72522.945	72618.950
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + nitrate + phosphate + phosphate_star	39063.6018	39105.0214	3.4399817	2000-2009 –> 2010-2019	6.574390	73464.676	73546.966
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	38595.0824	38643.4053	3.3317049	2000-2009 –> 2010-2019	6.358008	72526.921	72622.927
mod	Indo-Pacific	(27.25,27.5]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	38603.8126	38645.2322	3.3341357	2000-2009 –> 2010-2019	6.365115	72554.380	72636.671
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	42308.1886	42357.3288	3.1238761	2000-2009 –> 2010-2019	5.964878	80455.830	80553.650
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	47639.3504	47688.7654	3.8598228	2000-2009 –> 2010-2019	7.825872	90907.253	91005.340
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	42353.2772	42402.4174	3.1324066	2000-2009 –> 2010-2019	5.976106	80515.611	80613.431
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + silicate + phosphate	47723.8686	47766.2243	3.8792916	2000-2009 –> 2010-2019	7.881480	91130.058	91214.133
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + aou + silicate + phosphate_star	43031.8626	43073.9828	3.2640357	2000-2009 –> 2010-2019	6.222655	81805.547	81889.393
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	42324.0324	42373.1726	3.1268710	2000-2009 –> 2010-2019	5.965730	80459.990	80557.811
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	43106.0046	43148.1248	3.2787052	2000-2009 –> 2010-2019	6.248798	81939.611	82023.457
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	48372.3044	48421.7194	4.0278986	2000-2009 –> 2010-2019	8.031486	91785.900	91883.988
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	48403.9553	48453.3703	4.0353191	2000-2009 –> 2010-2019	8.048285	91853.738	91951.826
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate	47197.2492	47246.6642	3.7618533	2000-2009 –> 2010-2019	7.578263	89870.324	89968.412
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	42838.5570	42887.6972	3.2257065	2000-2009 –> 2010-2019	6.162816	81501.271	81599.091
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	47250.7909	47300.2059	3.7735845	2000-2009 –> 2010-2019	7.602180	89973.169	90071.257
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate	42580.6974	42629.8376	3.1757896	2000-2009 –> 2010-2019	6.099952	81175.012	81272.833
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	42760.8188	42809.9590	3.2105758	2000-2009 –> 2010-2019	6.115323	81252.000	81349.820
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	47263.4082	47312.8232	3.7763543	2000-2009 –> 2010-2019	7.606430	89991.763	90089.850
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + silicate + phosphate	42798.0632	42840.1833	3.2182053	2000-2009 –> 2010-2019	6.168804	81529.718	81613.564
mod	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	42795.7355	42844.8757	3.2173631	2000-2009 –> 2010-2019	6.146913	81418.548	81516.368
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	47243.9182	47293.3332	3.7720766	2000-2009 –> 2010-2019	7.599384	89961.091	90059.179
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ sal + temp + silicate + phosphate_star	47541.9235	47584.2792	3.8384625	2000-2009 –> 2010-2019	7.698913	90390.452	90474.527
obs	Indo-Pacific	(27.5,27.75]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	48428.0989	48477.5139	4.0409888	2000-2009 –> 2010-2019	8.053157	91874.809	91972.896
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	15008.5248	15051.0189	2.5210824	2000-2009 –> 2010-2019	4.812698	28480.223	28564.745
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + aou + nitrate	14990.3443	15026.7678	2.5147147	2000-2009 –> 2010-2019	4.811002	28472.235	28544.683
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + aou + nitrate + phosphate_star	14726.9406	14769.4347	2.4125326	2000-2009 –> 2010-2019	4.606783	27938.747	28023.269
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + aou + nitrate + silicate	14980.3666	15022.8608	2.5100113	2000-2009 –> 2010-2019	4.797667	28441.713	28526.235
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + aou + phosphate_star	14921.0932	14957.5167	2.4876425	2000-2009 –> 2010-2019	4.739486	28285.990	28358.437
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	14890.4862	14932.9803	2.4749967	2000-2009 –> 2010-2019	4.725437	28253.652	28338.174
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	24207.5352	24252.4908	3.4606132	2000-2009 –> 2010-2019	6.971167	45731.571	45820.621
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star	14743.9867	14786.4808	2.4189688	2000-2009 –> 2010-2019	4.619201	27972.091	28056.614
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	15019.0288	15061.5229	2.5252249	2000-2009 –> 2010-2019	4.819420	28497.463	28581.985
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	24195.8463	24240.8018	3.4561680	2000-2009 –> 2010-2019	6.959797	45704.004	45793.055
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + phosphate + phosphate_star	14925.2611	14961.6846	2.4892636	2000-2009 –> 2010-2019	4.743232	28295.803	28368.250
mod	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	14898.3411	14940.8352	2.4780371	2000-2009 –> 2010-2019	4.730981	28268.163	28352.686
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	24205.3607	24250.3163	3.4597859	2000-2009 –> 2010-2019	6.969583	45727.664	45816.715
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ sal + temp + silicate + phosphate_star	24227.8668	24266.4002	3.4691217	2000-2009 –> 2010-2019	6.989424	45772.204	45848.533
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	24252.8963	24297.8519	3.4779180	2000-2009 –> 2010-2019	7.001173	45805.977	45895.028
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	24206.7320	24251.6875	3.4603076	2000-2009 –> 2010-2019	6.969010	45726.526	45815.577
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + aou + silicate + phosphate_star	24223.1629	24261.6962	3.4673277	2000-2009 –> 2010-2019	6.986567	45765.072	45841.401
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	24207.4555	24252.4111	3.4605829	2000-2009 –> 2010-2019	6.969523	45727.794	45816.845
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	24205.8627	24244.3960	3.4607379	2000-2009 –> 2010-2019	6.969978	45724.890	45801.219
obs	Indo-Pacific	(27.75,27.85]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	24220.2916	24258.8250	3.4662332	2000-2009 –> 2010-2019	6.984460	45759.886	45836.215
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou	17446.2537	17477.6951	2.1666842	2000-2009 –> 2010-2019	4.215707	33234.977	33297.489
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou + phosphate_star	17299.1343	17336.8640	2.1264422	2000-2009 –> 2010-2019	4.109251	32847.238	32922.252
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate	17448.1673	17485.8970	2.1666606	2000-2009 –> 2010-2019	4.210206	33219.120	33294.134
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate + phosphate	33496.3446	33543.8004	3.1804501	2000-2009 –> 2010-2019	6.333720	62721.978	62815.957
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	17301.0670	17345.0850	2.1264242	2000-2009 –> 2010-2019	4.102839	32827.314	32914.831
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	33497.3445	33544.8003	3.1806948	2000-2009 –> 2010-2019	6.340970	62748.227	62842.206
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + nitrate + silicate	17593.1740	17630.9037	2.2065225	2000-2009 –> 2010-2019	4.280765	33474.251	33549.265
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	17485.9711	17529.9890	2.1764355	2000-2009 –> 2010-2019	4.198333	33180.262	33267.780
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	33493.7833	33541.2391	3.1798235	2000-2009 –> 2010-2019	6.333461	62720.743	62814.723
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + phosphate + phosphate_star	17306.8084	17344.5381	2.1284948	2000-2009 –> 2010-2019	4.113707	32863.858	32938.873
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + silicate + phosphate	17566.4766	17604.2063	2.1991288	2000-2009 –> 2010-2019	4.273883	33449.376	33524.391
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	17308.8059	17352.8239	2.1284941	2000-2009 –> 2010-2019	4.106814	32842.169	32929.686
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	33498.3225	33545.7783	3.1809341	2000-2009 –> 2010-2019	6.341849	62751.507	62845.486
mod	Indo-Pacific	(27.85,27.95]	cstar_tref ~ sal + temp + silicate + phosphate_star	17488.0473	17525.7770	2.1775512	2000-2009 –> 2010-2019	4.200042	33182.509	33257.524
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	33486.4402	33533.8961	3.1780276	2000-2009 –> 2010-2019	6.322901	62681.734	62775.713
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	33493.3313	33540.7871	3.1797129	2000-2009 –> 2010-2019	6.333884	62722.215	62816.195
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + aou + silicate + phosphate	33496.5274	33537.2038	3.1809843	2000-2009 –> 2010-2019	6.335399	62724.289	62804.843
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + aou + silicate + phosphate_star	33500.6258	33541.3022	3.1819874	2000-2009 –> 2010-2019	6.344148	62756.295	62836.849
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	33494.0045	33541.4603	3.1798776	2000-2009 –> 2010-2019	6.334007	62722.737	62816.717
obs	Indo-Pacific	(27.85,27.95]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	33495.1409	33535.8173	3.1806450	2000-2009 –> 2010-2019	6.334816	62722.025	62802.579
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + aou	21157.1576	21189.8631	1.9075321	2000-2009 –> 2010-2019	3.751221	39264.816	39329.532
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + aou + phosphate_star	20905.0869	20944.3335	1.8607945	2000-2009 –> 2010-2019	3.578573	38384.345	38462.004
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate	21093.3721	21132.6187	1.8953191	2000-2009 –> 2010-2019	3.732828	39173.109	39250.768
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate + phosphate	56846.9333	56898.1621	3.1008687	2000-2009 –> 2010-2019	6.090162	104988.484	105089.879
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	20730.2709	20776.0586	1.8289458	2000-2009 –> 2010-2019	3.513551	38037.735	38128.337
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	56996.5658	57047.7945	3.1217623	2000-2009 –> 2010-2019	6.147833	105372.206	105473.601
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + nitrate + phosphate_star	21206.3004	21245.5470	1.9163325	2000-2009 –> 2010-2019	3.703702	39039.482	39117.141
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + nitrate + silicate	21217.7930	21257.0396	1.9184840	2000-2009 –> 2010-2019	3.810217	39556.716	39634.375
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	20826.1765	20871.9642	1.8461524	2000-2009 –> 2010-2019	3.562164	38298.259	38388.861
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	56949.6442	57000.8729	3.1151954	2000-2009 –> 2010-2019	6.131539	105263.647	105365.042
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + phosphate + phosphate_star	20910.6126	20949.8592	1.8617987	2000-2009 –> 2010-2019	3.580880	38396.624	38474.283
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	20734.1003	20779.8881	1.8296298	2000-2009 –> 2010-2019	3.515326	38047.337	38137.938
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	56995.7674	57046.9961	3.1216505	2000-2009 –> 2010-2019	6.147066	105367.264	105468.659
mod	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + silicate + phosphate_star	20835.1672	20874.4138	1.8481346	2000-2009 –> 2010-2019	3.578398	38378.963	38456.622
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ sal + temp + silicate + phosphate_star	56994.5664	57038.4767	3.1217624	2000-2009 –> 2010-2019	6.147984	105369.161	105456.071
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate	56875.8848	56927.1135	3.1049003	2000-2009 –> 2010-2019	6.094785	105021.217	105122.612
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	57026.2087	57077.4374	3.1259181	2000-2009 –> 2010-2019	6.149358	105385.198	105486.593
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + aou + silicate + phosphate	56988.7222	57032.6326	3.1209437	2000-2009 –> 2010-2019	6.135458	105289.112	105376.022
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	57017.4883	57068.7170	3.1246950	2000-2009 –> 2010-2019	6.145463	105359.554	105460.948
obs	Indo-Pacific	(27.95,28.05]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	57033.3318	57077.2421	3.1271983	2000-2009 –> 2010-2019	6.154509	105414.814	105501.724
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	48481.7949	48532.0606	2.9361098	2000-2009 –> 2010-2019	5.802653	86594.960	86693.874
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	48756.7531	48807.0187	2.9779808	2000-2009 –> 2010-2019	5.902011	87175.939	87274.853
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + aou + silicate + phosphate	48516.4802	48559.5650	2.9416621	2000-2009 –> 2010-2019	5.820926	86695.886	86780.669
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	48658.3669	48708.6326	2.9629302	2000-2009 –> 2010-2019	5.865533	86964.194	87063.108
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	48733.7132	48776.7981	2.9747558	2000-2009 –> 2010-2019	5.902074	87168.219	87253.002
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + aou + silicate + phosphate	48499.5945	48549.8602	2.9388024	2000-2009 –> 2010-2019	5.810209	86638.884	86737.797
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	10975.6533	11017.5517	1.5629222	2000-2009 –> 2010-2019	3.097909	20622.311	20705.272
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	48661.9859	48712.2515	2.9634825	2000-2009 –> 2010-2019	5.863923	86956.329	87055.243
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + nitrate + silicate	11339.1846	11375.0975	1.6632355	2000-2009 –> 2010-2019	3.387095	21588.897	21660.005
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	10961.8927	11003.7911	1.5592664	2000-2009 –> 2010-2019	3.093625	20606.418	20689.378
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	48827.7163	48877.9819	2.9888837	2000-2009 –> 2010-2019	5.915471	87260.373	87359.287
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + silicate + phosphate	11307.9674	11343.8803	1.6544227	2000-2009 –> 2010-2019	3.355322	21487.767	21558.876
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	10974.7857	11016.6840	1.5626914	2000-2009 –> 2010-2019	3.097496	20620.824	20703.784
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	48632.1354	48682.4010	2.9589304	2000-2009 –> 2010-2019	5.854673	86901.495	87000.409
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + silicate + phosphate_star	10995.9417	11031.8546	1.5688618	2000-2009 –> 2010-2019	3.106269	20648.814	20719.923
obs	Indo-Pacific	(28.05,28.1]	cstar_tref ~ sal + temp + silicate + phosphate_star	48881.1678	48924.2526	2.9974312	2000-2009 –> 2010-2019	5.925420	87319.203	87403.987
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + aou + silicate + phosphate_star	11348.1262	11384.0391	1.6657683	2000-2009 –> 2010-2019	3.293387	21298.306	21369.415
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + nitrate + silicate + phosphate_star	11351.9818	11387.8947	1.6668617	2000-2009 –> 2010-2019	3.297550	21311.984	21383.093
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + silicate + phosphate + phosphate_star	11345.9440	11381.8569	1.6651498	2000-2009 –> 2010-2019	3.292196	21294.288	21365.396
mod	Indo-Pacific	(28.05,28.1]	cstar_tref ~ temp + silicate + phosphate_star	11373.7388	11403.6663	1.6736146	2000-2009 –> 2010-2019	3.306770	21339.626	21398.883
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate	35560.7170	35608.9634	2.7837502	2000-2009 –> 2010-2019	5.488958	67639.022	67734.876
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	35711.8067	35760.0530	2.8128038	2000-2009 –> 2010-2019	5.573254	68058.609	68154.463
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + aou + silicate + phosphate	24673.8154	24714.8745	1.4350273	2000-2009 –> 2010-2019	2.936110	47047.250	47128.631
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + aou + silicate + phosphate	35578.0094	35619.3634	2.7874433	2000-2009 –> 2010-2019	5.499311	67686.972	67769.132
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + aou + silicate + phosphate_star	35728.1359	35769.4899	2.8163489	2000-2009 –> 2010-2019	5.582442	68100.062	68182.222
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	35650.5548	35698.8012	2.8009891	2000-2009 –> 2010-2019	5.539478	67891.267	67987.121
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	35680.8393	35722.1934	2.8072102	2000-2009 –> 2010-2019	5.555815	67968.526	68050.686
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + aou + nitrate + silicate	23386.3418	23434.2441	1.3074883	2000-2009 –> 2010-2019	2.765539	45405.409	45500.353
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + aou + silicate	24362.8753	24403.9344	1.4031782	2000-2009 –> 2010-2019	2.955240	47145.497	47226.878
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + aou + silicate + phosphate	22897.6812	22945.5835	1.2621742	2000-2009 –> 2010-2019	2.605990	43917.143	44012.087
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + aou + silicate + phosphate	35572.7929	35621.0392	2.7860613	2000-2009 –> 2010-2019	5.495450	67671.612	67767.466
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	22912.4166	22960.3189	1.2635174	2000-2009 –> 2010-2019	2.608746	43944.753	44039.697
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	35635.2020	35683.4484	2.7980355	2000-2009 –> 2010-2019	5.538574	67885.851	67981.704
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + nitrate + silicate	24030.5242	24071.5833	1.3699171	2000-2009 –> 2010-2019	2.896898	46613.538	46694.919
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	22949.1972	22997.0994	1.2668763	2000-2009 –> 2010-2019	2.622839	44078.909	44173.853
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	35752.9811	35801.2275	2.8207738	2000-2009 –> 2010-2019	5.601592	68197.425	68293.278
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + silicate + phosphate	23846.3776	23887.4367	1.3518287	2000-2009 –> 2010-2019	2.838996	46105.704	46187.085
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	22909.3788	22957.2810	1.2632403	2000-2009 –> 2010-2019	2.608231	43939.546	44034.491
obs	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	35618.2145	35666.4609	2.7947711	2000-2009 –> 2010-2019	5.527802	67832.426	67928.279
mod	Indo-Pacific	(28.1,28.2]	cstar_tref ~ sal + temp + silicate + phosphate_star	23099.2871	23140.3462	1.2808608	2000-2009 –> 2010-2019	2.641533	44269.480	44350.861
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star	12047.5849	12088.8243	2.2960054	2000-2009 –> 2010-2019	4.731846	22360.752	22441.959
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ aou + silicate + phosphate + phosphate_star	12052.9786	12088.3266	2.2991819	2000-2009 –> 2010-2019	4.737846	22369.313	22438.919
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star	12036.2459	12077.4852	2.2911425	2000-2009 –> 2010-2019	4.729617	22354.233	22435.441
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + aou + phosphate + phosphate_star	12037.6160	12072.9640	2.2925868	2000-2009 –> 2010-2019	4.733294	22357.685	22427.292
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + aou + silicate + phosphate	44025.9556	44072.6571	0.8365994	2000-2009 –> 2010-2019	1.668898	79532.750	79624.887
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star	12033.2108	12074.4501	2.2898426	2000-2009 –> 2010-2019	4.726283	22347.476	22428.684
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + aou + silicate + phosphate_star	45600.1212	45646.8226	0.8745531	2000-2009 –> 2010-2019	1.756114	82759.056	82851.193
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + nitrate + silicate + phosphate_star	46264.8661	46311.5675	0.8910930	2000-2009 –> 2010-2019	1.801635	84353.170	84445.307
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + silicate + phosphate + phosphate_star	45156.7706	45203.4720	0.8636929	2000-2009 –> 2010-2019	1.731355	81859.108	81951.246
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	12039.2258	12080.4651	2.2924195	2000-2009 –> 2010-2019	4.734175	22363.211	22444.418
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + temp + aou + nitrate + silicate	45987.8884	46042.3735	0.8841139	2000-2009 –> 2010-2019	1.815500	84728.518	84836.011
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + temp + aou + phosphate + phosphate_star	12036.1970	12077.4363	2.2911215	2000-2009 –> 2010-2019	4.729422	22353.867	22435.075
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + temp + aou + silicate + phosphate	42890.6116	42945.0967	0.8102067	2000-2009 –> 2010-2019	1.616963	77503.819	77611.312
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + temp + aou + silicate + phosphate	12046.9967	12088.2360	2.2957529	2000-2009 –> 2010-2019	4.740574	22376.577	22457.785
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	42911.9400	42966.4250	0.8106939	2000-2009 –> 2010-2019	1.617687	77533.594	77641.087
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	43301.9898	43356.4749	0.8196555	2000-2009 –> 2010-2019	1.640883	78426.033	78533.527
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	42909.2393	42963.7243	0.8106322	2000-2009 –> 2010-2019	1.617596	77529.832	77637.326
mod	Indo-Pacific	(28.2,28.3]	cstar_tref ~ sal + temp + silicate + phosphate_star	43437.9075	43484.6090	0.8228478	2000-2009 –> 2010-2019	1.648202	78703.938	78796.075
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star	12051.4033	12092.6426	2.2976453	2000-2009 –> 2010-2019	4.735933	22369.049	22450.257
obs	Indo-Pacific	(28.2,28.3]	cstar_tref ~ temp + aou + phosphate + phosphate_star	12055.2545	12090.6024	2.3001605	2000-2009 –> 2010-2019	4.740363	22374.402	22444.008
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ aou + nitrate + phosphate + phosphate_star	22070.1073	22111.1759	1.1861590	2000-2009 –> 2010-2019	2.285019	40437.369	40518.699
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ aou + nitrate + silicate + phosphate_star	22387.0144	22428.0830	1.2135609	2000-2009 –> 2010-2019	2.330143	40948.337	41029.668
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ aou + nitrate + silicate + phosphate_star	4211.9549	4240.6887	2.5751323	2000-2009 –> 2010-2019	5.156003	7722.477	7778.843
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ aou + silicate + phosphate_star	4222.9300	4246.8748	2.5940146	2000-2009 –> 2010-2019	5.178712	7733.643	7780.614
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + aou + nitrate + phosphate	22588.9020	22629.9707	1.2313465	2000-2009 –> 2010-2019	2.385118	41547.654	41628.984
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + aou + nitrate + phosphate_star	22751.4175	22792.4861	1.2458528	2000-2009 –> 2010-2019	2.379631	41498.131	41579.461
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star	4213.5141	4247.0369	2.5744933	2000-2009 –> 2010-2019	5.143587	7719.277	7785.037
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + aou + silicate + phosphate_star	4224.8911	4253.6249	2.5939579	2000-2009 –> 2010-2019	5.166261	7730.499	7786.865
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star	4219.1847	4252.7075	2.5827265	2000-2009 –> 2010-2019	5.168789	7734.678	7800.438
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + aou + nitrate + phosphate	21822.6582	21870.5716	1.1650260	2000-2009 –> 2010-2019	2.243428	39963.955	40058.841
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + aou + silicate + phosphate_star	4224.2493	4257.7721	2.5901022	2000-2009 –> 2010-2019	5.161366	7731.260	7797.020
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + nitrate + phosphate	21823.0587	21864.1274	1.1652275	2000-2009 –> 2010-2019	2.246200	39991.233	40072.563
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star	22643.9069	22691.8203	1.2360591	2000-2009 –> 2010-2019	2.387352	41578.575	41673.461
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star	4225.8411	4259.3639	2.5924246	2000-2009 –> 2010-2019	5.177773	7740.926	7806.686
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ temp + aou + nitrate + phosphate	22046.5002	22087.5688	1.1841428	2000-2009 –> 2010-2019	2.280422	40385.235	40466.565
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ temp + aou + nitrate + phosphate_star	22091.7382	22132.8068	1.1880096	2000-2009 –> 2010-2019	2.288235	40474.062	40555.392
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star	4212.1391	4245.6619	2.5725009	2000-2009 –> 2010-2019	5.150599	7723.073	7788.833
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ temp + aou + silicate + phosphate_star	4222.3771	4251.1109	2.5902886	2000-2009 –> 2010-2019	5.172416	7733.619	7789.985
mod	Indo-Pacific	(28.3,28.4]	cstar_tref ~ temp + nitrate + phosphate + phosphate_star	22051.3155	22092.3841	1.1845538	2000-2009 –> 2010-2019	2.281379	40396.094	40477.424
obs	Indo-Pacific	(28.3,28.4]	cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star	4216.3666	4249.8894	2.5786316	2000-2009 –> 2010-2019	5.169820	7734.786	7800.546

6.3 Target variable coefficients

A data frame to map the target variable is prepared.

# create table with two era belonging to one eras
eras_forward <- GLODAP_glanced %>%
  arrange(era) %>% 
  group_by(basin, gamma_slab, model) %>% 
  mutate(eras = paste(era, lead(era), sep = " --> ")) %>% 
  ungroup() %>% 
  select(era, eras) %>% 
  unique()

eras_backward <- GLODAP_glanced %>%
  arrange(era) %>% 
  group_by(basin, gamma_slab, model) %>% 
  mutate(eras = paste(lag(era), era, sep = " --> ")) %>% 
  ungroup() %>% 
  select(era, eras) %>% 
  unique()

eras_era <- full_join(eras_backward, eras_forward) %>% 
  filter(str_detect(eras, "NA") == FALSE)

# extend best model selection from eras to era
lm_best <- full_join(
  lm_best_eras %>% select(data_source, basin, gamma_slab, model, eras),
  eras_era) %>% 
  drop_na() %>% 
  arrange(era, basin, gamma_slab, data_source)

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

rm(eras_era, eras_forward, eras_backward)

6.4 Cant coeffcients

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

# subtract coefficients of adjacent era  
lm_best_cant <- lm_best %>%
  arrange(era) %>%
  group_by(data_source, basin, gamma_slab, eras, model, term) %>%
  mutate(delta_coeff = estimate - lag(estimate)) %>%
  ungroup() %>%
  arrange(basin, gamma_slab, model, term, eras) %>%
  drop_na() %>%
  select(-c(era,estimate))

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

6.5 Write files

# create table of target varaible coefficients in wide format
lm_best_target <- lm_best %>% 
  pivot_wider(names_from = "term",
              names_prefix = "coeff_",
              values_from = "estimate",
              values_fill = 0
              )


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

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

7 Model diagnotics

7.1 Selection criterion vs predictors

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

7.1.1 All models

GLODAP_glanced %>%
  group_split(era) %>%
  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
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
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

7.1.2 Best models

left_join(lm_best_target %>% select(basin, gamma_slab, era, model),
          GLODAP_glanced) %>%
  group_split(era) %>%
  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
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
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

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

7.2.1 All models

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

lm_all_eras %>%
  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
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)

7.2.2 Best models

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

lm_best_eras %>%
  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
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)

7.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_cant %>% 
  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	10	7	5	10	5
(26,26.5]	7	7	6	7	7	7	9
(26.5,26.75]	6	4	7	7	4	9	10
(26.75,27]	8	4	8	5	5	10	7
(27,27.25]	8	6	4	4	7	10	8
(27.25,27.5]	5	4	4	6	10	8	8
(27.5,27.75]	6	1	6	6	10	6	8
(27.75,27.85]	5	4	4	6	10	10	6
(27.85,27.95]	6	2	5	6	10	7	9
(27.95,28.05]	4	2	4	8	10	5	10
(28.05,28.1]	9	1	4	5	6	6	10
(28.1,28.15]	10	2	4	4	6	5	10
(28.15,28.2]	5	7	2	6	6	5	10
(28.2, Inf]	5	10	NA	5	8	6	6
total	89.00	59.00	—	82.00	104.00	104.00	116.00
Indo-Pacific - mod
(-Inf,26]	9	6	6	6	10	4	7
(26,26.5]	7	6	7	6	6	7	10
(26.5,26.75]	6	5	6	6	6	9	10
(26.75,27]	5	10	5	7	5	5	10
(27,27.25]	6	10	6	7	4	4	10
(27.25,27.5]	4	10	5	9	5	6	8
(27.5,27.75]	4	5	6	7	10	10	5
(27.75,27.85]	6	6	3	8	10	5	9
(27.85,27.95]	4	2	3	6	10	7	10
(27.95,28.05]	4	3	2	7	10	6	10
(28.05,28.1]	2	3	3	8	6	10	10
(28.1,28.2]	5	3	4	4	10	10	9
(28.2,28.3]	5	3	4	7	10	10	6
(28.3,28.4]	7	10	6	6	5	2	6
total	74.00	82.00	66.00	94.00	107.00	95.00	120.00
Atlantic - obs
(-Inf,26]	8	10	6	6	6	5	8
(26,26.5]	7	9	7	6	4	4	10
(26.5,26.75]	8	8	6	6	4	5	10
(26.75,27]	4	4	4	9	5	10	10
(27,27.25]	5	6	5	9	5	8	10
(27.25,27.5]	6	4	7	7	9	4	10
(27.5,27.75]	7	4	7	6	4	10	9
(27.75,27.85]	6	5	8	6	7	10	6
(27.85,27.95]	5	4	7	6	10	10	5
(27.95,28.05]	5	5	5	8	10	10	4
(28.05,28.1]	5	5	3	10	10	6	7
(28.1,28.15]	5	4	2	10	7	7	10
(28.15,28.2]	7	1	7	6	6	6	10
(28.2, Inf]	6	4	7	7	10	4	9
total	84.00	73.00	81.00	102.00	97.00	99.00	118.00
Indo-Pacific - obs
(-Inf,26]	8	6	6	6	10	5	8
(26,26.5]	5	5	5	7	10	5	10
(26.5,26.75]	6	7	7	5	8	4	10
(26.75,27]	7	10	6	5	4	7	8
(27,27.25]	6	5	6	7	10	9	4
(27.25,27.5]	6	4	6	7	10	10	4
(27.5,27.75]	7	5	6	5	9	8	8
(27.75,27.85]	4	5	4	9	4	10	10
(27.85,27.95]	6	5	5	7	4	10	10
(27.95,28.05]	5	5	5	7	5	10	10
(28.05,28.1]	5	4	6	7	10	10	5
(28.1,28.2]	6	4	6	7	10	10	4
(28.2,28.3]	10	4	10	8	6	4	5
(28.3,28.4]	7	5	3	10	5	10	5
total	88.00	74.00	81.00	97.00	105.00	112.00	101.00

7.4 RMSE alternatives

7.4.1 AIC

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

lm_all_eras %>% 
  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
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_eras %>% 
  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
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

7.4.2 AIC vs BIC

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

lm_all_eras %>% 
  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
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

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

7.6 Residual patterns

7.6.1 Fitted vs actual models

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

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

# calculate equal axis limits and binwidth
axis_lims <- GLODAP_augmented %>%
  summarise(
    max_value = max(
      c(max(.fitted, max(!!sym(params_local$MLR_target))))
      ),
    min_value = min(
      c(min(.fitted, min(!!sym(params_local$MLR_target))))
    )
  )

i_binwidth <- 1
# binwidth_value <- (axis_lims$max_value - axis_lims$min_value) / 40
axis_lims <- c(axis_lims$min_value, axis_lims$max_value)

GLODAP_augmented %>%
  ggplot(aes(cstar_tref, .fitted)) +
  geom_bin2d(binwidth = i_binwidth) +
  scale_fill_viridis_c() +
  geom_abline(slope = 1,
              col = 'red') +
  coord_equal(xlim = axis_lims,
              ylim = axis_lims) +
  labs(title = "All models") +
  facet_wrap(~data_source)

Past versions of fitted_vs_actual_all-1.png

Version	Author	Date
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
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

GLODAP_augmented_best %>%
  ggplot(aes(cstar_tref, .fitted)) +
  geom_bin2d(binwidth = i_binwidth) +
  scale_fill_viridis_c() +
  geom_abline(slope = 1,
              col = 'red') +
  coord_equal(xlim = axis_lims,
              ylim = axis_lims) +
  labs(title = "Selected models") +
  facet_wrap(~data_source)

Past versions of fitted_vs_actual_best-1.png

Version	Author	Date
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

rm(binwidth_value, axis_lims)

7.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 = "white") +
  scale_fill_viridis_c() +
  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
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 = "white") +
  scale_fill_viridis_c() +
  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
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

7.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, era) %>% 
  map(~ ggplot(data = .x, aes(lat, .resid)) +
      geom_bin2d(binwidth = i_binwidth) +
      geom_hline(yintercept = 0, col = "black") +
      scale_fill_viridis_c() +
      facet_grid(gamma_slab ~ data_source) +
      coord_cartesian(ylim = i_ylim) +
      labs(
        title = paste(
          "Target:",
          params_local$MLR_target,
          "| best models | basin:",
          unique(.x$basin),
          "| era:",
          unique(.x$era)
        )
      ))

[[1]]

Past versions of residual_patterns_by_model_domain-1.png

Version	Author	Date
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
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
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
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

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

7.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(basin, gamma_slab, lat_grid) %>% 
  mutate(eras = paste(lag(era), era, sep = " --> "),
         .resid_mean_offset = .resid_mean - lag(.resid_mean)
         ) %>% 
  ungroup() %>% 
  select(-era) %>% 
  drop_na() %>% 
  filter(eras != paste(unique(lat_residual$era)[1],
                       unique(lat_residual$era)[3],
                       sep = " --> "))

lat_residual_offset %>%
  ggplot(aes(lat_grid, .resid_mean_offset, col=eras)) +
  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
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 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
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
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] gt_0.2.2         corrr_0.4.3      broom_0.7.5      kableExtra_1.3.1
 [5] knitr_1.30       olsrr_0.5.3      GGally_2.0.0     lubridate_1.7.9 
 [9] metR_0.9.0       scico_1.2.0      patchwork_1.1.1  collapse_1.5.0  
[13] forcats_0.5.0    stringr_1.4.0    dplyr_1.0.5      purrr_0.3.4     
[17] readr_1.4.0      tidyr_1.1.2      tibble_3.0.4     ggplot2_3.3.3   
[21] tidyverse_1.3.0  workflowr_1.6.2 

loaded via a namespace (and not attached):
 [1] fs_1.5.0                 webshot_0.5.2            RColorBrewer_1.1-2      
 [4] httr_1.4.2               rprojroot_2.0.2          tools_4.0.3             
 [7] backports_1.1.10         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.13          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            plyr_1.8.6              
[55] grid_4.0.3               blob_1.2.1               parallel_4.0.3          
[58] promises_1.1.1           crayon_1.3.4             lattice_0.20-41         
[61] haven_2.3.1              hms_0.5.3                pillar_1.4.7            
[64] reprex_0.3.0             glue_1.4.2               evaluate_0.14           
[67] RcppArmadillo_0.10.1.2.0 data.table_1.13.2        modelr_0.1.8            
[70] vctrs_0.3.5              httpuv_1.5.4             cellranger_1.1.0        
[73] gtable_0.3.0             reshape_0.8.8            assertthat_0.2.1        
[76] xfun_0.18                openxlsx_4.2.3           RcppEigen_0.3.3.7.0     
[79] later_1.1.0.1            viridisLite_0.3.0        ellipsis_0.3.1          
[82] here_0.1