Last updated: 2021-04-13

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: 80e1ee3

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

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


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

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

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

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

File Version Author Date Message
html 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:

  • 50

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 ~ nitrate + silicate + phosphate + phosphate_star 174.46194 186.05290 1.1898481 2000-2009 –> 2010-2019 5.093228 820.4801 848.4882
mod Atlantic (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate_star 123.96210 135.55305 0.7252267 2000-2009 –> 2010-2019 5.112557 796.6284 824.6366
mod Atlantic (-Inf,26] cstar_tref ~ sal + nitrate + phosphate + phosphate_star 136.58932 148.18027 0.8208010 2000-2009 –> 2010-2019 4.984959 797.3525 825.3607
mod Atlantic (-Inf,26] cstar_tref ~ sal + nitrate + silicate + phosphate_star 87.37399 98.96494 0.5066282 2000-2009 –> 2010-2019 4.903610 760.5414 788.5495
mod Atlantic (-Inf,26] cstar_tref ~ sal + silicate + phosphate + phosphate_star 140.49803 152.08898 0.8528650 2000-2009 –> 2010-2019 5.040258 802.5299 830.5380
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate 759.32253 775.63352 6.8027215 2000-2009 –> 2010-2019 14.507992 1124.3285 1152.2305
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 758.05638 777.08587 6.7042506 2000-2009 –> 2010-2019 14.361261 1124.4215 1156.9738
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 758.08494 777.11443 6.7051053 2000-2009 –> 2010-2019 14.358472 1124.4015 1156.9538
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 750.41454 769.44403 6.4793901 2000-2009 –> 2010-2019 14.135191 1116.7636 1149.3158
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + silicate 112.65883 124.24978 0.6491527 2000-2009 –> 2010-2019 5.255886 796.4511 824.4593
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 111.71440 125.23718 0.6306816 2000-2009 –> 2010-2019 4.507887 758.1985 790.8747
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 751.53707 770.56656 6.5119416 2000-2009 –> 2010-2019 15.043120 1128.9290 1161.4813
obs Atlantic (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 758.23390 777.26339 6.7095657 2000-2009 –> 2010-2019 14.391773 1124.9341 1157.4864
mod Atlantic (-Inf,26] cstar_tref ~ sal + temp + nitrate + silicate 82.51548 94.10644 0.4830620 2000-2009 –> 2010-2019 5.489807 785.2927 813.3008
mod Atlantic (-Inf,26] cstar_tref ~ silicate + phosphate + phosphate_star 178.86516 188.52429 1.2669371 2000-2009 –> 2010-2019 5.675746 850.6450 873.9851
obs Atlantic (-Inf,26] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 751.48267 770.51216 6.5103605 2000-2009 –> 2010-2019 14.699080 1124.6957 1157.2479
obs Atlantic (-Inf,26] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 750.89289 769.92238 6.4932416 2000-2009 –> 2010-2019 14.675563 1124.0262 1156.5784
mod Atlantic (-Inf,26] cstar_tref ~ temp + aou + silicate + phosphate 175.95081 187.54176 1.2073433 2000-2009 –> 2010-2019 5.089615 820.7327 848.7408
obs Atlantic (-Inf,26] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 762.31520 778.62619 6.8942167 2000-2009 –> 2010-2019 15.056154 1133.1941 1161.0960
obs Atlantic (-Inf,26] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 748.91840 767.94789 6.4362572 2000-2009 –> 2010-2019 14.591774 1121.7170 1154.2693
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate 3542.30555 3568.28388 5.6265775 2000-2009 –> 2010-2019 11.803187 6470.0428 6520.6766
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 2064.97869 2092.52294 3.6473331 2000-2009 –> 2010-2019 7.934961 3407.0773 3458.7184
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 2065.21848 2092.76274 3.6484902 2000-2009 –> 2010-2019 8.005897 3414.7753 3466.4165
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + silicate 3544.02034 3574.32839 5.6251474 2000-2009 –> 2010-2019 11.709968 6460.2825 6519.3553
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + phosphate 3645.95271 3671.93104 6.1711086 2000-2009 –> 2010-2019 11.861832 6499.9510 6550.5849
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + phosphate_star 3649.18983 3675.16815 6.1889388 2000-2009 –> 2010-2019 11.982227 6519.2645 6569.8983
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 2029.41957 2056.96383 3.4797497 2000-2009 –> 2010-2019 8.218174 3417.7047 3469.3459
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 3646.63166 3676.93971 6.1638470 2000-2009 –> 2010-2019 11.845560 6501.2039 6560.2767
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2030.07551 2057.61977 3.4827702 2000-2009 –> 2010-2019 8.338476 3429.6565 3481.2976
obs Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 3649.70183 3680.00987 6.1807364 2000-2009 –> 2010-2019 11.953226 6518.5396 6577.6124
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 2066.12427 2093.66853 3.6528642 2000-2009 –> 2010-2019 8.082841 3423.3120 3474.9532
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 2020.11027 2047.65453 3.4371632 2000-2009 –> 2010-2019 8.598076 3447.8543 3499.4955
mod Atlantic (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 2030.45050 2057.99476 3.4844982 2000-2009 –> 2010-2019 8.377159 3433.5343 3485.1755
obs Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + phosphate 3605.76249 3631.74081 5.9539707 2000-2009 –> 2010-2019 11.601877 6452.9636 6503.5974
obs Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + phosphate_star 3606.82029 3632.79861 5.9595866 2000-2009 –> 2010-2019 11.686630 6466.5444 6517.1782
obs Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 3606.51994 3636.82799 5.9473807 2000-2009 –> 2010-2019 11.511604 6442.2862 6501.3590
obs Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 3607.55226 3637.86030 5.9528552 2000-2009 –> 2010-2019 11.574870 6452.6180 6511.6908
mod Atlantic (26,26.5] cstar_tref ~ temp + aou + silicate + phosphate 2108.24061 2131.84998 3.8723692 2000-2009 –> 2010-2019 8.640725 3497.4349 3541.6988
mod Atlantic (26,26.5] cstar_tref ~ temp + aou + silicate + phosphate_star 2105.45576 2129.06513 3.8581309 2000-2009 –> 2010-2019 8.725754 3504.1692 3548.4331
mod Atlantic (26,26.5] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 2138.71235 2162.32171 4.0316392 2000-2009 –> 2010-2019 8.718289 3519.9217 3564.1856
obs Atlantic (26.5,26.75] cstar_tref ~ aou + phosphate 5221.58680 5240.43376 5.7681014 2000-2009 –> 2010-2019 13.228730 9179.3090 9215.5804
obs Atlantic (26.5,26.75] cstar_tref ~ aou + silicate + phosphate 5212.18379 5235.74249 5.7282315 2000-2009 –> 2010-2019 13.045000 9149.4754 9194.8146
obs Atlantic (26.5,26.75] cstar_tref ~ aou + silicate + phosphate_star 5229.16297 5252.72167 5.7876990 2000-2009 –> 2010-2019 13.720246 9259.5424 9304.8816
obs Atlantic (26.5,26.75] cstar_tref ~ sal + aou + nitrate 5178.52762 5202.08632 5.6121546 2000-2009 –> 2010-2019 13.670654 9227.0547 9272.3939
obs Atlantic (26.5,26.75] cstar_tref ~ sal + aou + phosphate_star 5243.57169 5267.13039 5.8386478 2000-2009 –> 2010-2019 12.225883 9024.3442 9069.6834
obs Atlantic (26.5,26.75] cstar_tref ~ sal + nitrate + phosphate_star 5236.93439 5260.49309 5.8151230 2000-2009 –> 2010-2019 12.102600 8999.5726 9044.9118
obs Atlantic (26.5,26.75] cstar_tref ~ sal + nitrate + silicate + phosphate_star 5238.25652 5266.52696 5.8127258 2000-2009 –> 2010-2019 11.984487 8981.4956 9035.9027
obs Atlantic (26.5,26.75] cstar_tref ~ sal + phosphate_star 5246.24278 5265.08974 5.8552607 2000-2009 –> 2010-2019 13.057215 9163.3140 9199.5854
obs Atlantic (26.5,26.75] cstar_tref ~ sal + silicate + phosphate_star 5236.49315 5260.05185 5.8135625 2000-2009 –> 2010-2019 12.311653 9037.0879 9082.4271
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 2690.47339 2719.66067 3.9777639 2000-2009 –> 2010-2019 7.145938 4592.3904 4648.9152
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 2693.58533 2722.77261 3.9907333 2000-2009 –> 2010-2019 7.049442 4569.6929 4626.2177
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 2761.10287 2790.29014 4.2827701 2000-2009 –> 2010-2019 7.135396 4586.0118 4642.5366
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2767.31203 2796.49930 4.3106770 2000-2009 –> 2010-2019 7.109310 4578.1951 4634.7199
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 2701.31558 2730.50286 4.0231333 2000-2009 –> 2010-2019 6.948433 4544.6899 4601.2147
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 2792.63200 2821.81927 4.4263720 2000-2009 –> 2010-2019 7.351880 4636.0586 4692.5834
mod Atlantic (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 2768.34174 2797.52901 4.3153226 2000-2009 –> 2010-2019 7.110017 4578.1913 4634.7161
mod Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + phosphate_star 2724.82787 2749.84554 4.1319420 2000-2009 –> 2010-2019 7.349003 4635.9843 4684.4342
mod Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 2729.83630 2754.85396 4.1536458 2000-2009 –> 2010-2019 7.218559 4605.4314 4653.8812
obs Atlantic (26.5,26.75] cstar_tref ~ temp + phosphate_star 5196.19354 5215.04050 5.6796917 2000-2009 –> 2010-2019 11.994152 8961.7651 8998.0365
mod Atlantic (26.5,26.75] cstar_tref ~ temp + silicate + phosphate + phosphate_star 2773.35190 2798.36956 4.3470823 2000-2009 –> 2010-2019 7.359998 4636.3873 4684.8372
obs Atlantic (26.75,27] cstar_tref ~ aou + phosphate 6172.01976 6191.73392 4.9514495 2000-2009 –> 2010-2019 12.475734 10682.5871 10720.2398
obs Atlantic (26.75,27] cstar_tref ~ aou + silicate + phosphate 6150.37816 6175.02085 4.8944538 2000-2009 –> 2010-2019 12.405104 10660.5695 10707.6354
mod Atlantic (26.75,27] cstar_tref ~ sal + aou + phosphate 5330.91506 5355.00968 4.4311774 2000-2009 –> 2010-2019 7.984348 8704.7623 8751.0537
obs Atlantic (26.75,27] cstar_tref ~ sal + aou + phosphate_star 6237.05686 6261.69955 5.1066858 2000-2009 –> 2010-2019 12.885287 10793.1696 10840.2355
mod Atlantic (26.75,27] cstar_tref ~ sal + aou + silicate + phosphate 4557.65080 4586.56435 2.9009076 2000-2009 –> 2010-2019 5.436797 7511.2037 7566.7533
mod Atlantic (26.75,27] cstar_tref ~ sal + aou + silicate + phosphate_star 4651.36321 4680.27675 3.0533294 2000-2009 –> 2010-2019 6.018610 7800.7626 7856.3123
obs Atlantic (26.75,27] cstar_tref ~ sal + nitrate + phosphate_star 6244.46365 6269.10634 5.1252425 2000-2009 –> 2010-2019 12.931251 10805.1839 10852.2498
mod Atlantic (26.75,27] cstar_tref ~ sal + nitrate + silicate + phosphate_star 4827.08935 4856.00290 3.3610648 2000-2009 –> 2010-2019 6.486461 8042.3310 8097.8807
obs Atlantic (26.75,27] cstar_tref ~ sal + nitrate + silicate + phosphate_star 6219.96120 6249.53243 5.0591535 2000-2009 –> 2010-2019 12.838749 10778.2413 10834.7204
obs Atlantic (26.75,27] cstar_tref ~ sal + phosphate + phosphate_star 6230.80706 6255.44975 5.0910800 2000-2009 –> 2010-2019 12.843763 10782.5476 10829.6135
mod Atlantic (26.75,27] cstar_tref ~ sal + silicate + phosphate + phosphate_star 4691.30913 4720.22267 3.1207114 2000-2009 –> 2010-2019 6.244877 7906.0577 7961.6073
obs Atlantic (26.75,27] cstar_tref ~ sal + silicate + phosphate + phosphate_star 6222.15136 6251.72259 5.0645826 2000-2009 –> 2010-2019 12.813118 10775.1908 10831.6699
obs Atlantic (26.75,27] cstar_tref ~ sal + silicate + phosphate_star 6220.22964 6244.87233 5.0647768 2000-2009 –> 2010-2019 12.849686 10777.4043 10824.4701
mod Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + phosphate 5323.98779 5352.90134 4.4096135 2000-2009 –> 2010-2019 7.959598 8698.7119 8754.2616
obs Atlantic (26.75,27] cstar_tref ~ silicate + phosphate + phosphate_star 6226.03843 6250.68112 5.0792049 2000-2009 –> 2010-2019 12.930271 10794.2986 10841.3644
mod Atlantic (26.75,27] cstar_tref ~ temp + aou + silicate + phosphate 4677.98472 4706.89826 3.0980717 2000-2009 –> 2010-2019 5.488690 7557.6790 7613.2286
mod Atlantic (26.75,27] cstar_tref ~ temp + aou + silicate + phosphate_star 4697.64956 4726.56311 3.1315425 2000-2009 –> 2010-2019 5.585263 7609.9628 7665.5125
mod Atlantic (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate_star 4949.07018 4977.98373 3.5927360 2000-2009 –> 2010-2019 6.395333 8027.8255 8083.3751
obs Atlantic (26.75,27] cstar_tref ~ temp + phosphate_star 6130.69301 6150.40716 4.8522474 2000-2009 –> 2010-2019 11.541379 10487.1365 10524.7892
mod Atlantic (26.75,27] cstar_tref ~ temp + silicate + phosphate + phosphate_star 4705.03561 4733.94915 3.1442073 2000-2009 –> 2010-2019 5.629189 7633.1989 7688.7486
mod Atlantic (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate_star 5865.81736 5896.65532 2.4649283 2000-2009 –> 2010-2019 5.197406 10139.4682 10198.9789
mod Atlantic (27,27.25] cstar_tref ~ sal + nitrate + silicate + phosphate_star 6001.50404 6032.34200 2.6011768 2000-2009 –> 2010-2019 5.538160 10402.0372 10461.5479
mod Atlantic (27,27.25] cstar_tref ~ sal + silicate + phosphate + phosphate_star 5935.91154 5966.74950 2.5343972 2000-2009 –> 2010-2019 5.363741 10270.8045 10330.3152
mod Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + phosphate 6401.31293 6432.15089 3.0480208 2000-2009 –> 2010-2019 6.346066 11005.6807 11065.1914
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + phosphate_star 11085.68847 11118.36097 6.1423012 2000-2009 –> 2010-2019 12.053172 19425.7653 19489.4724
mod Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5747.00124 5782.97886 2.3496299 2000-2009 –> 2010-2019 4.774287 9812.5385 9881.9677
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + nitrate + phosphate_star 11185.98466 11218.65716 6.3248832 2000-2009 –> 2010-2019 12.440664 19614.8719 19678.5789
mod Atlantic (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 5896.23451 5932.21213 2.4928594 2000-2009 –> 2010-2019 5.142403 10117.7015 10187.1306
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 10953.86382 10991.98174 5.9068645 2000-2009 –> 2010-2019 11.774852 19276.9653 19351.2901
mod Atlantic (27,27.25] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5766.10178 5802.07940 2.3674926 2000-2009 –> 2010-2019 4.806491 9842.0069 9911.4360
mod Atlantic (27,27.25] cstar_tref ~ sal + temp + silicate + phosphate_star 6291.81781 6322.65577 2.9185195 2000-2009 –> 2010-2019 6.144655 10857.4307 10916.9414
obs Atlantic (27,27.25] cstar_tref ~ sal + temp + silicate + phosphate_star 10957.05590 10989.72841 5.9158284 2000-2009 –> 2010-2019 11.972814 19360.7688 19424.4759
obs Atlantic (27,27.25] cstar_tref ~ temp + aou + phosphate_star 11090.71035 11117.93744 6.1549106 2000-2009 –> 2010-2019 12.233073 19501.5186 19554.6078
mod Atlantic (27,27.25] cstar_tref ~ temp + aou + silicate + phosphate_star 6488.26687 6519.10483 3.1549436 2000-2009 –> 2010-2019 6.123696 10907.7143 10967.2250
obs Atlantic (27,27.25] cstar_tref ~ temp + aou + silicate + phosphate_star 10942.40090 10975.07340 5.8905622 2000-2009 –> 2010-2019 11.860092 19308.2121 19371.9191
obs Atlantic (27,27.25] cstar_tref ~ temp + nitrate + silicate + phosphate_star 10957.77492 10990.44742 5.9170708 2000-2009 –> 2010-2019 11.994651 19370.3336 19434.0407
obs Atlantic (27,27.25] cstar_tref ~ temp + phosphate + phosphate_star 11090.70853 11117.93562 6.1549073 2000-2009 –> 2010-2019 12.249818 19508.6879 19561.7771
mod Atlantic (27,27.25] cstar_tref ~ temp + silicate + phosphate + phosphate_star 6497.06219 6527.90015 3.1659655 2000-2009 –> 2010-2019 6.157129 10929.7306 10989.2413
obs Atlantic (27,27.25] cstar_tref ~ temp + silicate + phosphate + phosphate_star 10942.52895 10975.20145 5.8907825 2000-2009 –> 2010-2019 11.871918 19313.4016 19377.1086
obs Atlantic (27,27.25] cstar_tref ~ temp + silicate + phosphate_star 10965.45815 10992.68523 5.9338282 2000-2009 –> 2010-2019 12.210959 19460.2074 19513.2966
mod Atlantic (27.25,27.5] cstar_tref ~ sal + aou + phosphate 11143.78791 11171.87671 3.7358703 2000-2009 –> 2010-2019 7.586793 21353.8760 21409.5605
mod Atlantic (27.25,27.5] cstar_tref ~ sal + nitrate + silicate + phosphate_star 10964.25619 10997.96275 3.5728246 2000-2009 –> 2010-2019 7.709551 21440.2324 21507.0539
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + phosphate 11139.43908 11173.14564 3.7300444 2000-2009 –> 2010-2019 7.566337 21337.4969 21404.3184
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + phosphate 13218.85102 13252.62796 5.9877938 2000-2009 –> 2010-2019 10.521159 22221.3276 22287.1183
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + phosphate_star 11049.27730 11082.98385 3.6482824 2000-2009 –> 2010-2019 7.359512 21125.1697 21191.9912
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + phosphate_star 13231.52478 13265.30172 6.0062595 2000-2009 –> 2010-2019 10.528432 22226.4175 22292.2082
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 13208.14657 13247.55300 5.9693404 2000-2009 –> 2010-2019 10.501337 22211.6966 22288.4524
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 10814.51811 10853.84242 3.4420117 2000-2009 –> 2010-2019 7.106895 20846.0889 20924.0473
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 13222.16412 13261.57055 5.9897044 2000-2009 –> 2010-2019 10.511724 22218.9531 22295.7090
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + phosphate_star 11029.86720 11063.57376 3.6309164 2000-2009 –> 2010-2019 7.192958 20954.5258 21021.3473
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 10793.29727 10832.62159 3.4241031 2000-2009 –> 2010-2019 6.962146 20695.0375 20772.9958
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + phosphate + phosphate_star 11050.44531 11084.15187 3.6493300 2000-2009 –> 2010-2019 7.363320 21129.0779 21195.8994
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + phosphate + phosphate_star 13224.53570 13258.31264 5.9960694 2000-2009 –> 2010-2019 10.518930 22219.8946 22285.6853
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 10815.55554 10854.87986 3.4428896 2000-2009 –> 2010-2019 7.110266 20849.6336 20927.5920
obs Atlantic (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 13215.04366 13254.45008 5.9793515 2000-2009 –> 2010-2019 10.502193 22212.3900 22289.1459
mod Atlantic (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate_star 10895.21415 10928.92070 3.5126984 2000-2009 –> 2010-2019 7.307209 21052.9072 21119.7287
obs Atlantic (27.25,27.5] cstar_tref ~ temp + aou + phosphate 13224.72472 13252.87216 5.9992591 2000-2009 –> 2010-2019 10.700268 22336.6083 22391.4339
obs Atlantic (27.25,27.5] cstar_tref ~ temp + aou + silicate + phosphate 13224.62279 13258.39973 5.9961963 2000-2009 –> 2010-2019 10.683454 22329.5190 22395.3098
obs Atlantic (27.25,27.5] cstar_tref ~ temp + aou + silicate + phosphate_star 13241.86832 13275.64525 6.0213723 2000-2009 –> 2010-2019 10.698796 22340.3210 22406.1117
obs Atlantic (27.25,27.5] cstar_tref ~ temp + silicate + phosphate + phosphate_star 13233.35853 13267.13547 6.0089360 2000-2009 –> 2010-2019 10.686747 22332.0651 22397.8558
obs Atlantic (27.5,27.75] cstar_tref ~ aou + silicate + phosphate 14610.76633 14639.65320 5.1589645 2000-2009 –> 2010-2019 10.745994 26123.4778 26179.9304
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + phosphate 18972.26853 19003.86590 2.4396203 2000-2009 –> 2010-2019 4.973560 34986.4050 35048.6703
mod Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate 18299.35839 18337.27523 2.2470011 2000-2009 –> 2010-2019 4.761774 34263.7561 34338.4745
obs Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate 14550.58481 14585.24905 5.0921768 2000-2009 –> 2010-2019 10.536461 25970.4667 26038.2099
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + phosphate 18971.18823 19009.10507 2.4387047 2000-2009 –> 2010-2019 4.945705 34914.4928 34989.2111
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + phosphate 14814.41175 14849.07600 5.3816329 2000-2009 –> 2010-2019 10.723595 26164.7755 26232.5188
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + phosphate_star 18858.94299 18896.85983 2.4055742 2000-2009 –> 2010-2019 4.867933 34679.8200 34754.5384
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 18290.59970 18334.83602 2.2440571 2000-2009 –> 2010-2019 4.738568 34201.8748 34289.0462
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 14542.77887 14583.22049 5.0817236 2000-2009 –> 2010-2019 10.280056 25795.2795 25874.3132
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 18083.32663 18127.56295 2.1880849 2000-2009 –> 2010-2019 4.639353 33875.4423 33962.6137
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 14569.15595 14609.59756 5.1098905 2000-2009 –> 2010-2019 10.317062 25827.8813 25906.9151
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + nitrate + phosphate_star 19061.07106 19098.98790 2.4655634 2000-2009 –> 2010-2019 4.939850 34914.8756 34989.5940
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 18181.84211 18226.07842 2.2145118 2000-2009 –> 2010-2019 4.669754 33984.9975 34072.1689
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + phosphate + phosphate_star 18856.45274 18894.36958 2.4048443 2000-2009 –> 2010-2019 4.869478 34683.6238 34758.3422
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + phosphate + phosphate_star 14831.04913 14865.71337 5.4004285 2000-2009 –> 2010-2019 10.768787 26199.4731 26267.2163
mod Atlantic (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 18089.35385 18133.59016 2.1896926 2000-2009 –> 2010-2019 4.643085 33887.3746 33974.5461
obs Atlantic (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 14549.26822 14589.70984 5.0886389 2000-2009 –> 2010-2019 10.287221 25801.9448 25880.9785
obs Atlantic (27.5,27.75] cstar_tref ~ temp + aou + silicate + phosphate 14541.92975 14576.59400 5.0829494 2000-2009 –> 2010-2019 10.364514 25850.6315 25918.3747
obs Atlantic (27.5,27.75] cstar_tref ~ temp + aou + silicate + phosphate_star 14568.68140 14603.34564 5.1115242 2000-2009 –> 2010-2019 10.409532 25888.7725 25956.5158
obs Atlantic (27.5,27.75] cstar_tref ~ temp + silicate + phosphate + phosphate_star 14548.53965 14583.20390 5.0899948 2000-2009 –> 2010-2019 10.375481 25859.9604 25927.7036
obs Atlantic (27.75,27.85] cstar_tref ~ aou + silicate + phosphate 9121.87617 9148.25549 5.6630340 2000-2009 –> 2010-2019 10.922805 16570.8029 16622.6659
mod Atlantic (27.75,27.85] cstar_tref ~ sal + aou + nitrate + silicate 5054.89144 5086.78674 1.2937625 2000-2009 –> 2010-2019 2.539160 9356.2059 9419.1633
mod Atlantic (27.75,27.85] cstar_tref ~ sal + aou + silicate + phosphate 4812.38721 4844.28251 1.1935534 2000-2009 –> 2010-2019 2.396039 9021.9029 9084.8604
obs Atlantic (27.75,27.85] cstar_tref ~ sal + aou + silicate + phosphate_star 9091.96080 9123.61599 5.6008390 2000-2009 –> 2010-2019 10.856499 16540.9985 16603.2340
mod Atlantic (27.75,27.85] cstar_tref ~ sal + nitrate + silicate + phosphate_star 5446.89757 5478.79287 1.4738467 2000-2009 –> 2010-2019 2.856069 10021.1056 10084.0630
mod Atlantic (27.75,27.85] cstar_tref ~ sal + silicate + phosphate + phosphate_star 5436.67463 5468.56993 1.4688463 2000-2009 –> 2010-2019 2.858058 10024.0870 10087.0444
obs Atlantic (27.75,27.85] cstar_tref ~ sal + silicate + phosphate + phosphate_star 9087.45224 9119.10743 5.5921082 2000-2009 –> 2010-2019 10.750281 16491.2548 16553.4903
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + silicate 5010.42043 5047.63161 1.2739286 2000-2009 –> 2010-2019 2.516164 9307.0783 9380.5286
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 4685.02571 4722.23690 1.1433117 2000-2009 –> 2010-2019 2.340681 8885.3784 8958.8288
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4817.32302 4854.53420 1.1947188 2000-2009 –> 2010-2019 2.432787 9105.1850 9178.6353
obs Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 9080.34097 9117.27202 5.5745058 2000-2009 –> 2010-2019 10.724770 16482.4364 16555.0445
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 4824.44441 4861.65559 1.1975507 2000-2009 –> 2010-2019 2.433277 9107.3483 9180.7987
obs Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate 9136.71073 9168.36592 5.6882396 2000-2009 –> 2010-2019 11.259311 16726.5572 16788.7927
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4810.54116 4847.75234 1.1920282 2000-2009 –> 2010-2019 2.425936 9089.5891 9163.0394
obs Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 9087.60686 9124.53791 5.5885386 2000-2009 –> 2010-2019 10.729664 16485.4117 16558.0198
mod Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate_star 4920.94340 4952.83870 1.2374144 2000-2009 –> 2010-2019 2.558817 9377.3456 9440.3031
obs Atlantic (27.75,27.85] cstar_tref ~ silicate + phosphate + phosphate_star 9135.45951 9161.83884 5.6897136 2000-2009 –> 2010-2019 11.047372 16628.9364 16680.7993
obs Atlantic (27.75,27.85] cstar_tref ~ temp + aou + silicate + phosphate_star 9097.63181 9129.28699 5.6118403 2000-2009 –> 2010-2019 10.765997 16499.5528 16561.7883
obs Atlantic (27.75,27.85] cstar_tref ~ temp + silicate + phosphate 9154.63964 9181.01896 5.7276002 2000-2009 –> 2010-2019 11.304006 16744.7984 16796.6613
obs Atlantic (27.75,27.85] cstar_tref ~ temp + silicate + phosphate + phosphate_star 9104.58792 9136.24311 5.6253640 2000-2009 –> 2010-2019 10.772106 16503.0307 16565.2662
mod Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate 5971.52433 6003.30281 1.8243782 2000-2009 –> 2010-2019 4.732364 12517.7812 12580.6447
mod Atlantic (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate 5743.65682 5775.43530 1.6887625 2000-2009 –> 2010-2019 4.413245 12118.6155 12181.4790
obs Atlantic (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate_star 18659.77671 18696.06752 4.7630659 2000-2009 –> 2010-2019 10.165265 35548.3930 35620.1275
mod Atlantic (27.85,27.95] cstar_tref ~ sal + nitrate + silicate + phosphate_star 6223.31955 6255.09802 1.9869355 2000-2009 –> 2010-2019 4.817687 12698.8352 12761.6986
obs Atlantic (27.85,27.95] cstar_tref ~ sal + silicate + phosphate 19363.46142 19393.70376 5.3316287 2000-2009 –> 2010-2019 11.092434 36599.3277 36659.1065
mod Atlantic (27.85,27.95] cstar_tref ~ sal + silicate + phosphate + phosphate_star 6207.19079 6238.96927 1.9761019 2000-2009 –> 2010-2019 4.759991 12638.8362 12701.6997
obs Atlantic (27.85,27.95] cstar_tref ~ sal + silicate + phosphate + phosphate_star 19062.70657 19098.99738 5.0798306 2000-2009 –> 2010-2019 10.639425 36107.3820 36179.1165
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + silicate 5813.07902 5850.15391 1.7278035 2000-2009 –> 2010-2019 4.623019 12349.7696 12423.1103
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + phosphate 6385.69914 6417.47762 2.0993703 2000-2009 –> 2010-2019 4.840760 12776.9157 12839.7792
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 5435.28533 5472.36022 1.5201140 2000-2009 –> 2010-2019 4.244380 11812.0346 11885.3754
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5545.72357 5582.79846 1.5781006 2000-2009 –> 2010-2019 4.327110 11946.2340 12019.5747
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 18506.81689 18549.15617 4.6465717 2000-2009 –> 2010-2019 10.004724 35352.9453 35436.6355
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 5553.86696 5590.94185 1.5824630 2000-2009 –> 2010-2019 4.375365 11996.0069 12069.3476
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate 18560.92728 18597.21809 4.6884211 2000-2009 –> 2010-2019 10.108136 35467.1338 35538.8683
obs Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 18538.30705 18580.64633 4.6700121 2000-2009 –> 2010-2019 10.037131 35393.5220 35477.2122
mod Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate_star 5659.79066 5691.56914 1.6414283 2000-2009 –> 2010-2019 4.662716 12306.4967 12369.3601
obs Atlantic (27.85,27.95] cstar_tref ~ silicate + phosphate + phosphate_star 19722.41074 19752.65308 5.6463835 2000-2009 –> 2010-2019 11.461257 37009.0396 37068.8183
obs Atlantic (27.85,27.95] cstar_tref ~ temp + aou + silicate + phosphate_star 19631.39103 19667.68184 5.5630756 2000-2009 –> 2010-2019 11.332877 36877.7358 36949.4703
obs Atlantic (27.85,27.95] cstar_tref ~ temp + silicate + phosphate 19743.04391 19773.28625 5.6650308 2000-2009 –> 2010-2019 11.484333 37033.8099 37093.5887
obs Atlantic (27.85,27.95] cstar_tref ~ temp + silicate + phosphate + phosphate_star 19722.61991 19758.91072 5.6447679 2000-2009 –> 2010-2019 11.450601 37002.7945 37074.5290
mod Atlantic (27.95,28.05] cstar_tref ~ sal + aou + phosphate 9602.50852 9629.86818 3.7036033 2000-2009 –> 2010-2019 7.743635 17663.9551 17717.6419
obs Atlantic (27.95,28.05] cstar_tref ~ sal + aou + phosphate_star 21177.73803 21208.28416 5.8367840 2000-2009 –> 2010-2019 12.726008 40015.2365 40075.4891
mod Atlantic (27.95,28.05] cstar_tref ~ sal + aou + silicate + phosphate 9551.19731 9584.02890 3.6478711 2000-2009 –> 2010-2019 7.643201 17582.8220 17647.2461
obs Atlantic (27.95,28.05] cstar_tref ~ sal + aou + silicate + phosphate 21042.12619 21078.78154 5.7172416 2000-2009 –> 2010-2019 12.565831 39848.3663 39920.6694
obs Atlantic (27.95,28.05] cstar_tref ~ sal + aou + silicate + phosphate_star 20635.75890 20672.41425 5.3783336 2000-2009 –> 2010-2019 12.094359 39332.1111 39404.4142
obs Atlantic (27.95,28.05] cstar_tref ~ sal + nitrate + silicate + phosphate_star 20878.25946 20914.91481 5.5780812 2000-2009 –> 2010-2019 12.507027 39750.0805 39822.3836
obs Atlantic (27.95,28.05] cstar_tref ~ sal + phosphate + phosphate_star 21164.08344 21194.62956 5.8248115 2000-2009 –> 2010-2019 12.686041 39978.6910 40038.9436
obs Atlantic (27.95,28.05] cstar_tref ~ sal + silicate + phosphate + phosphate_star 20798.34605 20835.00140 5.5114502 2000-2009 –> 2010-2019 12.258039 39520.2248 39592.5279
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate 9456.17578 9489.00737 3.5506060 2000-2009 –> 2010-2019 7.375632 17363.2168 17427.6410
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate 21208.06517 21244.72052 5.8617001 2000-2009 –> 2010-2019 12.732439 40032.4586 40104.7617
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate_star 9486.47589 9519.30748 3.5813365 2000-2009 –> 2010-2019 7.404166 17391.8741 17456.2982
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate_star 21139.77194 21176.42729 5.8018105 2000-2009 –> 2010-2019 12.641580 39938.7672 40011.0704
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 9299.02647 9337.33000 3.3934733 2000-2009 –> 2010-2019 7.218437 17208.0208 17283.1823
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 9310.33109 9348.63461 3.4044016 2000-2009 –> 2010-2019 7.227096 17217.6279 17292.7894
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + phosphate + phosphate_star 21158.76258 21195.41793 5.8184026 2000-2009 –> 2010-2019 12.666475 39964.5783 40036.8815
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + phosphate_star 9567.07911 9594.43877 3.6664709 2000-2009 –> 2010-2019 7.496053 17475.5247 17529.2115
obs Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate 20830.15614 20866.81149 5.5378773 2000-2009 –> 2010-2019 12.476708 39709.9918 39782.2949
mod Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate_star 9310.07745 9342.90904 3.4060929 2000-2009 –> 2010-2019 7.231439 17217.3576 17281.7818
mod Atlantic (27.95,28.05] cstar_tref ~ temp + aou + silicate + phosphate 9599.71137 9632.54297 3.6985536 2000-2009 –> 2010-2019 7.753350 17673.5909 17738.0151
mod Atlantic (27.95,28.05] cstar_tref ~ temp + aou + silicate + phosphate_star 9629.57074 9662.40233 3.7300971 2000-2009 –> 2010-2019 7.806915 17718.9411 17783.3653
obs Atlantic (28.05,28.1] cstar_tref ~ sal + aou + phosphate 17673.99055 17703.63799 5.8143292 2000-2009 –> 2010-2019 13.326656 31185.4418 31243.0055
mod Atlantic (28.05,28.1] cstar_tref ~ sal + aou + phosphate_star 4065.71127 4090.35885 1.7599135 2000-2009 –> 2010-2019 3.341207 6966.5708 7014.4503
obs Atlantic (28.05,28.1] cstar_tref ~ sal + aou + phosphate_star 17539.67269 17569.32012 5.6754516 2000-2009 –> 2010-2019 13.142895 31027.5734 31085.1370
mod Atlantic (28.05,28.1] cstar_tref ~ sal + aou + silicate + phosphate_star 4049.12370 4078.70080 1.7439818 2000-2009 –> 2010-2019 3.268688 6895.8637 6953.3191
obs Atlantic (28.05,28.1] cstar_tref ~ sal + nitrate + phosphate_star 17508.09951 17537.74695 5.6432911 2000-2009 –> 2010-2019 13.143299 31013.1008 31070.6644
mod Atlantic (28.05,28.1] cstar_tref ~ sal + phosphate 3888.21296 3907.93102 1.6151122 2000-2009 –> 2010-2019 3.405518 6978.3395 7016.6431
obs Atlantic (28.05,28.1] cstar_tref ~ sal + phosphate + phosphate_star 17546.54237 17576.18980 5.6824734 2000-2009 –> 2010-2019 13.151519 31035.2862 31092.8498
mod Atlantic (28.05,28.1] cstar_tref ~ sal + silicate + phosphate 3843.23537 3867.88296 1.5784151 2000-2009 –> 2010-2019 3.278292 6855.4565 6903.3360
obs Atlantic (28.05,28.1] cstar_tref ~ sal + silicate + phosphate_star 17599.66652 17629.31395 5.7370674 2000-2009 –> 2010-2019 13.279261 31126.7111 31184.2748
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate 17671.36485 17706.94177 5.8094905 2000-2009 –> 2010-2019 13.321719 31184.7646 31253.8410
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate_star 17501.70037 17537.27729 5.6347665 2000-2009 –> 2010-2019 13.053737 30966.0072 31035.0836
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate 3973.57863 3998.22621 1.6823472 2000-2009 –> 2010-2019 3.496923 7086.3554 7134.2349
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + phosphate_star 17478.73464 17514.31156 5.6115232 2000-2009 –> 2010-2019 13.038654 30947.3618 31016.4382
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + phosphate 3878.48341 3903.13099 1.6058703 2000-2009 –> 2010-2019 3.374542 6951.8013 6999.6809
mod Atlantic (28.05,28.1] cstar_tref ~ sal + temp + phosphate_star 3755.25218 3779.89976 1.5119143 2000-2009 –> 2010-2019 2.748465 6277.3998 6325.2793
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + phosphate_star 17534.78379 17564.43122 5.6704598 2000-2009 –> 2010-2019 13.143183 31025.4619 31083.0256
obs Atlantic (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate_star 17528.34337 17563.92029 5.6618521 2000-2009 –> 2010-2019 13.105833 31005.8764 31074.9528
mod Atlantic (28.05,28.1] cstar_tref ~ temp + phosphate_star 3756.22418 3775.94225 1.5141143 2000-2009 –> 2010-2019 2.758304 6285.8559 6324.1595
mod Atlantic (28.05,28.1] cstar_tref ~ temp + silicate + phosphate 3861.04196 3885.68955 1.5922257 2000-2009 –> 2010-2019 3.435879 6998.3012 7046.1807
mod Atlantic (28.05,28.1] cstar_tref ~ temp + silicate + phosphate_star 3639.09903 3663.74661 1.4283933 2000-2009 –> 2010-2019 2.647323 6139.1426 6187.0222
obs Atlantic (28.1,28.15] cstar_tref ~ sal + aou + phosphate_star 14205.11202 14233.46416 6.6292303 2000-2009 –> 2010-2019 15.881098 23601.4480 23655.6044
mod Atlantic (28.1,28.15] cstar_tref ~ sal + nitrate 5290.41776 5311.64389 1.4243063 2000-2009 –> 2010-2019 2.973393 10062.4634 10104.3173
mod Atlantic (28.1,28.15] cstar_tref ~ sal + nitrate + silicate 5277.25369 5303.78635 1.4170769 2000-2009 –> 2010-2019 2.958687 10038.8857 10091.2032
mod Atlantic (28.1,28.15] cstar_tref ~ sal + phosphate 5282.22241 5303.44854 1.4203947 2000-2009 –> 2010-2019 2.950031 10021.8455 10063.6995
obs Atlantic (28.1,28.15] cstar_tref ~ sal + phosphate 14920.09824 14942.77995 7.8357478 2000-2009 –> 2010-2019 18.320899 24636.7810 24680.1061
obs Atlantic (28.1,28.15] cstar_tref ~ sal + phosphate_star 14207.21772 14229.89943 6.6355808 2000-2009 –> 2010-2019 16.028214 23640.4520 23683.7771
mod Atlantic (28.1,28.15] cstar_tref ~ sal + silicate + phosphate 5022.40993 5048.94258 1.3009285 2000-2009 –> 2010-2019 2.703448 9541.4070 9593.7244
obs Atlantic (28.1,28.15] cstar_tref ~ sal + silicate + phosphate_star 14207.43602 14235.78816 6.6328242 2000-2009 –> 2010-2019 15.923228 23614.4794 23668.6357
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate 5278.40958 5304.94224 1.4176267 2000-2009 –> 2010-2019 2.950937 10026.1888 10078.5062
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + silicate 5273.38148 5305.22067 1.4142873 2000-2009 –> 2010-2019 2.921333 9978.8267 10041.6076
mod Atlantic (28.1,28.15] cstar_tref ~ sal + temp + phosphate 5157.40954 5183.94219 1.3612183 2000-2009 –> 2010-2019 2.887277 9893.0239 9945.3413
obs Atlantic (28.1,28.15] cstar_tref ~ temp + aou + phosphate_star 13888.73686 13917.08900 6.1577238 2000-2009 –> 2010-2019 15.037372 23179.2934 23233.4497
obs Atlantic (28.1,28.15] cstar_tref ~ temp + nitrate 14919.73087 14942.41258 7.8350765 2000-2009 –> 2010-2019 18.761241 24742.5454 24785.8705
obs Atlantic (28.1,28.15] cstar_tref ~ temp + nitrate + phosphate_star 13909.73080 13938.08294 6.1879457 2000-2009 –> 2010-2019 15.080315 23203.9750 23258.1313
mod Atlantic (28.1,28.15] cstar_tref ~ temp + nitrate + silicate 5411.31754 5437.85020 1.4822840 2000-2009 –> 2010-2019 3.178411 10418.0529 10470.3704
mod Atlantic (28.1,28.15] cstar_tref ~ temp + phosphate 5226.06793 5247.29405 1.3938796 2000-2009 –> 2010-2019 2.993128 10079.8880 10121.7419
obs Atlantic (28.1,28.15] cstar_tref ~ temp + phosphate 14777.00618 14799.68790 7.5785809 2000-2009 –> 2010-2019 18.086744 24499.3365 24542.6616
obs Atlantic (28.1,28.15] cstar_tref ~ temp + phosphate_star 13920.16562 13942.84733 6.2059163 2000-2009 –> 2010-2019 15.161251 23230.5859 23273.9110
mod Atlantic (28.1,28.15] cstar_tref ~ temp + silicate + phosphate 5107.47189 5134.00455 1.3385976 2000-2009 –> 2010-2019 2.875662 9861.5251 9913.8425
obs Atlantic (28.1,28.15] cstar_tref ~ temp + silicate + phosphate_star 13880.71893 13909.07107 6.1462205 2000-2009 –> 2010-2019 15.006034 23165.5149 23219.6713
mod Atlantic (28.15,28.2] cstar_tref ~ aou + nitrate + silicate 8687.63236 8717.46356 1.0904052 2000-2009 –> 2010-2019 2.293078 15663.1722 15721.4204
obs Atlantic (28.15,28.2] cstar_tref ~ aou + phosphate 4476.46303 4494.98576 4.6116995 2000-2009 –> 2010-2019 9.538857 8179.2603 8215.4564
mod Atlantic (28.15,28.2] cstar_tref ~ aou + silicate + phosphate 8970.27695 9000.10815 1.1452072 2000-2009 –> 2010-2019 2.336377 15904.0663 15962.3145
obs Atlantic (28.15,28.2] cstar_tref ~ nitrate + phosphate_star 4569.75619 4588.27892 4.9044133 2000-2009 –> 2010-2019 10.271773 8377.4669 8413.6631
obs Atlantic (28.15,28.2] cstar_tref ~ sal + aou + phosphate_star 4168.07490 4191.22832 3.7578799 2000-2009 –> 2010-2019 8.473858 7819.1660 7864.4113
mod Atlantic (28.15,28.2] cstar_tref ~ sal + nitrate 8536.79160 8560.65655 1.0626087 2000-2009 –> 2010-2019 2.329607 15736.6697 15783.2683
mod Atlantic (28.15,28.2] cstar_tref ~ sal + nitrate + phosphate_star 8471.16785 8500.99905 1.0502149 2000-2009 –> 2010-2019 2.214224 15304.7612 15363.0094
obs Atlantic (28.15,28.2] cstar_tref ~ sal + phosphate 4434.67322 4453.19595 4.4863103 2000-2009 –> 2010-2019 10.058886 8288.3851 8324.5813
obs Atlantic (28.15,28.2] cstar_tref ~ sal + phosphate_star 4303.60129 4322.12403 4.1147232 2000-2009 –> 2010-2019 9.860039 8194.7399 8230.9361
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + phosphate_star 7956.00220 7991.79964 0.9600897 2000-2009 –> 2010-2019 2.066032 14569.3034 14639.2013
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + nitrate 8454.43787 8484.26907 1.0471711 2000-2009 –> 2010-2019 2.313760 15654.9116 15713.1599
mod Atlantic (28.15,28.2] cstar_tref ~ sal + temp + phosphate 8545.01931 8574.85051 1.0637574 2000-2009 –> 2010-2019 2.303239 15651.5158 15709.7641
obs Atlantic (28.15,28.2] cstar_tref ~ sal + temp + phosphate 4306.16226 4329.31568 4.1162460 2000-2009 –> 2010-2019 9.601212 8142.4463 8187.6916
mod Atlantic (28.15,28.2] cstar_tref ~ temp + aou + phosphate_star 8769.26625 8799.09745 1.1059582 2000-2009 –> 2010-2019 2.272374 15611.8329 15670.0811
obs Atlantic (28.15,28.2] cstar_tref ~ temp + aou + phosphate_star 3981.40471 4004.55813 3.3225135 2000-2009 –> 2010-2019 7.419638 7460.0334 7505.2786
mod Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + silicate 8419.44870 8449.27989 1.0408337 2000-2009 –> 2010-2019 2.305189 15612.2562 15670.5044
obs Atlantic (28.15,28.2] cstar_tref ~ temp + phosphate 4316.50908 4335.03181 4.1499070 2000-2009 –> 2010-2019 9.780862 8182.9979 8219.1941
obs Atlantic (28.15,28.2] cstar_tref ~ temp + phosphate_star 4060.59995 4079.12269 3.5053152 2000-2009 –> 2010-2019 7.665947 7556.0863 7592.2825
mod Atlantic (28.15,28.2] cstar_tref ~ temp + silicate + phosphate 8640.72602 8670.55721 1.0815677 2000-2009 –> 2010-2019 2.316458 15731.1004 15789.3486
obs Atlantic (28.15,28.2] cstar_tref ~ temp + silicate + phosphate_star 4045.85017 4069.00358 3.4667992 2000-2009 –> 2010-2019 7.598431 7534.7611 7580.0064
mod Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate 8863.87155 8896.12681 0.6231225 2000-2009 –> 2010-2019 1.377264 14750.9383 14812.4810
obs Atlantic (28.2, Inf] cstar_tref ~ sal + aou + phosphate_star 9462.43888 9489.86072 3.4426479 2000-2009 –> 2010-2019 6.924695 16549.4863 16602.8615
mod Atlantic (28.2, Inf] cstar_tref ~ sal + nitrate + phosphate_star 11658.74818 11691.00345 0.8399490 2000-2009 –> 2010-2019 1.861966 19117.2516 19178.7942
obs Atlantic (28.2, Inf] cstar_tref ~ sal + nitrate + phosphate_star 9490.30507 9517.72692 3.4697013 2000-2009 –> 2010-2019 7.041820 16645.1313 16698.5066
obs Atlantic (28.2, Inf] cstar_tref ~ sal + phosphate + phosphate_star 9330.87723 9358.29908 3.3177452 2000-2009 –> 2010-2019 6.725702 16360.8436 16414.2188
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate 8775.39376 8814.10008 0.6171282 2000-2009 –> 2010-2019 1.366971 14634.9041 14708.7553
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + phosphate 12433.15875 12471.86507 0.9122040 2000-2009 –> 2010-2019 1.864688 19529.3883 19603.2395
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + phosphate_star 12507.29452 12546.00084 0.9194578 2000-2009 –> 2010-2019 1.871366 19600.3950 19674.2462
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + phosphate_star 9286.01551 9318.92172 3.2743586 2000-2009 –> 2010-2019 6.649888 16292.6081 16356.6583
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + nitrate + phosphate_star 11239.86270 11278.56903 0.8030161 2000-2009 –> 2010-2019 1.628101 17593.7444 17667.5956
obs Atlantic (28.2, Inf] cstar_tref ~ sal + temp + nitrate + phosphate_star 9378.08367 9410.98988 3.3601440 2000-2009 –> 2010-2019 6.891356 16504.3421 16568.3924
mod Atlantic (28.2, Inf] cstar_tref ~ sal + temp + phosphate + phosphate_star 12481.86136 12520.56768 0.9169628 2000-2009 –> 2010-2019 1.869098 19576.1926 19650.0438
mod Atlantic (28.2, Inf] cstar_tref ~ temp + aou + phosphate 12659.67690 12691.93217 0.9347489 2000-2009 –> 2010-2019 1.925131 19955.6237 20017.1664
obs Atlantic (28.2, Inf] cstar_tref ~ temp + aou + phosphate 9548.62427 9576.04611 3.5270094 2000-2009 –> 2010-2019 7.062583 16676.1592 16729.5344
obs Atlantic (28.2, Inf] cstar_tref ~ temp + aou + phosphate_star 9615.61539 9643.03723 3.5940081 2000-2009 –> 2010-2019 7.055075 16686.6235 16739.9987
obs Atlantic (28.2, Inf] cstar_tref ~ temp + aou + silicate + phosphate 9541.41427 9574.32048 3.5178965 2000-2009 –> 2010-2019 7.052984 16670.5841 16734.6344
obs Atlantic (28.2, Inf] cstar_tref ~ temp + aou + silicate + phosphate_star 9593.84906 9626.75527 3.5700946 2000-2009 –> 2010-2019 7.029312 16665.4385 16729.4887
mod Atlantic (28.2, Inf] cstar_tref ~ temp + nitrate + phosphate_star 12460.59203 12492.84730 0.9150770 2000-2009 –> 2010-2019 1.937922 19923.2865 19984.8291
mod Atlantic (28.2, Inf] cstar_tref ~ temp + phosphate + phosphate_star 12734.26082 12766.51609 0.9422271 2000-2009 –> 2010-2019 1.936505 20050.5065 20112.0491
obs Atlantic (28.2, Inf] cstar_tref ~ temp + phosphate + phosphate_star 9596.99445 9624.41629 3.5752583 2000-2009 –> 2010-2019 7.069066 16692.9904 16746.3656
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate 19161.75120 19198.42817 4.2646832 2000-2009 –> 2010-2019 10.902085 35566.4636 35638.0246
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate 25487.83507 25525.90352 4.9934803 2000-2009 –> 2010-2019 11.450673 46367.2172 46441.6678
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 19113.25814 19156.04793 4.2325399 2000-2009 –> 2010-2019 10.048917 34866.3575 34949.8453
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 25410.11656 25454.52976 4.9464041 2000-2009 –> 2010-2019 11.248333 46136.8499 46223.7090
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 19125.52487 19168.31466 4.2403265 2000-2009 –> 2010-2019 10.159485 34965.3327 35048.8204
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 25432.58601 25476.99921 4.9596279 2000-2009 –> 2010-2019 11.271446 46169.2825 46256.1415
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 18649.14938 18691.93917 3.9482106 2000-2009 –> 2010-2019 10.138990 34711.0467 34794.5345
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 25486.77540 25531.18860 4.9916652 2000-2009 –> 2010-2019 11.445628 46364.9790 46451.8380
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + phosphate 21506.93600 21543.61297 6.0602988 2000-2009 –> 2010-2019 12.434432 37711.3102 37782.8711
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 21035.89291 21078.68270 5.6456233 2000-2009 –> 2010-2019 11.742475 37022.1124 37105.6002
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 26687.74902 26732.16222 5.7573126 2000-2009 –> 2010-2019 12.372954 47723.1653 47810.0243
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 21156.78270 21199.57249 5.7488173 2000-2009 –> 2010-2019 12.022239 37284.3223 37367.8101
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 26776.69299 26821.10619 5.8184811 2000-2009 –> 2010-2019 12.469549 47846.0439 47932.9030
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 19286.01827 19328.80807 4.3435319 2000-2009 –> 2010-2019 10.354688 35202.1687 35285.6565
mod Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 21245.19456 21287.98436 5.8254797 2000-2009 –> 2010-2019 12.135312 37401.3800 37484.8677
obs Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 26832.03990 26876.45310 5.8568718 2000-2009 –> 2010-2019 12.523343 47916.0893 48002.9483
obs Indo-Pacific (-Inf,26] cstar_tref ~ temp + aou + nitrate + silicate 26258.53734 26296.60579 5.4723545 2000-2009 –> 2010-2019 12.614060 47778.1289 47852.5796
obs Indo-Pacific (-Inf,26] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 26254.44552 26298.85872 5.4683948 2000-2009 –> 2010-2019 12.589568 47757.7435 47844.6025
obs Indo-Pacific (-Inf,26] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 26320.88653 26358.95499 5.5130465 2000-2009 –> 2010-2019 12.599409 47791.0478 47865.4984
mod Indo-Pacific (-Inf,26] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 20354.56982 20397.35961 5.0977254 2000-2009 –> 2010-2019 12.497066 37299.2005 37382.6883
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 12723.96327 12763.91900 4.2036361 2000-2009 –> 2010-2019 8.644646 26943.7062 27024.2631
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 12553.00786 12592.78524 4.3502520 2000-2009 –> 2010-2019 8.986445 27130.3755 27210.8253
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + phosphate_star 12855.34403 12889.43893 4.6662574 2000-2009 –> 2010-2019 9.589741 27727.2371 27796.1941
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 11126.59148 11166.54721 2.9363027 2000-2009 –> 2010-2019 7.184252 25129.3519 25209.9087
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 12631.64030 12671.41768 4.4297886 2000-2009 –> 2010-2019 9.177516 27326.2547 27406.7046
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 12790.81965 12830.77539 4.2672392 2000-2009 –> 2010-2019 8.749457 27055.6533 27136.2102
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 12560.54721 12600.32459 4.3578158 2000-2009 –> 2010-2019 8.999447 27143.6967 27224.1466
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate_star 12674.55994 12708.65484 4.4758756 2000-2009 –> 2010-2019 9.295041 27440.8319 27509.7889
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 12416.89638 12456.85211 3.9234725 2000-2009 –> 2010-2019 8.726295 27019.0067 27099.5636
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 12582.61771 12622.39509 4.3800333 2000-2009 –> 2010-2019 9.135622 27285.3913 27365.8411
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + phosphate + phosphate_star 12852.35867 12886.45357 4.6630487 2000-2009 –> 2010-2019 9.586196 27723.9152 27792.8722
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + phosphate_star 12951.29211 12979.70452 4.7727660 2000-2009 –> 2010-2019 9.703064 27828.0073 27885.4715
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 11186.91831 11226.87405 2.9763619 2000-2009 –> 2010-2019 7.256605 25226.6520 25307.2089
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 12632.43504 12672.21242 4.4305998 2000-2009 –> 2010-2019 9.180263 27329.0596 27409.5094
mod Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate_star 12877.35630 12911.60407 4.3529509 2000-2009 –> 2010-2019 9.155890 27477.5856 27546.6344
obs Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate_star 12632.23780 12666.33270 4.4324406 2000-2009 –> 2010-2019 9.205311 27350.9014 27419.8584
mod Indo-Pacific (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 11721.90646 11761.86219 3.3564026 2000-2009 –> 2010-2019 7.644718 25770.8387 25851.3956
mod Indo-Pacific (26,26.5] cstar_tref ~ temp + aou + silicate + phosphate_star 11916.65042 11950.89819 3.5080559 2000-2009 –> 2010-2019 8.217599 26421.0108 26490.0596
mod Indo-Pacific (26,26.5] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 11762.13148 11802.08722 3.3868660 2000-2009 –> 2010-2019 7.692737 25831.0093 25911.5662
mod Indo-Pacific (26,26.5] cstar_tref ~ temp + silicate + phosphate + phosphate_star 11945.71729 11979.96506 3.5310347 2000-2009 –> 2010-2019 8.261530 26471.7490 26540.7978
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate_star 11803.12704 11836.71743 4.6470432 2000-2009 –> 2010-2019 9.423296 25057.8478 25125.6698
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate_star 11115.33684 11148.47124 4.8723984 2000-2009 –> 2010-2019 10.073469 27540.0832 27608.5678
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 10566.93577 10606.12456 3.4072568 2000-2009 –> 2010-2019 8.020902 23669.8637 23748.9893
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 11110.36773 11149.02453 4.8632247 2000-2009 –> 2010-2019 9.996038 27466.4377 27546.3364
mod Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate_star 10930.23980 10963.83020 3.7339342 2000-2009 –> 2010-2019 8.366656 24049.4889 24117.3109
obs Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate_star 11311.43484 11344.56924 5.1377458 2000-2009 –> 2010-2019 10.739656 28133.3816 28201.8662
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + nitrate + phosphate_star 11979.52258 12013.11298 4.8570952 2000-2009 –> 2010-2019 9.821456 25405.7532 25473.5752
obs Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + nitrate + phosphate_star 11249.50024 11282.63464 5.0524147 2000-2009 –> 2010-2019 10.510098 27931.9016 28000.3863
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 10411.17781 10450.36660 3.2768102 2000-2009 –> 2010-2019 7.804291 23430.4004 23509.5261
obs Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 11196.97805 11235.63485 4.9784699 2000-2009 –> 2010-2019 10.414850 27860.4563 27940.3550
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + silicate + phosphate_star 10411.56708 10445.15748 3.2787730 2000-2009 –> 2010-2019 8.148801 23752.6169 23820.4388
obs Indo-Pacific (26.5,26.75] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 11251.69097 11284.82537 5.0554087 2000-2009 –> 2010-2019 10.518518 27939.4086 28007.8933
obs Indo-Pacific (26.5,26.75] cstar_tref ~ temp + nitrate + phosphate_star 11262.92341 11290.53541 5.0735307 2000-2009 –> 2010-2019 10.586320 27997.0731 28054.1436
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 10427.38889 10466.57768 3.2901507 2000-2009 –> 2010-2019 7.834393 23463.0185 23542.1441
obs Indo-Pacific (26.5,26.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 11201.86714 11240.52394 4.9850562 2000-2009 –> 2010-2019 10.424788 27868.6425 27948.5413
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + nitrate + silicate + phosphate_star 10807.59559 10841.18598 3.6209069 2000-2009 –> 2010-2019 8.623419 24267.8173 24335.6393
obs Indo-Pacific (26.5,26.75] cstar_tref ~ temp + nitrate + silicate + phosphate_star 11243.33973 11276.47413 5.0440049 2000-2009 –> 2010-2019 10.511781 27935.6262 28004.1109
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + silicate + phosphate + phosphate_star 10428.65116 10462.24156 3.2928419 2000-2009 –> 2010-2019 8.172236 23778.2314 23846.0534
obs Indo-Pacific (26.5,26.75] cstar_tref ~ temp + silicate + phosphate + phosphate_star 11356.59075 11389.72515 5.2008669 2000-2009 –> 2010-2019 10.914907 28284.5679 28353.0525
mod Indo-Pacific (26.5,26.75] cstar_tref ~ temp + silicate + phosphate_star 11210.65936 11238.65135 4.0078072 2000-2009 –> 2010-2019 9.037154 24692.6342 24749.1525
obs Indo-Pacific (26.75,27] cstar_tref ~ aou + silicate + phosphate 18996.90955 19027.11178 5.1526948 2000-2009 –> 2010-2019 10.466451 43893.8077 43955.5151
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + silicate + phosphate + phosphate_star 19049.89358 19086.13626 5.1951862 2000-2009 –> 2010-2019 10.614167 44106.7592 44180.8081
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + phosphate_star 13633.21638 13668.15105 3.7049528 2000-2009 –> 2010-2019 8.075701 32738.9888 32810.5318
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + phosphate_star 18652.34932 18688.59200 4.8729283 2000-2009 –> 2010-2019 10.198137 43568.5904 43642.6393
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 13151.06369 13191.82080 3.3625002 2000-2009 –> 2010-2019 7.719261 32237.6867 32321.1535
obs Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 18583.84743 18626.13055 4.8179014 2000-2009 –> 2010-2019 9.894450 43116.8481 43203.2384
mod Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + silicate + phosphate_star 13492.46959 13527.40426 3.6019525 2000-2009 –> 2010-2019 8.607506 33493.0218 33564.5647
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + phosphate_star 14198.44266 14233.37733 4.1491228 2000-2009 –> 2010-2019 8.664461 33518.9534 33590.4963
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + phosphate_star 18580.91494 18617.15762 4.8171778 2000-2009 –> 2010-2019 10.075831 43395.8359 43469.8848
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 13292.07558 13332.83270 3.4588369 2000-2009 –> 2010-2019 7.971766 32611.0651 32694.5318
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + phosphate_star 18976.30164 19006.50388 5.1356184 2000-2009 –> 2010-2019 10.624969 44135.1071 44196.8145
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + silicate + phosphate_star 13725.98514 13760.91981 3.7744475 2000-2009 –> 2010-2019 8.647154 33549.0609 33620.6038
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + silicate + phosphate_star 18922.63314 18958.87582 5.0897720 2000-2009 –> 2010-2019 10.334626 43716.3864 43790.4352
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 14202.26716 14237.20183 4.1523028 2000-2009 –> 2010-2019 8.655494 33505.0041 33576.5470
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + phosphate_star 18654.97183 18685.17407 4.8765581 2000-2009 –> 2010-2019 10.206894 43576.9665 43638.6739
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 13302.19806 13342.95517 3.4658576 2000-2009 –> 2010-2019 7.967657 32604.8952 32688.3620
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate_star 13461.08011 13496.01477 3.5793746 2000-2009 –> 2010-2019 8.163528 32881.3865 32952.9294
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate_star 18599.47445 18635.71713 4.8316008 2000-2009 –> 2010-2019 9.947087 43192.0308 43266.0796
obs Indo-Pacific (26.75,27] cstar_tref ~ temp + phosphate + phosphate_star 18957.80906 18988.01129 5.1203430 2000-2009 –> 2010-2019 10.594837 44094.7816 44156.4890
mod Indo-Pacific (26.75,27] cstar_tref ~ temp + silicate + phosphate + phosphate_star 13725.77996 13760.71463 3.7742923 2000-2009 –> 2010-2019 8.618077 33509.5769 33581.1198
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate_star 25282.80953 25321.05777 4.4599200 2000-2009 –> 2010-2019 9.328728 57582.0386 57659.8272
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + silicate + phosphate + phosphate_star 25220.42177 25258.67001 4.4279498 2000-2009 –> 2010-2019 9.258790 57435.4435 57513.2322
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + phosphate_star 25329.94483 25368.19307 4.4842272 2000-2009 –> 2010-2019 9.581579 58122.5731 58200.3618
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 25110.29900 25154.92195 4.3710681 2000-2009 –> 2010-2019 9.173072 57262.9248 57353.6783
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + nitrate + phosphate_star 16995.43384 17032.55096 2.5748300 2000-2009 –> 2010-2019 6.031537 39648.3309 39723.5859
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + nitrate + phosphate_star 25306.41564 25344.66389 4.4720769 2000-2009 –> 2010-2019 9.568883 58097.8912 58175.6799
mod Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 16874.16706 16917.47036 2.5310146 2000-2009 –> 2010-2019 5.772601 38982.0084 39069.8060
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 25157.43805 25202.06101 4.3948929 2000-2009 –> 2010-2019 9.234029 57392.9150 57483.6684
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + phosphate + phosphate_star 25328.71699 25366.96523 4.4835923 2000-2009 –> 2010-2019 9.580989 58121.4381 58199.2268
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + phosphate_star 25431.14029 25463.01383 4.5379076 2000-2009 –> 2010-2019 9.635341 58221.9393 58286.7632
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 25100.82097 25145.44392 4.3662933 2000-2009 –> 2010-2019 9.157580 57229.4129 57320.1663
obs Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + silicate + phosphate_star 25394.37686 25432.62510 4.5176686 2000-2009 –> 2010-2019 9.597670 58150.3314 58228.1200
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + phosphate_star 17117.05599 17154.17310 2.6188043 2000-2009 –> 2010-2019 5.929978 39403.7351 39478.9902
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 17113.70916 17157.01246 2.6168554 2000-2009 –> 2010-2019 5.855213 39213.0671 39300.8646
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + silicate + phosphate_star 20448.97259 20486.08970 4.1646983 2000-2009 –> 2010-2019 8.609095 45241.6729 45316.9279
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 17122.29032 17159.40744 2.6207136 2000-2009 –> 2010-2019 5.938072 39424.8568 39500.1118
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + nitrate + phosphate_star 17199.37073 17230.30166 2.6497295 2000-2009 –> 2010-2019 6.118351 39879.5635 39942.2761
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 17118.35301 17161.65631 2.6185480 2000-2009 –> 2010-2019 5.860664 39227.5856 39315.3831
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + nitrate + silicate + phosphate_star 17152.42761 17189.54473 2.6317338 2000-2009 –> 2010-2019 5.893652 39311.5040 39386.7591
mod Indo-Pacific (27,27.25] cstar_tref ~ temp + silicate + phosphate + phosphate_star 20411.24969 20448.36681 4.1428809 2000-2009 –> 2010-2019 8.563140 45157.5838 45232.8388
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 23769.63320 23814.81149 3.0388733 2000-2009 –> 2010-2019 6.318524 52196.3968 52287.7987
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + silicate + phosphate_star 26932.69797 26971.66766 3.7950177 2000-2009 –> 2010-2019 8.276787 59788.7946 59867.5752
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 23568.42443 23613.60272 2.9744355 2000-2009 –> 2010-2019 6.189163 51777.2496 51868.6514
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + silicate + phosphate + phosphate_star 26801.84695 26840.81664 3.7445806 2000-2009 –> 2010-2019 8.144473 59450.4825 59529.2631
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 23130.89678 23176.07506 2.8389926 2000-2009 –> 2010-2019 6.093520 51473.8381 51565.2400
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + phosphate 27342.22504 27381.19472 3.9573041 2000-2009 –> 2010-2019 8.538257 60444.6184 60523.3990
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + phosphate_star 27378.02246 27416.99214 3.9718154 2000-2009 –> 2010-2019 8.623692 60653.2889 60732.0695
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 26507.15605 26552.62068 3.6326889 2000-2009 –> 2010-2019 7.890879 58789.4463 58881.3570
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 23140.33306 23185.51134 2.8418476 2000-2009 –> 2010-2019 6.102422 51503.5084 51594.9103
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + phosphate_star 23389.96389 23428.68813 2.9190490 2000-2009 –> 2010-2019 6.346869 52296.3707 52374.7152
mod Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 22779.33444 22824.51273 2.7346439 2000-2009 –> 2010-2019 5.785773 50418.8412 50510.2431
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 26619.29568 26664.76031 3.6745819 2000-2009 –> 2010-2019 7.932790 58901.6347 58993.5454
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + phosphate + phosphate_star 27379.92247 27418.89215 3.9725871 2000-2009 –> 2010-2019 8.628337 60664.5527 60743.3334
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 26502.68474 26548.14937 3.6310285 2000-2009 –> 2010-2019 7.880416 58761.6923 58853.6030
obs Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate_star 26918.49193 26957.46162 3.7895092 2000-2009 –> 2010-2019 8.239731 59695.1070 59773.8876
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + phosphate_star 23142.09409 23180.81833 2.8429864 2000-2009 –> 2010-2019 6.179614 51754.5959 51832.9404
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 23153.07809 23191.80234 2.8463146 2000-2009 –> 2010-2019 6.191648 51793.9823 51872.3267
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + nitrate + phosphate_star 23673.95832 23706.22852 3.0093428 2000-2009 –> 2010-2019 6.623727 53156.0319 53221.3190
mod Indo-Pacific (27.25,27.5] cstar_tref ~ temp + nitrate + silicate + phosphate_star 22848.38782 22887.11207 2.7554197 2000-2009 –> 2010-2019 5.978125 51082.2301 51160.5746
obs Indo-Pacific (27.25,27.5] cstar_tref ~ temp + nitrate + silicate + phosphate_star 27644.57041 27683.54009 4.0815535 2000-2009 –> 2010-2019 8.602193 60597.8354 60676.6160
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 24584.24608 24630.09100 2.6130441 2000-2009 –> 2010-2019 5.526730 56752.8021 56846.0618
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate 29914.53674 29954.28441 3.5493243 2000-2009 –> 2010-2019 7.520556 65142.3524 65222.5832
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate_star 25051.53349 25090.82913 2.7345394 2000-2009 –> 2010-2019 5.774765 57767.4392 57847.3761
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 24621.21126 24667.05617 2.6224151 2000-2009 –> 2010-2019 5.531758 56770.4887 56863.7484
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + nitrate + silicate + phosphate_star 26030.61535 26069.91099 3.0065107 2000-2009 –> 2010-2019 6.040083 58718.2133 58798.1502
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + silicate + phosphate + phosphate_star 25170.50389 25209.79953 2.7662275 2000-2009 –> 2010-2019 5.817860 57934.8007 58014.7376
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 25805.61189 25851.45680 2.9411377 2000-2009 –> 2010-2019 6.157118 59249.7389 59342.9986
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + phosphate 31097.36278 31137.11045 3.9471458 2000-2009 –> 2010-2019 7.695602 65598.5600 65678.7909
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + phosphate_star 31043.38110 31083.12877 3.9280549 2000-2009 –> 2010-2019 7.717288 65680.7535 65760.9843
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 29658.73758 29705.10986 3.4680865 2000-2009 –> 2010-2019 7.130734 63870.4704 63964.0731
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + phosphate_star 25854.73506 25894.03070 2.9557352 2000-2009 –> 2010-2019 6.228033 59521.1920 59601.1289
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 24420.72294 24466.56786 2.5719897 2000-2009 –> 2010-2019 5.577846 56991.7163 57084.9760
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 25271.40082 25317.24573 2.7928482 2000-2009 –> 2010-2019 5.755981 57657.4079 57750.6676
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 29821.37908 29867.75136 3.5191189 2000-2009 –> 2010-2019 7.171671 63998.3713 64091.9739
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + phosphate + phosphate_star 31038.48062 31078.22828 3.9263264 2000-2009 –> 2010-2019 7.720495 65692.2351 65772.4659
mod Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate 24642.14915 24681.44479 2.6282469 2000-2009 –> 2010-2019 5.664883 57342.7882 57422.7251
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 29668.37797 29714.75026 3.4710907 2000-2009 –> 2010-2019 7.132855 63877.0746 63970.6773
obs Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate_star 30058.03066 30097.77832 3.5953637 2000-2009 –> 2010-2019 7.260851 64277.5162 64357.7470
obs Indo-Pacific (27.5,27.75] cstar_tref ~ temp + aou + silicate + phosphate_star 31031.42042 31071.16809 3.9238375 2000-2009 –> 2010-2019 7.744157 65771.6251 65851.8559
obs Indo-Pacific (27.5,27.75] cstar_tref ~ temp + silicate + phosphate + phosphate_star 31031.60372 31071.35139 3.9239021 2000-2009 –> 2010-2019 7.747099 65781.2835 65861.5143
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + aou 8376.38227 8398.69894 2.0526372 2000-2009 –> 2010-2019 4.429912 19808.8493 19854.4622
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + aou + phosphate_star 8366.34755 8394.24339 2.0463355 2000-2009 –> 2010-2019 4.409354 19770.7387 19827.7548
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + silicate 8336.16132 8358.47799 2.0316520 2000-2009 –> 2010-2019 4.412190 19775.4868 19821.0997
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + silicate + phosphate_star 8305.86038 8333.75622 2.0149544 2000-2009 –> 2010-2019 4.370207 19693.7926 19750.8087
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou 8363.78967 8391.68551 2.0449986 2000-2009 –> 2010-2019 4.414592 19782.0733 19839.0893
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + phosphate_star 8250.75490 8284.22990 1.9857695 2000-2009 –> 2010-2019 4.274929 19498.3719 19566.7912
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + phosphate_star 16099.60044 16135.61259 3.5751159 2000-2009 –> 2010-2019 6.687742 32051.0906 32123.3762
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 15916.48064 15958.49482 3.4660306 2000-2009 –> 2010-2019 6.545918 31803.9901 31888.3234
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + nitrate + phosphate_star 8398.48190 8431.95690 2.0621513 2000-2009 –> 2010-2019 4.425537 19805.6501 19874.0694
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 15905.29467 15947.30885 3.4595467 2000-2009 –> 2010-2019 6.547354 31808.8290 31893.1623
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + phosphate + phosphate_star 8251.60586 8285.08087 1.9862013 2000-2009 –> 2010-2019 4.277407 19503.6886 19572.1079
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + phosphate + phosphate_star 16093.83451 16129.84666 3.5716670 2000-2009 –> 2010-2019 6.685932 32048.6092 32120.8948
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + silicate 8337.03636 8364.93220 2.0310681 2000-2009 –> 2010-2019 4.411600 19778.3493 19835.3654
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 15915.72338 15957.73755 3.4655913 2000-2009 –> 2010-2019 6.545827 31803.9379 31888.2711
mod Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate_star 8279.37791 8312.85292 2.0003447 2000-2009 –> 2010-2019 4.324300 19602.4231 19670.8423
obs Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate_star 15915.26437 15951.27652 3.4664853 2000-2009 –> 2010-2019 6.584061 31876.6684 31948.9541
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + silicate + phosphate_star 15964.71331 16000.72546 3.4952978 2000-2009 –> 2010-2019 6.575263 31850.3790 31922.6647
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate_star 15940.63027 15976.64242 3.4812355 2000-2009 –> 2010-2019 6.569316 31842.7163 31915.0019
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + silicate + phosphate + phosphate_star 15961.78219 15997.79434 3.4935832 2000-2009 –> 2010-2019 6.573919 31848.1985 31920.4842
obs Indo-Pacific (27.75,27.85] cstar_tref ~ temp + silicate + phosphate_star 15994.61771 16024.62783 3.5140145 2000-2009 –> 2010-2019 6.644425 31979.6577 32039.8958
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + aou + phosphate_star 10166.83433 10195.93853 1.8569764 2000-2009 –> 2010-2019 4.060152 23381.2703 23440.3880
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + silicate 9924.20791 9947.49128 1.7694518 2000-2009 –> 2010-2019 3.930257 23020.5583 23067.8525
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + silicate + phosphate_star 9905.16826 9934.27246 1.7619979 2000-2009 –> 2010-2019 3.918737 22992.2588 23051.3765
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou 9932.90086 9962.00507 1.7718296 2000-2009 –> 2010-2019 3.896054 22929.1821 22988.2997
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + phosphate 22759.07202 22797.31335 3.3438996 2000-2009 –> 2010-2019 6.108400 43986.8754 44063.3912
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + phosphate_star 9536.34910 9571.27415 1.6356601 2000-2009 –> 2010-2019 3.694339 22347.2659 22418.2071
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 22504.16117 22548.77605 3.2461774 2000-2009 –> 2010-2019 5.999316 43698.0936 43787.3620
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + nitrate + phosphate_star 9839.33047 9874.25552 1.7381777 2000-2009 –> 2010-2019 3.903704 22952.7278 23023.6690
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 22492.29044 22536.90532 3.2417318 2000-2009 –> 2010-2019 6.002560 43710.5176 43799.7860
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + phosphate 10117.51054 10146.61475 1.8386896 2000-2009 –> 2010-2019 4.036550 23317.5072 23376.6249
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + phosphate + phosphate_star 9537.35775 9572.28279 1.6359911 2000-2009 –> 2010-2019 3.697624 22356.8459 22427.7871
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + silicate 9903.95327 9933.05748 1.7615685 2000-2009 –> 2010-2019 3.906768 22958.9721 23018.0898
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 22503.61372 22548.22859 3.2459723 2000-2009 –> 2010-2019 6.000564 43702.1428 43791.4112
mod Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate_star 9660.40142 9695.32646 1.6768827 2000-2009 –> 2010-2019 3.772091 22576.4612 22647.4024
obs Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate_star 22504.45694 22542.69826 3.2470379 2000-2009 –> 2010-2019 6.034511 43804.3414 43880.8572
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + silicate + phosphate 22551.84578 22590.08710 3.2648508 2000-2009 –> 2010-2019 6.008412 43713.4262 43789.9420
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + silicate + phosphate_star 22529.01490 22567.25622 3.2562568 2000-2009 –> 2010-2019 6.009469 43721.1832 43797.6990
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + nitrate + silicate + phosphate_star 22505.88900 22544.13033 3.2475748 2000-2009 –> 2010-2019 6.008404 43722.1181 43798.6339
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + silicate + phosphate + phosphate_star 22527.96376 22566.20508 3.2558617 2000-2009 –> 2010-2019 6.010483 43724.5867 43801.1025
obs Indo-Pacific (27.85,27.95] cstar_tref ~ temp + silicate + phosphate_star 22567.90103 22599.76880 3.2716632 2000-2009 –> 2010-2019 6.075430 43916.5505 43980.3136
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + silicate 13668.79431 13693.43296 1.7062154 2000-2009 –> 2010-2019 3.807242 27734.3111 27783.2979
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + silicate + phosphate_star 13658.73834 13689.53665 1.7032768 2000-2009 –> 2010-2019 3.785759 27668.5914 27729.8248
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou 13375.74324 13406.54154 1.6357336 2000-2009 –> 2010-2019 3.638606 27132.0841 27193.3175
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate 37789.15708 37830.51358 3.2413920 2000-2009 –> 2010-2019 6.143289 75671.8749 75754.8642
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate_star 13198.58104 13235.53900 1.5943639 2000-2009 –> 2010-2019 3.425305 26373.1716 26446.6518
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate_star 37925.22012 37966.57662 3.2718290 2000-2009 –> 2010-2019 6.226263 76081.4543 76164.4436
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + nitrate 13734.60586 13765.40417 1.7218537 2000-2009 –> 2010-2019 3.879520 27975.1367 28036.3701
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + nitrate + phosphate_star 13348.01279 13384.97076 1.6287951 2000-2009 –> 2010-2019 3.541408 26806.4381 26879.9182
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + phosphate 13661.96027 13692.75857 1.7040617 2000-2009 –> 2010-2019 3.836708 27826.6310 27887.8644
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + phosphate + phosphate_star 13197.13329 13234.09126 1.5940339 2000-2009 –> 2010-2019 3.426993 26378.8893 26452.3694
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + silicate 13275.67601 13306.47432 1.6124968 2000-2009 –> 2010-2019 3.662944 27184.7012 27245.9346
mod Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate_star 12944.94333 12981.90130 1.5375801 2000-2009 –> 2010-2019 3.396145 26216.9284 26290.4085
obs Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate_star 37471.68672 37513.04321 3.1714713 2000-2009 –> 2010-2019 6.052959 75246.8153 75329.8046
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + phosphate 37790.33870 37824.80245 3.2421005 2000-2009 –> 2010-2019 6.144441 75673.3852 75742.5429
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + phosphate_star 37938.66870 37973.13245 3.2753028 2000-2009 –> 2010-2019 6.232335 76106.3019 76175.4597
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + silicate + phosphate 37597.89184 37639.24834 3.1990847 2000-2009 –> 2010-2019 6.065333 75292.1843 75375.1736
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + silicate + phosphate_star 37633.41831 37674.77480 3.2069010 2000-2009 –> 2010-2019 6.108996 75517.1783 75600.1676
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + nitrate + silicate + phosphate_star 37553.87768 37595.23417 3.1894273 2000-2009 –> 2010-2019 6.065757 75301.6944 75384.6836
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + silicate + phosphate + phosphate_star 37638.97284 37680.32933 3.2081249 2000-2009 –> 2010-2019 6.114623 75545.8438 75628.8330
obs Indo-Pacific (27.95,28.05] cstar_tref ~ temp + silicate + phosphate_star 37761.57104 37796.03479 3.2357002 2000-2009 –> 2010-2019 6.179464 75860.6498 75929.8075
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + aou + phosphate 37523.97188 37558.50656 3.0698409 2000-2009 –> 2010-2019 5.818798 65466.2314 65534.0482
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + aou + silicate + phosphate_star 37222.72297 37264.16458 3.0074386 2000-2009 –> 2010-2019 5.740188 65099.0149 65180.3951
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + silicate + phosphate + phosphate_star 37192.52248 37233.96409 3.0012939 2000-2009 –> 2010-2019 5.718858 65004.7685 65086.1487
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou 7342.98810 7370.96001 1.5315841 2000-2009 –> 2010-2019 3.474434 15364.5689 15420.3490
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate 37514.87994 37556.32155 3.0675357 2000-2009 –> 2010-2019 5.765017 65241.8666 65323.2468
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate_star 7295.37528 7328.94157 1.5125821 2000-2009 –> 2010-2019 3.232707 14850.6720 14917.6081
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 37177.21269 37225.56123 2.9977776 2000-2009 –> 2010-2019 5.652825 64723.9508 64818.8943
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + nitrate + phosphate_star 7387.13444 7420.70072 1.5479138 2000-2009 –> 2010-2019 3.298273 15009.4446 15076.3808
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 37313.09510 37361.44365 3.0254916 2000-2009 –> 2010-2019 5.687589 64890.3160 64985.2596
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + phosphate + phosphate_star 7294.99487 7328.56116 1.5124373 2000-2009 –> 2010-2019 3.233877 14853.2290 14920.1652
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate 7167.18558 7195.15748 1.4653064 2000-2009 –> 2010-2019 3.356302 15084.7228 15140.5029
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 37163.55406 37211.90261 2.9950059 2000-2009 –> 2010-2019 5.647715 64700.1644 64795.1080
mod Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate_star 6898.55671 6932.12300 1.3688412 2000-2009 –> 2010-2019 2.951705 14134.0143 14200.9504
obs Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate_star 37330.35218 37371.79380 3.0294399 2000-2009 –> 2010-2019 5.691864 64906.9826 64988.3628
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + phosphate_star 7576.55649 7604.52839 1.6242996 2000-2009 –> 2010-2019 3.358445 15161.0729 15216.8530
obs Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate 37571.11075 37612.55236 3.0792396 2000-2009 –> 2010-2019 5.758388 65219.7169 65301.0971
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + nitrate + phosphate_star 7632.11713 7660.08904 1.6471685 2000-2009 –> 2010-2019 3.401192 15260.4701 15316.2503
obs Indo-Pacific (28.05,28.1] cstar_tref ~ temp + nitrate + silicate + phosphate_star 37621.95074 37663.39235 3.0898598 2000-2009 –> 2010-2019 5.794621 65379.9185 65461.2987
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + phosphate + phosphate_star 7579.04212 7607.01402 1.6253158 2000-2009 –> 2010-2019 3.361165 15167.3372 15223.1174
mod Indo-Pacific (28.05,28.1] cstar_tref ~ temp + silicate + phosphate_star 7388.55528 7416.52719 1.5492468 2000-2009 –> 2010-2019 3.158384 14685.3268 14741.1069
obs Indo-Pacific (28.1, Inf] cstar_tref ~ sal + aou + nitrate + phosphate 31874.05068 31914.86491 2.6556896 2000-2009 –> 2010-2019 5.740481 70734.2872 70816.7413
obs Indo-Pacific (28.1, Inf] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 31867.74608 31915.36268 2.6540319 2000-2009 –> 2010-2019 5.723018 70651.5886 70747.7851
obs Indo-Pacific (28.1, Inf] cstar_tref ~ sal + aou + silicate + phosphate 32030.27265 32071.08688 2.6870673 2000-2009 –> 2010-2019 5.800742 71032.7489 71115.2030
mod Indo-Pacific (28.1, Inf] cstar_tref ~ sal + aou + silicate + phosphate_star 77754.42502 77802.22914 1.4983651 2000-2009 –> 2010-2019 3.187749 158806.0200 158901.4950
mod Indo-Pacific (28.1, Inf] cstar_tref ~ sal + nitrate + silicate + phosphate_star 78004.75473 78052.55885 1.5071880 2000-2009 –> 2010-2019 3.220066 159632.2963 159727.7714
mod Indo-Pacific (28.1, Inf] cstar_tref ~ sal + silicate + phosphate + phosphate_star 77204.91895 77252.72308 1.4791785 2000-2009 –> 2010-2019 3.143430 157631.4792 157726.9542
mod Indo-Pacific (28.1, Inf] cstar_tref ~ sal + silicate + phosphate_star 79065.57895 79105.41572 1.5452293 2000-2009 –> 2010-2019 3.299434 161685.3205 161764.8830
obs Indo-Pacific (28.1, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate 31866.87866 31914.49527 2.6538588 2000-2009 –> 2010-2019 5.634313 70203.9792 70300.1756
obs Indo-Pacific (28.1, Inf] cstar_tref ~ sal + temp + aou + phosphate 32097.52400 32138.33823 2.7006889 2000-2009 –> 2010-2019 5.710494 70582.1847 70664.6387
mod Indo-Pacific (28.1, Inf] cstar_tref ~ sal + temp + aou + silicate 80040.87449 80088.67861 1.5809079 2000-2009 –> 2010-2019 3.405870 164311.6695 164407.1446
obs Indo-Pacific (28.1, Inf] cstar_tref ~ sal + temp + aou + silicate + phosphate 32021.44222 32069.05882 2.6848801 2000-2009 –> 2010-2019 5.691975 70494.3525 70590.5490
mod Indo-Pacific (28.1, Inf] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 74596.90535 74652.67683 1.3913487 2000-2009 –> 2010-2019 2.923689 151581.7254 151693.1129
obs Indo-Pacific (28.1, Inf] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 32179.71534 32227.33195 2.7170217 2000-2009 –> 2010-2019 5.807288 71069.0186 71165.2151
obs Indo-Pacific (28.1, Inf] cstar_tref ~ sal + temp + nitrate + phosphate_star 32444.38351 32485.19774 2.7720484 2000-2009 –> 2010-2019 5.880879 71423.0965 71505.5506
mod Indo-Pacific (28.1, Inf] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 74867.60771 74923.37918 1.4002103 2000-2009 –> 2010-2019 2.939807 152049.4540 152160.8416
obs Indo-Pacific (28.1, Inf] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 32397.18296 32444.79957 2.7618127 2000-2009 –> 2010-2019 5.863738 71343.9556 71440.1521
mod Indo-Pacific (28.1, Inf] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 74596.07084 74651.84232 1.3913215 2000-2009 –> 2010-2019 2.923644 151580.4069 151691.7944
mod Indo-Pacific (28.1, Inf] cstar_tref ~ sal + temp + silicate + phosphate_star 75045.63125 75093.43537 1.4061346 2000-2009 –> 2010-2019 2.945953 152231.4677 152326.9428
obs Indo-Pacific (28.1, Inf] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 32688.71103 32736.32764 2.8230184 2000-2009 –> 2010-2019 5.858549 71305.2646 71401.4611
mod Indo-Pacific (28.1, Inf] cstar_tref ~ temp + nitrate + silicate + phosphate_star 80993.68223 81041.48636 1.6166332 2000-2009 –> 2010-2019 3.384223 163932.4258 164027.9009

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

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

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

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

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

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

Version Author Date
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] 4 4 6 6 7 9 4
(26,26.5] 6 5 6 7 7 6 10
(26.5,26.75] 5 6 6 8 7 5 10
(26.75,27] 6 2 6 6 6 8 5
(27,27.25] 4 2 4 9 8 9 7
(27.25,27.5] 4 3 4 8 10 5 8
(27.5,27.75] 6 2 6 6 10 5 8
(27.75,27.85] 5 4 4 6 10 10 6
(27.85,27.95] 6 4 4 5 10 9 6
(27.95,28.05] 8 NA 5 5 8 6 8
(28.05,28.1] 2 1 4 5 7 4 6
(28.1,28.15] NA 5 5 NA 7 5 6
(28.15,28.2] 4 5 3 3 5 4 6
(28.2, Inf] 5 5 4 6 7 NA 8
total — — 67.00 — 109.00 — 98.00
Indo-Pacific - mod
(-Inf,26] 7 6 6 5 9 5 10
(26,26.5] 4 5 4 10 6 8 10
(26.5,26.75] 3 6 2 10 3 8 10
(26.75,27] 3 7 3 10 3 7 10
(27,27.25] 3 8 3 10 2 6 10
(27.25,27.5] 3 10 3 10 6 4 8
(27.5,27.75] 3 7 4 8 10 8 5
(27.75,27.85] 4 1 1 6 10 4 6
(27.85,27.95] 3 1 2 6 10 4 7
(27.95,28.05] 2 2 2 5 10 4 8
(28.05,28.1] 3 2 2 8 6 3 10
(28.1, Inf] 3 3 2 9 9 10 6
total 41.00 58.00 34.00 97.00 84.00 71.00 100.00
Atlantic - obs
(-Inf,26] 7 9 6 5 6 5 10
(26,26.5] 10 6 4 4 6 5 10
(26.5,26.75] 5 3 2 7 6 4 1
(26.75,27] 3 2 5 8 6 5 1
(27,27.25] 3 3 2 10 4 6 10
(27.25,27.5] 7 NA 7 6 6 6 10
(27.5,27.75] 7 NA 8 5 6 8 8
(27.75,27.85] 4 NA 7 7 5 10 6
(27.85,27.95] 3 NA 7 7 6 10 6
(27.95,28.05] 5 1 6 7 10 5 4
(28.05,28.1] 4 2 3 8 10 2 5
(28.1,28.15] 2 2 2 7 4 2 6
(28.15,28.2] 3 1 4 6 4 1 5
(28.2, Inf] 6 2 4 8 5 2 7
total 69.00 — 67.00 95.00 84.00 71.00 89.00
Indo-Pacific - obs
(-Inf,26] 8 7 4 5 7 6 10
(26,26.5] 3 4 3 10 10 4 10
(26.5,26.75] 2 8 3 10 3 6 10
(26.75,27] 4 5 3 9 3 5 8
(27,27.25] 3 2 3 10 10 6 8
(27.25,27.5] 4 2 4 9 9 7 8
(27.5,27.75] 5 1 5 8 8 7 9
(27.75,27.85] 3 2 3 10 6 8 10
(27.85,27.95] 4 2 4 8 5 9 10
(27.95,28.05] 6 1 4 7 3 6 10
(28.05,28.1] 5 2 5 7 8 8 7
(28.1, Inf] 7 7 8 3 9 5 7
total 54.00 43.00 49.00 96.00 81.00 77.00 107.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)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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