rating.period: Load estimation with log log rating curve based on monthly or...

Description Usage Arguments Value Author(s) References See Also Examples

View source: R/rating.period.R

Description

Estimate the load using a monthly or annual based log-log rating curve between values of concentration and river flow at the time of sampling (Quilbe' et al., 2006). Multiple regression analyses are performed. The estimation is performed on the time period spanned by flow records.

Usage

1
rating.period(db, ncomp, period)

Arguments

db

An input data frame with at least three column. A column with DateTime in standard format (yyyy-mm-dd HH:MM:SS) labeled "datetime", a column with flow records in cubic meter per second labeled "flow", and columns with scattered concentration data in milligram per litre. Alternatively, it can be used the output matrix of the function db.union.

ncomp

Number of compounds in the input data frame for which the load must be estimated.

period

A mandatory argument that specify the time period in which the regression relationship must be calculated. It can be "month" for a monthly estimation of load, or it can be "year" for annual estimation of load.

Value

A matrix with the estimated load for the different compounds. If the optional argument is equal to "month", the load values are reported by month in different rows. If the optional argument is equal to "year", the load values are reported by year in different rows. The load is expressed in gram per period of estimation.

Author(s)

Veronica Nava

References

Quilbe', R., Rousseau, A. N., Duchemin, M., Poulin, A., Gangbazo, G., & Villeneuve, J. P. (2006). Selecting a calculation method to estimate sediment and nutrient loads in streams: Application to the Beaurivage River (Quebec, Canada). Journal of Hydrology, 326(1-4), 295-310. https://doi.org/10.1016/j.jhydrol.2005.11.008

See Also

db.union rating rsquared.period

Examples

1
2
3
4
data("flow.data2","conc.data2")
union<-db.union(flow.data2, conc.data1)
reg.periodM<-rating.period(union, 1, "month")
reg.periodY<-rating.period(union, 1, "year")

RiverLoad documentation built on Jan. 4, 2022, 5:07 p.m.