[R] Datetime in doubly y-axis plot

Ogbos Okike g||ted|||e2014 @end|ng |rom gm@||@com
Thu Nov 21 09:49:06 CET 2019


Dear Members,
I have some hourly data. Usin:
dta$year <- with( dta, ifelse(year < 50, year + 2000, year + 1900))
dta$datetime <- with( dta, as.POSIXct(ISOdatetime(year, month,day,hour,0,0)))
I converted the hour to time format and stored in Year.
The data consists of two different observations and I wish to plot
same on one graph (two y-axis and common x-axis).

Part of my script is of the form:
library(plotrix)
plot(Year,Li,pch=16,axes=FALSE,xlab="",ylab="",type="l",col="black",ylim=c(1,22909))
axis(2, ylim=c(1,22909),col="black",las=1)
mtext("Lightning count/day", side=2, line=4)
box()
par(new=TRUE)
plot(Year,CR, pch=15,
xlab="",ylab="",ylim=c(7302,9983),axes=FALSE,type="l",col="red")
mtext("GCR count/day",side=4,col="red",line=3)
axis(4, ylim=c(7302,9983), col="red",col.axis="red",las=1)
The two y-axis worked and looked fine.

But I have difficulties with x-axis.
I tried things like in an attempt to draw the x-axis.
## Draw the time axis
#axis(1,pretty(range(Year),10))
axis(side=1,at=c(as.datetime("2005-01-01 01:00:00
GMT"),as.datetime("2007-01-01 01:00:00 GMT"),as.datetime("2009-01-01
01:00:00 GMT"),labels=c("2005","2007","2009"))
#axis(side=1, at =Year)
#axis(1,Year)

But none worked.

However, when I plotted the two graphs separately, the x-axis was
correctly labelled from 2005 to 2010 (the time range of my data).

Small portion of the data is:
04 12 31 10  8637 4992 0.310676228913214 49.7113438132607
04 12 31 11  8634 4183 0.275834035016405 25.4492290005748
04 12 31 12  8652 4160 0.484887198397259 24.7594531777172
04 12 31 13  8626 3411 0.182921517958247 2.2967535550946
04 12 31 15  8618 3904 0.09000900090009 17.0819483667808
04 12 31 16  8620 3428 0.113237130164629 2.80658785894585
04 12 31 17  8603 3041 -0.0842019685839552 -8.79964011696198
04 12 31 18  8579 2340 -0.362939519758427 -29.8228075875341
04 12 31 19  8588 3006 -0.258412938068 -9.8492989778322
04 12 31 20  8604 3280 -0.0725879039516855 -1.63196960987679
04 12 31 21  8568 3696 -0.490694230713394 10.8439757078949
04 12 31 22  8567 2307 -0.502308295345664 -30.8124859420688
04 12 31 23  8570 2726 -0.467466101448855 -18.2465698647939
05 01 01 00  8579 2879 -0.568878758318289 8.46755203918244
05 01 01 01  8581 1586 -0.545698639133843 -40.2467740416313
05 01 01 02  8562 1183 -0.765909771386082 -55.4299708015447
05 01 01 03  8602 1491 -0.302307387697159 -43.8259395309409
05 01 01 04  8576 1742 -0.603648937094958 -34.369407553923
05 01 01 05  8622 2001 -0.0705061958526974 -24.6114721672789
05 01 01 06  8613 2675 -0.174816732182705 0.781765093717623
05 01 01 07  8622 3293 -0.0705061958526974 24.0651784873316
05 01 01 08  8628 4333 -0.00096583829935895 63.2476217387209
05 01 01 09  8669 16 0.474226604981787 -99.3971931807479
05 01 01 10  8698 3671 0.810338333156256 38.3064895921635
05 01 01 11  8671 3856 0.497406724166233 45.2764434397664
05 01 01 12  8679 3411 0.590127200904018 28.5108787793162
05 01 01 13  8647 3428 0.219245293952879 29.1513610247716
05 01 01 14  8651 4596 0.265605532321772 73.156258830178
05 01 01 15  8638 3154 0.114934757622872 18.8282942450786
05 01 01 16  8635 2365 0.0801645788462025 -10.8976170292926
05 01 01 17  8628 2697 -0.00096583829935895 1.61062447018932
05 01 01 18  8653 3231 0.288785651506218 21.7293020627296
05 01 01 19  8657 3548 0.33514588987511 33.6724121691627
05 01 01 20  8650 2645 0.254015472729549 -0.348497692380145
05 01 01 21  8605 2153 -0.26753720892049 -18.8848073843835
05 01 01 22  8600 1759 -0.325487506881605 -33.7289253084676
Where the first four columns are year, month, day and hour. The last
four columns are the parameters of interest which I am trying to plot,
taking two at a time.
I will remain indebted for your assistance.
Best Regards
Ogbos



More information about the R-help mailing list