lmr <- lmom.ub(c(123,34,4,654,37,78))
g <- pargam(lmr)
quagam(0.5,g)
# generate 50 random samples from this fitted parent
Qsim <- rlmomco(5000,g)
# compute the apparent gamma parameter for this parent
gsim <- pargam(lmoms(Qsim))
Run the code above in your browser using DataLab