我在11個月內收集了來自9次調查的許多物種的調查數據。調查跨度從12月開始到10月結束。 5月或6月沒有任何物種的數據。以下是一個簡化的例子,雖然我的數據有很多豐富的數據列。如何指定一個破損的時間序列?
as.yearmon
使用從封裝動物園
require(zoo)
month <- as.yearmon(month, "%b%y")
一些假丰度數據
meanAbund <- rlnorm(9, 1, 1.2)
meanAbund2 <- meanAbund + 0.5 # to have two series at least
日期
month <- c("Dec10", "Jan11", "Feb11", "Mar11", "Apr11", "Jul11", "Aug11", "Sep11", "Oct11")
編纂所有接合在一起
df <- cbind.data.frame(meanAbund, meanAbund2, month)
儘管我很欣賞這個時間序列不夠長,無法進行時間序列分析,但我想指定這些數據是一個(斷開的)時間序列對象,以便圖表將打破這個序列而不是插入或忽略間隙。這應該怎麼做?
我懷疑這應該很容易,但我還沒有破解解決方案。我想要製作的東西可以在Shah,Zeileis和Grothendieck的quick reference guide的第4頁看到,但我看不到代碼告訴他們在哪個步驟中他們的例子。
我可以讀取它作爲一個時間序列對象
df.ts <- zoo(df[ , 1:2], df$month])
並繪製系列(最高到63,從我的實際數據集,雖然輸出是很難考慮你可能會懷疑),
plot.zoo(df.ts)
但該系列忽略了差距。
我應該如何指定圖形或讀入這些數據,以便繪圖函數「知道」打破系列?
請提供您的繪圖代碼 – HubertL
@HubertL我現在編輯了我的帖子,提供了我已經制作的情節的步驟(儘管按照Pierre Lafortune的例子整理)。我可以使用'plot.ts'或多個不使用'plot.zoo'斷點的情節來獲得單個情節或多個空白情節的休息時間。現在看來是一個不同的問題或問題。 – dhd