# Return the model specification and starting values for a lognormal Poisson, then run the model using run.jags:
data <- rpois(100, rlnorm(3, 0.2))
strings <- run.model(model="LP", data=data, call.jags=FALSE)
modelstring <- strings[[1]]
datastring <- strings[[2]]
initial.values <- strings[[3]]
monitors <- strings[[4]]
run.jags(model=modelstring, inits=initial.values, check.conv=TRUE, data=datastring, monitor=monitors)
# This is equivalent to using:
run.model(model="LP", data=data, call.jags=TRUE, check.conv=TRUE)
Run the code above in your browser using DataLab