cumQdate: Cumulative flow calculation

View source: R/cumQdate.R

cumQdateR Documentation

Cumulative flow calculation


This function computes the first day of the calendar year at which a specific fraction of the cumulative flow for that year has been exceeded. Typically one looks for the point where half the cumulative flow has happened (fract = 0.5). The portion of the year being considered is set by paStart and paLong. The matrix returned has 2 columns: the first is the year (integer when the period of analysis ends), the second is the day of the year when the fraction has been exceeded. None of the rows will have any NA values.


cumQdate(eList, paStart = 10, paLong = 12, fract = 0.5)



named list with at least the Sample and INFO dataframes


numeric integer specifying the starting month for the period of analysis, 1 <= paStart <= 12, default is 10


numeric integer specifying the length of the period of analysis, in months, 1 <= paLong <= 12, default is 12


numeric fraction of the flow


It is common to use this type of analysis on the snowmelt period of the year. If (for example) we assume that snowmelt starts with the month of March and ends in July then we would set paStart = 3 and paLong = 5


annualSeries an integer matrix of two columns. The first column is the calendar year for the end of the period The second column is day of the year when the flow has exceeded the specified fraction of the entire period being considered


eList <- Choptank_eList
annualFlow <- cumQdate(eList)
mod1 <- lm(annualFlow[,2] ~ annualFlow[,1])

EGRET documentation built on April 18, 2023, 5:09 p.m.