Dear Rui Barradas,
Many thanks for your time.
Though I have not used ggplot before, I quickly installed all the required
packages and attempted your code. It produced the plot and solved the
initial problem of the wrong labels on x-axis.
But I noticed several other problems with the plot (some of which I have
circumvented by adjusting the script I displayed earlier). As I am not
quite familiar with ggplot, it was very difficult for me to interact with
your code. I decided to go back to my own code. I tried to use axis(1 ...,
labels=c()...) to assign Thu, Fri .... It worked for me before I got a
third response (very helpful) from the list.
Best regards
Ogbos
On Fri, May 1, 2020 at 6:35 PM Rui Barradas <ruipbarra...@sapo.pt> wrote:
> Hello,
>
> I have just noticed that the facets' labels are on the right, in your
> code example they are on the left. Change the following:
>
>
> facet_grid(rows = vars(variable), switch = "y")
>
>
> Hope this helps,
>
> Rui Barradas
>
> Às 18:29 de 01/05/20, Rui Barradas escreveu:
> > Hello,
> >
> > If you don't mind a ggplot2 solution, here it is.
> > As usual with ggplot, it's better to have the data in long format so I
> > also load packages dplyr and tidyr. The reshaping accounts for half the
> > code, the code for the graph itself is not very complicated.
> >
> >
> >
> > library(ggplot2)
> > library(scales)
> > library(dplyr)
> > library(tidyr)
> >
> > dta %>%
> > select(datetime, B, BX, BZ) %>%
> > rename(x = datetime, B1 = B, B2 = BX, B3 = BZ) %>%
> > pivot_longer(
> > cols = starts_with("B"),
> > names_to = "variable",
> > values_to = "value"
> > ) %>%
> > ggplot(aes(x, value)) +
> > geom_line() +
> > xlab("") + ylab("") +
> > scale_x_datetime(
> > breaks = seq(min(x), max(x), by = "days"),
> > labels = time_format("%a")) +
> > facet_grid(rows = vars(variable)) +
> > theme_bw()
> >
> >
> >
> > Hope this helps,
> >
> > Rui Barradas
> >
> > Às 17:46 de 01/05/20, Ogbos Okike escreveu:
> >> Dear Contributors,
> >> I am trying to do a plot with multiple y-axis on a common x-axis. I have
> >> made some progress but I am having difficulties with datetime
> >> appearing on
> >> x-axis. Instead of as date, it appears as large numbers.
> >>
> >> Sample of my data:
> >> 78 09 28 0 6.7 -40.4 -3.5 2.3 -3.6 278036 5.8 612
> >> 78 09 28 1 5.7 7.3 -4.4 1.4 0.6 261169 4.9 623
> >> 78 09 28 2 5.6 -20.9 -3.6 3.1 -1.8 269323 4.8 625
> >> 78 09 28 3 5.8 -46.9 -3.4 0.2 -3.6 254221 4.7 626
> >> 78 09 28 4 6.2 -41.4 -3.9 1.2 -3.6 248567 4.2 618
> >> 78 09 28 5 6.1 -54.6 -3 1 -4.4 252669 4.4 629
> >> 78 09 28 6 NA NA NA NA NA 220480 4.3 606
> >> 78 09 28 7 7.1 26 -4.8 3.6 2.9 216887 4.4 613
> >> 78 09 28 8 7.1 -21.4 -4.7 3.5 -2.3 249928 4 643
> >> 78 09 28 9 7 -13.6 -2.7 5.2 -1.4 268806 4 656
> >> 78 09 28 10 6.6 -18.5 -3.7 3.6 -1.7 285608 4.2 651
> >> 78 09 28 11 6.7 -6.6 -2.1 4.7 -0.6 280446 4.2 661
> >> 78 09 28 12 6.7 -29.7 -1.9 4 -2.5 281928 4.2 659
> >> 78 09 28 13 6.2 -8.7 -2.1 4.7 -0.8 273465 4.4 651
> >> 78 09 28 14 6.1 31.8 -4.7 0.9 3 212374 4.2 631
> >> 78 09 28 15 5.4 9.3 -4.5 -0.9 0.8 219662 3.9 636
> >> 78 09 28 16 4.9 6.3 -4.5 -0.5 0.5 220610 3.8 628
> >> 78 09 28 17 4.7 4.5 -4.3 -0.6 0.3 214432 3.7 628
> >> 78 09 28 18 4.9 3.1 -4.5 0.3 0.2 202199 3.5 623
> >> 78 09 28 19 5 -13 -4.6 -0.3 -1.1 192859 3.2 619
> >> 78 09 28 20 5.1 -26 -3.5 2 -2 193868 3.1 627
> >> 78 09 28 21 10.1 -37.5 -5.1 4.9 -5.4 284122 5.9 683
> >> 78 09 28 22 8.4 -7.3 -3.6 5.5 -0.8 367499 6.2 694
> >> 78 09 28 23 8.2 -17.7 -4.7 4.4 -2.1 346644 4.9 689
> >> 78 09 29 0 8 6.3 -4.3 5.8 0.8 269569 4.7 708
> >> 78 09 29 1 8.6 35 -3 6.2 4.8 187132 3.2 709
> >> 78 09 29 2 8.1 24.8 -4.1 6 3.3 166644 3.5 689
> >> 78 09 29 3 15.9 29.6 -2.4 9.6 5.6 720902 6.6 866
> >> 78 09 29 4 14.9 18.9 -2.8 13.8 4.8 1587324 8.1 912
> >> 78 09 29 5 14.1 8.2 -10.2 8.3 1.9 1509336 8.1 871
> >> 78 09 29 6 15.6 -2.9 -6.1 11.5 -0.7 968890 7.8 878
> >> 78 09 29 7 19.6 -49.7 1.2 12.1 -14.4 574978 5.4 906
> >> 78 09 29 8 23.5 -44.2 -1.4 11.7 -11.5 287721 4.1 865
> >> 78 09 29 9 25.3 -71.8 -2.4 7.5 -23.9 58521 5.1 859
> >> 78 09 29 10 24.8 -63.1 -2.3 10.9 -22 74956 1.5 821
> >> 78 09 29 11 23.8 -50 -2 15.1 -18.1 64510 3 807
> >> 78 09 29 12 21.4 -34.3 0.1 17.6 -12 50247 2.3 795
> >> 78 09 29 13 18.4 -20 -2 17 -6.3 59939 2.9 802
> >> 78 09 29 14 16.4 -13.9 -2.8 15.6 -3.9 56499 1.4 799
> >> 78 09 29 15 15.3 -7.9 -2 14.9 -2.1 49346 1.3 759
> >> 78 09 29 16 14.1 -0.5 -2.9 13.8 -0.1 57732 0.5 730
> >> 78 09 29 17 12.6 14.7 -3.1 11.5 3.1 44351 1.7 722
> >> 78 09 29 18 11.2 74.3 0.6 2.9 10.4 68224 0.4 712
> >> 78 09 29 19 10.5 60.5 3.2 4 9.1 48082 0.2 706
> >> 78 09 29 20 10 54.7 3.9 4.2 8.1 56404 0.2 710
> >> 78 09 29 21 9.3 61.1 3.4 2.9 8.1 68824 0.2 694
> >> 78 09 29 22 8.6 64.4 2.6 2.7 7.8 44579 0.2 683
> >> 78 09 29 23 7.9 63.3 2 3 7.1 39760 0.2 661
> >> 78 09 30 0 7.3 59.8 2.3 2.8 6.3 NA NA NA
> >> 78 09 30 1 6.6 65.6 1.7 2.1 6 57382 0.5 664
> >> 78 09 30 2 6 70.8 1.3 1.4 5.5 63540 0.5 654
> >> 78 09 30 3 4.4 45.1 1 2.7 2.9 60856 2.5 635
> >> 78 09 30 4 3.6 -28 -1 -0.3 -0.6 80967 1.9 633
> >> 78 09 30 5 4 -47.6 -2.2 -0.1 -2.4 48391 1.2 642
> >> 78 09 30 6 4.1 -28.9 -3.1 -1.4 -1.9 45012 1.1 643
> >> 78 09 30 7 4.2 -42.7 -2.7 1.2 -2.7 45562 1.1 633
> >> and part of my code is:
> >> Sys.setenv( TZ="GMT" )
> >> dta<-read.table("IFEDIFIG1D",col.names=c("year", "month", "day",
> >> "hour","B","LAT","BX","BZ","BY","SWT","SWD","SW"))
> >> dta$year <- with( dta, ifelse(year < 50, year + 2000, year + 1900))
> >> dta$datetime <- with( dta, as.POSIXct(ISOdatetime(year,
> >> month,day,hour,0,0)))
> >> x = dta$datetime
> >> B=dta$B
> >> LAT=dta$LAT
> >> BX=dta$BX
> >> BZ=dta$BZ
> >> par(mfcol = c(3, 1), mar = numeric(4), oma = c(4, 4, 1, 1),
> >> mgp = c(2, 0.5, 0))#0.5 here controls x-axis margins
> >> plot(x, B,type="l", axes = FALSE)
> >> abline(v=-3)
> >> axis(2L)
> >> plot(x, LAT,type="l", axes = FALSE)
> >> abline(v=-3)
> >> axis(2L)
> >> plot(x, BX, type="l",axes = FALSE)
> >> abline(v=-3)
> >> axis(1L)
> >> axis(2L)
> >> mtext("A1", side = 1, outer = TRUE, line = 2.2)
> >> mtext("B1", side = 2, outer = TRUE, line = 2.2,at =0.2)
> >> mtext("B2", side = 2, outer = TRUE, line = 2.2,at =0.5)
> >> mtext("B2", side = 2, outer = TRUE, line = 2.2,at =0.8)
> >>
> >> If I do a simple plot of x and B, the x-axis will be fine, appearing as
> >> Thu, Fri, Sat, Sun.
> >> Please let me know what I am doing wrong with the multiple plot code
> >> above.
> >> I want the same Thu, Fri, Sat, Sun to appear on the common x-axis.
> >>
> >> I am most grateful for your kind response.
> >>
> >> Warmest regards
> >> Ogbos
> >>
> >> [[alternative HTML version deleted]]
> >>
> >> ______________________________________________
> >> 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.
>
[[alternative HTML version deleted]]
______________________________________________
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.
