[Rd] problem with cut.Date/date plotting in ggplot2

From: Ben Bolker <bolker_at_zoology.ufl.edu>
Date: Fri, 10 Oct 2008 22:15:49 -0400

   I've hit a problem in ggplot2 which I can trace back to cut.Date , which is either a bug or (??) ggplot2 trying to do something it shouldn't with cut.Date (although its use of cut.Date seems OK).

  Apparently any (?) call of the form

cut(as.Date("2008-07-07"),"weeks")

where the date *begins the week*, gives the error

Error in 1:(1 + max(which(breaks < maxx))) :   result would be too long a vector
In addition: Warning message:
In max(which(breaks < maxx)) :
  no non-missing arguments to max; returning -Inf

  In particular, try:

d <- seq.Date(as.Date("2008-07-05"),as.Date("2008-07-08"),by="days") weekdays(d)

cut(d[1],"weeks") ## OK
cut(d[2],"weeks") ## OK
cut(d[3],"weeks") ## Monday -- breaks
cut(d[4],"weeks") ## OK

cut(d[1],"weeks",start.on.monday=FALSE) ## OK
cut(d[2],"weeks",start.on.monday=FALSE) ## Sunday -- breaks
cut(d[3],"weeks",start.on.monday=FALSE) ## OK cut(d[4],"weeks",start.on.monday=FALSE) ## OK

  The particular line that fails within cut.Date is

breaks <- breaks[1:(1 + max(which(breaks < maxx)))]

  but I haven't figured out enough about what's going on to be able to suggest a fix ...

  cheers
    Ben Bolker

session info:

> sessionInfo()

R version 2.7.2 (2008-08-25)
i486-pc-linux-gnu

locale:
LC_CTYPE=en_US.UTF-8;LC_NUMERIC=C;LC_TIME=en_US.UTF-8;LC_COLLATE=en_US.UTF-8;LC_MONETARY=C;LC_MESSAGES=en_US.UTF-8;LC_PAPER=en_US.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=en_US.UTF-8;LC_IDENTIFICATION=C

attached base packages:
[1] splines grid stats graphics grDevices utils datasets [8] methods base

other attached packages:
[1] ggplot2_0.7 MASS_7.2-44 RColorBrewer_1.0-2 proto_0.3-8

[5] reshape_0.8.1 plyr_0.1

 Here's the context in which I encountered the problem:

library(ggplot2)
tdat <- data.frame(date=seq.Date(as.Date("2008-07-07"),

                     as.Date("2008-10-10"),by="day"),
                   x=1:96)

q2 = qplot(date,x,data=tdat,geom="smooth") print(q2+geom_point())

  I can track the problem down to a call like

floor_date(as.Date("2008-07-07"),"weeks")

and digging further in, I can narrow it down to the issue above.



R-devel_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel Received on Sat 11 Oct 2008 - 14:57:46 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Sat 11 Oct 2008 - 19:31:20 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-devel. Please read the posting guide before posting to the list.

list of date sections of archive