次の形式の数十年のデータを含む毎日の動物園 (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 のサブセットを見つけたので、次を使用してすべての 12 月を抽出できると考えました。
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 か月をサブセット化する方法を知っている人はいますか? 助けてくれてありがとう！