Hi casperyc, While Jim Holtman's solution is quite neat, I thought I would add a multiple histogram to the discussion in case that was what you wanted: library(plotrix) barp(t(x.m[,2:5]),names.arg=x.m[,1],col=rainbow(4)) legend(20,2500,paste("V",2:5,sep=""),fill=rainbow(4)) Jim