Convert a Vector of Parameters to a Parameter Object of a Distribution


This function converts a vector of parameters to a parameter object of a distribution. The type of distribution is specified in the argument list: aep4, cau, exp, gam, gep, gev, glo, gno, gpa, gum, kap, kur, lap, lmrq, ln3, nor, pe3, ray, revgum, rice, st3, texp, wak, and wei. These abbreviations and only these are used in routing logic within lmomco. There is no provision for fuzzy matching. However, if the distribution type is not identified, then the function issues a warning, but goes ahead and creates the parameter list and of course can not check for the validity of the parameters. If one has a need to determine on-the-fly the number of parameters in a distribution as supported in lmomco, then see the dist.list function.


vec2par(vec, type, nowarn=FALSE, paracheck=TRUE, ...)



A vector of parameter values for the distribution specified by type.


Three character distribution type (for example, type='gev').


A logical switch on warning suppression. If TRUE then options(warn=-1) is made and restored on return. This switch is to permit calls in which warnings are not desired as the user knows how to handle the returned value—say in an optimization algorithm.


A logical controlling whether the parameters and checked for validity. Overriding of this check might be extremely important and needed for use of the distribution quantile function in the context of TL-moments with nonzero trimming.


Additional arguments for the are.par.valid call that is made internally.


If the distribution is a Reverse Gumbel (type=revgum) or Generalized Pareto (type=gpa), which are 2-parameter or 3-parameter distributions, the third or fourth value in the vector is the ζ of the distribution. ζ represents the fraction of the sample that is noncensored, or number of observed (noncensored) values divided by the sample size. The ζ represents censoring on the right, that is there are unknown observations above a threshold or the largest observed sample. Consultation of parrevgum or pargpaRC should elucidate the censoring discussion.


An R list is returned. This list should contain at least the following items, but some distributions such as the revgum have extra.


The type of distribution in three character format.


The parameters of the distribution.


Attribute specifying source of the parameters—“vec2par”.


If the type is not amongst the official list given above, then the type given is loaded into the type element of the returned list and an other element isuser = TRUE is also added. There is no isuser created if the distribution is supported by lmomco. This is an attempt to given some level of flexibility so that others can create their own distributions or conduct research on derivative code from lmomco.


W.H. Asquith

See Also

lmom2par, par2vec


para <- vec2par(c(12,123,0.5),'gev')
Q <- quagev(0.5,para)

my.custom <- vec2par(c(2,2), type='myowndist') # Think about making your own

Want to suggest features or report bugs for Use the GitHub issue tracker. Vote for new features on Trello.

comments powered by Disqus