問題タブ [hierarchical-bayesian]

以下のようなデータセットがあります (最初の 5 行が表示されています)。CPA は、さまざまな広告フライトでの実験 (処理) から観察された結果です。フライトはキャンペーンで階層的にグループ化されます。

次のようなモデルに適合します。

私の知る限り、このモデルは単純に次のように述べています。

• 治療用スロープをご用意しております
• グローバルインターセプトがあります
• キャンペーンごとにインターセプトもあります

そのため、そのような変数ごとに1 つの 事後分布が得られます。

ただし、以下の結果を見ると、次の変数の事後分布も得られます1|campaign_uid_offset。それは何を表していますか？

モデルとプロットを適合させるためのコード:

JAGS+R で打ち切りデータのカウント データ用に階層的なポアソン モデルをコーディングしています。

A は行列で、行はさまざまな場所、列は雨の日を数えるさまざまな時間間隔です。共変量として、一連の X_k$ 行列があります。すべての場所と時間間隔について、平均気温 (X_1 に保存)、風の強い日の数 (X_2)、平均湿度 (X_3) などの共変量があります。

検閲されたデータを検閲されていないデータから分離する必要がありますか? マトリックスでそれを行うにはどうすればよいですか？

ご協力いただきありがとうございます！

アップデート：

私はこれを最後から2番目の登録された時間間隔までのサイクルで持っています（機器の故障による検閲があるかもしれない場所ごとに）

mu[i,j]<- a[1]*x[i,1]+a[2]*x[i,2]+ b[1,j]*varx[i,j]+b[2, j]*varx[i,j]
N[i,j] <- dpois(lambda[i,j])
log(lambda[i,j]) <- mu[i,j] + alpha[i]
alpha[ i]<- G0[潜在[i]]
潜在[i]~ dcat(確率[])

最後の時間間隔でこれを繰り返しました。
最後の時間間隔ビットにcensoor[i]<-step(-censored[i])を追加しました。censored[i]
は、機器に障害があるかどうかを示す vec です。

私は新しいのですが、うまくいきません。ありがとう

MCMCglmm で混合モデルを実行する際に問題が発生しています。線形回帰として実行しようとすると、うまく機能します。

しかし、ロジスティック回帰または順序回帰を実行するようには見えません。

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

sort.list(y) のエラー: 'x' は 'sort.list' のアトミックでなければなりません リストで 'sort' を呼び出しましたか? さらに: 警告メッセージ: 不明または初期化されていない列: '家族'。

前もって感謝します！