我想在我的带有4个内核的Windows计算机上并行运行jags模型,但是我们无法弄清楚为什么我的模型不会运行.我在网上广泛搜索了这些帖子:
http://andrewgelman.com/2011/07/23/parallel-jags-rngs/
http://users.soe.ucsc.edu/~draper/eBay-Google-2013-parallel-rjags-example.txt
当我运行一个简单的例子(见下面的代码)时%do%
,模型运行正常(当然是连续的).当我使用时%dopar%
,我收到错误:
Error in { : task 1 failed - "Symbol table is empty"
library(rjags)
library(coda)
library(foreach)
library(doParallel)
library(random)
load.module("lecuyer")
### Data generation
y <- rnorm(100)
n <- length(y)
win.data <- list(y=y, n=n)
# Define model
sink("model.txt")
cat("
model {
# Priors
mu ~ dnorm(0, 0.001)
tau <- 1 / (sigma * sigma)
sigma ~ dunif(0, 10)
# Likelihood
for (i in 1:n) {
y[i] ~ dnorm(mu, tau)
}
}
",fill=TRUE)
sink()
inits <- …
Run Code Online (Sandbox Code Playgroud)