3

我有一个日常动物园（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"))

然后为一月和二月做类似的事情，但我收到以下错误： 

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

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

4

1 回答 1

7

试试这个：

z.winter <- z[months(time(z), TRUE) %in% c("Dec", "Jan", "Feb")]
plot(z.winter)
于 2011-06-15T14:07:02.623 回答