R code used to produce the graph shown below.

myplot(summary(res), lcol="blue", level=0.95)