問題タブ [r2winbugs]

R 内から R2WinBUGS パッケージを使用して、WinBUGS (100 万回の反復) で長い計算を実行しました。

しかし、Rはクラッシュしました。winbugs を再度実行せずに bugs.object を R に再ロードするにはどうすればよいですか? 私はこれを試しました（私は3つのチェーンを持っています）：

しかし、outデータ構造は bugs オブジェクトとはまったく異なります (そのままでは使用できません)。私はそれを変換しようとしましたas.bugs.array:

しかし、うまくいきません。助けてください。ありがとう。

この例を使用して R2WinBUGS を使用しようとしています:

コード

(次の部分だけ考えてください: ### 5.4. WinBUGS を使用した分析)

次のエラー メッセージが表示されます。

これが正しい機能にとって重要かどうかはわかりません (他のすべては問題ないように見えます)。これを取り除く方法はありますか？

ありがとう。

キリスト教徒

PS:

これはRコードです：

誰かがパッケージR2WinBUGSを使用してロジスティック回帰を実行するRコードを持っているかどうか疑問に思っています。理想的には、シミュレートされたデータを使用して「真実」と2つの連続共変量を生成します。

ありがとう。

キリスト教徒

PS：

人工データ（1次元の場合）を生成し、r2winbugsを介してwinbugsを実行するための潜在的なコード（まだ機能していません）。

R（R2OpenBUGS）を介してOpenBUGSを使用して、観測「時間」を確率変数としてモデル化しようとしています。すべての観測時間が利用可能である（NAがない）場合、すべてが機能しますが、いずれかの時間をNAに設定すると、何も起こりません。WinBUGSで同じコードをテストしたところ、トラップエラー「NIL逆参照（読み取り）」が発生しました。だから私の質問は、私のコードに本当に何か問題があるのか​​、それとも私のモデルがバグに対してあまりにも奇妙なのかということです。

私のモデルは次のようなものです。

そして、Rコードは次のようになります。

R から WinBUGS を呼び出そうとしていますが、うまくいきません。

Rでこのコードを入力しました：

次に、WinBUGS を呼び出します。

そして、WinBUGS で次のエラーが発生します。

誰かが私を助けてくれれば幸いです。私のファイルは同じにあります

R2OpenBUGS を使っている人はいますか? r2winbugs を使用する必要がありますか? ...

私は、(単一の) 中間 (3 か月) の転帰を持つ患者のサンプルに対して、最終 (2 年間) の治療転帰 (成功、死亡、デフォルトまたは失敗など) をモデル化しようとしています。

R2OpenBUGS は多項式ノードで奇妙な事後分布を示しています。結果のうち 2 つは一定で、他の 2 つの結果は等しく、結果の総数はコホート サイズよりも大きくなっています。

私は何を間違っていますか？よろしくお願いします！コードアウトの出力は以下のとおりです。

出力の一部を次に示します。

rjagsR ライブラリを使用しています。関数coda.samplesは、mcmc.listたとえば (からexample(coda.samples)) を生成します。

ただし、オブジェクトを入力としてplot.bugs必要とする関数を使用したいと思います。bugs

mcmc.listオブジェクトをからbugsオブジェクトに変換することは可能plot.bugs(LINE.out)ですか?

1 か月以上未回答のstats.SE に関する同様の質問があることに注意してください。その質問には報奨金があり、2012 年 8 月 29 日に終了しました。

その他のヒント:

R2WinBUGS パッケージに "as.bugs.array" 関数があることを発見しましたが、この関数を mcmc.list に適用する方法が明確ではありません。

「R2Winbugs」パッケージを使用して、R でベイジアン分析を行っています。R 出力では、推定パラメーターが丸められます。私の質問は、どうすればこれを制御できますか?

これが私のコマンドです：

ここに私の出力の一部があります

私が望むのは、beta1 を 0.30000 として表示することです。digits=5 を試してみましたが、うまくいきませんでした。何をすべきか知っている人はいますか？

どんな助けでも大歓迎です

乾杯

グナル

私は R 2.15.2 の Win-7 ユーザーです

次のモデルが単純なロジット モデルの推定値にうまく収束しないのはなぜですか?

編集済み

尤度推定のモデル コード:

R2WinBUGS を使用したモデル コード: (編集済み)

そして、従来の方法で/ベイジアン法を使用せずに

以前の変更モデルまたは構文の異なるモデルを使用する必要があるかどうかを提案してください...

R から winBUGS を呼び出して、ロジスティック回帰を推定しようとしています。私は次のコードを使用しています：

WinBUGS は開きますが、エラーが発生します。「互換性のないコピー」。私はこのエラーを探していて、多くの人がこれを経験したことがわかりました. 私は提案された解決策を試しました.model.bugのファイルパスを短くする、作業ディレクトリに「データ」ファイルがこれ以上ないことを確認する、データ変数名の名前のリストを使用するなど.

しかし、私はまだこれを機能させることができません。私は何を間違っていますか？