問題タブ [openbugs]

私は最近、OpenBUGS と R を使用してベイジアン モデリングを試み始めましR2OpenBUGSた。

私は David Eagles によるワインと OpenBUGS をインストールするための素晴らしいリンクをたどり、そこに投稿された学校の例はうまくいきました。

次に、持っていたコードを実行しようとしましたがthe model is syntactically correct、データが正しく読み込まれていないことを伝えた後、エラーが発生し続けました。

数日間のトラブルシューティングの後、データファイルが特定の長さを超えると、トップビット (listデータをロードするために OpenBUGS で従来ハイライトする場所) が失われ、OpenBUGS でデータをロードできないようです。

リストセクションが切り取られた、読み込まれたデータを含む OpenBUGS プログラムのスクリーンショット

list{この図は、data.txt ファイルの上部に、モデルを実行しようとするすべての試みを基本的に停止している がなくなったことを示しています。

これを確認するために、ポイント数を 1000 に設定してから 100 に設定する単純な線形回帰の例を実行しました。

の場合、モデルはうまく動作linedata <- list(Y = c(1:100), x = c(1:100), N = 100, xbar = 3)しますが、 の場合、データの読み込み段階で崩壊しlinedata <- list(Y = c(1:1000), x = c(1:1000), N = 1000, xbar = 3)ます。

wine と OpenBUGS のデータ ファイル サイズの制限について Google 検索を使い果たしましたが、役立つものは何も見つかりませんでした。

何を試すか、どこから始めるか、これを以前に経験したことについて何か提案はありますか?

OpenBUGS で次のモデルを初期化する際に問題が発生しています

これは私のデータです

これが私のイニシャルのリストです

モデルのチェック、データの読み込み、コンパイルの手順は問題ありません。イニシャルをロードするとき、OpenBUGS は、初期値はロードされているが、チェーンには初期化されていない変数が含まれていると言います。次に、theta0 も初期化しようとしましたが、問題は解決しません。誰かがこれに関して私を助けてくれますか? ありがとうハリド

私は R2OpenBUGS を初めて使用し、非常に謎めいたエラーにかなりイライラしています。

私は非常に単純なモデルを実行しようとしています。以前に同様のモデルを実行することに成功しました。私の問題は、私が 2 次元配列 (行列) を持っているという事実によるものですか?

モデルを単純化しようとしましたが成功しませんでした。

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

そして、ここに私が使用するコードがあります

データファイル:

私は JAGS/OpenBUGS をまったく初めて使用するので、モデルの指定に関して正しい方向に進んでいただければ幸いです。私は、103 か国が 15 年間にわたってまとめたアンバランスな縦断データを使用しています。この場合は 12 年が選択されています。DV はジニ係数であり、対数正規ではなく、おそらくベータでモデル化する必要がありますが、現在は JAGS でモデルをコンパイルする方法を理解することに重点が置かれています。当分の間、固定効果モデルを使用しています。

私が実行しているデータとコード:

そして、私はこの実行時エラーを受け取り続けます:

私が別の方法で何をすべきかについての提案は、非常に高く評価されます! 不均衡なデータに適用できる 3 つの「トリック」があることは知っていますが、これらすべてがどのように機能するか、JAGS がデータ入力を読み取る方法についてはまだ少し混乱しています。

乾杯

J

R2OpenBUGS ライブラリを使用して、R から階層モデルを評価しようとしています。私はすでに同様の質問への回答を見てきましたが、何もうまくいかなかったようです。誰かが助けてくれたらとてもうれしいです。

モデルは以下。

ログで私が得る：

問題は、a[i] と d[i] に事前手段のベクトルを与えようとしたときに始まります。これらの手段を「0」に設定すると、コードが機能します。

x1b1+x2b2+……+xpbpを計算したい。ただし、pは不明で、事前に書き留めることはできません。だから私は次のことをしたい。しかし、バグでは、複数の定義エラーが表示されます。

バグモデルでそれを行う方法はありますか?

OpenBUGS でマルコフ モデルを実行して、遷移確率を推定しようとしています。しかし、「配列インデックスは整数ではありません」というエラーが表示されます。これが私のモデル コードです。

そして、ここにRコードを含む私のデータがあります:

誰でも親切にエラーを解決するのを手伝ってもらえますか?

ベイジアン推論用に R の BRugs パッケージを使用してモデルを実行しています。モデルは構文的に正しく、データは読み込まれますが、modelCompile 関数は永久に実行されます。15 時間後に停止しましたが、何が問題なのかわかりません。

データファイルはかなり大きいですが、リクエストに応じて提供できます...

私は OpenBUGS を初めて使用し、logit() 関数を使用してモデルを適合させる際に問題が発生しました。

周りを読んで、これに対する1つの可能な解決策は、WinBUGS独自のロジット関数を使用せずにロジット関数を明示的に指定することであることがわかりました:

より複雑なモデルでは、WinBUGS 独自の logit 関数を使用するときに、たとえば収束を達成する際に問題が発生することがよくあります (実際には、かなり単純なモデルでも問題が発生する可能性があります)。したがって、その変換を logit.p[i] <- log(p[i] / (1 – p[i])), p[i] <- exp(logit.p[i] で明示的に指定する方がよい場合がよくあります。 ]) / (1 + exp(logit.p[i])) または p[i] <- 1 / (1 + exp(- logit.p[i]))。

(詳細はこちらhttp://www.mbr-pwrc.usgs.gov/software/kerybook/AppendixA_list_of_WinBUGS_tricks.pdfのポイント 14.)。

問題は、それを行う方法がわからないことです。WinBUGS 統合ロジット関数を使用した元の尤度関数が次のとおりだったとします。

どのように私はそれを明示的に書きますか？

どうもありがとうございました。

ヴィンチェンツォ