Last updated: 2021-03-20

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: 6c5618b

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

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


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

Unstaged changes:
    Modified:   code/Workflowr_project_managment.R

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

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

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

i_gamma_slab <- unique(GLODAP_basin_era$gamma_slab)[1]
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(basin %in% unique(GLODAP$basin)[1],
  #        era %in% unique(GLODAP$era)[c(1,2)],
  #        gamma_slab %in% unique(GLODAP$gamma_slab)[c(5,6)]) %>%
  filter_all(any_vars(is.na(.)))

GLODAP_nested <- GLODAP %>% 
  # filter(basin %in% unique(GLODAP$basin)[1],
  #        era %in% unique(GLODAP$era)[c(1,2)],
  #        gamma_slab %in% unique(GLODAP$gamma_slab)[c(5,6)]) %>%
  group_by(gamma_slab, era, basin) %>% 
  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)
  )

# 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(era, basin, gamma_slab, model, n) %>% 
  summarise(vif_max = max(VIF)) %>% 
  ungroup()

# calculate RMSE from augmented output
GLODAP_glanced_rmse <- GLODAP_augmented %>%
  group_by(gamma_slab, era, basin, 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:

  • 10

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(era, basin, gamma_slab) %>% 
  count() %>% 
  filter(n < params_local$MLR_number)
# A tibble: 2 x 4
# Groups:   era, basin, gamma_slab [2]
  era       basin    gamma_slab       n
  <chr>     <chr>    <chr>        <int>
1 1982-1999 Atlantic (28.15,28.2]     6
2 2011-2019 Atlantic (28.15,28.2]     6

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(basin, gamma_slab, model, era, AIC, BIC, rmse) %>% 
  arrange(era) %>% 
  group_by(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(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(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(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")
basin gamma_slab model AIC BIC rmse eras rmse_sum aic_sum bic_sum
Atlantic (-Inf,26] cstar_tref ~ aou + nitrate + silicate 292.2512 301.1721 5.980242 1982-1999 –> 2000-2010 13.760946 906.1330 927.3835
Atlantic (-Inf,26] cstar_tref ~ sal + aou + nitrate 295.1746 304.0956 6.182248 1982-1999 –> 2000-2010 13.714746 903.4154 924.6658
Atlantic (-Inf,26] cstar_tref ~ sal + aou + nitrate + silicate 294.1381 304.8433 5.972564 1982-1999 –> 2000-2010 13.462331 903.3890 928.8895
Atlantic (-Inf,26] cstar_tref ~ temp + aou + nitrate + silicate 288.6759 299.3811 5.613117 1982-1999 –> 2000-2010 13.357669 903.7474 929.2480
Atlantic (-Inf,26] cstar_tref ~ temp + aou + phosphate 299.2832 308.2041 6.477729 1982-1999 –> 2000-2010 12.921248 880.3535 901.6040
Atlantic (-Inf,26] cstar_tref ~ temp + aou + phosphate_star 299.5667 308.4876 6.498632 1982-1999 –> 2000-2010 12.900270 879.5023 900.7528
Atlantic (-Inf,26] cstar_tref ~ temp + aou + silicate + phosphate 299.9221 310.6272 6.378312 1982-1999 –> 2000-2010 12.455355 872.8036 898.3041
Atlantic (-Inf,26] cstar_tref ~ temp + aou + silicate + phosphate_star 300.2683 310.9734 6.403455 1982-1999 –> 2000-2010 12.458091 872.5070 898.0076
Atlantic (-Inf,26] cstar_tref ~ temp + phosphate + phosphate_star 299.7412 308.6622 6.511537 1982-1999 –> 2000-2010 12.872127 878.5576 899.8081
Atlantic (-Inf,26] cstar_tref ~ temp + silicate + phosphate + phosphate_star 300.3981 311.1032 6.412904 1982-1999 –> 2000-2010 12.438192 871.7913 897.2919
Atlantic (-Inf,26] cstar_tref ~ aou + nitrate + phosphate 438.6515 450.4987 3.647512 1982-1999 –> 2011-2019 10.467536 1029.6029 1053.7797
Atlantic (-Inf,26] cstar_tref ~ aou + nitrate + phosphate_star 448.7454 460.5926 3.888141 1982-1999 –> 2011-2019 11.069289 1048.6746 1072.8514
Atlantic (-Inf,26] cstar_tref ~ sal + nitrate + phosphate 492.9752 504.8224 5.144182 1982-1999 –> 2011-2019 15.523998 1157.0067 1181.1834
Atlantic (-Inf,26] cstar_tref ~ sal + nitrate + silicate + phosphate 486.4207 500.6374 4.873071 1982-1999 –> 2011-2019 15.011644 1148.3605 1177.3726
Atlantic (-Inf,26] cstar_tref ~ temp + aou + nitrate + phosphate 422.8677 437.0844 3.259227 1982-1999 –> 2011-2019 9.380887 997.0220 1026.0341
Atlantic (-Inf,26] cstar_tref ~ aou + silicate + phosphate 449.9971 461.8443 3.919065 2000-2010 –> 2011-2019 10.882822 755.6470 776.4152
Atlantic (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate 447.0213 461.2380 3.797568 2000-2010 –> 2011-2019 10.542134 751.8568 776.7786
Atlantic (-Inf,26] cstar_tref ~ temp + aou + nitrate 472.4565 484.3037 4.517691 2000-2010 –> 2011-2019 10.142419 759.3142 780.0824
Atlantic (-Inf,26] cstar_tref ~ temp + aou + phosphate 447.1772 459.0244 3.849741 2000-2010 –> 2011-2019 10.327470 746.4604 767.2285
Atlantic (-Inf,26] cstar_tref ~ temp + aou + phosphate_star 445.9748 457.8220 3.820554 2000-2010 –> 2011-2019 10.319185 745.5414 766.3096
Atlantic (-Inf,26] cstar_tref ~ temp + aou + silicate + phosphate 438.9488 453.1655 3.608417 2000-2010 –> 2011-2019 9.986729 738.8709 763.7927
Atlantic (-Inf,26] cstar_tref ~ temp + aou + silicate + phosphate_star 438.2512 452.4679 3.592521 2000-2010 –> 2011-2019 9.995976 738.5196 763.4414
Atlantic (-Inf,26] cstar_tref ~ temp + nitrate + silicate + phosphate_star 458.7462 472.9629 4.090098 2000-2010 –> 2011-2019 10.792703 763.0324 787.9543
Atlantic (-Inf,26] cstar_tref ~ temp + phosphate + phosphate_star 445.2525 457.0998 3.803129 2000-2010 –> 2011-2019 10.314667 744.9938 765.7619
Atlantic (-Inf,26] cstar_tref ~ temp + silicate + phosphate + phosphate_star 437.5548 451.7715 3.576720 2000-2010 –> 2011-2019 9.989624 737.9528 762.8747
Atlantic (26,26.5] cstar_tref ~ aou + phosphate 3408.9885 3425.8709 7.111874 1982-1999 –> 2000-2010 12.122435 6598.9993 6632.9352
Atlantic (26,26.5] cstar_tref ~ aou + phosphate_star 3459.1532 3476.0355 7.475502 1982-1999 –> 2000-2010 12.754517 6703.9653 6737.9012
Atlantic (26,26.5] cstar_tref ~ aou + silicate + phosphate 3409.5531 3430.6560 7.101733 1982-1999 –> 2000-2010 12.093442 6597.6056 6640.0256
Atlantic (26,26.5] cstar_tref ~ aou + silicate + phosphate_star 3461.0422 3482.1452 7.474677 1982-1999 –> 2000-2010 12.717450 6700.6206 6743.0406
Atlantic (26,26.5] cstar_tref ~ nitrate + phosphate_star 3502.0430 3518.9254 7.801104 1982-1999 –> 2000-2010 13.166226 6763.8436 6797.7795
Atlantic (26,26.5] cstar_tref ~ nitrate + silicate + phosphate_star 3496.0837 3517.1866 7.739626 1982-1999 –> 2000-2010 13.098838 6758.7269 6801.1468
Atlantic (26,26.5] cstar_tref ~ phosphate + phosphate_star 3482.2075 3499.0899 7.648795 1982-1999 –> 2000-2010 13.080451 6756.9493 6790.8853
Atlantic (26,26.5] cstar_tref ~ sal + aou + silicate 3574.8158 3595.9187 8.369682 1982-1999 –> 2000-2010 14.488487 6976.6364 7019.0564
Atlantic (26,26.5] cstar_tref ~ silicate + phosphate + phosphate_star 3482.7381 3503.8411 7.637631 1982-1999 –> 2000-2010 13.056820 6757.0672 6799.4871
Atlantic (26,26.5] cstar_tref ~ silicate + phosphate_star 3496.5761 3513.4585 7.758825 1982-1999 –> 2000-2010 13.654008 6857.3039 6891.2399
Atlantic (26,26.5] cstar_tref ~ aou + nitrate 2023.3475 2038.3705 5.870281 2000-2010 –> 2011-2019 12.643137 5383.2002 5415.1055
Atlantic (26,26.5] cstar_tref ~ aou + phosphate 1883.1589 1898.1819 4.702457 2000-2010 –> 2011-2019 11.814331 5292.1475 5324.0528
Atlantic (26,26.5] cstar_tref ~ aou + phosphate_star 1941.5008 1956.5238 5.157223 2000-2010 –> 2011-2019 12.632725 5400.6540 5432.5593
Atlantic (26,26.5] cstar_tref ~ aou + silicate + phosphate 1863.9071 1882.6858 4.546960 2000-2010 –> 2011-2019 11.648693 5273.4602 5313.3419
Atlantic (26,26.5] cstar_tref ~ aou + silicate + phosphate_star 1909.2256 1928.0043 4.884981 2000-2010 –> 2011-2019 12.359658 5370.2678 5410.1495
Atlantic (26,26.5] cstar_tref ~ nitrate + phosphate_star 1975.3447 1990.3677 5.440923 2000-2010 –> 2011-2019 13.242026 5477.3877 5509.2930
Atlantic (26,26.5] cstar_tref ~ nitrate + silicate + phosphate_star 1971.8817 1990.6604 5.394094 2000-2010 –> 2011-2019 13.133721 5467.9654 5507.8471
Atlantic (26,26.5] cstar_tref ~ phosphate + phosphate_star 1966.4535 1981.4764 5.364914 2000-2010 –> 2011-2019 13.013708 5448.6610 5480.5663
Atlantic (26,26.5] cstar_tref ~ sal + aou + silicate 1990.4374 2009.2162 5.554814 2000-2010 –> 2011-2019 13.924496 5565.2532 5605.1349
Atlantic (26,26.5] cstar_tref ~ silicate + phosphate_star 2009.7742 2024.7972 5.745550 2000-2010 –> 2011-2019 13.504376 5506.3503 5538.2556
Atlantic (26.5,26.75] cstar_tref ~ aou + phosphate 4306.1932 4323.9951 7.214976 1982-1999 –> 2000-2010 11.809993 8806.6016 8842.9525
Atlantic (26.5,26.75] cstar_tref ~ aou + phosphate_star 4395.0200 4412.8219 7.739384 1982-1999 –> 2000-2010 12.613973 8985.5586 9021.9095
Atlantic (26.5,26.75] cstar_tref ~ nitrate + phosphate_star 4401.8498 4419.6516 7.781249 1982-1999 –> 2000-2010 12.624618 8982.5837 9018.9346
Atlantic (26.5,26.75] cstar_tref ~ phosphate + phosphate_star 4403.8469 4421.6487 7.793534 1982-1999 –> 2000-2010 12.724154 9011.8265 9048.1774
Atlantic (26.5,26.75] cstar_tref ~ sal + aou 4577.5941 4595.3960 8.939999 1982-1999 –> 2000-2010 14.554103 9383.6746 9420.0255
Atlantic (26.5,26.75] cstar_tref ~ sal + nitrate 4569.0258 4586.8277 8.879697 1982-1999 –> 2000-2010 14.508524 9379.1029 9415.4538
Atlantic (26.5,26.75] cstar_tref ~ sal + phosphate 4556.3213 4574.1232 8.791034 1982-1999 –> 2000-2010 14.432589 9369.8454 9406.1963
Atlantic (26.5,26.75] cstar_tref ~ sal + silicate 4576.3753 4594.1771 8.931396 1982-1999 –> 2000-2010 14.573160 9389.9556 9426.3066
Atlantic (26.5,26.75] cstar_tref ~ silicate + phosphate_star 4385.4830 4403.2849 7.681301 1982-1999 –> 2000-2010 12.736446 9031.5235 9067.8744
Atlantic (26.5,26.75] cstar_tref ~ temp + phosphate 4591.6524 4609.4543 9.039827 1982-1999 –> 2000-2010 15.143911 9525.4226 9561.7735
Atlantic (26.5,26.75] cstar_tref ~ aou + nitrate 3098.5268 3115.5804 4.592481 2000-2010 –> 2011-2019 12.137923 7461.4177 7496.2732
Atlantic (26.5,26.75] cstar_tref ~ aou + phosphate 3012.3900 3029.4436 4.230775 2000-2010 –> 2011-2019 11.445751 7318.5832 7353.4387
Atlantic (26.5,26.75] cstar_tref ~ aou + phosphate_star 3045.1473 3062.2009 4.364845 2000-2010 –> 2011-2019 12.104230 7440.1673 7475.0228
Atlantic (26.5,26.75] cstar_tref ~ nitrate + phosphate_star 3042.5982 3059.6517 4.354261 2000-2010 –> 2011-2019 12.135510 7444.4479 7479.3034
Atlantic (26.5,26.75] cstar_tref ~ phosphate + phosphate_star 3044.0817 3061.1353 4.360418 2000-2010 –> 2011-2019 12.153951 7447.9286 7482.7840
Atlantic (26.5,26.75] cstar_tref ~ sal + aou 3184.6403 3201.6939 4.984999 2000-2010 –> 2011-2019 13.924998 7762.2344 7797.0899
Atlantic (26.5,26.75] cstar_tref ~ sal + nitrate 3144.7785 3161.8321 4.799297 2000-2010 –> 2011-2019 13.678995 7713.8043 7748.6598
Atlantic (26.5,26.75] cstar_tref ~ sal + phosphate 3144.2815 3161.3351 4.797026 2000-2010 –> 2011-2019 13.588060 7700.6028 7735.4583
Atlantic (26.5,26.75] cstar_tref ~ sal + silicate 3107.6688 3124.7224 4.632640 2000-2010 –> 2011-2019 13.564036 7684.0440 7718.8995
Atlantic (26.5,26.75] cstar_tref ~ silicate + phosphate_star 3032.1721 3049.2257 4.311239 2000-2010 –> 2011-2019 11.992540 7417.6551 7452.5106
Atlantic (26.75,27] cstar_tref ~ aou + nitrate 5362.4666 5381.0314 7.973664 1982-1999 –> 2000-2010 15.209033 11588.6601 11626.5005
Atlantic (26.75,27] cstar_tref ~ aou + phosphate 5109.3130 5127.8777 6.759171 1982-1999 –> 2000-2010 11.971573 10735.3753 10773.2157
Atlantic (26.75,27] cstar_tref ~ aou + phosphate_star 5185.4638 5204.0285 7.103638 1982-1999 –> 2000-2010 12.786062 10969.5230 11007.3634
Atlantic (26.75,27] cstar_tref ~ nitrate + phosphate_star 5185.3757 5203.9405 7.103230 1982-1999 –> 2000-2010 12.768210 10963.8087 11001.6491
Atlantic (26.75,27] cstar_tref ~ phosphate + phosphate_star 5176.9117 5195.4764 7.064094 1982-1999 –> 2000-2010 12.770937 10968.8181 11006.6585
Atlantic (26.75,27] cstar_tref ~ sal + aou 5399.8113 5418.3761 8.170422 1982-1999 –> 2000-2010 15.107946 11549.0780 11586.9184
Atlantic (26.75,27] cstar_tref ~ sal + nitrate 5361.1556 5379.7203 7.966844 1982-1999 –> 2000-2010 14.902964 11510.0521 11547.8925
Atlantic (26.75,27] cstar_tref ~ sal + phosphate 5340.8566 5359.4214 7.861979 1982-1999 –> 2000-2010 14.775816 11483.8643 11521.7047
Atlantic (26.75,27] cstar_tref ~ sal + silicate 5377.9769 5396.5416 8.054801 1982-1999 –> 2000-2010 14.977650 11523.3685 11561.2089
Atlantic (26.75,27] cstar_tref ~ silicate + phosphate_star 5194.1318 5212.6965 7.143944 1982-1999 –> 2000-2010 12.837072 10981.6350 11019.4755
Atlantic (26.75,27] cstar_tref ~ aou + nitrate 4453.0029 4471.0913 6.356297 2000-2010 –> 2011-2019 14.329961 9815.4696 9852.1227
Atlantic (26.75,27] cstar_tref ~ aou + phosphate 3966.9860 3985.0744 4.446337 2000-2010 –> 2011-2019 11.205508 9076.2990 9112.9521
Atlantic (26.75,27] cstar_tref ~ aou + phosphate_star 4057.9140 4076.0024 4.753777 2000-2010 –> 2011-2019 11.857415 9243.3778 9280.0309
Atlantic (26.75,27] cstar_tref ~ nitrate + phosphate_star 4065.5172 4083.6055 4.780428 2000-2010 –> 2011-2019 11.883657 9250.8929 9287.5460
Atlantic (26.75,27] cstar_tref ~ phosphate + phosphate_star 4051.6063 4069.6946 4.731780 2000-2010 –> 2011-2019 11.795873 9228.5180 9265.1711
Atlantic (26.75,27] cstar_tref ~ sal + aou 4338.5304 4356.6188 5.843179 2000-2010 –> 2011-2019 14.013601 9738.3418 9774.9949
Atlantic (26.75,27] cstar_tref ~ sal + nitrate 4318.9258 4337.0141 5.759553 2000-2010 –> 2011-2019 13.726396 9680.0814 9716.7345
Atlantic (26.75,27] cstar_tref ~ sal + phosphate 4290.1636 4308.2519 5.639025 2000-2010 –> 2011-2019 13.501004 9631.0202 9667.6733
Atlantic (26.75,27] cstar_tref ~ sal + silicate 4265.2889 4283.3773 5.536823 2000-2010 –> 2011-2019 13.591624 9643.2658 9679.9189
Atlantic (26.75,27] cstar_tref ~ silicate + phosphate_star 4032.0811 4050.1694 4.664332 2000-2010 –> 2011-2019 11.808276 9226.2128 9262.8659
Atlantic (27,27.25] cstar_tref ~ aou + nitrate 11140.2878 11161.6796 8.716215 1982-1999 –> 2000-2010 18.142744 24640.8146 24684.2809
Atlantic (27,27.25] cstar_tref ~ aou + phosphate 10082.4650 10103.8568 6.200383 1982-1999 –> 2000-2010 12.930511 22341.7272 22385.1934
Atlantic (27,27.25] cstar_tref ~ aou + phosphate_star 10345.1147 10366.5065 6.747507 1982-1999 –> 2000-2010 14.086861 22923.6206 22967.0868
Atlantic (27,27.25] cstar_tref ~ nitrate + phosphate_star 10413.1202 10434.5119 6.896872 1982-1999 –> 2000-2010 14.362372 23054.4072 23097.8734
Atlantic (27,27.25] cstar_tref ~ phosphate + phosphate_star 10270.3772 10291.7690 6.587084 1982-1999 –> 2000-2010 13.841727 22806.1145 22849.5807
Atlantic (27,27.25] cstar_tref ~ sal + nitrate 11087.4700 11108.8618 8.569248 1982-1999 –> 2000-2010 17.649907 24450.2842 24493.7504
Atlantic (27,27.25] cstar_tref ~ sal + phosphate 10896.8382 10918.2300 8.059122 1982-1999 –> 2000-2010 16.827765 24130.8425 24174.3087
Atlantic (27,27.25] cstar_tref ~ sal + silicate 11186.9484 11208.3402 8.848144 1982-1999 –> 2000-2010 17.765096 24482.7409 24526.2071
Atlantic (27,27.25] cstar_tref ~ silicate + phosphate_star 10574.8220 10596.2137 7.265441 1982-1999 –> 2000-2010 14.729558 23215.4261 23258.8924
Atlantic (27,27.25] cstar_tref ~ temp + phosphate 11135.5852 11156.9770 8.703028 1982-1999 –> 2000-2010 18.338647 24716.9335 24760.3997
Atlantic (27,27.25] cstar_tref ~ aou + phosphate 8077.7009 8098.2659 5.904516 2000-2010 –> 2011-2019 12.104898 18160.1659 18202.1227
Atlantic (27,27.25] cstar_tref ~ aou + phosphate_star 8281.5169 8302.0819 6.400683 2000-2010 –> 2011-2019 13.148190 18626.6317 18668.5884
Atlantic (27,27.25] cstar_tref ~ nitrate + phosphate_star 8393.0132 8413.5782 6.689533 2000-2010 –> 2011-2019 13.586405 18806.1334 18848.0902
Atlantic (27,27.25] cstar_tref ~ nitrate + silicate + phosphate_star 8275.4017 8301.1079 6.380153 2000-2010 –> 2011-2019 13.215570 18662.7221 18715.1680
Atlantic (27,27.25] cstar_tref ~ phosphate + phosphate_star 8232.5379 8253.1029 6.277769 2000-2010 –> 2011-2019 12.864854 18502.9152 18544.8719
Atlantic (27,27.25] cstar_tref ~ sal + aou + silicate 8783.1728 8808.8790 7.800672 2000-2010 –> 2011-2019 16.601955 19955.6276 20008.0735
Atlantic (27,27.25] cstar_tref ~ sal + nitrate + silicate 8772.7067 8798.4130 7.768418 2000-2010 –> 2011-2019 16.242754 19827.5832 19880.0291
Atlantic (27,27.25] cstar_tref ~ sal + phosphate 8763.1444 8783.7094 7.745195 2000-2010 –> 2011-2019 15.804317 19659.9826 19701.9394
Atlantic (27,27.25] cstar_tref ~ sal + silicate 8781.3234 8801.8884 7.801137 2000-2010 –> 2011-2019 16.649281 19968.2718 20010.2286
Atlantic (27,27.25] cstar_tref ~ silicate + phosphate_star 8315.2910 8335.8559 6.486838 2000-2010 –> 2011-2019 13.752280 18890.1129 18932.0697
Atlantic (27.25,27.5] cstar_tref ~ aou + nitrate 11422.8627 11444.5713 7.216146 1982-1999 –> 2000-2010 14.081831 26646.0801 26690.7059
Atlantic (27.25,27.5] cstar_tref ~ aou + phosphate 10033.2793 10054.9879 4.773122 1982-1999 –> 2000-2010 9.831122 23866.7890 23911.4147
Atlantic (27.25,27.5] cstar_tref ~ aou + phosphate_star 10534.1649 10555.8734 5.539947 1982-1999 –> 2000-2010 11.471321 25092.1006 25136.7263
Atlantic (27.25,27.5] cstar_tref ~ aou + silicate + phosphate_star 10466.6252 10493.7609 5.426536 1982-1999 –> 2000-2010 11.144864 24860.1981 24915.9803
Atlantic (27.25,27.5] cstar_tref ~ nitrate + phosphate_star 10607.3980 10629.1066 5.661946 1982-1999 –> 2000-2010 11.522016 25110.3297 25154.9554
Atlantic (27.25,27.5] cstar_tref ~ nitrate + silicate + phosphate_star 10577.6313 10604.7670 5.608699 1982-1999 –> 2000-2010 11.299191 24949.0107 25004.7929
Atlantic (27.25,27.5] cstar_tref ~ phosphate + phosphate_star 10317.5783 10339.2869 5.194306 1982-1999 –> 2000-2010 10.869546 24674.7516 24719.3774
Atlantic (27.25,27.5] cstar_tref ~ sal + phosphate 11194.8901 11216.5987 6.743050 1982-1999 –> 2000-2010 14.658345 27065.1109 27109.7367
Atlantic (27.25,27.5] cstar_tref ~ sal + silicate + phosphate 11175.2641 11202.3999 6.699815 1982-1999 –> 2000-2010 14.612631 27046.0606 27101.8428
Atlantic (27.25,27.5] cstar_tref ~ silicate + phosphate + phosphate_star 10239.1674 10266.3032 5.071545 1982-1999 –> 2000-2010 10.490892 24388.5080 24444.2902
Atlantic (27.25,27.5] cstar_tref ~ aou + nitrate 9195.3126 9216.2550 6.623422 2000-2010 –> 2011-2019 13.839568 20618.1752 20660.8263
Atlantic (27.25,27.5] cstar_tref ~ aou + phosphate 8316.9977 8337.9401 4.826957 2000-2010 –> 2011-2019 9.600079 18350.2770 18392.9280
Atlantic (27.25,27.5] cstar_tref ~ aou + phosphate_star 8707.6778 8728.6202 5.556404 2000-2010 –> 2011-2019 11.096351 19241.8426 19284.4937
Atlantic (27.25,27.5] cstar_tref ~ aou + silicate + phosphate_star 8704.8906 8731.0687 5.546830 2000-2010 –> 2011-2019 10.973366 19171.5158 19224.8296
Atlantic (27.25,27.5] cstar_tref ~ nitrate + phosphate_star 8751.4928 8772.4353 5.644799 2000-2010 –> 2011-2019 11.306745 19358.8908 19401.5418
Atlantic (27.25,27.5] cstar_tref ~ nitrate + silicate + phosphate_star 8752.5179 8778.6960 5.642817 2000-2010 –> 2011-2019 11.251517 19330.1492 19383.4630
Atlantic (27.25,27.5] cstar_tref ~ phosphate + phosphate_star 8517.8625 8538.8050 5.189171 2000-2010 –> 2011-2019 10.383477 18835.4408 18878.0919
Atlantic (27.25,27.5] cstar_tref ~ sal + phosphate 9166.0576 9187.0001 6.553987 2000-2010 –> 2011-2019 13.297038 20360.9477 20403.5987
Atlantic (27.25,27.5] cstar_tref ~ sal + silicate + phosphate 9090.8817 9117.0598 6.374288 2000-2010 –> 2011-2019 13.074104 20266.1459 20319.4597
Atlantic (27.25,27.5] cstar_tref ~ silicate + phosphate + phosphate_star 8515.0875 8541.2656 5.180252 2000-2010 –> 2011-2019 10.251797 18754.2550 18807.5688
Atlantic (27.5,27.75] cstar_tref ~ aou + phosphate 11699.2903 11721.5431 5.033842 1982-1999 –> 2000-2010 10.000093 24917.7416 24962.7537
Atlantic (27.5,27.75] cstar_tref ~ aou + phosphate_star 12586.7598 12609.0126 6.338077 1982-1999 –> 2000-2010 12.328948 26625.2744 26670.2865
Atlantic (27.5,27.75] cstar_tref ~ aou + silicate + phosphate_star 12497.8237 12525.6397 6.190203 1982-1999 –> 2000-2010 11.969196 26380.9141 26437.1793
Atlantic (27.5,27.75] cstar_tref ~ nitrate + phosphate_star 12531.2682 12553.5210 6.247426 1982-1999 –> 2000-2010 12.149222 26504.2900 26549.3021
Atlantic (27.5,27.75] cstar_tref ~ nitrate + silicate + phosphate_star 12381.7520 12409.5680 6.006457 1982-1999 –> 2000-2010 11.628960 26144.8200 26201.0851
Atlantic (27.5,27.75] cstar_tref ~ phosphate + phosphate_star 11974.3169 11996.5697 5.406392 1982-1999 –> 2000-2010 10.836313 25583.0101 25628.0223
Atlantic (27.5,27.75] cstar_tref ~ sal + phosphate 12316.3937 12338.6465 5.908470 1982-1999 –> 2000-2010 12.271285 26618.2513 26663.2634
Atlantic (27.5,27.75] cstar_tref ~ sal + silicate + phosphate 12294.5236 12322.3396 5.871970 1982-1999 –> 2000-2010 12.230935 26595.7355 26652.0007
Atlantic (27.5,27.75] cstar_tref ~ silicate + phosphate + phosphate_star 11726.0758 11753.8918 5.066337 1982-1999 –> 2000-2010 10.184112 25077.9260 25134.1912
Atlantic (27.5,27.75] cstar_tref ~ temp + silicate + phosphate 12426.8279 12454.6439 6.077157 1982-1999 –> 2000-2010 12.860835 27010.7062 27066.9713
Atlantic (27.5,27.75] cstar_tref ~ aou + phosphate 10061.5344 10083.2477 4.796331 2000-2010 –> 2011-2019 9.830173 21760.8247 21804.7908
Atlantic (27.5,27.75] cstar_tref ~ aou + silicate + phosphate_star 10471.0722 10498.2139 5.413664 2000-2010 –> 2011-2019 11.603867 22968.8959 23023.8536
Atlantic (27.5,27.75] cstar_tref ~ nitrate + silicate + phosphate_star 10280.5126 10307.6543 5.115693 2000-2010 –> 2011-2019 11.122151 22662.2646 22717.2223
Atlantic (27.5,27.75] cstar_tref ~ phosphate + phosphate_star 10152.4308 10174.1441 4.927617 2000-2010 –> 2011-2019 10.334009 22126.7476 22170.7138
Atlantic (27.5,27.75] cstar_tref ~ sal + aou + phosphate_star 10529.0297 10556.1714 5.507686 2000-2010 –> 2011-2019 11.512753 22909.8899 22964.8475
Atlantic (27.5,27.75] cstar_tref ~ sal + phosphate 10257.4851 10279.1984 5.083835 2000-2010 –> 2011-2019 10.992305 22573.8788 22617.8449
Atlantic (27.5,27.75] cstar_tref ~ sal + silicate + phosphate 10196.5454 10223.6870 4.989657 2000-2010 –> 2011-2019 10.861627 22491.0690 22546.0267
Atlantic (27.5,27.75] cstar_tref ~ silicate + phosphate + phosphate_star 9850.4713 9877.6129 4.502139 2000-2010 –> 2011-2019 9.568476 21576.5471 21631.5048
Atlantic (27.5,27.75] cstar_tref ~ temp + phosphate 10339.7122 10361.4256 5.209556 2000-2010 –> 2011-2019 11.371584 22817.9630 22861.9291
Atlantic (27.5,27.75] cstar_tref ~ temp + silicate + phosphate 10260.5930 10287.7347 5.085509 2000-2010 –> 2011-2019 11.162666 22687.4208 22742.3785
Atlantic (27.75,27.85] cstar_tref ~ aou + nitrate 8753.9784 8774.6804 6.868072 1982-1999 –> 2000-2010 13.408648 15652.6328 15693.1418
Atlantic (27.75,27.85] cstar_tref ~ aou + phosphate 7985.4907 8006.1927 5.118668 1982-1999 –> 2000-2010 10.450344 14457.0337 14497.5428
Atlantic (27.75,27.85] cstar_tref ~ aou + phosphate_star 8208.1926 8228.8945 5.573872 1982-1999 –> 2000-2010 11.151428 14773.9632 14814.4722
Atlantic (27.75,27.85] cstar_tref ~ nitrate + phosphate_star 8498.5422 8519.2442 6.228682 1982-1999 –> 2000-2010 12.287291 15237.2171 15277.7262
Atlantic (27.75,27.85] cstar_tref ~ phosphate + phosphate_star 8020.8608 8041.5628 5.188399 1982-1999 –> 2000-2010 10.554718 14505.9393 14546.4484
Atlantic (27.75,27.85] cstar_tref ~ sal + aou + silicate 9135.4260 9161.3034 7.941031 1982-1999 –> 2000-2010 15.024794 16202.8207 16253.4570
Atlantic (27.75,27.85] cstar_tref ~ sal + nitrate 8713.9747 8734.6767 6.763765 1982-1999 –> 2000-2010 13.107005 15548.6005 15589.1095
Atlantic (27.75,27.85] cstar_tref ~ sal + phosphate 8070.7623 8091.4643 5.288398 1982-1999 –> 2000-2010 10.684670 14567.4743 14607.9834
Atlantic (27.75,27.85] cstar_tref ~ temp + nitrate 8735.2700 8755.9720 6.819092 1982-1999 –> 2000-2010 13.213478 15586.6803 15627.1893
Atlantic (27.75,27.85] cstar_tref ~ temp + phosphate 8081.9391 8102.6411 5.311058 1982-1999 –> 2000-2010 10.725412 14585.6428 14626.1518
Atlantic (27.75,27.85] cstar_tref ~ aou + phosphate_star 6806.0786 6825.8357 6.518891 2000-2010 –> 2011-2019 12.092763 15014.2712 15054.7302
Atlantic (27.75,27.85] cstar_tref ~ nitrate + phosphate_star 6959.7457 6979.5028 7.022753 2000-2010 –> 2011-2019 13.251436 15458.2880 15498.7469
Atlantic (27.75,27.85] cstar_tref ~ phosphate + phosphate_star 6721.8419 6741.5989 6.258195 2000-2010 –> 2011-2019 11.446594 14742.7027 14783.1617
Atlantic (27.75,27.85] cstar_tref ~ sal + aou + phosphate_star 6783.4081 6808.1044 6.441436 2000-2010 –> 2011-2019 11.796580 14888.9563 14939.5300
Atlantic (27.75,27.85] cstar_tref ~ sal + aou + silicate 7134.2843 7158.9806 7.635051 2000-2010 –> 2011-2019 15.576082 16269.7103 16320.2840
Atlantic (27.75,27.85] cstar_tref ~ sal + nitrate 6959.5528 6979.3098 7.022097 2000-2010 –> 2011-2019 13.785862 15673.5274 15713.9864
Atlantic (27.75,27.85] cstar_tref ~ sal + phosphate 6751.1710 6770.9280 6.347758 2000-2010 –> 2011-2019 11.636156 14821.9333 14862.3923
Atlantic (27.75,27.85] cstar_tref ~ temp + aou + silicate 7146.7331 7171.4294 7.681240 2000-2010 –> 2011-2019 15.692809 16305.2762 16355.8499
Atlantic (27.75,27.85] cstar_tref ~ temp + nitrate 6958.6400 6978.3970 7.018992 2000-2010 –> 2011-2019 13.838084 15693.9100 15734.3690
Atlantic (27.75,27.85] cstar_tref ~ temp + phosphate 6724.0858 6743.8429 6.265003 2000-2010 –> 2011-2019 11.576061 14806.0250 14846.4839
Atlantic (27.85,27.95] cstar_tref ~ aou + phosphate_star 18632.3296 18656.2838 5.702707 1982-1999 –> 2000-2010 11.204380 32779.5120 32826.3640
Atlantic (27.85,27.95] cstar_tref ~ nitrate + phosphate_star 19564.1948 19588.1489 6.679512 1982-1999 –> 2000-2010 12.886349 34257.2091 34304.0610
Atlantic (27.85,27.95] cstar_tref ~ phosphate + phosphate_star 18315.9748 18339.9289 5.404689 1982-1999 –> 2000-2010 10.671270 32265.5031 32312.3551
Atlantic (27.85,27.95] cstar_tref ~ sal + aou 20232.0070 20255.9611 7.480867 1982-1999 –> 2000-2010 13.817839 35018.9343 35065.7862
Atlantic (27.85,27.95] cstar_tref ~ sal + nitrate 19513.4754 19537.4296 6.622280 1982-1999 –> 2000-2010 12.780826 34171.1383 34217.9903
Atlantic (27.85,27.95] cstar_tref ~ sal + phosphate 17992.2441 18016.1982 5.115839 1982-1999 –> 2000-2010 10.152855 31740.0597 31786.9116
Atlantic (27.85,27.95] cstar_tref ~ temp + aou 19995.5422 20019.4963 7.186679 1982-1999 –> 2000-2010 13.417137 34705.7485 34752.6005
Atlantic (27.85,27.95] cstar_tref ~ temp + nitrate 19574.5993 19598.5535 6.691314 1982-1999 –> 2000-2010 12.889619 34261.3880 34308.2400
Atlantic (27.85,27.95] cstar_tref ~ temp + phosphate 18162.9631 18186.9173 5.266186 1982-1999 –> 2000-2010 10.398821 31995.8915 32042.7435
Atlantic (27.85,27.95] cstar_tref ~ temp + phosphate_star 20552.1968 20576.1510 7.898504 1982-1999 –> 2000-2010 13.619642 34876.4155 34923.2675
Atlantic (27.85,27.95] cstar_tref ~ sal + temp 20567.0631 20590.1288 18.855321 1982-1999 –> 2011-2019 33.426166 39122.4352 39168.3987
Atlantic (27.85,27.95] cstar_tref ~ aou + phosphate_star 15385.1300 15408.1956 6.289792 2000-2010 –> 2011-2019 11.992499 34017.4596 34064.4794
Atlantic (27.85,27.95] cstar_tref ~ aou + silicate 15613.2056 15636.2713 6.601184 2000-2010 –> 2011-2019 13.897982 35698.3746 35745.3945
Atlantic (27.85,27.95] cstar_tref ~ nitrate + phosphate_star 15820.6559 15843.7215 6.897786 2000-2010 –> 2011-2019 13.577298 35384.8506 35431.8705
Atlantic (27.85,27.95] cstar_tref ~ phosphate + phosphate_star 15206.4640 15229.5297 6.056154 2000-2010 –> 2011-2019 11.460843 33522.4388 33569.4586
Atlantic (27.85,27.95] cstar_tref ~ sal + nitrate 15419.2286 15442.2942 6.335396 2000-2010 –> 2011-2019 12.957676 34932.7040 34979.7238
Atlantic (27.85,27.95] cstar_tref ~ sal + phosphate 14859.3821 14882.4478 5.626799 2000-2010 –> 2011-2019 10.742638 32851.6262 32898.6460
Atlantic (27.85,27.95] cstar_tref ~ temp + aou 15582.8423 15605.9080 6.558856 2000-2010 –> 2011-2019 13.745535 35578.3845 35625.4043
Atlantic (27.85,27.95] cstar_tref ~ temp + aou + phosphate_star 15101.6052 15130.4373 5.920586 2000-2010 –> 2011-2019 11.056168 33118.5516 33177.3264
Atlantic (27.85,27.95] cstar_tref ~ temp + nitrate 15598.9399 15622.0055 6.581263 2000-2010 –> 2011-2019 13.272577 35173.5392 35220.5590
Atlantic (27.85,27.95] cstar_tref ~ temp + phosphate 15141.7473 15164.8130 5.973684 2000-2010 –> 2011-2019 11.239870 33304.7104 33351.7303
Atlantic (27.95,28.05] cstar_tref ~ aou + phosphate_star 20596.9835 20621.2395 6.174083 1982-1999 –> 2000-2010 11.907784 40413.4302 40461.8801
Atlantic (27.95,28.05] cstar_tref ~ phosphate + phosphate_star 20237.5521 20261.8081 5.834627 1982-1999 –> 2000-2010 11.223655 39666.0259 39714.4758
Atlantic (27.95,28.05] cstar_tref ~ sal + nitrate 21026.2730 21050.5291 6.605490 1982-1999 –> 2000-2010 12.458712 40971.8276 41020.2775
Atlantic (27.95,28.05] cstar_tref ~ sal + phosphate 19708.1014 19732.3575 5.368299 1982-1999 –> 2000-2010 10.193059 38444.4157 38492.8656
Atlantic (27.95,28.05] cstar_tref ~ sal + phosphate_star 19988.5170 20012.7730 5.610441 1982-1999 –> 2000-2010 10.150466 38344.1663 38392.6162
Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate 20947.0769 20977.3969 6.521644 1982-1999 –> 2000-2010 12.369015 40888.3737 40948.9360
Atlantic (27.95,28.05] cstar_tref ~ sal + temp + phosphate 19595.2283 19625.5483 5.272148 1982-1999 –> 2000-2010 10.072270 38301.5026 38362.0650
Atlantic (27.95,28.05] cstar_tref ~ temp + aou + phosphate_star 19521.9876 19552.3076 5.211745 1982-1999 –> 2000-2010 10.014245 38231.3622 38291.9246
Atlantic (27.95,28.05] cstar_tref ~ temp + nitrate 21156.1322 21180.3882 6.741835 1982-1999 –> 2000-2010 12.691242 41203.6883 41252.1382
Atlantic (27.95,28.05] cstar_tref ~ temp + phosphate 20038.3026 20062.5586 5.654559 1982-1999 –> 2000-2010 10.587856 38913.8351 38962.2850
Atlantic (27.95,28.05] cstar_tref ~ sal + temp + phosphate_star 14296.2863 14324.8395 5.938131 1982-1999 –> 2011-2019 10.451749 32617.4295 32676.2252
Atlantic (27.95,28.05] cstar_tref ~ aou + phosphate_star 15458.5767 15481.4193 7.707641 2000-2010 –> 2011-2019 13.881724 36055.5602 36102.6588
Atlantic (27.95,28.05] cstar_tref ~ nitrate + phosphate_star 15388.6096 15411.4522 7.587776 2000-2010 –> 2011-2019 14.177392 36399.5901 36446.6888
Atlantic (27.95,28.05] cstar_tref ~ sal + nitrate 14899.9838 14922.8264 6.801066 2000-2010 –> 2011-2019 13.406556 35926.2568 35973.3554
Atlantic (27.95,28.05] cstar_tref ~ sal + phosphate 14488.0540 14510.8966 6.201562 2000-2010 –> 2011-2019 11.569861 34196.1555 34243.2541
Atlantic (27.95,28.05] cstar_tref ~ sal + phosphate_star 14537.4675 14560.3101 6.270590 2000-2010 –> 2011-2019 11.881031 34525.9845 34573.0832
Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate 14812.6729 14841.2261 6.666350 2000-2010 –> 2011-2019 13.187993 35759.7498 35818.6231
Atlantic (27.95,28.05] cstar_tref ~ sal + temp + phosphate 14416.5570 14445.1103 6.100293 2000-2010 –> 2011-2019 11.372441 34011.7853 34070.6586
Atlantic (27.95,28.05] cstar_tref ~ temp + nitrate 15229.6975 15252.5401 7.322413 2000-2010 –> 2011-2019 14.064248 36385.8296 36432.9283
Atlantic (27.95,28.05] cstar_tref ~ temp + phosphate 14980.8562 15003.6988 6.925401 2000-2010 –> 2011-2019 12.579960 35019.1588 35066.2574
Atlantic (27.95,28.05] cstar_tref ~ temp + phosphate_star 14964.8001 14987.6428 6.900536 2000-2010 –> 2011-2019 13.196029 35685.5583 35732.6569
Atlantic (28.05,28.1] cstar_tref ~ aou + phosphate_star 15550.8961 15574.0339 6.141827 1982-1999 –> 2000-2010 11.867891 29227.2330 29273.0823
Atlantic (28.05,28.1] cstar_tref ~ sal + nitrate 15711.9679 15735.1058 6.351157 1982-1999 –> 2000-2010 12.231056 29502.8332 29548.6826
Atlantic (28.05,28.1] cstar_tref ~ sal + phosphate 14943.4411 14966.5790 5.412587 1982-1999 –> 2000-2010 10.313662 27947.6998 27993.5491
Atlantic (28.05,28.1] cstar_tref ~ sal + phosphate_star 14172.0414 14195.1793 4.609961 1982-1999 –> 2000-2010 8.540177 26222.6156 26268.4650
Atlantic (28.05,28.1] cstar_tref ~ sal + silicate 16343.0044 16366.1423 7.242300 1982-1999 –> 2000-2010 13.369965 30312.1743 30358.0237
Atlantic (28.05,28.1] cstar_tref ~ temp + aou 16497.9771 16521.1150 7.479639 1982-1999 –> 2000-2010 14.171720 30847.7872 30893.6366
Atlantic (28.05,28.1] cstar_tref ~ temp + nitrate 15772.0936 15795.2315 6.431113 1982-1999 –> 2000-2010 12.297439 29552.9756 29598.8250
Atlantic (28.05,28.1] cstar_tref ~ temp + phosphate 15119.0582 15142.1961 5.614027 1982-1999 –> 2000-2010 10.596813 28194.7456 28240.5950
Atlantic (28.05,28.1] cstar_tref ~ temp + phosphate_star 13950.0565 13973.1944 4.401874 1982-1999 –> 2000-2010 8.331661 26000.1593 26046.0086
Atlantic (28.05,28.1] cstar_tref ~ temp + silicate 16364.9715 16388.1094 7.275478 1982-1999 –> 2000-2010 13.307414 30266.1194 30311.9688
Atlantic (28.05,28.1] cstar_tref ~ aou + phosphate_star 11419.3313 11441.0303 7.267167 2000-2010 –> 2011-2019 13.408994 26970.2273 27015.0643
Atlantic (28.05,28.1] cstar_tref ~ nitrate + phosphate_star 11356.2156 11377.9147 7.131692 2000-2010 –> 2011-2019 13.497043 27078.9127 27123.7496
Atlantic (28.05,28.1] cstar_tref ~ sal + nitrate 11139.3746 11161.0736 6.685206 2000-2010 –> 2011-2019 13.036363 26851.3424 26896.1794
Atlantic (28.05,28.1] cstar_tref ~ sal + phosphate 10867.5572 10889.2563 6.164791 2000-2010 –> 2011-2019 11.577378 25810.9983 25855.8353
Atlantic (28.05,28.1] cstar_tref ~ sal + phosphate_star 10403.4460 10425.1450 5.368126 2000-2010 –> 2011-2019 9.978087 24575.4874 24620.3243
Atlantic (28.05,28.1] cstar_tref ~ sal + silicate 11645.7232 11667.4222 7.774628 2000-2010 –> 2011-2019 15.016928 27988.7276 28033.5646
Atlantic (28.05,28.1] cstar_tref ~ silicate + phosphate_star 11696.4082 11718.1073 7.893009 2000-2010 –> 2011-2019 14.723551 27758.0942 27802.9311
Atlantic (28.05,28.1] cstar_tref ~ temp + nitrate 11202.2560 11223.9551 6.811724 2000-2010 –> 2011-2019 13.242836 26974.3497 27019.1866
Atlantic (28.05,28.1] cstar_tref ~ temp + phosphate 10999.3624 11021.0614 6.411878 2000-2010 –> 2011-2019 12.025905 26118.4206 26163.2575
Atlantic (28.05,28.1] cstar_tref ~ temp + phosphate_star 10422.5552 10444.2542 5.398798 2000-2010 –> 2011-2019 9.800672 24372.6117 24417.4487
Atlantic (28.1,28.15] cstar_tref ~ aou + phosphate_star 12058.3245 12080.3924 6.406972 1982-1999 –> 2000-2010 11.695549 21264.2669 21307.5636
Atlantic (28.1,28.15] cstar_tref ~ sal + aou 12841.2981 12863.3660 7.926940 1982-1999 –> 2000-2010 14.163038 22538.6904 22581.9871
Atlantic (28.1,28.15] cstar_tref ~ sal + nitrate 12047.4741 12069.5420 6.388099 1982-1999 –> 2000-2010 11.766761 21303.7841 21347.0808
Atlantic (28.1,28.15] cstar_tref ~ sal + phosphate 11761.6260 11783.6939 5.910429 1982-1999 –> 2000-2010 10.379009 20465.1659 20508.4626
Atlantic (28.1,28.15] cstar_tref ~ sal + silicate 12556.2442 12578.3121 7.335786 1982-1999 –> 2000-2010 12.621655 21760.6589 21803.9556
Atlantic (28.1,28.15] cstar_tref ~ temp + aou 12959.4928 12981.5607 8.185814 1982-1999 –> 2000-2010 14.419865 22655.9064 22699.2031
Atlantic (28.1,28.15] cstar_tref ~ temp + nitrate 12029.2859 12051.3538 6.356587 1982-1999 –> 2000-2010 11.679379 21254.4588 21297.7555
Atlantic (28.1,28.15] cstar_tref ~ temp + phosphate 11741.6517 11763.7196 5.878418 1982-1999 –> 2000-2010 10.403408 20482.5999 20525.8966
Atlantic (28.1,28.15] cstar_tref ~ temp + phosphate_star 10365.2864 10387.3543 4.043343 1982-1999 –> 2000-2010 7.381951 18199.5222 18242.8190
Atlantic (28.1,28.15] cstar_tref ~ temp + silicate 12549.8447 12571.9126 7.323034 1982-1999 –> 2000-2010 12.573544 21734.2453 21777.5420
Atlantic (28.1,28.15] cstar_tref ~ aou + phosphate_star 7914.3158 7934.4109 8.176283 2000-2010 –> 2011-2019 14.583255 19972.6403 20014.8033
Atlantic (28.1,28.15] cstar_tref ~ sal + aou 8072.6643 8092.7593 8.773538 2000-2010 –> 2011-2019 16.700477 20913.9623 20956.1253
Atlantic (28.1,28.15] cstar_tref ~ sal + nitrate 7604.9281 7625.0231 7.124127 2000-2010 –> 2011-2019 13.512225 19652.4022 19694.5651
Atlantic (28.1,28.15] cstar_tref ~ sal + phosphate 7317.4005 7337.4956 6.268077 2000-2010 –> 2011-2019 12.178505 19079.0265 19121.1895
Atlantic (28.1,28.15] cstar_tref ~ sal + silicate 7852.4358 7872.5309 7.954092 2000-2010 –> 2011-2019 15.289878 20408.6800 20450.8430
Atlantic (28.1,28.15] cstar_tref ~ temp + aou 8032.3317 8052.4267 8.617393 2000-2010 –> 2011-2019 16.803207 20991.8245 21033.9874
Atlantic (28.1,28.15] cstar_tref ~ temp + nitrate 7521.0360 7541.1310 6.862936 2000-2010 –> 2011-2019 13.219523 19550.3219 19592.4848
Atlantic (28.1,28.15] cstar_tref ~ temp + phosphate 7257.0426 7277.1376 6.101875 2000-2010 –> 2011-2019 11.980292 18998.6943 19040.8572
Atlantic (28.1,28.15] cstar_tref ~ temp + phosphate_star 6647.8481 6667.9431 4.652295 2000-2010 –> 2011-2019 8.695638 17013.1344 17055.2974
Atlantic (28.1,28.15] cstar_tref ~ temp + silicate 7739.7571 7759.8522 7.564891 2000-2010 –> 2011-2019 14.887924 20289.6018 20331.7648
Atlantic (28.15,28.2] cstar_tref ~ sal + temp 6284.1044 6303.3318 7.786241 1982-1999 –> 2000-2010 14.600475 9816.9821 9853.2858
Atlantic (28.15,28.2] cstar_tref ~ temp + aou 6358.2438 6377.4711 8.112162 1982-1999 –> 2000-2010 14.193949 9771.0381 9807.3418
Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate 5698.1013 5717.3286 5.630743 1982-1999 –> 2000-2010 10.691956 8916.9164 8953.2201
Atlantic (28.15,28.2] cstar_tref ~ temp + phosphate 5496.2576 5515.4850 5.035950 1982-1999 –> 2000-2010 9.235796 8518.0672 8554.3709
Atlantic (28.15,28.2] cstar_tref ~ temp + phosphate_star 5008.3677 5027.5950 3.844915 1982-1999 –> 2000-2010 7.328979 7832.8652 7869.1689
Atlantic (28.15,28.2] cstar_tref ~ temp + silicate 6020.2127 6039.4400 6.728824 1982-1999 –> 2000-2010 11.979622 9277.8612 9314.1649
Atlantic (28.15,28.2] cstar_tref ~ sal + temp 3007.0652 3023.3757 7.541012 2000-2010 –> 2011-2019 15.327253 9291.1696 9326.7075
Atlantic (28.15,28.2] cstar_tref ~ temp + aou 2998.8701 3015.1807 7.470474 2000-2010 –> 2011-2019 15.582636 9357.1140 9392.6518
Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate 2664.5940 2680.9046 5.091713 2000-2010 –> 2011-2019 10.722456 8362.6953 8398.2332
Atlantic (28.15,28.2] cstar_tref ~ temp + phosphate 2560.5015 2576.8121 4.518780 2000-2010 –> 2011-2019 9.554730 8056.7591 8092.2970
Atlantic (28.15,28.2] cstar_tref ~ temp + phosphate_star 2436.9772 2453.2878 3.921937 2000-2010 –> 2011-2019 7.766853 7445.3449 7480.8828
Atlantic (28.15,28.2] cstar_tref ~ temp + silicate 2862.6826 2878.9932 6.390293 2000-2010 –> 2011-2019 13.119116 8882.8953 8918.4332
Atlantic (28.2, Inf] cstar_tref ~ aou + phosphate 10995.8218 11018.3823 3.395207 1982-1999 –> 2000-2010 6.858716 16917.0546 16959.6671
Atlantic (28.2, Inf] cstar_tref ~ aou + phosphate_star 10974.7866 10997.3471 3.378083 1982-1999 –> 2000-2010 7.323913 17185.7176 17228.3302
Atlantic (28.2, Inf] cstar_tref ~ nitrate + phosphate_star 11172.5654 11195.1259 3.542566 1982-1999 –> 2000-2010 7.535567 17409.9018 17452.5144
Atlantic (28.2, Inf] cstar_tref ~ phosphate + phosphate_star 10941.7341 10964.2946 3.351349 1982-1999 –> 2000-2010 6.996125 16976.3168 17018.9294
Atlantic (28.2, Inf] cstar_tref ~ sal + phosphate 11765.8539 11788.4144 4.085600 1982-1999 –> 2000-2010 7.956753 17934.3286 17976.9412
Atlantic (28.2, Inf] cstar_tref ~ silicate + phosphate_star 11410.1169 11432.6774 3.750747 1982-1999 –> 2000-2010 8.006714 17789.1702 17831.7827
Atlantic (28.2, Inf] cstar_tref ~ temp + aou + phosphate 10985.1077 11013.3083 3.384846 1982-1999 –> 2000-2010 6.431378 16623.3054 16676.5710
Atlantic (28.2, Inf] cstar_tref ~ temp + aou + phosphate_star 10898.4492 10926.6498 3.315065 1982-1999 –> 2000-2010 6.382138 16551.5789 16604.8446
Atlantic (28.2, Inf] cstar_tref ~ temp + phosphate 11514.9482 11537.5087 3.846466 1982-1999 –> 2000-2010 7.717034 17683.0872 17725.6998
Atlantic (28.2, Inf] cstar_tref ~ temp + phosphate_star 11732.3526 11754.9131 4.052830 1982-1999 –> 2000-2010 7.170279 17419.6812 17462.2937
Atlantic (28.2, Inf] cstar_tref ~ aou + phosphate 5495.6936 5515.2110 4.071246 2000-2010 –> 2011-2019 7.466453 16491.5154 16533.5933
Atlantic (28.2, Inf] cstar_tref ~ aou + phosphate_star 5678.8099 5698.3274 4.473382 2000-2010 –> 2011-2019 7.851465 16653.5965 16695.6745
Atlantic (28.2, Inf] cstar_tref ~ nitrate + phosphate_star 5700.3723 5719.8897 4.523276 2000-2010 –> 2011-2019 8.065842 16872.9377 16915.0156
Atlantic (28.2, Inf] cstar_tref ~ phosphate + phosphate_star 5591.4722 5610.9896 4.276855 2000-2010 –> 2011-2019 7.628204 16533.2063 16575.2842
Atlantic (28.2, Inf] cstar_tref ~ silicate + phosphate_star 5897.3350 5916.8524 5.005587 2000-2010 –> 2011-2019 8.756334 17307.4519 17349.5298
Atlantic (28.2, Inf] cstar_tref ~ temp + aou + nitrate 5830.9019 5855.2987 4.832445 2000-2010 –> 2011-2019 8.989565 17670.9480 17723.5454
Atlantic (28.2, Inf] cstar_tref ~ temp + aou + phosphate 5403.8617 5428.2584 3.879405 2000-2010 –> 2011-2019 7.264251 16388.9693 16441.5667
Atlantic (28.2, Inf] cstar_tref ~ temp + aou + phosphate_star 5444.9455 5469.3423 3.962263 2000-2010 –> 2011-2019 7.277328 16343.3947 16395.9921
Atlantic (28.2, Inf] cstar_tref ~ temp + phosphate 5905.1372 5924.6547 5.025718 2000-2010 –> 2011-2019 8.872184 17420.0854 17462.1633
Atlantic (28.2, Inf] cstar_tref ~ temp + phosphate_star 5487.8481 5507.3655 4.054849 2000-2010 –> 2011-2019 8.107679 17220.2008 17262.2787
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + phosphate_star 23823.9339 23861.0074 6.826033 1982-1999 –> 2000-2010 15.536394 32995.2366 33063.2284
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 23448.2206 23491.4730 6.473833 1982-1999 –> 2000-2010 14.569690 32434.5241 32513.8479
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + nitrate 24486.3573 24517.2519 7.492711 1982-1999 –> 2000-2010 16.446606 33726.1429 33782.8027
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate_star 23792.8541 23829.9276 6.796343 1982-1999 –> 2000-2010 15.270177 32893.7897 32961.7815
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + nitrate + silicate 24454.1574 24491.2309 7.456857 1982-1999 –> 2000-2010 16.269475 33655.2919 33723.2837
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + phosphate 24479.0441 24509.9387 7.485029 1982-1999 –> 2000-2010 16.526431 33743.6880 33800.3479
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + phosphate_star 23831.1007 23861.9953 6.834814 1982-1999 –> 2000-2010 15.678728 33039.2961 33095.9559
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate 24393.1598 24430.2333 7.393335 1982-1999 –> 2000-2010 16.122599 33570.0033 33637.9951
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate_star 23675.6956 23712.7691 6.685579 1982-1999 –> 2000-2010 14.838179 32677.8514 32745.8432
Indo-Pacific (-Inf,26] cstar_tref ~ temp + aou + silicate + phosphate_star 24677.0863 24714.1598 7.693689 1982-1999 –> 2000-2010 16.485233 33872.1013 33940.0931
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + phosphate_star 24393.0090 24430.4782 5.943804 2000-2010 –> 2011-2019 12.769837 48216.9429 48291.4856
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate 25243.3907 25280.8598 6.645943 2000-2010 –> 2011-2019 14.362627 49941.7551 50016.2978
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 24048.5769 24092.2909 5.679493 2000-2010 –> 2011-2019 12.153327 47496.7975 47583.7639
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + nitrate 25358.7441 25389.9684 6.749143 2000-2010 –> 2011-2019 14.241853 49845.1014 49907.2203
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + nitrate + silicate 25292.6211 25330.0903 6.689043 2000-2010 –> 2011-2019 14.145900 49746.7785 49821.3212
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + phosphate 25550.7430 25581.9673 6.921451 2000-2010 –> 2011-2019 14.406480 50029.7871 50091.9060
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + phosphate_star 24844.5764 24875.8007 6.308538 2000-2010 –> 2011-2019 13.143352 48675.6771 48737.7960
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate 25351.2322 25388.7013 6.740719 2000-2010 –> 2011-2019 14.134054 49744.3919 49818.9346
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate_star 24822.9189 24860.3881 6.288972 2000-2010 –> 2011-2019 12.974552 48498.6145 48573.1572
Indo-Pacific (-Inf,26] cstar_tref ~ temp + aou + silicate + phosphate_star 24511.6585 24549.1276 6.037128 2000-2010 –> 2011-2019 13.730817 49188.7448 49263.2874
Indo-Pacific (26,26.5] cstar_tref ~ aou + silicate + phosphate_star 17993.9139 18023.4742 6.519636 1982-1999 –> 2000-2010 14.820826 28546.2899 28602.3862
Indo-Pacific (26,26.5] cstar_tref ~ phosphate + phosphate_star 17905.9239 17929.5721 6.417762 1982-1999 –> 2000-2010 14.663462 28436.2995 28481.1766
Indo-Pacific (26,26.5] cstar_tref ~ sal + aou + phosphate_star 17996.6381 18026.1984 6.522890 1982-1999 –> 2000-2010 14.819226 28547.2702 28603.3665
Indo-Pacific (26,26.5] cstar_tref ~ sal + aou + silicate + phosphate_star 17994.9618 18030.4341 6.518499 1982-1999 –> 2000-2010 14.664377 28493.0171 28560.3327
Indo-Pacific (26,26.5] cstar_tref ~ sal + nitrate + phosphate_star 18046.6880 18076.2483 6.582958 1982-1999 –> 2000-2010 14.871455 28594.5010 28650.5973
Indo-Pacific (26,26.5] cstar_tref ~ sal + nitrate + silicate + phosphate_star 18047.3129 18082.7852 6.581300 1982-1999 –> 2000-2010 14.760386 28557.5002 28624.8157
Indo-Pacific (26,26.5] cstar_tref ~ sal + silicate + phosphate 17977.0920 18006.6522 6.499581 1982-1999 –> 2000-2010 14.779263 28521.7319 28577.8282
Indo-Pacific (26,26.5] cstar_tref ~ silicate + phosphate + phosphate_star 17900.4002 17929.9604 6.408925 1982-1999 –> 2000-2010 14.564904 28400.1512 28456.2475
Indo-Pacific (26,26.5] cstar_tref ~ temp + phosphate_star 17400.8047 17424.4529 5.850674 1982-1999 –> 2000-2010 13.446549 27686.3987 27731.2758
Indo-Pacific (26,26.5] cstar_tref ~ temp + silicate + phosphate_star 17229.7035 17259.2637 5.668096 1982-1999 –> 2000-2010 13.037069 27426.8619 27482.9582
Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou 12284.7260 12312.5264 5.915670 1982-1999 –> 2011-2019 14.919227 23079.3029 23133.6393
Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate 12272.7705 12306.1310 5.894210 1982-1999 –> 2011-2019 14.613971 22973.8400 23039.0437
Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate 12300.5014 12328.3018 5.940022 1982-1999 –> 2011-2019 15.149840 23162.6214 23216.9578
Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + silicate 12295.1277 12328.4882 5.928627 1982-1999 –> 2011-2019 14.917895 23086.9682 23152.1719
Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + phosphate 11732.1081 11759.9085 5.122762 1982-1999 –> 2011-2019 13.229098 22213.6530 22267.9894
Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate 13406.6965 13434.4969 7.923116 1982-1999 –> 2011-2019 17.878382 24500.9104 24555.2468
Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate 11722.6589 11756.0194 5.107511 1982-1999 –> 2011-2019 13.104065 22165.5434 22230.7471
Indo-Pacific (26,26.5] cstar_tref ~ temp + nitrate + phosphate_star 11511.0078 11538.8082 4.836134 1982-1999 –> 2011-2019 12.309463 21750.1000 21804.4364
Indo-Pacific (26,26.5] cstar_tref ~ temp + nitrate + silicate + phosphate_star 11323.9239 11357.2843 4.603767 1982-1999 –> 2011-2019 11.966747 21520.6563 21585.8600
Indo-Pacific (26,26.5] cstar_tref ~ aou + phosphate_star 11913.7518 11935.9921 5.373705 2000-2010 –> 2011-2019 11.900375 29911.5527 29957.4413
Indo-Pacific (26,26.5] cstar_tref ~ aou + silicate + phosphate_star 11901.9776 11929.7780 5.354464 2000-2010 –> 2011-2019 11.874100 29895.8915 29953.2522
Indo-Pacific (26,26.5] cstar_tref ~ phosphate + phosphate_star 11842.1755 11864.4159 5.274469 2000-2010 –> 2011-2019 11.692231 29748.0994 29793.9879
Indo-Pacific (26,26.5] cstar_tref ~ sal + aou + phosphate_star 11899.4199 11927.2203 5.350899 2000-2010 –> 2011-2019 11.873789 29896.0580 29953.4187
Indo-Pacific (26,26.5] cstar_tref ~ sal + aou + silicate + phosphate_star 11897.2359 11930.5964 5.345072 2000-2010 –> 2011-2019 11.863571 29892.1977 29961.0305
Indo-Pacific (26,26.5] cstar_tref ~ sal + phosphate 11825.4275 11847.6678 5.251515 2000-2010 –> 2011-2019 11.751177 29800.5880 29846.4766
Indo-Pacific (26,26.5] cstar_tref ~ sal + silicate + phosphate 11825.3597 11853.1601 5.248688 2000-2010 –> 2011-2019 11.748268 29802.4517 29859.8124
Indo-Pacific (26,26.5] cstar_tref ~ silicate + phosphate + phosphate_star 11821.5933 11849.3937 5.243542 2000-2010 –> 2011-2019 11.652467 29721.9935 29779.3542
Indo-Pacific (26,26.5] cstar_tref ~ temp + phosphate_star 11650.6653 11672.9056 5.017870 2000-2010 –> 2011-2019 10.868544 29051.4700 29097.3585
Indo-Pacific (26,26.5] cstar_tref ~ temp + silicate + phosphate_star 11325.4676 11353.2680 4.608018 2000-2010 –> 2011-2019 10.276114 28555.1711 28612.5317
Indo-Pacific (26.5,26.75] cstar_tref ~ aou + phosphate_star 19033.2147 19057.0963 6.475797 1982-1999 –> 2000-2010 12.986472 29675.6078 29721.0378
Indo-Pacific (26.5,26.75] cstar_tref ~ aou + silicate + phosphate_star 18729.4846 18759.3365 6.142613 1982-1999 –> 2000-2010 12.643290 29368.9134 29425.7008
Indo-Pacific (26.5,26.75] cstar_tref ~ nitrate + phosphate_star 19216.4593 19240.3409 6.684096 1982-1999 –> 2000-2010 13.377598 29948.3040 29993.7339
Indo-Pacific (26.5,26.75] cstar_tref ~ nitrate + silicate + phosphate_star 18840.9132 18870.7651 6.262014 1982-1999 –> 2000-2010 12.925350 29560.1680 29616.9554
Indo-Pacific (26.5,26.75] cstar_tref ~ phosphate + phosphate_star 18911.0276 18934.9092 6.340523 1982-1999 –> 2000-2010 12.717053 29486.1749 29531.6049
Indo-Pacific (26.5,26.75] cstar_tref ~ sal + silicate + phosphate 19255.9236 19285.7756 6.727501 1982-1999 –> 2000-2010 13.610522 30079.9512 30136.7386
Indo-Pacific (26.5,26.75] cstar_tref ~ silicate + phosphate + phosphate_star 18645.1649 18675.0169 6.053776 1982-1999 –> 2000-2010 12.419940 29217.0571 29273.8445
Indo-Pacific (26.5,26.75] cstar_tref ~ temp + phosphate_star 18896.9069 18920.7885 6.325073 1982-1999 –> 2000-2010 12.232262 29225.1079 29270.5378
Indo-Pacific (26.5,26.75] cstar_tref ~ temp + silicate + phosphate 19348.6601 19378.5121 6.836159 1982-1999 –> 2000-2010 13.753862 30188.9222 30245.7096
Indo-Pacific (26.5,26.75] cstar_tref ~ temp + silicate + phosphate_star 18803.3527 18833.2047 6.221509 1982-1999 –> 2000-2010 12.040317 29084.8622 29141.6496
Indo-Pacific (26.5,26.75] cstar_tref ~ sal + aou + phosphate_star 10829.1781 10856.3934 5.744645 1982-1999 –> 2011-2019 11.953275 21320.1364 21374.2873
Indo-Pacific (26.5,26.75] cstar_tref ~ sal + aou + silicate + phosphate_star 10826.9743 10859.6328 5.737580 1982-1999 –> 2011-2019 11.944423 21319.0032 21383.9842
Indo-Pacific (26.5,26.75] cstar_tref ~ sal + nitrate + phosphate_star 10856.9572 10884.1726 5.791552 1982-1999 –> 2011-2019 11.831431 21258.9090 21313.0598
Indo-Pacific (26.5,26.75] cstar_tref ~ sal + nitrate + silicate + phosphate_star 10858.4715 10891.1300 5.790728 1982-1999 –> 2011-2019 11.821497 21257.5476 21322.5286
Indo-Pacific (26.5,26.75] cstar_tref ~ aou + phosphate_star 10905.7579 10927.5303 5.878324 2000-2010 –> 2011-2019 12.354120 29938.9726 29984.6265
Indo-Pacific (26.5,26.75] cstar_tref ~ aou + silicate + phosphate_star 10857.4890 10884.7044 5.792453 2000-2010 –> 2011-2019 11.935066 29586.9735 29644.0409
Indo-Pacific (26.5,26.75] cstar_tref ~ nitrate + phosphate_star 11063.7951 11085.5674 6.156666 2000-2010 –> 2011-2019 12.840763 30280.2544 30325.9083
Indo-Pacific (26.5,26.75] cstar_tref ~ nitrate + silicate + phosphate_star 10993.1944 11020.4097 6.027198 2000-2010 –> 2011-2019 12.289213 29834.1075 29891.1749
Indo-Pacific (26.5,26.75] cstar_tref ~ phosphate + phosphate_star 10816.0010 10837.7733 5.725880 2000-2010 –> 2011-2019 12.066402 29727.0286 29772.6825
Indo-Pacific (26.5,26.75] cstar_tref ~ sal + silicate + phosphate 10960.0077 10987.2231 5.968927 2000-2010 –> 2011-2019 12.696428 30215.9313 30272.9987
Indo-Pacific (26.5,26.75] cstar_tref ~ silicate + phosphate + phosphate_star 10782.8783 10810.0937 5.667309 2000-2010 –> 2011-2019 11.721085 29428.0433 29485.1106
Indo-Pacific (26.5,26.75] cstar_tref ~ temp + phosphate_star 10748.4570 10770.2293 5.613775 2000-2010 –> 2011-2019 11.938848 29645.3639 29691.0178
Indo-Pacific (26.5,26.75] cstar_tref ~ temp + silicate + phosphate 11015.9649 11043.1803 6.067509 2000-2010 –> 2011-2019 12.903668 30364.6250 30421.6923
Indo-Pacific (26.5,26.75] cstar_tref ~ temp + silicate + phosphate_star 10444.6884 10471.9038 5.133115 2000-2010 –> 2011-2019 11.354624 29248.0411 29305.1085
Indo-Pacific (26.75,27] cstar_tref ~ aou + phosphate_star 29406.5940 29432.1304 6.954337 1982-1999 –> 2000-2010 13.544151 47886.5458 47935.8246
Indo-Pacific (26.75,27] cstar_tref ~ aou + silicate + phosphate_star 27535.5109 27567.4315 5.614748 1982-1999 –> 2000-2010 11.632279 45509.6223 45571.2208
Indo-Pacific (26.75,27] cstar_tref ~ nitrate + phosphate_star 29264.4969 29290.0334 6.842364 1982-1999 –> 2000-2010 13.436645 47748.2370 47797.5159
Indo-Pacific (26.75,27] cstar_tref ~ nitrate + silicate + phosphate_star 27652.7625 27684.6831 5.690458 1982-1999 –> 2000-2010 11.852324 45759.3709 45820.9694
Indo-Pacific (26.75,27] cstar_tref ~ phosphate + phosphate_star 29170.0306 29195.5671 6.768924 1982-1999 –> 2000-2010 13.237498 47546.1797 47595.4586
Indo-Pacific (26.75,27] cstar_tref ~ sal + aou + phosphate_star 28173.4868 28205.4074 6.039220 1982-1999 –> 2000-2010 12.054481 46145.4878 46207.0863
Indo-Pacific (26.75,27] cstar_tref ~ sal + silicate + phosphate 28533.5179 28565.4385 6.292778 1982-1999 –> 2000-2010 12.893438 47024.6633 47086.2618
Indo-Pacific (26.75,27] cstar_tref ~ sal + silicate + phosphate_star 29479.6997 29511.6203 7.011055 1982-1999 –> 2000-2010 13.395253 47784.4527 47846.0512
Indo-Pacific (26.75,27] cstar_tref ~ silicate + phosphate + phosphate_star 27344.0828 27376.0034 5.493300 1982-1999 –> 2000-2010 11.395359 45209.8830 45271.4815
Indo-Pacific (26.75,27] cstar_tref ~ temp + silicate + phosphate 28265.9851 28297.9057 6.103372 1982-1999 –> 2000-2010 12.452342 46539.8075 46601.4060
Indo-Pacific (26.75,27] cstar_tref ~ sal + nitrate + phosphate_star 17563.5744 17593.3198 5.360812 1982-1999 –> 2011-2019 11.237398 35405.1957 35464.6191
Indo-Pacific (26.75,27] cstar_tref ~ sal + phosphate + phosphate_star 17505.4951 17535.2405 5.306142 1982-1999 –> 2011-2019 11.297564 35455.2992 35514.7226
Indo-Pacific (26.75,27] cstar_tref ~ temp + phosphate_star 18366.0280 18389.8244 6.178618 1982-1999 –> 2011-2019 12.233052 36372.3155 36419.8542
Indo-Pacific (26.75,27] cstar_tref ~ aou + silicate + phosphate_star 17453.9159 17483.6614 5.258058 2000-2010 –> 2011-2019 10.872806 44989.4268 45051.0928
Indo-Pacific (26.75,27] cstar_tref ~ nitrate + phosphate_star 18515.2282 18539.0245 6.343478 2000-2010 –> 2011-2019 13.185842 47779.7251 47829.0579
Indo-Pacific (26.75,27] cstar_tref ~ nitrate + silicate + phosphate_star 17603.3893 17633.1348 5.398616 2000-2010 –> 2011-2019 11.089074 45256.1518 45317.8179
Indo-Pacific (26.75,27] cstar_tref ~ phosphate + phosphate_star 18312.3048 18336.1012 6.120311 2000-2010 –> 2011-2019 12.889235 47482.3354 47531.6682
Indo-Pacific (26.75,27] cstar_tref ~ sal + aou + phosphate_star 17600.0187 17629.7641 5.395405 2000-2010 –> 2011-2019 11.434625 45773.5054 45835.1715
Indo-Pacific (26.75,27] cstar_tref ~ sal + nitrate + silicate 18219.3320 18249.0775 6.018579 2000-2010 –> 2011-2019 12.837460 47455.7339 47517.4000
Indo-Pacific (26.75,27] cstar_tref ~ sal + silicate + phosphate 17716.2943 17746.0398 5.507271 2000-2010 –> 2011-2019 11.800049 46249.8122 46311.4783
Indo-Pacific (26.75,27] cstar_tref ~ silicate + phosphate + phosphate_star 17313.1994 17342.9449 5.129081 2000-2010 –> 2011-2019 10.622381 44657.2823 44718.9483
Indo-Pacific (26.75,27] cstar_tref ~ temp + nitrate + silicate 18156.4051 18186.1505 5.952106 2000-2010 –> 2011-2019 12.691014 47289.5320 47351.1980
Indo-Pacific (26.75,27] cstar_tref ~ temp + silicate + phosphate 17588.3978 17618.1432 5.384350 2000-2010 –> 2011-2019 11.487722 45854.3829 45916.0489
Indo-Pacific (27,27.25] cstar_tref ~ aou + silicate + phosphate_star 37091.1297 37124.5407 5.613244 1982-1999 –> 2000-2010 12.336203 59973.4618 60037.5889
Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + phosphate_star 37110.8717 37144.2827 5.622648 1982-1999 –> 2000-2010 12.304375 59950.8793 60015.0064
Indo-Pacific (27,27.25] cstar_tref ~ sal + phosphate + phosphate_star 37123.1505 37156.5615 5.628505 1982-1999 –> 2000-2010 12.328125 59981.5577 60045.6849
Indo-Pacific (27,27.25] cstar_tref ~ sal + silicate + phosphate_star 36871.7783 36905.1893 5.509810 1982-1999 –> 2000-2010 11.843215 59343.4408 59407.5679
Indo-Pacific (27,27.25] cstar_tref ~ silicate + phosphate + phosphate_star 36771.7972 36805.2082 5.463299 1982-1999 –> 2000-2010 11.958971 59417.5100 59481.6371
Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + phosphate_star 36933.6260 36967.0370 5.538780 1982-1999 –> 2000-2010 11.991717 59533.9264 59598.0535
Indo-Pacific (27,27.25] cstar_tref ~ temp + nitrate + phosphate_star 36901.0465 36934.4575 5.523501 1982-1999 –> 2000-2010 12.040414 59569.2214 59633.3485
Indo-Pacific (27,27.25] cstar_tref ~ temp + phosphate + phosphate_star 36932.0472 36965.4582 5.538038 1982-1999 –> 2000-2010 11.995498 59537.1678 59601.2950
Indo-Pacific (27,27.25] cstar_tref ~ temp + phosphate_star 36950.7165 36977.4453 5.547752 1982-1999 –> 2000-2010 12.082470 59635.6620 59686.9637
Indo-Pacific (27,27.25] cstar_tref ~ temp + silicate + phosphate_star 36949.0270 36982.4380 5.546017 1982-1999 –> 2000-2010 11.833879 59371.0372 59435.1644
Indo-Pacific (27,27.25] cstar_tref ~ sal + nitrate + silicate + phosphate_star 22611.2884 22648.7984 4.610250 1982-1999 –> 2011-2019 10.347605 44404.9336 44479.3029
Indo-Pacific (27,27.25] cstar_tref ~ aou + silicate + phosphate_star 23423.7466 23455.0049 5.126881 2000-2010 –> 2011-2019 10.740124 60514.8763 60579.5456
Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + phosphate_star 23414.5816 23445.8399 5.120756 2000-2010 –> 2011-2019 10.743404 60525.4533 60590.1226
Indo-Pacific (27,27.25] cstar_tref ~ sal + phosphate + phosphate_star 23408.5034 23439.7618 5.116699 2000-2010 –> 2011-2019 10.745203 60531.6539 60596.3232
Indo-Pacific (27,27.25] cstar_tref ~ sal + silicate + phosphate_star 23331.7628 23363.0211 5.065747 2000-2010 –> 2011-2019 10.575557 60203.5411 60268.2104
Indo-Pacific (27,27.25] cstar_tref ~ silicate + phosphate + phosphate_star 23202.5398 23233.7981 4.981093 2000-2010 –> 2011-2019 10.444392 59974.3369 60039.0062
Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + phosphate_star 22746.6346 22777.8929 4.693571 2000-2010 –> 2011-2019 10.232351 59680.2606 59744.9299
Indo-Pacific (27,27.25] cstar_tref ~ temp + nitrate + phosphate_star 22652.2127 22683.4710 4.636130 2000-2010 –> 2011-2019 10.159631 59553.2593 59617.9286
Indo-Pacific (27,27.25] cstar_tref ~ temp + phosphate + phosphate_star 22741.8807 22773.1391 4.690663 2000-2010 –> 2011-2019 10.228701 59673.9280 59738.5973
Indo-Pacific (27,27.25] cstar_tref ~ temp + phosphate_star 23102.9400 23127.9467 4.918094 2000-2010 –> 2011-2019 10.465847 60053.6565 60105.3919
Indo-Pacific (27,27.25] cstar_tref ~ temp + silicate + phosphate_star 22879.3457 22910.6040 4.775511 2000-2010 –> 2011-2019 10.321528 59828.3727 59893.0420
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + nitrate + silicate + phosphate_star 36016.7206 36056.9734 4.728983 1982-1999 –> 2000-2010 10.641525 59621.6992 59699.2340
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + silicate + phosphate 37960.5708 37994.1148 5.553147 1982-1999 –> 2000-2010 12.356639 62599.6869 62664.2992
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + silicate + phosphate_star 38002.8465 38036.3905 5.572564 1982-1999 –> 2000-2010 12.015794 62240.3385 62304.9508
Indo-Pacific (27.25,27.5] cstar_tref ~ silicate + phosphate + phosphate_star 37533.1174 37566.6614 5.360585 1982-1999 –> 2000-2010 11.925115 61908.2898 61972.9021
Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + phosphate_star 36537.8097 36571.3538 4.937691 1982-1999 –> 2000-2010 10.638228 59871.2310 59935.8432
Indo-Pacific (27.25,27.5] cstar_tref ~ temp + nitrate + phosphate_star 36722.5460 36756.0900 5.013579 1982-1999 –> 2000-2010 10.898430 60290.8692 60355.4815
Indo-Pacific (27.25,27.5] cstar_tref ~ temp + phosphate + phosphate_star 36546.7973 36580.3413 4.941356 1982-1999 –> 2000-2010 10.659125 59902.4987 59967.1110
Indo-Pacific (27.25,27.5] cstar_tref ~ temp + phosphate_star 37311.3025 37338.1377 5.264176 1982-1999 –> 2000-2010 11.605904 61429.5776 61481.2674
Indo-Pacific (27.25,27.5] cstar_tref ~ temp + silicate + phosphate 37952.6501 37986.1941 5.549517 1982-1999 –> 2000-2010 12.310038 62544.9946 62609.6069
Indo-Pacific (27.25,27.5] cstar_tref ~ temp + silicate + phosphate_star 36182.5722 36216.1163 4.794975 1982-1999 –> 2000-2010 10.291472 59246.9248 59311.5371
Indo-Pacific (27.25,27.5] cstar_tref ~ aou + silicate + phosphate_star 27483.2661 27515.3019 5.196362 1982-1999 –> 2011-2019 12.070755 52198.9141 52262.0182
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + silicate 28016.6120 28048.6477 5.515141 1982-1999 –> 2011-2019 13.402950 53747.3961 53810.5002
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + silicate + phosphate_star 24189.5151 24227.9580 3.596810 1982-1999 –> 2011-2019 9.471877 47747.5565 47823.2813
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + silicate + phosphate + phosphate_star 24041.6085 24080.0515 3.537910 1982-1999 –> 2011-2019 9.317758 47479.0267 47554.7516
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + nitrate + silicate + phosphate_star 24872.1811 24910.6240 3.881628 2000-2010 –> 2011-2019 8.610611 60888.9017 60967.5974
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + phosphate_star 24392.4557 24430.8986 3.679224 2000-2010 –> 2011-2019 8.464123 60551.5499 60630.2457
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + phosphate_star 24487.8729 24526.3158 3.718624 2000-2010 –> 2011-2019 8.589060 60861.5745 60940.2702
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + phosphate + phosphate_star 24393.8612 24432.3041 3.679802 2000-2010 –> 2011-2019 8.468381 60562.2681 60640.9639
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate 24975.3686 25013.8116 3.926599 2000-2010 –> 2011-2019 8.893558 61586.7589 61665.4547
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate_star 24005.1109 24043.5538 3.523525 2000-2010 –> 2011-2019 8.098257 59620.1724 59698.8681
Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + phosphate_star 25000.7731 25032.8089 3.938630 2000-2010 –> 2011-2019 8.876321 61538.5829 61604.1627
Indo-Pacific (27.25,27.5] cstar_tref ~ temp + nitrate + phosphate_star 25091.6700 25123.7058 3.978799 2000-2010 –> 2011-2019 8.992378 61814.2160 61879.7958
Indo-Pacific (27.25,27.5] cstar_tref ~ temp + phosphate + phosphate_star 25004.2863 25036.3220 3.940175 2000-2010 –> 2011-2019 8.881531 61551.0836 61616.6634
Indo-Pacific (27.25,27.5] cstar_tref ~ temp + silicate + phosphate_star 24867.6420 24899.6778 3.880528 2000-2010 –> 2011-2019 8.675503 61050.2142 61115.7940
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + nitrate + silicate + phosphate_star 39801.1685 39842.1086 4.529506 1982-1999 –> 2000-2010 10.653189 65448.7332 65527.3879
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + silicate + phosphate + phosphate_star 39003.4124 39044.3525 4.271123 1982-1999 –> 2000-2010 10.162595 64344.2665 64422.9212
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + phosphate_star 37711.1375 37752.0776 3.883476 1982-1999 –> 2000-2010 8.786293 61594.5492 61673.2039
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + phosphate_star 38120.2429 38161.1830 4.002230 1982-1999 –> 2000-2010 9.040933 62220.5596 62299.2143
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + phosphate + phosphate_star 37726.6581 37767.5982 3.887916 1982-1999 –> 2000-2010 8.809859 61640.9587 61719.6134
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate_star 37139.2333 37180.1734 3.723347 1982-1999 –> 2000-2010 8.692946 61129.9850 61208.6397
Indo-Pacific (27.5,27.75] cstar_tref ~ temp + aou + phosphate_star 37964.9243 37999.0411 3.957305 1982-1999 –> 2000-2010 8.874791 61870.0370 61935.5826
Indo-Pacific (27.5,27.75] cstar_tref ~ temp + nitrate + phosphate_star 38330.8361 38364.9528 4.065368 1982-1999 –> 2000-2010 9.150777 62502.3579 62567.9035
Indo-Pacific (27.5,27.75] cstar_tref ~ temp + phosphate + phosphate_star 37983.6859 38017.8027 3.962776 1982-1999 –> 2000-2010 8.899797 61920.2550 61985.8006
Indo-Pacific (27.5,27.75] cstar_tref ~ temp + silicate + phosphate_star 37824.9817 37859.0984 3.916740 1982-1999 –> 2000-2010 8.892145 61822.9967 61888.5423
Indo-Pacific (27.5,27.75] cstar_tref ~ aou + silicate 32209.7516 32235.8697 5.822883 1982-1999 –> 2011-2019 13.297904 59435.3819 59486.6430
Indo-Pacific (27.5,27.75] cstar_tref ~ aou + silicate + phosphate_star 31381.6461 31414.2936 5.364492 1982-1999 –> 2011-2019 12.567128 58314.7675 58378.8439
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + silicate 30232.1889 30264.8365 4.788723 1982-1999 –> 2011-2019 12.119914 57305.6695 57369.7459
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate_star 27404.6417 27443.8188 3.621089 1982-1999 –> 2011-2019 9.659325 52940.7197 53017.6114
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + nitrate + silicate + phosphate_star 28077.3880 28116.5651 3.869887 2000-2010 –> 2011-2019 8.399394 67878.5565 67958.6738
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + silicate + phosphate + phosphate_star 27126.1068 27165.2839 3.522822 2000-2010 –> 2011-2019 7.793945 66129.5192 66209.6364
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + phosphate_star 27497.1278 27536.3049 3.654320 2000-2010 –> 2011-2019 7.537796 65208.2653 65288.3826
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + phosphate_star 27564.6353 27603.8124 3.678769 2000-2010 –> 2011-2019 7.680999 65684.8781 65764.9954
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + phosphate + phosphate_star 27489.2134 27528.3905 3.651465 2000-2010 –> 2011-2019 7.539381 65215.8715 65295.9888
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate_star 26668.9420 26708.1191 3.367282 2000-2010 –> 2011-2019 7.090629 63808.1753 63888.2925
Indo-Pacific (27.5,27.75] cstar_tref ~ temp + aou + phosphate_star 27803.7865 27836.4340 3.767448 2000-2010 –> 2011-2019 7.724754 65768.7108 65835.4751
Indo-Pacific (27.5,27.75] cstar_tref ~ temp + nitrate + phosphate_star 27861.2521 27893.8997 3.788894 2000-2010 –> 2011-2019 7.854262 66192.0882 66258.8525
Indo-Pacific (27.5,27.75] cstar_tref ~ temp + phosphate + phosphate_star 27800.7758 27833.4234 3.766328 2000-2010 –> 2011-2019 7.729104 65784.4618 65851.2261
Indo-Pacific (27.5,27.75] cstar_tref ~ temp + silicate + phosphate_star 27618.8987 27651.5463 3.699270 2000-2010 –> 2011-2019 7.616011 65443.8803 65510.6447
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + aou + phosphate_star 19977.9431 20008.6301 4.483264 1982-1999 –> 2000-2010 9.835116 31481.7101 31540.0279
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + nitrate + phosphate_star 20126.3121 20156.9991 4.581575 1982-1999 –> 2000-2010 10.063253 31719.0509 31777.3688
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + silicate + phosphate_star 19282.2848 19312.9717 4.049717 1982-1999 –> 2000-2010 9.070775 30549.2188 30607.5367
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + phosphate_star 17768.0247 17804.8491 3.244533 1982-1999 –> 2000-2010 7.316736 28259.4553 28329.4368
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + nitrate + phosphate_star 17945.5955 17982.4199 3.329866 1982-1999 –> 2000-2010 7.455995 28485.8596 28555.8410
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate_star 17743.9036 17780.7279 3.233112 1982-1999 –> 2000-2010 7.492612 28402.2555 28472.2370
Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + phosphate_star 17775.3400 17806.0270 3.248955 1982-1999 –> 2000-2010 7.341634 28283.3896 28341.7074
Indo-Pacific (27.75,27.85] cstar_tref ~ temp + nitrate + phosphate_star 17968.2811 17998.9681 3.341906 1982-1999 –> 2000-2010 7.553756 28582.8739 28641.1917
Indo-Pacific (27.75,27.85] cstar_tref ~ temp + phosphate + phosphate_star 17777.5893 17808.2763 3.250024 1982-1999 –> 2000-2010 7.349047 28291.3885 28349.7064
Indo-Pacific (27.75,27.85] cstar_tref ~ temp + silicate + phosphate_star 17771.0716 17801.7586 3.246928 1982-1999 –> 2000-2010 7.508771 28429.4643 28487.7821
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + aou + phosphate_star 15577.0595 15606.5162 4.445447 2000-2010 –> 2011-2019 8.928711 35555.0026 35615.1463
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + silicate + phosphate_star 14841.1207 14870.5774 3.873933 2000-2010 –> 2011-2019 7.923650 34123.4055 34183.5491
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + phosphate_star 14087.0469 14122.3949 3.363208 2000-2010 –> 2011-2019 6.607741 31855.0716 31927.2440
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + nitrate + phosphate_star 14071.9793 14107.3272 3.353746 2000-2010 –> 2011-2019 6.683612 32017.5747 32089.7471
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + phosphate + phosphate_star 14082.3474 14117.6954 3.360254 2000-2010 –> 2011-2019 6.605924 31852.7696 31924.9419
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate_star 13888.2407 13923.5886 3.240479 2000-2010 –> 2011-2019 6.473591 31632.1442 31704.3166
Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + phosphate_star 14085.8351 14115.2917 3.363704 2000-2010 –> 2011-2019 6.612658 31861.1751 31921.3187
Indo-Pacific (27.75,27.85] cstar_tref ~ temp + nitrate + phosphate_star 14069.9983 14099.4550 3.353757 2000-2010 –> 2011-2019 6.695663 32038.2794 32098.4231
Indo-Pacific (27.75,27.85] cstar_tref ~ temp + phosphate + phosphate_star 14081.1673 14110.6239 3.360769 2000-2010 –> 2011-2019 6.610792 31858.7566 31918.9002
Indo-Pacific (27.75,27.85] cstar_tref ~ temp + silicate + phosphate_star 13979.5669 14009.0236 3.297524 2000-2010 –> 2011-2019 6.544452 31750.6385 31810.7822
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + aou + phosphate_star 27376.6614 27408.9656 4.377123 1982-1999 –> 2000-2010 9.645697 43418.6147 43480.2391
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + silicate + phosphate_star 26473.0787 26505.3828 3.978062 1982-1999 –> 2000-2010 8.998494 42263.9724 42325.5968
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + phosphate_star 23498.4759 23537.2409 2.903384 1982-1999 –> 2000-2010 7.021256 38260.0437 38333.9930
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + nitrate + phosphate_star 23746.8109 23785.5759 2.980676 1982-1999 –> 2000-2010 7.069426 38471.4437 38545.3929
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + phosphate + phosphate_star 23511.1887 23549.9537 2.907292 1982-1999 –> 2000-2010 7.029901 38278.7393 38352.6885
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate_star 23509.5018 23548.2669 2.906773 1982-1999 –> 2000-2010 7.168983 38450.3537 38524.3029
Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + phosphate_star 23543.6499 23575.9540 2.917911 1982-1999 –> 2000-2010 7.035820 38303.2644 38364.8888
Indo-Pacific (27.85,27.95] cstar_tref ~ temp + nitrate + phosphate_star 23853.3155 23885.6197 3.015090 1982-1999 –> 2000-2010 7.125140 38602.9883 38664.6126
Indo-Pacific (27.85,27.95] cstar_tref ~ temp + phosphate + phosphate_star 23556.6653 23588.9695 2.921931 1982-1999 –> 2000-2010 7.044569 38322.2516 38383.8759
Indo-Pacific (27.85,27.95] cstar_tref ~ temp + silicate + phosphate_star 23551.7164 23584.0206 2.920402 1982-1999 –> 2000-2010 7.225538 38542.7176 38604.3419
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + aou + phosphate_star 22374.3585 22405.7469 4.148655 2000-2010 –> 2011-2019 8.525779 49751.0199 49814.7124
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + silicate + phosphate_star 21561.6590 21593.0474 3.741620 2000-2010 –> 2011-2019 7.719682 48034.7377 48098.4302
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + phosphate_star 20402.9147 20440.5807 3.228537 2000-2010 –> 2011-2019 6.131921 43901.3905 43977.8216
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + nitrate + phosphate_star 20365.1789 20402.8449 3.213093 2000-2010 –> 2011-2019 6.193770 44111.9898 44188.4208
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + phosphate + phosphate_star 20413.9696 20451.6356 3.233075 2000-2010 –> 2011-2019 6.140367 43925.1583 44001.5893
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate_star 20076.3575 20114.0235 3.097313 2000-2010 –> 2011-2019 6.004086 43585.8594 43662.2904
Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + phosphate_star 20417.6225 20449.0108 3.235398 2000-2010 –> 2011-2019 6.153309 43961.2723 44024.9648
Indo-Pacific (27.85,27.95] cstar_tref ~ temp + nitrate + phosphate_star 20389.5218 20420.9102 3.223866 2000-2010 –> 2011-2019 6.238957 44242.8373 44306.5298
Indo-Pacific (27.85,27.95] cstar_tref ~ temp + phosphate + phosphate_star 20428.3879 20459.7762 3.239827 2000-2010 –> 2011-2019 6.161758 43985.0532 44048.7457
Indo-Pacific (27.85,27.95] cstar_tref ~ temp + silicate + phosphate_star 20151.9368 20183.3251 3.127997 2000-2010 –> 2011-2019 6.048399 43703.6532 43767.3457
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + aou + phosphate_star 45910.8221 45945.8699 4.000392 1982-1999 –> 2000-2010 8.763298 76410.0999 76477.8484
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + silicate + phosphate_star 44569.6902 44604.7380 3.685573 1982-1999 –> 2000-2010 8.218118 74561.7230 74629.4715
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate_star 40842.5901 40884.6475 2.934432 1982-1999 –> 2000-2010 6.976602 69665.0396 69746.3378
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + nitrate + phosphate_star 41174.0867 41216.1441 2.994490 1982-1999 –> 2000-2010 7.026206 69970.0387 70051.3368
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + phosphate + phosphate_star 40918.1979 40960.2553 2.948023 1982-1999 –> 2000-2010 7.011971 69795.6278 69876.9260
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate_star 40495.7845 40537.8419 2.872888 1982-1999 –> 2000-2010 6.912914 69312.8049 69394.1031
Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + phosphate_star 40861.5966 40896.6445 2.938201 1982-1999 –> 2000-2010 7.029940 69806.7568 69874.5053
Indo-Pacific (27.95,28.05] cstar_tref ~ temp + nitrate + phosphate_star 41256.0618 41291.1096 3.009898 1982-1999 –> 2000-2010 7.051259 70074.4614 70142.2099
Indo-Pacific (27.95,28.05] cstar_tref ~ temp + phosphate + phosphate_star 40939.0939 40974.1417 2.952151 1982-1999 –> 2000-2010 7.066724 69941.1955 70008.9439
Indo-Pacific (27.95,28.05] cstar_tref ~ temp + silicate + phosphate_star 40814.2670 40849.3148 2.929714 1982-1999 –> 2000-2010 7.145881 70065.9407 70133.6892
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + aou + phosphate_star 36727.2884 36761.3370 3.749291 2000-2010 –> 2011-2019 7.749682 82638.1105 82707.2069
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + silicate + phosphate_star 35806.3246 35840.3732 3.500227 2000-2010 –> 2011-2019 7.185800 80376.0148 80445.1112
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate_star 34762.9970 34803.8553 3.237517 2000-2010 –> 2011-2019 6.171948 75605.5871 75688.5028
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + nitrate + phosphate_star 34525.6871 34566.5454 3.180678 2000-2010 –> 2011-2019 6.175167 75699.7738 75782.6895
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + phosphate + phosphate_star 34816.2340 34857.0923 3.250407 2000-2010 –> 2011-2019 6.198429 75734.4319 75817.3476
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate_star 34282.1688 34323.0270 3.123389 2000-2010 –> 2011-2019 5.996277 74777.9533 74860.8690
Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + phosphate_star 34769.5883 34803.6369 3.239593 2000-2010 –> 2011-2019 6.177795 75631.1850 75700.2814
Indo-Pacific (27.95,28.05] cstar_tref ~ temp + nitrate + phosphate_star 34552.3392 34586.3878 3.187487 2000-2010 –> 2011-2019 6.197385 75808.4010 75877.4974
Indo-Pacific (27.95,28.05] cstar_tref ~ temp + phosphate + phosphate_star 34824.3061 34858.3546 3.252851 2000-2010 –> 2011-2019 6.205002 75763.3999 75832.4963
Indo-Pacific (27.95,28.05] cstar_tref ~ temp + silicate + phosphate_star 34552.7583 34586.8069 3.187587 2000-2010 –> 2011-2019 6.117301 75367.0253 75436.1217
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + aou + phosphate_star 31133.7409 31167.3302 3.088222 1982-1999 –> 2000-2010 7.399004 51269.5837 51333.9669
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + nitrate + phosphate_star 31075.6932 31109.2824 3.073590 1982-1999 –> 2000-2010 7.458417 51330.5487 51394.9319
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + phosphate + phosphate_star 31200.2630 31233.8522 3.105077 1982-1999 –> 2000-2010 7.456203 51401.2033 51465.5865
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + silicate + phosphate_star 30731.9897 30765.5790 2.988359 1982-1999 –> 2000-2010 7.314570 50892.8004 50957.1836
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + phosphate_star 31678.4989 31712.0881 3.228984 1982-1999 –> 2000-2010 8.023608 52557.6977 52622.0809
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate_star 29588.1568 29628.4639 2.720928 1982-1999 –> 2000-2010 6.819402 49373.0760 49450.3359
Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + phosphate_star 30246.0613 30279.6506 2.871878 1982-1999 –> 2000-2010 7.023488 50118.9954 50183.3787
Indo-Pacific (28.05,28.1] cstar_tref ~ temp + nitrate + phosphate_star 30145.0318 30178.6210 2.848236 1982-1999 –> 2000-2010 7.043853 50091.6492 50156.0324
Indo-Pacific (28.05,28.1] cstar_tref ~ temp + phosphate + phosphate_star 30375.7688 30409.3581 2.902518 1982-1999 –> 2000-2010 7.097381 50321.1298 50385.5130
Indo-Pacific (28.05,28.1] cstar_tref ~ temp + silicate + phosphate_star 29774.1895 29807.7787 2.763112 1982-1999 –> 2000-2010 7.000011 49789.2275 49853.6108
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate_star 35690.7278 35731.8695 3.068522 1982-1999 –> 2011-2019 7.098319 55357.5587 55435.6532
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + nitrate 37293.2203 37327.5051 3.439847 1982-1999 –> 2011-2019 9.795917 60142.4410 60207.5198
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + nitrate + phosphate_star 35511.1068 35552.2484 3.029531 1982-1999 –> 2011-2019 7.173079 55372.4571 55450.5516
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + phosphate 37473.5684 37507.8531 3.484299 1982-1999 –> 2011-2019 9.780550 60256.7115 60321.7902
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + phosphate + phosphate_star 35741.2032 35782.3448 3.079568 1982-1999 –> 2011-2019 7.132233 55447.5746 55525.6691
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + aou + phosphate_star 35691.2116 35725.4963 3.069064 2000-2010 –> 2011-2019 6.157287 66824.9526 66892.8265
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + nitrate + phosphate_star 35535.1080 35569.3928 3.035144 2000-2010 –> 2011-2019 6.108734 66610.8012 66678.6752
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + phosphate + phosphate_star 35739.8869 35774.1716 3.079718 2000-2010 –> 2011-2019 6.184795 66940.1498 67008.0238
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + silicate + phosphate_star 35369.8356 35404.1204 2.999640 2000-2010 –> 2011-2019 5.988000 66101.8254 66169.6993
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + phosphate_star 35978.2472 36012.5320 3.132427 2000-2010 –> 2011-2019 6.361411 67656.7462 67724.6201
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate_star 35325.9603 35367.1020 2.989859 2000-2010 –> 2011-2019 5.710787 64914.1171 64995.5659
Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + phosphate_star 36091.9221 36126.2068 3.157881 2000-2010 –> 2011-2019 6.029758 66337.9834 66405.8574
Indo-Pacific (28.05,28.1] cstar_tref ~ temp + nitrate + phosphate_star 35748.8512 35783.1360 3.081684 2000-2010 –> 2011-2019 5.929920 65893.8830 65961.7570
Indo-Pacific (28.05,28.1] cstar_tref ~ temp + phosphate + phosphate_star 36193.5070 36227.7917 3.180802 2000-2010 –> 2011-2019 6.083320 66569.2758 66637.1498
Indo-Pacific (28.05,28.1] cstar_tref ~ temp + silicate + phosphate_star 35871.9904 35906.2751 3.108820 2000-2010 –> 2011-2019 5.871932 65646.1798 65714.0538
Indo-Pacific (28.1, Inf] cstar_tref ~ aou + phosphate 43530.6574 43558.5218 3.893162 1982-1999 –> 2000-2010 8.470916 85759.1679 85814.5485
Indo-Pacific (28.1, Inf] cstar_tref ~ aou + silicate + phosphate 43444.3603 43479.1908 3.871281 1982-1999 –> 2000-2010 8.449033 85674.8660 85744.0918
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + aou + phosphate 42961.1124 42995.9429 3.753687 1982-1999 –> 2000-2010 8.319643 85154.5671 85223.7928
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + aou + silicate + phosphate 42934.0742 42975.8708 3.746736 1982-1999 –> 2000-2010 8.312477 85128.8546 85211.9255
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + phosphate + phosphate_star 43457.2482 43492.0787 3.874467 1982-1999 –> 2000-2010 8.562838 86030.6275 86099.8533
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + silicate + phosphate + phosphate_star 43458.8905 43500.6871 3.874379 1982-1999 –> 2000-2010 8.561482 86030.3887 86113.4596
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + temp + phosphate_star 42798.4146 42833.2451 3.714905 1982-1999 –> 2000-2010 7.855687 83588.2762 83657.5020
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + temp + silicate + phosphate_star 42799.7526 42841.5492 3.714748 1982-1999 –> 2000-2010 7.792239 83370.4293 83453.5003
Indo-Pacific (28.1, Inf] cstar_tref ~ temp + phosphate_star 43041.8687 43069.7331 3.773568 1982-1999 –> 2000-2010 7.915856 83834.9513 83890.3319
Indo-Pacific (28.1, Inf] cstar_tref ~ temp + silicate + phosphate_star 43043.3917 43078.2222 3.773453 1982-1999 –> 2000-2010 7.855964 83629.7359 83698.9617
Indo-Pacific (28.1, Inf] cstar_tref ~ aou + phosphate 33707.0051 33734.2103 3.057175 2000-2010 –> 2011-2019 6.950336 77237.6625 77292.7322
Indo-Pacific (28.1, Inf] cstar_tref ~ aou + silicate + phosphate 33693.9837 33727.9903 3.053720 2000-2010 –> 2011-2019 6.925001 77138.3440 77207.1811
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + aou + phosphate_star 32853.1871 32887.1937 2.866455 2000-2010 –> 2011-2019 6.775199 76448.4199 76517.2570
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + nitrate + phosphate_star 32780.6773 32814.6839 2.850854 2000-2010 –> 2011-2019 6.835958 76679.0063 76747.8434
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + nitrate + silicate + phosphate_star 32691.6599 32732.4678 2.831390 2000-2010 –> 2011-2019 6.744829 76307.6996 76390.3041
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + phosphate + phosphate_star 32885.9399 32919.9465 2.873530 2000-2010 –> 2011-2019 6.747997 76343.1881 76412.0252
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + phosphate_star 32909.3657 32936.5709 2.879035 2000-2010 –> 2011-2019 6.911309 76990.0389 77045.1086
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + silicate + phosphate + phosphate_star 32523.4825 32564.2904 2.795776 2000-2010 –> 2011-2019 6.670154 75982.3729 76064.9775
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + silicate + phosphate_star 32701.2405 32735.2471 2.833859 2000-2010 –> 2011-2019 6.756823 76353.3605 76422.1976
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + temp + phosphate_star 32906.0255 32940.0321 2.877878 2000-2010 –> 2011-2019 6.592783 75704.4400 75773.2771

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(basin, gamma_slab, model, eras),
  eras_era)

lm_best <- left_join(
  lm_best,
  GLODAP_tidy %>% select(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(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 %>%
  ggplot(aes(as.factor(n_predictors),
             !!sym(params_local$MLR_criterion),
             col = basin)) +
  geom_hline(yintercept = c(0,10)) +
  geom_boxplot() +
  facet_grid(gamma_slab~era) +
  scale_color_brewer(palette = "Set1") +
  ylim(c(0,NA)) +
  labs(x="Number of predictors")

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

7.1.2 Best models

left_join(lm_best_target %>% select(basin, gamma_slab, era, model),
          GLODAP_glanced) %>%
  ggplot(aes("",
             !!sym(params_local$MLR_criterion),
             col = basin)) +
  geom_hline(yintercept = c(0, 10)) +
  geom_boxplot() +
  facet_grid(gamma_slab ~ era) +
  scale_color_brewer(palette = "Set1") +
  ylim(c(0, NA)) +
  labs(x = "Number of predictors pooled")

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

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(eras ~ basin)

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

Version Author Date
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(basin, eras, 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", "eras")) %>% 
  summary_rows(
    groups = TRUE,
    fns = list(total = "sum")
  )
aou nitrate phosphate phosphate_star sal silicate temp
Atlantic - 1982-1999 --> 2000-2010
(-Inf,26] 8 4 4 4 2 6 7
(26,26.5] 5 2 4 7 1 6 NA
(26.5,26.75] 3 2 4 4 4 2 1
(26.75,27] 4 3 3 4 4 2 NA
(27,27.25] 3 3 4 4 3 2 1
(27.25,27.5] 4 3 5 6 2 4 NA
(27.5,27.75] 3 2 6 6 2 5 1
(27.75,27.85] 4 4 4 3 3 1 2
(27.85,27.95] 3 3 3 4 3 NA 4
(27.95,28.05] 2 3 4 4 5 NA 5
(28.05,28.1] 2 2 2 3 4 2 5
(28.1,28.15] 3 2 2 2 4 2 5
(28.15,28.2] 1 1 1 1 1 1 6
(28.2, Inf] 4 1 5 6 1 1 4
total 49.00 35.00 51.00 58.00 39.00 — —
Atlantic - 2000-2010 --> 2011-2019
(-Inf,26] 7 2 6 5 1 6 8
(26,26.5] 6 3 3 6 1 5 NA
(26.5,26.75] 4 3 3 4 4 2 NA
(26.75,27] 4 3 3 4 4 2 NA
(27,27.25] 3 3 3 5 4 5 NA
(27.25,27.5] 4 3 5 6 2 4 NA
(27.5,27.75] 3 1 7 5 3 5 2
(27.75,27.85] 4 3 3 4 4 2 3
(27.85,27.95] 4 3 3 4 2 1 4
(27.95,28.05] 1 4 3 4 5 NA 5
(28.05,28.1] 1 3 2 5 4 2 3
(28.1,28.15] 3 2 2 2 4 2 5
(28.15,28.2] 1 1 1 1 1 1 6
(28.2, Inf] 5 2 4 6 NA 1 5
total 50.00 36.00 48.00 61.00 — — —
Indo-Pacific - 1982-1999 --> 2000-2010
(-Inf,26] 3 3 2 6 9 5 10
(26,26.5] 3 2 3 9 5 6 2
(26.5,26.75] 2 2 4 8 1 6 3
(26.75,27] 3 2 4 8 3 6 1
(27,27.25] 3 1 3 10 3 4 5
(27.25,27.5] 1 2 4 8 3 6 6
(27.5,27.75] 2 3 3 10 6 4 8
(27.75,27.85] 3 3 1 10 6 3 7
(27.85,27.95] 3 2 2 10 6 3 8
(27.95,28.05] 3 2 2 10 6 3 8
(28.05,28.1] 2 2 2 10 6 3 6
(28.1, Inf] 4 NA 6 6 6 5 4
total 32.00 — 36.00 105.00 60.00 54.00 68.00
Indo-Pacific - 2000-2010 --> 2011-2019
(-Inf,26] 4 2 2 5 9 6 10
(26,26.5] 4 NA 4 8 4 5 2
(26.5,26.75] 2 2 4 8 1 6 3
(26.75,27] 2 4 4 6 3 7 2
(27,27.25] 3 1 3 10 3 4 5
(27.25,27.5] 2 3 3 9 6 4 9
(27.5,27.75] 2 3 3 10 6 4 8
(27.75,27.85] 3 2 2 10 6 3 8
(27.85,27.95] 3 2 2 10 6 3 8
(27.95,28.05] 3 2 2 10 6 3 8
(28.05,28.1] 2 2 2 10 6 3 6
(28.1, Inf] 3 2 4 8 8 4 1
total 33.00 — 35.00 104.00 64.00 52.00 70.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(eras~basin)

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

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

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

lm_all_eras %>% 
  ggplot(aes(rmse, BIC, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  facet_grid(eras~basin)

Version Author Date
108d6b2 jens-daniel-mueller 2021-03-20 
lm_best_eras %>% 
  ggplot(aes(rmse, BIC, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  facet_grid(eras~basin)

Version Author Date
108d6b2 jens-daniel-mueller 2021-03-20

7.4.3 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(eras~basin)

Version Author Date
108d6b2 jens-daniel-mueller 2021-03-20 
lm_best_eras %>% 
  ggplot(aes(AIC, BIC, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  facet_grid(eras~basin)

Version Author Date
108d6b2 jens-daniel-mueller 2021-03-20

7.5 RMSE vs VIF

GLODAP_glanced %>% 
  ggplot(aes(rmse, log10(vif_max))) +
  geom_hline(yintercept = 1) +
  geom_point() +
  scale_color_viridis_d() +
  facet_grid(gamma_slab~basin)

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

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

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

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

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

p_residuals <- function(df){
    ggplot(data = df, aes(lat, .resid)) +
      geom_bin2d(binwidth = i_binwidth) +
      geom_hline(yintercept = 0, col = "black") +
      scale_fill_viridis_c() +
      facet_grid(gamma_slab ~ era) +
      coord_cartesian(ylim = i_ylim) +
      labs(
        title = paste(
          "Target variable:",
          params_local$MLR_target,
          "| selected best models | basin:",
          unique(df$basin)
        )
      )
}

GLODAP_augmented_best %>% 
  group_split(basin) %>% 
  map(p_residuals)
[[1]]

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

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(basin, gamma_slab, era, lat_grid) %>% 
  summarise(.resid_mean = mean(.resid)) %>% 
  ungroup()

lat_residual %>%
  ggplot(aes(lat_grid, .resid_mean, col=era)) +
  geom_line() +
  geom_point() +
  geom_hline(yintercept = 0, col = "black") +
  facet_grid(gamma_slab ~ basin)

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

7.6.5 Latitudinal offset

# calculate residual offset for adjacent eras
lat_residual_offset <- lat_residual  %>%
  select(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
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))) + 
  geom_hline(yintercept = log10(params_local$vif_max)) +
  geom_violin() +
  facet_grid(gamma_slab~basin)

Version Author Date
108d6b2 jens-daniel-mueller 2021-03-20
330dcd0 jens-daniel-mueller 2021-03-20 
GLODAP_glanced_vif %>% 
  ggplot(aes(Variables, log10(VIF))) + 
  geom_hline(yintercept = log10(params_local$vif_max)) +
  geom_violin() +
  # geom_point(shape = 21) +
  facet_grid(gamma_slab~basin) +
  coord_flip()

Version Author Date
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.2      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.2      purrr_0.3.4     
[17] readr_1.4.0      tidyr_1.1.2      tibble_3.0.4     ggplot2_3.3.2   
[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.9              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_0.2.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