pumps.sim <- bugs(data, inits,
     model.file="pumpsmodel.txt",
     parameters=c("theta","S"),
     n.chains=3,n.iter=20000,
     codaPkg=T)
codaobject <- read.bugs(pumps.sim)
plot(codaobject)
codamenu()