| 1 |
justin vas escriure el dia dt, 09 abr 2013: |
| 2 |
|
| 3 |
> Could you provide me a sample script which should run jags in parallel? |
| 4 |
|
| 5 |
This file generates data and samples (call it desc-000.bug) |
| 6 |
|
| 7 |
-----8<--------------- |
| 8 |
data { |
| 9 |
Alpha <- 3 |
| 10 |
Beta <- 5 |
| 11 |
Sigma <- 0.5 |
| 12 |
N <- 500 |
| 13 |
for (i in 1:500) { |
| 14 |
epsilon[i] ~ dnorm(0, 1/Sigma^2) |
| 15 |
x[i] ~ dunif(0, 10) |
| 16 |
y[i] <- Alpha + Beta*x[i] + epsilon[i] |
| 17 |
} |
| 18 |
} |
| 19 |
model { |
| 20 |
for (i in 1:500) { |
| 21 |
y[i] ~ dnorm(mu[i], tau) |
| 22 |
mu[i] <- alpha + beta * x[i] |
| 23 |
} |
| 24 |
alpha ~ dnorm(0, 0.001) |
| 25 |
beta ~ dnorm(0, 0.001) |
| 26 |
tau <- pow(sigma, -2) |
| 27 |
sigma ~ dgamma(0.001, 0.001) |
| 28 |
} |
| 29 |
-----8<--------------- |
| 30 |
|
| 31 |
The sample script with 3 chains (m-000.jmd) that calls the previous file. |
| 32 |
-----8<--------------- |
| 33 |
model in "desc-000.bug" |
| 34 |
compile, nchains(3) |
| 35 |
initialize |
| 36 |
update 10000 |
| 37 |
monitor alpha |
| 38 |
monitor beta |
| 39 |
monitor sigma |
| 40 |
update 10000 |
| 41 |
coda *, stem(coda_m-000-) |
| 42 |
parameters to "outs-000" |
| 43 |
exit |
| 44 |
-----8<--------------- |
| 45 |
|
| 46 |
HTH, |
| 47 |
|
| 48 |
-- |
| 49 |
- Xavier - |
Archives