| 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)