Was using xeval for er, now only x and er, reverted to data quantiles
for er but bounded off min(x) and max(x) by
sqrt(.Machine$double.eps). I believe that for long tailed
distributions with sharp peaks there was insufficient local
information.