Dear Contributors, I have a data of the form: 4 8 10 8590 12516 4 8 11 8641 98143 4 8 12 8705 98916 4 8 13 8750 89911 4 8 14 8685 104835 4 8 15 8629 121963 4 8 16 8676 77655 4 8 17 8577 81081 4 8 18 8593 83385 4 8 19 8642 112164 4 8 20 8708 103684 4 8 21 8622 83982 4 8 22 8593 75944 4 8 23 8600 97036 4 8 24 8650 104911 4 8 25 8730 114098 4 8 26 8731 99421 4 8 27 8715 85707 4 8 28 8717 81273 4 8 29 8739 106462 4 8 30 8684 110635 4 8 31 8713 105214 4 9 1 8771 92456 4 9 2 8759 109270 4 9 3 8762 99150 4 9 4 8730 77306 4 9 5 8780 86324 4 9 6 8804 90214 4 9 7 8797 99894 4 9 8 8863 95177 4 9 9 8873 95910 4 9 10 8827 108511 4 9 11 8806 115636 4 9 12 8869 85542 4 9 13 8854 111018 4 9 14 8571 93247 4 9 15 8533 85105 4 9 16 8553 114725 4 9 17 8561 122195 4 9 18 8532 100945 4 9 19 8560 108552 4 9 20 8634 108707 4 9 21 8646 117420 4 9 22 8633 113823 4 9 23 8680 82763 4 9 24 8765 121072 4 9 25 8756 89835 4 9 26 8750 104578 4 9 27 8790 88429 4 9 28 8824 84022 4 9 29 8843 80413 4 9 30 8795 118462 4 10 1 8795 139761 4 10 2 8805 103049 4 10 3 8826 104996 4 10 4 8841 104496 4 10 5 8882 75603 4 10 6 8865 89768 4 10 7 8875 93353 4 10 8 8898 98410 4 10 9 8920 114540 4 10 10 8937 95220 4 10 11 8926 104083 4 10 12 8912 106089 4 10 13 8898 116228 4 10 14 8926 149610 4 10 15 8954 109594 4 10 16 8941 75008 4 10 17 8952 119182 4 10 18 9003 151011 4 10 19 9002 146797 4 10 20 8995 145769 4 10 21 8969 121248 4 10 22 8954 107991 4 10 23 8992 128 4 10 24 8981 23036 4 10 25 8920 137485 4 10 26 8899 131756 4 10 27 8913 108729 4 10 28 8874 109478 4 10 29 8846 119627 4 10 30 8867 89999 4 10 31 8868 64833 4 11 1 9004 95864 4 11 2 9028 82322 4 11 3 8969 95591 4 11 4 8932 69378 4 11 5 8929 74281 4 11 6 8916 103261 4 11 7 8807 92473 4 11 8 8449 84344 4 11 9 8484 127415 4 11 10 8148 123826 4 11 11 8282 100029 4 11 12 8305 76205 4 11 13 8380 105162 4 11 14 8530 119533 4 11 15 8642 106490 4 11 16 8780 114771 4 11 17 8890 55593 4 11 18 8962 227 4 11 19 8949 109699 4 11 20 8974 86004 4 11 21 8956 74496 4 11 22 8881 109350 4 11 23 8872 134020 4 11 24 8847 105212 4 11 25 8868 91512 where columns 1= year 2004, 2 = month, 3= day, 4= CR data, 5= Lightning data.
The data range is 2004/8/10 to 2008/8/22. With the code below: data <- read.table("CRandWWLLNremzro", col.names = c("year", "month", "day", "CR","WW")) new.century <- data$year < 50 data$year <- ifelse(new.century, data$year + 2000, data$year + 1900) data$date <- as.Date(ISOdate(data$year, data$month, data$day)) x = data$date CR = data$CR WWLLN=data$WW Year<-x Li<-WWLLN CR<-CR setEPS() postscript("twoinone2.eps") par(mar = c(5, 4, 4, 4) + 2) # Leave space for z axis library(plotrix) plot(Year,Li,pch=16,axes=F,xlab="",ylab="",type="l",col="black") axis(2, col="black",las=1) ## las=1 makes horizontal labels mtext("Lightning Stroke/day", side=2, line=4) #box() par(new=TRUE) plot(Year,CR, pch=15, xlab="",ylab="",axes=FALSE,type="l",col="red") mtext("GCR count/day",side=4,col="red",line=3) axis(4, ylim=c(-7499,9684), col="red",col.axis="red",las=1) axis(side=1) ## mtext("YEAR",side=1,col="black",line=2.5) mtext("Long-term Variation betwenn WWLLN and GCRs",side=3,col="black",line=2.5) legend("topleft",col=c("red","black"),lty=1,legend=c("GCRs","WWLLN")) dev.off() I got the correct axes except that the x-axis is not what I want. I want the date to appear on the x-axis. i actually got that by using: axis(side=1,at=c(1,400,800,1200,1500),labels=c("2004","2005","2006","2007","2008")) But there is a little time shift, indicating that the above line may not have been assigned correctly. I then fiddled with axis(side=1) to see if the correct date can naturally appear. I could not succeed. I can attach the plot I generated. Instead of date appearing on the x-axis, ordinary numbers appeared. I guess the best way of getting the correct date on x-axis is to allow the system to fix the date as data points are many. But I don't know how to do that. Thank you in advance for your usual kind help. Best regards Ogbos ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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.