y=c(15,11,12,5,12)
lambda=50
a=1
b=1
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)
}
sink("DensidadeMarginal.txt")
cat("
model{
  phi~dbeta(a,b)
  for(i in 1:n){
    N[i]~dpois(lambda)
    y[i]~dbin(phi,N[i])
  }
}",fill=T)
sink()
library(R2OpenBUGS)
fit=bugs(data,inits,params,model.file="DensidadeMarginal.txt",debug=F,n.chains=1,n.iter=100000,n.burnin=50000,n.thin=5, OpenBUGS.pgm='/home/dvd/bugs/bin/OpenBUGS')
phi=fit$sims.matrix[,"phi"]
 

eT1jKDE1LDExLDEyLDUsMTIpCmxhbWJkYT01MAphPTEKYj0xCmRhdGE9bGlzdChuPWxlbmd0aCh5KSx5PXksbGFtYmRhPWxhbWJkYSxhPWEsYj1iKQpwYXJhbXM9bGlzdCgicGhpIiwiTiIpCmluaXRzPWZ1bmN0aW9uKCl7CiAgcGhpPXJiZXRhKDEsYSxiKQogIE49cnBvaXMobGVuZ3RoKHkpLGxhbWJkYSkKICBsaXN0KHBoaT1waGksTj1OKQp9CnNpbmsoIkRlbnNpZGFkZU1hcmdpbmFsLnR4dCIpCmNhdCgiCm1vZGVsewogIHBoaX5kYmV0YShhLGIpCiAgZm9yKGkgaW4gMTpuKXsKICAgIE5baV1+ZHBvaXMobGFtYmRhKQogICAgeVtpXX5kYmluKHBoaSxOW2ldKQogIH0KfSIsZmlsbD1UKQpzaW5rKCkKbGlicmFyeShSMk9wZW5CVUdTKQpmaXQ9YnVncyhkYXRhLGluaXRzLHBhcmFtcyxtb2RlbC5maWxlPSJEZW5zaWRhZGVNYXJnaW5hbC50eHQiLGRlYnVnPUYsbi5jaGFpbnM9MSxuLml0ZXI9MTAwMDAwLG4uYnVybmluPTUwMDAwLG4udGhpbj01LCBPcGVuQlVHUy5wZ209Jy9ob21lL2R2ZC9idWdzL2Jpbi9PcGVuQlVHUycpCnBoaT1maXQkc2ltcy5tYXRyaXhbLCJwaGkiXQo=