Andreas Christoffersen wrote:
Hi there,
I would like to make a 'bumps chart' like the ones described e.g.
here: http://junkcharts.typepad.com/junk_charts/bumps_chart/
Purpose: I'd like to plot the proportion of people in select countries
living for less then one USD pr day in 1994 and 2004 respectively. I
have already constructed a barplot - but I think a bumps chart would
be better
# The barplot and data
countries <- c("U-lande", "Afrika syd for sahara", "Europa og
Centralasien", "Lantinamerika og Caribien","Mellemøstenog Nordafrika",
"Sydasien","ØStasien og stillehaveet", "Kina", "Brasilien")
poor_1990 <- c(28.7,46.7,0.5,10.2,2.3,43,29.8,33,14)
poor_2004 <- c(18.1,41.1,0.9,8.6,1.5,30.8,9.1,9.9,7.5)
poor <- cbind(poor_1990,poor_2004)
rownames(poor) <- countries
oldpar <- par(no.readonly=T)
par <- par(mar=c(15,5,5,1))
png("poor.png")
par <- par(mar=c(15,5,5,1))
barplot(t(poor[order(poor[,2]),]),beside=T,col=c(1,2),las=3,ylab="%
poor",main="Percent living for < 1 USD per day (1993
prices)",ylim=c(0,50))
legend("topleft",c("1990","2004"),fill=c(1,2),bty="n")
par(oldpar)
dev.off()
I Guess I need to start with an normal plot? Something like the below
- but there is a loong way to go...
# A meager start - how to finish my bumps chart
plot(c(rep(1,9),rep(2,9)),c(fattig_1990,fattig_2004),type="b",ann=F)
Hi Andreas,
Not too hard. Try this:
bump.plot<-function (y1,y2=NULL,top.labels=NULL,left.labels=NULL,
right.labels=NULL,rank=TRUE,descending=TRUE,linecol=par("fg"),
mar=c(2,2,4,2),...) {
if(missing(y1))
stop("Usage: spread.labels(y1,y2,labels,...)")
ydim<-dim(y1)
if(is.null(ydim) && is.null(y2))
stop("y1 must be a matrix or data frame if y2 is NULL")
oldmar<-par("mar")
par(mar=mar)
if(is.null(y2) && ydim[2] > 1) {
y2<-y1[,2]
y1<-y1[,1]
}
if(rank) {
left.labels<-left.labels[order(y1)]
right.labels<-right.labels[order(y1)]
y1<-rank(y1)
y2<-rank(y2)
if(descending) {
left.labels<-rev(left.labels)
right.labels<-rev(right.labels)
y1<-rev(y1)
y2<-rev(y2)
}
}
ny<-length(y1)
plot(c(rep(1,ny),rep(2,ny)),c(rev(y1),rev(y2)),xlim=c(0,3),xlab="",ylab="",
axes=FALSE,...)
segments(rep(1,ny),y1,rep(2,ny),y2,col=linecol)
text(0.8,rev(y1),left.labels,adj=1)
text(2.2,rev(y2),right.labels,adj=0)
par(mar=oldmar)
}
# use the above data
bump.plot(poor_1990,poor_2004,c("1990","2004"),countries,countries,
linecol=1:20,pch=16,main="Bumps plot")
Jim
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.