freq.curve.all {lmomco} | R Documentation |
This function is dispatcher on top of the suite of freq.curve.CCC
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
, and wak
. The Cauchy distribution
is not called 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).
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,...)
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 deviations. 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 , and pe3 . |
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 TRUE . |
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 . |
... |
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
freq.curve.exp
,
freq.curve.gam
,
freq.curve.gev
,
freq.curve.gld
,
freq.curve.glo
,
freq.curve.gno
,
freq.curve.gpa
,
freq.curve.gum
,
freq.curve.kap
,
freq.curve.nor
,
freq.curve.pe3
, and
freq.curve.wak
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)