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

Passing existing matrix into JAGS model using R package rjags