library("arm") library(R2WinBUGS) J=2 y=c(1125,1125) prob744.data=list("J","y") prob744.inits=function() list (lambda1 = dexp(0.01)) prob744.parameters=c("lambda1") prob744.sim <- bugs (prob744.data, prob744.inits, prob744.parameters, "prob744-bug.txt", n.chains=3, n.iter=1000) attach.bugs(prob744.sim) hist(lambda1)