rugarch (version 1.0-3)

rgarchdist: Distribution: rugarch distribution functions

Description

Density, distribution function, quantile function, random generation and fitting from the univariate distributions implemented in the rugarch package, with functions for skewness and excess kurtosis given density skew and shape parameters. ll{ rgarchdist rugarch univariate distributions, fitdist MLE parameter fit for the rugarch univariate distributions, }

Usage

ddist(distribution = "norm", y, mu = 0, sigma = 1, lambda = -0.5, skew = 1, 
shape = 5)
pdist(distribution = "norm", q, mu = 0, sigma = 1, lambda = -0.5, skew = 1, 
shape = 5)
qdist(distribution = "norm", p, mu = 0, sigma = 1, lambda = -0.5, skew = 1, 
shape = 5)
rdist(distribution = "norm", n, mu = 0, sigma = 1, lambda = -0.5, skew = 1, 
shape = 5)
fitdist(distribution = "norm", x, control=list())
dskewness(distribution = "norm", skew = 1, shape = 5, lambda = -0.5)
dkurtosis(distribution = "norm", skew = 1, shape = 5, lambda = -0.5)

Arguments

distribution

The distribution name. Valid choices are norm, snorm, std, sstd, ged, sged, nig, jsu.

mu, sigma, skew, shape

location, scale and skewness and shape parameters (see details).

lambda

The additional shape parameter for the Generalized Hyperbolic and NIG distributions.

n

The number of observations.

p

A numeric vector of probabilities.

y, q

A numeric vector of quantiles.

x

A univariate dataset (for fitting routine).

control

Control parameters passed to the solnp solver.

Value

• d* returns the density, p* returns the distribution function, q* returns the quantile function, and r* generates random deviates, all values are numeric vectors. fitdist returns a list with the following components:
• par
  The best set of parameters found.
• value
  The likelihood values of the optimization (vector whose length represents the number of major iterations).
• convergence
  An integer code. 0 indicates successful convergence.
• lagrange
  The lagrange multiplier value at convergence.
• h
  The hessian at the solution.
• xineq0
  The value of the inequality constraint multiplier (NULL for the distribution fit problems).
• dskewness returns the skewness of the distribution. dkurtosis returns the excess kurtosis of the distribution.

Details

For the dQuote{nig} and ghyp distributions, the shape, skew and lambda are transformed from the zeta-rho to the alpha-beta parametrization and then scaled by the mean and standard deviation. The fitting routines use the solnp solver and minimize the negative of the log-likelihood. The dskewness and dkurtosis functions take as inputs the distribution name, skew and shape parameters and return the skewneness and excess kurtosis of the distribution. The functions are not at present vectorized.

References

Johnson, N. L. 1954, Systems of frequency curves derived from the first law of Laplace, Trabajos de Estadistica, 5, 283--291. Barndorff-Nielsen, O. E. 1995, Normal inverse Gaussian processes and the modeling of stock returns, mimeo: Univ.of Aarhus Denmark. Fernandez C. and Steel, M.F.J. 1998, On Bayesian Modelling of Fat Tails and Skewness, Journal of the American Statistical Association, 359--371.