x<-seq(0,1,length=100)
par(mfrow=c(2,2),cex.main=3)
plot(x,dbeta(x,3,1),type="l",lwd=4,ylim=c(0,3))
lines(x,dbeta(x,7,7),lwd=4,col="red")
abline(v=.4,col="blue",lwd=4)
title(expression(paste(alpha==3,"and", beta==1)))
legend(0,3,legend=c("prior","mle","post"),col=c("black","blue","red"),lwd=4)
plot(x,dbeta(x,2,2),type="l",lwd=4,ylim=c(0,3))
lines(x,dbeta(x,6,8),lwd=4,col="red")
abline(v=.4,col="blue",lwd=4)
title(expression(paste(alpha==2,"and", beta==2)))
legend(0,3,legend=c("prior","mle","post"),col=c("black","blue","red"),lwd=4)
plot(x,dbeta(x,75,25),type="l",lwd=4,ylim=c(0,3))
lines(x,dbeta(x,79,31),lwd=4,col="red")
abline(v=.4,col="blue",lwd=4)
title(expression(paste(alpha==75,"and", beta==25)))
legend(0,3,legend=c("prior","mle","post"),col=c("black","blue","red"),lwd=4)
plot(x,dbeta(x,1,1),type="l",lwd=4,ylim=c(0,3))
lines(x,dbeta(x,5,7),lwd=4,col="red")
abline(v=.4,col="blue",lwd=4)
title(expression(paste(alpha==1,"and", beta==1)))
legend(0,3,legend=c("prior","mle","post"),col=c("black","blue","red"),lwd=4)

#par(mfrow=c(1,1))
#plot(x,dbeta(x,2,2),type="l",lwd=4,ylim=c(0,3))
#lines(x,dbeta(x,6,8),lwd=4,col="red")
#abline(v=.4,col="blue",lwd=4)
#title(expression(paste(alpha==2,"and", beta==2)))