問題タブ [winbugs]

R で R2WinBUGS を実行するとエラーが発生します。再現性のコードを次に示します。

だから今、私は次のように見える2つのシミュレートされたデータセットを読みました

私の次のステップは、これらのデータセットで多変量ロジット正規モデルを実行し、そのようなシミュレートされたデータセットに対してこのプロセスを自動化することです。モデルステートメントは次のとおりです。

次のコードを使用して R2WinBUGS で実行しようとしています

どんな提案も役に立ちます。ありがとう！

したがって、次のコードは OpenBUGS で機能すると思っていましたが、代わりに「ノード Z の複数の定義」エラーが表示されます。

に置き換えZ <- round(X)ても、Z <- X引き続き同じエラーが発生します。この事実から、エラーは監視可能な変数の論理割り当ての使用に起因するものであり、特にエラーはround()操作によるものではないと推測できます。

バグがこれを許可しないのはなぜですか? また、この場合の良い回避策は何ですか? これは、私が実装したいより一般的なバージョンです。これは本質的に、壁 (切り捨て) を持つ離散ガウスをモデル化しています。

基本的に、Z「壁」(切り捨て) を持つ離散ガウス分布のようなもので分散したいと考えており、mu上のデータから推定したいと考えていZます。Zカテゴリ変数にしてパラメーターを推定しようとすると思いますが、これは理論的には苦痛のようです。意図したモデルを取得するために使用できるバグのトリックはありますか?

OpenBUGS モデルで 2 つのベクトルを差し引くのに非常に苦労しています。以下のコードの最後の行は、「予想される右括弧エラー」というエラーを出し続けます。

その行を に変更するとK[i,j+1] <- sum(ve[,i]) - sum(ve[,j])、モデルは正常に機能しますが、それは私がやりたいことではありません。要素ごとに減算したいと思います。

SO で OpenBUGS を検索しましたが、無関係なトピックがいくつかあります。

OpenBUGS - 変数が定義されていません

OpenBUGS: ベルヌーイ分布に値がありません

Stats Stack Exchange には、近いこの投稿がありますが、モデルにこれを実装する方法をまだ作成できませんでした:

https://stats.stackexchange.com/questions/20653/vector-multiplication-in-bugs-and-jags/20739#20739

for ループを書かなければならないことは理解していますが、これは確かに頭の痛い問題です。:)

その行を次のように変更してみました：

ここで、「p」は各「ve」の行数です。現在、「ノード temp[1] の複数の定義」というエラーが発生し続けています。

私は間違いなくいくつかの助けを使うことができました。それは非常に高く評価されます。

よろしくお願いします。

PS: この質問に「OpenBUGS」というタグを追加したかったのですが、残念ながらそれは新しいタグであり、十分な評判がないためできませんでした。代わりに「winbugs」を追加しました。

R2OpenBUGS の呼び出しに問題があります。おそらく簡単に修正できると思いますが、解決策を見つけることができません。

私の呼び出しは次のようになります。

返されるエラーは次のとおりです。

現在、OpenBUGS は開いていないため、ログ ファイルは書き込まれません。これが私が続けなければならないすべてであるため、何が起こっているのかわかりません。OpenBUGS 実行可能ファイルを見つけており、モデル ファイルを見つけています。これは、以前に取得して修正したエラーであるためです。

アイデアはありますか？

これのチュートリアル/例を見つけるのに苦労したので、質問したかったのです:i回測定される変数Xiがあります。追加の測定ごとに、Xの分布の予測がより厳密になることを示したかったのです。もちろん、モデルを 1:2 1:3 1:4 などで再実行し続けることもできますが、これは面倒です。私が気付いていなかった段階的なコーディングがいくつかあることを望んでいました。

その時点で利用可能なデータに基づいて、各時点で pred.x を推定するモデルを指定する方法を知っている人はいますか?

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 つのプロセスを結合するにはどうすればよいですか?

動物が取り除かれ、検出が時間と空間で変化する複数の観測期間にわたって n.sites から動物を検出する確率を推定しようとしています。5 つの観察期間で次のようなことをするとうまくいきます。

時間 2 の確率は時間 1 の確率に依存し、時間 3 の確率は時間 1 と 2 の確率に依存します。これを 5 期間だけ行う場合、書くのは大したことではありません。このアウト。しかし、10、15、20 以上の期間を取得すると、書き出すのはかなり面倒です。各ステップを入力せずにこのループを記述する方法があるはずだと思いますが、その方法が思いつきません。おそらく、追加の索引付けまたは他の制御ステートメントまたは累乗関数です。p[i] が各 j 番目の観測で同じである場合 (つまり、p[i,1] = p[i,2] = p[i,3] など)、次のようになります。

どんな提案でも大歓迎です。

これはバグ言語コードです。私は R で作業しており、コードを rjags パッケージ経由で JAGS に送信しました。BUGS、R、または疑似コードが私の目的に適しています。

問題をシミュレートする R コードを次に示します。

助けてくれてありがとう。ダン

アタッチされたモデルで 3 つのチェーンの初期化をロードすると、「このチェーンには初期化されていない変数が含まれています」というメッセージが表示されます。続行するために「gen init」を押しました。それは正しいことですか？固定モデルでは発生しませんが、ランダム モデルでは発生する傾向があります。お知らせ下さい。

R パッケージ「BRugs」を使用して Gibbs サンプラーを実装しようとしていますが、samplesStats() などの要約事後統計を生成する関数は、平均値と中央値のみを返します。事後モードを抽出することは可能ですか?

WinBUGS/OpenBUGS は初めてで、コードのデバッグに苦労しています。

バグモデルの潜在的なエラーメッセージとその意味のリストを平易な英語で知っている人はいますか?