问题描述
|
我有一个每日动物园(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
我只想绘制冬季数据(十二月,一月和二月的整个月)。我找到了xts的子集,所以我想我可以使用以下命令提取所有的十二月:
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\")
,但是当我将它们结合使用时,它们将不起作用。有谁知道细分这三个月的时间序列的方法吗?感谢您的帮助!
解决方法
尝试这个:
z.winter <- z[months(time(z),TRUE) %in% c(\"Dec\",\"Jan\",\"Feb\")]
plot(z.winter)