Run linear model on station and parameters
pr_model_data(df)
A dataframe with the model attribute added
df <- planktonr::pr_get_EOVs(Survey = "NRS") %>%
dplyr::filter(StationCode != "VBM") %>%
dplyr::filter(Parameters == "PigmentChla_mgm3") %>%
pr_remove_outliers(2) %>%
pr_model_data()