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 です。
私は新しいのですが、うまくいきません。ありがとう