# TODO: Add comment
# 
# Author: E.Korsching 10.9.2009
###############################################################################



my.barplot<- function(x,y,col=c("black","blue"),width=1,xlab=x,density=-1,tick=T,adj=0.5,srt=0,horiz=F)
{
	# plot bars

	for(i in 1:length(x)){
		at.x <- x[i]
		par(adj = 1)
	    if(!horiz){
		    xb <- c(x[i]-width/2, x[i]+width/2, x[i]+width/2, x[i]-width/2)
            yb <- c(0, 0, y[i], y[i])
        }else{
		    yb <- c(x[i]-width/2, x[i]-width/2, x[i]+width/2, x[i]+width/2)
		    xb <- c(0, y[i], y[i], 0)
    	}
		
		polygon(xb, yb, density=density, col=col[2], border=col[1])
		
		if(!horiz){
			if(tick)axis(1, at=at.x, labels=xlab[i], tick=T)		    # x, bars vertical
		}else{
			if(tick)axis(2, at=at.x, labels=xlab[i], tick=T, adj=1)	    # y, bars horizontal
		}
	}
}