#---------------------------------------------------- # Model involving a Poisson distribution with # a Gamma distributed 'lambda' # # GdA 23/2/2015 #---------------------------------------------------- library(rjags) #------------ model ---------------------------- model = "temp_model.bug" write(" model { lambda ~ dgamma(10, 5) x ~ dpois(lambda) y <- 2*x^2 } ", model) jm <- jags.model(model, n.chains = 3) # define the model chain <- coda.samples(jm, c('lambda', 'x', 'y'), n.iter=10000) # sampling print(summary(chain)) plot(chain)