Re: [gentoo-science] sci-libs/atlas from science overlay - is it multithreaded?

["Xavier Fernández i Marín" <xfim.ll@gmail.com>]

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  -