# First calculate sum of deaths and sum of waiting times (a.k.a. "exposed to
# risk" in actuarial terms)
d<-3
v<-2.75

# Now calculate range of values over which to plot log-likelihood function:
mu<-seq(0.05, 4, 0.01)

# Now calculate log-likelihood function
logL<--mu*v+d*log(mu)

# Plot
par(mfrow=c(1, 2), mar=c(4, 4, 2, 0))
plot(x=mu, y=logL, type="n", bty="n", font.lab=2, font.axis=2, main="Log-likelihood for Example 5.1", cex=0.8)
lines(x=mu, y=logL, lwd=2)
lines(x=c(d/v, d/v), y=c(min(logL), max(logL)), lty=3)
text(x=d/v+.1, y=min(logL), expression(paste(hat(mu), "=1.09")))

plot(x=mu, y=logL, type="n", bty="n", font.lab=2, font.axis=2, main="...with 2-support interval")
lines(x=mu, y=logL, lwd=2)
lines(x=c(d/v, d/v), y=c(min(logL), max(logL)), lty=3)
text(x=d/v+.1, y=min(logL), expression(paste(hat(mu), "=1.09")))

lines(x=c(0.26, 2.88), y=c(max(logL)-2, max(logL)-2))
lines(x=c(0.26, 0.26), y=c(min(logL), max(logL)-2), lty=3)
lines(x=c(2.88, 2.88), y=c(min(logL), max(logL)-2), lty=3)