seasAdjflow: Create Daily Seasonally-adjusted Log Flow Residuals

In leppott/baytrends033: Long Term Water Quality Trend Analysis

Create Daily Seasonally-adjusted Log Flow Residuals

Description

Create Daily Seasonally-adjusted Log Flow Residuals. The procedure to compute daily seasonally-adjusted log flow residuals is the following: 1) Check to make sure that the raw flow data are in the data set and that seasonally adjusted values have not already been computed. If so, no additional computations are performed. If not, then proceed with remaining steps. 2) Add date features if not already in the data set. 3) Compute and store Log (ln) flow as 'LogQ...' 4) Compute GAM model and store seasonally adjusted LogQ flow residuals as 'sa0LogQ...' 5) Smooth seasonally adjusted Log flow residuals by an averaged values based on dvAvgWin, dvAvgWgt, and dvAvgSides and store as "saxLogQ..."

Usage

seasAdjflow(
  dvFlow = dvFlow,
  siteNumber = NULL,
  dvAvgWin = c(7, 31),
  dvAvgWgt = "weighted",
  dvAvgSides = 1,
  plotResid = c(1)
)

Arguments

dvFlow	data frame with daily flow data. The flow and data qualifier code for each site is organized into two columns (e.g., "q01594440", "q01594440cd" for USGS gage 01594440)
siteNumber	a single USGS gage ID
dvAvgWin	Averaging window (days) for smoothing the residuals of the seasonally adjusted daily flow values.
dvAvgWgt	Averaging method ("uniform", "weighted" [default], or "centered") for creating weights. If using "weighted" then use dvAvgSides=1. If using "centered" then use dvAvgSides=2.
dvAvgSides	If dvAvgSides=1 only past values are used, if dvAvgSides=2 then values are centered around lag 0.
plotResid	plot residuals for selected averaging windows.

Value

return data frame of flow data with additional seasonally adjusted values

Examples

#Set Retrieval Parameters
yearStart   <- 1983
yearEnd     <- 2015
siteNumbers <- c("01578310")

# Regular Retrieval (default usage)
df <- getUSGSflow(siteNumbers, yearStart, yearEnd, fill=TRUE)
# Apply default smoothing
df <- seasAdjflow(df,"01578310")

leppott/baytrends033 documentation built on Feb. 17, 2024, 9:27 a.m.