使用R包rjags将现有矩阵传递到JAGS模型中

Question

很简单的问题，但是我不知何故找不到解决方案。 如何使用rjags将现有值（矩阵，向量）传递到JAGS模型中？

这是一个示例代码：

model{
      A = inverse(B)
     }

我想以某种方式将B = diag(100)传递到上面的模型中，该怎么办？

Answer 1

您将数据传递给模型，并将data arg传递给jags.model ：

例如：

library(rjags)

M <- 'model {
  A <- inverse(B)
}'

j <- jags.model(textConnection(M), data=list(B=diag(10)), n.chains=3)
jags.samples(j, 'A', 1000)