#R code for Fig. 4 in lecture 18
plot(0:20-.1, dbinom(0:20,size=76, p=0.003203005), type='h', col='skyblue', lwd=4, xlab='x', ylab=expression(P(X==x)))
points(0:20+.1, dbinom(0:20,size=76, p=0.091849523), type='h', col='pink', lwd=4)
abline(v=2,lty=2)
arrows(2,.27,4,.27, code=2, angle=45, length=.07, col=4, lwd=2)
text(2,.31, expression(P(X >= 2) == .025), cex=1, col=4, pos=4)
arrows(0,.15,2,.15,code=1,angle=45, length=.07, col=2, lwd=2)
text(.1,.18,expression(P(X <= 2) == .025), cex=1, col=2, pos=4)
legend('topright',c(expression(paste("X ~ binomial(",n==76,", ",p[L],")")), expression(paste("X ~ binomial(",n==76,", ",p[U],")"))), col=c('skyblue','pink'), pch=15, pt.cex=1.1, cex=1, bty='n', text.col=c(4,2))