我正在编写JAGS + R中的分层泊松模型,用于检查数据的计数数据。
A是一个矩阵,行是不同的地方,列是我计算下雨天的不同时间间隔。作为协变量,我有一组X_k $矩阵。对于每个地点和时间间隔,我都有协变量,如平均温度(存储在X_1中),刮风天数(X_2),平均湿度(X_3)。
我应该将审查后的数据与未审查的数据分开吗?我如何用矩阵做到这一点?
感谢您的帮助!
更新
我有一个循环直到最后但是一个注册的时间间隔(对于每个地方可能由于设备故障而进行审查)
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 [latent [i]]
潜伏[i] ~dcat(prob [])
在最后一个时间间隔重复此操作。在最后一个时间间隔位中添加了
censoor [i]< -step(-censored [i])
,被删除[i]是指示设备是否存在故障的vec
我是新人&它不起作用,有什么帮助吗?感谢