Many thanks, the solution is simple and elegant, as usual ;). R is a great tool, the best I know, unfortunately I need more time to practice ...
regards Tomas On Sat, 2009-09-12 at 22:52 -0400, Gabor Grothendieck wrote: > Try this: > > library(zoo) > x <- zoo(rnorm(1:6000), as.Date("1992-11-11")+c(1:6000)) > plot(x, xaxt = "n") > # cut time into years, remove duplicate years and convert back to Date yy <- as.Date(unique(cut(time(x), "year"))) > axis(1, yy, format(yy, "%y")) > > There are a number of examples of custom X axes in the > examples section of ?plot.zoo > > On Sat, Sep 12, 2009 at 10:33 PM, <lanc...@fns.uniba.sk> wrote: > > On Sun, 2009-07-26 at 13:17 +0200, Poersching wrote: > >> stvienna wiener schrieb: > >> > Hi all, > >> > > >> > I am plotting a financial time series, but I need a more detailed X-Axis. > >> > > >> > Example: > >> > x <- zoo(rnorm(1:6000), as.Date("1992-11-11")+c(1:6000)) > >> > plot(x) > >> > > >> > The X-Axis is labeled "1995", "2000" and "2005". > >> > I would need either "1995", "1997", etc. or maybe yearly > >> > > >> > > >> > I used google first, then look at "?plot.zoo" but could't get it working. > >> > > >> > > >> > Regards, > >> > Steve > >> > > >> > ______________________________________________ > >> > R-help@r-project.org mailing list > >> > 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. > >> > > >> > > >> Hey, > >> > >> try something like the following: > >> > >> plot(x, y, xaxt="n") > >> axis.Date(1, at=seq(as.Date("1960-01-01"), max(as.Date(x)), "years"), labels = FALSE, tcl = -0.2) > >> axis.Date(1, at=seq(as.Date("1960-01-01"), max(as.Date(x)), "5 years"), > >> labels = TRUE, las=3, tcl = -0.2) > >> > >> Regards, > >> Christian > > > > Hello, > > > > I'm actually struggling with the similar problem. I applied Your script > > into mine like this: > > - > > library(Rdbi) > > library(RdbiPgSQL) > > conn <- dbConnect(PgSQL(), host="localhost", dbname="BVS", > > user="postgres", password = "******") > > > > query_duj_kal <- dbSendQuery(conn, "select zdroj as well, cas as date, fe2, fe3, mn, nh4, no2, no3, o2, teplota as temperature from analystats > > where zdroj like 'Dunaj Kalinkovo' order by date") > > watchem_duj_kal <- dbGetResult(query_duj_kal) > > > > date <- (watchem_duj_kal$date) > > date_p <- strptime(date, "%m-%d-%Y") > > no3 <- (watchem_duj_kal$no3) > > nh4 <- (watchem_duj_kal$nh4) > > > > par(mfrow=c(3,1), ylog = TRUE, yaxp = c(0.01, 100, 3)) > > maxy <- 100 > > miny <- 0.005 > > plot(date_p, no3,log = "y", type = "l", ylim = c(miny, maxy), col = "darkred", main = "Dunaj Kalinkovo", xlab = "time", ylab = "log(NO3-,NH4 > > +)", xaxt="n") > > axis.Date(1, at=seq(as.Date("1972-01-01"),max(as.Date(date_p))), "years"), labels = FALSE, tcl = -0.2) > > axis.Date(1, at=seq(as.Date("1970-01-01"),max(as.Date(date_p))), "5 years"), labels = TRUE, las=3, tcl = -0.2) > > - > > > > Everything goes well, except that I cannot get no tickles nor labels at > > the moment and I cannot find what I did wrong. > > > > Many thanks for any advice > > > > Tomas > > > > ______________________________________________ > > R-help@r-project.org mailing list > > 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 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.