Bumps chart in R
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