inst/doc/intro.R

In USGS-R/EflowStats: Hydrologic Indicator and Alteration Stats

## ----message=FALSE, echo=TRUE, eval=FALSE--------------------------------
#  library(dataRetrieval)
#  library(EflowStats)
#  
#  #Get some data
#  dailyQ <- readNWISdv(siteNumber = "04085427",
#                       parameterCd = "00060",
#                       startDate = "2000-10-01",
#                       endDate = "2012-9-30")
#  #Check data for completeness
#  dailyQClean <- validate_data(dailyQ[c("Date","X_00060_00003")],yearType="water")
#  
#  #Get drainage area
#  siteInfo <- readNWISsite(siteNumber = "04085427")
#  drainageArea <- siteInfo$drain_area_va
#  
#  #Get peak flows
#  peakFlows <- readNWISpeak(siteNumber = "04085427",
#                            startDate = "2000-10-01",
#                            endDate = "2012-9-30")
#  

## ----message=FALSE, echo=TRUE, eval=FALSE--------------------------------
#  
#  #Check data for completeness
#  dailyQClean <- dataCheck(dailyQ[c("Date", "X_00060_00003")], yearType="water")
#  
#  #Get flood recurence threshold
#  
#  floodThresh <- get_peakThreshold(dailyQClean[c("date","discharge")],
#  peakFlows[c("peak_dt","peak_va")])
#  
#  #Calculate all hit stats
#  calc_allHITOut <- calc_allHIT(dailyQClean,
#  drainArea=drainageArea,
#  floodThreshold=floodThresh)
#  
#  #Calculate mag7 stats
#  magnifStatsOut <- calc_magnifSeven(dailyQClean,yearType="water",digits=3)
#  
#