#
# R2WinBugs example.  Needs the R2WinBUGS package installed
#
library(R2WinBUGS)
n <- 12
y <- 9
mu <- 0.5
prec <- 1
npred <- 10

data <- list("y","n","mu","prec","npred")
inits <- function(){list(theta=0.5,ypred=5)}
#
# You need to change the directories to the directories where you have the .bug file installed and WinBUGS installed 
# respectively.
#
binomialnc.sim <- bugs(data,inits,model.file="d:/teaching/bayes/doctoral/2014/binomialnc.bug",
                    parameters=c("theta","ypred"),n.chains = 1,n.iter=10100,n.burnin=100,n.thin=1,
                    bugs.directory="C:/Program Files/winbugs14/WinBUGS14",codaPkg=FALSE)
print(binomialnc.sim)
theta = binomialnc.sim$sims.array[,1,1]
mean(theta)
hist(theta,freq=FALSE,xlab='theta',ylab='f',xlim=c(0,1),main='')
ftheta <- density(theta,from=0,to=1)
lines(ftheta,lwd=5,col='red')

ypred = binomialnc.sim$sims.array[,1,2]
freq <- table(ypred)/length(ypred)
plot(freq,xlab='y',ylab='p',col='red')