## JAGS
require(R2jags); invisible(runif(1))
Data <- list(y=31, n=39)
model <- function() {
  theta ~ unif(0,1)
  y ~ dbin(theta, n)
}
fit <- jags(Data, model=model, param="theta", n.iter=10000, n.chain=1, n.thin=1, n.burn=0)
fit
plot(fit)

## coda; more later
mcfit <- as.mcmc(fit)
summary(mcfit)

## BUGS
require(R2OpenBUGS)
fit <- bugs(Data, inits=NULL, model=model, param=c("theta"), n.iter=10000, n.chains=1, n.burn=0, DIC=FALSE)
print(fit, digits=3)
mcfit <- as.mcmc.list(fit)
summary(mcfit)

## Output to a folder
fit <- bugs(Data, inits=NULL, model=model, param="theta", n.iter=5, working.directory="premature")