在R中,为什么我在POSIXct和POSIXlt之间得到一毫秒的差异?

时间:2017-09-12 14:39:39

标签: r posixct strptime posixlt

此代码段

options(digits.secs=3)
s<-"12:00:00.188"
fmt<-"%I:%M:%OS"
print(strptime(s,fmt))
print(as.POSIXct(strptime(s,fmt)))

给出了这个文本输出:

[1] "2017-09-12 00:00:00.188 CEST"
[1] "2017-09-12 00:00:00.187 CEST"

虽然我预计上述结果是相同的。我错过了什么?

我的会话信息:

print(sessionInfo())

给出:

R version 3.4.0 (2017-04-21)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows 7 x64 (build 7601) Service Pack 1

Matrix products: default

locale:
[1] LC_COLLATE=Italian_Italy.1252  LC_CTYPE=Italian_Italy.1252    LC_MONETARY=Italian_Italy.1252 LC_NUMERIC=C                  
[5] LC_TIME=Italian_Italy.1252    

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

loaded via a namespace (and not attached):
[1] compiler_3.4.0 tools_3.4.0   

Linux中的结果相同:

R version 3.3.3 (2017-03-06)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Debian GNU/Linux 9 (stretch)

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C               LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8     LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8   
 [7] 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] stats     graphics  grDevices utils     datasets  methods   base     

loaded via a namespace (and not attached):
[1] tools_3.3.3

修改(在Roland's comment之后)

也许我错了,但在我看来,0.002在浮点数上无法表示,仍然有0.002,POSIXct和POSIXlt之间没有区别:

options(digits.secs=3)
s<-"12:00:00.002"
fmt<-"%I:%M:%OS"
print(strptime(s,fmt))
print(as.POSIXct(strptime(s,fmt)))

给出:

[1] "2017-09-12 00:00:00.002 CEST"
[1] "2017-09-12 00:00:00.002 CEST"

1 个答案:

答案 0 :(得分:1)

您可以在日期时间的文档中阅读相关内容:https://stat.ethz.ch/R-manual/R-devel/library/base/html/DateTimeClasses.html

特别是:

  

类“POSIXct”表示自1970年初(UTC时区)以来的(带符号)秒数作为数字向量。

strptime和另一个班级POSIXlt以不同方式存储日期时间。

因此亚秒精度存在问题:

  

亚秒级准确度

     

类“POSIXct”和“POSIXlt”能够表示一秒的分数。 (两种形式之间的分数转换可能不准确,但精度要高于微秒。)

因此,您会看到POSIXltstrptime准确打印这些内容:

strptime(s,fmt)
as.POSIXlt(strptime(s,fmt), format = "%Y-%m-%d %H:%M:%OS")

但是,因为POSIXct进行计算以将日期表示为数字,所以由于浮点精度,闰秒等原因可能会出现不准确的情况。