On 06/15/2011 09:51 PM, bjmjarrett wrote:
Hi all,
I have scoured the archives of this forum but nothing quite seems to fit the
bill...
I would like to plot a graph displaying two variables (y axes) that share
date as the x axis. However, the date values for each variable are not the
same - for example, some parasitoids were not released on days that
collections from the trap took place, whilst sometimes releases did occur on
the same day. I would like to align them. My code is:
####
CollectionDate<-as.Date(CollectionDate,"%m/%d/%Y") # days that had
collections
release.date.Total<-as.Date(release.date.Total,"%m/%d/%Y") # days that had
releases
par(mar=c(5,4,4,4)+0.1)
plot(CollectionDate[data$Trap=="DPAU1"],TotalP[data$Trap=="DPAU1"],type="n",ylim=c(0,80),xlab="Time",ylab="Number
of egg masses") # I am using one trap of many
lines(CollectionDate[data$Trap=="DPAU1"],TotalP[data$Trap=="DPAU1"],lwd=2)
par(new=T)
plot(release.date.Total[data$Trap=="DPAU1"],parasitoid.total[data$Trap=="DPAU1"],axes=F,type="h",lty=1,col="red",xlab="",ylab="",lwd=1.2)
axis(4,las=1)
mtext(side=4,line=3,"Total parasitoids released")
###
As in the above code, I suppressed the axes for the second plot, but if I
don't the x axis (date) has two tick marks for each year.
Is there a way to align the date objects to begin with without compromising
the graphs themselves?
Hi Ben,
This might do what you want:
collection_dates<-as.Date(paste(sort(sample(1:30,15)),"05","2011",sep="/"),
"%d/%m/%Y")
release_dates<-as.Date(paste(sort(sample(1:30,15)),"05","2011",sep="/"),
"%d/%m/%Y")
release_n<-sample(50:100,15)
collect_n<-sample(50:100,15)
xtickpos<-as.numeric(as.Date(paste(c(1,10,20,30),"05","2011",sep="/"),"%d/%m/%Y"))
xticklab<-format(as.Date(paste(c(1,10,20,30),"05","2011",sep="/"),"%d/%m/%Y"))
library(plotrix)
twoord.plot(sort(collection_dates),collect_n,sort(release_dates),release_n,
lylim=c(45,105),rylim=c(45,105),xlab="Date",
ylab="Parasitoids collected",main="Collect/Release",
rylab="Parasitoids released",xtickpos=xtickpos,xticklab=xticklab)
legend(as.numeric(as.Date("1/5/2011","%d/%m/%Y")),55,
c("Collection","Release"),lty=1,col=1:2,pch=1:2)
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.
