Dear Bert and Jim, Here is the code: data <- read.table("CLMX.txt", col.names = c("year", "month", "day", "CLMX"))
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)) x1 = data$date Li<-data$CLMX ########### data <- read.table("DPRV.txt", col.names = c("year", "month", "day", "DPRV")) 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)) x2 = data$date CR<-data$DPRV ################### date_x1<-seq(as.Date("1953-01-02"),as.Date("2006-11-15"),length.out=1575) date_x2<-seq(as.Date("1957-07-26"),as.Date("1994-07-17"),length.out=1203) pdf("CLMXandDPRV.pdf") par(mar = c(5, 4, 4, 4) + 2) library(plotrix) plot(date_x1,Li,pch=16,axes=FALSE,xlab="",ylab="",type="l",col="black",ylim=c(-0.5001663,-10.27402)) points(date_x1,Li,col="black") axis(2, ylim=c(-0.5001663,-10.27402),col="black",las=1) mtext("FD (%)", side=2, line=4) text(as.Date("1992-02-21"),3523,"FD1") text(as.Date("1972-02-28"),-8,"Total FDs = 1203",col="red") text(as.Date("1972-03-02"),-9,"Total FDs = 1575") #box() par(new=TRUE) plot(date_x2,CR, pch=16, xlab="",ylab="",ylim=c(-0.5030138,-14.39884),axes=FALSE,type="l",col="red") points(date_x2,CR,col="red") mtext("FD (%)",side=4,col="red",line=3) axis(4, ylim=c(-0.5030138,-14.39884), col="red",col.axis="red",las=1) mtext("Date",side=1,col="black",line=2.5) legend("topleft",col=c("red","black"),lty=1,legend=c("DPRV","CLMX")) dev.off() Thanks again. Best wishes Ogbos On Sun, Aug 25, 2019 at 4:52 PM Ogbos Okike <giftedlife2...@gmail.com> wrote: > > Dear Jim and Bert, > Thank your for looking into this for me. I will try to reproduce a > part of my data and the plot I have. > > Data 1: > 53 01 02 -1.28792560381641 > 53 01 06 -1.1854773963453 > 53 01 08 -1.55920165458006 > 53 01 15 -1.29196482429683 > 53 01 20 -1.06082194329819 > 53 01 22 -1.15430411152234 > 53 01 24 -1.14775155345262 > 53 01 26 -1.19116423468105 > 53 01 28 -0.924542121164923 > 53 01 31 -0.869543637175866 > 53 02 03 -0.764466677206437 > 53 02 05 -0.707705051599105 > 53 02 07 -0.750908466521414 > 53 02 16 -1.6848903649715 > 53 02 20 -0.760876361940755 > 53 02 22 -0.803816591115474 > 53 02 26 -0.818160891300922 > 53 02 28 -0.831013020625678 > 53 03 04 -1.06661138281899 > till the year 2006. The length is 1575 > > Data B: > 57 07 26 -3.62351759422703 > 57 07 29 -4.51297753038541 > 57 08 02 -3.29214621128099 > 57 08 05 -6.44144118698436 > 57 08 16 -3.09828319252649 > 57 08 22 -1.58619965424759 > 57 08 28 -4.3638763217976 > 57 08 30 -9.93305062749234 > 57 09 03 -9.29132498735799 > 57 09 05 -8.50042591710719 > 57 09 13 -5.58659819576993 > 57 09 19 -2.4834966827445 > 57 09 23 -5.69132573226493 > 57 09 26 -3.6296458026372 > to the year 1994. The length is 1203. > > The plot is attached. How to label the x-axis is the problem. > I will send the code in another email so that the moderator will let it go. > Thank again > Best > Ogbos > > On Sun, Aug 25, 2019 at 3:58 PM Bert Gunter <bgunter.4...@gmail.com> wrote: > > > > I think you will need to show us your code to get useful help (see the > > posting guide). A small subset of your data and a reprex may also be > > needed. R has several different graphics systems, and the answer depends on > > which you are using. In base graphics, the axis() function (see its Help > > page) may be what you need, but that's just a guess. Others may be able to > > be more specific. > > > > Cheers, > > Bert > > > > Bert Gunter > > > > "The trouble with having an open mind is that people keep coming along and > > sticking things into it." > > -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) > > > > > > On Sat, Aug 24, 2019 at 9:49 PM Ogbos Okike <giftedlife2...@gmail.com> > > wrote: > >> > >> Dear Contributors, > >> I have two dataset of different lengths, each containing year, month, > >> day and counts. > >> After converting the date using as.Date function, I plotted the two > >> dateset on one graph. That was fine. > >> > >> I am, however, having problems with the axis 1 where I am trying to > >> put the dates. > >> Since the two dates are not exactly the same, I stored the first date > >> as x1 and the second x2. x1 runs from 1953-01-02 to 2006-11-15 while > >> the range of x2 is between 1957-07-26 and 1994-07-17. > >> I tired a number of approaches but the one that seems to come close to > >> what I am looking for is: axis.Date(1, at=seq(min(x1), max(x1), by="50 > >> mon"), format="%Y"). > >> > >> Unfortunately, it only labeled the axis up to 1994 even when I tried > >> to replace x1 with x2 in the code. Since one of the dates is up to > >> 2006, I I wish to display the minimum (1953) and maximum (2006) dates > >> and some possible intermediary dates. > >> I am always indebted to you. > >> > >> Best wishes > >> 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. ______________________________________________ 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.