rgarchdist {rugarch} | R Documentation |
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.
rgarchdist | rugarch univariate distributions, |
fitdist | MLE parameter fit for the rugarch univariate distributions, |
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)
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 |
For the dQuotenig 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.
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.
Diethelm Wuertz for the Rmetrics R-port of the “norm”, “snorm”,
“std”, “sstd”, “ged”, “sged” and “nig”
distrbutions.
Rigby, R. A. and Stasinopoulos D. M for the JSU distribution in the
gamlss
package.
Alexios Ghalanos for rugarch implementation and higher moment distribution
functions.
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.