如何对1985-1994年的xts日度序列按跨年度区间（11月至次年3月）求和？

我们需要对1985-01-01至1994-12-31的日度xts时间序列，计算每个当年11月1日到次年3月31日区间内的数值总和，比如1985-11-01至1986-03-31、1986-11-01至1987-03-31这类跨年度区间的总和。

实现代码 #

# 生成原始时间序列（用户提供的示例代码）
dates <- seq(from = as.Date("1985-01-01"), to = as.Date("1994-12-31"), by = 1)
data <- rnorm(length(dates))
my.ts <- xts(x=data, order.by=dates)

# 确定所有需要计算的区间起始年份（1985到1993，对应最后一个区间1993-11至1994-3）
start_years <- 1985:1993

# 生成每个区间的起止日期，并计算总和
sum_results <- sapply(start_years, function(y) {
  start_date <- as.Date(paste0(y, "-11-01"))
  end_date <- as.Date(paste0(y + 1, "-03-31"))
  # 提取区间数据并求和
  sum(my.ts[paste(start_date, end_date, sep = "/")])
})

# 将结果转为xts对象，用区间结束日期作为索引
my.ts2 <- xts(sum_results, order.by = as.Date(paste0(start_years + 1, "-03-31")))

代码说明 #

• 首先通过start_years确定所有需要计算的区间起始年份，确保覆盖到1993-11至1994-3的最后一个区间。
• 利用sapply遍历每个起始年份，生成对应区间的起止日期，再通过xts的切片语法paste(start_date, end_date, sep="/")提取区间内的数据并求和。
• 最后将求和结果转为xts对象，用每个区间的结束日期作为时间索引，方便后续的时间序列操作。

由于原始序列中没有2月29日的数据，无需额外处理闰年的日期匹配问题。

内容的提问来源于stack exchange，提问作者San-Diego-WX