使用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)