lubridate::floor_date毫秒级向下取整结果不准确，求替代方案

我在处理数据库表时，部分表存储的时间戳精度为毫秒。出于技术要求，需要将这些时间戳向下取整到完整毫秒（即floor操作）。

尝试用lubridate::floor_date()并设置unit = ".001s"，但结果不准确：偶尔会生成带有毫秒以下精度的时间戳，甚至出现跳过毫秒值的情况。

示例代码 #

options(digits.secs=7)
for (i in seq_len(100)) {
  print(lubridate::floor_date(Sys.time(), unit = ".001s"))
}

观察结果（去重后） #

[1] "2026-06-15 11:17:51.882999 CEST"
[1] "2026-06-15 11:17:51.884 CEST"
[1] "2026-06-15 11:17:51.884999 CEST"
[1] "2026-06-15 11:17:51.885999 CEST"
[1] "2026-06-15 11:17:51.887 CEST"
[1] "2026-06-15 11:17:51.888 CEST"
[1] "2026-06-15 11:17:51.888999 CEST"
[1] "2026-06-15 11:17:51.89 CEST"

预期结果 #

• 所有输出必须是精确的毫秒级时间戳，例如51.884999应变为51.885
• 不能跳过任何毫秒值（部分结果因略低于精确毫秒值导致跳过，即使设置options(digits.secs=3)也存在该问题）

推测问题源于数值精度以及R存储秒以下时间戳的方式，但正是为了避免下游数据库的精度问题，才需要严格向下取整到完整毫秒。虽然?lubridate::floor_date说明支持秒的分数单位，但实际处理毫秒场景时表现不稳定。

需求：寻找无需转换字符串的替代方法，实现POSIXct类型时间的精确向下取整到毫秒。

解决方案：基于数值运算的精确取整 #

由于POSIXct本质是存储从1970-01-01以来的秒数（含小数），可以直接对数值进行floor运算，再转回POSIXct类型，完全避免lubridate的精度问题：

floor_ms <- function(x) {
  as.POSIXct(floor(as.numeric(x) * 1000) / 1000, origin = "1970-01-01", tz = attr(x, "tzone"))
}

# 测试
options(digits.secs=7)
for (i in seq_len(100)) {
  print(floor_ms(Sys.time()))
}

原理说明 #

1. as.numeric(x)将POSIXct时间转为秒级数值（含小数，比如1718438271.884999）
2. 乘以1000后得到毫秒级数值（1718438271884.999），用floor()取整得到整数毫秒数1718438271884
3. 除以1000转回秒级数值（1718438271.884）
4. as.POSIXct()将数值转回带时区的POSIXct类型，保留原时间对象的时区信息

这种方法直接操作底层数值，完全规避了lubridate在处理分数秒单位时的精度误差，输出结果严格符合毫秒级要求，不会出现跳过或多余小数位的情况。

内容的提问来源于stack exchange，提问作者Falk Mielke