squall44 wrote:
Hello,
I have the following problem: I created an ecdf and a barplot. Unfortunatly,
the bars are not where I would like them to be (please see picture below).
http://www.nabble.com/file/p12877530/problem.gif
...
Can anyone tell me how I can adjust the bars to the left?
Hi Tobias (and anyone else who is interested),
This is a rather intriguing request. If there are people out there who
would like to place the bars in a barplot at arbitrary positions, I
could rather easily modify "barp" in the plotrix package to do this.
(I've already tried it with the "xpos" argument as a test, and it
produces nice overlapped bars, if you like that sort of thing.)
However, as a one-off, I would probably not include it in plotrix. I
attach the perhaps buggy altered code for Tobias if it helps. Use it
instead of barplot like this:
library(plotrix)
# add your path to the barp.R file to the following
source("barp.R")
barp(height,width=0.2,col="gray",xlim=c(-1.2,4.2),xpos=0:2)
If I get a few requests, I can thrash the code a bit and it will turn up
in the next version of plotrix.
Jim
barp<-function(height,width=0.4,names.arg=NULL,legend.lab=NULL,legend.pos="e",
col=NULL,border=par("fg"),main=NULL,xlab="",ylab="",xlim=NULL,ylim=NULL,
xpos=NULL,staxx=FALSE,staxy=FALSE,height.at=NULL,height.lab=NULL,
cex.axis=par("cex.axis"),cylindrical=FALSE,shadow=FALSE) {
if(is.data.frame(height)) its_ok<-is.numeric(unlist(height))
else its_ok<-is.numeric(height)
if(!its_ok) stop("barp can only display bars with numeric heights")
hdim<-dim(height)
if(is.null(hdim)) {
ngroups<-length(height)
if(is.null(xpos)) xpos<-1:ngroups
barcol=col
barpinfo<-list(x=xpos,y=height)
}
else {
ngroups<-hdim[2]
if(is.null(xpos)) xpos<-1:ngroups
if(!is.matrix(col) && length(col)==hdim[1])
barcol<-matrix(rep(col,each=ngroups),nrow=hdim[1],byrow=TRUE)
else barcol<-col
barpinfo<-list(x=matrix(rep(xpos,each=hdim[1]),ncol=hdim[2]),
y=as.matrix(height))
}
negy<-any(height<0)
if(is.null(xlim)) xlim<-c(xpos[1]-0.6,xpos[ngroups]+0.6)
if(is.null(ylim)) {
if(negy) miny<-min(height)*1.05
else miny<-0
ylim<-c(miny,max(height)*1.05)
}
else miny<-ylim[1]
plot(0,type="n",main=main,xlab=xlab,ylab=ylab,axes=FALSE,xlim=xlim,ylim=ylim,
xaxs="i",yaxs="i")
if(negy) abline(h=0)
if(is.null(names.arg)) names.arg<-xpos
if(staxx) {
axis(1,at=xpos,labels=rep("",ngroups),cex.axis=cex.axis)
staxlab(1,at=xpos,labels=names.arg,cex=cex.axis)
}
else axis(1,at=xpos,labels=names.arg)
if(is.null(height.at)) height.at<-pretty(ylim)
if(is.null(height.lab)) height.lab<-pretty(ylim)
if(staxy) {
axis(2,at=height.at,labels=rep("",length(height.lab)),cex.axis=cex.axis)
staxlab(2,at=height.at,labels=height.lab,cex=cex.axis)
}
else axis(2,at=height.at,labels=height.lab,cex.axis=cex.axis)
bottoms<-ifelse(negy,0,miny)
if(is.null(hdim)) {
if(shadow) {
for(bar in xpos)
polygon.shadow(c(bar-width,bar-width,bar+width,bar+width),
c(bottoms,height[bar],height[bar],bottoms),
offset=c(0.2*width,0.05*(height[bar]-ylim[2])))
}
if(cylindrical)
cylindrect(xpos-width,bottoms,xpos+width,height,col=barcol)
else rect(xpos-width,bottoms,xpos+width,height,col=barcol)
}
else {
bottoms<-matrix(bottoms,nrow=hdim[1],ncol=hdim[2])
for(subgroup in 1:hdim[1]) {
if(shadow) {
for(bar in xpos[subgroup,]) {
barleft<-bar-width+(subgroup-1)*2*width/hdim[1]
barright<-barleft+2*width/hdim[1]
polygon.shadow(c(barleft,barleft,barright,barright),
c(bottoms[bar],height[subgroup,bar],height[subgroup,bar],bottoms[bar]),
offset=c(0.2*width,0.05*(height[subgroup,bar]-ylim[2])))
}
}
if(cylindrical)
cylindrect(xpos[subgroup,]-width+(subgroup-1)*2*width/hdim[1],
bottoms[subgroup,],
xpos[subgroup,]-width+(subgroup)*2*width/hdim[1],
height[subgroup,],
col=barcol[subgroup,])
else rect(xpos[subgroup,]-width+(subgroup-1)*2*width/hdim[1],
bottoms[subgroup,],
xpos[subgroup,]-width+(subgroup)*2*width/hdim[1],
height[subgroup,],
col=barcol[subgroup,])
}
}
if(!is.null(legend.lab)) {
if(is.na(legend.pos[1])) {
cat("Click at the lower left corner of the legend\n")
legend.pos<-locator(1)
xjust<-yjust<-0
}
if(legend.pos[1] == "e") {
legend.pos<-emptyspace(barpinfo,bars=TRUE)
xjust<-yjust<-0.5
}
legend(legend.pos,legend=legend.lab,fill=col,xjust=xjust,yjust=yjust)
}
box()
return(barpinfo)
}
______________________________________________
[email protected] 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.
