diff --git a/R/readNWISdata.r b/R/readNWISdata.r
index c66e7f7f53375a93c5c36fd3f55ef5599fe6f0b2..9b3e5a76c53d4fdf4cb729191fb454ffd58ae8b9 100644
--- a/R/readNWISdata.r
+++ b/R/readNWISdata.r
@@ -46,11 +46,13 @@
 #' #Not empty:
 #' multiSite <- readNWISdata(sites=c("04025500","040263491"), service="iv", parameterCd="00060")
 #' bBoxEx <- readNWISdata(bBox=c(-83,36.5,-81,38.5), parameterCd="00010")
+#' 
 #' startDate <- as.Date("2013-10-01")
 #' endDate <- as.Date("2014-09-30")
 #' waterYear <- readNWISdata(bBox=c(-83,36.5,-81,38.5), parameterCd="00010", 
 #'                   service="dv", startDate=startDate, endDate=endDate)
-#' siteInfo <- readNWISdata(stateCd="WI", parameterCd="00010",hasDataTypeCd="iv", service="site")
+#' siteInfo <- readNWISdata(stateCd="WI", parameterCd="00010",
+#'                   hasDataTypeCd="iv", service="site")
 #' }
 readNWISdata <- function(service="dv", ...){
   
@@ -94,9 +96,7 @@ readNWISdata <- function(service="dv", ...){
   } else {
     retval <- importWaterML1(urlCall, asDateTime = ("iv" == service))
     if("dv" == service){
-      retval$dateTime <- as.Date(retval$dateTime)
-      retval$tz_cd <- NULL      
-      names(retval)[names(retval) == "dateTime"] <- "Date"
+      retval$dateTime <- as.POSIXct(retval$dateTime)
     }
   }
   
diff --git a/R/readNWISqw.r b/R/readNWISqw.r
index 5f3679e411a0d907dcfb9277a59381442eb509ad..81c210d534294b9278d5b84e0f8cf09deb5596de 100644
--- a/R/readNWISqw.r
+++ b/R/readNWISqw.r
@@ -96,9 +96,6 @@ readNWISqw <- function (siteNumbers,parameterCd,startDate="",endDate="",
   
   pgrp <- c("INF", "PHY", "INM", "INN", "NUT", "MBI", "BIO", "IMM", "IMN", "TOX",
                            "OPE", "OPC", "OOT", "RAD", "XXX", "SED", "POP")
-  
-  
-  
 
   if(any(parameterCd == "all") | any(parameterCd == "All") ){
     
@@ -128,8 +125,6 @@ readNWISqw <- function (siteNumbers,parameterCd,startDate="",endDate="",
                             endDate,"qw",expanded=expanded)    
   }
 
-
-  
   data <- importRDB1(url,asDateTime=TRUE, qw=TRUE, tz = tz)
   originalHeader <- comment(data)
   parameterCd <- unique(data$parm_cd)
diff --git a/R/whatNWISsites.R b/R/whatNWISsites.R
index 82c847f4684a33816cd5ac41f794a365e96b057b..b271d21270169cfc800c9014a95c7f587728120b 100644
--- a/R/whatNWISsites.R
+++ b/R/whatNWISsites.R
@@ -33,11 +33,12 @@
 whatNWISsites <- function(...){
   
   matchReturn <- list(...)
-
-  values <- sapply(matchReturn, function(x) URLencode(as.character(paste(eval(x),collapse="",sep=""))))
+  values <- sapply(matchReturn, function(x) URLencode(as.character(paste(eval(x),collapse=",",sep=""))))
   
   urlCall <- paste(paste(names(values),values,sep="="),collapse="&")
   
+  names(values)[names(values) == "siteNumber"] <- "sites"
+  names(values)[names(values) == "siteNumbers"] <- "sites"
   
   baseURL <- "http://waterservices.usgs.gov/nwis/site/?format=mapper&"
   urlCall <- paste(baseURL,
diff --git a/inst/doc/dataRetrieval.R b/inst/doc/dataRetrieval.R
index 488c76630dd3fe1ec8de8c7d0f83608860878007..8e96c27ce3807879c04e1173b51fa20a088f1e58 100644
--- a/inst/doc/dataRetrieval.R
+++ b/inst/doc/dataRetrieval.R
@@ -73,83 +73,82 @@ print(xtable(data.df,label="tab:stat",
       )
 
 
-## ----getSite, echo=TRUE-----------------------------------
-siteNumbers <- c("01491000","01645000") 
-siteINFO <- readNWISsite(siteNumbers)
-
-## ----siteNames2, echo=TRUE--------------------------------
-siteINFO$station_nm
+## ----getSite, echo=TRUE, eval=FALSE-----------------------
+#  siteNumbers <- c("01491000","01645000")
+#  siteINFO <- readNWISsite(siteNumbers)
 
 ## ----siteNames3, echo=TRUE, eval=FALSE--------------------
 #  comment(siteINFO)
 
-## ----getSiteExtended, echo=TRUE---------------------------
-# Continuing from the previous example:
-# This pulls out just the daily, mean data:
-
-dailyDataAvailable <- whatNWISdata(siteNumbers,
-                    service="dv", statCd="00003")
-
-
-
-## ----tablegda, echo=FALSE,results='asis'------------------
-tableData <- with(dailyDataAvailable, 
-      data.frame( 
-      siteNumber= site_no,
-      srsname=srsname, 
-      startDate=as.character(begin_date), 
-      endDate=as.character(end_date), 
-      count=as.character(count_nu),
-      units=parameter_units,
-#       statCd = stat_cd,
-      stringsAsFactors=FALSE)
-      )
-
-tableData$units[which(tableData$units == "ft3/s")] <- "ft$^3$/s"
-tableData$units[which(tableData$units == "uS/cm @25C")] <- "$\\mu$S/cm @25C"
-
-
-print(xtable(tableData,label="tab:gda",
-    caption="Reformatted version of output from \\texttt{whatNWISdata} function for the Choptank River near Greensboro, MD, and from Seneca Creek at Dawsonville, MD from the daily values service [Some columns deleted for space considerations]"),
-       caption.placement="top",
-       size = "\\footnotesize",
-       latex.environment=NULL,
-       sanitize.text.function = function(x) {x},
-       sanitize.colnames.function =  bold.colHeaders,
-       sanitize.rownames.function = addSpace
-      )
-
-
-## ----label=getPCodeInfo, echo=TRUE------------------------
-# Using defaults:
-parameterCd <- "00618" 
-parameterINFO <- readNWISpCode(parameterCd)
-colnames(parameterINFO)
-
-## ----siteNames, echo=TRUE---------------------------------
-parameterINFO$parameter_nm
-
-## ----label=getNWISDaily, echo=TRUE, eval=TRUE-------------
+## ----getSiteExtended, echo=TRUE, eval=FALSE---------------
+#  # Continuing from the previous example:
+#  # This pulls out just the daily, mean data:
+#  
+#  dailyDataAvailable <- whatNWISdata(siteNumbers,
+#                      service="dv", statCd="00003")
+#  
+#  
 
-# Choptank River near Greensboro, MD:
-siteNumber <- "01491000"
-parameterCd <- "00060"  # Discharge
-startDate <- "2009-10-01"  
-endDate <- "2012-09-30" 
+## ----tablegda, echo=FALSE,eval=FALSE----------------------
+#  tableData <- with(dailyDataAvailable,
+#        data.frame(
+#        siteNumber= site_no,
+#        srsname=srsname,
+#        startDate=as.character(begin_date),
+#        endDate=as.character(end_date),
+#        count=as.character(count_nu),
+#        units=parameter_units,
+#  #       statCd = stat_cd,
+#        stringsAsFactors=FALSE)
+#        )
+#  
+#  tableData$units[which(tableData$units == "ft3/s")] <- "ft$^3$/s"
+#  tableData$units[which(tableData$units == "uS/cm @25C")] <- "$\\mu$S/cm @25C"
+#  
+#  
+#  print(xtable(tableData,label="tab:gda",
+#      caption="Reformatted version of output from \\texttt{whatNWISdata} function for the Choptank River near Greensboro, MD, and from Seneca Creek at Dawsonville, MD from the daily values service [Some columns deleted for space considerations]"),
+#         caption.placement="top",
+#         size = "\\footnotesize",
+#         latex.environment=NULL,
+#         sanitize.text.function = function(x) {x},
+#         sanitize.colnames.function =  bold.colHeaders,
+#         sanitize.rownames.function = addSpace
+#        )
+#  
 
-discharge <- readNWISdv(siteNumber, 
-                    parameterCd, startDate, endDate)
-names(discharge)
+## ----label=getPCodeInfo, echo=TRUE, eval=FALSE------------
+#  # Using defaults:
+#  parameterCd <- "00618"
+#  parameterINFO <- readNWISpCode(parameterCd)
 
-## ----label=getNWIStemperature, echo=TRUE------------------
+## ----label=getNWISDaily, echo=TRUE, eval=FALSE------------
+#  
+#  # Choptank River near Greensboro, MD:
+#  siteNumber <- "01491000"
+#  parameterCd <- "00060"  # Discharge
+#  startDate <- "2009-10-01"
+#  endDate <- "2012-09-30"
+#  
+#  discharge <- readNWISdv(siteNumber,
+#                      parameterCd, startDate, endDate)
 
-parameterCd <- c("00010","00060")  # Temperature and discharge
-statCd <- c("00001","00003")  # Mean and maximum
-startDate <- "2012-01-01"
-endDate <- "2012-05-01"
+## ----label=getNWIStemperature, echo=TRUE, eval=FALSE------
+#  siteNumber <- "01491000"
+#  parameterCd <- c("00010","00060")  # Temperature and discharge
+#  statCd <- c("00001","00003")  # Mean and maximum
+#  startDate <- "2012-01-01"
+#  endDate <- "2012-05-01"
+#  
+#  temperatureAndFlow <- readNWISdv(siteNumber, parameterCd,
+#          startDate, endDate, statCd=statCd)
+#  
 
-temperatureAndFlow <- readNWISdv(siteNumber, parameterCd, 
-        startDate, endDate, statCd=statCd)
+## ----label=getNWIStemperature2, echo=FALSE, eval=TRUE-----
+filePath <- system.file("extdata", package="dataRetrieval")
+fileName <- "temperatureAndFlow.RData"
+fullPath <- file.path(filePath, fileName)
+load(fullPath)
 
 
 ## ----label=renameColumns, echo=TRUE-----------------------
@@ -186,77 +185,61 @@ title(paste(siteInfo$station_nm,"2012"))
 legend("topleft", variableInfo$param_units, 
        col=c("black","red"),lty=c(NA,1),pch=c(1,NA))
 
-## ----label=readNWISuv, echo=TRUE--------------------------
-
-parameterCd <- "00060"  # Discharge
-startDate <- "2012-05-12" 
-endDate <- "2012-05-13" 
-dischargeUnit <- readNWISuv(siteNumber, parameterCd, 
-        startDate, endDate)
-dischargeUnit <- renameNWISColumns(dischargeUnit)
-
-## ----dischargeData, echo=TRUE-----------------------------
-head(dischargeUnit)
-
-## ----label=getQW, echo=TRUE, eval=TRUE--------------------
- 
-# Dissolved Nitrate parameter codes:
-parameterCd <- c("00618","71851")
-startDate <- "1985-10-01"
-endDate <- "2012-09-30"
-
-dfLong <- readNWISqw(siteNumber, parameterCd, 
-      startDate, endDate)
-
-# Or the wide return:
-# dfWide <- readNWISqw(siteNumber, parameterCd, 
-#       startDate, endDate, reshape=TRUE)
-
+## ----label=readNWISuv, eval=FALSE-------------------------
+#  
+#  parameterCd <- "00060"  # Discharge
+#  startDate <- "2012-05-12"
+#  endDate <- "2012-05-13"
+#  dischargeUnit <- readNWISuv(siteNumber, parameterCd,
+#          startDate, endDate)
+#  dischargeUnit <- renameNWISColumns(dischargeUnit)
+
+## ----label=getQW, echo=TRUE, eval=FALSE-------------------
+#  
+#  # Dissolved Nitrate parameter codes:
+#  parameterCd <- c("00618","71851")
+#  startDate <- "1985-10-01"
+#  endDate <- "2012-09-30"
+#  
+#  dfLong <- readNWISqw(siteNumber, parameterCd,
+#        startDate, endDate)
+#  
+#  # Or the wide return:
+#  # dfWide <- readNWISqw(siteNumber, parameterCd,
+#  #       startDate, endDate, reshape=TRUE)
+#  
 
 ## ----qwmeta, echo=TRUE, eval=FALSE------------------------
 #  
 #  comment(dfLong)
 #  
 
-## ----gwlexample, echo=TRUE, eval=TRUE---------------------
-siteNumber <- "434400121275801"
-groundWater <- readNWISgwl(siteNumber)
-
-names(groundWater)
-
+## ----gwlexample, echo=TRUE, eval=FALSE--------------------
+#  siteNumber <- "434400121275801"
+#  groundWater <- readNWISgwl(siteNumber)
 
-## ----peakexample, echo=TRUE, eval=TRUE--------------------
-siteNumber <- '01594440'
-peakData <- readNWISpeak(siteNumber)
-
-
-names(peakData)
-
-
-## ----ratingexample, echo=TRUE, eval=TRUE------------------
-ratingData <- readNWISrating(siteNumber, "base")
-attr(ratingData, "RATING")
-
-names(ratingData)
-
-
-## ----surfexample, echo=TRUE, eval=TRUE--------------------
-surfaceData <- readNWISmeas(siteNumber)
+## ----peakexample, echo=TRUE, eval=FALSE-------------------
+#  siteNumber <- '01594440'
+#  peakData <- readNWISpeak(siteNumber)
+#  
 
-names(surfaceData)
+## ----ratingexample, echo=TRUE, eval=FALSE-----------------
+#  ratingData <- readNWISrating(siteNumber, "base")
+#  attr(ratingData, "RATING")
+#  
 
+## ----surfexample, echo=TRUE, eval=FALSE-------------------
+#  surfaceData <- readNWISmeas(siteNumber)
+#  
 
 ## ----label=getQWData, echo=TRUE, eval=FALSE---------------
 #  specificCond <- readWQPqw('WIDNR_WQX-10032762',
 #                  'Specific conductance','2011-05-01','2011-09-30')
 
-## ----siteSearch-------------------------------------------
-sites <- whatNWISsites(bBox="-83.0,36.5,-81.0,38.5", 
-                      parameterCd="00010,00060",
-                      hasDataTypeCd="dv")
-
-names(sites)
-nrow(sites)
+## ----siteSearch, eval=FALSE-------------------------------
+#  sites <- whatNWISsites(bBox=c(-83.0,36.5,-81.0,38.5),
+#                        parameterCd=c("00010","00060"),
+#                        hasDataTypeCd="dv")
 
 ## ----dataExample------------------------------------------
 dischargeWI <- readNWISdata(service="dv",
@@ -306,16 +289,12 @@ variableInfo <- attr(dischargeWI, "variableInfo")
 
 
 
-## ----meta5, eval=TRUE, eval=FALSE-------------------------
+## ----meta5, eval=FALSE------------------------------------
 #  comment(peakData)
 #  
 #  #Which is equivalent to:
 #  # attr(peakData, "comment")
 
-## ----meta6, eval=TRUE, eval=TRUE--------------------------
-comment(peakData)[c(1:15,58:66)]
-
-
 ## ----helpFunc,eval = FALSE--------------------------------
 #  ?readNWISpCode
 
@@ -328,19 +307,18 @@ comment(peakData)[c(1:15,58:66)]
 ## ----openLibraryTest, eval=FALSE--------------------------
 #  library(dataRetrieval)
 
-## ----label=getSiteApp, echo=TRUE--------------------------
-availableData <- whatNWISdata(siteNumber, "dv")
-dailyData <- availableData["00003" == availableData$stat_cd,]
-
-tableData <- with(dailyData, 
-      data.frame(
-        shortName=srsname, 
-        Start=begin_date, 
-        End=end_date, 
-        Count=count_nu,
-        Units=parameter_units)
-      )
-tableData
+## ----label=getSiteApp, echo=TRUE, eval=FALSE--------------
+#  availableData <- whatNWISdata(siteNumber, "dv")
+#  dailyData <- availableData["00003" == availableData$stat_cd,]
+#  
+#  tableData <- with(dailyData,
+#        data.frame(
+#          shortName=srsname,
+#          Start=begin_date,
+#          End=end_date,
+#          Count=count_nu,
+#          Units=parameter_units)
+#        )
 
 ## ----label=saveData, echo=TRUE, eval=FALSE----------------
 #  write.table(tableData, file="tableData.tsv",sep="\t",
diff --git a/inst/doc/dataRetrieval.Rnw b/inst/doc/dataRetrieval.Rnw
index 40390f31350c030097b6860cf987b404139a1cdb..bcdf370041f96531fb2e995897892266f7ad4f04 100644
--- a/inst/doc/dataRetrieval.Rnw
+++ b/inst/doc/dataRetrieval.Rnw
@@ -354,17 +354,13 @@ Examples for using these siteNumbers, parameter codes, and stat codes will be pr
 Use the \texttt{readNWISsite} function to obtain all of the information available for a particular USGS site (or sites) such as full station name, drainage area, latitude, and longitude. \texttt{readNWISsite} can also access information about multiple sites with a vector input.
 
 
-<<getSite, echo=TRUE>>=
+<<getSite, echo=TRUE, eval=FALSE>>=
 siteNumbers <- c("01491000","01645000") 
 siteINFO <- readNWISsite(siteNumbers)
 @
 
-A specific example piece of information can be retrieved, in this case a station name, as follows:
-
-<<siteNames2, echo=TRUE>>=
-siteINFO$station_nm
-@
-Site information is obtained from \url{http://waterservices.usgs.gov/rest/Site-Test-Tool.html}
+Site information is obtained from:
+\url{http://waterservices.usgs.gov/rest/Site-Test-Tool.html}
 
 Information on the returned data can be found with the \texttt{comment} function as described in section \ref{sec:metadata}.
 
@@ -385,7 +381,7 @@ To discover what data is available for a particular USGS site, including measure
 In the following example, we limit the retrieved data to only daily data. The default for \texttt{"}service\texttt{"} is \enquote{all}, which returns all of the available data for that site. Likewise, there are arguments for parameter code (\texttt{parameterCd}) and statistic code (\texttt{statCd}) to filter the results. The default for both is to return all possible values (\enquote{all}). The returned \texttt{"}count\_nu\texttt{"} for \texttt{"}uv\texttt{"} data is the count of days with returned data, not the actual count of returned values.
 
 
-<<getSiteExtended, echo=TRUE>>=
+<<getSiteExtended, echo=TRUE, eval=FALSE>>=
 # Continuing from the previous example:
 # This pulls out just the daily, mean data:
 
@@ -395,7 +391,7 @@ dailyDataAvailable <- whatNWISdata(siteNumbers,
 
 @
 
-<<tablegda, echo=FALSE,results='asis'>>=
+<<tablegda, echo=FALSE,eval=FALSE>>=
 tableData <- with(dailyDataAvailable, 
       data.frame( 
       siteNumber= site_no,
@@ -424,7 +420,24 @@ print(xtable(tableData,label="tab:gda",
 
 @
 
-
+\begin{table}[ht]
+\caption{Reformatted version of output from \texttt{whatNWISdata} function for the Choptank River near Greensboro, MD, and from Seneca Creek at Dawsonville, MD from the daily values service [Some columns deleted for space considerations]} 
+\label{tab:gda}
+{\footnotesize
+\begin{tabular}{rllllll}
+  \hline
+ & \multicolumn{1}{c}{\textbf{\textsf{siteNumber}}} & \multicolumn{1}{c}{\textbf{\textsf{srsname}}} & \multicolumn{1}{c}{\textbf{\textsf{startDate}}} & \multicolumn{1}{c}{\textbf{\textsf{endDate}}} & \multicolumn{1}{c}{\textbf{\textsf{count}}} & \multicolumn{1}{c}{\textbf{\textsf{units}}} \\ 
+  \hline
+ & 01491000 & Temperature, water & 2010-10-01 & 2012-05-09 & 529 & deg C \\ 
+  [5pt] & 01645000 & Stream flow, mean. daily & 1930-09-26 & 2015-02-19 & 30828 & ft$^3$/s \\ 
+  [5pt] & 01491000 & Stream flow, mean. daily & 1948-01-01 & 2015-02-19 & 24522 & ft$^3$/s \\ 
+  [5pt] & 01491000 & Specific conductance & 2010-10-01 & 2012-05-09 & 527 & $\mu$S/cm @25C \\ 
+  [5pt] & 01491000 & Suspended sediment concentration (SSC) & 1980-10-01 & 1991-09-30 & 3651 & mg/l \\ 
+  [5pt] & 01491000 & Suspended sediment discharge & 1980-10-01 & 1991-09-30 & 3652 & tons/day \\ 
+   \hline
+\end{tabular}
+}
+\end{table}
 
 See Section \ref{app:createWordTable} for instructions on converting an R data frame to a table in Microsoft\textregistered\ software Excel or Word to display a data availability table similar to Table \ref{tab:gda}. Excel, Microsoft, PowerPoint, Windows, and Word are registered trademarks of Microsoft Corporation in the United States and other countries.
 
@@ -436,19 +449,13 @@ See Section \ref{app:createWordTable} for instructions on converting an R data f
 %------------------------------------------------------------
 To obtain all of the available information concerning a measured parameter (or multiple parameters), use the \texttt{readNWISpCode} function:
 
-<<label=getPCodeInfo, echo=TRUE>>=
+<<label=getPCodeInfo, echo=TRUE, eval=FALSE>>=
 # Using defaults:
 parameterCd <- "00618" 
 parameterINFO <- readNWISpCode(parameterCd)
-colnames(parameterINFO)
 @
 
-A specific example piece of information, in this case parameter name, can be obtained as follows:
-
-<<siteNames, echo=TRUE>>=
-parameterINFO$parameter_nm
-@
-Parameter information can obtained from:
+More information on parameter codes can obtained from:
 
 \url{http://help.waterdata.usgs.gov/codes-and-parameters/parameters}
 
@@ -462,7 +469,7 @@ To obtain daily records of USGS data, use the \texttt{readNWISdv} function. The
 
 The dates (start and end) must be in the format \texttt{"}YYYY-MM-DD\texttt{"} (note: the user must include the quotes).  Setting the start date to \texttt{"}\texttt{"} (no space) will prompt the program to ask for the earliest date, and setting the end date to \texttt{"}\texttt{"} (no space) will prompt for the latest available date.
 
-<<label=getNWISDaily, echo=TRUE, eval=TRUE>>=
+<<label=getNWISDaily, echo=TRUE, eval=FALSE>>=
 
 # Choptank River near Greensboro, MD:
 siteNumber <- "01491000"
@@ -472,7 +479,6 @@ endDate <- "2012-09-30"
 
 discharge <- readNWISdv(siteNumber, 
                     parameterCd, startDate, endDate)
-names(discharge)
 @
 
 The column \texttt{"}datetime\texttt{"} in the returned data frame is automatically imported as a variable of class \texttt{"}Date\texttt{"} in R. Each requested parameter has a value and remark code column.  The names of these columns depend on the requested parameter and stat code combinations. USGS daily value qualification codes are often \texttt{"}A\texttt{"} (approved for publication) or \texttt{"}P\texttt{"} (provisional data subject to revision). A more complete list of daily value qualification codes can be found here:
@@ -481,8 +487,8 @@ The column \texttt{"}datetime\texttt{"} in the returned data frame is automatica
 
 Another example would be a request for mean and maximum daily temperature and discharge in early 2012:
 
-<<label=getNWIStemperature, echo=TRUE>>=
-
+<<label=getNWIStemperature, echo=TRUE, eval=FALSE>>=
+siteNumber <- "01491000"
 parameterCd <- c("00010","00060")  # Temperature and discharge
 statCd <- c("00001","00003")  # Mean and maximum
 startDate <- "2012-01-01"
@@ -493,6 +499,13 @@ temperatureAndFlow <- readNWISdv(siteNumber, parameterCd,
 
 @
 
+<<label=getNWIStemperature2, echo=FALSE, eval=TRUE>>=
+filePath <- system.file("extdata", package="dataRetrieval")
+fileName <- "temperatureAndFlow.RData"
+fullPath <- file.path(filePath, fileName)
+load(fullPath)
+
+@
 
 The column names can be shortened and simplified using the \texttt{renameNWISColumns} function.  This is not necessary, but may streamline subsequent data analysis and presentation. Site information, daily statistic information, and measured parameter information is attached to the data frame as attributes. This is discused further in section \ref{sec:metadata}.
 
@@ -549,7 +562,7 @@ There are occasions where NWIS values are not reported as numbers, instead there
 %------------------------------------------------------------
 Any data collected at regular time intervals (such as 15-minute or hourly) are known as \enquote{unit values.} Many of these are delivered on a real time basis and very recent data (even less than an hour old in many cases) are available through the function \texttt{readNWISuv}.  Some of these unit values are available for many years, and some are only available for a recent time period such as 120 days.  Here is an example of a retrieval of such data.  
 
-<<label=readNWISuv, echo=TRUE>>=
+<<label=readNWISuv, eval=FALSE>>=
 
 parameterCd <- "00060"  # Discharge
 startDate <- "2012-05-12" 
@@ -561,11 +574,6 @@ dischargeUnit <- renameNWISColumns(dischargeUnit)
 
 The retrieval produces a data frame that contains 96 rows (one for every 15 minute period in the day).  They include all data collected from the startDate through the endDate (starting and ending with midnight locally-collected time). The dateTime column is converted to \enquote{UTC} (Coordinated Universal Time), so midnight EST will be 5 hours earlier in the dateTime column (the previous day, at 7pm).
 
-
-<<dischargeData, echo=TRUE>>=
-head(dischargeUnit)
-@
-
 To override the UTC timezone, specify a valid timezone in the tz argument. Default is \texttt{""}, which will keep the dateTime column in UTC. Other valid timezones are:
 
 \begin{verbatim}
@@ -594,7 +602,7 @@ Data are retrieved from \url{http://waterservices.usgs.gov/rest/IV-Test-Tool.htm
 %------------------------------------------------------------
 To get USGS water quality data from water samples collected at the streamgage or other monitoring site (as distinct from unit values collected through some type of automatic monitor) we can use the function \texttt{readNWISqw}, with the input arguments: siteNumber, parameterCd, startDate, and endDate. Additionally, the argument \texttt{"}expanded\texttt{"} is a logical input that allows the user to choose between a simple return of datetimes/qualifier/values (expanded=FALSE), or a more complete and verbose output (expanded=TRUE). Expanded = TRUE includes such columns as remark codes, value qualifying text, and detection level for each parameter code. There also includes an argument \texttt{"}reshape\texttt{"}, that converts the expanded dataset to a \texttt{"}wide\texttt{"} format (each requested parameter code gets individual columns). The defaults are expanded=TRUE, and reshape=FALSE. 
 
-<<label=getQW, echo=TRUE, eval=TRUE>>=
+<<label=getQW, echo=TRUE, eval=FALSE>>=
  
 # Dissolved Nitrate parameter codes:
 parameterCd <- c("00618","71851")
@@ -626,12 +634,9 @@ comment(dfLong)
 %------------------------------------------------------------
 Groundwater level measurements can be obtained with the \texttt{readNWISgwl} function. Information on the returned data can be found with the \texttt{comment} function, and attached attributes as described in section \ref{sec:metadata}.
 
-<<gwlexample, echo=TRUE, eval=TRUE>>=
+<<gwlexample, echo=TRUE, eval=FALSE>>=
 siteNumber <- "434400121275801"
 groundWater <- readNWISgwl(siteNumber)
-
-names(groundWater)
-
 @
 
 %------------------------------------------------------------
@@ -641,13 +646,10 @@ names(groundWater)
 
 Peak flow data are instantaneous discharge or stage data that record the maximum values of these variables during a flood event.  They include the annual peak flood event but can also include records of other peaks that are lower than the annual maximum. Peak discharge measurements can be obtained with the \texttt{readNWISpeak} function. Information on the returned data can be found with the \texttt{comment} function and attached attributes as described in section \ref{sec:metadata}.
 
-<<peakexample, echo=TRUE, eval=TRUE>>=
+<<peakexample, echo=TRUE, eval=FALSE>>=
 siteNumber <- '01594440'
 peakData <- readNWISpeak(siteNumber)
 
-
-names(peakData)
-
 @
 
 
@@ -659,12 +661,10 @@ Rating curves are the calibration curves that are used to convert measurements o
 
 Rating curves can be obtained with the \texttt{readNWISrating} function.
 
-<<ratingexample, echo=TRUE, eval=TRUE>>=
+<<ratingexample, echo=TRUE, eval=FALSE>>=
 ratingData <- readNWISrating(siteNumber, "base")
 attr(ratingData, "RATING")
 
-names(ratingData)
-
 @
 
 
@@ -677,11 +677,9 @@ These data are the discrete measurements of discharge that are made for the purp
 
 Surface-water measurement data can be obtained with the \texttt{readNWISmeas} function.
 
-<<surfexample, echo=TRUE, eval=TRUE>>=
+<<surfexample, echo=TRUE, eval=FALSE>>=
 surfaceData <- readNWISmeas(siteNumber)
 
-names(surfaceData)
-
 @
 
 
@@ -697,7 +695,9 @@ specificCond <- readWQPqw('WIDNR_WQX-10032762',
                 'Specific conductance','2011-05-01','2011-09-30')
 @
 
-A tool for finding NWIS characteristic names can be found at: \url{http://www.waterqualitydata.us/public_srsnames.jsp}. 
+A tool for finding NWIS characteristic names can be found at: 
+
+\url{http://www.waterqualitydata.us/public_srsnames.jsp}
 
 \FloatBarrier
 
@@ -721,13 +721,10 @@ to discover many options for searching for NWIS sites. For example, you may want
 
 The following dataRetrieval code can be used to get those sites:
 
-<<siteSearch>>=
-sites <- whatNWISsites(bBox="-83.0,36.5,-81.0,38.5", 
-                      parameterCd="00010,00060",
+<<siteSearch, eval=FALSE>>=
+sites <- whatNWISsites(bBox=c(-83.0,36.5,-81.0,38.5), 
+                      parameterCd=c("00010","00060"),
                       hasDataTypeCd="dv")
-
-names(sites)
-nrow(sites)
 @
 
 
@@ -850,19 +847,13 @@ variableInfo <- attr(dischargeWI, "variableInfo")
 
 Data obtained from \texttt{readNWISpeak}, \texttt{readNWISmeas}, and \texttt{readNWISrating}, the \texttt{comment} attribute is useful.
 
-<<meta5, eval=TRUE, eval=FALSE>>=
+<<meta5, eval=FALSE>>=
 comment(peakData)
 
 #Which is equivalent to:
 # attr(peakData, "comment")
 @
 
-A subset (due to space considerations) of the \texttt{comment} metatdata is shown here:
-
-<<meta6, eval=TRUE, eval=TRUE>>=
-comment(peakData)[c(1:15,58:66)]
-
-@
 
 
 %------------------------------------------------------------ 
@@ -921,7 +912,7 @@ library(dataRetrieval)
 %------------------------------------------------------------
 There are a few steps that are required in order to create a table in Microsoft\textregistered\ software (Excel, Word, PowerPoint, etc.) from an R data frame. There are certainly a variety of good methods, one of which is detailed here. The example we will step through here will be to create a table in Microsoft Excel based on the data frame tableData:
 
-<<label=getSiteApp, echo=TRUE>>=
+<<label=getSiteApp, echo=TRUE, eval=FALSE>>=
 availableData <- whatNWISdata(siteNumber, "dv")
 dailyData <- availableData["00003" == availableData$stat_cd,]
 
@@ -933,7 +924,6 @@ tableData <- with(dailyData,
         Count=count_nu,
         Units=parameter_units)
       )
-tableData
 @
 
 First, save the data frame as a tab delimited file (you don't want to use comma delimited because there are commas in some of the data elements):
diff --git a/inst/doc/dataRetrieval.pdf b/inst/doc/dataRetrieval.pdf
index ab667863f3c7b1d7dc1461293661f11ca2655827..74f8cb0eb3d64fe70e4996df3121cd654a39a3e1 100644
Binary files a/inst/doc/dataRetrieval.pdf and b/inst/doc/dataRetrieval.pdf differ
diff --git a/inst/extdata/temperatureAndFlow.RData b/inst/extdata/temperatureAndFlow.RData
new file mode 100644
index 0000000000000000000000000000000000000000..80062797ae53223fb15f7152d3d1f663dd465bec
Binary files /dev/null and b/inst/extdata/temperatureAndFlow.RData differ
diff --git a/man/readNWISdata.Rd b/man/readNWISdata.Rd
index 0973fbecfb02319caf1232fa5e565576c1bac84e..b9977bc462ecf6810d05e508d951ed6275c97875 100644
--- a/man/readNWISdata.Rd
+++ b/man/readNWISdata.Rd
@@ -54,11 +54,13 @@ multiSite <- readNWISdata(sites=c("04025000","04072150"), service="iv", paramete
 #Not empty:
 multiSite <- readNWISdata(sites=c("04025500","040263491"), service="iv", parameterCd="00060")
 bBoxEx <- readNWISdata(bBox=c(-83,36.5,-81,38.5), parameterCd="00010")
+
 startDate <- as.Date("2013-10-01")
 endDate <- as.Date("2014-09-30")
 waterYear <- readNWISdata(bBox=c(-83,36.5,-81,38.5), parameterCd="00010",
                   service="dv", startDate=startDate, endDate=endDate)
-siteInfo <- readNWISdata(stateCd="WI", parameterCd="00010",hasDataTypeCd="iv", service="site")
+siteInfo <- readNWISdata(stateCd="WI", parameterCd="00010",
+                  hasDataTypeCd="iv", service="site")
 }
 }
 \seealso{
diff --git a/tests/testthat/tests_general.R b/tests/testthat/tests_general.R
index 979d4aeac62c45257d1cbcb904f6c94986ee4f60..c714f9ba24f5842a06127d5395c419804256af71 100644
--- a/tests/testthat/tests_general.R
+++ b/tests/testthat/tests_general.R
@@ -14,7 +14,7 @@ test_that("General NWIS retrievals working", {
   endDate <- as.Date("2014-09-30")
   waterYear <- readNWISdata(bBox=c(-83,36.5,-81,38.5), parameterCd="00010", 
                    service="dv", startDate=startDate, endDate=endDate)
-  expect_is(waterYear$Date, 'Date')
+  expect_is(waterYear$dateTime, 'POSIXct')
   
   siteInfo <- readNWISdata(stateCd="WI", parameterCd="00010",hasDataTypeCd="iv", 
                            service="site")
diff --git a/vignettes/dataRetrieval.Rnw b/vignettes/dataRetrieval.Rnw
index 40390f31350c030097b6860cf987b404139a1cdb..bcdf370041f96531fb2e995897892266f7ad4f04 100644
--- a/vignettes/dataRetrieval.Rnw
+++ b/vignettes/dataRetrieval.Rnw
@@ -354,17 +354,13 @@ Examples for using these siteNumbers, parameter codes, and stat codes will be pr
 Use the \texttt{readNWISsite} function to obtain all of the information available for a particular USGS site (or sites) such as full station name, drainage area, latitude, and longitude. \texttt{readNWISsite} can also access information about multiple sites with a vector input.
 
 
-<<getSite, echo=TRUE>>=
+<<getSite, echo=TRUE, eval=FALSE>>=
 siteNumbers <- c("01491000","01645000") 
 siteINFO <- readNWISsite(siteNumbers)
 @
 
-A specific example piece of information can be retrieved, in this case a station name, as follows:
-
-<<siteNames2, echo=TRUE>>=
-siteINFO$station_nm
-@
-Site information is obtained from \url{http://waterservices.usgs.gov/rest/Site-Test-Tool.html}
+Site information is obtained from:
+\url{http://waterservices.usgs.gov/rest/Site-Test-Tool.html}
 
 Information on the returned data can be found with the \texttt{comment} function as described in section \ref{sec:metadata}.
 
@@ -385,7 +381,7 @@ To discover what data is available for a particular USGS site, including measure
 In the following example, we limit the retrieved data to only daily data. The default for \texttt{"}service\texttt{"} is \enquote{all}, which returns all of the available data for that site. Likewise, there are arguments for parameter code (\texttt{parameterCd}) and statistic code (\texttt{statCd}) to filter the results. The default for both is to return all possible values (\enquote{all}). The returned \texttt{"}count\_nu\texttt{"} for \texttt{"}uv\texttt{"} data is the count of days with returned data, not the actual count of returned values.
 
 
-<<getSiteExtended, echo=TRUE>>=
+<<getSiteExtended, echo=TRUE, eval=FALSE>>=
 # Continuing from the previous example:
 # This pulls out just the daily, mean data:
 
@@ -395,7 +391,7 @@ dailyDataAvailable <- whatNWISdata(siteNumbers,
 
 @
 
-<<tablegda, echo=FALSE,results='asis'>>=
+<<tablegda, echo=FALSE,eval=FALSE>>=
 tableData <- with(dailyDataAvailable, 
       data.frame( 
       siteNumber= site_no,
@@ -424,7 +420,24 @@ print(xtable(tableData,label="tab:gda",
 
 @
 
-
+\begin{table}[ht]
+\caption{Reformatted version of output from \texttt{whatNWISdata} function for the Choptank River near Greensboro, MD, and from Seneca Creek at Dawsonville, MD from the daily values service [Some columns deleted for space considerations]} 
+\label{tab:gda}
+{\footnotesize
+\begin{tabular}{rllllll}
+  \hline
+ & \multicolumn{1}{c}{\textbf{\textsf{siteNumber}}} & \multicolumn{1}{c}{\textbf{\textsf{srsname}}} & \multicolumn{1}{c}{\textbf{\textsf{startDate}}} & \multicolumn{1}{c}{\textbf{\textsf{endDate}}} & \multicolumn{1}{c}{\textbf{\textsf{count}}} & \multicolumn{1}{c}{\textbf{\textsf{units}}} \\ 
+  \hline
+ & 01491000 & Temperature, water & 2010-10-01 & 2012-05-09 & 529 & deg C \\ 
+  [5pt] & 01645000 & Stream flow, mean. daily & 1930-09-26 & 2015-02-19 & 30828 & ft$^3$/s \\ 
+  [5pt] & 01491000 & Stream flow, mean. daily & 1948-01-01 & 2015-02-19 & 24522 & ft$^3$/s \\ 
+  [5pt] & 01491000 & Specific conductance & 2010-10-01 & 2012-05-09 & 527 & $\mu$S/cm @25C \\ 
+  [5pt] & 01491000 & Suspended sediment concentration (SSC) & 1980-10-01 & 1991-09-30 & 3651 & mg/l \\ 
+  [5pt] & 01491000 & Suspended sediment discharge & 1980-10-01 & 1991-09-30 & 3652 & tons/day \\ 
+   \hline
+\end{tabular}
+}
+\end{table}
 
 See Section \ref{app:createWordTable} for instructions on converting an R data frame to a table in Microsoft\textregistered\ software Excel or Word to display a data availability table similar to Table \ref{tab:gda}. Excel, Microsoft, PowerPoint, Windows, and Word are registered trademarks of Microsoft Corporation in the United States and other countries.
 
@@ -436,19 +449,13 @@ See Section \ref{app:createWordTable} for instructions on converting an R data f
 %------------------------------------------------------------
 To obtain all of the available information concerning a measured parameter (or multiple parameters), use the \texttt{readNWISpCode} function:
 
-<<label=getPCodeInfo, echo=TRUE>>=
+<<label=getPCodeInfo, echo=TRUE, eval=FALSE>>=
 # Using defaults:
 parameterCd <- "00618" 
 parameterINFO <- readNWISpCode(parameterCd)
-colnames(parameterINFO)
 @
 
-A specific example piece of information, in this case parameter name, can be obtained as follows:
-
-<<siteNames, echo=TRUE>>=
-parameterINFO$parameter_nm
-@
-Parameter information can obtained from:
+More information on parameter codes can obtained from:
 
 \url{http://help.waterdata.usgs.gov/codes-and-parameters/parameters}
 
@@ -462,7 +469,7 @@ To obtain daily records of USGS data, use the \texttt{readNWISdv} function. The
 
 The dates (start and end) must be in the format \texttt{"}YYYY-MM-DD\texttt{"} (note: the user must include the quotes).  Setting the start date to \texttt{"}\texttt{"} (no space) will prompt the program to ask for the earliest date, and setting the end date to \texttt{"}\texttt{"} (no space) will prompt for the latest available date.
 
-<<label=getNWISDaily, echo=TRUE, eval=TRUE>>=
+<<label=getNWISDaily, echo=TRUE, eval=FALSE>>=
 
 # Choptank River near Greensboro, MD:
 siteNumber <- "01491000"
@@ -472,7 +479,6 @@ endDate <- "2012-09-30"
 
 discharge <- readNWISdv(siteNumber, 
                     parameterCd, startDate, endDate)
-names(discharge)
 @
 
 The column \texttt{"}datetime\texttt{"} in the returned data frame is automatically imported as a variable of class \texttt{"}Date\texttt{"} in R. Each requested parameter has a value and remark code column.  The names of these columns depend on the requested parameter and stat code combinations. USGS daily value qualification codes are often \texttt{"}A\texttt{"} (approved for publication) or \texttt{"}P\texttt{"} (provisional data subject to revision). A more complete list of daily value qualification codes can be found here:
@@ -481,8 +487,8 @@ The column \texttt{"}datetime\texttt{"} in the returned data frame is automatica
 
 Another example would be a request for mean and maximum daily temperature and discharge in early 2012:
 
-<<label=getNWIStemperature, echo=TRUE>>=
-
+<<label=getNWIStemperature, echo=TRUE, eval=FALSE>>=
+siteNumber <- "01491000"
 parameterCd <- c("00010","00060")  # Temperature and discharge
 statCd <- c("00001","00003")  # Mean and maximum
 startDate <- "2012-01-01"
@@ -493,6 +499,13 @@ temperatureAndFlow <- readNWISdv(siteNumber, parameterCd,
 
 @
 
+<<label=getNWIStemperature2, echo=FALSE, eval=TRUE>>=
+filePath <- system.file("extdata", package="dataRetrieval")
+fileName <- "temperatureAndFlow.RData"
+fullPath <- file.path(filePath, fileName)
+load(fullPath)
+
+@
 
 The column names can be shortened and simplified using the \texttt{renameNWISColumns} function.  This is not necessary, but may streamline subsequent data analysis and presentation. Site information, daily statistic information, and measured parameter information is attached to the data frame as attributes. This is discused further in section \ref{sec:metadata}.
 
@@ -549,7 +562,7 @@ There are occasions where NWIS values are not reported as numbers, instead there
 %------------------------------------------------------------
 Any data collected at regular time intervals (such as 15-minute or hourly) are known as \enquote{unit values.} Many of these are delivered on a real time basis and very recent data (even less than an hour old in many cases) are available through the function \texttt{readNWISuv}.  Some of these unit values are available for many years, and some are only available for a recent time period such as 120 days.  Here is an example of a retrieval of such data.  
 
-<<label=readNWISuv, echo=TRUE>>=
+<<label=readNWISuv, eval=FALSE>>=
 
 parameterCd <- "00060"  # Discharge
 startDate <- "2012-05-12" 
@@ -561,11 +574,6 @@ dischargeUnit <- renameNWISColumns(dischargeUnit)
 
 The retrieval produces a data frame that contains 96 rows (one for every 15 minute period in the day).  They include all data collected from the startDate through the endDate (starting and ending with midnight locally-collected time). The dateTime column is converted to \enquote{UTC} (Coordinated Universal Time), so midnight EST will be 5 hours earlier in the dateTime column (the previous day, at 7pm).
 
-
-<<dischargeData, echo=TRUE>>=
-head(dischargeUnit)
-@
-
 To override the UTC timezone, specify a valid timezone in the tz argument. Default is \texttt{""}, which will keep the dateTime column in UTC. Other valid timezones are:
 
 \begin{verbatim}
@@ -594,7 +602,7 @@ Data are retrieved from \url{http://waterservices.usgs.gov/rest/IV-Test-Tool.htm
 %------------------------------------------------------------
 To get USGS water quality data from water samples collected at the streamgage or other monitoring site (as distinct from unit values collected through some type of automatic monitor) we can use the function \texttt{readNWISqw}, with the input arguments: siteNumber, parameterCd, startDate, and endDate. Additionally, the argument \texttt{"}expanded\texttt{"} is a logical input that allows the user to choose between a simple return of datetimes/qualifier/values (expanded=FALSE), or a more complete and verbose output (expanded=TRUE). Expanded = TRUE includes such columns as remark codes, value qualifying text, and detection level for each parameter code. There also includes an argument \texttt{"}reshape\texttt{"}, that converts the expanded dataset to a \texttt{"}wide\texttt{"} format (each requested parameter code gets individual columns). The defaults are expanded=TRUE, and reshape=FALSE. 
 
-<<label=getQW, echo=TRUE, eval=TRUE>>=
+<<label=getQW, echo=TRUE, eval=FALSE>>=
  
 # Dissolved Nitrate parameter codes:
 parameterCd <- c("00618","71851")
@@ -626,12 +634,9 @@ comment(dfLong)
 %------------------------------------------------------------
 Groundwater level measurements can be obtained with the \texttt{readNWISgwl} function. Information on the returned data can be found with the \texttt{comment} function, and attached attributes as described in section \ref{sec:metadata}.
 
-<<gwlexample, echo=TRUE, eval=TRUE>>=
+<<gwlexample, echo=TRUE, eval=FALSE>>=
 siteNumber <- "434400121275801"
 groundWater <- readNWISgwl(siteNumber)
-
-names(groundWater)
-
 @
 
 %------------------------------------------------------------
@@ -641,13 +646,10 @@ names(groundWater)
 
 Peak flow data are instantaneous discharge or stage data that record the maximum values of these variables during a flood event.  They include the annual peak flood event but can also include records of other peaks that are lower than the annual maximum. Peak discharge measurements can be obtained with the \texttt{readNWISpeak} function. Information on the returned data can be found with the \texttt{comment} function and attached attributes as described in section \ref{sec:metadata}.
 
-<<peakexample, echo=TRUE, eval=TRUE>>=
+<<peakexample, echo=TRUE, eval=FALSE>>=
 siteNumber <- '01594440'
 peakData <- readNWISpeak(siteNumber)
 
-
-names(peakData)
-
 @
 
 
@@ -659,12 +661,10 @@ Rating curves are the calibration curves that are used to convert measurements o
 
 Rating curves can be obtained with the \texttt{readNWISrating} function.
 
-<<ratingexample, echo=TRUE, eval=TRUE>>=
+<<ratingexample, echo=TRUE, eval=FALSE>>=
 ratingData <- readNWISrating(siteNumber, "base")
 attr(ratingData, "RATING")
 
-names(ratingData)
-
 @
 
 
@@ -677,11 +677,9 @@ These data are the discrete measurements of discharge that are made for the purp
 
 Surface-water measurement data can be obtained with the \texttt{readNWISmeas} function.
 
-<<surfexample, echo=TRUE, eval=TRUE>>=
+<<surfexample, echo=TRUE, eval=FALSE>>=
 surfaceData <- readNWISmeas(siteNumber)
 
-names(surfaceData)
-
 @
 
 
@@ -697,7 +695,9 @@ specificCond <- readWQPqw('WIDNR_WQX-10032762',
                 'Specific conductance','2011-05-01','2011-09-30')
 @
 
-A tool for finding NWIS characteristic names can be found at: \url{http://www.waterqualitydata.us/public_srsnames.jsp}. 
+A tool for finding NWIS characteristic names can be found at: 
+
+\url{http://www.waterqualitydata.us/public_srsnames.jsp}
 
 \FloatBarrier
 
@@ -721,13 +721,10 @@ to discover many options for searching for NWIS sites. For example, you may want
 
 The following dataRetrieval code can be used to get those sites:
 
-<<siteSearch>>=
-sites <- whatNWISsites(bBox="-83.0,36.5,-81.0,38.5", 
-                      parameterCd="00010,00060",
+<<siteSearch, eval=FALSE>>=
+sites <- whatNWISsites(bBox=c(-83.0,36.5,-81.0,38.5), 
+                      parameterCd=c("00010","00060"),
                       hasDataTypeCd="dv")
-
-names(sites)
-nrow(sites)
 @
 
 
@@ -850,19 +847,13 @@ variableInfo <- attr(dischargeWI, "variableInfo")
 
 Data obtained from \texttt{readNWISpeak}, \texttt{readNWISmeas}, and \texttt{readNWISrating}, the \texttt{comment} attribute is useful.
 
-<<meta5, eval=TRUE, eval=FALSE>>=
+<<meta5, eval=FALSE>>=
 comment(peakData)
 
 #Which is equivalent to:
 # attr(peakData, "comment")
 @
 
-A subset (due to space considerations) of the \texttt{comment} metatdata is shown here:
-
-<<meta6, eval=TRUE, eval=TRUE>>=
-comment(peakData)[c(1:15,58:66)]
-
-@
 
 
 %------------------------------------------------------------ 
@@ -921,7 +912,7 @@ library(dataRetrieval)
 %------------------------------------------------------------
 There are a few steps that are required in order to create a table in Microsoft\textregistered\ software (Excel, Word, PowerPoint, etc.) from an R data frame. There are certainly a variety of good methods, one of which is detailed here. The example we will step through here will be to create a table in Microsoft Excel based on the data frame tableData:
 
-<<label=getSiteApp, echo=TRUE>>=
+<<label=getSiteApp, echo=TRUE, eval=FALSE>>=
 availableData <- whatNWISdata(siteNumber, "dv")
 dailyData <- availableData["00003" == availableData$stat_cd,]
 
@@ -933,7 +924,6 @@ tableData <- with(dailyData,
         Count=count_nu,
         Units=parameter_units)
       )
-tableData
 @
 
 First, save the data frame as a tab delimited file (you don't want to use comma delimited because there are commas in some of the data elements):
diff --git a/vignettes/figure/getNWIStemperaturePlot-1.pdf b/vignettes/figure/getNWIStemperaturePlot-1.pdf
index 0f18e7d998023d7fb221005eac8d6fce5289aaaf..f13c842b9a2021bc9db04f76595d427fbb129f04 100644
Binary files a/vignettes/figure/getNWIStemperaturePlot-1.pdf and b/vignettes/figure/getNWIStemperaturePlot-1.pdf differ