# rm(list=ls()) ## Clear workspace. NOT in vignette.
library(rjags) # ? Is this needed?
library(R2jags)
y<-c(5.308, 7.24, 9.638, 12.866, 17.069, 23.192, 31.443, 38.558, 50.156, 62.948,
         75.995, 91.972)
t<-1:12

 tfit_jags <- jags(model="../BUGS/weeds_bugs.txt",
      data=list(N=length(y),t=t, y=y),
      parameters.to.save=c("b1","b2","b3"),
      n.iter=10000,
      n.burnin=1000,
      n.thin=10,
      n.chains=3,
      inits=list(list(b1=1,b2=1,b3=1),list(b1=2,b2=2,b3=2),list(b1=0,b2=0,b3=0))
      )
# wpred <- .... stuff to create predicted model
cat("output of jags() run\n")
 tfit_jags

