ineq (version 0.2-13)

Pen: Pen's Parade

Description

plots Pen's Parade of a vector x

Usage

Pen(x, n = rep(1, length(x)), group = NULL, scaled = TRUE, abline = TRUE, add = FALSE, segments = NULL,   main = "Pen's Parade", ylab = NULL, xlab = NULL,  col = NULL, lwd = NULL, las = 1, fill = NULL, ...)

Arguments

x

a vector containing non-negative elements.

n

a vector of frequencies or weights, must be same length as x.

group

a factor coding different groups, must be same length as x. See also details.

scaled

logical. Should Pen's parade be divided by mean(x)?

abline

logical. Should a horizontal line for the mean be drawn?

add

logical. Should the plot be added to an existing plot?

segments

logical. Should histogram-like segments be drawn?

col

a (vector of) color(s) for drawing the curve.

fill

a (vector of) color(s) for filling the area under the curve.

xlab,ylab

axis labels. Suitable defaults depending on scaled and n are chosen.

main, lwd, las, ...

further high-level plot parameters.

Details

Pen's Parade is basically the inverse distribution function (standardized by mean(x)).

Pen allows for fine control of the layout---the graphical parameters col and fill can be vectorized if histogram-like segments are drawn (segments = TRUE)---but implements several heuristics in choosing its default plotting parameters. If a grouping factor group is given, the default is to draw segments with a grey-shaded filling. If no fill color is used, the default is to draw a thick blue curve. But as all of these are just defaults, they can of course easily be changed. See also the examples.

References

F A Cowell: Measurement of Inequality, 2000, in A B Atkinson / F Bourguignon (Eds): Handbook of Income Distribution, Amsterdam,

F A Cowell: Measuring Inequality, 1995 Prentice Hall/Harvester Wheatshef,

J Pen: Income Distribution, 1971, Harmondsworth: Allen Lane.

See Also

Lc, plot.Lc

Examples

# load and attach Philippine income data
data(Ilocos)
attach(Ilocos)
# plot Pen's Parade of income
Pen(income)
Pen(income, fill = hsv(0.1, 0.3, 1))

# 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)
Pen(x, n = n)
# create artificial grouping variable
myfac <- factor(c(1, 1, 1, 2, 2, 2, 3, 3, 3, 3))
Pen(x, n = n, group = myfac)