截尾模型JAGS计数数据

时间:2017-07-11 13:28:01

标签: hierarchy bayesian poisson jags hierarchical-bayesian

我正在编写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

我是新人&它不起作用,有什么帮助吗?感谢

0 个答案:

没有答案