# NOT RUN {
# Let us parametize some "income" distribution.
A <- vec2par(c(123, 264, 2.11), type="gov")
riglmomco(0.5, A)
# }
# NOT RUN {
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