data(skewhypParam)
### Testing the accuracy of skewhypMean
for (i in 1:nrow(skewhypSmallParam)) {
param <- skewhypSmallParam[i, ]
x <- rskewhyp(1000, param = param)
sampleMean <- mean(x)
distMean <- skewhypMean(param = param)
difference <- abs(sampleMean - distMean)
print(difference)
}
Run the code above in your browser using DataLab