library(dplyr)
# data to model
tomod <- rawdat %>%
filter(station %in% 34) %>%
filter(param %in% 'chl') %>%
filter(yr > 2015)
mod <- anlz_gam(tomod, trans = 'log10')
show_prd3d(mod, ylab = 'Chlorophyll-a (ug/L)')
Run the code above in your browser using DataLab