r - R:xts 时间戳与真实数据时间戳相差 1 毫秒
问题描述
所以我有以下数据。
tt <- structure(list(Timestamp = c("2018-03-01 09:51:59.969", "2018-03-01 09:51:59.969",
"2018-03-01 09:51:59.970", "2018-03-01 09:51:59.971", "2018-03-01 09:51:59.987",
"2018-03-01 09:51:59.988"), Mid_Px = c(30755.5, 30755, 30755.5,
30756, 30756.5, 30756.5)), .Names = c("Timestamp", "Mid_Px"), class = "data.frame", row.names = 85774:85779)
看起来像这样:
Timestamp Mid_Px
85774 2018-03-01 09:51:59.969 30755.5
85775 2018-03-01 09:51:59.969 30755.0
85776 2018-03-01 09:51:59.970 30755.5
85777 2018-03-01 09:51:59.971 30756.0
85778 2018-03-01 09:51:59.987 30756.5
85779 2018-03-01 09:51:59.988 30756.5
当我尝试使用下面的代码从中创建一个 xts 对象时,事情开始变得糟糕。
tt_ts <- strptime(tt[,1],"%Y-%m-%d %H:%M:%OS")
tt_ts
[1] "2018-03-01 09:51:59.969 CST" "2018-03-01 09:51:59.969 CST" "2018-03-01 09:51:59.970 CST" "2018-03-01 09:51:59.971 CST" "2018-03-01 09:51:59.987 CST"
[6] "2018-03-01 09:51:59.988 CST"
xts(x=tt[,c(-1)], order.by=tt_ts)
[,1]
2018-03-01 09:51:59.969 30755.5
2018-03-01 09:51:59.969 30755.0
2018-03-01 09:51:59.970 30755.5
2018-03-01 09:51:59.970 30756.0
2018-03-01 09:51:59.986 30756.5
2018-03-01 09:51:59.987 30756.5
请注意,第 4,5 和 6 行中的毫秒数不正确。
我在这里做错了什么?如何修复它以显示正确的时间戳?
解决方案
这类似于带有舍入毫秒的 R 问题。一种简单的解决方案是按照那里的建议添加 0.5 ms:
tt_ts <- strptime(tt[,1],"%Y-%m-%d %H:%M:%OS") + 0.0005
xts::xts(x=tt[,c(-1)], order.by=tt_ts)
# [,1]
# 2018-03-01 09:51:59.969 30755.5
# 2018-03-01 09:51:59.969 30755.0
# 2018-03-01 09:51:59.970 30755.5
# 2018-03-01 09:51:59.971 30756.0
# 2018-03-01 09:51:59.987 30756.5
# 2018-03-01 09:51:59.988 30756.5
我们可以从一个简单的例子中看出这一点:
st <- strptime("2018-03-01 09:51:59.971", "%Y-%m-%d %H:%M:%OS")
format(st, "%Y-%m-%d %H:%M:%OS3")
#> [1] "2018-03-01 09:51:59.971"
pt <- as.POSIXct(st)
format(pt, "%Y-%m-%d %H:%M:%OS3")
#> [1] "2018-03-01 09:51:59.970"
转换为POSIXct
ms后是错误的。提高输出精度,我们看到用于表示时间的浮点数刚好低于要求的值,但是 R 截断了数字而不是四舍五入:
format(pt, "%Y-%m-%d %H:%M:%OS6")
#> [1] "2018-03-01 09:51:59.970999"
移动所需精度的一半可以解决此问题。
format(pt + 0.0005, "%Y-%m-%d %H:%M:%OS3")
#> [1] "2018-03-01 09:51:59.971"
通常,如果x是 3 个十进制数字,则开放范围 ( x - 0.0005, x + 0.0005) 内的任何数字都将四舍五入为x。在截断时,这仍然适用于 [ x , x + 0.0005) 内的那些。但是如您所见, ( x - 0.0005, x ) 内的那些将由x - 0.001 表示。如果我们在截断之前将相关数字移动 0.0005,我们说的是范围(x,x + 0.001)。所有这些数字都将根据需要截断为x。
我将点x ± 0.0005 排除在外,因为对它们进行四舍五入的规则不同,并且表示时间点的实际浮点数将比这更接近所需值。
编辑:关于评论中关于采取差异的问题:如果将它添加到两个点,是否添加半毫秒都无关紧要。需要自行调整的时间点的示例:
st1 <- strptime("2018-03-01 09:51:59.971", "%Y-%m-%d %H:%M:%OS")
format(st1, "%Y-%m-%d %H:%M:%OS3")
#> [1] "2018-03-01 09:51:59.970"
pt1 <- as.POSIXct(st1)
format(pt1, "%Y-%m-%d %H:%M:%OS3")
#> [1] "2018-03-01 09:51:59.970"
format(pt1 + 0.0005, "%Y-%m-%d %H:%M:%OS3")
#> [1] "2018-03-01 09:51:59.971"
还有一个不需要调整的时间点:
st2 <- strptime("2018-03-01 09:51:59.969", "%Y-%m-%d %H:%M:%OS")
format(st2, "%Y-%m-%d %H:%M:%OS3")
#> [1] "2018-03-01 09:51:59.969"
pt2 <- as.POSIXct(st2)
format(pt2, "%Y-%m-%d %H:%M:%OS3")
#> [1] "2018-03-01 09:51:59.969"
format(pt2 + 0.0005, "%Y-%m-%d %H:%M:%OS3")
#> [1] "2018-03-01 09:51:59.969"
差异是相同的,与任何调整无关:
difftime(pt1, pt2, "secs")
#> Time difference of 0.001999855 secs
difftime(pt1 + 0.0005, pt2 + 0.0005, "secs")
#> Time difference of 0.001999855 secs
推荐阅读
- compiler-construction - 在为项目开发编译器时,我遇到了某种类型的生产问题
- php - 如何修复此错误“htmlspecialchars() 期望参数 1 为字符串,给定数组 (View:\frontend\passenger.blade.php)”
- laravel - 我可以在 laravel 和 mongodb 中更新我的密码吗?
- sql-server - CTE 基于输入创建递归动态行和列
- python - 如何在sqlalchemy中选择列存在
- java - 从 java 应用程序发送时没有在服务器中获取正确的 json
- python - Plotly-Dash:如何打印条形图的 hoverData?
- javascript - 如何在 jquery 中创建一个按钮,其 href 属性等于 ajax 响应
- php - 使用 php mailer 邮件从 localhost xamp 成功发送,但未能在生产 Windows 托管中发送
- angular - 如何在不重新加载整个站点的情况下更新子组件