2011-06-15 17 views
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的子集,所以我想我可以使用提取所有的12月份:

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.yearlylast(almorol, "1 month"),但是當我將它們組合時,它不起作用。有沒有人知道對這3個月的時間序列進行子集化的方法?感謝您的幫助!

+0

請使用'dput(頭(almorol))'提供數據的樣本。也可以輸入(頭(almorol,15)) – Marek 2011-06-15 13:41:41

+0

dput(頭(almorol,15)) 結構(c(183.9,208.2,153.7,84.8,52.5,35.5,49.2,77.7, 73,104.9,124.8 ,103,75.9,108.7,138.2),指數=結構(c(1370, 1371,1372,1373,1374,1375,1376,1377,1378,1379,1380,1381, 1382,1383,1384),等級=「Date」),class =「zoo」) – sbg 2011-06-15 13:47:21

回答

6

試試這個:

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

謝謝!我不知道'months()'技巧! – sbg 2011-06-15 14:21:36

相關問題