每日时间序列(动物园)的冬季(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的子集,所以我想我可以使用以下方法提取所有的Decembers:

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

然后在Jan和Feb做类似的事情,但是我收到以下错误:

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)