問題タブ [jags]

Rで次のjagsモデルを呼び出そうとしています:

これは階層モデルで、ficon は 0 ～ 60 の連続変数であり、国によって平均または分布が異なる場合があります。N = 総観測数 (2244) および J = 国の数 (34)。このモデルを実行すると、次のエラー メッセージが表示され続けます。

このコードは以前は機能していましたが、現在は機能していません。問題は 34 か国しかないことだと思います。そのため、i=35 で行き詰まっていますが、問題を解決する方法がわかりません。どんなアドバイスでも大歓迎です！

モデルを呼び出すために使用する R コード:

私はベイズと JAGS の両方に慣れていないので、私の無知を許してください。

JAGS コードを使用して作成された R スクリプトが (同僚から) 送られてきました。

このコードの作成者は、コーダ サンプルのセットを次のように定義しています。

私は以下のものを手に入れたいと思っていますが、限られた成功しか収めていません:

Gelman 診断: 単一のシミュレーションに適した "show(gelman.diag(codaSamples))" を使用しました。しかし、対象のすべてのシミュレーションについて、パラメーターごとに各ゲルマン診断をファイルに出力するにはどうすればよいでしょうか? さらに興味深いのは、Rhat 値が 1.1 を超えるシミュレーションの割合を記録することは可能ですか?

密度プロット: 「show(densplot(codaSamples))」を使用しました。ただし、これにより、各プロットが個別のプロットに生成されます (モデルには 96 個のパラメーターがあります)。ページごとに複数のプロットを配置する「autocorr.plot」と同等のものはありますか?

分位数: 「show (summary(codaSamples))」を使用しましたが、これにより各パラメーターの平均、SD、および特定の百分位数が得られましたが (これは私が望んでいたものです)、MCMC マトリックスも得られました。各パラメータの基本的な統計プロパティを指定する方法はありますか?

事後分布: 各パラメータの特定の値 (0 など) が下/上にあるセンタイルを計算する方法はありますか? それでは、すべてのシミュレーションを要約しますか?

あなたが提供できる助けを前もって感謝します。

Wabersich と Vandekerckhove (2013) の 26 ページにあるグラフィカル モデルを実装しようとしています。jags モデル スクリプトは次のように記述されています。

次に、次の初期リストを使用しました。

ただし、実行jags.modelするとエラーが発生します。

私のデータ構造は

このエラーは、初期化のためにジャグで非常に一般的であるようです。初期値の変更や削除を試みましたが、また失敗しました。エラーを削除するとdwiener、エラーが解決されます。

誰かが私を助けてくれれば幸いです。

私は本当に立ち往生していて、コーディングエラーが見つからないので、誰かがこれで私を助けてくれることを願っています!

JAGS（R2Jagsを使用）でゼロ膨張ポアソン/負の二項GLM（ランダム効果なし）をフィッティングしていますが、パラメーター推定、事前確率、初期値、およびチェーン収束ですべて問題ありません。すべての結果は、たとえば、モデルのピアソン残差の計算を含む、pscl-パッケージからの推定値と完全に一致しています...

私が作業できない唯一のことは、モデルから新しいサンプルをサンプリングして、モデルの適合性を評価するためのベイジアン p 値を取得することです。私が当てはめた「通常の」ポアソンおよび負の二項モデルはすべて、予想される複製サンプルを提供し、問題は発生しませんでした。

これまでのコードは次のとおりですが、重要な部分は「#New Samples」です。

大きな問題は、うまくいく/うまくいくはずだと思うすべての組み合わせと変更を実際に試したことですが、シミュレートされたサンプルを見ると、次のようになります。

さらに、Fit と FitNew の手段を使用する場合は、非常に奇妙になります。

私が間違っていることを誰かが手がかりを持っていますか? どんな助けでも大歓迎です！

よろしく、ウルフ

パッケージを使用してR で MCMC を実行していますが、後で別の R セッションで使用するためにrjags関数の出力を保存したいと考えています。jags.model

正規分布の平均の簡単な例を次に示します。

mu次のようなサンプルを生成できます。

ここで、後で新しい R セッションで使用するためにモデル オブジェクトを保存しjagsて、マルコフ連鎖を再度初期化して焼き付ける必要がないようにします。

ただし、新しい R セッションを開始してモデルをリロードすると、JAGS モデルを再コンパイルする必要があるというエラー メッセージが表示されます。

何故ですか？後で使用するためにオブジェクトを Rに保存するにはどうすればよいですか?jags

注：質問は以前に尋ねられましたが、OPは問題についてあまり具体的ではありませんでした。

JAGS を使用して、13 のサイト (9 は鳥類、4 つは潜在的な生息地) の林冠被覆率の平均と標準偏差を推定したいと思います。データが 0 と 1 でバインドされているという事実を説明するために、ベータ分布を使用しています。

他の分布 (ポアソン分布と対数正規分布) で完全に機能するモデル ステートメントのコードがあり、そのコードを適応させようとしましたが、惨めに失敗しました。

以下は、R コード、モデル ステートメント、およびデータです。Windows Vista で R 3.1.1 を使用しています。モデルステートメントを見て、私をまっすぐにしてくれれば、とても感謝しています.

ありがとう、

ジェフ

すべてのプロットを含む runjags オブジェクトが大きすぎるため、 を使用して、結果のオブジェクトをファイルに保存し、新しい R セッションで (as として) 復元してから、プロットを生成しようとしrun.jagsましたplot=FALSErunjagsout

(このトリックについては、こちらのディスカッションを参照してください: https://stackoverflow.com/a/21859618/684229 )

ただし、理由は不明ですが、これはモデルを再コンパイルして再度適応させます! を設定してもsample = 0, adapt = 0！

グラフをプロットするだけでもかなり時間がかかり、かなり面倒です。プロットを使用して runjags オブジェクトを計算し、それらを削除して runjags オブジェクトを小さく保存しようとすると、同じことが起こります。

この問題を修正する方法に関するヒントはありますか (独自のプロット関数を作成することは別として)?

警告:extend.jags同じ runjags オブジェクトで関数を2 回実行すると、すでに高速です。しかし、runjags オブジェクトを保存して新しいセッションで再度ロードすると、再びextend.jags遅くなります。または JAGS が何かをキャッシュしているようrunjagsです (ただし、元の runjags オブジェクト内ではありません)。

13 のグループ (植生の有無にかかわらず 20 からのポイント数として記録) と合計 223 のサンプルについて、中央値と、5、25、75、および 95 パーセンタイルのキャノピー カバーを推定したいと思います。以前、ベータ版配布を想定してこれを投稿しましたが、それは正しくありませんでした。これは提出期限が過ぎた原稿のためのもので、これが最後に欠落している部分です。誰かが私を完成させるのを手伝ってくれたら本当にありがたいです（コードが機能するまで）。私はそれに近いと思いますが、微調整が必​​要です-私は思います。

（私は編集して2つの反対票を修正しましたが、何が明確でないのかわかりません）。

大変感謝します！

以下は、私のモデル ステートメント、R コード、およびデータです。そのまま、私が得るエラーは

ただし、以下のスペースを削除したことに注意してください。エラーは可能性のステートメントを参照しています。

R コード:

データ：

MCMC と R の理解を深めるために改善しようとしている添付モデル (Rasmus Bååth の厚意による) があります。モデルにフィードされるのは、各行に (HomeTeam、AwayTeam、Season、HomeGoals、AwayGoals、MatchResult を含む) サッカーの試合です。 、ホームフォーム、アウェイフォーム)。このモデルは、ポアソン分布としてモデル化された HomeGoals と AwayGoals を使用して、2 つのチーム間の試合でサッカーの結果を予測するために使用されます。私がここ数日達成しようとしてきたことは、モデル内のチームの現在のフォームに関するこのデータを使用することですが、それを実現する方法が思いつきません。どのようにすればよいか、いくつかの提案を期待しています。

いくつかの説明:

• シーズンは、ゲームがプレイされたシーズンです
• MatchResult は、アウェイ チームが勝った場合は -1、引き分けの場合は 0、ホーム チームが勝った場合は 1 です。
• HomeForm と AwayForm は、それぞれホーム チームとアウェイ チームの現在のフォームを定量化する値です (-1 対 1、5 のウィンドウでのチームの最後のゲームのローリング平均試合結果、チームがアウェイ サイドでプレーした場合は負の試合結果)。 .
• テキスト内の列は、モデル (因子) に渡される前に整数に変換されます。

ありがとうございました！

モデル：

データ (の一部):

JAGS/WinBUGS でベクトル リテラルを記述する方法はありますか? JAGS でこれらを試してみましたが、すべてエラーになります。

それらを1つずつ割り当てたくありません。

ここまたはマニュアルで実際に役立つものは何も見つかりませんでした。