Last updated: 2021-01-31

Checks: 7 0

Knit directory: emlr_mod_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: 9781315

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 9781315. 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 a618965 Donghe-Zhu 2021-01-31 Build site.
html 59e006e Donghe-Zhu 2021-01-31 Build site.
html a1c8f87 Donghe-Zhu 2021-01-31 Build site.
html ae5c18f Donghe-Zhu 2021-01-31 Build site.
html b50fe52 Donghe-Zhu 2021-01-31 Build site.
Rmd ac99ae5 jens-daniel-mueller 2021-01-29 code review
html ac99ae5 jens-daniel-mueller 2021-01-29 code review
html b5bdcaf Donghe-Zhu 2021-01-29 Build site.
Rmd b234505 Donghe-Zhu 2021-01-29 MLR approach across all depth
html 442010d Donghe-Zhu 2021-01-29 Build site.
Rmd e67e7dd Donghe-Zhu 2021-01-29 surface equilibrium approach across all latitudes irrespective of gamma
html 372adf5 Donghe-Zhu 2021-01-29 Build site.
html af8788e Donghe-Zhu 2021-01-29 Build site.
html 21c91c9 Donghe-Zhu 2021-01-29 Build site.
html eded038 Donghe-Zhu 2021-01-29 Build site.
html 541d4dd Donghe-Zhu 2021-01-29 Build site.
html 6a75576 Donghe-Zhu 2021-01-28 Build site.
html 16fba40 Donghe-Zhu 2021-01-28 Build site.
Rmd aecf0c6 Donghe-Zhu 2021-01-28 diagnostic ploting
Rmd a02684e Donghe-Zhu 2021-01-28 error
Rmd 4584be9 Donghe-Zhu 2021-01-28 latest
html 12bc567 Donghe-Zhu 2021-01-27 Build site.
html ceed31b Donghe-Zhu 2021-01-27 Build site.
html 342402d Donghe-Zhu 2021-01-27 Build site.
html 5bad5c2 Donghe-Zhu 2021-01-27 Build site.
Rmd c2c9529 Donghe-Zhu 2021-01-27 random subsetting based on lat
html 61efb56 Donghe-Zhu 2021-01-25 Build site.
html 48f638e Donghe-Zhu 2021-01-25 Build site.
html c1cec47 Donghe-Zhu 2021-01-25 Build site.
html 05ffb0c Donghe-Zhu 2021-01-25 Build site.
html 8b97165 Donghe-Zhu 2021-01-25 Build site.
html c569946 Donghe-Zhu 2021-01-24 Build site.
html a2f0d56 Donghe-Zhu 2021-01-23 Build site.
html 28509fc Donghe-Zhu 2021-01-23 Build site.
html 4c28e4a Donghe-Zhu 2021-01-22 Build site.
html 24cc264 jens-daniel-mueller 2021-01-22 cleaned /docs before creating copies
html 88eb28f Donghe-Zhu 2021-01-21 Build site.
html 2679490 Donghe-Zhu 2021-01-21 Build site.
html 7891955 Donghe-Zhu 2021-01-21 Build site.
html d4cf1cb Donghe-Zhu 2021-01-21 Build site.
Rmd 167eeec Donghe-Zhu 2021-01-21 surface DIC calculation with atmospheric equilibrium option
html 1f3e5b6 jens-daniel-mueller 2021-01-20 Build site.
html 0e7bdf1 jens-daniel-mueller 2021-01-15 cleaning template repository
html 73cbef3 jens-daniel-mueller 2021-01-15 Build site.
html 4571843 jens-daniel-mueller 2021-01-14 revision and html deleted for template copying
html 23151cd jens-daniel-mueller 2021-01-14 Build site.
html b3564aa jens-daniel-mueller 2021-01-14 Build site.
html 8d032c3 jens-daniel-mueller 2021-01-14 Build site.
html 022871c Donghe-Zhu 2021-01-13 Build site.
Rmd d44f36f Donghe-Zhu 2021-01-13 reorder analysis final
html 17dee1d jens-daniel-mueller 2021-01-13 Build site.
Rmd 9e04fd7 jens-daniel-mueller 2021-01-13 local rebuild after revision
html a076226 Donghe-Zhu 2021-01-11 Build site.
Rmd 52eff18 Donghe-Zhu 2021-01-09 Implemet model_run and subsetting
html 7cdea0c jens-daniel-mueller 2021-01-06 Build site.
Rmd b5934dd jens-daniel-mueller 2021-01-06 local rebuild after revision
html fa85b93 jens-daniel-mueller 2021-01-06 Build site.
html e5cb81a Donghe-Zhu 2021-01-05 Build site.
Rmd 608cc45 Donghe-Zhu 2021-01-05 modification of analysis
html a499f10 Donghe-Zhu 2021-01-05 Build site.
Rmd 715bdb4 Donghe-Zhu 2021-01-02 model modification
html fb8a752 Donghe-Zhu 2020-12-23 Build site.
Rmd 82e3c9c Donghe-Zhu 2020-12-23 first build after creating model template
html 8fae0b2 Donghe-Zhu 2020-12-21 Build site.
Rmd 00a1322 Donghe-Zhu 2020-12-21 first build after creating model template
Rmd d73ae35 Donghe-Zhu 2020-12-21 first version with lm error
html c8b76b3 jens-daniel-mueller 2020-12-19 Build site.
Rmd b5fedce jens-daniel-mueller 2020-12-19 first build after creating model template
Rmd 8e8abf5 Jens Müller 2020-12-18 Initial commit

1 Required data

Required are:

  • cleaned and prepared GLODAP-based synthetic model subsetting 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)

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

# remove predictor combinations with nitrate and phosphate
lm_all <- lm_all %>%
  mutate(lm_coeff_filter = str_remove(lm_coeff, "phosphate_star")) %>%
  filter(!(
    str_detect(lm_coeff_filter, "nitrate") &
      str_detect(lm_coeff_filter, "phosphate")
  )) %>%
  select(-lm_coeff_filter)

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

  • 86

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.

# loop across all basins, era, gamma slabs, and MLRs
# fit all MLR models
for (i_basin in unique(GLODAP$basin)) {
  for (i_era in unique(GLODAP$era)) {
    # i_basin <- unique(GLODAP$basin)[1]
    # i_era   <- unique(GLODAP$era)[1]
    print(i_basin)
    print(i_era)
    
    GLODAP_basin_era <- GLODAP %>%
      filter(basin == i_basin, era == i_era)
    
    for (i_gamma_slab in unique(GLODAP_basin_era$gamma_slab)) {
      # i_gamma_slab <- unique(GLODAP_basin_era$gamma_slab)[1]
      print(i_gamma_slab)
      
      GLODAP_basin_era_slab <- GLODAP_basin_era %>%
        filter(gamma_slab == i_gamma_slab)
      
      for (i_predictors in unique(lm_all$predictors)) {
        # i_predictors <- unique(lm_all$predictors)[110]
        
        # extract one model definition
        i_lm <- lm_all %>%
          filter(predictors == i_predictors) %>%
          select(lm_coeff) %>%
          pull()
        
        # extract number of predictors
        i_n_predictors <- lm_all %>%
          filter(predictors == i_predictors) %>%
          select(n) %>%
          pull()
        
        # fit model
        i_lm_fit <- lm(as.formula(i_lm),
                       data = GLODAP_basin_era_slab)
        
        # find max predictor correlation
        i_cor_max <- GLODAP_basin_era_slab %>%
          select(!!!syms(str_split(i_predictors, " ",
                                   simplify = TRUE))) %>%
          correlate(quiet = TRUE) %>%
          select(-term) %>%
          abs() %>%
          max(na.rm = TRUE)
        
        # calculate root mean squared error
        i_rmse <- sqrt(c(crossprod(i_lm_fit$residuals)) /
                         length(i_lm_fit$residuals))
        
        # calculate maximum residual
        i_resid_max <- max(abs(i_lm_fit$residuals))
        
        # calculate Akaike information criterion aic
        i_aic <- AIC(i_lm_fit)
        
        # collect model coefficients and diagnostics
        coefficients <- tidy(i_lm_fit)
        
        coefficients <- coefficients %>%
          mutate(
            basin = i_basin,
            era = i_era,
            gamma_slab = i_gamma_slab,
            model = i_lm,
            rmse = i_rmse,
            aic = i_aic,
            resid_max = i_resid_max,
            n_predictors = i_n_predictors,
            na_predictor = anyNA(coefficients$estimate),
            cor_max = i_cor_max
          )
        
        if (exists("lm_all_fitted")) {
          lm_all_fitted <- bind_rows(lm_all_fitted, coefficients)
        }
        
        if (!exists("lm_all_fitted")) {
          lm_all_fitted <- coefficients
        }
      }
      
      
    }
  }
  
}

rm(
  i_lm_fit,
  coefficients,
  i_rmse,
  GLODAP_basin_era,
  GLODAP_basin_era_slab,
  i_lm,
  i_basin,
  i_era,
  i_gamma_slab,
  i_predictors,
  i_aic,
  i_n_predictors,
  i_resid_max
)

5 Prepare coeffcients

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

5.1 Formatting

# select relevant columns
lm_all_fitted <- lm_all_fitted %>%
  select(
    basin,
    gamma_slab,
    era,
    model,
    n_predictors,
    term,
    estimate,
    rmse,
    aic,
    resid_max,
    na_predictor,
    cor_max
  )

# set coefficient to zero if not fitted (=NA)
lm_all_fitted <- lm_all_fitted %>%
  mutate(estimate = if_else(is.na(estimate), 0, estimate))

# Prepare model coefficients for mapping of target variable
lm_all_fitted_wide <- lm_all_fitted %>%
  pivot_wider(
    values_from = estimate,
    names_from = term,
    names_prefix = "coeff_",
    values_fill = 0
  )

5.2 Predictor selection

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

  • 10

The criterion used to select the best models was:

  • rmse

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

Please note, that currently the lm() function produces NAs for some predictors. It is not yet entirely clear when this happens, but presumably it is caused by some form of collinearity between predictors, such that including another predictor does not help to explain the target variable any better. The issues also expresses as exactly identical rmse values of different models. As an interim solution, models with fitted NA predictors were not included.

# remove models with predictors fitted as NA

lm_all_fitted_wide <- lm_all_fitted_wide %>%
  filter(na_predictor == FALSE)
# calculate RMSE sum for adjacent eras
lm_all_fitted_wide_eras <- lm_all_fitted_wide  %>%
  select(basin, gamma_slab, model, era, rmse, aic, resid_max) %>%
  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 <- lm_all_fitted_wide_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 <- lm_all_fitted_wide_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 %>%
  kable() %>%
  add_header_above() %>%
  kable_styling() %>%
  scroll_box(width = "100%", height = "400px")
basin gamma_slab model rmse aic resid_max eras rmse_sum aic_sum
Atlantic (-Inf,26] cstar_tref ~ sal + aou + phosphate + phosphate_star 4.5828295 1782.63560 47.2761834 1982-1999 –> 2000-2012 8.273915 3140.716
Atlantic (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 3.9145177 1689.74566 34.3406149 1982-1999 –> 2000-2012 7.266315 3002.193
Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 4.2975616 1745.94577 39.0790802 1982-1999 –> 2000-2012 7.689503 3064.274
Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 3.0746466 1544.36136 11.9637284 1982-1999 –> 2000-2012 6.683866 2893.362
Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 4.1907786 1730.79869 39.5959733 1982-1999 –> 2000-2012 7.496552 3036.415
Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4.2899390 1744.87705 40.9596911 1982-1999 –> 2000-2012 7.586689 3049.144
Atlantic (-Inf,26] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 4.3077335 1747.36895 41.0322273 1982-1999 –> 2000-2012 7.584080 3048.569
Atlantic (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4.3052878 1747.02707 40.7779269 1982-1999 –> 2000-2012 7.600182 3051.015
Atlantic (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate_star 4.4128160 1759.87784 34.5411394 1982-1999 –> 2000-2012 7.920060 3092.720
Atlantic (-Inf,26] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 4.9989890 1836.96089 28.8835408 1982-1999 –> 2000-2012 8.353833 3149.857
Atlantic (-Inf,26] cstar_tref ~ sal + aou + phosphate + phosphate_star 0.0000000 -Inf 0.0000000 2000-2012 –> 2013-2019 4.582829 -Inf
Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate 0.0000000 -Inf 0.0000000 2000-2012 –> 2013-2019 4.654563 -Inf
Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + phosphate 0.0000000 -Inf 0.0000000 2000-2012 –> 2013-2019 4.920736 -Inf
Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + phosphate_star 0.0000000 -Inf 0.0000000 2000-2012 –> 2013-2019 4.972483 -Inf
Atlantic (-Inf,26] cstar_tref ~ sal + temp + aou + silicate 0.0000000 -Inf 0.0000000 2000-2012 –> 2013-2019 4.717846 -Inf
Atlantic (-Inf,26] cstar_tref ~ sal + temp + nitrate + silicate 0.0000000 -Inf 0.0000000 2000-2012 –> 2013-2019 4.745230 -Inf
Atlantic (-Inf,26] cstar_tref ~ sal + temp + phosphate + phosphate_star 0.0000000 -Inf 0.0000000 2000-2012 –> 2013-2019 4.995094 -Inf
Atlantic (-Inf,26] cstar_tref ~ sal + temp + silicate 0.0290845 -11.18612 0.0467729 2000-2012 –> 2013-2019 4.816118 1795.693
Atlantic (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate 0.0000000 -Inf 0.0000000 2000-2012 –> 2013-2019 4.784756 -Inf
Atlantic (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate_star 0.0000000 -Inf 0.0000000 2000-2012 –> 2013-2019 4.412816 -Inf
Atlantic (26,26.5] cstar_tref ~ aou + silicate + phosphate + phosphate_star 5.2331502 4241.75789 16.5200105 1982-1999 –> 2000-2012 9.269188 7298.724
Atlantic (26,26.5] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 4.8510735 4139.43865 18.1541948 1982-1999 –> 2000-2012 8.571802 7110.391
Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 5.0903881 4205.69862 18.0523733 1982-1999 –> 2000-2012 9.086848 7254.003
Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 5.0441888 4193.15334 23.3381326 1982-1999 –> 2000-2012 9.106402 7259.114
Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5.0968271 4207.43807 24.6248004 1982-1999 –> 2000-2012 9.281845 7305.625
Atlantic (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5.1087642 4210.65699 24.8169049 1982-1999 –> 2000-2012 9.326619 7317.300
Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + phosphate_star 5.3305337 4267.12851 13.9285460 1982-1999 –> 2000-2012 9.486430 7355.759
Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5.2612574 4251.12862 18.1981171 1982-1999 –> 2000-2012 9.380623 7332.207
Atlantic (26,26.5] cstar_tref ~ temp + aou + silicate + phosphate 5.2596682 4248.71292 19.3666605 1982-1999 –> 2000-2012 9.399526 7333.160
Atlantic (26,26.5] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 5.2283982 4242.50784 15.5787758 1982-1999 –> 2000-2012 9.231396 7292.580
Atlantic (26,26.5] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2.1609414 1008.02506 9.2801010 2000-2012 –> 2013-2019 7.012015 5147.464
Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate 1.6684215 888.59068 10.1335049 2000-2012 –> 2013-2019 7.225472 5212.983
Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 1.6384153 882.35128 9.7158579 2000-2012 –> 2013-2019 6.728803 5088.050
Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + nitrate + silicate 1.6320250 880.57708 9.8067953 2000-2012 –> 2013-2019 6.974952 5152.901
Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 2.1721100 1010.36546 8.9767754 2000-2012 –> 2013-2019 7.216299 5203.519
Atlantic (26,26.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2.1992179 1015.99631 8.7240489 2000-2012 –> 2013-2019 7.296045 5223.434
Atlantic (26,26.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 2.2037978 1016.94078 8.6874260 2000-2012 –> 2013-2019 7.312562 5227.598
Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + phosphate_star 1.6413580 881.16597 9.6469308 2000-2012 –> 2013-2019 6.971892 5148.294
Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate 1.6323682 878.67255 9.8549118 2000-2012 –> 2013-2019 7.307070 5231.893
Atlantic (26,26.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1.6279676 879.44699 9.6673167 2000-2012 –> 2013-2019 6.889225 5130.576
Atlantic (26.5,27.25] cstar_tref ~ aou + silicate + phosphate + phosphate_star 4.1304796 18392.23713 16.4822877 1982-1999 –> 2000-2012 6.920251 30159.308
Atlantic (26.5,27.25] cstar_tref ~ sal + aou + silicate + phosphate 3.9856786 18161.06352 19.3336082 1982-1999 –> 2000-2012 7.610436 31187.001
Atlantic (26.5,27.25] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 3.7044507 17689.05142 16.3967046 1982-1999 –> 2000-2012 6.347888 29199.069
Atlantic (26.5,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate 3.8662359 17965.96266 17.3406698 1982-1999 –> 2000-2012 6.780649 29945.186
Atlantic (26.5,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4.0999901 18346.24173 18.9870846 1982-1999 –> 2000-2012 7.357874 30861.120
Atlantic (26.5,27.25] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4.1695086 18455.16048 20.5068020 1982-1999 –> 2000-2012 7.552913 31151.803
Atlantic (26.5,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 3.8982696 18019.41498 15.6575524 1982-1999 –> 2000-2012 6.951070 30221.684
Atlantic (26.5,27.25] cstar_tref ~ temp + aou + silicate + phosphate 4.4629117 18893.68535 17.3249451 1982-1999 –> 2000-2012 7.490025 31053.329
Atlantic (26.5,27.25] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 4.0903640 18331.01467 16.0482877 1982-1999 –> 2000-2012 6.870925 30084.185
Atlantic (26.5,27.25] cstar_tref ~ temp + aou + silicate + phosphate_star 4.6621925 19176.67299 17.5224437 1982-1999 –> 2000-2012 7.990535 31792.426
Atlantic (26.5,27.25] cstar_tref ~ aou + silicate + phosphate + phosphate_star 3.2698587 3355.13193 10.4947186 2000-2012 –> 2013-2019 7.400338 21747.369
Atlantic (26.5,27.25] cstar_tref ~ sal + aou + silicate + phosphate 3.0960299 3284.99194 11.9477735 2000-2012 –> 2013-2019 7.081709 21446.055
Atlantic (26.5,27.25] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2.9608157 3229.65389 10.2512761 2000-2012 –> 2013-2019 6.665266 20918.705
Atlantic (26.5,27.25] cstar_tref ~ sal + temp + aou + nitrate + silicate 3.1414755 3305.70238 12.2050116 2000-2012 –> 2013-2019 7.530317 22092.971
Atlantic (26.5,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate 3.0683805 3275.47354 11.1287415 2000-2012 –> 2013-2019 6.934616 21241.436
Atlantic (26.5,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 3.2658555 3355.55901 15.0533681 2000-2012 –> 2013-2019 7.365846 21701.801
Atlantic (26.5,27.25] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 3.3163044 3375.24181 16.3331194 2000-2012 –> 2013-2019 7.485813 21830.402
Atlantic (26.5,27.25] cstar_tref ~ temp + aou + nitrate + silicate 3.1560458 3309.64385 12.2578048 2000-2012 –> 2013-2019 7.638891 22232.199
Atlantic (26.5,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2.9565206 3227.78988 10.8428371 2000-2012 –> 2013-2019 6.854790 21247.205
Atlantic (26.5,27.25] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 3.1563265 3311.75803 11.7658020 2000-2012 –> 2013-2019 7.246690 21642.773
Atlantic (27.25,27.5] cstar_tref ~ sal + aou + silicate + phosphate 3.4367243 12764.51582 19.3549171 1982-1999 –> 2000-2012 5.136796 18262.716
Atlantic (27.25,27.5] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 3.4336790 12762.25534 19.7316052 1982-1999 –> 2000-2012 4.995095 18023.047
Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + silicate 3.3082064 12583.34669 12.2714795 1982-1999 –> 2000-2012 5.213635 18404.445
Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 3.4301570 12757.32314 18.9279349 1982-1999 –> 2000-2012 4.886845 17822.746
Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 3.3266989 12610.13688 20.2900432 1982-1999 –> 2000-2012 4.757197 17624.505
Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 3.3029495 12575.70366 19.6013031 1982-1999 –> 2000-2012 4.803986 17725.520
Atlantic (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 3.3274602 12611.23657 20.2641140 1982-1999 –> 2000-2012 4.767211 17643.749
Atlantic (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate_star 3.3477053 12638.38880 20.9516770 1982-1999 –> 2000-2012 5.572217 18893.180
Atlantic (27.25,27.5] cstar_tref ~ temp + aou + nitrate + silicate 3.5205583 12880.34426 10.8326642 1982-1999 –> 2000-2012 5.451764 18737.256
Atlantic (27.25,27.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 3.5072449 12864.13540 10.7297312 1982-1999 –> 2000-2012 5.218080 18382.094
Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate 3.0354968 2212.50286 10.5112938 2000-2012 –> 2013-2019 6.434650 14924.189
Atlantic (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 3.0228627 2210.87426 10.8919133 2000-2012 –> 2013-2019 6.366292 14845.121
Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + nitrate + silicate 2.7106315 2116.02454 10.6030217 2000-2012 –> 2013-2019 6.018838 14699.371
Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 3.2475216 2273.24267 11.3852910 2000-2012 –> 2013-2019 6.677679 15030.566
Atlantic (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 3.0762144 2226.09532 9.7338204 2000-2012 –> 2013-2019 6.402913 14836.232
Atlantic (27.25,27.5] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 3.0833217 2228.10303 9.4137639 2000-2012 –> 2013-2019 6.386271 14803.807
Atlantic (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 3.0730664 2225.20455 9.7809944 2000-2012 –> 2013-2019 6.400527 14836.441
Atlantic (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate_star 3.0848784 2226.54217 9.1879408 2000-2012 –> 2013-2019 6.432584 14864.931
Atlantic (27.25,27.5] cstar_tref ~ temp + aou + nitrate + silicate 2.7675418 2132.10129 10.8403267 2000-2012 –> 2013-2019 6.288100 15012.446
Atlantic (27.25,27.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2.6954446 2111.13648 11.1195188 2000-2012 –> 2013-2019 6.202689 14975.272
Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2.7987832 19878.07239 19.6550981 1982-1999 –> 2000-2012 4.525050 30828.703
Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 3.1726120 20895.32839 24.7977813 1982-1999 –> 2000-2012 4.544697 30567.839
Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + phosphate 3.1378414 20803.91113 18.4929507 1982-1999 –> 2000-2012 4.916489 31918.921
Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 3.0683222 20624.12349 15.4298356 1982-1999 –> 2000-2012 4.747908 31422.168
Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 2.8181068 19933.90127 20.7472313 1982-1999 –> 2000-2012 4.400938 30401.704
Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2.8954211 20153.50859 23.0449360 1982-1999 –> 2000-2012 4.350971 30154.706
Atlantic (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 2.9025072 20173.34207 22.6874895 1982-1999 –> 2000-2012 4.352538 30153.395
Atlantic (27.5,27.75] cstar_tref ~ sal + temp + phosphate + phosphate_star 3.2449829 21076.33839 20.2462241 1982-1999 –> 2000-2012 4.936900 31913.099
Atlantic (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 2.8873189 20130.77164 22.7690006 1982-1999 –> 2000-2012 4.341917 30128.328
Atlantic (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate_star 3.0194470 20491.83530 25.9368483 1982-1999 –> 2000-2012 4.482099 30518.125
Atlantic (27.5,27.75] cstar_tref ~ sal + aou + phosphate + phosphate_star 3.4938225 4342.63446 19.7011501 2000-2012 –> 2013-2019 6.590496 25039.386
Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate 3.1061778 4151.88171 20.3499569 2000-2012 –> 2013-2019 6.005333 24313.847
Atlantic (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 3.0994585 4150.36918 19.5951821 2000-2012 –> 2013-2019 5.898242 24028.442
Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + phosphate 3.5220210 4355.67302 20.5841897 2000-2012 –> 2013-2019 6.659862 25159.584
Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 3.4515248 4324.87803 18.5862785 2000-2012 –> 2013-2019 6.519847 24949.002
Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 3.1015372 4151.45661 19.8633878 2000-2012 –> 2013-2019 5.919644 24085.358
Atlantic (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 3.1373540 4170.08030 20.9981844 2000-2012 –> 2013-2019 6.032775 24323.589
Atlantic (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 3.1577835 4180.60801 21.4193134 2000-2012 –> 2013-2019 6.060291 24353.950
Atlantic (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 3.1326269 4167.63457 20.8525678 2000-2012 –> 2013-2019 6.019946 24298.406
Atlantic (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate_star 3.2215170 4211.01884 22.6113428 2000-2012 –> 2013-2019 6.240964 24702.854
Atlantic (27.75,27.85] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1.9881768 7170.72001 10.2642278 1982-1999 –> 2000-2012 2.869421 9684.450
Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 2.2197531 7545.10471 11.7512077 1982-1999 –> 2000-2012 3.165189 10194.817
Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 1.4600325 6121.55569 12.0498550 1982-1999 –> 2000-2012 2.264755 8459.607
Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 2.1855949 7492.40875 11.6234970 1982-1999 –> 2000-2012 3.127160 10134.187
Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2.3044964 7672.41483 12.0327175 1982-1999 –> 2000-2012 3.318332 10457.217
Atlantic (27.75,27.85] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 2.3154634 7688.54732 12.1574357 1982-1999 –> 2000-2012 3.329695 10474.105
Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 2.2949054 7658.24329 11.9571080 1982-1999 –> 2000-2012 3.308354 10442.306
Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate_star 2.3661490 7760.12728 12.6327354 1982-1999 –> 2000-2012 3.380697 10544.288
Atlantic (27.75,27.85] cstar_tref ~ temp + aou + phosphate + phosphate_star 1.8355053 6897.22582 14.0440158 1982-1999 –> 2000-2012 2.910585 9793.464
Atlantic (27.75,27.85] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1.7574653 6751.59228 12.7257741 1982-1999 –> 2000-2012 2.809668 9608.233
Atlantic (27.75,27.85] cstar_tref ~ aou + silicate + phosphate + phosphate_star 1.6718699 1121.47381 5.6482917 2000-2012 –> 2013-2019 3.778278 8486.484
Atlantic (27.75,27.85] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1.6107605 1102.10021 5.8523820 2000-2012 –> 2013-2019 3.598937 8272.820
Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 1.0763727 870.71527 4.0182245 2000-2012 –> 2013-2019 3.296126 8415.820
Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 1.0905803 878.24221 6.8582071 2000-2012 –> 2013-2019 2.550613 6999.798
Atlantic (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 1.7992538 1165.62225 5.7939819 2000-2012 –> 2013-2019 3.984849 8658.031
Atlantic (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 1.8656360 1186.41826 5.8537907 2000-2012 –> 2013-2019 4.160541 8844.662
Atlantic (27.75,27.85] cstar_tref ~ temp + aou + nitrate + phosphate_star 1.4705731 1047.83507 5.0367269 2000-2012 –> 2013-2019 3.879236 8868.474
Atlantic (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1.4653949 1047.81032 4.9242659 2000-2012 –> 2013-2019 3.819790 8793.017
Atlantic (27.75,27.85] cstar_tref ~ temp + aou + phosphate + phosphate_star 1.3401873 994.54329 8.3302122 2000-2012 –> 2013-2019 3.175693 7891.769
Atlantic (27.75,27.85] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1.2571481 959.82817 6.1086443 2000-2012 –> 2013-2019 3.014613 7711.420
Atlantic (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2.6193544 7826.52187 12.9954728 1982-1999 –> 2000-2012 3.682022 10888.746
Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 2.9577916 8225.09110 14.4898183 1982-1999 –> 2000-2012 3.972594 11192.375
Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 2.3899509 7525.89324 12.8731042 1982-1999 –> 2000-2012 3.332096 10340.139
Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 2.6827902 7905.01062 13.1713047 1982-1999 –> 2000-2012 3.786378 11045.072
Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2.7319420 7964.56014 13.1805752 1982-1999 –> 2000-2012 3.884267 11193.645
Atlantic (27.85,27.95] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 2.7450889 7980.30657 13.4070570 1982-1999 –> 2000-2012 3.898684 11211.659
Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 2.7211968 7951.63385 13.1104360 1982-1999 –> 2000-2012 3.872859 11179.532
Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate_star 2.8864126 8142.96564 15.1603385 1982-1999 –> 2000-2012 4.040818 11373.766
Atlantic (27.85,27.95] cstar_tref ~ temp + aou + phosphate + phosphate_star 2.5101775 7684.87776 12.4806275 1982-1999 –> 2000-2012 3.629024 10851.226
Atlantic (27.85,27.95] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 2.5044142 7679.33829 12.2635275 1982-1999 –> 2000-2012 3.620115 10841.886
Atlantic (27.85,27.95] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1.2102817 777.04127 5.5865242 2000-2012 –> 2013-2019 3.829636 8603.563
Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 0.9630503 668.73096 6.2708895 2000-2012 –> 2013-2019 3.920842 8893.822
Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 0.9439067 659.21382 7.1568716 2000-2012 –> 2013-2019 3.333858 8185.107
Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate 1.2678996 799.08631 5.6777445 2000-2012 –> 2013-2019 3.950690 8704.097
Atlantic (27.85,27.95] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1.3206689 818.41445 5.8719880 2000-2012 –> 2013-2019 4.052611 8782.975
Atlantic (27.85,27.95] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 1.3490776 828.50247 6.0846501 2000-2012 –> 2013-2019 4.094167 8808.809
Atlantic (27.85,27.95] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 1.3101869 814.63736 5.7999938 2000-2012 –> 2013-2019 4.031384 8766.271
Atlantic (27.85,27.95] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1.2848730 805.38963 8.1583309 2000-2012 –> 2013-2019 4.217350 9002.288
Atlantic (27.85,27.95] cstar_tref ~ temp + aou + phosphate + phosphate_star 1.2670816 796.78039 7.8560905 2000-2012 –> 2013-2019 3.777259 8481.658
Atlantic (27.85,27.95] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1.2668847 798.70674 7.8103392 2000-2012 –> 2013-2019 3.771299 8478.045
Atlantic (27.95,28.05] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 4.7725520 11523.82405 19.5727243 1982-1999 –> 2000-2012 8.484198 18544.066
Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 3.6165121 10453.17157 23.1399681 1982-1999 –> 2000-2012 5.831518 16148.792
Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 4.7821515 11531.58029 19.7543927 1982-1999 –> 2000-2012 8.632037 18645.656
Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 4.7763718 11526.91225 19.5233482 1982-1999 –> 2000-2012 8.485114 18545.146
Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 4.7806912 11530.40135 19.4911998 1982-1999 –> 2000-2012 8.490193 18549.161
Atlantic (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 4.8012609 11546.97399 19.4136588 1982-1999 –> 2000-2012 8.316739 18427.883
Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 4.7801524 11529.96634 19.4991041 1982-1999 –> 2000-2012 8.490141 18549.063
Atlantic (27.95,28.05] cstar_tref ~ temp + aou + nitrate + phosphate_star 3.7435727 10584.45879 23.0305784 1982-1999 –> 2000-2012 6.056346 16388.911
Atlantic (27.95,28.05] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 3.0985202 9856.47183 28.1316256 1982-1999 –> 2000-2012 4.739309 14782.093
Atlantic (27.95,28.05] cstar_tref ~ temp + nitrate + silicate + phosphate_star 5.0088477 11708.35745 20.0312896 1982-1999 –> 2000-2012 8.638294 18669.133
Atlantic (27.95,28.05] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 0.7052324 583.08605 2.6380008 2000-2012 –> 2013-2019 5.477784 12106.910
Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 0.6090338 505.06630 2.7410760 2000-2012 –> 2013-2019 4.225546 10958.238
Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 0.4282354 317.69557 1.7580876 2000-2012 –> 2013-2019 5.210387 11849.276
Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 0.7717527 631.03884 3.2511944 2000-2012 –> 2013-2019 5.548125 12157.951
Atlantic (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 0.8233219 665.45022 3.5360274 2000-2012 –> 2013-2019 5.604013 12195.852
Atlantic (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 0.8140349 659.41520 3.5104036 2000-2012 –> 2013-2019 5.594187 12189.382
Atlantic (27.95,28.05] cstar_tref ~ temp + aou + nitrate + phosphate_star 1.0018614 767.86467 3.4924624 2000-2012 –> 2013-2019 4.745434 11352.323
Atlantic (27.95,28.05] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 0.8248152 666.41425 4.0059800 2000-2012 –> 2013-2019 3.923335 10522.886
Atlantic (27.95,28.05] cstar_tref ~ temp + aou + phosphate + phosphate_star 0.7129087 586.84546 3.8376580 2000-2012 –> 2013-2019 5.649299 12238.957
Atlantic (27.95,28.05] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 0.7112986 587.64258 3.8532384 2000-2012 –> 2013-2019 5.646352 12240.708
Atlantic (28.05,28.1] cstar_tref ~ sal + aou + phosphate + phosphate_star 1.4547274 3947.50296 10.6472542 1982-1999 –> 2000-2012 2.169645 5530.362
Atlantic (28.05,28.1] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1.4335612 3917.34586 10.0452931 1982-1999 –> 2000-2012 2.147715 5500.653
Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate 1.4778810 3982.14796 10.3635143 1982-1999 –> 2000-2012 2.200587 5580.717
Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 1.3824916 3837.76007 11.3003457 1982-1999 –> 2000-2012 2.060432 5345.613
Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 1.4625528 3961.27353 9.8019103 1982-1999 –> 2000-2012 2.184382 5560.082
Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1.4832548 3992.11118 9.7392278 1982-1999 –> 2000-2012 2.217653 5615.952
Atlantic (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 1.4746403 3979.33171 9.6428618 1982-1999 –> 2000-2012 2.205474 5596.117
Atlantic (28.05,28.1] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1.4126694 3885.13662 8.6567917 1982-1999 –> 2000-2012 2.213167 5633.941
Atlantic (28.05,28.1] cstar_tref ~ temp + aou + phosphate + phosphate_star 1.3827446 3836.16148 11.3339728 1982-1999 –> 2000-2012 2.079381 5381.459
Atlantic (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1.2354526 3591.04476 9.9114589 1982-1999 –> 2000-2012 1.865486 4992.632
Atlantic (28.05,28.1] cstar_tref ~ sal + aou + phosphate + phosphate_star 0.4291476 225.15010 1.7846102 2000-2012 –> 2013-2019 1.883875 4172.653
Atlantic (28.05,28.1] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 0.4195171 218.70700 1.7112354 2000-2012 –> 2013-2019 1.853078 4136.053
Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 0.4499877 244.79012 1.7810074 2000-2012 –> 2013-2019 1.943204 4251.587
Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate 0.4659034 255.72009 1.9728635 2000-2012 –> 2013-2019 1.943784 4237.868
Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 0.3366016 136.79097 1.5570948 2000-2012 –> 2013-2019 1.719093 3974.551
Atlantic (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 0.4484156 243.48821 1.8624228 2000-2012 –> 2013-2019 1.910968 4204.762
Atlantic (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 0.4681852 259.53759 1.9383973 2000-2012 –> 2013-2019 1.942826 4238.869
Atlantic (28.05,28.1] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 0.4733482 263.61742 2.3906660 2000-2012 –> 2013-2019 1.886018 4148.754
Atlantic (28.05,28.1] cstar_tref ~ temp + aou + phosphate + phosphate_star 0.4580099 249.36351 2.6788983 2000-2012 –> 2013-2019 1.840754 4085.525
Atlantic (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 0.3431583 143.96758 2.4199022 2000-2012 –> 2013-2019 1.578611 3735.012
Atlantic (28.1,28.15] cstar_tref ~ aou + silicate + phosphate + phosphate_star 1.3132228 5302.77062 11.7129137 1982-1999 –> 2000-2012 1.997526 7075.828
Atlantic (28.1,28.15] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1.2955374 5262.35892 10.2036101 1982-1999 –> 2000-2012 1.909983 6855.003
Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 1.3197918 5320.37831 12.5607717 1982-1999 –> 2000-2012 1.983473 7043.599
Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + silicate + phosphate 1.3083688 5293.18720 10.0207472 1982-1999 –> 2000-2012 1.932937 6913.511
Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1.3226012 5327.02985 9.8715304 1982-1999 –> 2000-2012 1.961659 6986.207
Atlantic (28.1,28.15] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 1.3126104 5303.31145 10.0247145 1982-1999 –> 2000-2012 1.946252 6948.068
Atlantic (28.1,28.15] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1.2032252 5031.13734 9.5167410 1982-1999 –> 2000-2012 1.874728 6774.206
Atlantic (28.1,28.15] cstar_tref ~ temp + aou + phosphate + phosphate_star 1.3262336 5333.60889 11.2785823 1982-1999 –> 2000-2012 1.995514 7069.062
Atlantic (28.1,28.15] cstar_tref ~ temp + aou + silicate + phosphate 1.3344027 5352.81701 11.8338778 1982-1999 –> 2000-2012 2.052718 7208.045
Atlantic (28.1,28.15] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1.1245324 4819.56476 13.6240895 1982-1999 –> 2000-2012 1.629827 6080.899
Atlantic (28.1,28.15] cstar_tref ~ aou + silicate + phosphate + phosphate_star 0.4895828 408.06216 1.7039831 2000-2012 –> 2013-2019 1.802806 5710.833
Atlantic (28.1,28.15] cstar_tref ~ sal + aou + phosphate + phosphate_star 0.4037117 299.67895 1.7095319 2000-2012 –> 2013-2019 1.781627 5752.867
Atlantic (28.1,28.15] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 0.3856113 275.89937 1.4310718 2000-2012 –> 2013-2019 1.681149 5538.258
Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 0.3922726 285.52476 1.6867747 2000-2012 –> 2013-2019 1.712064 5605.903
Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + silicate + phosphate 0.4043632 302.58522 1.4758757 2000-2012 –> 2013-2019 1.712732 5595.772
Atlantic (28.1,28.15] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 0.4200192 323.93378 1.5314873 2000-2012 –> 2013-2019 1.742620 5650.964
Atlantic (28.1,28.15] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 0.4178209 320.98475 1.5136772 2000-2012 –> 2013-2019 1.730431 5624.296
Atlantic (28.1,28.15] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 0.5341289 459.00314 1.6821036 2000-2012 –> 2013-2019 1.737354 5490.140
Atlantic (28.1,28.15] cstar_tref ~ temp + aou + phosphate + phosphate_star 0.4126692 312.01225 1.6827518 2000-2012 –> 2013-2019 1.738903 5645.621
Atlantic (28.1,28.15] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 0.2771681 90.32370 1.1114077 2000-2012 –> 2013-2019 1.401701 4909.888
Atlantic (28.15,28.2] cstar_tref ~ sal + aou + nitrate + phosphate_star 1.0048762 8438.06575 11.5293844 1982-1999 –> 2000-2012 1.652106 11389.914
Atlantic (28.15,28.2] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 0.9937731 8374.31211 12.3479982 1982-1999 –> 2000-2012 1.640860 11327.501
Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate 0.9600630 8168.08193 11.7993448 1982-1999 –> 2000-2012 1.614205 11151.671
Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 0.9555765 8142.36133 12.5573844 1982-1999 –> 2000-2012 1.601177 11088.679
Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate + silicate 0.9583722 8159.65070 12.3221925 1982-1999 –> 2000-2012 1.610870 11137.718
Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 0.7654124 6829.15714 25.0683926 1982-1999 –> 2000-2012 1.549113 10354.687
Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + phosphate_star 0.9576736 8153.33518 13.3268283 1982-1999 –> 2000-2012 1.669683 11390.207
Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate 0.9653568 8200.62438 13.8491542 1982-1999 –> 2000-2012 1.691710 11497.093
Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 0.9572758 8152.87609 13.2107909 1982-1999 –> 2000-2012 1.668265 11387.463
Atlantic (28.15,28.2] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 0.8182781 7224.40533 13.5294387 1982-1999 –> 2000-2012 1.621210 10822.370
Atlantic (28.15,28.2] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 0.6358764 860.38018 2.3162524 2000-2012 –> 2013-2019 1.596021 9030.969
Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate 0.6284098 848.03313 3.0550910 2000-2012 –> 2013-2019 1.588473 9016.115
Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 0.5673141 760.43649 2.1629643 2000-2012 –> 2013-2019 1.522891 8902.798
Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + nitrate + silicate 0.5172736 679.54551 2.2423627 2000-2012 –> 2013-2019 1.475646 8839.196
Atlantic (28.15,28.2] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 0.4842522 621.75921 1.6393860 2000-2012 –> 2013-2019 1.249665 7450.916
Atlantic (28.15,28.2] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 0.6682826 903.92347 2.0565213 2000-2012 –> 2013-2019 1.627562 9069.174
Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + phosphate_star 0.5689706 760.99054 2.2087785 2000-2012 –> 2013-2019 1.526644 8914.326
Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate 0.5200793 682.28413 2.2713444 2000-2012 –> 2013-2019 1.485436 8882.909
Atlantic (28.15,28.2] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 0.5159584 677.31541 2.3130401 2000-2012 –> 2013-2019 1.473234 8830.191
Atlantic (28.15,28.2] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 0.7477977 1002.40453 1.7226763 2000-2012 –> 2013-2019 1.566076 8226.810
Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate 0.6790836 12178.48882 4.5921194 1982-1999 –> 2000-2012 1.226230 16168.451
Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + phosphate_star 0.6784036 12168.67486 4.5156779 1982-1999 –> 2000-2012 1.196928 15898.543
Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate 0.6615632 11872.26055 4.3245696 1982-1999 –> 2000-2012 1.196177 15751.188
Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 0.6380153 11446.88016 4.6598571 1982-1999 –> 2000-2012 1.145340 15072.234
Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate 0.6786319 12172.64357 4.4942970 1982-1999 –> 2000-2012 1.220970 16121.548
Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 0.6783961 12170.54493 4.5269884 1982-1999 –> 2000-2012 1.189030 15827.608
Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + silicate 0.6572459 11797.05538 4.5957727 1982-1999 –> 2000-2012 1.191738 15676.874
Atlantic (28.2, Inf] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 0.7749638 13739.88315 4.4228123 1982-1999 –> 2000-2012 1.459062 18823.517
Atlantic (28.2, Inf] cstar_tref ~ temp + aou + nitrate + silicate 0.8108738 14272.01531 11.1403370 1982-1999 –> 2000-2012 1.532409 19613.543
Atlantic (28.2, Inf] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 0.7986185 14094.43368 12.1243589 1982-1999 –> 2000-2012 1.510631 19373.154
Atlantic (28.2, Inf] cstar_tref ~ aou + nitrate + silicate + phosphate_star 0.4501223 993.95159 1.5648889 2000-2012 –> 2013-2019 1.258437 15228.692
Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate 0.4547557 1008.15284 2.0609851 2000-2012 –> 2013-2019 1.133839 13186.642
Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + phosphate_star 0.4444906 974.03342 1.9381796 2000-2012 –> 2013-2019 1.122894 13142.708
Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate 0.4449066 975.51344 1.8686116 2000-2012 –> 2013-2019 1.106470 12847.774
Atlantic (28.2, Inf] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 0.4267565 911.62170 1.5887800 2000-2012 –> 2013-2019 1.064772 12358.502
Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate 0.4542792 1008.49415 1.9878249 2000-2012 –> 2013-2019 1.132911 13181.138
Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 0.4378294 952.14581 2.3570728 2000-2012 –> 2013-2019 1.116226 13122.691
Atlantic (28.2, Inf] cstar_tref ~ sal + temp + aou + nitrate + silicate 0.4448646 977.36411 1.8488413 2000-2012 –> 2013-2019 1.102110 12774.419
Atlantic (28.2, Inf] cstar_tref ~ temp + aou + nitrate + silicate 0.4698081 1061.66895 1.5717936 2000-2012 –> 2013-2019 1.280682 15333.684
Atlantic (28.2, Inf] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 0.4182040 879.59547 1.6653982 2000-2012 –> 2013-2019 1.216822 14974.029
Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + phosphate + phosphate_star 7.3104369 31490.51912 39.1711843 1982-1999 –> 2000-2012 13.404252 47059.318
Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 7.0929655 31213.59635 40.5228658 1982-1999 –> 2000-2012 13.040247 46667.027
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 7.5464568 31785.99411 29.0631319 1982-1999 –> 2000-2012 13.440101 47195.739
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + phosphate 7.6411626 31899.18186 32.2400366 1982-1999 –> 2000-2012 13.461879 47246.887
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 7.3103981 31492.47009 39.2906784 1982-1999 –> 2000-2012 13.055907 46779.466
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + phosphate_star 7.7734389 32057.69838 31.9620946 1982-1999 –> 2000-2012 13.734788 47520.521
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 7.3877289 31589.65724 30.9613066 1982-1999 –> 2000-2012 12.897975 46675.049
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 7.4936485 31721.13559 30.7246034 1982-1999 –> 2000-2012 13.093212 46884.062
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 7.7874131 32076.28693 30.8033176 1982-1999 –> 2000-2012 13.745540 47538.503
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 7.5337672 31770.45041 30.5102780 1982-1999 –> 2000-2012 13.163029 46958.883
Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + phosphate + phosphate_star 6.6326824 24320.97907 36.5419814 2000-2012 –> 2013-2019 13.943119 55811.498
Indo-Pacific (-Inf,26] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 6.0966097 23704.22114 33.2248005 2000-2012 –> 2013-2019 13.189575 54917.817
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 6.4719095 24142.82023 35.9956123 2000-2012 –> 2013-2019 14.018366 55928.814
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + nitrate + silicate 6.0000623 23587.02095 46.9934950 2000-2012 –> 2013-2019 13.743833 55611.402
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + phosphate 6.9709903 24686.22874 36.3244617 2000-2012 –> 2013-2019 14.612153 56585.411
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 6.6147557 24303.10838 36.6037799 2000-2012 –> 2013-2019 13.925154 55795.578
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + phosphate_star 7.1024493 24823.39497 36.4083814 2000-2012 –> 2013-2019 14.875888 56881.093
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate 6.4685351 24138.99118 33.0116315 2000-2012 –> 2013-2019 13.856264 55728.648
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 6.5688037 24251.92630 32.9760048 2000-2012 –> 2013-2019 14.062452 55973.062
Indo-Pacific (-Inf,26] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 6.6473334 24339.17899 32.9971359 2000-2012 –> 2013-2019 14.181101 56109.629
Indo-Pacific (26,27] cstar_tref ~ sal + aou + phosphate + phosphate_star 5.9990836 61545.31045 28.1502555 1982-1999 –> 2000-2012 10.460331 91225.567
Indo-Pacific (26,27] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 5.5849096 60176.20385 27.6492836 1982-1999 –> 2000-2012 9.700523 89037.538
Indo-Pacific (26,27] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 6.0345447 61660.26866 34.8318819 1982-1999 –> 2000-2012 10.411459 91148.243
Indo-Pacific (26,27] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 5.9779330 61479.61867 29.3817196 1982-1999 –> 2000-2012 10.327438 90903.645
Indo-Pacific (26,27] cstar_tref ~ sal + temp + aou + silicate + phosphate 5.6458883 60384.33337 31.5029359 1982-1999 –> 2000-2012 9.923785 89639.365
Indo-Pacific (26,27] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5.8028002 60909.73179 31.8362794 1982-1999 –> 2000-2012 10.540927 91204.960
Indo-Pacific (26,27] cstar_tref ~ sal + temp + nitrate + phosphate_star 6.0358861 61662.52863 35.0894963 1982-1999 –> 2000-2012 10.695322 91785.268
Indo-Pacific (26,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 5.7521299 60741.63831 32.6247644 1982-1999 –> 2000-2012 10.384569 90807.222
Indo-Pacific (26,27] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5.7995541 60899.00735 31.7889577 1982-1999 –> 2000-2012 10.555842 91233.181
Indo-Pacific (26,27] cstar_tref ~ sal + temp + silicate + phosphate_star 5.8047303 60914.10545 31.8798539 1982-1999 –> 2000-2012 10.580810 91288.552
Indo-Pacific (26,27] cstar_tref ~ aou + silicate + phosphate + phosphate_star 5.4333408 32421.30831 22.8404576 2000-2012 –> 2013-2019 12.017365 95749.812
Indo-Pacific (26,27] cstar_tref ~ sal + aou + phosphate + phosphate_star 5.8492026 33189.50015 24.6874156 2000-2012 –> 2013-2019 11.848286 94734.811
Indo-Pacific (26,27] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 4.9604904 31474.93572 18.2435343 2000-2012 –> 2013-2019 10.545400 91651.140
Indo-Pacific (26,27] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 5.7652639 33040.94302 24.5362181 2000-2012 –> 2013-2019 11.743197 94520.562
Indo-Pacific (26,27] cstar_tref ~ sal + temp + aou + silicate + phosphate 5.1251300 31821.13763 20.2265627 2000-2012 –> 2013-2019 10.771018 92205.471
Indo-Pacific (26,27] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 5.2258404 32017.72406 20.6713269 2000-2012 –> 2013-2019 11.028641 92927.456
Indo-Pacific (26,27] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 5.1503194 31866.09958 21.5800772 2000-2012 –> 2013-2019 10.902449 92607.738
Indo-Pacific (26,27] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5.2320507 32030.09495 20.7703259 2000-2012 –> 2013-2019 11.031605 92929.102
Indo-Pacific (26,27] cstar_tref ~ sal + temp + silicate + phosphate_star 5.2496504 32063.07371 21.1640102 2000-2012 –> 2013-2019 11.054381 92977.179
Indo-Pacific (26,27] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 5.4332741 32423.18036 22.8037658 2000-2012 –> 2013-2019 11.974337 95628.214
Indo-Pacific (27,27.25] cstar_tref ~ aou + phosphate + phosphate_star 5.2371007 27995.97919 32.2351312 1982-1999 –> 2000-2012 8.173612 39004.063
Indo-Pacific (27,27.25] cstar_tref ~ aou + silicate + phosphate + phosphate_star 4.9933460 27564.16134 27.5072204 1982-1999 –> 2000-2012 7.929622 38573.893
Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + phosphate + phosphate_star 5.1987480 27931.07759 33.6435577 1982-1999 –> 2000-2012 8.050886 38812.712
Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 4.9517171 27489.96100 28.9064215 1982-1999 –> 2000-2012 7.803850 38373.588
Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 5.5502532 28528.58245 33.6306411 1982-1999 –> 2000-2012 8.093087 38906.454
Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 4.2520766 26103.48135 25.8470726 1982-1999 –> 2000-2012 7.059038 36916.768
Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + phosphate_star 5.5535806 28532.03757 34.1239670 1982-1999 –> 2000-2012 8.136183 38976.283
Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5.1784214 27897.41987 29.9525855 1982-1999 –> 2000-2012 7.754068 38331.783
Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + phosphate + phosphate_star 4.2956149 26194.20566 24.4429067 1982-1999 –> 2000-2012 7.160764 37095.894
Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 3.6892995 24811.26247 16.5259999 1982-1999 –> 2000-2012 6.543890 35698.684
Indo-Pacific (27,27.25] cstar_tref ~ aou + phosphate + phosphate_star 5.1106716 15999.51199 20.5011904 2000-2012 –> 2013-2019 10.347772 43995.491
Indo-Pacific (27,27.25] cstar_tref ~ aou + silicate + phosphate + phosphate_star 4.9936597 15880.09789 19.3017367 2000-2012 –> 2013-2019 9.987006 43444.259
Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + phosphate + phosphate_star 5.0186900 15906.30742 21.7121993 2000-2012 –> 2013-2019 10.217438 43837.385
Indo-Pacific (27,27.25] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 4.9323641 15817.35591 20.4566521 2000-2012 –> 2013-2019 9.884081 43307.317
Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 4.4494231 15277.19939 19.5881654 2000-2012 –> 2013-2019 8.701500 41380.681
Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + aou + silicate + phosphate 4.9909008 15879.20105 21.0434228 2000-2012 –> 2013-2019 10.218575 43862.782
Indo-Pacific (27,27.25] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 5.0068556 15895.93181 21.1755589 2000-2012 –> 2013-2019 10.362665 44099.922
Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 5.0972789 15989.75709 18.2443059 2000-2012 –> 2013-2019 10.275700 43887.177
Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + phosphate + phosphate_star 4.5159050 15352.94440 18.5650212 2000-2012 –> 2013-2019 8.811520 41547.150
Indo-Pacific (27,27.25] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 4.0705633 14810.69785 16.0636754 2000-2012 –> 2013-2019 7.759863 39621.960
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 3.3849884 30533.72744 11.0754723 1982-1999 –> 2000-2012 5.952622 43618.610
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + silicate + phosphate 3.7888309 31835.52222 16.1173962 1982-1999 –> 2000-2012 6.437956 45091.314
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 3.5834239 31192.73625 11.8582925 1982-1999 –> 2000-2012 6.202082 44386.513
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + silicate + phosphate_star 3.5906502 31214.04044 11.7127181 1982-1999 –> 2000-2012 6.353720 44702.886
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + silicate + phosphate + phosphate_star 3.6490848 31400.78377 12.8107579 1982-1999 –> 2000-2012 6.374190 44813.064
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 3.5263556 31007.02574 11.3701775 1982-1999 –> 2000-2012 6.173519 44260.717
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 3.5680431 31142.97700 11.8185943 1982-1999 –> 2000-2012 6.284558 44539.784
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 3.5446542 31066.89780 11.5804435 1982-1999 –> 2000-2012 6.246542 44433.826
Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 3.9568587 32339.49177 15.0901136 1982-1999 –> 2000-2012 6.296866 44910.650
Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 3.1123768 29562.43421 18.0064898 1982-1999 –> 2000-2012 5.407434 42026.253
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 3.4559699 18499.65775 12.1867657 2000-2012 –> 2013-2019 6.840958 49033.385
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + silicate + phosphate 3.9438430 19415.68669 13.6147060 2000-2012 –> 2013-2019 7.732674 51251.209
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 3.5505311 18687.32037 14.1874309 2000-2012 –> 2013-2019 7.133955 49880.057
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + aou + silicate + phosphate_star 3.6516395 18880.52595 12.8887909 2000-2012 –> 2013-2019 7.242290 50094.566
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + silicate + phosphate + phosphate_star 3.8388284 19228.06317 12.4331989 2000-2012 –> 2013-2019 7.487913 50628.847
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate 3.4958077 18579.33668 13.7420053 2000-2012 –> 2013-2019 7.022163 49586.362
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 3.5193671 18626.03147 14.0873266 2000-2012 –> 2013-2019 7.087410 49769.008
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate 3.5199729 18625.22792 15.3312305 2000-2012 –> 2013-2019 7.503877 51041.520
Indo-Pacific (27.25,27.5] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 3.4988333 18585.35098 13.8701905 2000-2012 –> 2013-2019 7.043487 49652.249
Indo-Pacific (27.25,27.5] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 3.4147628 18416.26765 16.3339820 2000-2012 –> 2013-2019 6.527140 47978.702
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 3.2112424 35141.91301 16.4034407 1982-1999 –> 2000-2012 5.275493 48571.226
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 3.0250790 34330.54803 11.0114608 1982-1999 –> 2000-2012 5.292166 48346.413
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate_star 3.2504188 35304.65632 15.4799081 1982-1999 –> 2000-2012 5.536903 49371.835
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + nitrate + silicate 3.2004156 35096.03017 12.5321673 1982-1999 –> 2000-2012 5.560490 49363.449
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 2.9851108 34149.84916 10.3732193 1982-1999 –> 2000-2012 5.237249 48124.312
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2.9944106 34192.10925 10.5059002 1982-1999 –> 2000-2012 5.280726 48260.826
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 3.4398779 36076.33199 20.5231614 1982-1999 –> 2000-2012 5.527895 49577.283
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate 2.9856060 34150.10289 10.4099816 1982-1999 –> 2000-2012 5.362578 48460.169
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 2.9853761 34151.05676 10.3885141 1982-1999 –> 2000-2012 5.265905 48203.914
Indo-Pacific (27.5,27.75] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 3.2392001 35259.68360 26.8858334 1982-1999 –> 2000-2012 5.473127 49183.337
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 3.5326504 21231.33076 12.1447892 2000-2012 –> 2013-2019 6.743893 56373.244
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2.8161137 19437.30143 12.9190676 2000-2012 –> 2013-2019 5.841193 53767.849
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + aou + silicate + phosphate_star 3.5789089 21332.28839 11.9951356 2000-2012 –> 2013-2019 6.829328 56636.945
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + nitrate + silicate 2.8253648 19463.25657 13.8804937 2000-2012 –> 2013-2019 6.025780 54559.287
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate 2.8440045 19513.29591 14.4603188 2000-2012 –> 2013-2019 6.103330 54855.128
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate 2.8006015 19393.58751 12.3250577 2000-2012 –> 2013-2019 5.785712 53543.437
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2.8014773 19396.06200 12.3871862 2000-2012 –> 2013-2019 5.795888 53588.171
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate 3.1329458 20279.05940 9.0077237 2000-2012 –> 2013-2019 6.118552 54429.162
Indo-Pacific (27.5,27.75] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 2.7987433 19388.33491 12.2484203 2000-2012 –> 2013-2019 5.784119 53539.392
Indo-Pacific (27.5,27.75] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 3.5994976 21379.68552 13.8983258 2000-2012 –> 2013-2019 6.838698 56639.369
Indo-Pacific (27.75,27.85] cstar_tref ~ aou + silicate + phosphate + phosphate_star 3.0485778 13754.32087 21.2149937 1982-1999 –> 2000-2012 4.853759 18959.127
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2.9864839 13644.70331 21.7452820 1982-1999 –> 2000-2012 4.788392 18846.820
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 3.0709187 13795.92461 17.4217391 1982-1999 –> 2000-2012 4.775258 18854.192
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + silicate 3.0312326 13725.37236 19.8760952 1982-1999 –> 2000-2012 4.866540 18974.946
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 3.2096685 14035.61594 20.6600951 1982-1999 –> 2000-2012 4.807580 18927.268
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 2.9846650 13641.39879 21.6202029 1982-1999 –> 2000-2012 4.839816 18918.760
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2.9844286 13640.96925 21.5584275 1982-1999 –> 2000-2012 4.862650 18950.266
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 2.9844334 13640.97800 21.5247014 1982-1999 –> 2000-2012 4.867761 18957.291
Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 3.0254827 13715.07388 20.0994248 1982-1999 –> 2000-2012 4.763970 18824.602
Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 3.0460606 13751.84039 20.9920953 1982-1999 –> 2000-2012 4.657167 18664.742
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 2.5277333 7470.41912 8.4417890 2000-2012 –> 2013-2019 5.514217 21115.122
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + nitrate + silicate 2.8386769 7839.11498 10.9077148 2000-2012 –> 2013-2019 5.869910 21564.487
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate 2.5612931 7512.33454 7.5788589 2000-2012 –> 2013-2019 5.545958 21153.733
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 2.5700955 7523.23770 7.5816681 2000-2012 –> 2013-2019 5.554524 21164.207
Indo-Pacific (27.75,27.85] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 2.5772559 7532.07947 7.5699017 2000-2012 –> 2013-2019 5.561689 21173.057
Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2.9503252 7961.71377 9.3984757 2000-2012 –> 2013-2019 5.975808 21676.788
Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + silicate + phosphate 2.9565757 7966.43951 9.5367896 2000-2012 –> 2013-2019 6.004606 21719.787
Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 2.9562916 7968.13415 9.5578677 2000-2012 –> 2013-2019 6.002352 21719.975
Indo-Pacific (27.75,27.85] cstar_tref ~ temp + aou + silicate + phosphate_star 2.9570397 7966.93829 9.5252980 2000-2012 –> 2013-2019 6.004238 21718.804
Indo-Pacific (27.75,27.85] cstar_tref ~ temp + silicate + phosphate + phosphate_star 2.9576726 7967.61831 9.5120200 2000-2012 –> 2013-2019 6.004479 21718.787
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate 2.9763726 16635.88090 12.8739590 1982-1999 –> 2000-2012 4.834131 22432.602
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 2.9753738 16635.65767 12.7056086 1982-1999 –> 2000-2012 4.743923 22294.718
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate 3.0119732 16714.64077 12.4248859 1982-1999 –> 2000-2012 4.801556 22405.255
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 3.0094926 16711.18311 12.2412154 1982-1999 –> 2000-2012 4.798584 22403.017
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + silicate 2.9762974 16637.71349 12.8420797 1982-1999 –> 2000-2012 4.717297 22252.217
Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate 3.0214813 16733.51844 12.1756810 1982-1999 –> 2000-2012 4.895962 22553.672
Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate + phosphate_star 3.0197956 16731.82184 12.0306931 1982-1999 –> 2000-2012 4.890372 22548.057
Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate + silicate 2.9903551 16666.92636 12.4996275 1982-1999 –> 2000-2012 4.819031 22418.868
Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2.9897076 16667.49202 12.3975301 1982-1999 –> 2000-2012 4.818041 22420.902
Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 3.2266639 17172.72536 19.0200577 1982-1999 –> 2000-2012 4.869583 22622.668
Indo-Pacific (27.85,27.95] cstar_tref ~ aou + nitrate + silicate + phosphate_star 2.9540891 9485.05992 10.2048048 2000-2012 –> 2013-2019 5.956014 26177.566
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + aou + nitrate 2.9190648 9437.90416 11.0080109 2000-2012 –> 2013-2019 5.934044 26157.153
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + aou + nitrate + phosphate_star 2.8453506 9343.06963 14.2670774 2000-2012 –> 2013-2019 5.860329 26064.317
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate 2.9148297 9434.40728 10.8783798 2000-2012 –> 2013-2019 5.891202 26070.288
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 2.8201826 9311.43225 14.5544280 2000-2012 –> 2013-2019 5.795556 25947.090
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate 2.8538540 9354.36730 11.6329621 2000-2012 –> 2013-2019 5.865827 26069.008
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 2.8373081 9334.35312 13.4974823 2000-2012 –> 2013-2019 5.846801 26045.536
Indo-Pacific (27.85,27.95] cstar_tref ~ sal + temp + aou + nitrate + silicate 2.7860076 9265.27328 11.4499589 2000-2012 –> 2013-2019 5.762305 25902.987
Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate + silicate 2.9519650 9482.33668 10.3208135 2000-2012 –> 2013-2019 5.942320 26149.263
Indo-Pacific (27.85,27.95] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 2.9406968 9469.85715 11.9349933 2000-2012 –> 2013-2019 5.930404 26137.349
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1.7874554 15895.86831 7.8151143 1982-1999 –> 2000-2012 3.638642 23051.895
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1.7286768 15630.31324 7.7259247 1982-1999 –> 2000-2012 3.657598 22930.721
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate 1.7180393 15579.29062 8.8610580 1982-1999 –> 2000-2012 3.635996 22857.688
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 1.7131939 15558.86034 8.4375330 1982-1999 –> 2000-2012 3.624809 22827.633
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate + silicate 1.6667019 15340.35497 8.1214416 1982-1999 –> 2000-2012 3.490450 22443.967
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + phosphate + phosphate_star 1.7505669 15730.25071 8.8790256 1982-1999 –> 2000-2012 3.652754 22981.669
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 1.7394473 15679.64225 7.5925193 1982-1999 –> 2000-2012 3.686232 23012.406
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1.7423080 15692.69315 7.5165669 1982-1999 –> 2000-2012 3.693224 23032.898
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 1.7432650 15697.05391 7.4945473 1982-1999 –> 2000-2012 3.694948 23038.638
Indo-Pacific (27.95,28.05] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1.7164828 15574.09216 8.1537405 1982-1999 –> 2000-2012 3.414315 22426.593
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1.5914479 8591.83603 5.6830274 2000-2012 –> 2013-2019 3.378903 24487.704
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1.6390817 8726.14212 6.3665459 2000-2012 –> 2013-2019 3.367758 24356.455
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate 1.7296062 8968.95424 8.8832430 2000-2012 –> 2013-2019 3.447646 24548.245
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate + phosphate_star 1.6881292 8860.41581 7.9824969 2000-2012 –> 2013-2019 3.401323 24419.276
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + nitrate + silicate 1.6136967 8655.06105 7.1018200 2000-2012 –> 2013-2019 3.280399 23995.416
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate 1.6404783 8730.02086 6.3193490 2000-2012 –> 2013-2019 3.379926 24409.663
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1.6406708 8730.55496 6.3060246 2000-2012 –> 2013-2019 3.382979 24423.248
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 1.6313849 8704.70716 5.9735729 2000-2012 –> 2013-2019 3.399475 24514.063
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 1.6406058 8730.37470 6.3030216 2000-2012 –> 2013-2019 3.383871 24427.429
Indo-Pacific (27.95,28.05] cstar_tref ~ sal + temp + silicate + phosphate_star 1.6431153 8735.33511 6.6113621 2000-2012 –> 2013-2019 3.422819 24594.686
Indo-Pacific (28.05,28.1] cstar_tref ~ aou + silicate + phosphate + phosphate_star 1.6447239 8678.46306 7.7207500 1982-1999 –> 2000-2012 3.219938 12789.308
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1.5915015 8531.71352 7.7621673 1982-1999 –> 2000-2012 3.145566 12614.983
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + silicate 1.6218520 8617.13784 8.3281132 1982-1999 –> 2000-2012 3.228276 12772.911
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 1.5956884 8543.59422 7.8212458 1982-1999 –> 2000-2012 3.159611 12640.699
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1.5967017 8546.46500 7.8401506 1982-1999 –> 2000-2012 3.164411 12648.861
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 1.5986906 8552.09413 7.9440822 1982-1999 –> 2000-2012 3.155713 12639.524
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 1.5964228 8545.67518 7.8335931 1982-1999 –> 2000-2012 3.163050 12646.560
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate_star 1.6131414 8590.78564 8.1419640 1982-1999 –> 2000-2012 3.219929 12745.054
Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1.6551150 8708.94237 7.6426135 1982-1999 –> 2000-2012 3.248519 12846.909
Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1.5817655 8503.96540 6.9776970 1982-1999 –> 2000-2012 3.026668 12427.878
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1.3027624 4273.06645 5.5980153 2000-2012 –> 2013-2019 2.894264 12804.780
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + nitrate + silicate 1.3130657 4292.99709 5.2892135 2000-2012 –> 2013-2019 2.934918 12910.135
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate 1.3058438 4279.04349 5.5563164 2000-2012 –> 2013-2019 2.901532 12822.638
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1.3066744 4280.65235 5.5448453 2000-2012 –> 2013-2019 2.903376 12827.117
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate 1.3519556 4364.84136 4.9525298 2000-2012 –> 2013-2019 3.045685 13176.072
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 1.3057827 4278.92521 5.5444741 2000-2012 –> 2013-2019 2.904473 12831.019
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate 1.3457180 4353.14158 4.9632050 2000-2012 –> 2013-2019 3.021116 13115.164
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 1.3064165 4280.15278 5.5391203 2000-2012 –> 2013-2019 2.902839 12825.828
Indo-Pacific (28.05,28.1] cstar_tref ~ sal + temp + silicate + phosphate_star 1.3263075 4316.38343 6.4976889 2000-2012 –> 2013-2019 2.939449 12907.169
Indo-Pacific (28.05,28.1] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1.3125511 4292.00522 6.4087236 2000-2012 –> 2013-2019 2.894317 12795.971
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1.5883171 89371.82695 9.7632536 1982-1999 –> 2000-2012 3.234197 135521.653
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + aou + silicate + phosphate 1.5868540 89326.06037 10.0894722 1982-1999 –> 2000-2012 3.005214 131893.792
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1.5449975 88058.59952 9.5563055 1982-1999 –> 2000-2012 2.860038 128808.291
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + temp + aou + silicate + phosphate 1.5451389 88062.94641 9.5573889 1982-1999 –> 2000-2012 2.859391 128798.203
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1.5452865 88067.48241 9.5583763 1982-1999 –> 2000-2012 2.858544 128784.508
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 1.5517841 88266.74970 9.2817857 1982-1999 –> 2000-2012 2.879350 129244.542
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 1.5452757 88067.15154 9.5586239 1982-1999 –> 2000-2012 2.858578 128785.003
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + temp + silicate + phosphate_star 1.5519083 88268.54903 9.2248545 1982-1999 –> 2000-2012 2.881166 129275.008
Indo-Pacific (28.1, Inf] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1.4895568 86323.14245 7.4525907 1982-1999 –> 2000-2012 2.806933 127115.521
Indo-Pacific (28.1, Inf] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1.5819058 89179.74622 9.4623053 1982-1999 –> 2000-2012 3.037218 132368.396
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + aou + nitrate + silicate + phosphate_star 1.4706717 53633.99233 8.6400885 2000-2012 –> 2013-2019 3.058989 143005.819
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + aou + silicate + phosphate 1.3484173 51053.35522 8.6511037 2000-2012 –> 2013-2019 2.935271 140379.416
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + aou + silicate + phosphate + phosphate_star 1.2967699 49894.95069 9.4879928 2000-2012 –> 2013-2019 2.841767 137953.550
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + temp + aou + silicate + phosphate 1.2971098 49902.73877 9.4884111 2000-2012 –> 2013-2019 2.842249 137965.685
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + temp + aou + silicate + phosphate_star 1.2975644 49913.14971 9.4916336 2000-2012 –> 2013-2019 2.842851 137980.632
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + temp + nitrate + silicate + phosphate_star 1.3092903 50180.44731 9.7188452 2000-2012 –> 2013-2019 2.861074 138447.197
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + temp + silicate + phosphate + phosphate_star 1.2975181 49912.08927 9.4905243 2000-2012 –> 2013-2019 2.842794 137979.241
Indo-Pacific (28.1, Inf] cstar_tref ~ sal + temp + silicate + phosphate_star 1.3255162 50544.40212 10.3689998 2000-2012 –> 2013-2019 2.877425 138812.951
Indo-Pacific (28.1, Inf] cstar_tref ~ temp + aou + nitrate + silicate + phosphate_star 1.4289536 52778.97514 11.7985403 2000-2012 –> 2013-2019 2.918510 139102.118
Indo-Pacific (28.1, Inf] cstar_tref ~ temp + aou + silicate + phosphate + phosphate_star 1.3513288 51119.43989 9.9488153 2000-2012 –> 2013-2019 2.933235 140299.186

5.3 Target variable coefficients

A data frame to map the target variable is prepared.

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

eras_backward <- lm_all_fitted_wide %>%
  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_target <- full_join(
  lm_best %>% select(basin, gamma_slab, model, eras),
  eras_era)

lm_best_target <- left_join(lm_best_target, lm_all_fitted_wide)

rm(eras_era, eras_forward, eras_backward,
   lm_all_fitted)

5.4 Plot selected model residuals

# plot model diagnostics, if activated
if (params_local$plot_all_figures == "y") {
  # mutate predictors column
  lm_best_plot <- lm_best_target %>%
    select(basin, gamma_slab, model, eras, era) %>%
    mutate(
      predictors = str_remove(model, paste(params_local$MLR_target, "~ ")),
      predictors = str_replace_all(predictors, "\\+ ", "")
    )
  
  # loop across all basins, era, gamma slabs, and MLRs
  # fit all MLR models
  for (i_basin in unique(GLODAP$basin)) {
    for (i_era in unique(GLODAP$era)) {
      #i_basin <- unique(GLODAP$basin)[1]
      #i_era   <- unique(GLODAP$era)[2]
      print(i_basin)
      print(i_era)
      
      GLODAP_basin_era <- GLODAP %>%
        filter(basin == i_basin, era == i_era)
      
      for (i_gamma_slab in unique(GLODAP_basin_era$gamma_slab)) {
        #i_gamma_slab <- unique(GLODAP_basin_era$gamma_slab)[1]
        print(i_gamma_slab)
        
        GLODAP_basin_era_slab <- GLODAP_basin_era %>%
          filter(gamma_slab == i_gamma_slab)
        
        lm_best_basin_era_slab <- lm_best_plot %>%
          filter(basin == i_basin, era == i_era, gamma_slab == i_gamma_slab)
        
        for (i_eras in unique(lm_best_basin_era_slab$eras)) {
          #i_eras <- unique(lm_best_basin_era_slab$eras)[1]
          print(i_eras)
          
          lm_best_basin_era_slab_eras <- lm_best_basin_era_slab %>%
            filter(eras == i_eras)
          
          for (i_predictors in unique(lm_best_basin_era_slab_eras$predictors)) {
            #i_predictors <- unique(lm_best_basin_era_slab$predictors)[1]
            print(i_predictors)
            # extract one model definition
            i_lm <- lm_all %>%
              filter(predictors == i_predictors) %>%
              select(lm_coeff) %>%
              pull()
            
            # fit model
            i_lm_fit <- lm(as.formula(i_lm),
                           data = GLODAP_basin_era_slab)
            
            p_model <- ggnostic(
              i_lm_fit,
              columnsY = c(params_local$MLR_target, ".fitted", ".resid"),
              title = paste(
                "era:",
                i_era,
                "| eras:",
                i_eras,
                "| basin:",
                i_basin,
                "| gamma slab:",
                i_gamma_slab,
                "| predictors:",
                i_predictors
              )
            )
            
            ggsave(
              plot = p_model,
              path = paste(path_version_figures, "eMLR_diagnostics/", sep = ""),
              filename = paste(
                "MLR_residuals",
                i_era,
                i_eras,
                i_basin,
                i_gamma_slab,
                i_predictors,
                "predictors.png",
                sep = "_"
              ),
              width = 14,
              height = 8
            )

            rm(p_model)
            
          }
          
        }
        
      }
    }
  }
  
  rm(
    lm_best_plot,
    lm_best_basin_era_slab,
    i_rmse,
    GLODAP_basin_era,
    GLODAP_basin_era_slab,
    i_lm,
    lm_all_fitted,
    i_basin,
    i_era,
    i_gamma_slab,
    i_predictors,
    lm_all,
    i_aic,
    i_n_predictors,
    i_resid_max
  )
  
}

Individual residual plots of the MLR models for each basin, era, eras and neutral density (gamma) slab are available at:

/nfs/kryo/work/jenmueller/emlr_cant/model/v_XXX/figures/eMLR_diagnostics/

5.5 Cant coeffcients

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

# pivot long format
lm_best_long <- lm_best_target %>%
  pivot_longer(cols = starts_with("coeff_"),
               names_to = "term",
               values_to = "estimate",
               names_prefix = "coeff_")

# subtract coefficients of adjacent era  
lm_best_long <- lm_best_long %>%
  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 back to wide format
lm_best_cant <- lm_best_long %>%
  pivot_wider(values_from = delta_coeff,
              names_from = term,
              names_prefix = "delta_coeff_",
              values_fill = 0)

5.6 Write files

lm_best_target %>%
  select(
    basin,
    gamma_slab,
    model,
    eras,
    era,
    starts_with("coeff_")
  ) %>%
  write_csv(paste(path_version_data,
                  "lm_best_target.csv",
                  sep = ""))

lm_best_cant %>%
  select(
    basin,
    gamma_slab,
    model,
    eras,
    starts_with("delta_coeff_")
  ) %>%
  write_csv(paste(path_version_data,
                  "lm_best_cant.csv",
                  sep = ""))

6 Model diagnotics

6.1 Selection criterion vs predictors

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

6.1.1 All models

lm_all_fitted_wide %>%
  ggplot(aes(as.factor(n_predictors),
             !!sym(params_local$MLR_criterion),
             col = basin)) +
  geom_hline(yintercept = 10) +
  geom_boxplot() +
  facet_grid(gamma_slab~era) +
  scale_color_brewer(palette = "Set1") +
  labs(x="Number of predictors")

Version Author Date
a618965 Donghe-Zhu 2021-01-31
59e006e Donghe-Zhu 2021-01-31
a1c8f87 Donghe-Zhu 2021-01-31
ae5c18f Donghe-Zhu 2021-01-31
b50fe52 Donghe-Zhu 2021-01-31
ac99ae5 jens-daniel-mueller 2021-01-29
b5bdcaf Donghe-Zhu 2021-01-29
442010d Donghe-Zhu 2021-01-29
372adf5 Donghe-Zhu 2021-01-29
af8788e Donghe-Zhu 2021-01-29
21c91c9 Donghe-Zhu 2021-01-29
eded038 Donghe-Zhu 2021-01-29
541d4dd Donghe-Zhu 2021-01-29
6a75576 Donghe-Zhu 2021-01-28
16fba40 Donghe-Zhu 2021-01-28
12bc567 Donghe-Zhu 2021-01-27
ceed31b Donghe-Zhu 2021-01-27
342402d Donghe-Zhu 2021-01-27
5bad5c2 Donghe-Zhu 2021-01-27
61efb56 Donghe-Zhu 2021-01-25
48f638e Donghe-Zhu 2021-01-25
c1cec47 Donghe-Zhu 2021-01-25
05ffb0c Donghe-Zhu 2021-01-25
8b97165 Donghe-Zhu 2021-01-25
c569946 Donghe-Zhu 2021-01-24
a2f0d56 Donghe-Zhu 2021-01-23
28509fc Donghe-Zhu 2021-01-23
4c28e4a Donghe-Zhu 2021-01-22
24cc264 jens-daniel-mueller 2021-01-22
7891955 Donghe-Zhu 2021-01-21
d4cf1cb Donghe-Zhu 2021-01-21
1f3e5b6 jens-daniel-mueller 2021-01-20
0e7bdf1 jens-daniel-mueller 2021-01-15
4571843 jens-daniel-mueller 2021-01-14
b3564aa jens-daniel-mueller 2021-01-14
8d032c3 jens-daniel-mueller 2021-01-14
17dee1d jens-daniel-mueller 2021-01-13
7cdea0c jens-daniel-mueller 2021-01-06
fa85b93 jens-daniel-mueller 2021-01-06
e5cb81a Donghe-Zhu 2021-01-05
a499f10 Donghe-Zhu 2021-01-05
fb8a752 Donghe-Zhu 2020-12-23
8fae0b2 Donghe-Zhu 2020-12-21
c8b76b3 jens-daniel-mueller 2020-12-19

6.1.2 Best models

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

Version Author Date
a618965 Donghe-Zhu 2021-01-31
59e006e Donghe-Zhu 2021-01-31
a1c8f87 Donghe-Zhu 2021-01-31
ae5c18f Donghe-Zhu 2021-01-31
b50fe52 Donghe-Zhu 2021-01-31
ac99ae5 jens-daniel-mueller 2021-01-29
b5bdcaf Donghe-Zhu 2021-01-29
442010d Donghe-Zhu 2021-01-29
372adf5 Donghe-Zhu 2021-01-29
af8788e Donghe-Zhu 2021-01-29
21c91c9 Donghe-Zhu 2021-01-29
eded038 Donghe-Zhu 2021-01-29
541d4dd Donghe-Zhu 2021-01-29
6a75576 Donghe-Zhu 2021-01-28
16fba40 Donghe-Zhu 2021-01-28
12bc567 Donghe-Zhu 2021-01-27
ceed31b Donghe-Zhu 2021-01-27
342402d Donghe-Zhu 2021-01-27
5bad5c2 Donghe-Zhu 2021-01-27
61efb56 Donghe-Zhu 2021-01-25
48f638e Donghe-Zhu 2021-01-25
c1cec47 Donghe-Zhu 2021-01-25
05ffb0c Donghe-Zhu 2021-01-25
8b97165 Donghe-Zhu 2021-01-25
c569946 Donghe-Zhu 2021-01-24
a2f0d56 Donghe-Zhu 2021-01-23
28509fc Donghe-Zhu 2021-01-23
4c28e4a Donghe-Zhu 2021-01-22
24cc264 jens-daniel-mueller 2021-01-22
7891955 Donghe-Zhu 2021-01-21
d4cf1cb Donghe-Zhu 2021-01-21
1f3e5b6 jens-daniel-mueller 2021-01-20
0e7bdf1 jens-daniel-mueller 2021-01-15
4571843 jens-daniel-mueller 2021-01-14
b3564aa jens-daniel-mueller 2021-01-14
8d032c3 jens-daniel-mueller 2021-01-14
17dee1d jens-daniel-mueller 2021-01-13
7cdea0c jens-daniel-mueller 2021-01-06
fa85b93 jens-daniel-mueller 2021-01-06
e5cb81a Donghe-Zhu 2021-01-05
a499f10 Donghe-Zhu 2021-01-05
fb8a752 Donghe-Zhu 2020-12-23
8fae0b2 Donghe-Zhu 2020-12-21
c8b76b3 jens-daniel-mueller 2020-12-19

6.2 RMSE correlation between eras

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

6.2.1 All models

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

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

Version Author Date
a618965 Donghe-Zhu 2021-01-31
59e006e Donghe-Zhu 2021-01-31
a1c8f87 Donghe-Zhu 2021-01-31
ae5c18f Donghe-Zhu 2021-01-31
b50fe52 Donghe-Zhu 2021-01-31
ac99ae5 jens-daniel-mueller 2021-01-29
b5bdcaf Donghe-Zhu 2021-01-29
442010d Donghe-Zhu 2021-01-29
372adf5 Donghe-Zhu 2021-01-29
af8788e Donghe-Zhu 2021-01-29
21c91c9 Donghe-Zhu 2021-01-29
eded038 Donghe-Zhu 2021-01-29
541d4dd Donghe-Zhu 2021-01-29
6a75576 Donghe-Zhu 2021-01-28
16fba40 Donghe-Zhu 2021-01-28
12bc567 Donghe-Zhu 2021-01-27
ceed31b Donghe-Zhu 2021-01-27
342402d Donghe-Zhu 2021-01-27
5bad5c2 Donghe-Zhu 2021-01-27
61efb56 Donghe-Zhu 2021-01-25
48f638e Donghe-Zhu 2021-01-25
c1cec47 Donghe-Zhu 2021-01-25
05ffb0c Donghe-Zhu 2021-01-25
8b97165 Donghe-Zhu 2021-01-25
c569946 Donghe-Zhu 2021-01-24
a2f0d56 Donghe-Zhu 2021-01-23
28509fc Donghe-Zhu 2021-01-23
4c28e4a Donghe-Zhu 2021-01-22
24cc264 jens-daniel-mueller 2021-01-22
7891955 Donghe-Zhu 2021-01-21
d4cf1cb Donghe-Zhu 2021-01-21
1f3e5b6 jens-daniel-mueller 2021-01-20
0e7bdf1 jens-daniel-mueller 2021-01-15
4571843 jens-daniel-mueller 2021-01-14
b3564aa jens-daniel-mueller 2021-01-14
8d032c3 jens-daniel-mueller 2021-01-14
17dee1d jens-daniel-mueller 2021-01-13
7cdea0c jens-daniel-mueller 2021-01-06
fa85b93 jens-daniel-mueller 2021-01-06
e5cb81a Donghe-Zhu 2021-01-05
a499f10 Donghe-Zhu 2021-01-05
fb8a752 Donghe-Zhu 2020-12-23
8fae0b2 Donghe-Zhu 2020-12-21
c8b76b3 jens-daniel-mueller 2020-12-19 
rm(max_rmse)

6.2.2 Best models

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

lm_best %>%
  ggplot(aes(rmse, rmse_sum - rmse, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  coord_equal(xlim = c(0,max_rmse),
              ylim = c(0,max_rmse)) +
  facet_grid(eras ~ basin)

Version Author Date
a618965 Donghe-Zhu 2021-01-31
59e006e Donghe-Zhu 2021-01-31
a1c8f87 Donghe-Zhu 2021-01-31
ae5c18f Donghe-Zhu 2021-01-31
b50fe52 Donghe-Zhu 2021-01-31
ac99ae5 jens-daniel-mueller 2021-01-29
b5bdcaf Donghe-Zhu 2021-01-29
442010d Donghe-Zhu 2021-01-29
372adf5 Donghe-Zhu 2021-01-29
af8788e Donghe-Zhu 2021-01-29
21c91c9 Donghe-Zhu 2021-01-29
eded038 Donghe-Zhu 2021-01-29
541d4dd Donghe-Zhu 2021-01-29
6a75576 Donghe-Zhu 2021-01-28
16fba40 Donghe-Zhu 2021-01-28
12bc567 Donghe-Zhu 2021-01-27
ceed31b Donghe-Zhu 2021-01-27
342402d Donghe-Zhu 2021-01-27
5bad5c2 Donghe-Zhu 2021-01-27
61efb56 Donghe-Zhu 2021-01-25
48f638e Donghe-Zhu 2021-01-25
c1cec47 Donghe-Zhu 2021-01-25
05ffb0c Donghe-Zhu 2021-01-25
8b97165 Donghe-Zhu 2021-01-25
c569946 Donghe-Zhu 2021-01-24
a2f0d56 Donghe-Zhu 2021-01-23
28509fc Donghe-Zhu 2021-01-23
4c28e4a Donghe-Zhu 2021-01-22
24cc264 jens-daniel-mueller 2021-01-22
7891955 Donghe-Zhu 2021-01-21
d4cf1cb Donghe-Zhu 2021-01-21
1f3e5b6 jens-daniel-mueller 2021-01-20
0e7bdf1 jens-daniel-mueller 2021-01-15
4571843 jens-daniel-mueller 2021-01-14
b3564aa jens-daniel-mueller 2021-01-14
8d032c3 jens-daniel-mueller 2021-01-14
17dee1d jens-daniel-mueller 2021-01-13
7cdea0c jens-daniel-mueller 2021-01-06
fa85b93 jens-daniel-mueller 2021-01-06
e5cb81a Donghe-Zhu 2021-01-05
a499f10 Donghe-Zhu 2021-01-05
fb8a752 Donghe-Zhu 2020-12-23
8fae0b2 Donghe-Zhu 2020-12-21
c8b76b3 jens-daniel-mueller 2020-12-19 
rm(max_rmse)

6.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_long %>% 
  filter(term != "(Intercept)",
         delta_coeff != 0) %>% 
  group_by(basin, eras, gamma_slab) %>% 
  count(term) %>% 
  ungroup() %>% 
  pivot_wider(values_from = n,
              names_from = term,
              values_fill = 0)

# 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 2 6 9 9 7 8
(26,26.5] 9 3 6 8 5 8 8
(26.5,27.25] 9 1 7 7 5 10 7
(27.25,27.5] 7 4 4 6 8 10 8
(27.5,27.75] 6 2 6 8 10 6 9
(27.75,27.85] 7 2 6 9 8 7 9
(27.85,27.95] 7 2 6 9 8 7 9
(27.95,28.05] 7 5 4 9 7 7 9
(28.05,28.1] 9 1 8 8 7 6 8
(28.1,28.15] 9 1 8 8 5 8 8
(28.15,28.2] 10 8 2 7 6 5 8
(28.2, Inf] 9 10 0 5 8 6 6
total 96.00 41.00 63.00 93.00 86.00 87.00 97.00
Atlantic - 2000-2012 --> 2013-2019
(-Inf,26] 5 2 4 4 10 5 9
(26,26.5] 9 6 3 6 7 7 9
(26.5,27.25] 9 3 6 6 6 10 7
(27.25,27.5] 7 6 2 6 8 10 8
(27.5,27.75] 7 1 7 7 10 7 7
(27.75,27.85] 9 3 7 9 5 6 8
(27.85,27.95] 8 3 6 9 7 7 9
(27.95,28.05] 9 3 6 9 6 6 9
(28.05,28.1] 9 2 8 8 7 5 8
(28.1,28.15] 9 1 8 9 6 7 7
(28.15,28.2] 9 7 3 7 6 6 9
(28.2, Inf] 10 10 0 5 7 6 5
total 100.00 47.00 60.00 85.00 85.00 82.00 95.00
Indo-Pacific - 1982-1999 --> 2000-2012
(-Inf,26] 8 2 6 8 10 5 8
(26,27] 6 3 5 9 10 6 8
(27,27.25] 10 3 7 10 4 4 6
(27.25,27.5] 8 2 6 8 8 10 5
(27.5,27.75] 7 3 5 7 9 10 7
(27.75,27.85] 9 3 6 8 7 8 8
(27.85,27.95] 10 9 1 5 5 6 8
(27.95,28.05] 9 4 5 7 9 7 8
(28.05,28.1] 7 3 5 8 7 10 8
(28.1, Inf] 7 3 5 8 8 10 7
total 81.00 35.00 51.00 78.00 77.00 76.00 73.00
Indo-Pacific - 2000-2012 --> 2013-2019
(-Inf,26] 9 2 6 7 10 5 8
(26,27] 7 1 7 9 8 8 7
(27,27.25] 9 1 9 9 5 6 6
(27.25,27.5] 7 1 7 7 9 10 5
(27.5,27.75] 8 2 5 6 9 10 7
(27.75,27.85] 8 2 6 7 5 10 9
(27.85,27.95] 10 10 0 5 7 6 5
(27.95,28.05] 7 5 3 7 10 8 8
(28.05,28.1] 5 3 5 6 9 10 9
(28.1, Inf] 7 3 5 8 8 10 7
total 77.00 30.00 53.00 71.00 80.00 83.00 71.00

6.4 RMSE alternatives

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

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

Version Author Date
a618965 Donghe-Zhu 2021-01-31
59e006e Donghe-Zhu 2021-01-31
a1c8f87 Donghe-Zhu 2021-01-31
ae5c18f Donghe-Zhu 2021-01-31
b50fe52 Donghe-Zhu 2021-01-31
ac99ae5 jens-daniel-mueller 2021-01-29
b5bdcaf Donghe-Zhu 2021-01-29
442010d Donghe-Zhu 2021-01-29
372adf5 Donghe-Zhu 2021-01-29
af8788e Donghe-Zhu 2021-01-29
21c91c9 Donghe-Zhu 2021-01-29
eded038 Donghe-Zhu 2021-01-29
541d4dd Donghe-Zhu 2021-01-29
6a75576 Donghe-Zhu 2021-01-28
16fba40 Donghe-Zhu 2021-01-28
12bc567 Donghe-Zhu 2021-01-27
ceed31b Donghe-Zhu 2021-01-27
342402d Donghe-Zhu 2021-01-27
5bad5c2 Donghe-Zhu 2021-01-27
61efb56 Donghe-Zhu 2021-01-25
48f638e Donghe-Zhu 2021-01-25
c1cec47 Donghe-Zhu 2021-01-25
05ffb0c Donghe-Zhu 2021-01-25
8b97165 Donghe-Zhu 2021-01-25
c569946 Donghe-Zhu 2021-01-24
a2f0d56 Donghe-Zhu 2021-01-23
28509fc Donghe-Zhu 2021-01-23
4c28e4a Donghe-Zhu 2021-01-22
24cc264 jens-daniel-mueller 2021-01-22
7891955 Donghe-Zhu 2021-01-21
d4cf1cb Donghe-Zhu 2021-01-21
1f3e5b6 jens-daniel-mueller 2021-01-20
0e7bdf1 jens-daniel-mueller 2021-01-15
4571843 jens-daniel-mueller 2021-01-14
b3564aa jens-daniel-mueller 2021-01-14
8d032c3 jens-daniel-mueller 2021-01-14
17dee1d jens-daniel-mueller 2021-01-13
7cdea0c jens-daniel-mueller 2021-01-06
fa85b93 jens-daniel-mueller 2021-01-06
e5cb81a Donghe-Zhu 2021-01-05
a499f10 Donghe-Zhu 2021-01-05
fb8a752 Donghe-Zhu 2020-12-23
8fae0b2 Donghe-Zhu 2020-12-21
c8b76b3 jens-daniel-mueller 2020-12-19 
lm_best %>% 
  ggplot(aes(rmse, aic, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  facet_grid(eras~basin)

Version Author Date
a618965 Donghe-Zhu 2021-01-31
59e006e Donghe-Zhu 2021-01-31
a1c8f87 Donghe-Zhu 2021-01-31
ae5c18f Donghe-Zhu 2021-01-31
b50fe52 Donghe-Zhu 2021-01-31
ac99ae5 jens-daniel-mueller 2021-01-29
b5bdcaf Donghe-Zhu 2021-01-29
442010d Donghe-Zhu 2021-01-29
372adf5 Donghe-Zhu 2021-01-29
af8788e Donghe-Zhu 2021-01-29
21c91c9 Donghe-Zhu 2021-01-29
eded038 Donghe-Zhu 2021-01-29
541d4dd Donghe-Zhu 2021-01-29
6a75576 Donghe-Zhu 2021-01-28
16fba40 Donghe-Zhu 2021-01-28
12bc567 Donghe-Zhu 2021-01-27
ceed31b Donghe-Zhu 2021-01-27
342402d Donghe-Zhu 2021-01-27
5bad5c2 Donghe-Zhu 2021-01-27
61efb56 Donghe-Zhu 2021-01-25
48f638e Donghe-Zhu 2021-01-25
c1cec47 Donghe-Zhu 2021-01-25
05ffb0c Donghe-Zhu 2021-01-25
8b97165 Donghe-Zhu 2021-01-25
c569946 Donghe-Zhu 2021-01-24
a2f0d56 Donghe-Zhu 2021-01-23
28509fc Donghe-Zhu 2021-01-23
4c28e4a Donghe-Zhu 2021-01-22
24cc264 jens-daniel-mueller 2021-01-22
7891955 Donghe-Zhu 2021-01-21
d4cf1cb Donghe-Zhu 2021-01-21
1f3e5b6 jens-daniel-mueller 2021-01-20
0e7bdf1 jens-daniel-mueller 2021-01-15
4571843 jens-daniel-mueller 2021-01-14
b3564aa jens-daniel-mueller 2021-01-14
8d032c3 jens-daniel-mueller 2021-01-14
17dee1d jens-daniel-mueller 2021-01-13
7cdea0c jens-daniel-mueller 2021-01-06
fa85b93 jens-daniel-mueller 2021-01-06
e5cb81a Donghe-Zhu 2021-01-05
a499f10 Donghe-Zhu 2021-01-05
fb8a752 Donghe-Zhu 2020-12-23
8fae0b2 Donghe-Zhu 2020-12-21
c8b76b3 jens-daniel-mueller 2020-12-19 

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.3      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.3   
[21] tidyverse_1.3.0  workflowr_1.6.2 

loaded via a namespace (and not attached):
 [1] fs_1.5.0                 webshot_0.5.2            RColorBrewer_1.1-2      
 [4] httr_1.4.2               rprojroot_2.0.2          tools_4.0.3             
 [7] backports_1.1.10         R6_2.5.0                 nortest_1.0-4           
[10] DBI_1.1.0                colorspace_2.0-0         withr_2.3.0             
[13] gridExtra_2.3            tidyselect_1.1.0         curl_4.3                
[16] compiler_4.0.3           git2r_0.27.1             cli_2.2.0               
[19] rvest_0.3.6              xml2_1.3.2               sass_0.2.0              
[22] labeling_0.4.2           scales_1.1.1             checkmate_2.0.0         
[25] goftest_1.2-2            digest_0.6.27            foreign_0.8-80          
[28] rmarkdown_2.5            rio_0.5.16               pkgconfig_2.0.3         
[31] htmltools_0.5.0          highr_0.8                dbplyr_1.4.4            
[34] rlang_0.4.10             readxl_1.3.1             rstudioapi_0.13         
[37] farver_2.0.3             generics_0.1.0           jsonlite_1.7.2          
[40] zip_2.1.1                car_3.0-10               magrittr_2.0.1          
[43] Matrix_1.2-18            Rcpp_1.0.5               munsell_0.5.0           
[46] fansi_0.4.1              abind_1.4-5              lifecycle_0.2.0         
[49] stringi_1.5.3            whisker_0.4              yaml_2.2.1              
[52] carData_3.0-4            plyr_1.8.6               grid_4.0.3              
[55] blob_1.2.1               parallel_4.0.3           promises_1.1.1          
[58] crayon_1.3.4             lattice_0.20-41          haven_2.3.1             
[61] hms_0.5.3                pillar_1.4.7             reprex_0.3.0            
[64] glue_1.4.2               evaluate_0.14            RcppArmadillo_0.10.1.2.2
[67] data.table_1.13.6        modelr_0.1.8             vctrs_0.3.6             
[70] httpuv_1.5.4             cellranger_1.1.0         gtable_0.3.0            
[73] reshape_0.8.8            assertthat_0.2.1         xfun_0.20               
[76] openxlsx_4.2.3           RcppEigen_0.3.3.9.1      later_1.1.0.1           
[79] viridisLite_0.3.0        ellipsis_0.3.1           here_1.0.1