## Univariate: Parametric
vg.dist <- VG(lambda = 1.1, mu = 10, sigma = 10, gamma = 2)
mean(vg.dist)
vcov(vg.dist)
ghyp.skewness(vg.dist)
ghyp.kurtosis(vg.dist)
## Univariate: Empirical
vg.sim <- rghyp(10000, vg.dist)
mean(vg.sim)
var(vg.sim)
## Multivariate: Parametric
vg.dist <- VG(lambda = 0.1, mu = c(55, 33), sigma = diag(c(22, 888)), gamma = 1:2)
mean(vg.dist)
vcov(vg.dist)
## Multivariate: Empirical
vg.sim <- rghyp(50000, vg.dist)
colMeans(vg.sim)
var(vg.sim)
Run the code above in your browser using DataLab