問題タブ [winbugs]

WINBUGS から次の Panel データ モデルを推定しています。

データを長方形形式で読み込もうとすると、次のエラー メッセージが表示されます: expected comma (cursor next to y[]) :

エラーを修正するにはどうすればよいですか?

助けてください！

Winbugs でローカル レベルのモデルを定義したいと考えています。モデルは構文的に正しいです。しかし、実行すると、「ノード y[1] の複数の定義」というエラーが発生しました。

この原因と解決方法を教えてください。エラーを検索しましたが、それぞれに固有のケースがあり、これが発生する理由に関する一般的な解決策や説明はありません。

私は、世論調査の結果に基づいて多党派の主要な投票意図を推定するために、JAGS で階層的ディリクレ多項過程隠れマルコフ モデルを開発しようとしています。また、予備選挙の推定値を使用して、オーストラリアの優先投票システムの下での 2 党の優先投票シェアを計算します。

dmulti() 多項分布が実行時エラー メッセージで失敗しました: 適切なサンプラーが見つかりません。一連の二項分布と sum-to-N 制約を使用する回避策があります。理論的には、これで同じ結果が得られるはずですが、モデル内のスペースと時間が非効率になります。

私の質問は、多項分布を機能させるために、以下のモデルの隠れた一時的な部分で何かできることがあるかどうかです。

モデル (および周囲の R コード) は次のとおりです。

6 か月にわたるモデルの入力データは次のとおりです。

出力の比較 (私が持っている他のモデルとの比較) は次のとおりです。次のチャートの赤い線は、上記から生成されたものです。

対数尤度式にゼロトリックを使用した以下のコードを実行しようとしています (phi は対数尤度です):

データ付き:

モデルは構文的に正しく、データはロードされていますが、次のエラーが発生しました:

誰でもその問題を解決するのを手伝ってもらえますか?

OpenBUGS を使用してモデルを適合させようとしています。コードは次のとおりです。

しかし、このコードを実行すると、「コンマが必要です」というエラーが表示されます。誰かがこの問題を解決するのを手伝ってくれるなら、それは素晴らしいことです。

私は現在、いくつかの臨床試験にわたる生存データのメタ分析に取り組んでいます。

これを行うために、同じ方法論を使用して公開された分析からのコードがあります。ただし、公開された分析のデータを使用してこのコードを実行すると、結果を再現できません。実際、結果はいかなる種類の合理的な推定値にも収束しません。

コード自体 (データを含まない) は、作成者から直接提供されているため、正しいはずです。問題は初期値またはサンプリングの実行方法のパラメーターに関係していると思いますが、多くの初期値、バーンインの長さ、間引きなどで遊んだ後...意味のある結果が得られませんでした。

これを実行する方法（初期値など）について、適切に実行するための提案をいただければ幸いです。あるいは、コードに問題がある場合、またはデータがコードと一致しない方法で設定されている場合は、それを知っておくと役立ちます。

補足として、私は R2WinBUG を使用して分析を行っていますが、WinBUG のみを使用しても同様の問題が発生しました。

メソッドの背景を少し追加します。

これが機能する方法は、変量効果を使用して、複数の研究にわたる処理間の再パラメータ化されたワイブル分布の形状およびスケール パラメータの違いを推定することです。

ワイブル分布は、ハザード率の対数が a+b*log(t) になるように再パラメーター化されます。ここで、a はスケール パラメーター、b は形状パラメーターです。このことから、一定期間内に指定された数の患者から指定された数の失敗の尤度関数を計算できます。

残念ながら、この記事は公開されていますが、アクセスできる場合は、次のリンクを参照してください 。

モデルに入力された変数の簡単な要約:

NT: 含まれる個別の治療の数。

N: メイン データセットの行数。NS: 研究数

s: データの行が対応する研究 (これは 1:6 の番号が付けられています)

r: この治療/研究の間隔で失敗した患者の数

n：この治療/研究の間隔の開始時にリスクのある患者の数

t: この行のデータが対応する処理 (1:3 の番号)

b: どの処理が他の処理と比較されるベースライン処理であるかを示します (各行で 1 に設定)。

bs: この研究の対照群である治療

bt: 本研究の研究部門である治療

WinBUGS コード (データを含む):

私は WinBUGS/R を使用して 3 つのレベルを持つベイジアン階層モデルを開発していますが、適切なコード例を見つけるのに苦労しています。誰か提案してもらえますか？私は WinBUGS を初めて使用しますが、マルチレベル モデリングではなく、この言語でモデルを指定するのに助けが必要です。

地域内にネストされた地方自治体の 4 つの反復測定データ ポイントがあります。したがって、地方自治体レベルと地域レベルで 2 つのランダム インターセプトを指定します。

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

次の論文のスクリプトを使用しています (Zipkin, EF, Royle, JA, Dawson, DK, Bates, S., 2010. 保存および管理アクションの影響を評価するための複数種の発生モデル。Biological Conservation 143, 479- 484) 鳥種の占有率を推定します。検出推定の変数の 1 つ (以下のコードの K ループ) は、レベル 1 ～ 6 のカテゴリ変数である風です。OpenBUGS で関数を使用しようとしましたdcatが、これは一様な事前 (beta(1,1)) であることを望んでいますが、OpenBUGS はエラーで失敗します:

これは、行を削除b3[i] ~ dcat(p[i,])#WIND Dataしても発生しません。

dcat を適切に指定する方法、または WinBUGS/OpenBUGS のカテゴリ変数をコーディングする方法についてアドバイスをいただければ幸いです。

誰かが反復測定 ANOVA の BUGS/JAGS モデルのコードを持っているかどうか疑問に思っていましたか? 基本的に、Time of day、Day、およびTreatmentに対してモデル化したい応答 ( y ) があります。また、治療 x 時間帯と治療 x 日という 2 つの交互作用用語も含めたいと思います。. この研究には約 20 人が参加し、約 1 週間にわたって 1 日 4 回測定されました。どこから始めればよいかよくわかりません。時刻の共変量も日の共変量内にネストする必要があるのではないかと心配しています。BUGS/JAGS モデルの尤度部分のコードを誰かが持っていれば、大歓迎です。私は事前の世話をすることができます。これでうまくいくようには見えません。

作業モデルを適応させ、WinBUGSR を使用して軽減しようとしているときにR2WinBUGS、いくつかのエラー メッセージが表示されました。の仕様に関連していると思いますがinits、問題を解決できませんでした。

最初はエラーメッセージが使用されていましたinits1：

Uwe Liggers による提案された修正「リストを含むリスト ソリューション」を読んだ後、inits1 をinits2次のように変更しました。

エラーを受け取りました：

この投稿"getting-winbugs-leuk-example-to-work-from-r-using-r2winbugs" で AndyCによって提案された修正も試みました。を次のように変更initsします。

エラーを受け取りました：

これは私の試みR2WinBUGSです。注意してください、私はinits動作しなかったコードにいくつかを含めました：

これは、ブリストル大学とレスター大学が開催したコースから採用された作業WinBUGSコードです。WinBUGS