lfstat (version 0.9.4)

rfa: Regional Frequency Analysis

Description

This function uses J.R.M. Hosking's package produce an object of class "rfd", containing the specification of the regional frequency distribution.

Usage

rfa(lflist, n = 7, event = 100, dist =  c("wei","gev","ln3","gum","pe3"))

Arguments

lflist

A list of "lfobj"s

n

MAM-n is used (e.g. n=7 means MAM7)

event

A value for T, e.g. event = 100 means the 100 years extreme low flow event

dist

A vector of distribution to fit, the names are according to Hostking's in his "lmom" package. Can be an of "wei","gev","ln3","gum","pe3".

References

Manual on Low-flow Estimation and Prediction, Operational Hydrology Report No. 50, Koblenz 2009

J. R. M. Hosking (2012). L-moments. R package, version 1.6. URL: http://CRAN.R-project.org/package=lmom.

See Also

lmom,lmomRFA

Examples

data(ngaruroro)

#Toy example to get some more "rivers"
seventies <- subset(ngaruroro, hyear %in% 1970:1979)
eighties <- subset(ngaruroro, hyear %in% 1980:1989)
nineties <- subset(ngaruroro, hyear %in% 1990:1999)

toyrfa <- rfa(list(seventies,eighties,nineties), n=3,dist = "gev")

# Now you can work on using Hoskings lmomRFA-package, e.g.
require(lmomRFA)
regquant(c(1/1000,1/100),toyrfa)
sitequant(1/100,toyrfa)