r - 从每日时间序列(动物园)中提取冬季(Dez、Jan、Feb)

Question

我有一个每日动物园 (xts)，其中包含以下格式的几十年数据:

head(almorol)
1973-10-02 1973-10-03 1973-10-04 1973-10-05 1973-10-06 1973-10-07
     183.9      208.2      153.7       84.8       52.5       35.5

我只想绘制冬季数据(12 月、1 月和 2 月的整个月份)。我找到了 xts 的子集，所以我想我可以使用以下方法提取所有十二月:

x<-apply.yearly(almorol, FUN=last(almorol, "1 month"))

然后对 1 月和 2 月执行类似的操作，但出现以下错误:

Error in get(as.character(FUN), mode = "function", envir = envir) :
object 'FUN' of mode 'function' was not found

我可以单独使用 apply.yearly 和 last(almorol, "1 month") 但是当我将它们组合时它不起作用。有谁知道对时间序列的这 3 个月进行子集化的方法？感谢您的帮助!

Answer 1

试试这个:

z.winter <- z[months(time(z), TRUE) %in% c("Dec", "Jan", "Feb")]
plot(z.winter)