Re: [Rd] as.Date function yields inconsistent results (PR#14166)

From: Tony Plate <tplate_at_acm.org>
Date: Mon, 28 Dec 2009 02:15:06 -0700

I think you're right that this is a timezone issue -- it seems to be a consequence of the behavior described about in the article by Gabor Grothendieck and Thomas Petzoldt: "R help desk: Date and time classes in R." R News, 4(1):29-32, June 2004.
<http://cran.r-project.org/doc/Rnews/Rnews_2004-1.pdf> http://cran.r-project.org/doc/Rnews/Rnews_2004-1.pdf

(there is a reference at
http://wiki.r-project.org/rwiki/doku.php?id=guides:times-dates)

BEGIN QUOTE from p31 http://cran.r-project.org/doc/Rnews/Rnews_2004-1.pdf

Regarding POSIX classes, the user should be aware of the following: ...
* POSIXlt. The tzone attribute on POSIXlt times are ignored so it is safer to use POSIXct than POSIXlt when performing arithmetic or other manipulations that may depend on time zones.

END QUOTE The ignoring of the tzone attribute appears to apply to as.Date() conversions too. The behavior is especially interesting because POSIXlt objects are converted to POSIXct objects by arithmetic, and the arithmetic operation preserves the tzone attribute. So, if you want as.Date() to pay attention to the tzone attribute of a POSIXlt object, you might be able to just add 0 to it (or just use as.POSIXct() on it). (Though I don't know if there are other pitfalls in this path to catch the unwary). The following demonstrates various aspects of the behavior:

> d.ct.CET <- as.POSIXct("1999-03-18", tz="CET")
> d.lt.CET <- as.POSIXlt("1999-03-18", tz="CET")
> d.ct.UTC <- as.POSIXct("1999-03-18", tz="UTC")
> d.lt.UTC <- as.POSIXlt("1999-03-18", tz="UTC")
> d.ct.EST <- as.POSIXct("1999-03-18", tz="EST")
> d.lt.EST <- as.POSIXlt("1999-03-18", tz="EST")
> # ii is hours to catch the date changes in the time-zones here
> ii <- c(0,1,19,20,23,24)*3600
> # all columns of x1 are as.Date() of POSIXct objects
> # each column follows its tzone attribute because as.Date() pays
> # attention to tzone attr of POSIXct object, and arithmetic converts
> # POSIXlt object to POSIXct object preserving the tzone attr
> (x1 <- data.frame(hour=ii, ct.UTC=as.Date(d.ct.UTC + ii),

lt.UTC=as.Date(d.lt.UTC + ii), ct.CET=as.Date(d.ct.CET + ii), 
lt.CET=as.Date(d.lt.CET + ii), ct.EST=as.Date(d.ct.EST + ii), 
lt.EST=as.Date(d.lt.EST + ii)))
   hour     ct.UTC     lt.UTC     ct.CET     lt.CET     ct.EST     lt.EST
1     0 1999-03-18 1999-03-18 1999-03-17 1999-03-17 1999-03-18 1999-03-18
2  3600 1999-03-18 1999-03-18 1999-03-18 1999-03-18 1999-03-18 1999-03-18
3 68400 1999-03-18 1999-03-18 1999-03-18 1999-03-18 1999-03-19 1999-03-19
4 72000 1999-03-18 1999-03-18 1999-03-18 1999-03-18 1999-03-19 1999-03-19
5 82800 1999-03-18 1999-03-18 1999-03-18 1999-03-18 1999-03-19 1999-03-19 6 86400 1999-03-19 1999-03-19 1999-03-18 1999-03-18 1999-03-19 1999-03-19
> all.equal(x1$lt.UTC, x1$ct.UTC)

[1] TRUE
> all.equal(x1$lt.CET, x1$ct.CET)

[1] TRUE
> all.equal(x1$lt.EST, x1$ct.EST)

[1] TRUE
> class(d.lt.EST)

[1] "POSIXt" "POSIXlt"
> class(d.lt.EST + ii)

[1] "POSIXt" "POSIXct"
> as.POSIXlt(d.lt.EST + ii)
[1] "1999-03-18 00:00:00 EST" "1999-03-18 01:00:00 EST"
[3] "1999-03-18 19:00:00 EST" "1999-03-18 20:00:00 EST"
[5] "1999-03-18 23:00:00 EST" "1999-03-19 00:00:00 EST"
 >

> # the lt.* columns of x2 are as.Date() of POSIXlt objects, and these
> # are all the same the ct.UTC column because as.Date() ignores the
> # tzone attribute of POSIXlt objects
> (x2 <- data.frame(hour=ii, ct.UTC=as.Date(d.ct.UTC + ii),
lt.UTC=as.Date(as.POSIXlt(d.lt.UTC + ii)), ct.CET=as.Date(d.ct.CET + ii), lt.CET=as.Date(as.POSIXlt(d.lt.CET + ii)), ct.EST=as.Date(d.ct.EST + ii), lt.EST=as.Date(as.POSIXlt(d.lt.EST + ii))))
   hour     ct.UTC     lt.UTC     ct.CET     lt.CET     ct.EST     lt.EST
1     0 1999-03-18 1999-03-18 1999-03-17 1999-03-18 1999-03-18 1999-03-18
2  3600 1999-03-18 1999-03-18 1999-03-18 1999-03-18 1999-03-18 1999-03-18
3 68400 1999-03-18 1999-03-18 1999-03-18 1999-03-18 1999-03-19 1999-03-18
4 72000 1999-03-18 1999-03-18 1999-03-18 1999-03-18 1999-03-19 1999-03-18
5 82800 1999-03-18 1999-03-18 1999-03-18 1999-03-18 1999-03-19 1999-03-18 6 86400 1999-03-19 1999-03-19 1999-03-18 1999-03-19 1999-03-19 1999-03-19
> all.equal(x2$lt.UTC, x2$ct.UTC)

[1] TRUE
> all.equal(x2$lt.CET, x2$ct.UTC)

[1] TRUE
> all.equal(x2$lt.EST, x2$ct.UTC)

[1] TRUE
 >
> sessionInfo()

R version 2.10.1 (2009-12-14)
i486-pc-linux-gnu

locale:

 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C             
 [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8   
 [5] LC_MONETARY=C              LC_MESSAGES=en_US.UTF-8  
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                
 [9] LC_ADDRESS=C               LC_TELEPHONE=C           
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C      

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

loaded via a namespace (and not attached): [1] tools_2.10.1
 >
 >

diverses_at_univecom.ch wrote:

> Full_Name: Mario Luoni
> Version: 2.10.0
> OS: Windows XP HE SP3
> Submission from: (NULL) (217.194.59.134)
>
>
> This piece of code:
>
> zzz1 <- as.POSIXct("1999-03-18", tz="CET")
> zzz2 <- as.POSIXlt("1999-03-18", tz="CET")
> zzz1 == zzz2
> as.Date(zzz1)
> as.Date(zzz2)
>
> yields TRUE for "zzz1==zzz2", but the two dates returned by as.Date are 
> different:
>
>   
>> as.Date(zzz1)
>>     
> [1] "1999-03-17"
>   
>> as.Date(zzz2)
>>     
> [1] "1999-03-18"
>
> For me this looks like a bug, even though it could be a problem with timezones,
> but I couldn't find documentation that would explain that behaviour.
>
>
>   
>> sessionInfo()
>>     
> R version 2.10.0 (2009-10-26) 
> i386-pc-mingw32 
>
> locale:
> [1] LC_COLLATE=English_United States.1252  LC_CTYPE=English_United States.1252  
>
> [3] LC_MONETARY=English_United States.1252 LC_NUMERIC=C                         
>
> [5] LC_TIME=English_United States.1252    
>
> attached base packages:
> [1] grDevices datasets  splines   graphics  stats     tcltk     utils    
> methods   base     
>
> other attached packages:
> [1] svSocket_0.9-48 TinnR_1.0.3     R2HTML_1.59-1   Hmisc_3.7-0    
> survival_2.35-7
>
> loaded via a namespace (and not attached):
> [1] cluster_1.12.1  grid_2.10.0     lattice_0.17-26 svMisc_0.9-56   tools_2.10.0
>
> ______________________________________________
> R-devel_at_r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-devel
>
>

______________________________________________
R-devel_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel Received on Mon 28 Dec 2009 - 09:25:57 GMT

This archive was generated by hypermail 2.2.0 : Mon 28 Dec 2009 - 12:51:15 GMT