lubridate包parse_date_time()解析时间偶现1微秒偏差问题咨询

parse_date_time()解析带小数秒的时间时出现的微秒位偏差，本质是双精度浮点数的固有精度限制：R的POSIXct时间类型用双精度浮点数存储从纪元起算的秒数，部分十进制小数秒没法被二进制浮点数精确表示，解析过程里进制转换的舍入就会导致最终值和输入差1微秒，parse_date_time()本身没有内置参数可以直接调高计算精度来规避这个问题。

• 手动拆分拼接时间绕开内置解析的浮点误差：先把时间字符串拆成整秒部分和小数秒部分，整秒部分用parse_date_time()正常解析，小数秒部分单独按十进制数值算好后加到解析完的整秒时间上，能把偏差控制到可忽略的范围。参考实现：

library(data.table)
library(lubridate)

parse_high_precise_time <- function(time_str, tz = "") {
  # 按小数点拆分整段时间和小数秒
  parts <- tstrsplit(time_str, "\\.", fixed = FALSE)
  # 解析整秒部分
  base_time <- parse_date_time(parts[[1]], orders = "Y-m-d H:M:S", tz = tz)
  # 统一补全6位微秒，避免小数位数不一致导致计算错误
  frac_pad <- paste0(parts[[2]], "000000")
  micro_sec <- as.numeric(substr(frac_pad, 1, 6)) / 1000000
  # 加回小数秒得到最终时间
  base_time + micro_sec
}

• 如果需要长期处理高精度时间数据，直接换存储类型：别用POSIXct类型存微秒/纳秒级时间，改用nanotime包提供的纳秒级时间类型，或者用64位整数直接存从纪元起算的微秒/纳秒时间戳，从根上避免双精度浮点数的精度损失。
• 注意：options(digits.secs)只控制时间打印时展示的小数位数，不会改底层存储的时间数值，调这个参数解决不了实际的舍入偏差。

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