freq.curve.all {lmomco} | R Documentation |
This function is dispatcher on top of the suite of quaCCC
functions that compute frequency curves for the L-moments. Frequency curves in hydrologic science is a term typically renaming the more conventional quantile function. The notation CCC
represents the three character notation for the distribution: exp
, gam
, gev
, gld
, glo
, gno
, gpa
, gum
, kap
, nor
, pe3
, wak
, and wei
. The Cauchy distribution is not used because of its dependency on trimmed L-moments and its general lack of use in applied research problems (at least those familiar to the author). The nonexceedance probabilities to construct the curves are derived from nonexceeds
.
freq.curve.all(lmom,aslog10=FALSE,asprob=TRUE, no2para=FALSE,no3para=FALSE, no4para=FALSE,no5para=FALSE, step=FALSE,show=FALSE, xmin=NULL,xmax=NULL,xlim=NULL, ymin=NULL,ymax=NULL,ylim=NULL, exp=TRUE,gam=TRUE,gev=TRUE,gld=FALSE, glo=TRUE,gno=TRUE,gpa=TRUE,gum=TRUE, kap=TRUE,nor=TRUE,pe3=TRUE,wak=TRUE, wei=TRUE,...)
lmom |
A L-moment object from lmom.ub or similar. |
aslog10 |
Compute log10 of quantiles—note that NaNs produced in: log(x, base)will be produced for less than zero values. |
asprob |
The qnorm function is used to convert nonexceedance probabilities, which are produced by nonexceeds , to standard normal deviates. The normal distribution will plot as straight line when this argument is TRUE |
no2para |
If TRUE , do not run the 2-parameter distributions: exp , gam , gum , and nor . |
no3para |
If TRUE , do not run the 3-parameter distributions: gev , glo , gno , gpa , pe3 , and wei . |
no4para |
If TRUE , do not run the 4-parameter distributions: kap and gld . |
no5para |
If TRUE , do not run the 5-parameter distributions: wak . |
step |
Shows incremental processing of each distribution. |
show |
Plots all the frequency curves in a simple (crowded) plot . |
xmin |
Minimum x-axis value to use instead of the automatic value determined from the nonexceedance probabilities. This argument is only used is show=TRUE . |
xmax |
Maximum x-axis value to use instead of the automatic value determined from the nonexceedance probabilities. This argument is only used is show=TRUE . |
xlim |
Both limits of the x-axis. This argument is only used is show=TRUE . |
ymin |
Minimum y-axis value to use instead of the automatic value determined from the nonexceedance probabilities. This argument is only used is show=TRUE . |
ymax |
Maximum y-axis value to use instead of the automatic value determined from the nonexceedance probabilities. This argument is only used is show=TRUE . |
ylim |
Both limits of the y-axis. This argument is only used is show=TRUE . |
exp |
A logical switch on computation of corresponding distribution—default is TRUE . |
gam |
A logical switch on computation of corresponding distribution—default is TRUE . |
gev |
A logical switch on computation of corresponding distribution—default is TRUE . |
gld |
A logical switch on computation of corresponding distribution—default is FALSE . |
glo |
A logical switch on computation of corresponding distribution—default is TRUE . |
gno |
A logical switch on computation of corresponding distribution—default is TRUE . |
gpa |
A logical switch on computation of corresponding distribution—default is TRUE . |
gum |
A logical switch on computation of corresponding distribution—default is TRUE . |
kap |
A logical switch on computation of corresponding distribution—default is TRUE . |
nor |
A logical switch on computation of corresponding distribution—default is TRUE . |
pe3 |
A logical switch on computation of corresponding distribution—default is TRUE . |
wak |
A logical switch on computation of corresponding distribution—default is TRUE . |
wei |
A logical switch on computation of corresponding distribution—default is TRUE . |
... |
Additional parameters are passed to the parameter estimation routines such as parexp . |
An extensive R data.frame
of frequency curves. The nonexceedance probability values, which are provided by nonexceeds
, are the first item in the data.frame
under the heading of nonexceeds
. If a particular distribution could not be fit to the L-moments of the data; this particular function returns zeros so that a data.frame
can be returned.
W.H. Asquith
quaexp
,
quagam
,
quagev
,
quagld
,
quaglo
,
quagno
,
quagpa
,
quagum
,
quakap
,
quanor
,
quape3
,
quawak
, and
quawei
.
L <- vec2lmom(c(35612,23593,0.48,0.21,0.11)) freq.curve.all(L,gld=FALSE) freq.curve.all(L,step=TRUE,no2para=TRUE,no4para=TRUE)