# Comparing GHD fBasics and ghyp functions
set.seed(4242)
plot(density(rhyp(100,2,0,1,0),bw = .5),col="black",ylim=c(0,.6), main="")
set.seed(4242)
lines(density(rghyp(100, ghyp.ad(lambda=1, alpha=2, delta=1, mu=0, beta=0)), bw=.5), col="red")
# Comparing random samples
set.seed(4242)
fb_ghd <- rhyp(100,2,0,1,0)
set.seed(4242)
gh_ghd <- rghyp(100, ghyp.ad(lambda=1, alpha=2, delta=1, mu=0, beta=0))
tail(cbind(fb_ghd, gh_ghd))