justin vas escriure el dia dt, 09 abr 2013:
> Could you provide me a sample script which should run jags in parallel?
This file generates data and samples (call it desc-000.bug)
-----8<---------------
data {
Alpha <- 3
Beta <- 5
Sigma <- 0.5
N <- 500
for (i in 1:500) {
epsilon[i] ~ dnorm(0, 1/Sigma^2)
x[i] ~ dunif(0, 10)
y[i] <- Alpha + Beta*x[i] + epsilon[i]
}
}
model {
for (i in 1:500) {
y[i] ~ dnorm(mu[i], tau)
mu[i] <- alpha + beta * x[i]
}
alpha ~ dnorm(0, 0.001)
beta ~ dnorm(0, 0.001)
tau <- pow(sigma, -2)
sigma ~ dgamma(0.001, 0.001)
}
-----8<---------------
The sample script with 3 chains (m-000.jmd) that calls the previous file.
-----8<---------------
model in "desc-000.bug"
compile, nchains(3)
initialize
update 10000
monitor alpha
monitor beta
monitor sigma
update 10000
coda *, stem(coda_m-000-)
parameters to "outs-000"
exit
-----8<---------------
HTH,
--
- Xavier -
