# Let us parametize some "income" distribution.
A <- vec2par(c(123, 264, 2.11), type="gov")
riglmomco(0.5, A)
if (FALSE) {
F <- nonexceeds(f01=TRUE)
plot(F, riglmomco(F,A), type="l",
xlab="NONEXCEEDANCE PROBABILITY", ylab="INCOME GAP RATIO")}
Run the code above in your browser using DataLab