########## R script: binomBeta ##########

# For MCMC analysis of a simple binomial/beta
# situation using `BRugs'.

library(BRugs)
source("summariseMCMC.r")

allData <- list(n=20,y=c(0,0,1,0,1,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0))

parInit <- list(list(p=mean(allData$y)))

BRugsFit(data=allData,inits=parInit,parametersToSave="p",nBurnin=100,
         nIter=1000,modelFile="binomBetaModel.txt",numChains=1)

summariseMCMC("p","probab. of success")

########## End of binomBeta ##########