data(vgParam)
## Testing the accuracy of vgMean
for (i in 1:nrow(smallParam)) {
param <- smallParam[i,]
x <- rvg(10000,param = param)
sampleMean <- mean(x)
funMean <- vgMean(param = param)
difference <- abs(sampleMean - funMean)
print(difference)
}
Run the code above in your browser using DataLab