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