Saat ini saya menggunakan proses berikut untuk bootstrap seri waktu multivarian di R:
- Tentukan ukuran blok - jalankan fungsi
b.star
dalamnp
paket yang menghasilkan ukuran blok untuk setiap seri - Pilih ukuran blok maksimum
- Jalankan
tsboot
pada seri apa saja menggunakan ukuran blok yang dipilih - Gunakan indeks dari output bootstrap untuk merekonstruksi rangkaian waktu multivarian
Seseorang menyarankan untuk menggunakan paket meboot sebagai alternatif dari blok bootstrap tetapi karena saya tidak menggunakan seluruh kumpulan data untuk memilih ukuran blok, saya tidak yakin bagaimana menjaga korelasi antara seri jika saya menggunakan indeks yang dibuat dengan menjalankan meboot
pada satu seri. Jika ada yang punya pengalaman dengan meboot dalam pengaturan multivarian, saya akan sangat menghargai saran tentang prosesnya.