ifelse() 从时间戳向量中剥离 POSIXct 属性？

Question

这很奇怪：R 的 ifelse() 似乎做了一些（不需要的）转换： 假设我有一个时间戳向量（可能是 NA），NA 值的处理方式应与现有日期不同，例如，只是忽略：

formatString = "%Y-%m-%d %H:%M:%OS"
timestamp = c(as.POSIXct(strptime("2000-01-01 12:00:00.000000", formatString)) + (1:3)*30, NA)

现在

timestamp
#[1] "2000-01-01 12:00:30 CET" "2000-01-01 12:01:00 CET" "2000-01-01 12:01:30 CET"
#[6] NA    

根据需要，但翻译 30 秒会导致

ifelse(is.na(timestamp), NA, timestamp+30)
#[1] 946724460 946724490 946724520        NA

请注意，timestamp+30 仍按预期工作，但假设我想用固定日期替换 NA 日期并将所有其他日期翻译 30 秒：

fixedDate = as.POSIXct(strptime("2000-01-01 12:00:00.000000", formatString))
ifelse(is.na(timestamp), fixedDate, timestamp+30)
#[1] 946724460 946724490 946724520 946724400

问题：这个解决方案有什么问题，为什么不能按预期工作？

编辑：所需的输出是一个转换为 30 秒的时间戳（不是整数）向量，并且 NA 被替换为任何内容......

Answer 1

如果您查看ifelse 的编写方式，它有一段代码如下所示：

ans <- test
ok <- !(nas <- is.na(test))
if (any(test[ok]))
  ans[test & ok] <- rep(yes, length.out = length(ans))[test & ok]

请注意，答案以逻辑向量开始，与测试相同。然后将具有test == TRUE 的元素分配给yes 的值。

这里的问题是，将逻辑向量的一个或多个元素分配为 POSIX.ct 类的日期会发生什么。如果你这样做，你可以看到会发生什么：

x <- c(TRUE, FALSE)
class(x)
# logical
x[1] <- Sys.time()
class(x)
# numeric

你可以通过写作来解决这个问题：

timestamp <- timestamp + 30
timestamp[is.na(timestamp)] <- fixedDate

你也可以这样做：

fixedDate = as.POSIXct(strptime("2000-01-01 12:00:00.000000", formatString))
unlist(ifelse(is.na(timestamp), as.list(fixedDate), as.list(timestamp+30)))

这利用了替换运算符[<- 处理右侧列表的方式。

你也可以像这样重新添加类属性：

x <- ifelse(is.na(timestamp), fixedDate, timestamp+30)
class(x) <- c("POSIXct", "POSIXt")

或者，如果您不顾一切地想像这样在一行中做到这一点：

`class<-`(ifelse(is.na(timestamp), fixedDate, timestamp+30), c("POSIXct", "POSIXt"))

或者通过复制fixedDate的属性：

x <- ifelse(is.na(timestamp), fixedDate, timestamp+30)
attributes(x) <- attributes(fixedDate)

最后一个版本还具有复制tzone 属性的优点。

从 dplyr 0.5.0 开始，您还可以使用 dplyr::if_else 在输出中保留类，并为真假参数强制使用相同的类。

Answer 2

正如 Henrik 所说，ifelse() 会剥离属性，这与简单的 for 循环不同。

在不悲伤的情况下填充 NA 的解决方法是更简单、更清晰的函数 zoo::na.fill

那么你会这样做：na.fill(timestamp, fixedDate)

另请参阅 na.locf, na.approx, na.spline ...，其他来自 zoo 的出色便利功能。