3
動的線形モデルのパラメーターの推定
私は(Rで)私は2つの持っているため、以下の非常に単純な動的な線形モデルを実装したい未知の時間変化のパラメータ(観測誤差の分散と状態誤差の分散ε 2 tで)。ε1tϵt1\epsilon^1_tε2tϵt2\epsilon^2_t Ytθt + 1==θt+ ϵ1tθt+ ϵ2tYt=θt+ϵt1θt+1=θt+ϵt2 \begin{matrix} Y_t & = & \theta_t + \epsilon^1_t\\ \theta_{t+1} & = & \theta_{t}+\epsilon^2_t \end{matrix} 先読みバイアスなしで、各時点でこれらのパラメーターを推定したいと思います。私が理解していることから、MCMC(先読みバイアスを回避するためにローリングウィンドウ上)またはパーティクルフィルター(またはシーケンシャルモンテカルロ-SMC)のいずれかを使用できます。 どちらの方法は、あなたが使用して、 これらの2つの方法の長所と短所は何ですか? おまけの質問:これらの方法では、パラメーターの変更速度をどのように選択しますか?大量のデータを使用してパラメーターを推定することと、少ないデータを使用してパラメーターの変更により迅速に反応することの間に取引があるため、ここに情報を入力する必要があると思いますか?
11
r
mcmc
dlm
particle-filter