r - Bursts パッケージを使用した R のバースト検出 (Kleinberg)

Question

> 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+ に高レベルのバーストがあることを示しており、元のプロットからはバーストが見られません。バーストプロットの読み方を理解してくれる人はいますか?

Answer 1

インデックスを x として、オフセットを y としてプロットしていますが、それらを切り替える必要があると思います。これを試して：

plot(offsets,1:length(offsets))