Plot Theoretical Lorenz Curves
Plotting method for objects of class "theorLc" (theoretical Lorenz
curves)
## S3 method for class 'theorLc' plot(x, parameter=NULL, xlab="p", ylab="L(p)", lwd=2, las=1, ...)
C Dagum: Income Distribution Models, 1983, in: Johnson / Kotz (Eds): Encyclopedia of Statistical Sciences Vol.4, 27-34.
J B McDonald: Some generalized functions for the size distribution of income, 1984, Econometrica 52, 647-664.
# income distribution of the USA in 1968 (in 10 classes) # x vector of class means, n vector of class frequencies x <- c(541, 1463, 2445, 3438, 4437, 5401, 6392, 8304, 11904, 22261) n <- c(482, 825, 722, 690, 661, 760, 745, 2140, 1911, 1024) # compute minimal Lorenz curve (= no inequality in each group) Lc.min <- Lc(x, n=n) # compute maximal Lorenz curve (limits of Mehran) Lc.max <- Lc.mehran(x,n) # plot both Lorenz curves in one plot plot(Lc.min) lines(Lc.max, col=4) # add the theoretic Lorenz curve of a Lognormal-distribution with variance 0.78 lines(Lc.lognorm, parameter=0.78) # add the theoretic Lorenz curve of a Dagum-distribution lines(Lc.dagum, parameter=c(3.4,2.6))
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.