問題タブ [jags]

JAGS を使用して、ユニット固有の時間傾向を含むモデルを推定しようとしています。ただし、問題は、これをモデル化する方法がわからないことであり、これまでのところ解決策を見つけることができませんでした。

例として、次のデータがあるとします。

通常の線形回帰を使用して、モデルを次のように推定します。

factor(ccode)*yearユニット固有の時間トレンドはどこにありますか。次に、JAGS を使用してモデルを推定します。そこで、インデックス作成用のパラメーターを作成します。

モデルを推定し、

次のモデルファイルを使用して、これが現時点でエラーが発生している場所です。

b0コードを使用したときにエラーが発生したことを考えると、定義の方法とインデックス付けが間違っていると確信していますCompilation error on line 7. Dimension mismatch taking subset of b0。しかし、これを解決する方法がわからないので、ここで誰かがこれについて提案を持っているかどうか疑問に思いましたか?

ifステートメントを使用せずにこの for ループに入らないようにする方法を見つけようとしています(JAGS にはありません)。

つまり、このループが逆の順序で実行されないようにする必要があります。つまり、Je[i] - 1が 2 未満の場合です。

の場合Je[i] - 1 == 1、ループは引き続き実行されますが、逆の順序で実行されるため、j最初は 2 に等しく、2 回目の反復jでは 1 に等しくなります。

これを修正する方法を知っている人はいますか？

と を使用%dopar%しforeach、出力を結合する必要があります。

並行して呼び出される関数には、呼び出しごとに一定の長さのリストが出力されます。ただし、このリストの要素の長さは常に一定ではありません。

結合後、各要素の元のリスト (反復) を識別できるようにしながら、結果をできるだけ単純化したいと考えています。

B/c の長さは、この深いレベルでは異なります。

データの例を次に示します。

do.call(Map, c(c, list(list1, list2, list3)))与えます：

そしてdo.call(Map, c(list, list(list1, list2, list3)))与えます：

編集、正解は次のようにする必要があります（RNGを許してください）：

また、リストの要素は必ずしも数値ではないことも指摘しておく必要があります。たとえば、 からのモデル出力である可能性がありますjags()。

Usingcは最初の部分を正しく (私が望むものと比較して) 取得し、 usinglistは最後の部分を正しく取得します。両方の長所を活かすにはどうすればよいですか?

rstan を介して Stan を学習しようとしています (R に精通しているため)。シンプルな混合パレート モデルとノーマル モデルを実行してみました。（私が知る限り）正常にコンパイルされますが、サンプリングに失敗し、エラーが発生します：

"(-2, 2) 間の初期化は、100 回試行した後に失敗しました。初期値を指定するか、制約された値の範囲を縮小するか、モデルを再パラメータ化してください。

サンプラーの呼び出し中にエラーが発生しました。サンプリングは行われていません」

パラメーター化するさまざまな方法を試し、初期値を設定しようとしましたが、すべて役に立ちませんでした。

私のR + rstanコードは以下の通りです：

この例は JAGS で問題なく動作し (そのため、JAGS にもタグを付けました)、そのコードを投稿できます。

ちなみに、パレート分布を追加の正規分布に変更すると、正常に動作します (もちろん、ナンセンスな答えが得られます)。

私が間違っていることについての提案は大歓迎です! どういうわけか、スタンではなくJAGSを考えているのではないかと心配していますが、パレートモデルをスタンに当てはめる人々の例を見つけることができなかったので、私のアプローチを相互検証するのは困難でした.

Rにパッケージ「rjags」をインストールしようとしていますが、次のエラーが発生しています:

Mac (Mavericks を実行) に JAGS をインストールし、R-3.1.2 を実行しています。このサイトdevtoolsで推奨されているように、パッケージを使用してURLからインストールしようとしました:

この他のサイトの手順も確認し、さまざまなバージョンの JAGS を試してみましたが、同じエラーが発生しました。次に何を試せばよいか途方に暮れていますが、解決策は簡単だと確信しています。誰にも洞察がありますか？本当にありがとう！

関数と同様jagsのR関数はありますrepか? 次のようなコードを使用して配列を作成したいと思います。

マニュアルを読みましたが、これを達成する方法を見つけることができませんでした。Stan がおそらくこれを許可することは理解していますが、離散パラメーターで推論を行う必要があるため、Stan を使用できませんでした。あなたの助けに本当に感謝します！

この質問は、JAGS ヘルプ フォーラムにも投稿されています。

COV既知の分布パラメーターからカバー率 ( ) の応答値を推定しようとしています。OpenBUGS で応答データを NA として指定することでこれを行うことができます (たとえば、以下のコード) が、JAGS はこれを許可しません。JAGSでこれを達成する方法を知っている人はいますか?

これは「教師なし統計学習」の範疇に入ると思います

R で coda (mcmc.list) オブジェクトを保存する方法がわかりません。他の人から同様の質問がありましたが、回答が特に明確ではないことがわかりました。理想的には、coda オブジェクトを R.data ファイルまたはテキスト ファイル (csv など) として保存し、モデルを再実行することなく再インポートして JAGS チェーンを分析できるようにしたいと考えています (これには約 30 分かかります)。私のコンピューター）。現在、私の coda オブジェクト "coda.samples" は次のようになっています。

ご覧のとおり、これは 1094 個のパラメーターの 3334 個の推定値を含む 3 つのマトリックスのリストです (つまり、長さ 3334 の 3 つのチェーン)。この coda オブジェクトを保存して、毎回モデルを再実行することなく R にコールバックできるようにしたいと考えています。また、独自のチェーンが 3 つあるという事実も維持したいと考えています。

次のコード (Gelman および Hill's Book に記載されているコードから適応) を使用して、Jags の可変係数/切片順序付きプロビット モデルを推定しようとしています。ただし、「初期化時に観察されていない親と矛盾する観察されたノード。適切な初期値を設定してください」というメッセージが表示されます。どこが間違っていますか？誰か助けてくれませんか？前もって感謝します ！！