text: Heres how these plots are done (Ive placed the method below in a wrapper function). The mle2 function used here is from Ben Bolkers bbmle package (I suspect that the method could be modified to use the more common mle function).