問題タブ [winbugs]

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 は初めてで、コードのデバッグに苦労しています。

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

誰かが OpenBUGS で非線形モデルを適合させるのを手伝ってくれることを願っています。この問題はやや適切ではありませんが、多くのデータがあります。

マッピングされた樹木の森林区画で落ち葉を収集するトラップのネットワークがあります。

次のモデルを当てはめることで、樹木の直径とごみの生産量のアロメトリック関係を特定しようとしています。

直径「DBH」cmの木から「z」メートル離れた場所で収集されるごみの予想量

ただし、トラップ データには複数のツリーからの入力が含まれています。

データは次のようにフォーマットされます。

私はこれについて進歩を遂げ、投稿を数回更新しました。Gibb のサンプラーを実行できましたが、満足のいくミキシングが得られず、不合理な結果が得られません。観察されたごみの値は数十グラムですが、私は <10^-7 g のごみの投入量を見積もっています。

すべてのパラメータが正の値を取ると予想されるため、事前ガンマを使用することにしました。問題は事前にあると思います。私の意図は無知な事前確率を使用することですが、パラメーターが次の値を取るという漠然とした考えがあります: アルファ = 2 : 20 ベータ = ~2 ガンマ = .2 デルタ = .02

OpenBUGS と R パッケージを使用して二項混合モデルに取り組んでいますR2OpenBUGS。より単純なモデルを正常に構築できましたが、不完全な検出のために別のレベルを追加すると、一貫してエラーが発生しますvariable X is not defined in model or in data set。データの構造を変更したり、データを OpenBUGS に直接入力したりするなど、さまざまなことを試しました。他の誰かがこのエラーを経験したことがあり、OpenBUGS が変数 X を認識していない理由を知っていることを願ってこれを投稿しています。

エラーも発生しました。expected the collection operator c error pos 8これは以前に発生したエラーではありませんが、同様に困惑しています。

モデルとデータ シミュレーション機能の両方は、Kery の「Introduction to WinBUGS for Ecoologists (2010)」に由来します。ここでのデータ セットは、同様の私自身のデータの代わりであることに注意してください。

データセットとモデルを構築する関数を含めています。長々と申し訳ありません。

そして、ここにモデルがあります：