In jsta/nlaR: Interface to the National Lakes Assessment

# rmarkdown::render("vignettes/summary_2012.Rmd", rmarkdown::html_document(pandoc_args = "--self-contained"))

library(dplyr)
library(magrittr)
library(purrr)

library(nlaR)

dt <- nla_load(2012)

tp <- select(dt$waterchem_wide, PTL_RESULT, UID) %>%
      left_join(select(dt$wide_siteinfo, AREA_HA, SITE_ID, UID)) %>%
      group_by(SITE_ID) %>%
      summarize(tp = median(PTL_RESULT), area = median(AREA_HA) / 100)

chl  <- select(dt$chla_wide, CHLX_RESULT, UID)  %>%
        left_join(select(dt$wide_siteinfo, SITE_ID, UID)) %>%
        group_by(SITE_ID) %>%
        summarize(chl = median(CHLX_RESULT))

secchi  <- select(dt$secchi, SECCHI, SITE_ID)  %>%
           group_by(SITE_ID) %>%
           summarize(secchi = median(SECCHI))

res <- reduce(list(tp, chl, secchi), left_join)
res <- left_join(res, select(dt$wide_phab, depth = DEPTH_AT_STATION, SITE_ID))


skimr::skim(res)

jsta/nlaR documentation built on March 23, 2023, 11:53 p.m.