library(fBasics)
test <- nig.parameter(0, 2, 5, 1)
random <- rnig(1000000, alpha = test$alpha, beta = test$beta,
mu = test$mu, delta = test$delta)
mean(random)
var(random)
kurtosis(random)
skewness(random)
Run the code above in your browser using DataLab