[R] Correct x-axis of two in one graph

Discussion:

Ogbos Okike

2018-11-29 13:33:15 UTC

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-***@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.

PIKAL Petr

2018-11-29 14:03:15 UTC

Permalink

Hi

If I understand correctly you want Li and CR appear in one plot with the same x axis. Although it is not usually recommended you could use twoord.plot from plotrix or undocumented code below.

plot.yy <- function (x, yright, yleft, yleftlim = NULL, yrightlim = NULL, xlab = NULL, yylab = list(NA, NA), pch = c(1, 2), col = c(1, 2), linky = F, smooth = 0, lwds = 1, length = 10, format = "%d/%m", rect = NULL, type = "p", ...)
{
par(mar = c(5, 4, 4, 2), oma = c(0, 0, 0, 3))
plot(x, yright, ylim = yrightlim, axes = F, ylab = "", xlab = xlab, pch = pch[1], col = col[1], type = type, ...)
if (!is.null(rect))
rect(x[rect[1]], rect[2], x[rect[3]], rect[4], col = "grey")
points(x, yright, ylim = yrightlim, ylab = "", xlab = xlab, pch = pch[1], col = col[1], ...)
axis(4, pretty(range(yright, na.rm = T), 10), col = col[1])
if (linky)
lines(x, yright, col = col[1], ...)
if (smooth != 0)
lines(supsmu(x, yright, span = smooth), col = col[1], lwd = lwds, ...)
if (is.na(yylab[[1]]))
mtext(deparse(substitute(yright)), side = 4, outer = T, line = 1, col = col[1], ...) else mtext(yylab[[1]], side = 4, outer = T, line = 1, col = col[1], ...)
par(new = T)
plot(x, yleft, ylim = yleftlim, ylab = "", axes = F, xlab = xlab,
pch = pch[2], col = col[2], ...)
box()
axis(2, pretty(range(yleft, na.rm = T), 10), col = col[2],
col.axis = col[2])
if (!inherits(x, c("Date", "POSIXt")))
axis(1, pretty(range(x, na.rm = T), 10))
else {
l <- length(x)
axis(1, at = x[seq(1, l, length = length)], labels = format(as.POSIXct(x[seq(1, l, length = length)]), format = format))
}
if (is.na(yylab[[2]]))
mtext(deparse(substitute(yleft)), side = 2, line = 2, col = col[2], ...)
else mtext(yylab[[2]], side = 2, line = 2, col = col[2], ...)
if (linky)
lines(x, yleft, col = col[2], lty = 2, ...)
if (smooth != 0)
lines(supsmu(x, yleft, span = smooth), col = col[2], lty = 2, lwd = lwds, ...)
}

something like

plot.yy(Year, Li, CR)

Cheers
Petr

-----Original Message-----
Sent: Thursday, November 29, 2018 2:33 PM
Subject: [R] Correct x-axis of two in one graph
Dear Contributors,
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.
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
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
______________________________________________
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.

Osobní údaje: Informace o zpracování a ochraně osobních údajů obchodních partnerů PRECHEZA a.s. jsou zveřejněny na: https://www.precheza.cz/zasady-ochrany-osobnich-udaju/ | Information about processing and protection of business partner’s personal data are available on website: https://www.precheza.cz/en/personal-data-protection-principles/
Důvěrnost: Tento e-mail a jakékoliv k němu připojené dokumenty jsou důvěrné a podléhají tomuto právně závaznému prohláąení o vyloučení odpovědnosti: https://www.precheza.cz/01-dovetek/ | This email and any documents attached to it may be confidential and are subject to the legally binding disclaimer: https://www.precheza.cz/en/01-disclaimer/

______________________________________________
R-***@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, reproducib

Ogbos Okike

2018-11-29 16:17:18 UTC

Permalink

Dear Petr,

Thank you so much for your contribution.

Let me show you what I have please.

I am attaching two plots. The first (twoineone) is my own while the
second (testrun) is plotted with the code you send.

The first is closer to my interest. The major problem I have is the
date. Those were set arbitrary with:
axis(side=1,at=c(1,400,800,1200,1500),labels=c("2004","2005","2006","2007","2008"))
and they seem not to be correct.

The second one does not show date and date is very important to me here.

Is there a way your code can implement the correct date on x-axis? I
would be very glad for further assistance.

Thank you again.
Best wishes
Ogbos

Post by PIKAL Petr
Hi
If I understand correctly you want Li and CR appear in one plot with the same x axis. Although it is not usually recommended you could use twoord.plot from plotrix or undocumented code below.
plot.yy <- function (x, yright, yleft, yleftlim = NULL, yrightlim = NULL, xlab = NULL, yylab = list(NA, NA), pch = c(1, 2), col = c(1, 2), linky = F, smooth = 0, lwds = 1, length = 10, format = "%d/%m", rect = NULL, type = "p", ...)
{
par(mar = c(5, 4, 4, 2), oma = c(0, 0, 0, 3))
plot(x, yright, ylim = yrightlim, axes = F, ylab = "", xlab = xlab, pch = pch[1], col = col[1], type = type, ...)
if (!is.null(rect))
rect(x[rect[1]], rect[2], x[rect[3]], rect[4], col = "grey")
points(x, yright, ylim = yrightlim, ylab = "", xlab = xlab, pch = pch[1], col = col[1], ...)
axis(4, pretty(range(yright, na.rm = T), 10), col = col[1])
if (linky)
lines(x, yright, col = col[1], ...)
if (smooth != 0)
lines(supsmu(x, yright, span = smooth), col = col[1], lwd = lwds, ...)
if (is.na(yylab[[1]]))
mtext(deparse(substitute(yright)), side = 4, outer = T, line = 1, col = col[1], ...) else mtext(yylab[[1]], side = 4, outer = T, line = 1, col = col[1], ...)
par(new = T)
plot(x, yleft, ylim = yleftlim, ylab = "", axes = F, xlab = xlab,
pch = pch[2], col = col[2], ...)
box()
axis(2, pretty(range(yleft, na.rm = T), 10), col = col[2],
col.axis = col[2])
if (!inherits(x, c("Date", "POSIXt")))
axis(1, pretty(range(x, na.rm = T), 10))
else {
l <- length(x)
axis(1, at = x[seq(1, l, length = length)], labels = format(as.POSIXct(x[seq(1, l, length = length)]), format = format))
}
if (is.na(yylab[[2]]))
mtext(deparse(substitute(yleft)), side = 2, line = 2, col = col[2], ...)
else mtext(yylab[[2]], side = 2, line = 2, col = col[2], ...)
if (linky)
lines(x, yleft, col = col[2], lty = 2, ...)
if (smooth != 0)
lines(supsmu(x, yleft, span = smooth), col = col[2], lty = 2, lwd = lwds, ...)
}
something like
plot.yy(Year, Li, CR)
Cheers
Petr

OsobnÃ Ãºdaje: Informace o zpracovÃ¡nÃ a ochranÄ osobnÃch ÃºdajÅ¯ obchodnÃch partnerÅ¯ PRECHEZA a.s. jsou zveÅejnÄny na: https://www.precheza.cz/zasady-ochrany-osobnich-udaju/ | Information about processing and protection of business partnerâs personal data are available on website: https://www.precheza.cz/en/personal-data-protection-principles/
DÅ¯vÄrnost: Tento e-mail a jakÃ©koliv k nÄmu pÅipojenÃ© dokumenty jsou dÅ¯vÄrnÃ© a podlÃ©hajÃ tomuto prÃ¡vnÄ zÃ¡vaznÃ©mu prohlÃ¡ÄenÃ o vylouÄenÃ odpovÄdnosti: https://www.precheza.cz/01-dovetek/ | This email and any documents attached to it may be confidential and are subject to the legally binding disclaimer: https://www.precheza.cz/en/01-disclaimer/

PIKAL Petr

2018-11-30 08:24:25 UTC

Permalink

Hi

You cannot expect any code to give you results precisely according to your wishes just out of the box.

You could modify x axis by changing format parameter. It is preset to = "%d/%m", so you could change it to "%Y", if you want to display only year.

format parameter goes to this line in my code.

axis(1, at = x[seq(1, l, length = length)], labels = format(as.POSIXct(x[seq(1, l, length = length)]), format = format))

The code actually has several other parameters which you could modify or you are free to modify the code itself.

Cheers
Petr

-----Original Message-----
Sent: Thursday, November 29, 2018 5:17 PM
Subject: Re: [R] Correct x-axis of two in one graph
Dear Petr,
Thank you so much for your contribution.
Let me show you what I have please.
I am attaching two plots. The first (twoineone) is my own while the second
(testrun) is plotted with the code you send.
The first is closer to my interest. The major problem I have is the date. Those
axis(side=1,at=c(1,400,800,1200,1500),labels=c("2004","2005","2006","2007",
"2008"))
and they seem not to be correct.
The second one does not show date and date is very important to me here.
Is there a way your code can implement the correct date on x-axis? I would be
very glad for further assistance.
Thank you again.
Best wishes
Ogbos

Post by PIKAL Petr
Hi
If I understand correctly you want Li and CR appear in one plot with the same

x axis. Although it is not usually recommended you could use twoord.plot from
plotrix or undocumented code below.

Post by PIKAL Petr
plot.yy <- function (x, yright, yleft, yleftlim = NULL, yrightlim = NULL, xlab =

NULL, yylab = list(NA, NA), pch = c(1, 2), col = c(1, 2), linky = F, smooth = 0,
lwds = 1, length = 10, format = "%d/%m", rect = NULL, type = "p", ...)

Post by PIKAL Petr
{
par(mar = c(5, 4, 4, 2), oma = c(0, 0, 0, 3))
plot(x, yright, ylim = yrightlim, axes = F, ylab = "", xlab = xlab, pch = pch[1],

col = col[1], type = type, ...)

Post by PIKAL Petr
if (!is.null(rect))
rect(x[rect[1]], rect[2], x[rect[3]], rect[4], col = "grey")
points(x, yright, ylim = yrightlim, ylab = "", xlab = xlab, pch = pch[1], col =

col[1], ...)

Post by PIKAL Petr
axis(4, pretty(range(yright, na.rm = T), 10), col = col[1])
if (linky)
lines(x, yright, col = col[1], ...)
if (smooth != 0)
lines(supsmu(x, yright, span = smooth), col = col[1], lwd = lwds, ...)
if (is.na(yylab[[1]]))
mtext(deparse(substitute(yright)), side = 4, outer = T, line = 1, col =

col[1], ...) else mtext(yylab[[1]], side = 4, outer = T, line = 1, col = col[1], ...)

Post by PIKAL Petr
par(new = T)
plot(x, yleft, ylim = yleftlim, ylab = "", axes = F, xlab = xlab,
pch = pch[2], col = col[2], ...)
box()
axis(2, pretty(range(yleft, na.rm = T), 10), col = col[2],
col.axis = col[2])
if (!inherits(x, c("Date", "POSIXt")))
axis(1, pretty(range(x, na.rm = T), 10))
else {
l <- length(x)
axis(1, at = x[seq(1, l, length = length)], labels =

format(as.POSIXct(x[seq(1, l, length = length)]), format = format))

Post by PIKAL Petr
}
if (is.na(yylab[[2]]))
mtext(deparse(substitute(yleft)), side = 2, line = 2, col = col[2], ...)
else mtext(yylab[[2]], side = 2, line = 2, col = col[2], ...)
if (linky)
lines(x, yleft, col = col[2], lty = 2, ...)
if (smooth != 0)
lines(supsmu(x, yleft, span = smooth), col = col[2], lty = 2, lwd = lwds,

...)

Post by PIKAL Petr
}
something like
plot.yy(Year, Li, CR)
Cheers
Petr

data.

Post by PIKAL Petr

The data range is 2004/8/10 to 2008/8/22.
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.
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
______________________________________________
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.

Osobní údaje: Informace o zpracování a ochraně osobních údajů
https://www.precheza.cz/zasady-ochrany-osobnich-udaju/ | Information
about processing and protection of business partner’s personal data
https://www.precheza.cz/en/personal-data-protection-principles/
Důvěrnost: Tento e-mail a jakékoliv k němu připojené dokumenty jsou
důvěrné a podléhají tomuto právně závaznému prohláąení o vyloučení
odpovědnosti: https://www.precheza.cz/01-dovetek/ | This email and any
documents attached to it may be confidential and are subject to the
legally binding disclaimer: https://www.precheza.cz/en/01-disclaimer/

Ogbos Okike

2018-11-30 13:15:41 UTC

Permalink

Dear Petr,

Great!!! It worked. The years are interestingly displayed.

Please one thing more. I want all the data points represented by lines
only and not points. I tried to change type = type to type="l". But it
did work.

Thanks for more assistance.

Ogbos

Post by PIKAL Petr
Hi
You cannot expect any code to give you results precisely according to your wishes just out of the box.
You could modify x axis by changing format parameter. It is preset to = "%d/%m", so you could change it to "%Y", if you want to display only year.
format parameter goes to this line in my code.
axis(1, at = x[seq(1, l, length = length)], labels = format(as.POSIXct(x[seq(1, l, length = length)]), format = format))
The code actually has several other parameters which you could modify or you are free to modify the code itself.
Cheers
Petr

Post by PIKAL Petr
Hi
If I understand correctly you want Li and CR appear in one plot with the same

x axis. Although it is not usually recommended you could use twoord.plot from
plotrix or undocumented code below.

Post by PIKAL Petr
plot.yy <- function (x, yright, yleft, yleftlim = NULL, yrightlim = NULL, xlab =

NULL, yylab = list(NA, NA), pch = c(1, 2), col = c(1, 2), linky = F, smooth = 0,
lwds = 1, length = 10, format = "%d/%m", rect = NULL, type = "p", ...)

Post by PIKAL Petr
{
par(mar = c(5, 4, 4, 2), oma = c(0, 0, 0, 3))
plot(x, yright, ylim = yrightlim, axes = F, ylab = "", xlab = xlab, pch = pch[1],

col = col[1], type = type, ...)

Post by PIKAL Petr
if (!is.null(rect))
rect(x[rect[1]], rect[2], x[rect[3]], rect[4], col = "grey")
points(x, yright, ylim = yrightlim, ylab = "", xlab = xlab, pch = pch[1], col =

col[1], ...)

col[1], ...) else mtext(yylab[[1]], side = 4, outer = T, line = 1, col = col[1], ...)

format(as.POSIXct(x[seq(1, l, length = length)]), format = format))

...)

Post by PIKAL Petr
}
something like
plot.yy(Year, Li, CR)
Cheers
Petr

data.

Post by PIKAL Petr

______________________________________________
R-***@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.

PIKAL Petr

2018-11-30 13:44:06 UTC

Permalink

Hi

If you wanted only lines you need to set parameter linky to TRUE (sorry, the code is around 15 year old and not mentioned to distribution) and it was extended rather chaotically during years.

so

plot.yy(..., linky=TRUE) add lines
plot.yy(..., linky=TRUE, pch=c(NA,NA)) add lines and remove points
plot.yy(..., linky=TRUE, pch=c(NA,NA), smooth = some number) add lines, remove points and add smoothed lines (smooth is parameter to supsmu smoother)

HTH.
Cheers
Petr

-----Original Message-----
Sent: Friday, November 30, 2018 2:16 PM
Subject: Re: [R] Correct x-axis of two in one graph
Dear Petr,
Great!!! It worked. The years are interestingly displayed.
Please one thing more. I want all the data points represented by lines only and
not points. I tried to change type = type to type="l". But it did work.
Thanks for more assistance.
Ogbos

Post by PIKAL Petr
Hi
You cannot expect any code to give you results precisely according to your

wishes just out of the box.

Post by PIKAL Petr
You could modify x axis by changing format parameter. It is preset to =

"%d/%m", so you could change it to "%Y", if you want to display only year.

Post by PIKAL Petr
format parameter goes to this line in my code.
axis(1, at = x[seq(1, l, length = length)], labels =
format(as.POSIXct(x[seq(1, l, length = length)]), format = format))
The code actually has several other parameters which you could modify or

you are free to modify the code itself.

Post by PIKAL Petr
Cheers
Petr

-----Original Message-----
Sent: Thursday, November 29, 2018 5:17 PM
Subject: Re: [R] Correct x-axis of two in one graph
Dear Petr,
Thank you so much for your contribution.
Let me show you what I have please.
I am attaching two plots. The first (twoineone) is my own while the second
(testrun) is plotted with the code you send.
The first is closer to my interest. The major problem I have is the
axis(side=1,at=c(1,400,800,1200,1500),labels=c("2004","2005","2006",
"2007",
"2008"))
and they seem not to be correct.
The second one does not show date and date is very important to me here.
Is there a way your code can implement the correct date on x-axis? I
would be very glad for further assistance.
Thank you again.
Best wishes
Ogbos

Post by PIKAL Petr
Hi
If I understand correctly you want Li and CR appear in one plot with the same

x axis. Although it is not usually recommended you could use
twoord.plot from plotrix or undocumented code below.

Post by PIKAL Petr
plot.yy <- function (x, yright, yleft, yleftlim = NULL, yrightlim = NULL, xlab

Post by PIKAL Petr

NULL, yylab = list(NA, NA), pch = c(1, 2), col = c(1, 2), linky = F, smooth = 0,
lwds = 1, length = 10, format = "%d/%m", rect = NULL, type = "p", ...)

Post by PIKAL Petr
{
par(mar = c(5, 4, 4, 2), oma = c(0, 0, 0, 3))
plot(x, yright, ylim = yrightlim, axes = F, ylab = "", xlab =
xlab, pch = pch[1],

col = col[1], type = type, ...)

Post by PIKAL Petr
if (!is.null(rect))
rect(x[rect[1]], rect[2], x[rect[3]], rect[4], col = "grey")
points(x, yright, ylim = yrightlim, ylab = "", xlab = xlab,
pch = pch[1], col =

col[1], ...)

col[1], ...) else mtext(yylab[[1]], side = 4, outer = T, line = 1, col = col[1],

...)

Post by PIKAL Petr

format(as.POSIXct(x[seq(1, l, length = length)]), format = format))

lwds,

Post by PIKAL Petr

...)

Post by PIKAL Petr
}
something like
plot.yy(Year, Li, CR)
Cheers
Petr

data.

Post by PIKAL Petr

The data range is 2004/8/10 to 2008/8/22.
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","WWL
LN")
)
dev.off()
I got the correct axes except that the x-axis is not what I want.
axis(side=1,at=c(1,400,800,1200,1500),labels=c("2004","2005","20
06",
"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
______________________________________________
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.

Osobní údaje: Informace o zpracování a ochraně osobních údajů
https://www.precheza.cz/zasady-ochrany-osobnich-udaju/ |
Information about processing and protection of business partner’s
https://www.precheza.cz/en/personal-data-protection-principles/
Důvěrnost: Tento e-mail a jakékoliv k němu připojené dokumenty
jsou důvěrné a podléhají tomuto právně závaznému prohláąení o
vyloučení
odpovědnosti: https://www.precheza.cz/01-dovetek/ | This email and
any documents attached to it may be confidential and are subject
https://www.precheza.cz/en/01-disclaimer/