The following is R code that would generate the plot in figure 1:

boxplot(count ~ spray, data = InsectSprays, col = 'lightgray')
rect(xleft=c(3,1), ytop=c(25,5),xright=c(4,2), ybottom=c(20,0))
title(main='A', cex=3)