如何在xts对象中提取多只股票每年最后一条月度观测值？

你说得对，to.yearly()确实主要是为OHLC（开盘/最高/最低/收盘）数据设计的，对于仅包含收盘价的单变量xts对象，它没法直接满足你的需求。下面给你推荐几个实用的方法，都能高效处理你那2000只股票的数据集：

方法1：用xts原生工具（最推荐，无需额外包） #

这是最贴合xts/zoo生态的方案，利用endpoints()定位每年的最后一个观测点，再用period.apply()提取对应数据：

# 生成每年最后一个观测值的位置索引
year_endpoints <- endpoints(myxtsfile, on = "years")

# 提取每个年份的最后一条数据
yearly_last_obs <- period.apply(myxtsfile, INDEX = year_endpoints, FUN = function(x) tail(x, 1))

endpoints()会帮你找出每个年份区间的结束位置，period.apply()则按这些位置拆分数据，对每个子区间执行tail(x,1)来取最后一行，完美匹配你的需求。

方法2：用zoo的聚合函数 #

因为xts是zoo的子类，所以可以直接用zoo的aggregate()函数按年份分组提取：

# 按年份分组，每组取最后一条观测值
yearly_last_obs <- aggregate(myxtsfile, as.integer(format(index(myxtsfile), "%Y")), tail, 1)

这里as.integer(format(index(myxtsfile), "%Y"))把xts的时间索引转换成年份整数，aggregate()以此为分组依据，对每组执行tail(x,1)获取最后一条数据，代码更简洁。

方法3：Tidyverse风格的处理（适合习惯dplyr的用户） #

如果你平时常用tidyverse工具链，可以结合timetk把xts转成 tibble，用更直观的分组逻辑处理：

# 先加载需要的包
library(dplyr)
library(timetk)
library(tidyr)
library(lubridate)

# 转换格式→分组提取→转回xts
yearly_last_obs <- myxtsfile %>%
  # 把xts转成带日期索引的宽格式tibble
  tk_tbl(preserve_index = TRUE, rename_index = "date") %>%
  # 转成长格式，方便按股票和年份分组
  pivot_longer(cols = -date, names_to = "symbol", values_to = "close") %>%
  # 提取年份
  mutate(year = year(date)) %>%
  # 按年份和股票分组，取每组最后一条
  group_by(year, symbol) %>%
  slice_tail(n = 1) %>%
  ungroup() %>%
  # 转回宽格式
  pivot_wider(names_from = "symbol", values_from = "close") %>%
  # 再转成xts对象（可选，根据你的需求决定）
  tk_xts(date_var = date)

这种方法逻辑更直观，如果你后续还要做其他数据清洗或分析，用tidy格式会更方便。

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