問題タブ [winbugs]

マルチグレイントピックモデルをモデル化するための Winbugs/Jags モデルを作成しようとしています (まさにこの論文 -> http://www.ryanmcd.com/papers/mg_lda.pdf )

ここでは、特定の値に基づいて別の分布を選択したいと思います。例：次のようなことをしたい

これは Winbugs/JAGS で実行できますか?

不均衡なデータを使用した入力需要のモデルがあります。従属変数はyijtであり、iは生産関数の入力を(i=1,2,3)意味し、j企業は企業(j=1,..,21)をt意味し、時間を意味する(ti=1,..,Ti)。では、パネルが 3 つの方程式でバランスが取れていないため、winbug で使用する正しい形式は何でしょうか?

インデックスが少し違うSTATAロングフォーマット、3*107の形で使ってみました。行 ( i) は入力を意味し、列 ( j) はパネルを形成しi*tます。したがって、明示的な時間インデックスはありません。列にネストされたインデックスを使用して、データをパネルとして形成しました。ちょっとしたコード:

ここa[i]で、 は入力固有の自由パラメータ、uは企業固有の望ましい誤差成分であり、企業全体で推定する必要があります (21)。このために、ネストされたインデックスを使用します。モデル内の u の式は列u[firm[j]]よりも優先され、 foruは会社よりも優先されますu[k]。

これは、ネストされたインデックス作成に関する質問です。モデルに時間依存の多くの共変量が含まれている場合、これらすべての共変量に対してネストされたインデックスを使用する必要があります。そうしないと、Winbugs はそれをパネルとして取りませんか?

すべての多次元変数に多次元配列を使用しようとしました。バランスをとるために、多くの NA が含まれていました。たとえば、従属変数 yijt はy[ , , ](list コマンドを使用して) 3 次元配列ですが、何らかの理由で Wi​​nbugs は NA を認識しませんでしたか? 主なポイントは、事後分布が収束しないことであり、データのフォーマットが理由の 1 つになる可能性があります

WinBUGSとから使用しようとしましたRが、コードは次のとおりです。BRugsR2WinBUGS

結果が出ませんでした、そして私はの一部を見つけlog.txtますWinBUGS:

であることは明らかですerrorがmultiple definitions of node y[1]、それはどういう意味ですか? y[1]を使用しているので、複数の定義があるとは思いませy[i]んyが、loop.

コードは次のとおりです。

理由はわかりませんが、 を に置き換えるtheta[j] <- b*x[j-1]とプログラムは正しくなりますtheta[j] <- 1*x[j-1]が、 を定義しb ~ dunif(-1, 1)ました。確かに、最終モデルを設定する必要がありますが、追加しようとtheta[j] <- a - b*x[j-1]すると間違っていることがわかります。誰でも問題がどこにあるかを見つけますか?ab

Excel から WinBugs にデータを読み込む必要がありますが、何らかの理由で Wi​​nBugs が認識しません。私が試したことがたくさんあります (すべて以下にリストされています) が、どれも機能しません。

1. Excel からコピーし、「形式を選択して貼り付け」を「プレーン テキスト」として WinBugs にコピーしようとしました。[]すべての変数の後、さらには最後に追加しまし'END'た。エラーメッセージは次のようになりました。

   LoadNumeric Dataモジュールの手順で問題が発生しましたBugRectData

2. コピーしてメモ帳に貼り付け、[]すべての変数の後に追加'END'し、最後に追加してから、行間のすべてのスペースと、各データ間のスペースを 1 つだけ削除しました。それを Winbugs にドラッグすると、行間に 2 つのスペースが分散し、データ間に大きなギャップが生じ、(データをロードしようとしたときに) エラー メッセージが表示されました。

   LoadNumeric Dataモジュールの手順で問題が発生しましたBugRectData

3. 私が最後に行ったのは、.txt ファイルを Winbugs で開くことでした。すべてがごちゃごちゃになり、データをロードしようとすると、代わりに、NA2.0あるべきデータが強調表示NA 2.0され、「数値または NA または終わり"。これらの間違いを修正した後、データは最後まで問題ないように見え'END'、同じことを で強調表示するだけでしEたEND。

誰かが洞察を与えることができれば、私は最も感謝しています。

WINBUGS コード (v. 1.4.3) のエラーを見つけるのに助けが必要です。

「モデル仕様」ステップでは、モデルは構文的に正しいように見えます。ただし、データをロードしようとすると、次のエラーが発生しました。

配列インデックスが phi3 の配列上限を超えています

誰か助けてくれませんか？私のモデルは以下に提供されています：

pnormRで標準正規分布に使用します。ベイジアン プロビットを使用するため、Winbugs でこれを行うためにどのコードを使用できるか疑問に思っています。

たとえば、p = Phi(Xbeta) を pnorm(Xbeta) として計算します。

私の質問は、Winbugs でこれを行う方法です。

BUGS モデルの問題が見つかりません。モデルの構文は正しく、データは正常にロードされていますが、init のロード中に次のエラー メッセージが表示されます。

ベクトルの長さと等しくない初期値の数

感謝のために私のコードを以下に示します。

ではごきげんよう！

• カミラ

  /li>

Would someone please let me know what the equals function in WinBUGS return if one of the elements in the equals brackets is NA?

Specifically, I want to do the transformation:

and want to preserve NA's, so records with NA in oldvar[i] also has NA in newvar[i], but I am afraid that this tranformation turns NA's into 0's.

Any and all advice is much appreciated.

Trang