xx <- seq(0, 1, len=199)
x <- 31
n <- 39
plot(xx, dbinom(x, prob=xx, size=n), type="l", col="blue", lwd=3, xlab=expression(theta), ylab=expression(L(theta*"|"*x)))
abline(v=31/39)
text(31/39, 0.17, label=expression(hat(theta)), xpd=TRUE)