## Not run:
# x1 <- readNWISstat(siteNumbers=c("02319394"),
# parameterCd=c("00060"),
# statReportType="annual")
#
# #all the annual mean discharge data for two sites
# x2 <- readNWISstat(siteNumbers=c("02319394","02171500"),
# parameterCd=c("00010","00060"),
# statReportType="annual")
#
# #Request p25, p75, and mean values for temperature and discharge for the 2000s
# #Note that p25 and p75 were not available for temperature, and return NAs
# x <- readNWISstat(siteNumbers=c("02171500"),
# parameterCd=c("00010","00060"),
# statReportType="daily",
# statType=c("mean","median"),
# startDate="2000",endDate="2010")
# ## End(Not run)
Run the code above in your browser using DataLab