library(R2OpenBUGS)
windows(record=T)
data <- list(k=10, Y=c(5,1,5,14,3,19,1,1,4,22),
      t=c(94.320,15.720,62.880,125.760,5.240,
      31.440,1.048,1.048,2.096,10.480))
inits <- function(){
    list(a=c(1,1,1,1,1,1,1,1,1,1), nu=2,S=2)
}
pumps.sim <- bugs(data, inits, model.file="pumpsmodel.txt",
    parameters=c("theta","S"), n.chains=3, n.iter=20000)
print(pumps.sim)
plot(pumps.sim)
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)