多変量時系列のブロックブートストラップの代替


11

私は現在、Rで多変量時系列をブートストラップするために次のプロセスを使用しています。

  1. ブロックサイズの決定- 各シリーズのブロックサイズを生成するパッケージb.star内の関数を実行しますnp
  2. 最大ブロックサイズを選択
  3. tsboot選択したブロックサイズを使用して任意のシリーズで実行
  4. ブートストラップ出力のインデックスを使用して、多変量時系列を再構築します

誰かがブロックブートストラップの代わりにmebootパッケージを使用することを提案しましたが、データセット全体を使用してブロックサイズを選択していないため、実行mebootして作成されたインデックスを使用する場合にシリーズ間の相関関係を保持する方法がわかりません1つのシリーズ。多変量設定でmebootの経験がある人がいれば、そのプロセスについてのアドバイスをいただければ幸いです。

回答:


8

まず、最大のエントロピーブートストラップ(meboot)を強くお勧めします。私はmebootを優先してブロックブートストラップを放棄しましたが、その結果には非常に満足しています。このアルゴリズムはブロッキングをまったく使用せず、定常性を必要としませんが、データの相関構造を組み込んでいます。かっこいいね。

次に、mebootを使用して多変量ブートストラップを実行したことがないことを認めますが、時系列データをパネルデータとして再キャストし、meboot.pdata.frame関数を使用して基本的に多変量ブートストラップを実行できます。


関数はmeboot.pdata.frameだと思います。
ProbablePattern 2011

うわぁ!PPさん、ありがとうございます。関数のスペルを間違えました。正しい名前はmeboot.pdata.frameです。ごめんなさい。
pteetor

私はこのスレッドに出会い、MEブートストラップに興味を持ちました。少し遊んでみたところ、時系列が非常に密接に複製されていることがわかりました。ブートストラップDGPの変動はそれほど多くありません。さらに掘り下げてみると、MEブートストラップの批評を含む、Davidsonによる次の論文を見つけました。www.monticini.eu/ wp / rdavidson.pdf。それに対する反論があったかどうかはわかりませんが、私自身のシミュレーション結果を考えると、私はそれを使用するのをためらいます。
generic_user

@generic_userそのリファレンスをありがとう!とても興味を持って読みました。私も、特定の状況下でのブートストラップ複製の変動性の欠如に言及しました。結果として、MEブートストラップが有効になることはないと私はまだ確信していません。しかし、紙は地味です。さらに調査が必要です。
pteetor 2016
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.