
These objects store different parameter sets of the generalized inverse Gaussian distribution as matrices for testing or demonstration purposes.
The parameter sets gigSmallParam
and gigLargeParam
give
combinations of values of the parameters gigSmallParam
,
the values of gigLargeParam
, the values of
gigSmallParam
gigLargeParam
gigSmallParam
: a 125 by 3 matrix;
gigLargeParam
: a 1100 by 3 matrix.
David Scott d.scott@auckland.ac.nz
data(gigParam)
## Check values of chi and psi
plot(gigLargeParam[, 1], gigLargeParam[, 2])
### Check all three parameters
pairs(gigLargeParam,
labels = c(expression(chi),expression(psi),expression(lambda)))
## Testing the accuracy of gigMean
for (i in 1:nrow(gigSmallParam)) {
param <- gigSmallParam[i, ]
x <- rgig(1000, param = param)
sampleMean <- mean(x)
funMean <- gigMean(param = param)
difference <- abs(sampleMean - funMean)
print(difference)
}
Run the code above in your browser using DataLab