問題タブ [r2jags]

JAGSとBUGSでモデルを実行する方法の詳細について簡単に質問します。

n.burnin=5000、、n.iter=5000およびを使用してモデルを実行するとしthin=2ます。これは、プログラムが次のことを行うことを意味しますか？

1. 5,000回の反復を実行し、結果を破棄します。その後
2. 1秒おきの結果のみを保持しながら、さらに10,000回の反復を実行しますか？

これらのシミュレーションをCODAオブジェクトとして保存すると、10,000個すべてが保存されますか、それとも間引きされた5,000個だけが保存されますか？ACFプロットを作成するためにどの反復セットが使用されているかを理解しようとしていますか？

R2jags と CODA を使用して、MCMC チェーンの診断統計を作成していますが、問題があります。次のようにMCMCを実行したい：

エラーは次のとおりです。

RStudio 0.97.551 と R バージョン 3.0.0 (2013-04-03) を使用しています。

どんな助けにも感謝します！

これは私のRスクリプトです：

JAGS で GARCH(1,1) モデルを構築しようとしています。簡単にするために、平均方程式が AR(1) プロセスに従うと仮定します。AR(1) および GARCH(1,1) プロセスに参加できるようにする 1 つの JAGS モデルを構築しようとしています。

今のところ、2 つの別個の JAGS モデルを構築することによってのみ結果を得ることができます (プレゼンテーションを明確にするために簡略化しています)。最初の JAGS モデルは、AR(1) プロセスのパラメーターを推定します。

パラメータの推定値を取得して、AR(1) プロセスのデータを生成し、残差と分散を取得します (何らかのウィンドウを想定)。

次のブロックは、JAGS の GARCH(1,1) 散文です。

依存している 2 つのプロセスを結合するにはどうすればよいですか?

jags.model() オブジェクトを RData または txt ファイルに保存する方法はありますか?

より優れたコンピューターで MCMC を実行するには、モデルを 1 台に保存し、それを新しいワークスペースで使用する必要があります。しかし、Rで「save()」と「load()」関数を使用するのは少し難しいです。アドバイスをありがとう。

追加した：

私は試した：

jags <- jags.model('regression.bug', data = my.data, n.chains = 4, n.adapt = 1000)

次に、「ジャグ」を保存したいと思います

保存 (jags , file="jags.RData")

保存されているようです。しかし、私がしようとすると：

ld.jags <- load( "jags.RData" )

ld.jags

[1]「ギザギザ」

そして、「ld.jags」を使用して分析を実行する方法がわかりません。

JAGS (または WinBUGS) を使用してベイジアン線形混合モデルを当てはめたとします。出力オブジェクトにはモデルの残差が含まれますか? どうすれば残差を見つけることができますか?

ありがとう！

ポーリング データを分析的に処理するだけでなく、ベイジアン推論を使用するように言われました。ただし、問題があります。政党の以前の分布に関する推測を含む小さなデータセットと、世論調査のデータがあります。Gibbs シミュレーションから限界値を取得するにはどうすればよいですか?

R に相当するギザギザは何is.na()ですか? r2jags を使用して R データ オブジェクトを jags に渡しています。私の jags モデルでは、ベクトルの要素に R の NA 値と同等の値が含まれているかどうかをテストしたいと考えています。Rが要素がNAであると言ったときにジャグが「見る」ものはわかりません。