如何在R中用变量获取XTS数据框最后X条观测值？含时间设置与报错处理

问题说明 #

我在R里尝试计算XTS数据集最后X分钟的平均值，目前不知道怎么把起始日期设为今日午夜，暂时用结束日期减86400秒代替，写了如下代码：

xts <- xts::xts(tempdat, order.by = data()$date)
enddate <- end(xts)
todaydate <- enddate - 86400 # 不知道怎么获取今日午夜的日期
todaymins <- difftime(enddate, todaydate, units="hours")
lastavg <- xts |> utils::tail(todayhrs) |> mean(na.rm =TRUE)

运行后报错：

Warning: Error in ==: comparison of these types is not implemented
  101: Summary.difftime
  100: tail.zoo
   97: ::

想解决两个问题：

• 怎么用变量获取XTS数据框的最后X条观测值？
• 如何把起始日期设置为今日午夜？

解决方法 #

1. 获取今日午夜的日期

直接用Sys.Date()配合as.POSIXct()就能生成今日午夜的时间戳——Sys.Date()返回当日日期，转成POSIXct类型时默认就是午夜时刻：

today_midnight <- as.POSIXct(Sys.Date())

如果你的XTS数据时区和系统时区不一致，记得指定时区参数，比如：

today_midnight <- as.POSIXct(Sys.Date(), tz = "UTC")

也可以匹配XTS数据本身的时区：

today_midnight <- as.POSIXct(Sys.Date(), tz = attr(index(xts), "tzone"))

2. 获取最后X条/最后X分钟数据并计算均值

你之前报错的原因是tail()接收的参数类型不对——你传了difftime类型的对象，但tail.zoo（xts继承自zoo）需要整数型的观测条数，或者直接用时间范围筛选更准确（毕竟最后X分钟的数据条数不一定固定）。

方式一：按时间范围筛选（推荐）

如果要获取今日午夜到当前的所有数据，直接用XTS的时间索引筛选：

# 筛选今日午夜之后的数据
today_data <- xts[paste0(today_midnight, "/")]
# 计算均值
today_avg <- mean(today_data, na.rm = TRUE)

如果要获取最后X分钟的数据，先算出X分钟前的时间点再筛选：

x_minutes <- 60 # 示例：取最后60分钟
start_time <- end(xts) - (x_minutes * 60) # 分钟转成秒数计算
last_x_data <- xts[paste0(start_time, "/")]
last_x_avg <- mean(last_x_data, na.rm = TRUE)

方式二：按固定条数筛选

如果确实需要固定条数的最后X条，把参数转成整数即可：

x_rows <- 100 # 示例：取最后100条
last_x_data <- tail(xts, n = x_rows)
last_x_avg <- mean(last_x_data, na.rm = TRUE)

修正后的完整示例代码

library(xts)

# 构建XTS对象（假设tempdat是数据列，data()$date是时间列）
xts_obj <- xts(tempdat, order.by = data()$date)
end_date <- end(xts_obj)

# 获取匹配数据时区的今日午夜时间
today_midnight <- as.POSIXct(Sys.Date(), tz = attr(index(xts_obj), "tzone"))

# 计算今日所有数据的均值
today_data <- xts_obj[paste0(today_midnight, "/")]
today_avg <- mean(today_data, na.rm = TRUE)

# 计算最后30分钟数据的均值
x_minutes <- 30
start_time <- end_date - (x_minutes * 60)
last_30min_data <- xts_obj[paste0(start_time, "/")]
last_30min_avg <- mean(last_30min_data, na.rm = TRUE)

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