fishSurvival.Bayes = function() { # assign data values y=c(15,11,12,5,12) lambda=50 a=1 b=1 # model specification in WinBUGS modelFilename = 'fishSurvival.txt' cat(' model { phi ~ dbeta(a,b) for (i in 1:n) { N[i] ~ dpois(lambda) y[i] ~ dbin(phi, N[i]) } } } ', fill=TRUE, file=modelFilename) # arguments for bugs() data = list(n=length(y), y=y, lambda=lambda, a=a, b=b) params = list('phi','N') inits = function() { phi = rbeta(1,a,b) N = rpois(length(y),lambda) list(phi=phi, N=N) } # call to bugs() library(R2WinBUGS) fit = bugs(data, inits, params, model.file=modelFilename, n.chains=1, n.iter=100000, n.burnin=50000, n.thin=5, bugs.seed=sample(1:9999,size=1), debug=FALSE, DIC=FALSE) phiPosterior = fit$sims.matrix[,'phi'] summary(phiPosterior) }