dataOverview.r

#' Data Overview for WRTDS
#'
#' Gives a summary of data to be used for WRTDS analysis
#'
#' @param localDaily dataframe
#' @param localSample dataframe
#' @keywords data import USGS WRTDS
#' @export
#' @seealso \code{\link{mergeReport}}
#' @examples
#' # These examples require an internet connection to run
#' exDaily <- getDVData('01594440','00060', '1985-01-01', '1985-03-31', interactive=FALSE)
#' exSample <- getSampleData('01594440','01075', '1985-01-01', '1985-03-31', interactive=FALSE)
#' dataOverview(localDaily = exDaily, localSample = exSample)
dataOverview <- function(localDaily = Daily, localSample = Sample){
  numDays<-length(localDaily$Date)
  numSamples<-length(localSample$Date)
  numYears<-round(numDays/365.25,digits=0)
  cat("\n Discharge Record is",numDays,"days long, which is",numYears,"years")
  cat("\n First day of the discharge record is", as.character(localDaily$Date[1]),"and last day is",as.character(localDaily$Date[numDays]))
  cat("\n The water quality record has",numSamples,"samples")
  cat("\n The first sample is from", as.character(localSample$Date[1]),"and the last sample is from",as.character(localSample$Date[numSamples]))
  if(localSample$Date[1]<localDaily$Date[1]) cat("\n WE HAVE A PROBLEM first sample is from before the first daily discharge")	
  if(localSample$Date[numSamples]>localDaily$Date[numDays]) cat("\n WE HAVE A PROBLEM last sample is from after the last daily discharge")
  Qmin<-signif(min(localDaily$Q),digits=3)
  Qmean<-signif(mean(localDaily$Q),digits=3)
  Qmax<-signif(max(localDaily$Q),digits=3)
  Cmin<-signif(min(localSample$ConcHigh),digits=2)
  Cmean<-signif(mean(localSample$ConcHigh),digits=2)
  Cmax<-signif(max(localSample$ConcHigh),digits=2)
  cat("\n Discharge: Minimum, mean and maximum",Qmin,Qmean,Qmax)
  cat("\n Concentration: Minimum, mean and maximum",Cmin,Cmean,Cmax)
  pct<-sum(localSample$Uncen)
  pct<-((numSamples-pct)/numSamples)*100
  cat("\n Percentage of the sample values that are censored is",signif(pct,digits=2),"%")
}