Last updated: 2021-03-15

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: 9ebb3f2

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

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


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

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

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

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

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

# GLODAP_NA <- GLODAP %>% 
#   filter_all(any_vars(is.na(.)))

# 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)]) %>%
  drop_na() %>% 
  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)
  )

# 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)) %>%
  unnest(glanced) %>% 
  rename(n_predictors = n)

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

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

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

# 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 to glanced output
GLODAP_glanced <- full_join(GLODAP_glanced, GLODAP_glanced_rmse)
rm(GLODAP_glanced_rmse)


# extract input data
GLODAP_data <- GLODAP_nested_lm_fit %>%
  select(-c(fit, tidied, glanced, augmented)) %>%
  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 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  %>%
  select(basin, gamma_slab, model, era, AIC, 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)
         ) %>% 
  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)
} else {
  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 rmse eras rmse_sum aic_sum
Atlantic (-Inf,26] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 551.3828 6.956692 1982-1999 –> 2000-2012 12.949861 1123.846
Atlantic (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 545.1070 6.689105 1982-1999 –> 2000-2012 12.490040 1111.898
Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 546.4730 6.746459 1982-1999 –> 2000-2012 12.528415 1112.693
Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 546.6136 6.752389 1982-1999 –> 2000-2012 12.531024 1112.734
Atlantic (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 547.3751 6.784603 1982-1999 –> 2000-2012 12.555208 1113.254
Atlantic (-Inf,26] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 550.8536 6.933722 1982-1999 –> 2000-2012 12.720052 1117.206
Atlantic (-Inf,26] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 550.8148 6.932037 1982-1999 –> 2000-2012 12.874994 1121.814
Atlantic (-Inf,26] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 550.6738 6.925934 1982-1999 –> 2000-2012 12.848828 1121.085
Atlantic (-Inf,26] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 548.9192 6.936565 1982-1999 –> 2000-2012 12.959869 1120.255
Atlantic (-Inf,26] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 550.2292 6.906712 1982-1999 –> 2000-2012 12.793594 1119.579
Atlantic (-Inf,26] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 571.2704 6.945624 2000-2012 –> 2013-2019 13.902316 1122.653
Atlantic (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 573.6845 7.047372 2000-2012 –> 2013-2019 13.736477 1118.791
Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 573.8101 7.052707 2000-2012 –> 2013-2019 13.799166 1120.283
Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 573.8171 7.053006 2000-2012 –> 2013-2019 13.805395 1120.431
Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 568.9683 6.849967 2000-2012 –> 2013-2019 13.687451 1117.586
Atlantic (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 573.8865 7.055954 2000-2012 –> 2013-2019 13.840556 1121.262
Atlantic (-Inf,26] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 574.6824 7.089865 2000-2012 –> 2013-2019 14.023587 1125.536
Atlantic (-Inf,26] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 570.7391 6.923432 2000-2012 –> 2013-2019 13.855469 1121.554
Atlantic (-Inf,26] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 570.4927 6.913161 2000-2012 –> 2013-2019 13.839094 1121.166
Atlantic (-Inf,26] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 569.7742 6.883305 2000-2012 –> 2013-2019 13.790017 1120.003
Atlantic (26,26.5] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 4926.3364 5.491807 1982-1999 –> 2000-2012 10.372023 8092.671
Atlantic (26,26.5] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 4928.1618 5.491198 1982-1999 –> 2000-2012 10.348647 8091.586
Atlantic (26,26.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 4926.8474 5.486614 1982-1999 –> 2000-2012 10.361696 8094.076
Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 4964.2662 5.618611 1982-1999 –> 2000-2012 10.517020 8136.508
Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 4984.4990 5.691301 1982-1999 –> 2000-2012 10.616843 8162.540
Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 4980.7762 5.677856 1982-1999 –> 2000-2012 10.475307 8131.151
Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + phosphate 4962.9640 5.621102 1982-1999 –> 2000-2012 10.519619 8133.228
Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 4912.2038 5.435807 1982-1999 –> 2000-2012 10.244055 8064.939
Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 4961.5712 5.608999 1982-1999 –> 2000-2012 10.486236 8129.265
Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 4979.8917 5.674666 1982-1999 –> 2000-2012 10.579132 8153.431
Atlantic (26,26.5] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 1522.4898 7.047624 2000-2012 –> 2013-2019 12.539431 6448.826
Atlantic (26,26.5] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 1524.1465 7.042224 2000-2012 –> 2013-2019 12.533423 6452.308
Atlantic (26,26.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 1502.9817 6.717267 2000-2012 –> 2013-2019 12.203881 6429.829
Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 1504.5924 6.741461 2000-2012 –> 2013-2019 12.360072 6468.859
Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 1504.6511 6.742344 2000-2012 –> 2013-2019 12.433645 6489.150
Atlantic (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 1505.8737 6.760769 2000-2012 –> 2013-2019 12.560753 6520.147
Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + phosphate 1522.4038 7.046271 2000-2012 –> 2013-2019 12.667374 6485.368
Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 1524.3855 7.045984 2000-2012 –> 2013-2019 12.481791 6436.589
Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 1524.1085 7.041627 2000-2012 –> 2013-2019 12.650626 6485.680
Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1524.1346 7.042037 2000-2012 –> 2013-2019 12.716703 6504.026
Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 7227.9952 5.308155 1982-1999 –> 2000-2012 9.434228 11587.152
Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 7234.8360 5.323722 1982-1999 –> 2000-2012 9.233727 11511.591
Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 7224.8631 5.301043 1982-1999 –> 2000-2012 9.440930 11589.141
Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 7230.5719 5.314013 1982-1999 –> 2000-2012 9.410586 11578.737
Atlantic (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 7228.6503 5.309644 1982-1999 –> 2000-2012 9.421884 11582.663
Atlantic (26.5,26.75] cstar_tref ~ sal + temp + phosphate + phosphate_star 7233.3884 5.324981 1982-1999 –> 2000-2012 9.450866 11590.476
Atlantic (26.5,26.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 7231.8368 5.316892 1982-1999 –> 2000-2012 9.402748 11575.989
Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 7244.8285 5.346544 1982-1999 –> 2000-2012 9.310085 11542.417
Atlantic (26.5,26.75] cstar_tref ~ temp + aou + phosphate + phosphate_star 7242.9579 5.346840 1982-1999 –> 2000-2012 9.323692 11543.683
Atlantic (26.5,26.75] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 7239.2155 5.333713 1982-1999 –> 2000-2012 9.303822 11539.341
Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 1621.2571 7.394655 2000-2012 –> 2013-2019 12.714114 8854.221
Atlantic (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 1621.0266 7.391030 2000-2012 –> 2013-2019 12.699185 8849.022
Atlantic (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 1620.3396 7.380234 2000-2012 –> 2013-2019 12.689878 8848.990
Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + phosphate 1626.2281 7.505149 2000-2012 –> 2013-2019 12.860342 8872.832
Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 1623.5926 7.431491 2000-2012 –> 2013-2019 12.778035 8868.421
Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + phosphate_star 1625.7438 7.497419 2000-2012 –> 2013-2019 12.845224 8869.123
Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 1628.2261 7.505117 2000-2012 –> 2013-2019 12.838612 8867.346
Atlantic (26.5,26.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1627.7259 7.497135 2000-2012 –> 2013-2019 12.841276 8871.504
Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 1625.1556 7.488042 2000-2012 –> 2013-2019 12.837266 8869.155
Atlantic (26.5,26.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 1627.1415 7.487818 2000-2012 –> 2013-2019 12.834547 8872.051
Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 8327.1290 5.319844 1982-1999 –> 2000-2012 10.135093 13915.843
Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate 8350.1359 5.365538 1982-1999 –> 2000-2012 10.172686 13935.711
Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 8322.7843 5.311259 1982-1999 –> 2000-2012 10.120039 13908.992
Atlantic (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 8324.8201 5.315280 1982-1999 –> 2000-2012 10.107040 13904.419
Atlantic (26.75,27] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 8321.9431 5.309598 1982-1999 –> 2000-2012 10.119607 13908.627
Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 8336.1698 5.337753 1982-1999 –> 2000-2012 10.141884 13920.575
Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 8353.7955 5.372843 1982-1999 –> 2000-2012 10.174778 13937.348
Atlantic (26.75,27] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 8313.6376 5.293230 1982-1999 –> 2000-2012 10.132355 13911.571
Atlantic (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 8354.6361 5.374522 1982-1999 –> 2000-2012 10.175708 13937.898
Atlantic (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate_star 8352.8863 5.375022 1982-1999 –> 2000-2012 10.187272 13938.439
Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 2310.5578 6.034717 2000-2012 –> 2013-2019 11.354561 10637.687
Atlantic (26.75,27] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2307.1669 6.006125 2000-2012 –> 2013-2019 11.317384 10629.951
Atlantic (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 2309.6708 6.027225 2000-2012 –> 2013-2019 11.342505 10634.491
Atlantic (26.75,27] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 2307.0254 6.004935 2000-2012 –> 2013-2019 11.314533 10628.968
Atlantic (26.75,27] cstar_tref ~ sal + temp + silicate + phosphate_star 2310.1436 6.048135 2000-2012 –> 2013-2019 11.379757 10641.222
Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 2309.7054 6.027517 2000-2012 –> 2013-2019 11.365271 10645.875
Atlantic (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2305.0339 5.988209 2000-2012 –> 2013-2019 11.361052 10658.829
Atlantic (26.75,27] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 2305.1347 5.989055 2000-2012 –> 2013-2019 11.282285 10618.772
Atlantic (26.75,27] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 2304.6268 5.984796 2000-2012 –> 2013-2019 11.359318 10659.263
Atlantic (26.75,27] cstar_tref ~ temp + silicate + phosphate + phosphate_star 2305.1328 6.005838 2000-2012 –> 2013-2019 11.382737 10658.958
Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate 15410.9036 5.367043 1982-1999 –> 2000-2012 10.723249 29232.196
Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 15350.2636 5.301931 1982-1999 –> 2000-2012 10.621326 29140.812
Atlantic (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 15386.3591 5.340593 1982-1999 –> 2000-2012 10.632514 29153.823
Atlantic (27,27.25] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 15355.4180 5.307435 1982-1999 –> 2000-2012 10.625160 29144.567
Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 15394.3114 5.349148 1982-1999 –> 2000-2012 10.811791 29303.323
Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 15394.7461 5.349616 1982-1999 –> 2000-2012 10.736864 29241.801
Atlantic (27,27.25] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 15358.8162 5.311066 1982-1999 –> 2000-2012 10.752991 29250.888
Atlantic (27,27.25] cstar_tref ~ temp + aou + silicate + phosphate_star 15419.5374 5.378544 1982-1999 –> 2000-2012 10.842479 29327.604
Atlantic (27,27.25] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 15408.1701 5.364091 1982-1999 –> 2000-2012 10.758373 29261.042
Atlantic (27,27.25] cstar_tref ~ temp + silicate + phosphate + phosphate_star 15425.4540 5.384953 1982-1999 –> 2000-2012 10.848241 29332.993
Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 4118.4123 7.441818 2000-2012 –> 2013-2019 12.846536 19564.067
Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate 4132.9960 7.532963 2000-2012 –> 2013-2019 12.900006 19543.900
Atlantic (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4130.6524 7.518241 2000-2012 –> 2013-2019 12.820172 19480.916
Atlantic (27,27.25] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 4119.8882 7.450992 2000-2012 –> 2013-2019 12.874139 19582.455
Atlantic (27,27.25] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 4128.1529 7.502572 2000-2012 –> 2013-2019 12.843164 19514.512
Atlantic (27,27.25] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4130.6275 7.518085 2000-2012 –> 2013-2019 12.825519 19486.045
Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 4119.4564 7.448307 2000-2012 –> 2013-2019 12.877990 19588.007
Atlantic (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 4113.3763 7.410601 2000-2012 –> 2013-2019 12.760217 19508.122
Atlantic (27,27.25] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 4130.8419 7.519431 2000-2012 –> 2013-2019 12.830497 19489.658
Atlantic (27,27.25] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 4115.0700 7.421085 2000-2012 –> 2013-2019 12.785176 19523.240
Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 17379.8800 4.462806 1982-1999 –> 2000-2012 9.079380 31788.903
Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 17296.4637 4.400759 1982-1999 –> 2000-2012 9.010134 31697.868
Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 17263.0730 4.376165 1982-1999 –> 2000-2012 9.023495 31704.513
Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + phosphate 17322.1757 4.421276 1982-1999 –> 2000-2012 9.126105 31821.648
Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 17302.7779 4.405426 1982-1999 –> 2000-2012 9.105065 31798.861
Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 17320.2314 4.418350 1982-1999 –> 2000-2012 9.121243 31819.694
Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 17302.6551 4.405335 1982-1999 –> 2000-2012 9.140648 31835.656
Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 17269.5907 4.380955 1982-1999 –> 2000-2012 9.009391 31691.142
Atlantic (27.25,27.5] cstar_tref ~ sal + temp + phosphate + phosphate_star 17302.3257 4.406570 1982-1999 –> 2000-2012 9.137346 31828.647
Atlantic (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 17303.2046 4.405741 1982-1999 –> 2000-2012 9.132257 31827.128
Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 4639.7057 7.482702 2000-2012 –> 2013-2019 11.945508 22019.586
Atlantic (27.25,27.5] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 4642.0247 7.495585 2000-2012 –> 2013-2019 11.968398 22035.248
Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 4646.3694 7.519783 2000-2012 –> 2013-2019 11.920542 21942.833
Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 4660.7070 7.600192 2000-2012 –> 2013-2019 11.976357 21923.780
Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + phosphate 4649.8670 7.550514 2000-2012 –> 2013-2019 11.971789 21972.043
Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 4644.0729 7.506983 2000-2012 –> 2013-2019 11.912409 21946.851
Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 4648.3541 7.530863 2000-2012 –> 2013-2019 11.949213 21968.586
Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 4654.7501 7.566680 2000-2012 –> 2013-2019 11.947635 21924.341
Atlantic (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4655.9811 7.573593 2000-2012 –> 2013-2019 11.979334 21959.186
Atlantic (27.25,27.5] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 4636.9118 7.467209 2000-2012 –> 2013-2019 11.992651 22099.831
Atlantic (27.5,27.75] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 20432.6923 4.642043 1982-1999 –> 2000-2012 9.233040 34520.288
Atlantic (27.5,27.75] cstar_tref ~ aou + silicate + phosphate + phosphate_star 20441.6431 4.649404 1982-1999 –> 2000-2012 9.275707 34563.874
Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 20393.4336 4.615752 1982-1999 –> 2000-2012 9.232471 34507.747
Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 20385.9247 4.610741 1982-1999 –> 2000-2012 9.228220 34501.026
Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 20422.2204 4.635016 1982-1999 –> 2000-2012 9.278531 34564.210
Atlantic (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 20396.6099 4.617874 1982-1999 –> 2000-2012 9.248663 34525.476
Atlantic (27.5,27.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 20433.5434 4.642615 1982-1999 –> 2000-2012 9.236198 34523.833
Atlantic (27.5,27.75] cstar_tref ~ temp + aou + silicate + phosphate 20443.6371 4.650745 1982-1999 –> 2000-2012 9.279979 34568.897
Atlantic (27.5,27.75] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 20443.3552 4.649210 1982-1999 –> 2000-2012 9.275477 34567.549
Atlantic (27.5,27.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 20447.2716 4.651845 1982-1999 –> 2000-2012 9.260534 34553.261
Atlantic (27.5,27.75] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 5424.8712 6.426161 2000-2012 –> 2013-2019 11.068205 25857.564
Atlantic (27.5,27.75] cstar_tref ~ aou + silicate + phosphate + phosphate_star 5431.1017 6.458296 2000-2012 –> 2013-2019 11.107700 25872.745
Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 5431.8832 6.453528 2000-2012 –> 2013-2019 11.069281 25825.317
Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 5431.7555 6.453029 2000-2012 –> 2013-2019 11.063770 25817.680
Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5435.1755 6.466418 2000-2012 –> 2013-2019 11.101434 25857.396
Atlantic (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5432.0894 6.454335 2000-2012 –> 2013-2019 11.072209 25828.699
Atlantic (27.5,27.75] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 5426.4751 6.432411 2000-2012 –> 2013-2019 11.075026 25860.019
Atlantic (27.5,27.75] cstar_tref ~ temp + aou + silicate + phosphate 5431.1191 6.458364 2000-2012 –> 2013-2019 11.109110 25874.756
Atlantic (27.5,27.75] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 5433.0052 6.457918 2000-2012 –> 2013-2019 11.107129 25876.360
Atlantic (27.5,27.75] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 5428.1592 6.438980 2000-2012 –> 2013-2019 11.090825 25875.431
Atlantic (27.75,27.85] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 12968.7355 4.658827 1982-1999 –> 2000-2012 9.372098 19717.155
Atlantic (27.75,27.85] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 12946.6528 4.635397 1982-1999 –> 2000-2012 9.305332 19674.122
Atlantic (27.75,27.85] cstar_tref ~ sal + aou + silicate + phosphate 12988.2669 4.681785 1982-1999 –> 2000-2012 9.377975 19726.452
Atlantic (27.75,27.85] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 12990.0510 4.681554 1982-1999 –> 2000-2012 9.376672 19729.718
Atlantic (27.75,27.85] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 12948.3154 4.637157 1982-1999 –> 2000-2012 9.315025 19679.634
Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 12990.0574 4.681561 1982-1999 –> 2000-2012 9.373960 19728.411
Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 12990.4450 4.681976 1982-1999 –> 2000-2012 9.384830 19733.846
Atlantic (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 12964.1778 4.653981 1982-1999 –> 2000-2012 9.368430 19713.164
Atlantic (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 12957.6699 4.647072 1982-1999 –> 2000-2012 9.372295 19711.834
Atlantic (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 12961.3246 4.650951 1982-1999 –> 2000-2012 9.368708 19711.902
Atlantic (27.75,27.85] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 3339.7700 7.788418 2000-2012 –> 2013-2019 12.447244 16308.506
Atlantic (27.75,27.85] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 3338.0443 7.774400 2000-2012 –> 2013-2019 12.409797 16284.697
Atlantic (27.75,27.85] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 3329.6820 7.706833 2000-2012 –> 2013-2019 12.456525 16383.021
Atlantic (27.75,27.85] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 3337.8066 7.772471 2000-2012 –> 2013-2019 12.454026 16327.858
Atlantic (27.75,27.85] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 3336.0194 7.757985 2000-2012 –> 2013-2019 12.395142 16284.335
Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 3337.3621 7.768866 2000-2012 –> 2013-2019 12.455188 16331.871
Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 3340.6678 7.795720 2000-2012 –> 2013-2019 12.477695 16331.113
Atlantic (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 3339.4428 7.785758 2000-2012 –> 2013-2019 12.439739 16303.621
Atlantic (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 3335.7672 7.755943 2000-2012 –> 2013-2019 12.403014 16293.437
Atlantic (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 3338.8527 7.780963 2000-2012 –> 2013-2019 12.431914 16300.177
Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 27061.8568 4.892557 1982-1999 –> 2000-2012 9.014203 40616.729
Atlantic (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 26939.4156 4.826418 1982-1999 –> 2000-2012 8.989562 40542.134
Atlantic (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate_star 27059.8601 4.892559 1982-1999 –> 2000-2012 9.056278 40661.238
Atlantic (27.85,27.95] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 27430.1663 5.097023 1982-1999 –> 2000-2012 9.284551 41060.777
Atlantic (27.85,27.95] cstar_tref ~ sal + silicate + phosphate + phosphate_star 27442.1688 5.104963 1982-1999 –> 2000-2012 9.318910 41100.817
Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 26988.5792 4.852867 1982-1999 –> 2000-2012 9.027603 40604.578
Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 26951.0319 4.832654 1982-1999 –> 2000-2012 8.996088 40554.083
Atlantic (27.85,27.95] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 27088.0490 4.906823 1982-1999 –> 2000-2012 9.102982 40728.493
Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate 27094.7128 4.911551 1982-1999 –> 2000-2012 9.239820 40881.204
Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 26980.9379 4.848746 1982-1999 –> 2000-2012 9.018271 40590.971
Atlantic (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 8460.3794 5.412077 2000-2012 –> 2013-2019 10.304635 35522.236
Atlantic (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 8437.7756 5.367255 2000-2012 –> 2013-2019 10.193673 35377.191
Atlantic (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate_star 8475.3291 5.445933 2000-2012 –> 2013-2019 10.338492 35535.189
Atlantic (27.85,27.95] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 8574.0050 5.643124 2000-2012 –> 2013-2019 10.740147 36004.171
Atlantic (27.85,27.95] cstar_tref ~ sal + silicate + phosphate + phosphate_star 8595.2487 5.691590 2000-2012 –> 2013-2019 10.796553 36037.417
Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 8442.3888 5.376373 2000-2012 –> 2013-2019 10.229240 35430.968
Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 8439.0110 5.369695 2000-2012 –> 2013-2019 10.202350 35390.043
Atlantic (27.85,27.95] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 8404.6131 5.302167 2000-2012 –> 2013-2019 10.208990 35492.662
Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate 8442.3584 5.380270 2000-2012 –> 2013-2019 10.291821 35537.071
Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 8442.5932 5.376777 2000-2012 –> 2013-2019 10.225523 35423.531
Atlantic (27.95,28.05] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 31113.3234 5.972570 1982-1999 –> 2000-2012 10.070028 49189.609
Atlantic (27.95,28.05] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 30731.1502 5.741830 1982-1999 –> 2000-2012 9.633991 48479.284
Atlantic (27.95,28.05] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 30755.0680 5.756005 1982-1999 –> 2000-2012 9.651009 48507.862
Atlantic (27.95,28.05] cstar_tref ~ sal + aou + silicate + phosphate_star 30758.0485 5.758962 1982-1999 –> 2000-2012 9.654128 48509.110
Atlantic (27.95,28.05] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 30837.5121 5.805136 1982-1999 –> 2000-2012 9.754426 48678.668
Atlantic (27.95,28.05] cstar_tref ~ sal + silicate + phosphate + phosphate_star 30867.7868 5.824484 1982-1999 –> 2000-2012 9.774370 48707.906
Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 30894.1339 5.839122 1982-1999 –> 2000-2012 9.933170 48965.104
Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 30814.9759 5.791665 1982-1999 –> 2000-2012 9.715586 48614.991
Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 30759.4505 5.758607 1982-1999 –> 2000-2012 9.653577 48512.190
Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 30791.5917 5.777720 1982-1999 –> 2000-2012 9.687720 48568.919
Atlantic (27.95,28.05] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 8629.5932 6.550409 2000-2012 –> 2013-2019 12.292239 39360.743
Atlantic (27.95,28.05] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 8636.6235 6.568063 2000-2012 –> 2013-2019 12.324069 39391.691
Atlantic (27.95,28.05] cstar_tref ~ sal + aou + silicate + phosphate_star 8658.1281 6.627434 2000-2012 –> 2013-2019 12.386396 39416.177
Atlantic (27.95,28.05] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 8683.4518 6.686878 2000-2012 –> 2013-2019 12.492015 39520.964
Atlantic (27.95,28.05] cstar_tref ~ sal + silicate + phosphate + phosphate_star 8715.2185 6.773884 2000-2012 –> 2013-2019 12.598368 39583.005
Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 8582.7697 6.434031 2000-2012 –> 2013-2019 12.273153 39476.904
Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 8677.2604 6.671047 2000-2012 –> 2013-2019 12.462712 39492.236
Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 8650.6494 6.603427 2000-2012 –> 2013-2019 12.362034 39410.100
Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate 8619.3942 6.524882 2000-2012 –> 2013-2019 12.517384 39765.037
Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 8671.4462 6.656214 2000-2012 –> 2013-2019 12.433933 39463.038
Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 23981.3340 6.128209 1982-1999 –> 2000-2012 10.054164 36350.888
Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 23723.6882 5.918960 1982-1999 –> 2000-2012 9.756201 35991.899
Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate_star 23981.7529 6.130209 1982-1999 –> 2000-2012 10.059846 36353.465
Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 23935.2677 6.090261 1982-1999 –> 2000-2012 10.019716 36308.774
Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 23987.2462 6.133097 1982-1999 –> 2000-2012 10.052836 36349.775
Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 23966.5722 6.116023 1982-1999 –> 2000-2012 10.042033 36336.188
Atlantic (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 23977.8897 6.125364 1982-1999 –> 2000-2012 10.054015 36350.489
Atlantic (28.05,28.1] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 23827.7283 6.002583 1982-1999 –> 2000-2012 9.944591 36215.376
Atlantic (28.05,28.1] cstar_tref ~ temp + aou + phosphate + phosphate_star 23825.7455 6.002597 1982-1999 –> 2000-2012 9.947380 36214.514
Atlantic (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 23710.9470 5.908799 1982-1999 –> 2000-2012 9.836881 36082.903
Atlantic (28.05,28.1] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 6976.3119 6.990116 2000-2012 –> 2013-2019 13.171076 31021.208
Atlantic (28.05,28.1] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 6988.4359 7.031178 2000-2012 –> 2013-2019 13.147921 30955.882
Atlantic (28.05,28.1] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 6988.6813 7.032011 2000-2012 –> 2013-2019 13.160731 30970.632
Atlantic (28.05,28.1] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 6983.0121 7.012778 2000-2012 –> 2013-2019 13.142129 30965.727
Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 6874.6228 6.655024 2000-2012 –> 2013-2019 12.573984 30598.311
Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 6981.4876 7.007615 2000-2012 –> 2013-2019 13.140712 30968.734
Atlantic (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 6988.0171 7.029755 2000-2012 –> 2013-2019 13.145778 30954.589
Atlantic (28.05,28.1] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 6933.9581 6.848546 2000-2012 –> 2013-2019 12.851129 30761.686
Atlantic (28.05,28.1] cstar_tref ~ temp + aou + phosphate + phosphate_star 6949.9622 6.908372 2000-2012 –> 2013-2019 12.910969 30775.708
Atlantic (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 6912.3373 6.777387 2000-2012 –> 2013-2019 12.686186 30623.284
Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 18122.3755 6.710673 1982-1999 –> 2000-2012 10.455233 26512.980
Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 18134.0548 6.725069 1982-1999 –> 2000-2012 10.441978 26501.994
Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 18146.3462 6.740253 1982-1999 –> 2000-2012 10.485556 26537.557
Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 18134.0260 6.725033 1982-1999 –> 2000-2012 10.451163 26509.542
Atlantic (28.1,28.15] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 18137.3398 6.729124 1982-1999 –> 2000-2012 10.446109 26505.342
Atlantic (28.1,28.15] cstar_tref ~ sal + temp + silicate + phosphate_star 18139.0293 6.733681 1982-1999 –> 2000-2012 10.462810 26515.006
Atlantic (28.1,28.15] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 18155.4901 6.751571 1982-1999 –> 2000-2012 10.483590 26535.835
Atlantic (28.1,28.15] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 18131.6128 6.722056 1982-1999 –> 2000-2012 10.454467 26512.279
Atlantic (28.1,28.15] cstar_tref ~ temp + aou + silicate + phosphate_star 18154.1467 6.752385 1982-1999 –> 2000-2012 10.485498 26533.389
Atlantic (28.1,28.15] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 18157.0222 6.753469 1982-1999 –> 2000-2012 10.484824 26536.823
Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 5034.9358 8.431169 2000-2012 –> 2013-2019 15.174695 23183.928
Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 4992.3304 8.180918 2000-2012 –> 2013-2019 14.891590 23114.706
Atlantic (28.1,28.15] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 5021.8013 8.353215 2000-2012 –> 2013-2019 15.093469 23168.148
Atlantic (28.1,28.15] cstar_tref ~ sal + temp + phosphate + phosphate_star 5035.4900 8.446413 2000-2012 –> 2013-2019 15.188967 23181.696
Atlantic (28.1,28.15] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5037.4748 8.446322 2000-2012 –> 2013-2019 15.175445 23174.815
Atlantic (28.1,28.15] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 4984.1941 8.133979 2000-2012 –> 2013-2019 14.884780 23139.062
Atlantic (28.1,28.15] cstar_tref ~ temp + aou + phosphate + phosphate_star 4990.5398 8.182129 2000-2012 –> 2013-2019 14.932943 23143.418
Atlantic (28.1,28.15] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 4982.4420 8.123906 2000-2012 –> 2013-2019 14.845962 23114.055
Atlantic (28.1,28.15] cstar_tref ~ temp + nitrate + phosphate + phosphate_star 5025.5781 8.387412 2000-2012 –> 2013-2019 15.162740 23198.212
Atlantic (28.1,28.15] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 5026.3756 8.380282 2000-2012 –> 2013-2019 15.133751 23183.398
Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate + phosphate 6370.1568 3.693912 1982-1999 –> 2000-2012 7.842771 9544.207
Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 6369.2041 3.692404 1982-1999 –> 2000-2012 7.861733 9548.728
Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate 6351.2686 3.664114 1982-1999 –> 2000-2012 7.815332 9525.951
Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 6365.4173 3.686413 1982-1999 –> 2000-2012 7.852697 9544.128
Atlantic (28.15,28.2] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 6359.4818 3.677042 1982-1999 –> 2000-2012 7.837180 9536.551
Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 6377.4801 3.705531 1982-1999 –> 2000-2012 7.845315 9549.095
Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 6361.8957 3.680850 1982-1999 –> 2000-2012 7.819058 9533.088
Atlantic (28.15,28.2] cstar_tref ~ temp + aou + silicate + phosphate 6360.4177 3.681674 1982-1999 –> 2000-2012 7.838346 9534.560
Atlantic (28.15,28.2] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 6361.2563 3.679841 1982-1999 –> 2000-2012 7.834949 9536.980
Atlantic (28.15,28.2] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 6375.8353 3.702918 1982-1999 –> 2000-2012 7.853122 9550.246
Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate + phosphate 1015.1912 2.781502 2000-2012 –> 2013-2019 6.475415 7385.348
Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 1015.8571 2.786024 2000-2012 –> 2013-2019 6.512187 7406.286
Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + phosphate 1013.2580 2.781956 2000-2012 –> 2013-2019 6.476281 7381.676
Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 1015.1012 2.780892 2000-2012 –> 2013-2019 6.473296 7384.305
Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + phosphate_star 1014.0305 2.787202 2000-2012 –> 2013-2019 6.513884 7402.784
Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate 1012.8285 2.765520 2000-2012 –> 2013-2019 6.429634 7364.097
Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1015.2607 2.781974 2000-2012 –> 2013-2019 6.468387 7380.678
Atlantic (28.15,28.2] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 1014.9989 2.780198 2000-2012 –> 2013-2019 6.498623 7400.580
Atlantic (28.15,28.2] cstar_tref ~ sal + temp + phosphate + phosphate_star 1013.1182 2.781007 2000-2012 –> 2013-2019 6.500966 7397.661
Atlantic (28.15,28.2] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 1013.5425 2.770340 2000-2012 –> 2013-2019 6.447382 7373.024
Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 14237.7674 3.309133 1982-1999 –> 2000-2012 6.690135 20401.510
Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 14305.3751 3.350530 1982-1999 –> 2000-2012 6.754193 20484.695
Atlantic (28.2, Inf] cstar_tref ~ sal + aou + phosphate + phosphate_star 14245.1575 3.314852 1982-1999 –> 2000-2012 6.710913 20417.263
Atlantic (28.2, Inf] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 14236.7928 3.308540 1982-1999 –> 2000-2012 6.600847 20338.542
Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate 14265.1728 3.325852 1982-1999 –> 2000-2012 6.735142 20448.345
Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + phosphate 14276.7743 3.334181 1982-1999 –> 2000-2012 6.758639 20468.299
Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 14233.2950 3.306413 1982-1999 –> 2000-2012 6.700203 20405.841
Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + silicate + phosphate 14275.7526 3.332329 1982-1999 –> 2000-2012 6.612985 20369.235
Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 14208.0305 3.291087 1982-1999 –> 2000-2012 6.608718 20327.648
Atlantic (28.2, Inf] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 14218.4793 3.297417 1982-1999 –> 2000-2012 6.592799 20322.406
Atlantic (28.2, Inf] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2017.7709 3.200287 2000-2012 –> 2013-2019 6.508827 16254.564
Atlantic (28.2, Inf] cstar_tref ~ sal + nitrate + silicate + phosphate + phosphate_star 2022.2780 3.218929 2000-2012 –> 2013-2019 6.527016 16258.326
Atlantic (28.2, Inf] cstar_tref ~ sal + silicate + phosphate + phosphate_star 2021.0618 3.222182 2000-2012 –> 2013-2019 6.540228 16271.457
Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 2005.9330 3.151837 2000-2012 –> 2013-2019 6.457143 16237.408
Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 2003.8191 3.143263 2000-2012 –> 2013-2019 6.449675 16237.114
Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + phosphate_star 2003.9511 3.151911 2000-2012 –> 2013-2019 6.459794 16237.665
Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1999.9462 3.127614 2000-2012 –> 2013-2019 6.418701 16207.977
Atlantic (28.2, Inf] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 2012.4504 3.178420 2000-2012 –> 2013-2019 6.482114 16241.273
Atlantic (28.2, Inf] cstar_tref ~ sal + temp + phosphate + phosphate_star 2010.5849 3.178971 2000-2012 –> 2013-2019 6.486337 16243.449
Atlantic (28.2, Inf] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 2005.5741 3.150379 2000-2012 –> 2013-2019 6.447796 16224.053
Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 29168.2678 5.954357 1982-1999 –> 2000-2012 17.762836 42066.292
Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 29642.9152 6.273101 1982-1999 –> 2000-2012 17.786860 42457.178
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 29244.7951 6.004631 1982-1999 –> 2000-2012 17.905852 42168.746
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 29259.3854 6.014264 1982-1999 –> 2000-2012 17.943918 42191.244
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 29397.2145 6.106029 1982-1999 –> 2000-2012 18.537608 42465.652
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 29819.9948 6.396339 1982-1999 –> 2000-2012 17.975202 42652.944
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 29876.7422 6.436342 1982-1999 –> 2000-2012 18.044124 42717.958
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 29335.9395 6.065061 1982-1999 –> 2000-2012 18.025167 42276.247
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 30186.1345 6.658885 1982-1999 –> 2000-2012 18.395757 43064.002
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 29902.1779 6.454354 1982-1999 –> 2000-2012 18.075543 42747.219
Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 16975.7498 4.538412 2000-2012 –> 2013-2019 10.492769 46144.018
Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 18071.3307 5.484508 2000-2012 –> 2013-2019 11.757609 47714.246
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate 17400.7481 4.886008 2000-2012 –> 2013-2019 10.993399 46797.992
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate 17282.2002 4.785264 2000-2012 –> 2013-2019 10.789895 46526.995
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 17307.6919 4.806394 2000-2012 –> 2013-2019 10.820657 46567.077
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 17394.6862 4.879205 2000-2012 –> 2013-2019 10.985234 46791.901
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 18311.7731 5.717224 2000-2012 –> 2013-2019 12.113563 48131.768
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 17470.1085 4.943224 2000-2012 –> 2013-2019 11.008285 46806.048
Indo-Pacific (-Inf,26] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 18026.9773 5.442627 2000-2012 –> 2013-2019 12.068540 48167.931
Indo-Pacific (-Inf,26] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 17978.7005 5.397404 2000-2012 –> 2013-2019 12.081109 48198.697
Indo-Pacific (26,26.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 18689.7484 4.849007 1982-1999 –> 2000-2012 12.324833 32618.315
Indo-Pacific (26,26.5] cstar_tref ~ sal + aou + phosphate + phosphate_star 19017.8778 5.112887 1982-1999 –> 2000-2012 12.963653 33142.932
Indo-Pacific (26,26.5] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 18735.5911 4.884819 1982-1999 –> 2000-2012 12.495236 32736.531
Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 18647.3511 4.816120 1982-1999 –> 2000-2012 12.436670 32653.688
Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 18666.8858 4.831245 1982-1999 –> 2000-2012 12.534220 32716.857
Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 18981.1284 5.081184 1982-1999 –> 2000-2012 12.914661 33099.239
Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 18772.2465 4.913645 1982-1999 –> 2000-2012 12.707695 32869.892
Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 18813.1659 4.946024 1982-1999 –> 2000-2012 12.831822 32958.277
Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 18672.5949 4.835674 1982-1999 –> 2000-2012 12.566669 32737.293
Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 18814.1970 4.946843 1982-1999 –> 2000-2012 12.852198 32969.355
Indo-Pacific (26,26.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 8982.8028 4.186732 2000-2012 –> 2013-2019 9.035739 27672.551
Indo-Pacific (26,26.5] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 8908.0090 4.088370 2000-2012 –> 2013-2019 8.973189 27643.600
Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 9037.5941 4.260288 2000-2012 –> 2013-2019 9.076408 27684.945
Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 9064.7680 4.297246 2000-2012 –> 2013-2019 9.128491 27731.654
Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 8989.4175 4.195545 2000-2012 –> 2013-2019 9.109189 27761.664
Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 9027.0885 4.246085 2000-2012 –> 2013-2019 9.192110 27840.254
Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 9069.1093 4.303180 2000-2012 –> 2013-2019 9.138854 27741.704
Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 8978.3858 4.180858 2000-2012 –> 2013-2019 9.099365 27756.794
Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 9026.6011 4.245427 2000-2012 –> 2013-2019 9.192270 27840.798
Indo-Pacific (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate_star 9025.7912 4.247034 2000-2012 –> 2013-2019 9.198671 27844.023
Indo-Pacific (26.5,26.75] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 20791.8548 5.803199 1982-1999 –> 2000-2012 12.068817 33898.015
Indo-Pacific (26.5,26.75] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 20786.0513 5.796279 1982-1999 –> 2000-2012 12.025528 33870.785
Indo-Pacific (26.5,26.75] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 20326.3631 5.402854 1982-1999 –> 2000-2012 11.248783 33155.532
Indo-Pacific (26.5,26.75] cstar_tref ~ sal + aou + phosphate + phosphate_star 20644.9292 5.674279 1982-1999 –> 2000-2012 12.089490 33846.033
Indo-Pacific (26.5,26.75] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 20646.2812 5.673717 1982-1999 –> 2000-2012 11.880293 33716.342
Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 20406.8564 5.469763 1982-1999 –> 2000-2012 11.533415 33383.170
Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 20460.7852 5.515053 1982-1999 –> 2000-2012 11.695405 33513.808
Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 20644.2775 5.671979 1982-1999 –> 2000-2012 12.058038 33829.054
Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 20467.3240 5.520570 1982-1999 –> 2000-2012 11.737478 33544.078
Indo-Pacific (26.5,26.75] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 20746.5129 5.761342 1982-1999 –> 2000-2012 11.946246 33802.498
Indo-Pacific (26.5,26.75] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 8306.2948 4.979050 2000-2012 –> 2013-2019 10.775328 29092.346
Indo-Pacific (26.5,26.75] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 8259.6354 4.895040 2000-2012 –> 2013-2019 10.297894 28585.998
Indo-Pacific (26.5,26.75] cstar_tref ~ sal + aou + phosphate + phosphate_star 8381.5421 5.121313 2000-2012 –> 2013-2019 10.795593 29026.471
Indo-Pacific (26.5,26.75] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 8337.0319 5.035177 2000-2012 –> 2013-2019 10.708895 28983.313
Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 8281.4896 4.934210 2000-2012 –> 2013-2019 10.403973 28688.346
Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 8300.0972 4.967808 2000-2012 –> 2013-2019 10.482861 28760.882
Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 8380.5012 5.115637 2000-2012 –> 2013-2019 10.787616 29024.779
Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 8301.0350 4.969508 2000-2012 –> 2013-2019 10.490077 28768.359
Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + phosphate_star 8303.5489 4.977695 2000-2012 –> 2013-2019 10.548298 28827.879
Indo-Pacific (26.5,26.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 8283.1036 4.937115 2000-2012 –> 2013-2019 10.452955 28744.822
Indo-Pacific (26.75,27] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 30930.6695 5.521128 1982-1999 –> 2000-2012 12.557209 52234.265
Indo-Pacific (26.75,27] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 30795.8581 5.445248 1982-1999 –> 2000-2012 12.447513 52071.016
Indo-Pacific (26.75,27] cstar_tref ~ aou + silicate + phosphate + phosphate_star 31022.8192 5.572832 1982-1999 –> 2000-2012 12.672678 52383.415
Indo-Pacific (26.75,27] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 30825.2347 5.461453 1982-1999 –> 2000-2012 12.432941 52072.562
Indo-Pacific (26.75,27] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 30967.3988 5.540558 1982-1999 –> 2000-2012 12.555416 52253.909
Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 30946.9493 5.529109 1982-1999 –> 2000-2012 12.547358 52236.513
Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 30755.8512 5.423256 1982-1999 –> 2000-2012 12.432799 52037.572
Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 30896.8589 5.501164 1982-1999 –> 2000-2012 12.640563 52294.554
Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + phosphate + phosphate_star 30949.9641 5.531914 1982-1999 –> 2000-2012 12.577797 52262.356
Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 30841.3084 5.470339 1982-1999 –> 2000-2012 12.431318 52079.104
Indo-Pacific (26.75,27] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 13601.7816 4.408879 2000-2012 –> 2013-2019 9.930007 44532.451
Indo-Pacific (26.75,27] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 13594.4071 4.400061 2000-2012 –> 2013-2019 9.845308 44390.265
Indo-Pacific (26.75,27] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 13563.6016 4.371205 2000-2012 –> 2013-2019 9.832658 44388.836
Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + nitrate + phosphate 13709.2279 4.509301 2000-2012 –> 2013-2019 10.131270 44820.833
Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 13768.4355 4.566687 2000-2012 –> 2013-2019 10.095795 44715.385
Indo-Pacific (26.75,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 13825.0487 4.622241 2000-2012 –> 2013-2019 10.142543 44756.239
Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 13541.6532 4.350762 2000-2012 –> 2013-2019 9.774018 44297.504
Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 13724.9285 4.524448 2000-2012 –> 2013-2019 10.025612 44621.787
Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + phosphate + phosphate_star 13767.6126 4.567835 2000-2012 –> 2013-2019 10.099749 44717.577
Indo-Pacific (26.75,27] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 13769.2507 4.567482 2000-2012 –> 2013-2019 10.037822 44610.559
Indo-Pacific (27,27.25] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 40417.2524 5.127728 1982-1999 –> 2000-2012 11.789096 65897.106
Indo-Pacific (27,27.25] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 40172.1648 5.032864 1982-1999 –> 2000-2012 11.692365 65651.864
Indo-Pacific (27,27.25] cstar_tref ~ aou + silicate + phosphate + phosphate_star 40311.2670 5.086820 1982-1999 –> 2000-2012 11.814468 65867.177
Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 40108.5859 5.008740 1982-1999 –> 2000-2012 11.669854 65590.147
Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 39815.6209 4.899062 1982-1999 –> 2000-2012 11.626630 65373.440
Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 39769.1375 4.881882 1982-1999 –> 2000-2012 11.428764 65117.817
Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate 39871.4009 4.919758 1982-1999 –> 2000-2012 11.786820 65586.876
Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 39791.8636 4.890274 1982-1999 –> 2000-2012 11.444860 65149.577
Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + phosphate + phosphate_star 40090.0335 5.002478 1982-1999 –> 2000-2012 11.567064 65457.457
Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 39604.6733 4.821579 1982-1999 –> 2000-2012 11.384086 64971.664
Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 18293.7823 4.179346 2000-2012 –> 2013-2019 9.188085 58402.368
Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 18424.6960 4.265498 2000-2012 –> 2013-2019 9.190467 58310.105
Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate 18424.2810 4.266552 2000-2012 –> 2013-2019 9.191525 58307.702
Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 18262.2201 4.158837 2000-2012 –> 2013-2019 9.057899 58077.841
Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 18240.6362 4.144870 2000-2012 –> 2013-2019 9.026751 58009.774
Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate 18302.8863 4.185280 2000-2012 –> 2013-2019 9.105038 58174.287
Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 18362.2061 4.224155 2000-2012 –> 2013-2019 9.201795 58388.376
Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 18355.8927 4.220001 2000-2012 –> 2013-2019 9.186942 58353.593
Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 18282.9195 4.172276 2000-2012 –> 2013-2019 9.062549 58074.783
Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 18260.0353 4.157421 2000-2012 –> 2013-2019 8.979000 57864.709
Indo-Pacific (27.25,27.5] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 41257.8873 4.624698 1982-1999 –> 2000-2012 10.006976 66874.211
Indo-Pacific (27.25,27.5] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 40990.4809 4.536429 1982-1999 –> 2000-2012 9.916879 66606.000
Indo-Pacific (27.25,27.5] cstar_tref ~ aou + silicate + phosphate + phosphate_star 41208.2783 4.608316 1982-1999 –> 2000-2012 10.015164 66862.195
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 40982.0591 4.533697 1982-1999 –> 2000-2012 9.898817 66574.026
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 40597.6238 4.410724 1982-1999 –> 2000-2012 9.803217 66231.596
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 40691.9051 4.440570 1982-1999 –> 2000-2012 9.589859 65944.963
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 40739.5914 4.455743 1982-1999 –> 2000-2012 10.081926 66723.726
Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 40498.9732 4.379709 1982-1999 –> 2000-2012 9.533398 65759.081
Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + phosphate + phosphate_star 40807.2076 4.477987 1982-1999 –> 2000-2012 9.636669 66073.308
Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 40146.2507 4.270589 1982-1999 –> 2000-2012 9.412224 65387.031
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 19866.2097 3.609297 2000-2012 –> 2013-2019 8.142995 60848.269
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 19739.0912 3.547378 2000-2012 –> 2013-2019 7.958102 60336.715
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate 19891.9270 3.621955 2000-2012 –> 2013-2019 8.242534 61139.356
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 19967.8397 3.659578 2000-2012 –> 2013-2019 8.100148 60659.745
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 19780.2285 3.567299 2000-2012 –> 2013-2019 8.023042 60519.820
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 19819.0219 3.586187 2000-2012 –> 2013-2019 8.092543 60716.516
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 19718.9772 3.537678 2000-2012 –> 2013-2019 8.005542 60496.545
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 19817.4305 3.585410 2000-2012 –> 2013-2019 8.084764 60693.186
Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 20109.4369 3.730802 2000-2012 –> 2013-2019 8.110511 60608.410
Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 19967.4158 3.659367 2000-2012 –> 2013-2019 7.929956 60113.666
Indo-Pacific (27.5,27.75] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 46210.8869 4.209858 1982-1999 –> 2000-2012 9.271791 73370.140
Indo-Pacific (27.5,27.75] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 46204.5955 4.207700 1982-1999 –> 2000-2012 9.269607 73365.801
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 45802.2700 4.104326 1982-1999 –> 2000-2012 9.153421 72940.853
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + phosphate + phosphate_star 46092.0923 4.179051 1982-1999 –> 2000-2012 9.236512 73243.455
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 45965.0368 4.145838 1982-1999 –> 2000-2012 9.198302 73109.573
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 45676.9295 4.072642 1982-1999 –> 2000-2012 9.281618 73093.834
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 46082.7497 4.176121 1982-1999 –> 2000-2012 9.060621 72925.438
Indo-Pacific (27.5,27.75] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 46142.2615 4.191515 1982-1999 –> 2000-2012 9.058099 72952.142
Indo-Pacific (27.5,27.75] cstar_tref ~ temp + aou + phosphate + phosphate_star 46469.3092 4.277662 1982-1999 –> 2000-2012 9.162308 73310.265
Indo-Pacific (27.5,27.75] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 46471.2765 4.277653 1982-1999 –> 2000-2012 9.162298 73314.229
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate 27030.3776 3.708175 2000-2012 –> 2013-2019 7.780817 72707.307
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 26854.3882 3.642825 2000-2012 –> 2013-2019 7.714315 72526.741
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 26728.8955 3.596931 2000-2012 –> 2013-2019 7.773052 72811.645
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 26992.9458 3.694178 2000-2012 –> 2013-2019 7.853671 73011.166
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 26835.4133 3.635849 2000-2012 –> 2013-2019 7.832189 72996.282
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + phosphate + phosphate_star 26836.7584 3.636343 2000-2012 –> 2013-2019 7.709955 72517.541
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 26882.9573 3.653355 2000-2012 –> 2013-2019 7.824900 72947.974
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + phosphate + phosphate_star 26862.6871 3.646618 2000-2012 –> 2013-2019 7.868434 73119.453
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 26830.6532 3.634101 2000-2012 –> 2013-2019 7.831629 72996.100
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate_star 26884.0437 3.654495 2000-2012 –> 2013-2019 7.857794 73069.715
Indo-Pacific (27.75,27.85] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 21736.6889 3.582942 1982-1999 –> 2000-2012 7.978984 33656.624
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 21755.1927 3.591177 1982-1999 –> 2000-2012 7.988789 33676.594
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 21681.5612 3.558519 1982-1999 –> 2000-2012 7.978296 33623.606
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 21639.1687 3.539852 1982-1999 –> 2000-2012 7.975763 33596.175
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 21622.7633 3.532654 1982-1999 –> 2000-2012 7.982714 33592.844
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 21637.1918 3.538984 1982-1999 –> 2000-2012 7.978887 33597.890
Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 21668.5837 3.552794 1982-1999 –> 2000-2012 7.974970 33612.856
Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 21632.8075 3.537059 1982-1999 –> 2000-2012 7.978306 33594.748
Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 21604.0629 3.524467 1982-1999 –> 2000-2012 7.922551 33525.904
Indo-Pacific (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 21633.3223 3.537285 1982-1999 –> 2000-2012 7.981613 33598.112
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 10533.7769 2.916430 2000-2012 –> 2013-2019 6.474949 32215.338
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 10518.2244 2.905717 2000-2012 –> 2013-2019 6.445569 32157.393
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 10519.2237 2.906404 2000-2012 –> 2013-2019 6.439058 32141.987
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 10518.4288 2.905857 2000-2012 –> 2013-2019 6.444841 32155.621
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate_star 10517.4085 2.906531 2000-2012 –> 2013-2019 6.457566 32179.999
Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 10583.7994 2.951157 2000-2012 –> 2013-2019 6.488216 32216.607
Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 10580.5180 2.948866 2000-2012 –> 2013-2019 6.473333 32184.581
Indo-Pacific (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 10584.4574 2.951616 2000-2012 –> 2013-2019 6.488901 32217.780
Indo-Pacific (27.75,27.85] cstar_tref ~ temp + nitrate + silicate + phosphate_star 10586.4683 2.954419 2000-2012 –> 2013-2019 6.491907 32218.252
Indo-Pacific (27.75,27.85] cstar_tref ~ temp + silicate + phosphate + phosphate_star 10595.9865 2.961081 2000-2012 –> 2013-2019 6.507480 32248.049
Indo-Pacific (27.85,27.95] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 29452.0123 3.204080 1982-1999 –> 2000-2012 7.673553 46413.019
Indo-Pacific (27.85,27.95] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 29302.4021 3.161736 1982-1999 –> 2000-2012 7.631093 46265.258
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 29406.8174 3.190839 1982-1999 –> 2000-2012 7.653800 46361.351
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 29322.6705 3.167365 1982-1999 –> 2000-2012 7.668471 46326.668
Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 29451.2896 3.203315 1982-1999 –> 2000-2012 7.672621 46414.079
Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 29296.5028 3.160100 1982-1999 –> 2000-2012 7.641585 46275.110
Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 29324.3651 3.167836 1982-1999 –> 2000-2012 7.663836 46321.766
Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 29327.5823 3.168730 1982-1999 –> 2000-2012 7.679647 46344.233
Indo-Pacific (27.85,27.95] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 29324.3806 3.167840 1982-1999 –> 2000-2012 7.665098 46323.407
Indo-Pacific (27.85,27.95] cstar_tref ~ temp + nitrate + silicate + phosphate_star 29322.4002 3.167845 1982-1999 –> 2000-2012 7.671457 46327.631
Indo-Pacific (27.85,27.95] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 15182.8373 2.954543 2000-2012 –> 2013-2019 6.116280 44485.239
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 15138.9160 2.933214 2000-2012 –> 2013-2019 6.101659 44465.472
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 15120.0523 2.924101 2000-2012 –> 2013-2019 6.098103 44466.579
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 15115.4262 2.921870 2000-2012 –> 2013-2019 6.089235 44438.097
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 15120.1034 2.924125 2000-2012 –> 2013-2019 6.098713 44468.732
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate_star 15118.1410 2.924143 2000-2012 –> 2013-2019 6.110157 44505.714
Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 15151.8042 2.939457 2000-2012 –> 2013-2019 6.099557 44448.307
Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 15143.1911 2.935283 2000-2012 –> 2013-2019 6.103119 44467.556
Indo-Pacific (27.85,27.95] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 15143.6006 2.935482 2000-2012 –> 2013-2019 6.103321 44467.981
Indo-Pacific (27.85,27.95] cstar_tref ~ temp + nitrate + silicate + phosphate_star 15154.4805 2.941725 2000-2012 –> 2013-2019 6.109571 44476.881
Indo-Pacific (27.95,28.05] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 48327.7274 3.003766 1982-1999 –> 2000-2012 7.496978 81194.988
Indo-Pacific (27.95,28.05] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 48221.5330 2.986872 1982-1999 –> 2000-2012 7.480026 81090.648
Indo-Pacific (27.95,28.05] cstar_tref ~ aou + silicate + phosphate + phosphate_star 48283.6428 2.996871 1982-1999 –> 2000-2012 7.506547 81192.036
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 48328.1549 3.003520 1982-1999 –> 2000-2012 7.469467 81128.965
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + aou + phosphate + phosphate_star 48353.5284 3.007809 1982-1999 –> 2000-2012 7.478186 81163.487
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 48203.9016 2.984128 1982-1999 –> 2000-2012 7.445229 80992.504
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 48310.6311 3.000778 1982-1999 –> 2000-2012 7.457818 81088.991
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 48353.6512 3.007515 1982-1999 –> 2000-2012 7.497419 81214.628
Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 48296.3433 2.998543 1982-1999 –> 2000-2012 7.488602 81157.708
Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 48280.6227 2.996087 1982-1999 –> 2000-2012 7.499129 81174.460
Indo-Pacific (27.95,28.05] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 27163.4148 3.072148 2000-2012 –> 2013-2019 6.075914 75491.142
Indo-Pacific (27.95,28.05] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 27164.1919 3.071796 2000-2012 –> 2013-2019 6.058668 75385.725
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 27150.9049 3.067979 2000-2012 –> 2013-2019 6.071499 75479.060
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 27167.2637 3.072680 2000-2012 –> 2013-2019 6.056807 75371.165
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 27123.3648 3.060080 2000-2012 –> 2013-2019 6.067595 75477.016
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 27120.9257 3.059382 2000-2012 –> 2013-2019 6.091598 75631.478
Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 27144.9173 3.066260 2000-2012 –> 2013-2019 6.064803 75441.261
Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 27117.1529 3.058302 2000-2012 –> 2013-2019 6.066197 75473.231
Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 27266.8544 3.101455 2000-2012 –> 2013-2019 6.097542 75547.477
Indo-Pacific (27.95,28.05] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 27126.9300 3.061102 2000-2012 –> 2013-2019 6.098326 75669.137
Indo-Pacific (28.05,28.1] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 38615.1438 2.933200 1982-1999 –> 2000-2012 7.075500 60593.172
Indo-Pacific (28.05,28.1] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 38608.9231 2.931642 1982-1999 –> 2000-2012 7.068374 60578.548
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 38555.1426 2.921469 1982-1999 –> 2000-2012 7.037439 60485.852
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + aou + phosphate + phosphate_star 38588.4922 2.928152 1982-1999 –> 2000-2012 7.044516 60517.943
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 38502.1850 2.911487 1982-1999 –> 2000-2012 7.027479 60432.936
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 38580.2006 2.926205 1982-1999 –> 2000-2012 7.040588 60507.928
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 38556.7160 2.921766 1982-1999 –> 2000-2012 7.052681 60515.457
Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 38614.2200 2.932646 1982-1999 –> 2000-2012 7.070671 60586.262
Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + phosphate + phosphate_star 38670.3420 2.943683 1982-1999 –> 2000-2012 7.089657 60655.226
Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 38663.3820 2.941979 1982-1999 –> 2000-2012 7.087094 60648.664
Indo-Pacific (28.05,28.1] cstar_tref ~ aou + nitrate + phosphate + phosphate_star 26829.4207 2.855005 2000-2012 –> 2013-2019 5.788205 65444.565
Indo-Pacific (28.05,28.1] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 26739.1554 2.830865 2000-2012 –> 2013-2019 5.762507 65348.078
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 26653.1965 2.808559 2000-2012 –> 2013-2019 5.730029 65208.339
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 26669.1551 2.812687 2000-2012 –> 2013-2019 5.740879 65259.856
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + aou + phosphate + phosphate_star 26785.8109 2.843570 2000-2012 –> 2013-2019 5.771722 65374.303
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + aou + silicate + phosphate 26807.0982 2.849146 2000-2012 –> 2013-2019 5.796368 65496.027
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 26778.9309 2.841247 2000-2012 –> 2013-2019 5.752734 65281.116
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 26752.6285 2.834378 2000-2012 –> 2013-2019 5.760582 65332.829
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 26809.0925 2.849145 2000-2012 –> 2013-2019 5.770911 65365.809
Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 26799.0464 2.846512 2000-2012 –> 2013-2019 5.779158 65413.266
Indo-Pacific (28.1, Inf] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 49698.3618 3.070012 1982-1999 –> 2000-2012 7.055519 91512.282
Indo-Pacific (28.1, Inf] cstar_tref ~ aou + silicate + phosphate + phosphate_star 50435.8817 3.188329 1982-1999 –> 2000-2012 7.186748 92296.060
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 51237.2780 3.321359 1982-1999 –> 2000-2012 7.357900 93241.032
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 50435.8195 3.187992 1982-1999 –> 2000-2012 7.184049 92289.180
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate 51439.3931 3.355864 1982-1999 –> 2000-2012 7.313721 93149.439
Indo-Pacific (28.1, Inf] cstar_tref ~ temp + aou + nitrate + phosphate 51717.5497 3.404286 1982-1999 –> 2000-2012 7.367167 93444.525
Indo-Pacific (28.1, Inf] cstar_tref ~ temp + aou + nitrate + phosphate + phosphate_star 51495.6925 3.365539 1982-1999 –> 2000-2012 7.317471 93183.388
Indo-Pacific (28.1, Inf] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 50962.3799 3.274998 1982-1999 –> 2000-2012 7.223446 92636.916
Indo-Pacific (28.1, Inf] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 50072.9439 3.129382 1982-1999 –> 2000-2012 7.069410 91715.627
Indo-Pacific (28.1, Inf] cstar_tref ~ temp + nitrate + silicate + phosphate + phosphate_star 51466.8050 3.360571 1982-1999 –> 2000-2012 7.392587 93453.822
Indo-Pacific (28.1, Inf] cstar_tref ~ aou + nitrate + silicate + phosphate + phosphate_star 22772.1003 2.718492 2000-2012 –> 2013-2019 5.788503 72470.462
Indo-Pacific (28.1, Inf] cstar_tref ~ aou + silicate + phosphate + phosphate_star 22821.2368 2.733308 2000-2012 –> 2013-2019 5.921637 73257.119
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + aou + nitrate + phosphate 22819.4088 2.732777 2000-2012 –> 2013-2019 6.113601 74401.717
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + aou + nitrate + phosphate + phosphate_star 22703.7931 2.698825 2000-2012 –> 2013-2019 6.020185 73941.071
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + aou + nitrate + silicate + phosphate 22818.8939 2.732047 2000-2012 –> 2013-2019 6.075269 74184.479
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 22584.2623 2.664753 2000-2012 –> 2013-2019 5.852746 73020.082
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate 22817.6347 2.731681 2000-2012 –> 2013-2019 6.087545 74257.028
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 22671.3415 2.689532 2000-2012 –> 2013-2019 6.135522 74629.003
Indo-Pacific (28.1, Inf] cstar_tref ~ temp + aou + nitrate + silicate + phosphate 23240.2784 2.857197 2000-2012 –> 2013-2019 6.132195 74202.658
Indo-Pacific (28.1, Inf] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 22669.4717 2.688998 2000-2012 –> 2013-2019 5.818380 72742.416

6.2 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.3 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.4 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
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
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
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
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-2012
(-Inf,26] 7 10 8 8 4 4 8
(26,26.5] 10 9 8 7 4 3 7
(26.5,26.75] 7 3 8 9 7 4 10
(26.75,27] 6 5 6 9 5 8 10
(27,27.25] 6 4 6 9 4 9 10
(27.25,27.5] 7 4 8 7 10 3 9
(27.5,27.75] 8 3 9 7 4 10 7
(27.75,27.85] 7 6 9 6 6 10 5
(27.85,27.95] 5 3 7 7 10 10 5
(27.95,28.05] 7 3 7 8 10 9 4
(28.05,28.1] 7 4 6 10 7 4 10
(28.1,28.15] 5 4 5 10 6 8 10
(28.15,28.2] 8 4 9 6 5 7 10
(28.2, Inf] 9 3 9 6 10 5 6
total 99.00 65.00 105.00 109.00 92.00 94.00 111.00
Atlantic - 2000-2012 --> 2013-2019
(-Inf,26] 8 10 7 7 5 5 8
(26,26.5] 9 10 8 7 4 3 7
(26.5,26.75] 7 10 7 7 3 3 10
(26.75,27] 5 4 6 10 5 8 10
(27,27.25] 6 6 6 8 6 8 10
(27.25,27.5] 8 5 9 7 9 3 8
(27.5,27.75] 8 3 9 7 4 10 7
(27.75,27.85] 7 7 7 8 6 10 5
(27.85,27.95] 5 3 7 7 10 10 5
(27.95,28.05] 6 3 7 8 10 9 5
(28.05,28.1] 7 6 8 10 7 5 6
(28.1,28.15] 5 5 9 10 5 3 10
(28.15,28.2] 7 3 7 7 10 3 10
(28.2, Inf] 5 3 7 10 10 5 7
total 93.00 78.00 104.00 113.00 94.00 85.00 108.00
Indo-Pacific - 1982-1999 --> 2000-2012
(-Inf,26] 7 6 6 7 10 6 8
(26,26.5] 8 4 8 8 10 4 7
(26.5,26.75] 9 7 9 9 7 2 5
(26.75,27] 10 5 10 9 3 5 5
(27,27.25] 10 4 10 9 4 5 5
(27.25,27.5] 10 4 10 9 4 5 5
(27.5,27.75] 10 5 10 9 5 3 5
(27.75,27.85] 7 5 7 8 5 10 8
(27.85,27.95] 7 9 7 9 2 7 7
(27.95,28.05] 10 4 10 9 5 5 4
(28.05,28.1] 10 4 10 9 5 4 5
(28.1, Inf] 9 7 10 7 3 6 6
total 107.00 64.00 107.00 102.00 63.00 62.00 70.00
Indo-Pacific - 2000-2012 --> 2013-2019
(-Inf,26] 8 8 7 6 8 4 8
(26,26.5] 6 5 6 8 10 6 8
(26.5,26.75] 7 7 7 9 9 3 6
(26.75,27] 9 7 9 8 4 4 7
(27,27.25] 9 3 9 7 8 7 6
(27.25,27.5] 8 4 8 8 8 6 8
(27.5,27.75] 5 4 6 8 10 5 10
(27.75,27.85] 4 4 5 9 5 10 10
(27.85,27.95] 5 6 5 8 5 10 9
(27.95,28.05] 8 7 9 8 4 7 6
(28.05,28.1] 10 5 10 7 7 5 3
(28.1, Inf] 10 6 10 6 6 6 4
total 89.00 66.00 91.00 92.00 84.00 73.00 85.00

7.4 RMSE alternatives

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
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
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.5 Residual patterns

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

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