library(fmx)
(d2 = fmx('GH', A = c(1,6), B = 2, g = c(0,.3), h = c(.2,0), w = c(1,2)))
set.seed(100); hist(x2 <- rfmx(n = 1e3L, dist = d2))
# equidistant in probabilities
(p1 = QLMDp())
# equidistant in quantiles
(p2 = QLMDp(equidistant = 'quantile', x = x2))
Run the code above in your browser using DataLab