2013-10-30 62 views
0
> library(bursts) 
> par(mfrow=c(2,1)) 
> offsets <- c(seq(0, 400, 100), seq(410, 450, 5), seq(451, 470, 2), 
+    seq(480, 600, 5), 700, seq(710, 800, 5), 900, 1000) 
> bursts <- kleinberg(offsets) 
> plot(offsets) 
> plot(bursts) 

遵循來自突發包的演示代碼。我在下面繪製了一張情節,我無法將爆發結果與原始數據關聯起來。爆裂圖表明在時間400+處有高電平爆發,而我從原始圖中看不到任何爆發。任何人都可以幫助我理解如何閱讀爆發情節?帶爆發包的R中的突發檢測包(Kleinberg)

enter image description here

+0

是可以適應您的時間戳的代碼? –

回答

1

您正在策劃指數爲x和在y軸偏移,但我想你應該切換他們。試試這個:

plot(offsets,1:length(offsets)) 

enter image description here

+0

我的不好..非常感謝解釋。 –