#  R CODE TO DRAW POOLED DATA LIKELIHOOD FROM THE JANICAK DATA,
#  PLUS VERTICAL BARS AT EACH INDIVIDUAL STUDY OBSERVATION

par(lwd=1.5,cex=1.5,mar=c(5,4,1,1))


x<-seq(0,1,length=401)
phat<- c(1, .4, .69, .53, .36, .783, .9, .778, .67)
n <- c(20, 10, 16, 19, 14, 46, 10, 9, 6)

plot(x,dbeta(x,107,45),type="l",lwd=2,xlab="p",ylab="density",cex=1.5)
lines(phat,n/12, type = "h", col = "red", lwd=5)