#source("c:/Work07/ShortCourse07/Programs/PoissonCompletion.txt",print.eval=TRUE)#
#This does the gibbs sampler for the poisson completion#
nsim<-500;lam<-array(313/360,dim=c(nsim,1));y<-array(0,dim=c(13,1));
for (j in 2:nsim)  {
for(i in 1:13){while(y[i] < 4) y[i] <- rpois(1,lam[j -1])};a<-313+sum(y);
lam[j]<-rgamma(1,a,scale=1/360);
y<-y*0;
}
den<-1:(nsim)
meanlam<-cumsum(lam)/den;
par(mfrow=c(1,2))
plot(meanlam,type="l",ylim=c(.9,1.1),xlab="iteration",ylab="estimate",col="red")
hist(lam,main="Mean",freq=F,col="green")