pulling hydrograph data and making a NWIS-web style graph

usgs-054279465_hydrograph.r

library(dataRetrieval)
library(dplyr)
library(gsplot)

middleton_q <- readNWISuv(siteNumbers = '054279465', 
                          parameterCd = '00060',
                          startDate = '2016-09-19')
middleton_q_median <- readNWISstat(siteNumbers = '054279465', 
                          parameterCd = '00060',
                          statType = 'median') %>% 
  filter(month_nu==9) %>% 
  filter(day_nu >= 19 & day_nu <= 22) %>% 
  mutate(date = as.POSIXct(paste(paste("2016", month_nu, day_nu, sep="-"), "12:00")))

middleton_attr <- attributes(middleton_q)

middleton_hydrograph <- gsplot(yaxs='r', ylog=TRUE) %>% 
  lines(x=middleton_q$dateTime, y=middleton_q$X_00060_00011, col="blue") %>% 
  points(x=middleton_q_median$date, y=middleton_q_median$p50_va, pch=2, 
         col="orange", legend.name="median") %>% 
  axis(side=2, n.minor=5) %>% 
  grid() %>% 
  legend(location="below") %>% 
  title(main=middleton_attr$siteInfo$siteName,
        ylab=middleton_attr$variableInfo$variableDescription)
middleton_hydrograph